Affiliation:
1. State University of New York at Binghamton
Abstract
A new technique for the development of finite difference schemes for diffusion equations is presented. The model equations are the one space variable advection diffusion equation and the two space variable diffusion equation, each with Dirichlet boundary conditions. A two-step hybrid technique, which combines perturbation methods, based on the parameter ρ = Δt / (Δx)2, with the Galerkin method, provides a systematic way to develop new finite difference methods, referred to as hybrid equations. The main contributions of this paper include: 1) recovery of classical explicit or implicit finite difference schemes using only the perturbation terms; 2) development of new finite difference schemes, referred to as hybrid equations, which have better stability properties than the classical finite difference equations, permitting the use of larger values of the parameter ρ; and 3) higher order accurate methods, with either O((Δx)4) or O((Δx)6) truncation error, formed by convex linear combinations of the classical and hybrid equations. The solution of the hybrid finite difference equations requires only a tridiagonal equation solver and, hence, does not lead to excessive computational effort.
Reference17 articles.
1. Patankar, S., 1980, Numerical Heat Transfer and Fluid Flow, Hemisphere, New York.
2. Roache, P., 1972, Computational Fluid Dynamics, Kenmore Publishers, Albuquerque, NM.
3. Fillo, J. and Geer, J., 1999, “A New Numerical Technique for the Solution to the Heat Equation Based on a Hybrid Perturbation Galerkin Method,” Proceedings of the 33rd National Heat Transfer Conference, Albuquerque, NM.
4. Anderson
C. M.
, and GeerJ. F., 1991, “Investigating a Hybrid Perturbation Galerkin Technique Using Computer Algebra,” J. Symb. Comp., Vol. 12, pp. 695–714.
5. Geer
J. F.
and AndersenC. M., 1989, “A Hybrid Perturbation Galerkin Technique for Differential Equations Containing a Parameter,” Appl. Mech. Rev., Vol. 42, no. 11, pp. 569–577.