CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > General Forums > Main CFD Forum

Discretization of the diffusion term

Register Blogs Community New Posts Updated Threads Search

Like Tree3Likes
  • 1 Post By selig5576
  • 2 Post By FMDenaro

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   August 29, 2017, 22:16
Default Discretization of the diffusion term
  #1
New Member
 
Join Date: Jul 2017
Posts: 6
Rep Power: 9
mseka is on a distinguished road
For incompressible NS equation in the finite volume discretization on a staggered grid, I have encountered a doubt.

We can put the diffusive term together with the advection term by factoring out gradient, for instance for x-direction for one term as below:
\frac{\partial }{\partial x}\left ( u^2-2 \mu \frac{\partial u}{\partial x} \right )

I am fine with the way to discretize the second term in staggered grid but for the first term in the code that I see someone has discretized the first term as
\frac{u^2_{i+1,j}-u^2_{i-1,j}}{2\Delta x }

Why is this correct? Even if it is the central difference, shouldn't it also include u_{ij}?
mseka is offline   Reply With Quote

Old   August 29, 2017, 22:29
Default Diffusion
  #2
Senior Member
 
Selig
Join Date: Jul 2016
Posts: 213
Rep Power: 11
selig5576 is on a distinguished road
The finite difference discretization below ensures proper conservation.
\frac{u_{i+1}^{2} - u_{i-1}^{2}}{2 \Delta x}

The advection term of the form
u_{i} \frac{u_{i+1} - u_{i-1}}{2 \Delta x}

is not conservative and can lead to the numerical solution undershooting. As a side note, central differencing like that of the advection term can lead spurious oscillations. QUICK scheme or a stable advection scheme (ENO, WENO) can remedy such a problem.
selig5576 is offline   Reply With Quote

Old   August 29, 2017, 22:36
Default
  #3
New Member
 
Join Date: Jul 2017
Posts: 6
Rep Power: 9
mseka is on a distinguished road
So you mean all is good with the conservative form of discretization? And no oscillations occur?
mseka is offline   Reply With Quote

Old   August 29, 2017, 23:02
Default Conservative schemes
  #4
Senior Member
 
Selig
Join Date: Jul 2016
Posts: 213
Rep Power: 11
selig5576 is on a distinguished road
Yes, you will want your scheme to be conservative. In terms of your second question, spurious oscillations can occur regardless if your scheme is conservative or not. A simple 1D example would be looking at the convection-diffusion equation and adjust the peclet number. You will find that upwind differencing will perform better than central differencing.
FMDenaro likes this.
selig5576 is offline   Reply With Quote

Old   August 30, 2017, 04:59
Default
  #5
Senior Member
 
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 6,882
Rep Power: 73
FMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura about
I agree, consider the Burgers equation in divergence and non-divergence (quasi-linear) form:

du/dt + d(u^2/2)/dx=d(mu*du/dx)/dx (divergence form)

du/dt + u*du/dx=mu*d2u/dx^2 (quasi linear)

No difference exists for the staggered or non staggered grid since the only involved variable is u.
Discretize with second order central difference a see what happens. The appearence of u(i,j) in the convective term:

d(u^2/2)/dx -> (u(i+1/2)^2-u(i-1/2)^2)/2/h
du/dx -> (u(i+1/2)-u(i-1/2))/h

Now you have two chance:
1) linear interpolation for u on the half-nodes
2) linear interpolation for u^2 on the half nodes

The first choice is congruent to a linear reconstruction of the flux function, the second one introduces aliasing in the terms.
FMDenaro is offline   Reply With Quote

Old   August 30, 2017, 11:41
Default
  #6
Senior Member
 
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 6,882
Rep Power: 73
FMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura about
Quote:
Originally Posted by mseka View Post
So you mean all is good with the conservative form of discretization? And no oscillations occur?
numerical oscillations are due to the lack of monotonicity of the scheme... Godunov theores states that monotonic linear scheme can be only first order accurate. Have also a look to the examples in the book of Peric & Ferziger
FMDenaro is offline   Reply With Quote

Reply

Tags
discretization scheme, navier stokes


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 Off
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
Why numerical diffusion is related to the advection term Houthuys Main CFD Forum 6 January 28, 2017 23:01
Moving mesh Niklas Wikstrom (Wikstrom) OpenFOAM Running, Solving & CFD 122 June 15, 2014 07:20
FVM discretization of diffusion term on crvlnr gr Michail Main CFD Forum 3 March 14, 2008 07:52
Discretization of the Diffusion Term Rui CFX 0 October 24, 2005 13:43


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