0

Questions Regarding Local Refinement and Mesh Size Limitations (number of elements)

1. Local Mesh Refinement

I am importing my mesh from ECLIPSE and not using the MeshMaker utility. I would like to apply local refinement, similar to the approach shown in the snippet below. However, when I write the refined mesh using the toughio library in Python, the connectivity between coarse and fine cells is lost, and these cells are treated as disconnected.

a. How can I apply local refinement and maintain a fully functional mesh without completely remeshing the domain?

b. Is there a tool or library that allows remeshing around a specific point or region while keeping the new cells connected to the rest of the original mesh?

 

2. TOUGH3 Large Mesh Size Limit

I am attempting to run TOUGH3 simulations with more than 500,000 elements, but the simulation fails to run—regardless of whether I use a PC with 64 GB RAM or with 512 GB RAM. Reducing the number of cells allows the simulation to run on both machines without any issues.

Is there a known limitation in TOUGH3 regarding the maximum number of grid blocks or connections?

If so, is this limitation imposed by the source code (e.g., fixed-size arrays) and can it be adjusted by modifying some parameters during compilation?

Any guidance you could provide on these issues would be greatly appreciated.

Kind regards,
Javad

2 replies

null
    • kenny
    • 2 days ago
    • Reported - view

    1. Local Mesh Refinement.

    (a) d1 and d2 are the distances from center of element 1 and 2 to their interface (perpendicular distance to the interface). Connection is the interface area shared by the two elements. 

    (b) MVIEW, developed by Geofirma,  is good for generating such a mesh.

    2. To my knowledge, there is no hard limitation on the element and connection number for TOUGH3.  Did you get any error messages?

      • Mohammadjavad_Shokriafra
      • 2 days ago
      • Reported - view

       Thank you very much for your response. I understand how the mesh works and I'd be able to manually add those values to the Mesh file for a small model; however it has been difficult to manually add those values for a bigger model. I will check the MVIEW, thanks for your suggestion.
      2. I did not get any error message in the output file, the code stopped running I think at shown below. 

      OUTPUT

Content aside

  • 2 days agoLast active
  • 2Replies
  • 35Views
  • 2 Following