CFD Online Logo CFD Online URL
Home > Forums > OpenFOAM Post-Processing

Transient Error Estimate with icoErrorEstimate ?

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

Like Tree1Likes
  • 1 Post By hjasak

LinkBack Thread Tools Display Modes
Old   June 25, 2009, 13:54
Default Transient Error Estimate with icoErrorEstimate ?
New Member
Christophe Kassiotis
Join Date: Mar 2009
Location: Paris
Posts: 17
Rep Power: 10
kassiotis is on a distinguished road
Hello everybody

I would like to estimate the errors in a transient computation (for instance flow around a cylinder bluff body) with low Reynolds number (around 200). So the idea was to use an error estimate like icoErrorEstimate.

I launch some computations and I observe that the error is at it biggest behind the bluff body,for few points in the Van Karman vortex street. With a rather coarse mesh, the error estimate give something around 1.8 in relative value (mag(e)/mag(U)) for some of this points. Elsewhere, the error is rather small (less than 10%) - but the flow is quasi static elsewhere.

My guess was that the transient term is not estimate corectly with icoErrorEstimate, but maybe i'm wrong.
Indeed, the equation solved is:

errorEstimate<vector> ee( resError::div(phi, U) - resError::laplacian(nu, U) == -fvc::grad(p) );

According to the PhD of H. Jasak, p 196-200, one should take care of the transient term, adding something like (Eq 4.117)

\frac{\rho\Phi-\rho\Phi^0}{\Delta t} to this equation

So do I need to add this term ? Is this normal that nothing like this appears in icoErrorEstimate ? Maybe I didn't get something...

Other questions : if I want a global error estimate of the computation error, should I consider the max/mean total/relative error ? Is it possible to get an error approximate for the pressure field ?

Thank you in advance.

kassiotis is offline   Reply With Quote

Old   June 28, 2009, 05:24
Senior Member
Hrvoje Jasak
Join Date: Mar 2009
Location: London, England
Posts: 1,802
Rep Power: 24
hjasak will become famous soon enough
Formally, you need to separate the temporal error from its spatial counterpart.

For the spatial error, consider the ddt as a fixed source term and use the residual estimate as before.

For the temporal component, you need some form of Richardson extrapolation, depending on how much reliable temporal data you have got.

Please note: you are NOT allowed to sum up the two errors for the "total error"


hua1015 likes this.
Hrvoje Jasak
Providing commercial FOAM/OpenFOAM and CFD Consulting:
hjasak 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
Transient Thermal Stress Distribution Ijaz CFX 0 February 19, 2009 17:27
How are a transient simulations performed?? Luk Main CFD Forum 0 October 19, 2007 10:09
Transient vs Steady State Adam CFX 1 April 12, 2007 11:34
transient simulation: natural convection problem? Basics CFX 3 September 25, 2002 09:42
difference between false and true transient mahesh prakash Main CFD Forum 1 January 21, 1999 14:45

All times are GMT -4. The time now is 11:44.