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

SimpleFoam convergen problem

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

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   August 10, 2010, 10:38
Default SimpleFoam convergen problem
  #1
Member
 
Maolong LIU
Join Date: Apr 2010
Location: USA
Posts: 31
Rep Power: 15
maolongliu is on a distinguished road
Hi, Now I am using simpleFoam to calculate a pipe flow.
Now I face a converge problem.I am using low-Re k-epsilon model. After about 10 000 times irritation, k ,epsilon nearly converged, but ux and uy an not converge. Because the mesh number is very huge, and I have tried many ways, but still can not solve this problem. Mesh check is OK and the max skewness is about 0.5
Hope someone can give a help.
Below is the output:
================================================== ===============
smoothSolver: Solving for Ux, Initial residual = 0.00828136, Final residual = 7.04211e-05, No Iterations 9
smoothSolver: Solving for Uy, Initial residual = 0.0071747, Final residual = 4.91001e-05, No Iterations 10
smoothSolver: Solving for Uz, Initial residual = 7.61274e-06, Final residual = 8.03272e-07, No Iterations 4
GAMG: Solving for p_rgh, Initial residual = 0.000732179, Final residual = 7.03776e-05, No Iterations 13
GAMG: Solving for p_rgh, Initial residual = 0.000109843, Final residual = 1.07213e-05, No Iterations 27
time step continuity errors : sum local = 1.73118e-09, global = 8.26877e-13, cumulative = -5.21255e-10
smoothSolver: Solving for epsilon, Initial residual = 1.02623e-05, Final residual = 4.30998e-07, No Iterations 3
smoothSolver: Solving for k, Initial residual = 4.14421e-09, Final residual = 4.14421e-09, No Iterations 0
================================================== ===============
U:
Boundary condition:
inlet
{
type fixedValue;
value uniform (0 0 0.501);
}
outlet
{
type zeroGradient;
}
wall
{
type fixedValue;
value uniform (0 0 0);
}
}
p:
outlet
{
type fixedValue;
value uniform 0;
}
inlet
{
type zeroGradient;
}
wall
{
type zeroGradient;
}

fvSchemes:
divSchemes
{
default no;
div(rho*phi,U) Gauss linearUpwindV cellLimited Gauss linear 1;
div(phi,alpha1) Gauss linearUpwind Gauss linear;
div(phirb,alphai1) Gauss linearUpwind Gauss linear;
div(phi,k) Gauss linearUpwind Gauss linear;
div(phi,epsilon) Gauss linearUpwind Gauss linear;
div(phi,R) Gauss linear;
}

laplacianSchemes
{
default Gauss linear corrected;

}

interpolationSchemes
{
default linear;
}

snGradSchemes
{
default corrected;
}

fvSolution:
p_rgh
{
solver GAMG;
preconditioner DIC;
tolerance 1e-06;
relTol 0.1;
minIter 5;
maxIter 50;
smoother GaussSeidel;
nPreSweeps 1;
nPostSweeps 3;
nFinestSweeps 3;
scaleCorrection true;
directSolveCoarsest false;
cacheAgglomeration on;
nCellsInCoarsestLevel 50;
agglomerator faceAreaPair;
mergeLevels 1;
}
k
{
solver smoothSolver;
smoother GaussSeidel;
tolerance 1e-06;
relTol 0.01;
nSweeps 1;
maxIter 20;
}

epsilon
{
solver smoothSolver;
smoother GaussSeidel;
tolerance 1e-06;
relTol 0.01;
nSweeps 1;
maxIter 20;
}
U
{
solver smoothSolver;
smoother GaussSeidel;
tolerance 1e-06;
relTol 0.01;
nSweeps 1;
maxIter 20;
}
maolongliu is offline   Reply With Quote

Old   August 10, 2010, 11:46
Default
  #2
Member
 
Rasoul
Join Date: Feb 2010
Posts: 32
Rep Power: 16
aut_iut is on a distinguished road
Have you tried with a more fine mesh?
Try refineMesh command and repeat simulation again.
aut_iut is offline   Reply With Quote

Old   August 10, 2010, 11:56
Default
  #3
Member
 
