|July 23, 2009, 08:19||
divergence with higher order scheme
Join Date: Jul 2009
Posts: 5Rep Power: 5
I am doing DNS and am in process of implementing higher order schemes. I am done with implementing it in momentum equation etc.
However, when it comes to calculate divergence in a cell, i.e. Div = du/dx+dv/dy+dw/dz, I am getting a bit confused whether I should implement higher order scheme here as well or I should keep it just second order accurate.
I am asking it because if 2nd order accurate,
du/dx = (u[i+1]-u[i])/(x[i+1]-x[i]) etc. on staggered grids. It gives a feeling of the rate of mass accumulation in the cell, a physically meaningful quantity.
However, if we change this to 4th order discretization,
du/dx = f(u[i+2], u[i+1], u[i], u[i-1]) etc. which doesn't look like giving mass accumulation in the cell in discrete manner.
Could anyone please point out which one should be preferred?
|July 23, 2009, 13:53||
Join Date: Mar 2009
Posts: 62Rep Power: 6
Usually you should keep the same order of discretization on all spatial derivatives if you want to have a pure 4 order scheme let's say. Otherwise your DNS will be only 2nd order.
|divergence, higher order|
|Thread||Thread Starter||Forum||Replies||Last Post|
|Higher order terms in the equations||Argyle||Main CFD Forum||7||October 31, 2008 08:41|
|CFL condition for higher order schemes||Shyam||Main CFD Forum||2||February 14, 2008 14:24|
|Higher order downwind scheme||jelmer||OpenFOAM Running, Solving & CFD||4||August 9, 2006 06:43|
|Higher order discretization on staggered grid||Chandra Shekhar||Main CFD Forum||9||January 27, 2005 16:31|
|Discretisation scheme in CFX-1st or 2nd order?||Pete||CFX||10||January 12, 2005 12:48|