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

SimpleFoam unable to reach steady state

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

Reply
 
LinkBack Thread Tools Display Modes
Old   January 28, 2009, 06:59
Default Hi Foamers, I'm playing wit
  #1
Senior Member
 
Francois Beaubert
Join Date: Mar 2009
Location: Lille, France
Posts: 138
Rep Power: 9
francois is on a distinguished road
Hi Foamers,

I'm playing with the simpleFoam solver for a simple open rotor/stator laminar flow configuration but I'm not able to reach the steady state of the flow.

The angular velocity on the rotor is around 300 tr/min so a Reynolds number (based on the rotor radius) around 2.10^5 is expected. According to the literature the flow is still laminar.

So I try simpleFoam with no turbulence models but even if convergence seems ok I'm not able to reach this steady state. Perhaps the flow is locally turbulent or in transition and a RANS model could be the way to go. Some experiences on those kind of flows ?

Here is one of my last iteration output in the log file:

Time = 95000

DILUPBiCG: Solving for Ux, Initial residual = 0.0114524, Final residual = 0.000374143, No Iterations 2
DILUPBiCG: Solving for Uy, Initial residual = 0.00827667, Final residual = 0.000238708, No Iterations 2
DILUPBiCG: Solving for Uz, Initial residual = 0.00574431, Final residual = 0.000142418, No Iterations 2
DICPCG: Solving for p, Initial residual = 0.0246974, Final residual = 0.000246175, No Iterations 54
time step continuity errors : sum local = 7.81603e-05, global = -2.14943e-06, cumulative = 0.07646
ExecutionTime = 182884 s ClockTime = 183024 s

Any advices to get better convergence to reach the steady state of the flow ?

My setup is as follow:

===========
:fvSchemes:
===========

ddtSchemes
{
default steadyState;
}

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

divSchemes
{
default none;
div(phi,U) Gauss upwind;
div(phi,k) Gauss upwind;
div(phi,epsilon) Gauss upwind;
div(phi,R) Gauss upwind;
div(R) Gauss linear;
div(phi,nuTilda) Gauss upwind;
div((nuEff*dev(grad(U).T()))) 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(DepsilonEff,epsilon) Gauss linear corrected;
laplacian(DREff,R) Gauss linear corrected;
laplacian(DnuTildaEff,nuTilda) Gauss linear corrected;
}

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

snGradSchemes
{
default corrected;
}

fluxRequired
{
default no;
p;
}

============
:fvSolution:
============

solvers
{
p PCG
{
preconditioner DIC;
tolerance 1e-06;
relTol 0.01;
};
U PBiCG
{
preconditioner DILU;
tolerance 1e-05;
relTol 0.1;
};
k PBiCG
{
preconditioner DILU;
tolerance 1e-05;
relTol 0.1;
};
epsilon PBiCG
{
preconditioner DILU;
tolerance 1e-05;
relTol 0.1;
};
R PBiCG
{
preconditioner DILU;
tolerance 1e-05;
relTol 0.1;
};
nuTilda PBiCG
{
preconditioner DILU;
tolerance 1e-05;
relTol 0.1;
};
}

SIMPLE
{
nNonOrthogonalCorrectors 0;
}

relaxationFactors
{
p 0.3;
U 0.7;
k 0.7;
epsilon 0.7;
R 0.7;
nuTilda 0.7;
}


===============
:RASProperties:
===============

RASModel laminar;

turbulence off;

printCoeffs on;

laminarCoeffs
{
}


Thank you very much.
Have a nice day

Francois
francois is offline   Reply With Quote

Old   February 25, 2009, 12:33
Default Hi all, Hopefully I have so
  #2
Senior Member
 
Francois Beaubert
Join Date: Mar 2009
Location: Lille, France
Posts: 138
Rep Power: 9
francois is on a distinguished road
Hi all,

Hopefully I have some time to go back on this case ...

My runs don't reach the steady state and convergence is still not achieved with residuals not so low.

I've attached some pictures to better describe my case. It's a simple rotor/stator configuration with a hole in the stator.

rotor_stator.jpg

I've also attached the distribution of the u,v,w velocity components on a radial line between the rotor and stator.

velocity_distribution.jpg