Maolong LIU
Join Date: Apr 2010
Location: USA
Posts: 31
Rep Power: 15
maolongliu is on a distinguished road
Thanks, I'll try that later.
Quote:
Originally Posted by aut_iut View Post
Have you tried with a more fine mesh?
Try refineMesh command and repeat simulation again.
maolongliu is offline   Reply With Quote

Old   August 10, 2010, 15:15
Default
  #4
New Member
 
Ricardo Flatschart
Join Date: Apr 2010
Posts: 12
Rep Power: 15
rflats is on a distinguished road
You can try to decrease your relTol parameter, as it is probably preventing the linear solver to reach the specified tolerance.
The values you are using seems to be a bit high.
rflats is offline   Reply With Quote

Old   August 11, 2010, 02:47
Default
  #5
Member
 
Maolong LIU
Join Date: Apr 2010
Location: USA
Posts: 31
Rep Power: 15
maolongliu is on a distinguished road
I tried to decrease reTol parameter of p to 1e-3. the others to 0, but still, the residual of ux, uy and p decrease at first to about 1e-4, and then start to increase.

Quote:
Originally Posted by rflats View Post
You can try to decrease your relTol parameter, as it is probably preventing the linear solver to reach the specified tolerance.
The values you are using seems to be a bit high.
maolongliu is offline   Reply With Quote

Old   August 12, 2010, 03:27
Default
  #6
Member
 
Hrushikesh Khadamkar
Join Date: Jul 2010
Location: Mumbai
Posts: 68
Rep Power: 15
Hrushi is on a distinguished road
Hi Maolong LIU,

If you are having convergence problems with Above solver settings then try with PCG solver for p and PBiCG solver for other variables. It may take more time but it seems that these solvers are more stable as compared to GAMG solver.

You can also try Gauss linearUpwind cellLimited Gauss linear 1 scheme for divScheme for div (phi, U). for other variables limitedLinear 1.

Regards

Hrushikesh
Hrushi is offline   Reply With Quote

Old   August 12, 2010, 03:32
Default
  #7
Member
 
Maolong LIU
Join Date: Apr 2010
Location: USA
Posts: 31
Rep Power: 15
maolongliu is on a distinguished road
Thanks, Hrushikesh
Now I am trying to start with a higher viscosity, and then gradually decrease it to the real value. This seems work, the calculation is still going on. If this fail, I will try your advice.
Thanks again.

Quote:
Originally Posted by Hrushi View Post
Hi Maolong LIU,

If you are having convergence problems with Above solver settings then try with PCG solver for p and PBiCG solver for other variables. It may take more time but it seems that these solvers are more stable as compared to GAMG solver.

You can also try Gauss linearUpwind cellLimited Gauss linear 1 scheme for divScheme for div (phi, U). for other variables limitedLinear 1.

Regards

Hrushikesh
maolongliu is offline   Reply With Quote

Old   August 13, 2010, 11:17
Default
  #8
Member
 
Maolong LIU
Join Date: Apr 2010
Location: USA
Posts: 31
Rep Power: 15
maolongliu is on a distinguished road
I think I find the reason why my calculation seen hard to cenverge.
Now I am doing a 3D simulation, and through the calculation results I found tat ux and uy is very small, about 1e-5.
I heard some said that the convergence criteria of OpenFOAM for this kond very small results may have some problem.
Is this true?
Thanks.

Quote:
Originally Posted by maolongliu View Post
Thanks, Hrushikesh
Now I am trying to start with a higher viscosity, and then gradually decrease it to the real value. This seems work, the calculation is still going on. If this fail, I will try your advice.
Thanks again.
maolongliu is offline   Reply With Quote

Reply

Thread Tools Search this Thread
Search this Thread:

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


Similar Threads
Thread Thread Starter Forum Replies Last Post
UDF compiling problem Wouter Fluent UDF and Scheme Programming 6 June 6, 2012 05:43
Incoherent problem table in hollow-fiber spinning Gianni FLUENT 0 April 5, 2008 11:33
natural convection problem for a CHT problem Se-Hee CFX 2 June 10, 2007 07:29
Adiabatic and Rotating wall (Convection problem) ParodDav CFX 5 April 29, 2007 20:13
convergence problem Trushar Phoenics 5 August 28, 2002 00:40


All times are GMT -4. The time now is 08:51.