A numerical framework for diffusion-controlled bimolecular-reactive systems to enforce maximum principles and the non-negative constraint

We present a novel computational framework for diffusive-reactive systems that satisfies the non-negative constraint and maximum principles on general computational grids. The governing equations for the concentration of reactants and product are written in terms of tensorial diffusion-reaction equations. We restrict our studies to fast irreversible bimolecular reactions. If one assumes that the reaction is diffusion-limited and all chemical species have the same diffusion coefficient, one can employ a linear transformation to rewrite the governing equations in terms of invariants, which are unaffected by the reaction. This results in two uncoupled tensorial diffusion equations in terms of these invariants, which are solved using a novel non-negative solver for tensorial diffusion-type equations. The concentrations of the reactants and the product are then calculated from invariants using algebraic manipulations. The novel aspect of the proposed computational framework is that it will always produce physically meaningful non-negative values for the concentrations of all chemical species. Several representative numerical examples are presented to illustrate the robustness, convergence, and the numerical performance of the proposed computational framework. We will also compare the proposed framework with other popular formulations. In particular, we will show that the Galerkin formulation (which is the standard single-field formulation) does not produce reliable solutions, and the reason can be attributed to the fact that the single-field formulation does not guarantee non-negative solutions. We will also show that the clipping procedure (which produces non-negative solutions but is considered as a variational crime) does not give accurate results when compared with the proposed computational framework.

[1]  Byron Goldstein,et al.  Diffusion Limited Reactions , 2007, SIAM J. Appl. Math..

[2]  Eckehard Schöll,et al.  Nonlinear Spatio-Temporal Dynamics and Chaos in Semiconductors , 2001 .

[3]  R J Barsotti,et al.  Laser photolysis of caged calcium: rates of calcium release by nitrophenyl-EGTA and DM-nitrophen. , 1996, Biophysical journal.

[4]  John P. Crimaldi,et al.  Reaction enhancement of isolated scalars by vortex stirring , 2008 .

[5]  Daniil Svyatskiy,et al.  Interpolation-free monotone finite volume method for diffusion equations on polygonal meshes , 2009, J. Comput. Phys..

[6]  G. M.,et al.  Partial Differential Equations I , 2023, Applied Mathematical Sciences.

[7]  J. Verwer,et al.  Numerical solution of time-dependent advection-diffusion-reaction equations , 2003 .

[8]  J. David Logan,et al.  An Introduction to Nonlinear Partial Differential Equations , 1994 .

[9]  M. Cross,et al.  Pattern formation outside of equilibrium , 1993 .

[10]  P. G. Ciarlet,et al.  Maximum principle and uniform convergence for the finite element method , 1973 .

[11]  Abdul-Majid Wazwaz,et al.  Partial differential equations : methods and applications , 2002 .

[12]  D. Z. Turner,et al.  A stabilized mixed finite element method for Darcy flow based on a multiscale decomposition of the solution , 2006 .

[13]  Noam Agmon,et al.  Diffusion-Limited Acid-Base Nonexponential Dynamics , 2001 .

[14]  Peter K. Kitanidis,et al.  Modeling microbial reactions at the plume fringe subject to transverse mixing in porous media: When can the rates of microbial reaction be assumed to be instantaneous? , 2005 .

[15]  Perry L. McCarty,et al.  Chemical and Biological Processes: The Need for Mixing , 2012 .

[16]  M. Farkas,et al.  Dynamical models in biology , 2001 .

[17]  Nicholas I. M. Gould,et al.  Preprocessing for quadratic programming , 2004, Math. Program..

[18]  R. M. Bowen Part I – Theory of Mixtures , 1976 .

[19]  P. Bedient,et al.  Transport of dissolved hydrocarbons influenced by oxygen-limited biodegradation , 1986 .

[20]  Joel Keizer,et al.  Nonequilibrium statistical thermodynamics and the effect of diffusion on chemical reaction rates , 1982 .

