CFD Online URL
[Sponsors]
Home > Forums > OpenFOAM Pre-Processing

Convergence crtieria

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

Reply
 
LinkBack Thread Tools Display Modes
Old   April 3, 2007, 13:31
Default Hi I am trying to test SIMP
  #1
lam
Member
 
Hoang-Lam
Join Date: Mar 2009
Location: Lausanne, Switzerland
Posts: 60
Rep Power: 7
lam is on a distinguished road
Hi

I am trying to test SIMPLEFOAM with a Bird-Carreau law and I have a very trivial question: I noticed that, when running the non linear computation, I have to specify the number of iterations a priori, (their number is being suggested by the proogramme itself).

Is it possible to run the programme in such a way that the computation interrupts automatically as soon as convergence is reached ? (comvergence criterium).

Thank you in davance,

Lam
lam is offline   Reply With Quote

Old   April 4, 2007, 06:07
Default Oh, THANKS a lot Mark, for rep
  #2
lam
Member
 
Hoang-Lam
Join Date: Mar 2009
Location: Lausanne, Switzerland
Posts: 60
Rep Power: 7
lam is on a distinguished road
Oh, THANKS a lot Mark, for replying me!

I am going give a go right now!

Cheers,

Lam
lam is offline   Reply With Quote

Old   April 5, 2007, 05:40
Default I guess I wasn't explicit enou
  #3
Senior Member
 
Mark Olesen
Join Date: Mar 2009
Location: http://olesenm.github.io/
Posts: 777
Rep Power: 17
olesen will become famous soon enough
I guess I wasn't explicit enough. You can either solve the equation and ignore the return value:

solve(UEqn() == -fvc::grad(p));

OR solve the equation and use the return value:

eqnResidual = solve(UEqn() == -fvc::grad(p)).initialResidual();

but not both.

It is also essential to capture the mass residual. ie. replace

pEqn.solve();

with
if (nonOrth == 0)
{
eqnResidual = pEqn().solve().initialResidual();
maxResidual = max(eqnResidual, maxResidual);
}
else
{
pEqn().solve();
}

The residuals from the turbulence equation are not directly accessible, but this might change in OpenFOAM 1.4
olesen is offline   Reply With Quote

Old   April 5, 2007, 07:46
Default It looks mostly reasonable exc
  #4
Senior Member
 
Mark Olesen
Join Date: Mar 2009
Location: http://olesenm.github.io/
Posts: 777
Rep Power: 17
olesen will become famous soon enough
It looks mostly reasonable except you don't need this line twice.

eqnResidual = solve(UEqn() == -fvc::grad(p)).initialResidual();
olesen is offline   Reply With Quote

Old   April 5, 2007, 10:20
Default Hi again Mark, I had remove
  #5
lam
Member
 
Hoang-Lam
Join Date: Mar 2009
Location: Lausanne, Switzerland
Posts: 60
Rep Power: 7
lam is on a distinguished road
Hi again Mark,

I had removed the line which is not needed as you said. And of course, I did the Allmake to update the simpleFoam.C file.
I made a computation, but the result has not changed.
You can take a look at my fvSolution file:

***********************

solvers
{

p ICCG 1e-03 0.01;
U BICCG 1e-03 0.1;
k BICCG 1e-03 0.1;
epsilon BICCG 1e-03 0.1;
R BICCG 1e-03 0.1;
nuTilda BICCG 1e-03 0.1;

}

***********************

I just want a 1e-03 solver tolerance.
But the computation still stopps at the endtime (=4000 itrations). I just have a solution wich seems to not menage to reach the convergence.

Can you give some advice to get it right?

Thanks in advance,

Lam
lam is offline   Reply With Quote

Old   April 5, 2007, 10:50
Default I don't wish to be impolite, b
  #6
Senior Member
 
Mark Olesen
Join Date: Mar 2009
Location: http://olesenm.github.io/
Posts: 777
Rep Power: 17
olesen will become famous soon enough
I don't wish to be impolite, but re-check the original suggestion:
# include "readSIMPLEControls.H"
...
if (simple.found("convergence"))
{
...
}

Check "readSIMPLEControls.H" file to see what is being read there and its relation to the following if (simple.found(..)) conditional.

Again, I really don't wish to be impolite, but I really don't have enough time to help much more than I already did. Sorry but you need to do this alone (or keep some basic commercial support).

/mark
olesen 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
URF and convergence in natural convergence Marie-Anne Main CFD Forum 11 September 11, 2009 11:07
convergence RMS anna CFX 7 October 3, 2006 05:50
convergence vijay FLUENT 6 February 1, 2006 04:04
convergence vijay Main CFD Forum 1 January 30, 2006 14:13
too bad convergence Davoche Main CFD Forum 2 November 20, 2005 06:08


All times are GMT -4. The time now is 05:26.