Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

[GeoMechanicsApplication] Remove swapping faces in condition process #12227

Merged
merged 24 commits into from
Mar 28, 2024

Conversation

rfaasse
Copy link
Contributor

@rfaasse rfaasse commented Mar 26, 2024

📝 Description
We have identified that certain 3D conditions in GeoMechanicsApp are not functioning correctly due to the presence of edge and face swapping. These lines were originally introduced as a workaround for a bug that has since been fixed. See issues #9818 #9422.

However, it has come to our attention that these workarounds are now causing issues with our 3D boundary conditions, therefore they are removed, test cases adjusted and new test cases added.

@rfaasse rfaasse self-assigned this Mar 26, 2024
@rfaasse rfaasse added the GeoMechanics Issues related to the GeoMechanicsApplication label Mar 26, 2024
@rfaasse rfaasse requested review from WPK4FEM and avdg81 March 28, 2024 10:00
@rfaasse rfaasse marked this pull request as ready for review March 28, 2024 10:00
Copy link
Contributor

@WPK4FEM WPK4FEM left a comment

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

And thank you for also updating the existing README with the added thermal tests.
Remaining comment is only on consistency in formatting of mentioned units.

- The X, Y and Z displacements in the bottom four nodes (1, 2, 3 and 4) are fixed.
- Material:
- The material is described using a linear elastic material with a `LinearElastic3DLaw`, a Young's modulus
of 3.0e7 kPa and a Poisson ratio of 0.0.
Copy link
Contributor

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Mostly we have put units between [ ].
See also line 19, 23 and 24

Copy link
Contributor Author

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Done

Copy link
Contributor

@avdg81 avdg81 left a comment

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

@mnabideltares and @rfaasse: thank you for taking out the work-around that was previously required. This helps to make the code more readable. And I'm pleased to see that you've added a few more tests. Great!

"kratos_module": "KratosMultiphysics.GeoMechanicsApplication",
"process_name": "ApplyScalarConstraintTableProcess",
"Parameters": {
"model_part_name": "PorousDomain.BottomFixed",
Copy link
Contributor

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Question: is it intentional that only the bottom nodes have prescribed water pressures? Or should the fixed value be applied to all nodes?

Copy link
Contributor Author

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

I think this should be enough (correct me if I'm wrong @WPK4FEM), since now we have a reference pressure (and also since ignore undrained is set to true)

Copy link
Contributor Author

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Changed to force the entire element to 0.0 water pressure

Copy link
Contributor Author

@rfaasse rfaasse left a comment

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

I processed/replied to all that I could.

- The X, Y and Z displacements in the bottom four nodes (1, 2, 3 and 4) are fixed.
- Material:
- The material is described using a linear elastic material with a `LinearElastic3DLaw`, a Young's modulus
of 3.0e7 kPa and a Poisson ratio of 0.0.
Copy link
Contributor Author

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Done

"kratos_module": "KratosMultiphysics.GeoMechanicsApplication",
"process_name": "ApplyScalarConstraintTableProcess",
"Parameters": {
"model_part_name": "PorousDomain.BottomFixed",
Copy link
Contributor Author

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

I think this should be enough (correct me if I'm wrong @WPK4FEM), since now we have a reference pressure (and also since ignore undrained is set to true)

WPK4FEM
WPK4FEM previously approved these changes Mar 28, 2024
Copy link
Contributor

@WPK4FEM WPK4FEM left a comment

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Thank you, I think this is good to go.

Copy link
Contributor

@avdg81 avdg81 left a comment

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Many thanks for the final changes. I believe we're good to go :-)

@rfaasse rfaasse enabled auto-merge (squash) March 28, 2024 16:04
@rfaasse rfaasse merged commit ca1560e into master Mar 28, 2024
10 of 11 checks passed
@rfaasse rfaasse deleted the geo/12026-remove-swapping-faces branch March 28, 2024 16:58
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
GeoMechanics Issues related to the GeoMechanicsApplication
Projects
None yet
Development

Successfully merging this pull request may close these issues.

[GeoMechanicsApplication] Remove swapping edges and faces from "FindNeighbourElementsOfConditionsProcess"
4 participants