SimVascular / SimVascular/svMultiPhysics

Testing Taylor-Hood 10-noded tetrahedra

Open
#414 1 comment 0 reactions 0 assignees View on GitHub

Nobody has claimed this yet.

bug
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

Image

pipe_RCR_3d - tet10 - taylor-hood

Image

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

Open the contributing guide

First steps

  1. Read the whole issue, then the project's contributing guide.
  2. Comment on the issue to say you are picking it up — it saves two people doing the same work.
  3. Fork the repository and make your change on a branch.
  4. 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

Get new issues in your inbox

A short digest of beginner-friendly GitHub issues.