# Transient Error Estimate with icoErrorEstimate ?

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

 June 25, 2009, 13:54 Transient Error Estimate with icoErrorEstimate ? #1 New Member   Christophe Kassiotis Join Date: Mar 2009 Location: Paris Posts: 17 Rep Power: 8 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 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. CK

 June 28, 2009, 05:24 #2 Senior Member   Hrvoje Jasak Join Date: Mar 2009 Location: London, England Posts: 1,763 Rep Power: 21 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" Enjoy, Hrv hua1015 likes this. __________________ Hrvoje Jasak Providing commercial FOAM/OpenFOAM and CFD Consulting: http://wikki.co.uk

 Thread Tools Display Modes Linear Mode

 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 OffTrackbacks are On Pingbacks are On Refbacks are On Forum Rules

 Similar Threads Thread Thread Starter Forum Replies Last Post Ijaz CFX 0 February 19, 2009 17:27 Luk Main CFD Forum 0 October 19, 2007 10:09 Adam CFX 1 April 12, 2007 11:34 Basics CFX 3 September 25, 2002 09:42 mahesh prakash Main CFD Forum 1 January 21, 1999 14:45

All times are GMT -4. The time now is 10:55.