[liguifan]

Dear all,

I want to simulate an artery with bifurcation. After using Gmsh to generate a mesh file and convert it to openfoam.

You can see the screenshot in attachment, there are one "inlet", one "straight_outlet", one "bent_outlet". And there are some defaultFaces at the interaction.

When I using icoFoam to run the file, the program will be interupted because courant number become to large.

I am just wondering it that because the defaultFaces at interaction block the fluids entering the bifurcation(daughter arteries)?

How can I get rid of the defaultFaces, I mean delete them( I can not delete it in gmsh, and delete them in polyMesh/boundary as well)?


Here are my p and U file.

U
dimensions [0 1 -1 0 0 0 0];

internalField uniform (0 0 0);

boundaryField
{
artery_wall
{
type fixedValue;
value uniform (0 0 0);
}
inlet
{
type fixedValue;
value uniform (0 0 -0.11);
}
straight_outlet
{
type zeroGradient;
}
bent_outlet
{
type zeroGradient;
}

}

p
dimensions [0 2 -2 0 0 0 0];
internalField uniform 0;
boundaryField
{
artery_wall
{
type zeroGradient;
}
inlet
{
type zeroGradient;
}
straight_outlet
{
type fixedValue;
value uniform 0;
}
bent_outlet
{
type fixedValue;
value uniform 0;
}
}


Many thanks!