CFD Online Logo CFD Online URL
Home > Forums > OpenFOAM Running, Solving & CFD


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

LinkBack Thread Tools Display Modes
Old   October 14, 2013, 06:17
Question simpleFoamResidual
Senior Member
M. Montero
Join Date: Mar 2009
Location: Madrid
Posts: 125
Rep Power: 10
be_inspired is on a distinguished road
Hi all,

My question is related to the utility simpleFoamResidual implemented into OF1.6ext. According to the code, it analyze the momentum equation. The perfect situation, I think, would be that uResidual is equal "0".

The attached convergence is OK but the maxium uResidual that I obtain is around 76 and the mean around 0.09. In other different mesh, the maximum uResidual obtained was around 2000.

Could anyone tell me how much is a high value of uResidual? Is possible to have a physical meaning about the uResidual? If I have higher values than 0.05, for example, does it tell me that it is require to refine the mesh?

I would appreciate any help. I would like to have a criteria to know if current mesh resolution is enough or not for my purposes without doing a mesh resolution analysis in detail. Thank you very much.

            Info<< "    Calculating uResidual" << endl;
            volScalarField uResidual
                    fvc::div(phi, U)
                  + fvc::div(turbulence->R())
                  + fvc::grad(p)
Time = 3000

DILUPBiCG:  Solving for Ux, Initial residual = 1.48121e-08, Final residual = 6.66887e-11, No Iterations 2
DILUPBiCG:  Solving for Uy, Initial residual = 1.09835e-05, Final residual = 1.32954e-08, No Iterations 2
DILUPBiCG:  Solving for Uz, Initial residual = 1.96974e-05, Final residual = 3.96294e-08, No Iterations 2
GAMG:  Solving for p, Initial residual = 0.00085226, Final residual = 4.08688e-05, No Iterations 8
time step continuity errors : sum local = 2.32144e-10, global = -6.39797e-11, cumulative = 2.45607e-07
DILUPBiCG:  Solving for epsilon, Initial residual = 2.62941e-07, Final residual = 2.54438e-10, No Iterations 1
DILUPBiCG:  Solving for k, Initial residual = 9.73059e-07, Final residual = 3.07595e-10, No Iterations 2
ExecutionTime = 7790.51 s  ClockTime = 7792 s
be_inspired is offline   Reply With Quote

Old   October 14, 2013, 06:26
Senior Member
akidess's Avatar
Anton Kidess
Join Date: May 2009
Location: Germany
Posts: 1,274
Rep Power: 23
akidess will become famous soon enough
I don't think this tool can replace a mesh sensitivity study. All it tells you is that a solution is converged, not that it is converged and accurate.
*On twitter @akidTwit
*Spend as much time formulating your questions as you expect people to spend on their answer.
akidess is offline   Reply With Quote

Old   October 14, 2013, 07:10
Senior Member
M. Montero
Join Date: Mar 2009
Location: Madrid
Posts: 125
Rep Power: 10
be_inspired is on a distinguished road
I am trying to link Dr.Jasak Thesis about error estimation with my case and to know if it is needed to refine the mesh once convegence has been "obtained" (residuals).

Chapter 5 is dedicated to the adaptative mesh refinement to reduce the numerical error up to the desired level. There should be a criteria ( I hope that) about when the error (uResidual) has been reduced enough to minimize the numerical error.
be_inspired is offline   Reply With Quote


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
Custom Boundary Condition with OpenFOAM trusharg OpenFOAM Programming & Development 44 Yesterday 23:39
simpleFoamResidual File OpenFOAM Running, Solving & CFD 2 October 4, 2013 10:50
simpleFoamResidual error val46 OpenFOAM 1 October 19, 2010 02:47

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