The modified characteristic finite difference fractional steps method for the coupled system of fluid dynamics in porous media and its analysis

For the coupled system of multilayer fluid dynamics in porous media, the modified characteristic finite difference fractional steps method applicable to parallel arithmetic is put forward and two‐dimensional and three‐dimensional schemes are used to form a complete set. Some techniques, such as calculus of variations, energy method, piecewise biquadratic interpolation, multiplicative commutation rule of difference operators, decomposition of high order difference operators and prior estimates are adopted. Optimal order estimates in L² norm are derived to determine the error in the approximate solution. This method has already been applied to the numerical simulation of multilayer fluid dynamics in porous media. © 2003 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 19: 665–681, 2003.