Numerical Methods for Stiff Reaction-diffusion Systems
نویسندگان
چکیده
In a previous study [21], a class of efficient semi-implicit schemes was developed for stiff reaction-diffusion systems. This method which treats linear diffusion terms exactly and nonlinear reaction terms implicitly has excellent stability properties, and its second-order version, with a name IIF2, is linearly unconditionally stable. In this paper, we present another linearly unconditionally stable method that approximates both diffusions and reactions implicitly using a second order Crank-Nicholson scheme. The nonlinear system resulted from the implicit approximation at each time step is solved using a multi-grid method. We compare this method (CN-MG) with IIF2 for their accuracy and efficiency. Numerical simulations demonstrate that both methods are accurate and robust with convergence using even very large size of time steps. IIF2 is found to be more accurate for systems with large diffusion while CN-MG is more efficient when the number of spatial grid points is large.
منابع مشابه
A hybrid method with optimal stability properties for the numerical solution of stiff differential systems
In this paper, we consider the construction of a new class of numerical methods based on the backward differentiation formulas (BDFs) that be equipped by including two off--step points. We represent these methods from general linear methods (GLMs) point of view which provides an easy process to improve their stability properties and implementation in a variable stepsize mode. These superioritie...
متن کاملA numerical treatment of a reaction-diffusion model of spatial pattern in the embryo
In this work the mathematical model of a spatial pattern in chemical and biological systems is investigated numerically. The proposed model considered as a nonlinear reaction-diffusion equation. A computational approach based on finite difference and RBF-collocation methods is conducted to solve the equation with respect to the appropriate initial and boundary conditions. The ability and robust...
متن کاملOperator splitting implicit integration factor methods for stiff reaction–diffusion–advection systems
For reaction–diffusion–advection equations, the stiffness from the reaction and diffusion terms often requires very restricted time step size, while the nonlinear advection term may lead to a sharp gradient in localized spatial regions. It is challenging to design numerical methods that can efficiently handle both difficulties. For reaction–diffusion systems with both stiff reaction and diffusi...
متن کاملOperator splitting implicit integration factor methods for stiff reaction-diffusion-advection systems
For reaction-diffusion-advection equations, the stiffness from the reaction and diffusion terms often requires very restricted time step size, while the nonlinear advection term may lead to a sharp gradient in localized spatial regions. It is challenging to design numerical methods that can efficiently handle both difficulties. For reaction-diffusion systems with both stiff reaction and diffusi...
متن کاملEfficient semi-implicit schemes for stiff systems
When explicit time discretization schemes are applied to stiff reaction–diffusion equations, the stability constraint on the time step depends on two terms: the diffusion and the reaction. The part of the stability constraint due to diffusion can be totally removed if the linear diffusions are treated exactly using integration factor (IF) or exponential time differencing (ETD) methods. For syst...
متن کاملPositivity-preserving nonstandard finite difference Schemes for simulation of advection-diffusion reaction equations
Systems in which reaction terms are coupled to diffusion and advection transports arise in a wide range of chemical engineering applications, physics, biology and environmental. In these cases, the components of the unknown can denote concentrations or population sizes which represent quantities and they need to remain positive. Classical finite difference schemes may produce numerical drawback...
متن کامل