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
3 replies
-
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?
-
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.
-
-
Dear Javad,
The mesh refinement you would like to use seems not to follow the requirements of IFDM implemented by TOUGH2/3. The code should be modified as described by Garcia and Pruess (2000, LBNL-45646).
It seems better to use the capabilities of the IFDM to build local grid refinements. The same topic was discussed in a previous post (How to create meshes with local grid refinements) a couple of years ago, which I think is still giving useful advices.
I see some issue in applying a local grid refinement within an existing grid without regridding from scratch. Local grid refinements can be applied in AMESH and in Petrasim by adding the additional nodes required, but the grid needs to be rebuilt.
Regards,
Alfredo