Higher order integration for Neumann BC
Open
@castelletto1 is already working on this.
Since Aug 7, 2020.
type: feature
- Dominant language
- C++
- Stars
- 287
- Forks
- 109
- Avg merge
- 4d 41m
- Merged PRs (30d)
- 5
Description
What is the requested feature?
Currently, we impose traction boundary condition as F*A/n, with n number of nodes for the specific face. We should use a proper integration. E.g.:
https://github.com/GEOSX/GEOSX/blob/ee9de4c0693462b2471567882aeb62f16d452c1c/src/coreComponents/physicsSolvers/solidMechanics/SolidMechanicsLagrangianFEM.cpp#L736
Is your request related to a specific problem?
For quadrilateral faces with generic angles, this is a rough approximation.
Describe the solution you'd like
To use a proper integration, from the new finite element library.
Contributor guide
First steps
- Read the whole issue, then the project's contributing guide.
- Comment on the issue to say you are picking it up — it saves two people doing the same work.
- Fork the repository and make your change on a branch.
- Open a pull request that references the issue number.
Assessment
This issue has not been assessed yet.