Abstract
The physical content of Parlange's method of solving the flow equation is explored. His first approximation, which was developed by Macey (1959), satisfies continuity in the integral sense, but the second and higher approximations do not. Approximations beyond the first lack any constraining link between the separate steps of ‘satisfying continuity’ and ‘satisfying Darcy's law’, which make up each iteration. There is, consequently, nothing in the procedure to ensure convergence.A detailed investigation establishes the nonconvergence of the method when applied to one-dimensional sorption. It is found that the first approximation is best, and that the higher approximations make oscillations of increasing magnitude about the exact solution. Two illustrative examples are given. There is no reason to expect the procedure to be any more useful in other cases.The utility of Parlange's method is thus simply the utility of the first approximation: the dependence of this on the shape of the diffusivity function and on the flux-concentration relation is discussed.Contrary to Parlange's claim, the method cannot be applied to two-and three-dimensional systems other than radially symmetrical ones. The physical content of Parlange's method of solving the flow equation is explored. His first approximation, which was developed by Macey (1959), satisfies continuity in the integral sense, but the second and higher approximations do not. Approximations beyond the first lack any constraining link between the separate steps of ‘satisfying continuity’ and ‘satisfying Darcy's law’, which make up each iteration. There is, consequently, nothing in the procedure to ensure convergence. A detailed investigation establishes the nonconvergence of the method when applied to one-dimensional sorption. It is found that the first approximation is best, and that the higher approximations make oscillations of increasing magnitude about the exact solution. Two illustrative examples are given. There is no reason to expect the procedure to be any more useful in other cases. The utility of Parlange's method is thus simply the utility of the first approximation: the dependence of this on the shape of the diffusivity function and on the flux-concentration relation is discussed. Contrary to Parlange's claim, the method cannot be applied to two-and three-dimensional systems other than radially symmetrical ones. © Williams & Wilkins 1973. All Rights Reserved.