Global spectral analysis for convection-diffusion-reaction equation in one and two-dimensions: Effects of numerical anti-diffusion and dispersion

Abstract Convection-diffusion-reaction (CDR) equation plays a central role in many disciplines of engineering, science and finances. As a consequence, importance of analysis of numerical methods for the accurate solution of CDR equation has motivated the present research. We have used the global spectral analysis to characterize all the three important physical processes in terms of the non-dimensional numerical parameters, namely, the non-dimensional wavenumber (kh); Courant-Friedrich-Lewy (CFL) number, N c ; the Peclet number (Pe) and the Damkohler number (Da). For the purpose of illustration, we have focused on two space-time discretization schemes known for accuracy and robustness. The basic properties relate to numerical issues arising for the numerical amplification factor, numerical diffusion coefficient, numerical phase speed and numerical group velocity. With the help of model one-dimensional (1D) and two-dimensional (2D) CDR equations, we have reported the numerical property charts for the cases: (i) When all the processes of convection, diffusion and reaction are of same order, with critical numerical behaviour enforcing low values of Da for the 1D CDR equation studied here. (ii) The 2D CDR equation considered is diffusion-reaction dominated, and as a consequence, this enforces Da to be larger. We have thoroughly analyzed these cases to identify the essential roles of anti-diffusion on the critical N c and Pe values, which in turn decides admissible space and time steps to be used with the discretization schemes. The property charts have been used to calibrate the analysis with two model equations, one of which has an exact solution for a 1D CDR equation, and the second case for the 2D CDR equation has numerical solution available in the literature. These cases help to identify the importance of such analysis in explaining the utility of the choice one can exercise in fixing the numerical parameters. This also identifies and explains some hitherto unknown numerical problems for CDR equation and their alleviation techniques.

[1]  D. Drikakis,et al.  Large eddy simulation of compressible turbulence using high‐resolution methods , 2005 .

[2]  Éloïse Comte,et al.  Optimal Buffer Zone for the Control of Groundwater Pollution from Agricultural Activities , 2019, Applied Mathematics & Optimization.

[3]  R. Takaki,et al.  Thermal convection of a fluid with temperature-dependent viscosity , 1992 .

[4]  L. Kramer,et al.  Travelling-stripe forcing of Turing patterns , 2004 .

[5]  L. Leboucher,et al.  Monotone Scheme and Boundary Conditions for Finite Volume Simulation of Magnetohydrodynamic Internal Flows at High Hartmann Number , 1999 .

[6]  P. Moin,et al.  Progress-variable approach for large-eddy simulation of non-premixed turbulent combustion , 2004, Journal of Fluid Mechanics.

[7]  Andrew J. Wathen,et al.  A moving grid finite element method applied to a model biological pattern generator , 2003 .

[8]  Jean-Philippe Nicot,et al.  Convection-diffusion-reaction of CO2-enriched brine in porous media: A pore-scale study , 2019, Comput. Geosci..

[9]  Tapan K. Sengupta,et al.  Further improvement and analysis of CCD scheme: Dissipation discretization and de-aliasing properties , 2009, J. Comput. Phys..

[10]  Tapan K. Sengupta,et al.  Analysis of central and upwind compact schemes , 2003 .

[11]  Z. Mei.  Numerical Bifurcation Analysis for Reaction-Diffusion Equations , 2000 .

[12]  Tony W. H. Sheu,et al.  An implicit scheme for solving the convection-diffusion-reaction equation in two dimensions , 2000 .

[13]  R. Codina.  Comparison of some finite element methods for solving the diffusion-convection-reaction equation , 1998 .

[14]  T. Hughes,et al.  Finite element methods for the Helmholtz equation in an exterior domain: model problems , 1991 .

[15]  T. K. Sengupta,et al.  Error dynamics: Beyond von Neumann analysis , 2007, J. Comput. Phys..

[16]  Tapan K. Sengupta,et al.  Error dynamics of diffusion equation: Effects of numerical diffusion and dispersive diffusion , 2014, J. Comput. Phys..

