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

convergence problem

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

Reply
 
LinkBack Thread Tools Display Modes
Old   June 17, 2011, 04:08
Post convergence problem
  #1
Member
 
Gregoire Junqua
Join Date: Jun 2011
Location: China
Posts: 58
Blog Entries: 1
Rep Power: 7
gregjunqua is on a distinguished road
Hi there.
I am quite new into modelling but i got a convergence problem as follow DICPCG ask a many iteration. How to make it less?
is improving the grid resolution or lowering the grid resolution could help?
or a better approximation of k and epsilon is the solution?

And all link about it are welcome


Quote:
Courant Number mean: 0.0750431 max: 0.291312
DILUPBiCG: Solving for Ux, Initial residual = 0.00135266, Final residual = 1.80396e-06, No Iterations 1
DILUPBiCG: Solving for Uy, Initial residual = 0.00114192, Final residual = 1.76496e-06, No Iterations 1
DICPCG: Solving for p, Initial residual = 0.047671, Final residual = 0.00470314, No Iterations 17
time step continuity errors : sum local = 7.64284e-08, global = -2.95465e-10, cumulative = 1.88697e-06
DICPCG: Solving for p, Initial residual = 0.0290166, Final residual = 9.92062e-07, No Iterations 777
time step continuity errors : sum local = 1.58005e-11, global = -3.76457e-14, cumulative = 1.88697e-06
DILUPBiCG: Solving for epsilon, Initial residual = 0.00315981, Final residual = 6.3086e-06, No Iterations 1
DILUPBiCG: Solving for k, Initial residual = 0.0141768, Final residual = 7.62168e-06, No Iterations 1
ExecutionTime = 655.38 s ClockTime = 658 s
gregjunqua is offline   Reply With Quote

Old   June 17, 2011, 04:34
Post hello...
  #2
Senior Member
 
Tushar Chourushi
Join Date: Jul 2009
Location: IIT-Indore, India
Posts: 319
Blog Entries: 1
Rep Power: 9
Tushar@cfd is on a distinguished road
I think the problem is because of the mesh resolution. I just want to know which type of grid have you used? Is it Unstructured?
If it's unstructured then it's usual you will observe the error, check out with structured if you get the same.


--------

Last edited by Tushar@cfd; June 17, 2011 at 06:27.
Tushar@cfd is offline   Reply With Quote

Old   June 18, 2011, 00:02
Default
  #3
Member
 
Alan Russell
Join Date: Aug 2009
Location: Boise, Idaho USA
Posts: 61
Rep Power: 8
AlanR is on a distinguished road
Check the mesh quality with checkMesh - you could have problems with skewed cells, etc. If you get good results from checkMesh, you are ready to start working on the solution. First, instead of refining the mesh, I would do the opposite. A coarse mesh lets you get the 'outline' of the solution quickly, and shows if your setup is valid. If you can get a coarse mesh solution to converge, then start refining your mesh.

Good luck,
Alan
AlanR is offline   Reply With Quote

Old   June 20, 2011, 21:14
Default
  #4
Member
 
Gregoire Junqua
Join Date: Jun 2011
Location: China
Posts: 58
Blog Entries: 1
Rep Power: 7
gregjunqua is on a distinguished road
Hi thanks for your answers.

What do you mean by unstructured?
My grid got a sinus form...
and my
my grid resolution is around 0.01m
and dT = 0.001;

This is the answer of the checkMesh
Quote:

Checking topology...
Boundary definition OK.
Point usage OK.
Upper triangular ordering OK.
Face vertices OK.
Number of regions: 1 (OK).

Checking patch topology for multiply connected surfaces ...
Patch Faces Points Surface topology
movingWall 1500 3002 ok (non-closed singly connected)
fixedWalls 1100 2202 ok (non-closed singly connected)
outlet 400 802 ok (non-closed singly connected)
frontAndBack 880000 883002 ok (non-closed singly connected)

