SimVascular / SimVascular/svMultiPhysics
Testing Taylor-Hood 10-noded tetrahedra
Nobody has claimed this yet.
- Dominant language
- C++
- Stars
- 45
- Forks
- 60
- Avg merge
- 5d 23h
- Merged PRs (30d)
- 11
Description
Description
I am testing Taylor-Hood 10-noded tetrahedra for fluids problems.
Running some simple tests does not produce good results.
Reproduction
I've converted the pipe_RCR_3d mesh into 10-noded tetrahedra.
Running a simulation using the NS solver converges
NS 1-1 1.915e+00 [0 1.000e+00 1.000e+00 4.158e-04] [4 -11 77]
NS 1-2 3.734e+00 [-41 8.252e-03 8.252e-03 4.477e-04] [4 -9 76]
NS 1-3 5.581e+00 [-87 4.314e-05 4.314e-05 9.197e-04] [4 -10 76]
NS 1-4 7.418e+00 [-126 4.745e-07 4.745e-07 6.957e-04] [4 -10 76]
NS 1-5 9.310e+00 [-168 3.874e-09 3.874e-09 7.596e-04] [4 -10 74]
NS 2-1 1.007e+01 [0 1.000e+00 2.414e+02 3.866e-04] [0 -79 43]
NS 2-2 1.232e+01 [-67 4.079e-04 9.847e-02 6.880e-04] [5 -12 81]
NS 2-3 1.363e+01 [-89 3.429e-05 8.277e-03 9.008e-04] [2 -9 67]
NS 2-4 1.571e+01 [-142 7.758e-08 1.873e-05 6.966e-04] [5 -10 79]
NS 2-5s 1.707e+01 [-162 7.394e-09 1.785e-06 8.178e-04] [2 -12 64]
Velocity results seem to match svMultiPhysics/tests/cases/fluid/pipe_RCR_3d but pressure does not
pipe_RCR_3d - tet4
pipe_RCR_3d - tet10 - taylor-hood
Pressure results match from a simulation using tet10 without taylor-hood.
Expected behavior
I would expect taylor-hood results to match tet4 and tet10 simulations.
Additional context
No response
Code of Conduct
- I agree to follow this project's Code of Conduct and Contributing Guidelines
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.
Research direction
Start with the svMultiPhysics/tests/cases/fluid/pipe_RCR_3d case and compare the tet4, tet10, and Taylor-Hood tet10 simulations using the NS solver. Check why velocity agrees while pressure differs; the issue is done when Taylor-Hood pressure results match the tet4 and non-Taylor-Hood tet10 cases.
Written by the indexing model from the issue text.
Assessment
- Tech stack
- cpp
- Domain
- hpc
- Issue type
- Bug
- Difficulty
- 4/5
- Estimated time
- 3-5 days
- Activity status
- Stale
- Clarity
- Mostly clear
- Newbie friendliness
- 38/100