# DIFFUSION free schemes

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

 May 12, 2004, 07:59 DIFFUSION free schemes #1 question Guest   Posts: n/a my lecturer says the scheme he uses if diffusion free. what does it mean by diffusion free, and which schemes are and which are not diffsion free? also he seemed to imply it had something to do with the type of mesh? can some one please explain. thanks

 May 12, 2004, 13:08 Re: DIFFUSION free schemes #2 centaur_ks Guest   Posts: n/a i guess by diffusion free he means Numerical diffusion free

 May 13, 2004, 01:11 Re: DIFFUSION free schemes #3 agg Guest   Posts: n/a The higher the grid resolution, the lower the diffusion. The higher the order of accuracy of discretization, the lower the diffusion.

 May 13, 2004, 09:39 Re: DIFFUSION free schemes #4 Li Yang Guest   Posts: n/a As far as I know, the central difference scheme is less diffusive and upwind schemes are quite diffusive. But is there any scheme which has no numerical diffusion ? I have never heard of. Li

 May 13, 2004, 21:00 Re: DIFFUSION free schemes #5 Chris Sideroff Guest   Posts: n/a In determining the numerical behavior of a scheme you must look at order of the truncated terms (from the discretization) There are two types of numerical error - dissipation (diffusion) and dispersion (oscillatory) Even order truncatation tends to be dissipative Odd order truncation tends to be dispersive First orders schemes (first order upwind, Lax scheme, etc) have second order truncation error thus they are dissipative. This is why obtaining solutions tend to be easy for first order schemes but are notoriously inaccurate unless you use VERY fine grids. Second order schemes (central difference, second order upwind) have third order truncation error so they tend to be disperive or oscillatory. Due to this solutions can be difficult to obtain. Usually modifications are made to damp out the oscillations: artificial viscosity for central differencing and TVD and ENO schemes for second order upwinding. The same holds true for higher order schemes. The difficulty with higher order schemes is maintaining a uniform accuracy on the boundaries (because the stencil gets large). Hope that clears thing up - Oh and to answer your inital question NO CFD algorithm is without numerical dissipation - even the best schemes contain a little

 May 17, 2004, 03:14 Is ideal DIFFUSION free schemes possible? #6 Zeng Guest   Posts: n/a I have two genernal questions: (1) Is really a diffusion free scheme possible? (2) What wiil be happen if we choose a high-order scheme with lower-order boundary ? I mean, for example, we use a 4-order scheme, but in boundary condition, it is only 2-order precision, what is the effect for the lower-order scheme in boundary? Thanks your responding in advance!

 May 17, 2004, 04:15 Re: Is ideal DIFFUSION free schemes possible? #7 Praveen Guest   Posts: n/a As others have written before some numerical diffusion is necessary if you want to get reasonable numerical solution especially for convection-dominated systems. It is enough if the boundary scheme is one order less accurate than the interior scheme. The overall accuracy will still be equal to that of the interior scheme. I cant recall where I have read this, probably in a paper by D. W. Zingg. Hope somebody else can point to some references.

 May 19, 2004, 05:02 Re: DIFFUSION free schemes #8 Onno Guest   Posts: n/a Be careful Chris; In answer to your answer: " ... NO CFD algorithm is without numerical dissipation - even the best schemes contain a little" The Hyper-C scheme of Leonard see 1991 Comp. Meth. in Appl. Mech. Eng. Vol. 88, p.17-74 is completely diffusion free, but pretty useless as it turns any moderate gradient into a step profile.

 May 24, 2004, 14:23 Re: Is ideal DIFFUSION free schemes possible? #9 QK Guest   Posts: n/a I'm not 100% sure.. but, Hirsh's book discusses the issue on the interior/boundary accuarcy. It should have a pointer to some papers. I kind of remember seeing a diff. free schemes in Lohner's book. It supposedly makes the dissipation term propotional to residual. I feel like an idiot.. Time to read/study more.

 May 25, 2004, 01:48 Re: Is ideal DIFFUSION free schemes possible? #10 Praveen Guest   Posts: n/a In A residual-based compact scheme for the compressible Navier-Stokes equations by Alain Lerat and Christophe Corre a centered scheme is modified by adding dissipation that depends on the residual.
... the process for getting a high accuracy as well as the construction of the numerical dissipation rely strongly on the residual vanishing at steady state. Starting from a basic centered scheme, the idea for increasing the accuracy order is not to cancel the leading term of the truncation error, but to add a new error term so that the resulting truncation error be expressed in terms of derivatives of the residual only.
Even though the residual vanishes at steady state, the scheme still has dissipation.
For simplicity and robustness, we use a first-order dissipation operator, which is compatible with any accuracy order. More precisely, by applying this operator, a centered scheme accurate at order 2p becomes a dissipative scheme accurate at order 2p-1 (dissipation at steady state comes from the truncation error of the dissipation operator).
I haven't come across any scheme which is totally free of dissipation. Would be interested to hear about such schemes. If you want the above paper then mail me at cpraveen[at]gmail.com

 June 1, 2004, 05:04 Re: Is ideal DIFFUSION free schemes possible? #11 Zeng Guest   Posts: n/a Thanks for all of your kindly reply. Zeng

 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 Niklas Wikstrom (Wikstrom) OpenFOAM Running, Solving & CFD 122 June 15, 2014 06:20 botp OpenFOAM 2 March 11, 2011 16:27 therockyy FLOW-3D 1 June 20, 2010 19:36 QU KUN Main CFD Forum 0 January 22, 2004 02:21 Viatcheslav Anissimov CFX 0 April 3, 2002 06:27

All times are GMT -4. The time now is 14:02.