CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > OpenFOAM Running, Solving & CFD

continuity error with simpleFOAM in duct flow

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

Reply
 
LinkBack Thread Tools Display Modes
Old   February 11, 2010, 09:56
Default continuity error with simpleFOAM in duct flow
  #1
New Member
 
Luc Bordier
Join Date: Feb 2010
Posts: 10
Rep Power: 7
lbordier is on a distinguished road
Hello,

I'm trying to set up a case for duct flow with simpleFOAM.
I know the inlet & outlet pressure, so I tried to fix BC in consequence but i guess i've mistaken somewhere because I've very high continuity errors.

Thank you for any comment / advice.

The setting of the bc and some logs are given below :

for velocity :
dimensions [0 1 -1 0 0 0 0];

internalField uniform (0 0 0);

boundaryField
{
Inlet
{
type pressureInletVelocity;
value uniform (0 0 0);
}

Outlet
{
type zeroGradient;
//value uniform (0 0 0);
}

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


for pressure :

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

internalField uniform 0;

boundaryField
{
Inlet
{
type fixedValue;
value uniform 0;
}

Outlet
{
type fixedValue;
value uniform -1659.751;
}

Wall
{
type zeroGradient;
}
}


some log for laminar run :

Time = 268

DILUPBiCG: Solving for Ux, Initial residual = 0.0250111, Final residual = 1.79843e-06, No Iterations 3
DILUPBiCG: Solving for Uy, Initial residual = 0.0247072, Final residual = 1.19593e-06, No Iterations 3
DILUPBiCG: Solving for Uz, Initial residual = 0.0285379, Final residual = 1.2351e-06, No Iterations 3
GAMG: Solving for p, Initial residual = 0.0515492, Final residual = 0.000424851, No Iterations 5
GAMG: Solving for p, Initial residual = 0.0101859, Final residual = 9.71891e-05, No Iterations 3
time step continuity errors : sum local = 259841, global = -10979.7, cumulative = -206963
ExecutionTime = 7931.91 s ClockTime = 7933 s

Time = 269

DILUPBiCG: Solving for Ux, Initial residual = 0.0250091, Final residual = 1.86461e-06, No Iterations 3
DILUPBiCG: Solving for Uy, Initial residual = 0.0247248, Final residual = 1.15719e-06, No Iterations 3
DILUPBiCG: Solving for Uz, Initial residual = 0.0285198, Final residual = 1.22932e-06, No Iterations 3
GAMG: Solving for p, Initial residual = 0.0516174, Final residual = 0.000425292, No Iterations 5
GAMG: Solving for p, Initial residual = 0.0101802, Final residual = 9.74468e-05, No Iterations 3
time step continuity errors : sum local = 272986, global = -11533.7, cumulative = -218496
ExecutionTime = 7962.34 s ClockTime = 7964 s

Mesh check :

Checking geometry...
Domain bounding box: (0.478596 0.0186977 0.111428) (0.683153 0.121314 0.475742)
Boundary openness (-3.99958e-20 -1.02563e-19 -8.1168e-20) OK.
Max cell openness = 2.13144e-16 OK.
Max aspect ratio = 6.79448 OK.
Minumum face area = 1.45022e-09. Maximum face area = 6.95171e-06. Face area magnitudes OK.
Min volume = 2.10867e-13. Max volume = 2.42119e-08. Total volume = 0.00117081. Cell volumes OK.
Mesh non-orthogonality Max: 71.8538 average: 13.8243
*Number of severely non-orthogonal faces: 1.
Non-orthogonality check OK.
<<Writing 1 non-orthogonal faces to set nonOrthoFaces
Face pyramids OK.
Max skewness = 2.20433 OK.
*There are 250 faces with concave angles between consecutive edges. Max concave angle = 89.457 degrees.
<<Writing 250 faces with concave angles to set concaveFaces
Face flatness (1 = flat, 0 = butterfly) : average = 0.990544 min = 0.755141
*There are 47 faces with ratio between projected and actual area < 0.8
Minimum ratio (minimum flatness, maximum warpage) = 0.755141
<<Writing 47 warped faces to set warpedFaces

Mesh OK.
lbordier is offline   Reply With Quote

Old   February 12, 2010, 06:10
Default
  #2
Senior Member
 
santos's Avatar
 
Jose Luis Santos
Join Date: Mar 2009
Location: Portugal
Posts: 213
Rep Power: 9
santos is on a distinguished road
Send a message via Skype™ to santos
Hi,

Please change the pressure inlet boundary condition to zeroGradient.

Regards,
Jose Santos
santos is offline   Reply With Quote

Old   February 15, 2010, 06:20
Default
  #3
New Member
 
Luc Bordier
Join Date: Feb 2010
Posts: 10
Rep Power: 7
lbordier is on a distinguished road
It seems that it converge better with the following settings for p :

boundaryField
{
Inlet
{
//type zeroGradient;
type fixedValue;
value uniform 0.0;
}

Outlet
{
//type zeroGradient;
type timeVaryingUniformFixedValue;
fileName "pressure_ramp.dat";
outOfBounds clamp;
//value uniform -2000.0;
}

Wall
{
type zeroGradient;
}
}

adding a pressure ramp makes the continuity errors to be reduced
lbordier is offline   Reply With Quote

Reply

Thread Tools
Display Modes

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 Off
Trackbacks are On
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
Turbulent Flow in a Square Duct using LES Hock Ming FLUENT 0 February 7, 2009 21:25
Convergence moving mesh lr103476 OpenFOAM Running, Solving & CFD 30 November 19, 2007 15:09
IcoFoam parallel woes msrinath80 OpenFOAM Running, Solving & CFD 9 July 22, 2007 02:58
Turbulence flow in a rectangular duct Watchapon Main CFD Forum 0 April 7, 2007 07:34
Straight duct flow in CFX 5.2 Bo Jensen Main CFD Forum 3 November 25, 1998 07:15


All times are GMT -4. The time now is 16:58.