[21]  I. Epstein,et al.  An Introduction to Nonlinear Chemical Dynamics , 1998 .

[22]  Gun-Young Park,et al.  Oxidation of pharmaceuticals during ozonation and advanced oxidation processes. , 2003, Environmental science & technology.

[23]  Zhiming Chen,et al.  A mixed multiscale finite element method for elliptic problems with oscillating coefficients , 2003, Math. Comput..

[24]  Harsha Nagarajan,et al.  Enforcing the non‐negativity constraint and maximum principles for diffusion with decay on general computational grids , 2010, ArXiv.

[25]  Katharina Krischer,et al.  Spatio-Temporal Pattern Formation , 1999 .

[26]  Albert J. Valocchi,et al.  Reply to comments on “Two-dimensional concentration distribution for mixing-controlled bioreactive transport in steady state” by H. Shao et al. , 2009 .

[27]  P E Arratia,et al.  Predicting the progress of diffusively limited chemical reactions in the presence of chaotic advection. , 2006, Physical review letters.

[28]  Peter K. Kitanidis,et al.  Concentration fluctuations and dilution in two-dimensionally periodic heterogeneous porous media , 1996 .

[29]  Stephen P. Boyd,et al.  Convex Optimization , 2004, Algorithms and Theory of Computation Handbook.

[30]  John Thuburn,et al.  A parameterization of mixdown time for atmospheric chemicals , 1997 .

[31]  Jérôme Droniou,et al.  Construction and Convergence Study of Schemes Preserving the Elliptic Local Maximum Principle , 2011, SIAM J. Numer. Anal..

[32]  Robert A. Alberty,et al.  Application of the Theory of Diffusion-controlled Reactions to Enzyme Kinetics , 1958 .

[33]  P. Raviart,et al.  A mixed finite element method for 2-nd order elliptic problems , 1977 .

[34]  Péter Érdi,et al.  Mathematical Models of Chemical Reactions: Theory and Applications of Deterministic and Stochastic Models , 1989 .

[35]  Albert J. Valocchi,et al.  Using dispersivity values to quantify the effects of pore-scale flow focusing on enhanced reaction along a transverse mixing zone , 2010 .

[36]  Eugene A. Kotomin,et al.  Modern aspects of diffusion-controlled reactions : cooperative phenomena in bimolecular processes , 1996 .

[37]  N. Rashevsky,et al.  Mathematical biology , 1961, Connecticut medicine.

[38]  Albert J. Valocchi,et al.  Non-negative mixed finite element formulations for a tensorial diffusion equation , 2008, J. Comput. Phys..

[39]  R. Glasow,et al.  Importance of the surface reaction OH + Cl − on sea salt aerosol for the chemistry of the marine boundary layer – a model study , 2006 .

[40]  Stephen A. Rice Diffusion-limited reactions , 1985 .

[41]  C. Kuo-chen,et al.  Studies on the rate of diffusion-controlled reactions of enzymes. Spatial factor and force field factor. , 1974, Scientia Sinica.

[42]  K. Chou,et al.  Role of the protein outside active site on the diffusion-controlled reaction of enzymes , 1982 .

[43]  Jørg E. Aarnes,et al.  On the Use of a Mixed Multiscale Finite Element Method for GreaterFlexibility and Increased Speed or Improved Accuracy in Reservoir Simulation , 2004, Multiscale Model. Simul..

[44]  M. Fixman,et al.  General theory of diffusion‐controlled reactions , 1973 .

[45]  Gianmarco Manzini,et al.  Analysis of the monotonicity conditions in the mimetic finite difference method for elliptic problems , 2011, J. Comput. Phys..

[46]  Harsha Nagarajan,et al.  A numerical methodology for enforcing maximum principles and the non-negative constraint for transient diffusion equations , 2012, ArXiv.

[47]  J. N. Reddy,et al.  On the performance of high‐order finite elements with respect to maximum principles and the nonnegative constraint for diffusion‐type equations , 2011, ArXiv.

