Feature/oduran/test fem stress consistency#3957
Draft
Conversation
…l surfaceGenerator
The previous UMFPACK_ORDERING_BEST setting triggered external METIS calls that caused stack corruption due to library mismatch.
…st_fem_stress_consistency
OmarDuran
force-pushed
the
feature/oduran/test_fem_stress_consistency
branch
from
850ac67 to
85fc4de
Compare
This reverts commit deef39c.
This was
linked to
issues
Feb 6, 2026
Update SolidMechanicsAugmentedLagrangianContact.cpp Update SolidMechanicsAugmentedLagrangianContact.cpp Update SolidMechanicsAugmentedLagrangianContact.cpp Update SolidMechanicsAugmentedLagrangianContact.cpp Update SolidMechanicsAugmentedLagrangianContact.cpp Update SolidMechanicsAugmentedLagrangianContact.cpp Update SolidMechanicsAugmentedLagrangianContact.cpp Update SolidMechanicsAugmentedLagrangianContact.cpp Update SolidMechanicsAugmentedLagrangianContact.cpp wip: solution Update SolidMechanicsAugmentedLagrangianContact.cpp
This file contains hidden or bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
Sign up for free
to join this conversation on GitHub.
Already have an account?
Sign in to comment
2 participants
Add this suggestion to a batch that can be applied as a single commit.This suggestion is invalid because no changes were made to the code.Suggestions cannot be applied while the pull request is closed.Suggestions cannot be applied while viewing a subset of changes.Only one suggestion per line can be applied in a batch.Add this suggestion to a batch that can be applied as a single commit.Applying suggestions on deleted lines is not supported.You must change the existing code in this line in order to create a valid suggestion.Outdated suggestions cannot be applied.This suggestion has been applied or marked resolved.Suggestions cannot be applied from pending reviews.Suggestions cannot be applied on multi-line comments.Suggestions cannot be applied while the pull request is queued to merge.Suggestion cannot be applied right now. Please check back later.
This PR adds a new integration test, to verify the consistency of the Finite Element Method implementation in GEOS. The test focuses on stress transmission across fractured domains using the SurfaceGenerator class.
The test ensures that the mechanical solvers preserve stress equilibrium and correctly handle contact/interface conditions across different mesh types and fracture configurations.
The physical setup is a pure compression case with no shear induced by the boundary conditions. As a result, the stress field is constant, which allows for direct verification: for volume elements, the averaged computed stress tensor is checked against the applied boundary tractions, and for fracture surfaces, the traction vector 𝑡 = 𝜎⋅𝑛 is projected from neighboring volume cells onto fracture faces, is computed and verified on both fracture sides against the expected traction.
The figures show (top-left) a structured hexahedral mesh, (top-right) a structured tetrahedral mesh, (bottom-left) a hexahedral mesh with curved fractures, and (bottom-right) a tetrahedral mesh with curved fractures.
Fracture leading orientations:
f1: normal aligned with the x-axis
f2: normal aligned with the y-axis
f3: normal aligned with the z-axis
The test covers all combinations of fractures:
f1 -meshes:
f2 -meshes:
f3 -meshes:
f1-f2 -meshes:
f2-f3 -meshes:
f1-f3 -meshes:
f1-f2-f3 -meshes: