# Problem with pressure boundary conditions

 Register Blogs Members List Search Today's Posts Mark Forums Read

June 20, 2013, 21:31
Problem with pressure boundary conditions
#1
New Member

Saša G.
Join Date: May 2013
Location: Zagreb
Posts: 1
Rep Power: 0
Greatings to everyone.
I hope someone can help me with a problem.
I have been learning to use openFoam for past couple of weeks so I could use it for a bit of analyzing of Smagorinsky model, which is topic of my thesis.
I have run already a couple simple tests like lid driven cavity flow, flow over an obstacle and so on, mainly simple 2D test problems.
My mentor sugested i try to model a blood flow through bifurcation with an aneurysm.
So i have basicly three pipes, of different diameters, one inlet, two outlets, that are connected at a junction and one of the pipes has a widening. The fluid is assumed to be incompressible Newtonian fluid.
At first try I have setup the case in a way that i prescribe time-varying Pouiselle velocity profile on inlet, fixedValue pressure and zeroGradient velocity on outlet.
So far so good, I solved the case with IcoFoam and with pimpleFoam with Smagorinsky model for comparisson.
Then my mentor suggested I try to prescribe time-varying pressure on inlet and fixedValue pressure on outlets. Thats where the trouble begins.
Again i try to solve it with icoFoam first. I used the folowing initial-B.C.:
Code:
```
0/p
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

dimensions      [0 2 -2 0 0 0 0];

internalField   uniform 0;

boundaryField
{

inlet
{
type            uniformFixedValue;
uniformValue    tableFile;
tableFileCoeffs
{
fileName    "pressuredrop.txt";
outOfBounds    repeat;
}
}

outlet1
{
type            fixedValue;
value           uniform 0;

}

outlet2
{
type            fixedValue;
value           uniform 0;

}
walls
{
}

UpDown
{
type            empty;
}
}

// ************************************************************************* //```
and following 0/U

Code:
```
0/U
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

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

internalField   uniform (0 0 0);

boundaryField
{
inlet
{
type            pressureInletOutletVelocity;
value           uniform (0 0 0);
}

outlet1
{
}

outlet2
{

}

walls
{
type    fixedValue;
value    uniform (0 0 0);
}

UpDown
{
type    empty;
}
}

// ************************************************************************* //```
No matter how small deltaT i use, sollution blows up at aproxx t = 0.206. When i look at the results just before blowup i can see that velocity tends to grow to physicaly impossible proportions, and it happens mainly in a small part of the inlet. Now I'm not sure is this due to icoFoam-s limitations, or due to bad mesh, or maybe bad choice of boundary conditions.
When i change inlet in 0/U to:
Code:
```inlet
{
type            pressureInletUniformVelocity;
value           uniform (0 0 0);
}```
the simulaton runs it's course whitout a problem because, that inlet condition "irons" the velocity so it never get the chance to blow up anywhere. Unfortunatly that constraint is unphysical.

I would think it's icoFoam thats causing the trouble but i'm not sure because very similar case with pulsating Pouiselle profile works like a charm.
I would think that it's the mesh but mesh looks ok and pasess checkMesh with reasonable results, screenshot is attached below
So i'm thinking it's boundary conditions, since the whole trouble starts at inlet.

I have tried using different B.C for velocity on outlet too, ( InletOutlet ) but to no avail. I have also tried to use different, nCorrection, nNonOrthCorr and so on but nothing helped.

Sorry if it's a stupid question, I'm very new to OF

Attached Images
 mesh.jpg (24.9 KB, 26 views)

 Thread Tools Display Modes Linear Mode

 Posting Rules You may not post new threads You may not post replies You may not post attachments You may not edit your posts BB code is On Smilies are On [IMG] code is On HTML code is OffTrackbacks are On Pingbacks are On Refbacks are On Forum Rules

 Similar Threads Thread Thread Starter Forum Replies Last Post volo87 CFX 5 June 14, 2013 17:44 Attesz CFX 7 January 5, 2013 04:32 mbgnfrx2 STAR-CD 2 November 26, 2012 00:55 M_Flodin FLUENT 0 April 1, 2011 09:32 lyang Main CFD Forum 0 September 19, 1999 18:29