Exponential and Power Law Schemes
From ThermalFluidsPedia
Computational methodologies for forced convection 
Since the exact solution of eq. (4.201) exists, one can reasonably expect that an accurate scheme can be derived if the result of the exact solution, eq. (4.210), is utilized. Equation (4.201) can be rewritten as

Defining the total flux of due to convection and diffusion

eq. (4.229) becomes

Integrating eq. (4.231) over the control volume P (shaded area in Fig. 4.17), yields

Instead of assuming piecewise linear distribution of as with central difference scheme or assuming at the face of the control volume is equal to the value of at the grid point on the upwind side in the upwind scheme, the distribution of between grid points can be taken as that obtained from the exact solution, eq. (4.210). Applying eq. (4.210) between grid points E and P, we have

Substituting eq. (4.233) into eq. (4.230) and evaluating the result at x = x_{e}, the total flux of at the face of control volume becomes

Similarly, the total flux at the west face of the control volume is

Substituting eqs. (4.234) and (4.235) into eq. (4.232) and rearranging the resulting equation yields

where



Equations (4.237) and (4.238) can be rewritten in a format similar to that of eqs. (4.225) – (4.228), i.e.,


The comparison of a_{E} / D_{e} for different schemes is shown in Fig. 1. It can be seen that the hybrid scheme can be viewed as an envelope of the exponential scheme. The hybrid scheme is a good approximation if the absolute value of the grid Peclet number is either very large or near zero.
While the exponential scheme is accurate, the computational time is much longer than for the central difference, upwind or hybrid schemes. Patankar (1981) proposed a power law scheme that has almost the same accuracy as the exponential scheme but a substantially shorter computational time. The coefficient of the neighbor grid point on the east side can be obtained by
which can be rewritten in the following compact form