Checking geometry...
Overall domain bounding box (0 -0.5 0) (11 10 0.1)
Mesh (non-empty, non-wedge) directions (1 1 0)
Mesh (non-empty) directions (1 1 0)
All edges aligned with or perpendicular to non-empty directions.
Boundary openness (-6.5805e-20 5.32246e-20 -7.18612e-17) OK.
Max cell openness = 1.78614e-16 OK.
Max aspect ratio = 4.0807 OK.
Minumum face area = 0.000237529. Maximum face area = 0.002625. Face area magnitudes OK.
Min volume = 2.37529e-05. Max volume = 2.62471e-05. Total volume = 10.9684. Cell volumes OK.
Mesh non-orthogonality Max: 57.1056 average: 28.6341
Non-orthogonality check OK.
Face pyramids OK.
Max skewness = 1.29563 OK.

Mesh OK.
Today i will do some test by improving the grid and the time resolution hope it would work
gregjunqua is offline   Reply With Quote

Old   June 21, 2011, 00:20
Default
  #5
Member
 
Alan Russell
Join Date: Aug 2009
Location: Boise, Idaho USA
Posts: 61
Rep Power: 8
AlanR is on a distinguished road
In your checkMesh, I see an outlet but no inlet, which will cause problems with the continuity equation (conservation of mass). Most of your setup is like the cavity tutorial, except for the outlet. Can you describe what you are trying to model? You will certainly need some kind of inlet if you have an outlet.

Alan
AlanR is offline   Reply With Quote

Old   June 21, 2011, 01:37
Default
  #6
Senior Member
 
Tushar Chourushi
Join Date: Jul 2009
Location: IIT-Indore, India
Posts: 319
Blog Entries: 1
Rep Power: 9
Tushar@cfd is on a distinguished road
Can you please upload your case file??
Tushar@cfd is offline   Reply With Quote

Old   June 21, 2011, 02:55
Default
  #7
Member
 
Gregoire Junqua
Join Date: Jun 2011
Location: China
Posts: 58
Blog Entries: 1
Rep Power: 7
gregjunqua is on a distinguished road
I don't specify the inlet, i just want a uniform wind field at the boundaries


Moving wall
+-------------------------------------------------------------------+
| .................................................. ........................................ |
|................................................. .......................................... |
Moving Wall .................................................. .....................outlet
| .................................................. .........................................|
| .................................................. .........................................|
_.-._.-._.-._.-._.-._.-._.-._.-._-._.-._.-._.-._.-._.-._.-._.-._.|

The ground is a sin function, and i want to view the turbulence around it.

[QUOTEFoamFile
{
version 2.0;
format ascii;
class volVectorField;
object U;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

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

internalField uniform (0 0 0);

boundaryField
{
movingWall
{
type fixedValue;
value uniform (1 0 0);
}

fixedWalls
{
type fixedValue;
value uniform (0 0 0);
}
outlet
{
type zeroGradient;
}
frontAndBack
{
type empty;
}
}

// ************************************************** *********************** //
][/QUOTE]

If you need more just ask
Thank you for your interest anyway
gregjunqua is offline   Reply With Quote

Old   June 21, 2011, 03:33
Default
  #8
Member
 
Gregoire Junqua
Join Date: Jun 2011
Location: China
Posts: 58
Blog Entries: 1
Rep Power: 7
gregjunqua is on a distinguished road
I had try a lower grid scale as around 0.05m and there is a good convergence

...
How to explain it? my first guess was more the grid scale is small more the convergence is fast ...
gregjunqua is offline   Reply With Quote

Reply

Tags
convergence, rans

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
convergence problem when use pisoFoam, LES for wind tunnel case Forrest_Lei OpenFOAM 3 July 19, 2011 06:00
convergence problem commonyue Main CFD Forum 1 December 1, 2009 04:54
Convergence of CFX field in FSI analysis nasdak CFX 2 June 29, 2009 01:17
3D Fluid Flow Convergence problem Emily FLUENT 2 March 21, 2007 23:18
Non Convergence of 3D Heat transfer cfd problem Balraj Main CFD Forum 3 December 9, 2004 01:24


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