search

OpenCMISS / iron

Source code repository for OpenCMISS-Iron
9 stars 62 forks source link

Navier Stokes depends on the mesh component ordering #76

Open chrispbradley opened 7 years ago

chrispbradley commented 7 years ago

For a transient NS problem it is common to have quadratic for velocity and linear for pressure. Currently NS works if the velocity (quadratic) mesh component is number 1 and the pressure (linear) mesh component is number 2. The NS code does not work, however, if you have the linear pressure component as mesh component 1 and the quadratic velocity component is mesh component 2. There is obviously some assumptions going on in the code which needs to be fixed.

dladd commented 7 years ago

I believe I have gotten to the bottom of this- the main culprit appears to have been an assumption of the mesh component ordering in the boundary face integration routines. The fix for this should come with the pull associated with issue #77