This is a read only copy of the old FEniCS QA forum. Please visit the new QA forum to ask questions

Help for solving Navier Stokes inside two materials

–2 votes

Hi all,

I am currently have a geometry with two materials, one is wall and the other is some fluid. The gmsh file for the geometry is as following,

Point(1) = {0.25, 0.25, 0, 0.4};
Point(2) = {9.75, 0.25, 0, 0.4};
Point(3) = {5.5, 1.0, 0, 0.4};
Point(4) = {5.75, 1.0, 0, 0.4};
Point(5) = {0.25, 5.75, 0, 0.4};
Point(6) = {5.5, 5.75, 0, 0.4};
Point(7) = {5.75, 5.75, 0, 0.4};
Point(8) = {9.75, 5.75, 0, 0.4};
Point(9) = {0, 0, 0, 1};
Point(10) = {10, 0, 0, 1};
Point(11) = {0, 6, 0, 1};
Point(12) = {10, 6, 0, 1};
Line(1) = {5, 1};
Line(2) = {1, 2};
Line(3) = {2, 8};
Line(4) = {8, 7};
Line(5) = {7, 4};
Line(6) = {4, 3};
Line(7) = {3, 6};
Line(8) = {6, 5};
Line(11) = {11, 9};
Line(12) = {9, 10};
Line(13) = {10, 12};
Line(14) = {12, 11};
Line Loop(10) = {8, 1, 2, 3, 4, 5, 6, 7};
Plane Surface(10) = {10};
Line Loop(16) = {8, 1, 2, 3, 4, 5, 6, 7, -14, -13, -12, -11};
Plane Surface(16) = {16};
Physical Surface(0) = {10};
Physical Surface(1) = {10,16};

Where the Surface 10 is the fluid and Surface 16 is the wall.
I am trying to get the fluids' velocity and pressure from the Navier Stokes Equation. So firstly, I just assume the hole geometry satisfies the Navier Stokes Equation, then I set the both the pressure and velocity inside the walls to be 0 as a boundary condition for the equations.
However, the result I get is not good. Could I ask for some help for my problem? Since there are some further uses, I cannot change a new geometry without walls.
Thanks for your reading and help!

asked Jan 28, 2015 by Hrunx FEniCS Novice (910 points)

Could you give some more detail on the mathematical formulation or a minimal working code sample ?

...