You can see that steady state is not reach here: one can observe large variation of u,v,w profiles between iteration 80000 and 90000. The behaviour is the same for a greater number of iteration (currently I'm at 120000 iterations).

Below you can see isocontours of the u component around the rotor/stator and one can note a "shear layer" near the rotor which oscillates with iterations (this trend is also observed on the uvw distribution figure just above)

rotor_stator_Ux.jpg

And to conclude below are the evolution of residuals with the iterations for Ux and Uy for example:

UxFinalRes.jpg

UyFinalRes.jpg

So I still have a few questions:

* what is the best procedure to obtain a steady state

* is my computational domain too small (Ddomain = 4rotor and Hdomain = 3Drotor)for this problem with a "feedback" of the outlet boundary condition
?

I will be very pleased to have some advices/ideas frome somone on those points and forgive me for this long post.

Thanks you very much
Regards

Francois

Last edited by francois; March 18, 2009 at 16:56. Reason: Add images
francois is offline   Reply With Quote

Old   February 25, 2009, 12:40
Default Oup's sorry for the images Wi
  #3
Senior Member
 
Francois Beaubert
Join Date: Mar 2009
Location: Lille, France
Posts: 138
Rep Power: 9
francois is on a distinguished road
Oup's sorry for the images
Will try again here :

Computationaldomain:



distribution of the u,v,w velocity components on a radial line between the rotor and stator:



u v w contours:





Residuals:










francois is offline   Reply With Quote

Old   February 25, 2009, 12:50
Default Humm ... sorry Don't work her
  #4
Senior Member
 
Francois Beaubert
Join Date: Mar 2009
Location: Lille, France
Posts: 138
Rep Power: 9
francois is on a distinguished road
Humm ... sorry
Don't work here with firefox 3.06 on ubuntu 8.10

What the best way to send/attach picture or files on the message board ?
I use:

\ image{foo} or \ attach{foo} with no space between \ and the after word

It seems that I always have to struggle until death to be able to send some files or images here.
francois is offline   Reply With Quote

Old   February 25, 2009, 23:46
Default I think your syntax is probabl
  #5
Member
 
Andrew King
Join Date: Mar 2009
Location: Perth, Western Australia, Australia
Posts: 82
Rep Power: 9
andersking is on a distinguished road
I think your syntax is probably fine, but there's a very small limit on image sizes, 50K I think. It doesn't stop you from posting, but won't show the images.

Less than 50K is fine.

Cheers
Andrew
__________________
Dr Andrew King
Fluid Dynamics Research Group
Curtin University
andersking is offline   Reply With Quote

Old   February 26, 2009, 06:18
Default Thanks Andrew for your advice
  #6
Senior Member
 
Francois Beaubert
Join Date: Mar 2009
Location: Lille, France
Posts: 138
Rep Power: 9
francois is on a distinguished road
Thanks Andrew for your advice

But my images sizes are always under 50k and I never saw the dialog for file upload

So the problem seems to be elsewhere ...

Cheers
Francois
francois is offline   Reply With Quote

Old   February 27, 2009, 01:19
Default Hmmm, hopefully it will all ge
  #7
Member
 
Andrew King
Join Date: Mar 2009
Location: Perth, Western Australia, Australia
Posts: 82
Rep Power: 9
andersking is on a distinguished road
Hmmm, hopefully it will all get betted with the new forum software.

Cheers
Andrew
__________________
Dr Andrew King
Fluid Dynamics Research Group
Curtin University
andersking is offline   Reply With Quote

Old   March 18, 2009, 16:59
Default
  #8
Senior Member
 
Francois Beaubert
Join Date: Mar 2009
Location: Lille, France
Posts: 138
Rep Power: 9
francois is on a distinguished road
I've just modified the first post and was able to include some of the images which could not be posted in the last forum...

Yeah ... the new forum is better
francois is offline   Reply With Quote

Old   November 19, 2009, 14:33
Default Same problem too...
  #9
New Member
 
Valerio Novaresio
Join Date: Mar 2009
Location: Polonghera, Cuneo, Italy
Posts: 27
Rep Power: 9
novyno is on a distinguished road
I have the same convergence problem... Can someone help us?
__________________
...The best way to acquire new knowledge is to share it...
novyno 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
Convergence problem using simpleFoam steady state vvqf OpenFOAM Running, Solving & CFD 12 May 18, 2011 07:51
time to reach steady state with transient simulati prabhu CFX 2 April 21, 2008 10:18
What if I use turbfoam instead of simplefoam what is steady turblence lakeat OpenFOAM Running, Solving & CFD 4 January 18, 2008 13:13
How to reach steady state solution Sally FLUENT 10 June 14, 2007 09:26
Time to reach steady state for cavity flow? Quarkz Main CFD Forum 8 July 15, 2005 21:42


All times are GMT -4. The time now is 19:42.