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

HDG formulation for linear elasticity #11

Merged
merged 9 commits into from
Feb 16, 2022

Conversation

amartinhuertas
Copy link
Member

In this PR I am implementing the HDG formulation for linear elasticity with symmetric stresses and k+1-order and k+2-order, k>=1, convergence for stresses and displacements resp.

This is WIP, not ready to be merged.

stresses and k+1-order and k+2-order, k>=1, convergence for stresses and
displacements resp.
@amartinhuertas
Copy link
Member Author

cc @ruizbaier @santiagobadia FYI ... Work in Progress, keep you posted when we have a validated implementation for linear elasticity. Hopefully soon.

@codecov-commenter
Copy link

codecov-commenter commented Feb 7, 2022

Codecov Report

Merging #11 (bcff4ac) into main (28f8cd5) will decrease coverage by 1.70%.
The diff coverage is 45.31%.

Impacted file tree graph

@@            Coverage Diff             @@
##             main      #11      +/-   ##
==========================================
- Coverage   81.77%   80.07%   -1.71%     
==========================================
  Files          11       11              
  Lines        1317     1380      +63     
==========================================
+ Hits         1077     1105      +28     
- Misses        240      275      +35     
Impacted Files Coverage Δ
src/GridapAPIExtensions.jl 66.82% <0.00%> (-9.65%) ⬇️
src/Skeleton.jl 77.52% <0.00%> (-1.60%) ⬇️
src/GridapTmpModifications.jl 87.17% <100.00%> (+4.03%) ⬆️
src/HybridAffineFEOperators.jl 96.17% <100.00%> (ø)

Continue to review full report at Codecov.

Legend - Click here to learn more
Δ = absolute <relative> (impact), ø = not affected, ? = missing data
Powered by Codecov. Last update 28f8cd5...bcff4ac. Read the comment docs.

  L2-projection of the bulk displacements on the skeleton (TO-DO)

* Found issues related to some of the terms in the weak formulation.
  (TO-INVESTIGATE)
…p/GridapHybrid.jl into linear_elasticity_symmetric_stresses
symmetric stresses HDG formulation for linear elasticity runs bottom top
without fatal errors.

Next step manofacture exact solution in FE space and check
consistency/correctness of the current implementation
* HDG formulation with symmetric stresses for linear elasticity,
  still without P_M projection, validated for solution in FE space
  and analtyic solution from paper.

* Next step: implement P_M projection.
Surprise, surprise. It crashes. To DEBUG.
* Added necessary fixes to support triangle meshes.
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

Successfully merging this pull request may close these issues.

None yet

2 participants