The coefficients of equation (?) must remain positive, since an increase in a neighboring value must result in an increase in to obtain physically realistic solutions. Thus, the inequalities and must be satisfied. The Peclet number is the ratio between convective strength and diffusive conductance. To achieve physically realistic solutions, the inequality
must be satisfied. The parameter is defined by the chosen scheme, depending on Equation (1). The various differencing schemes are:
so that satisfies Equation (1). Thus, the central differencing scheme is only numerically stable for a low values of .
This formulation can be derived from the exact solution, and thus, guarantees positive coefficients while not over-predicting the diffusive terms. However, the computation of exponentials is slow and therefore a faster scheme is generally used, especially in higher dimensions.
The hybrid scheme is formulated by allowing , and in the exponential scheme. The hybrid scheme is an improvement on the upwind scheme, however, it deviates from the exponential scheme at .
The power law scheme overcomes the inaccuracies of the hybrid scheme, while improving on the computational time for the exponential scheme.
VanLeerConvectionTerm not mentioned and no discussion of explicit forms.
All of the numerical schemes presented here are available in FiPy and can be selected by the user.