You signed in with another tab or window. Reload to refresh your session.You signed out in another tab or window. Reload to refresh your session.You switched accounts on another tab or window. Reload to refresh your session.Dismiss alert
I'm attempting to import a 3D Tetrahedral Mesh (.vtu) and use it for a single-phase flow simulation. I've successfully simulated using an internal mesh and a provided sample (egg.vtu) without encountering any errors. However, when I import my own mesh, I encounter the following error:
***** ERROR
***** LOCATION: /home/sachdeva/GEOS/codes/GEOS/src/coreComponents/mesh/CellElementSubRegion.cpp:397
***** Controlling expression (should be false): m_elementVolume[k] <= 0.0 ***** Rank 0: Negative volume for element 6020549 type C3D4 in subregion domain/MeshBodies/carteisanMesh/meshLevels/Level0/ElementRegions/elementRegionsGroup/reservoir(Mesh1_co2Injection.xml,l.56)/elementSubRegions/tetrahedra
I've also tried using mesh_doctor.py to investigate the root cause, and it seems that the volume of each element in my mesh is negative, unlike the sample mesh.
I suspect that the source code doesn't account for irregular tetrahedra, where the centroid lies outside the element, causing the negative volume.
Nevertheless, I'm unable to import my desired reservoir geometry. Could you please provide insights if you've encountered a similar problem, or suggest ideas on how I could tackle this error?
reacted with thumbs up emoji reacted with thumbs down emoji reacted with laugh emoji reacted with hooray emoji reacted with confused emoji reacted with heart emoji reacted with rocket emoji reacted with eyes emoji
-
I'm attempting to import a 3D Tetrahedral Mesh (.vtu) and use it for a single-phase flow simulation. I've successfully simulated using an internal mesh and a provided sample (egg.vtu) without encountering any errors. However, when I import my own mesh, I encounter the following error:
***** ERROR
***** LOCATION: /home/sachdeva/GEOS/codes/GEOS/src/coreComponents/mesh/CellElementSubRegion.cpp:397
***** Controlling expression (should be false): m_elementVolume[k] <= 0.0 ***** Rank 0: Negative volume for element 6020549 type C3D4 in subregion domain/MeshBodies/carteisanMesh/meshLevels/Level0/ElementRegions/elementRegionsGroup/reservoir(Mesh1_co2Injection.xml,l.56)/elementSubRegions/tetrahedra
I've also tried using mesh_doctor.py to investigate the root cause, and it seems that the volume of each element in my mesh is negative, unlike the sample mesh.
I suspect that the source code doesn't account for irregular tetrahedra, where the centroid lies outside the element, causing the negative volume.
Nevertheless, I'm unable to import my desired reservoir geometry. Could you please provide insights if you've encountered a similar problem, or suggest ideas on how I could tackle this error?
Beta Was this translation helpful? Give feedback.
All reactions