The Time Second Order Mass Conservative Characteristic FDM for Advection–Diffusion Equations in High Dimensions

We propose in this paper a time second order mass conservative algorithm for solving advection–diffusion equations. A conservative interpolation and a continuous discrete flux are coupled to the characteristic finite difference method, which enables using large time step size in computation. The advection–diffusion equations are first transformed to the characteristic form, for which the integration over the irregular tracking cells at previous time level is proposed to be computed using conservative interpolation. In order to get second order in time solution, we treat the diffusion terms by taking the average along the characteristics and use high order accurate discrete flux that are continuous at tracking cell boundaries to obtain mass conservative solution. We demonstrate the second order temporal and spatial accuracy, as well as mass conservation property by comparing results with exact solutions. Comparisons with standard characteristic finite difference methods show the excellent performance of our method that it can get much more stable and accurate solutions and avoid non-physical numerical oscillation.

[1]  Hongxing Rui,et al.  A conservative characteristic finite volume element method for solution of the advection–diffusion equation ☆ , 2008 .

[2]  T. F. Russell,et al.  An Eulerian-Lagrangian localized adjoint method for the advection-diffusion equation , 1990 .

[3]  Nigel Wood,et al.  SLICE: A Semi‐Lagrangian Inherently Conserving and Efficient scheme for transport problems , 2002 .

[4]  Andrew Staniforth,et al.  Application of the parabolic spline method (PSM) to a multi-dimensional conservative semi-Lagrangian transport scheme (SLICE) , 2007, J. Comput. Phys..

[5]  Philip John Binning,et al.  A Finite Volume Eulerian‐Lagrangian Localized Adjoint Method for Solution of the Contaminant Transport Equations in Two‐Dimensional Multiphase flow Systems , 1996 .

[6]  M. Zerroukat,et al.  A three‐dimensional monotone and conservative semi‐Lagrangian scheme (SLICE‐3D) for transport problems , 2012 .

[7]  P. Woodward,et al.  The Piecewise Parabolic Method (PPM) for Gas Dynamical Simulations , 1984 .

[8]  M. Memmesheimer,et al.  Modal aerosol dynamics model for Europe: development and first applications , 1998 .

[9]  Hongxing Rui,et al.  A Mass-Conservative Characteristic Finite Element Scheme for Convection-Diffusion Problems , 2008, J. Sci. Comput..

[10]  Hong Wang,et al.  A fractional step ELLAM approach to high-dimensional convection-diffusion problems with forward particle tracking , 2007, J. Comput. Phys..

[11]  Hongxing Rui,et al.  A second order characteristic finite element scheme for convection-diffusion problems , 2002, Numerische Mathematik.

[12]  Dong Liang,et al.  An efficient second‐order characteristic finite element method for non‐linear aerosol dynamic equations , 2009, International Journal for Numerical Methods in Engineering.

[13]  Peter Hansbo,et al.  The characteristic streamline diffusion method for convection-diffusion problems , 1992 .

[14]  K. Morton Numerical Solution of Convection-Diffusion Problems , 2019 .

[15]  Harley Flanders,et al.  Differentiation Under the Integral Sign , 1973 .

[16]  T. F. Russell,et al.  NUMERICAL METHODS FOR CONVECTION-DOMINATED DIFFUSION PROBLEMS BASED ON COMBINING THE METHOD OF CHARACTERISTICS WITH FINITE ELEMENT OR FINITE DIFFERENCE PROCEDURES* , 1982 .

[17]  F. Semazzi,et al.  Efficient Conservative Global Transport Schemes for Climate and Atmospheric Chemistry Models , 2002 .

[18]  Richard W. Healy,et al.  Solution of the advection-dispersion equation in two dimensions by a finite-volume Eulerian-Lagrangian localized adjoint method , 1998 .

[19]  J. Seinfeld,et al.  Atmospheric Chemistry and Physics: From Air Pollution to Climate Change , 1997 .

[20]  J. Bear Hydraulics of Groundwater , 1979 .

[21]  Mary F. Wheeler,et al.  Some improved error estimates for the modified method of characteristics , 1989 .

[22]  Rodolfo Bermejo,et al.  A Galerkin-characteristic algorithm for transport-diffusion equations , 1995 .

[23]  K. Aziz,et al.  Petroleum Reservoir Simulation , 1979 .