Abstract: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.

Group classification of systems of non-linear reaction-diffusion equations with general diffusion matrix. III. Triangular diffusion matrix

Convergence to equilibrium for a degenerate triangular reaction-diffusion system

Existence of solution of a triangular degenerate reaction-diffusion system

Convergence to equilibrium for a degenerate three species reaction-diffusion system

Designing reaction-cross-diffusion systems with Turing and wave instabilities

Global wellposedness for a class of reaction-advection-anisotropic-diffusion systems

A general class of invariant diffusion processes in one dimension

Asymptotics of solutions in nA+nB->C reaction Diffusion systems

Numerical studies of a class of reaction–diffusion equations with Stefan conditions

Global Existence to a Class of Triangular Block Matrix Cross Diffusion Systems and the Spectral Gap Condition

A class of fractional parabolic reaction-diffusion systems with control of total mass: theory and numerics

Stationary states of reaction-diffusion and Schrödinger systems with inhomogeneous or controlled diffusion

A class of fractional parabolic reaction–diffusion systems with control of total mass: theory and numerics

Complex spatiotemporal phenomena in multi-dimensional reaction–diffusion system

Symmetry and Numerical Solutions for Systems of Non-linear Reaction Diffusion Equations

New conditional symmetries and exact solutions of nonlinear reaction-diffusion-convection equations. I

A monotone finite element method for reaction-drift-diffusion equations with discontinuous reaction coefficients

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

The generalized method of separation of variables for diffusion-influenced reactions: Irreducible Cartesian tensors technique

Study of a class of triangular starvation driven cross-diffusion systems