Comparison of finite difference and control volume methods for solving differential equations

Abstract Comparisons are made between the finite difference method (FDM) and the control volume formulation (CVF). An analysis of truncation errors for the two methods is presented. Some rules-of-thumb related to the accuracy of the methods are included. It is shown that the truncation error is the same for both methods when the boundary conditions are of the Dirichlet type, the system equations are linear and represented in Cartesian coordinates. A technique to analyze the accuracy of the methods is presented. Two examples representing different physical situations are solved using the methods. The FDM failed to conserve mass for a small number of nodes when both boundary conditions include a derivative term (i.e. either a Robin or Neumann type boundary condition) whereas the CVF method did conserve mass for these cases. The FDM is more accurate than the CVF for problems with interfaces between adjacent regions. The CVF is (ΔX) order of accuracy for a Neumann type boundary condition whereas the FDM is (ΔX) 2 order.

[1]  Mikhail Shashkov,et al.  The sensitivity and accuracy of fourth order finite-difference schemes on nonuniform grids in one dimension , 1995 .

[2]  J. Newman.  Numerical Solution of Coupled, Ordinary Differential Equations , 1968 .

[3]  E. Christiansen,et al.  Handbook of Numerical Analysis , 1996 .

[4]  T. Fuller,et al.  Influence of rib spacing in proton-exchange membrane electrode assemblies , 1996 .

[5]  Claus-Dieter Munz,et al.  A finite-volume particle-in-cell method for the numerical treatment of Maxwell-Lorentz equations on boundary-fitted meshes , 1999 .

[6]  Richard H. Stark.  Rates of Convergence in Numerical Solution of the Diffusion Equation , 1956, JACM.

[7]  S. Patankar,et al.  Nonaxisymmetric Radiative Transfer in Cylindrical Enclosures , 1996, Heat Transfer: Volume 1 — Heat Transfer in Microgravity Systems; Radiative Heat Transfer and Radiative Heat Transfer in Low-Temperature Environments; Thermal Contact Conductance and Inverse Problems in Heat Transfer.

[8]  B. P. Leonard.  Order of accuracy of QUICK and related convection-diffusion schemes , 1995 .

[9]  Stefan T. Thynell,et al.  Discrete-ordinates method in radiative heat transfer , 1998 .

[10]  J. V. Zee,et al.  Extension of Newman’s Numerical Technique To Pentadiagonal Systems Of Equations , 1984 .

[11]  Lars Schiøtt Sørensen,et al.  An introduction to Computational Fluid Dynamics: The Finite Volume Method , 1999 .

[12]  Joel H. Ferziger,et al.  Introduction to Theoretical and Computational Fluid Dynamics , 1996 .

[13]  René A. J. Janssen,et al.  Electrochemical Society Proceedings , 2000 .

[14]  Anthony Ralston,et al.  Mathematical Methods for Digital Computers , 1960 .

[15]  Sam Qian,et al.  Wavelets and the numerical solution of boundary value problems , 1993 .

[16]  Yannis Kallinderis,et al.  A 3D Finite-Volume Scheme for the Euler Equations on Adaptive Tetrahedral Grids , 1994 .

[17]  B. P. Leonard.  Comparison of truncation error of finite-difference and finite-volume formulations of convection terms , 1994 .

[18]  C. Chan,et al.  Solution of incompressible flows with or without a free surface using the finite volume method on unstructured triangular meshes , 1999 .

[19]  Mark E. Davis,et al.  Numerical methods and modeling for chemical engineers , 1984 .

[20]  L. Davidson,et al.  Simulation of three-dimensional flow around a square cylinder at moderate Reynolds numbers , 1999 .

[21]  E. Cuthill,et al.  Numerical Methods for Nuclear Reactor Calculations , 1961 .

[22]  Hong-Koo Noh Noh,et al.  Temperature distribution of a low temperature heat pipe with multiple heaters for electronic cooling , 1998 .

[23]  H. Dwyer,et al.  Fluid Flow and Heat Transfer Over a Three-Dimensional Spherical Object in a Pipe , 1998 .

[24]  D. Fan,et al.  Optimization and extension of pentadiagonal band(J) solver to multiregion systems containing interior boundaries , 1991 .

[25]  M. D. Su,et al.  Numerical simulation of fluid flow and thermal performance of a dry-cooling tower under cross wind condition , 1999 .

[26]  Seung Wook Baek,et al.  Nonorthogonal finite-volume solutions of radiative heat transfer in a three-dimensional enclosure , 1998 .

[27]  A new finite volume method for the solution of convection–diffusion equations: analysis of stability and convergence , 1997 .

[28]  R. Varga.  Numerical Solution of the Two-Group Diffusion Equations in X-Y Geometry , 1957, IRE Transactions on Nuclear Science.

[29]  J. Gillis,et al.  Matrix Iterative Analysis , 1961 .

[30]  Robert J. MacKinnon,et al.  Equivalent versions of the quick scheme for finite‐difference and finite‐volume numerical methods , 1992 .

[31]  B. P. Leonard,et al.  Beyond first‐order upwinding: The ultra‐sharp alternative for non‐oscillatory steady‐state simulation of convection , 1990 .

[32]  James P. Braselton,et al.  The Maple V handbook , 1994 .