Central Difference Scheme
From ThermalFluidsPedia
Line 1:  Line 1:  
{{Comp Method for Forced Convection Category}}  {{Comp Method for Forced Convection Category}}  
  Integrating [[OneDimensional_SteadyState_Convection_and_Diffusion#equation_.282.29the governing equation]] over the control volume P (shaded area in  +  Integrating [[OneDimensional_SteadyState_Convection_and_Diffusion#equation_.282.29the governing equation]] over the control volume P (shaded area in the figure to the right), one obtains 
  [[Image:Fig4.17.pngthumb400 pxalt=Control volume for onedimensional problem   +  [[Image:Fig4.17.pngthumb400 pxalt=Control volume for onedimensional problem  Control volume for onedimensional problem.]] 
{ class="wikitable" border="0"  { class="wikitable" border="0"  
Line 12:  Line 12:  
{{EquationRef(1)}}  {{EquationRef(1)}}  
}  }  
  The righthand side of eq. (  +  The righthand side of eq. (1) can be obtained by assuming the distribution of <math>\varphi </math> between any two neighboring grid points is piecewise linear, i.e., 
<center><math>\left( \Gamma \frac{d\varphi }{dx} \right)_{e}=\Gamma _{e}\frac{\varphi _{E}\varphi _{P}}{(\delta x)_{e}}</math></center>  <center><math>\left( \Gamma \frac{d\varphi }{dx} \right)_{e}=\Gamma _{e}\frac{\varphi _{E}\varphi _{P}}{(\delta x)_{e}}</math></center>  
Line 19:  Line 19:  
<center><math>\left( \Gamma \frac{d\varphi }{dx} \right)_{w}=\Gamma _{w}\frac{\varphi _{P}\varphi _{W}}{(\delta x)_{w}}</math></center>  <center><math>\left( \Gamma \frac{d\varphi }{dx} \right)_{w}=\Gamma _{w}\frac{\varphi _{P}\varphi _{W}}{(\delta x)_{w}}</math></center>  
  where  +  where Γ<sub>e</sub> and Γ<sub>w</sub> are the diffusivities at the faces of the control volume. To ensure that the flux of <math>\varphi </math> across the faces of the control volume is continuous, the harmonic mean diffusivity at the faces should be used. To evaluate the left hand side of eq. (1), it is necessary to know the values of <math>\varphi </math> at the faces of the control volume. If the piecewise linear profile of <math>\varphi </math> is chosen, it follows that 
<center><math>\varphi _{e}=\frac{\varphi _{E}+\varphi _{P}}{2}</math></center>  <center><math>\varphi _{e}=\frac{\varphi _{E}+\varphi _{P}}{2}</math></center>  
Line 26:  Line 26:  
<center><math>\varphi _{w}=\frac{\varphi _{P}+\varphi _{W}}{2}</math></center>  <center><math>\varphi _{w}=\frac{\varphi _{P}+\varphi _{W}}{2}</math></center>  
  Therefore, eq. (  +  Therefore, eq. (1) becomes 
{ class="wikitable" border="0"  { class="wikitable" border="0" 
Revision as of 03:17, 21 July 2010
Computational methodologies for forced convection 
Integrating the governing equation over the control volume P (shaded area in the figure to the right), one obtains

The righthand side of eq. (1) can be obtained by assuming the distribution of between any two neighboring grid points is piecewise linear, i.e.,
where Γ_{e} and Γ_{w} are the diffusivities at the faces of the control volume. To ensure that the flux of across the faces of the control volume is continuous, the harmonic mean diffusivity at the faces should be used. To evaluate the left hand side of eq. (1), it is necessary to know the values of at the faces of the control volume. If the piecewise linear profile of is chosen, it follows that
Therefore, eq. (1) becomes

Defining the mass flux and diffusive conductance

eq. (4.212) can be rearranged as

where



This scheme is termed the central difference scheme because the values of at the faces of the control volume are taken as the averaged value between two grid points. The continuity equation (4.202) requires that F_{e} = F_{w} and therefore, eq. (4.217) reduces to
a_{P} = a_{W} + a_{E}
To evaluate the performance of the central difference scheme, let us consider the case of a uniform grid, i.e., (δx)_{e} = (δx)_{w} = δx, for which case eq. (4.212) can be rearranged as

where

is the Peclet number using grid size as the characteristic length, which is referred to as the grid Peclet number. The grid Pe is a ratio of the strength of convection over diffusion. To ensure stability of the discretization scheme, the value of should always fall between and , which requires that the coefficients, and , are positive, i.e.,

This is the criterion for stability of the central difference scheme. It can be demonstrated that the central difference becomes unstable if eq. (4.220) is violated (see Problem 4.22). The fact that the central difference scheme is stable under small grid Peclet number indicates that the central difference scheme is accurate only if the convection is not very significant.