[17]  Dominique Pelletier,et al.  Positivity Preservation and Adaptive Solution for the k-? Model of Turbulence , 1998 .

[18]  Tapan K. Sengupta,et al.  Analysis of anisotropy of numerical wave solutions by high accuracy finite difference methods , 2011, J. Comput. Phys..

[19]  Tapan K. Sengupta,et al.  Spurious waves in discrete computation of wave phenomena and flow problems , 2012, Appl. Math. Comput..

[20]  Vaibhav Joshi,et al.  A positivity preserving variational method for multi-dimensional convection-diffusion-reaction equation , 2017, J. Comput. Phys..

[21]  Tapan K. Sengupta,et al.  A new combined stable and dispersion relation preserving compact scheme for non-periodic problems , 2009, J. Comput. Phys..

[22]  V. K. Suman,et al.  Focusing Phenomenon in Numerical Solution of Two-Dimensional Navier–Stokes Equation , 2019, High-Performance Computing of Big Data for Turbulence and Combustion.

[23]  A. Turing.  The chemical basis of morphogenesis , 1952, Philosophical Transactions of the Royal Society of London. Series B, Biological Sciences.

[24]  B. Rothschild,et al.  Population-dynamic instability as a cause of patch structure , 1996 .

[25]  Enzo Tiezzi,et al.  Dynamics of pattern formation in biomimetic systems. , 2008, Journal of theoretical biology.

[26]  Gabriel R. Barrenechea,et al.  A stabilised finite element method for the convection–diffusion–reaction equation in mixed form , 2018, Computer Methods in Applied Mechanics and Engineering.

[27]  Anotida Madzvamuse,et al.  Pigmentation pattern formation in butterflies: experiments and models. , 2003, Comptes rendus biologies.

[28]  Tsing-Fa Lin,et al.  Visualization of roll patterns in Rayleigh–Bénard convection of air in a rectangular shallow cavity , 2001 .

[29]  Satoru Ikeuchi,et al.  Formation of dissipative structures in galaxies , 1984 .

[30]  F. Busse,et al.  Thermal convection in rotating spherical shells , 1997 .

[31]  Diego Alexander Garzón-Alvarado,et al.  Computational examples of reaction–convection–diffusion equations solution under the influence of fluid flow: First example , 2012 .

[32]  Tapan K. Sengupta,et al.  High Accuracy Computing Methods: Fluid Flows and Wave Phenomena , 2013 .

[33]  Clint Dawson,et al.  Stability of Nonlinear Convection–Diffusion–Reaction Systems in Discontinuous Galerkin Methods , 2017, J. Sci. Comput..

[34]  D. Drikakis,et al.  Large eddy simulation using high-resolution and high-order methods , 2009, Philosophical Transactions of the Royal Society A: Mathematical, Physical and Engineering Sciences.

[35]  C. Hirsch.  Numerical Computation of Internal and External Flows, Volume 2: Computational Methods for Inviscid and Viscous Flows , 1990 .

[36]  Junjie Wei,et al.  Bifurcation and spatiotemporal patterns in a homogeneous diffusive predator-prey system ✩ , 2009 .

[37]  Thilo Gross,et al.  Instabilities in spatially extended predator-prey systems: spatio-temporal patterns in the neighborhood of Turing-Hopf bifurcations. , 2007, Journal of theoretical biology.

[38]  L. Kramer,et al.  Theory of pattern forming systems under traveling-wave forcing , 2007 .

[39]  Tapan K. Sengupta,et al.  Spectral analysis of finite difference schemes for convection diffusion equation , 2017 .

[40]  D. White,et al.  The planforms and onset of convection with a temperature-dependent viscosity , 1988, Journal of Fluid Mechanics.

[41]  T. Poinsot,et al.  Theoretical and numerical combustion , 2001 .

[42]  P. Chu,et al.  A Three-Point Combined Compact Difference Scheme , 1998 .