[lam]

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]

Oh, THANKS a lot Mark, for replying me!

I am going give a go right now!

Cheers,

Lam

[olesen]

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]

It looks mostly reasonable except you don't need this line twice.

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

[lam]

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

[olesen]

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