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

Bad pressure residuals

Register Blogs Community New Posts Updated Threads Search

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   December 10, 2015, 09:14
Default Bad pressure residuals
  #1
New Member
 
A Z
Join Date: Dec 2015
Posts: 16
Rep Power: 10
don6alfonso is on a distinguished road
Hey,
I´m running a simpleFoam simulation. It´s a inducer. I use OF 2.3.1.
The mesh is splinted in 3 cell zones. Is just simulate the half and it is copied over cyclic.

here my output:

Time = 1998

smoothSolver: Solving for Ux, Initial residual = 0.000486577, Final residual = 3.79976e-05, No Iterations 2
smoothSolver: Solving for Uy, Initial residual = 0.000456892, Final residual = 3.53004e-05, No Iterations 2
smoothSolver: Solving for Uz, Initial residual = 0.000887549, Final residual = 7.22756e-05, No Iterations 2
GAMG: Solving for p, Initial residual = 0.0220465, Final residual = 0.0009725, No Iterations 5
time step continuity errors : sum local = 0.0180951, global = 0.000183385, cumulative = -0.970663
smoothSolver: Solving for omega, Initial residual = 1.58397e-06, Final residual = 7.70641e-08, No Iterations 2
bounding omega, min: -1687.16 max: 3.07084e+09 average: 376133
smoothSolver: Solving for k, Initial residual = 0.00104681, Final residual = 0.000103889, No Iterations 2
bounding k, min: -0.00542659 max: 0.977381 average: 0.0406476
ExecutionTime = 4503.4 s ClockTime = 4527 s

Time = 1999

smoothSolver: Solving for Ux, Initial residual = 0.000484606, Final residual = 3.78135e-05, No Iterations 2
smoothSolver: Solving for Uy, Initial residual = 0.000455814, Final residual = 3.52004e-05, No Iterations 2
smoothSolver: Solving for Uz, Initial residual = 0.000883571, Final residual = 7.17319e-05, No Iterations 2
GAMG: Solving for p, Initial residual = 0.0219943, Final residual = 0.00109084, No Iterations 4
time step continuity errors : sum local = 0.0203253, global = -0.000235497, cumulative = -0.970898
smoothSolver: Solving for omega, Initial residual = 1.60132e-06, Final residual = 7.71437e-08, No Iterations 2
smoothSolver: Solving for k, Initial residual = 0.00104701, Final residual = 0.000104054, No Iterations 2
bounding k, min: -9.39376e-07 max: 0.977339 average: 0.0406466
ExecutionTime = 4505.5 s ClockTime = 4529 s

Time = 2000

smoothSolver: Solving for Ux, Initial residual = 0.000483367, Final residual = 3.76716e-05, No Iterations 2
smoothSolver: Solving for Uy, Initial residual = 0.000454811, Final residual = 3.50667e-05, No Iterations 2
smoothSolver: Solving for Uz, Initial residual = 0.000885823, Final residual = 7.24752e-05, No Iterations 2
GAMG: Solving for p, Initial residual = 0.0219506, Final residual = 0.00097752, No Iterations 5
time step continuity errors : sum local = 0.0182286, global = 0.000121839, cumulative = -0.970776
smoothSolver: Solving for omega, Initial residual = 1.57142e-06, Final residual = 7.6524e-08, No Iterations 2
smoothSolver: Solving for k, Initial residual = 0.00104601, Final residual = 0.000103989, No Iterations 2
bounding k, min: -0.00553487 max: 0.977304 average: 0.040646
ExecutionTime = 4510.08 s ClockTime = 4534 s


Settings:

fvSchemes:

ddtSchemes
{
default steadyState;
}

gradSchemes
{
default Gauss linear;
grad(p) Gauss linear;
grad(U) Gauss linear;
}

divSchemes
{
default none;
div(phi,U) bounded Gauss limitedLinearV 1;
div(phi,k) bounded Gauss limitedLinear 1;
div(phi,omega) bounded Gauss limitedLinear 1;
div((nuEff*dev(T(grad(U))))) Gauss linear;
}

laplacianSchemes
{
default none;
laplacian(nuEff,U) Gauss linear corrected;
laplacian((1|A(U)),p) Gauss linear corrected;
laplacian(DkEff,k) Gauss linear corrected;
laplacian(DomegaEff,omega) Gauss linear corrected;
}

interpolationSchemes
{
default linear;
interpolate(U) linear;
}

snGradSchemes
{
default corrected;
}

fluxRequired
{
default no;
p ;
}



fvSolution:

p
{
solver GAMG;
tolerance 1e-08;
relTol 0.05;
smoother GaussSeidel;
cacheAgglomeration true;
nCellsInCoarsestLevel 20;
agglomerator faceAreaPair;
mergeLevels 1;
}
...
SIMPLE
{
nNonOrthogonalCorrectors 0;
pRefCell 0;
pRefValue 0;
}

relaxationFactors
{
fields
{
p 0.5;
}

equations
{
U 0.6;
k 0.6;
omega 0.5;
}


p:

internalField uniform 0;

boundaryField
{
INFLOW-IN
{
type zeroGradient;
}

OUTFLOW-OUT
{
type fixedValue;
value $internalField;
}

INDUCER-BLADE
{
type zeroGradient;
}

"(INFLOW-HUB|...|INDUCER-SHROUD)" //wall
{
type zeroGradient;
}

"(INDUCER-OUTFLOW-SIDE-1|...|INFLOW-INDUCER-SIDE-2)" //GGI
{
type cyclicAMI;
value $internalField;
}

"(O-P-S-1|O-P-S-2|I-P-S-1|I-P-S-2|If-P-S-1|If-P-S-2)" //periodic
{
type cyclic;
}

}


I hope someone can help me.

Best Regards
don6alfonso is offline   Reply With Quote

Reply

Tags
pressure, residuals


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


Similar Threads
Thread Thread Starter Forum Replies Last Post
Pressure Outlet Guage pressure Mohsin FLUENT 36 April 29, 2016 17:16
Pulsatile pressure inlet with pressure outlet a.lynchy FLUENT 3 March 23, 2012 13:45
UDF to define or adjust pressure??? engahmed FLUENT 0 July 6, 2010 17:19
Problems of Duns Codes! Martin J Main CFD Forum 8 August 14, 2003 23:19
pressure gradient term in low speed flow Atit Koonsrisuk Main CFD Forum 2 January 10, 2002 10:52


All times are GMT -4. The time now is 07:28.