[48]  Albert J. Valocchi,et al.  Two-dimensional concentration distribution for mixing-controlled bioreactive transport in steady state , 2007 .

[49]  Todd Arbogast,et al.  Numerical Subgrid Upscaling of Two-Phase Flow in Porous Media , 2000 .

[50]  Jean E. Roberts,et al.  Mixed and hybrid finite element methods , 1987 .

[51]  Andreas Englert,et al.  Mixing, spreading and reaction in heterogeneous media: a brief review. , 2011, Journal of contaminant hydrology.

[52]  P. Bedient,et al.  Transport of dissolved hydrocarbons influenced by oxygen‐limited biodegradation: 1. Theoretical development , 1986 .

[53]  Robert I. Cukier,et al.  Diffusion-influenced reactions , 1986 .

[54]  Christophe Le Potier,et al.  A nonlinear finite volume scheme satisfying maximum and minimum principles for diffusion operators , 2009 .

[55]  William P. Jencks,et al.  Diffusion-controlled and concerted base catalysis in the decomposition of hemithioacetals , 1969 .

[56]  R. Punnett,et al.  The Genetical Theory of Natural Selection , 1930, Nature.

[57]  Jacek Gondzio,et al.  Multiple centrality corrections in a primal-dual method for linear programming , 1996, Comput. Optim. Appl..

[58]  Daniil Svyatskiy,et al.  Monotone finite volume schemes for diffusion equations on unstructured triangular and shape-regular polygonal meshes , 2007, J. Comput. Phys..

[59]  A. Feinstein,et al.  Variational Methods for the Study of Nonlinear Operators , 1966 .

[60]  John P. Crimaldi,et al.  A proposed mechanism for turbulent enhancement of broadcast spawning efficiency , 2004 .

[61]  Yue-Kin Tsang,et al.  Predicting the evolution of fast chemical reactions in chaotic flows. , 2009, Physical review. E, Statistical, nonlinear, and soft matter physics.

[62]  Christophe Le Potier,et al.  Schéma volumes finis monotone pour des opérateurs de diffusion fortement anisotropes sur des maillages de triangles non structurés , 2005 .

[63]  Graham R. Fleming,et al.  Direct measurement of intrinsic proton transfer rates in diffusion-controlled reactions , 1997 .

[64]  Chia-Ven Pao,et al.  Nonlinear parabolic and elliptic equations , 1993 .

[65]  S. Benson,et al.  The Kinetics of Free Radical Polymerization under Conditions of Diffusion-controlled Termination , 1962 .

[66]  Knut-Andreas Lie,et al.  Mixed multiscale finite elements and streamline methods for reservoir simulation of large geomodels , 2005 .

[67]  A. Concha Theory of Mixtures , 2014 .

[68]  Charles J Werth,et al.  Evaluation of the effects of porous media structure on mixing-controlled reactions using pore-scale modeling and micromodel experiments. , 2008, Environmental science & technology.

[69]  Sanjay Mehrotra,et al.  On the Implementation of a Primal-Dual Interior Point Method , 1992, SIAM J. Optim..

[70]  Peter K. Kitanidis,et al.  Macrodispersion of sorbing solutes in heterogeneous porous formations with spatially periodic retardation factor and velocity field , 1992 .

[71]  D K Smith,et al.  Numerical Optimization , 2001, J. Oper. Res. Soc..

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

[73]  Thomas J. R. Hughes,et al.  A stabilized mixed finite element method for Darcy flow , 2002 .

[74]  P. Bassanini,et al.  Elliptic Partial Differential Equations of Second Order , 1997 .

[75]  K. Hjelmstad Fundamentals of Structural Mechanics , 1996 .

[76]  A. North,et al.  Diffusion-controlled reactions , 1966 .

[77]  Henning Prommer,et al.  Effects of hydrodynamic dispersion on plume lengths for instantaneous bimolecular reactions , 2004 .