Numerical simulation of reaction‐diffusion telegraph systems in unbounded domains

We present a boundary element method for computing numerical solutions of the reaction‐diffusion telegraph equation in unbounded domains. This technique does not need artificial boundary conditions at the computational domain and uses a new algorithm to compute the Fourier transform, the convolution theorem, and the fact that the exact solution of the telegraph equation can be written as an integral transform in terms of the fundamental solution. We use the logistic growth model to find how the population of an organism evolves according to its growth rate. © 2014 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 31: 326–335, 2015     

The attractor for a nonlinear reaction‐diffusion system in an unbounded domain

In this paper the quasi‐linear second‐order parabolic systems of reaction‐diffusion type in an unbounded domain are considered. Our aim is to study the long‐time behavior of parabolic systems for which the nonlinearity depends explicitly on the gradient of the unknown functions. To this end we give a systematic study of given parabolic systems and their attractors in weighted Sobolev spaces. Dependence of the Hausdorff dimension of attractors on the weight of the Sobolev spaces is considered. © 2001 John Wiley & Sons, Inc.     

An ETD Crank‐Nicolson method for reaction‐diffusion systems

A novel Exponential Time Differencing Crank‐Nicolson method is developed which is stable, second‐order convergent, and highly efficient. We prove stability and convergence for semilinear parabolic problems with smooth data. In the nonsmooth data case, we employ a positivity‐preserving initial damping scheme to recover the full rate of convergence. Numerical experiments are presented for a wide variety of examples, including chemotaxis and exotic options with transaction cost. © 2011Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2012     

Attractors of derivative complex Ginzburg-Landau equation in unbounded domains

The Ginzburg-Landau-type complex equations are simplified mathematical models for various pattern formation systems in mechanics, physics, and chemistry. In this paper, the derivative complex Ginzburg-Landau (DCGL) equation in an unbounded domain Ω ⊂ ℝ² is studied. We extend the Gagliardo-Nirenberg inequality to the weighted Sobolev spaces introduced by S. V. Zelik. Applied this Gagliardo-Nirenberg inequality of the weighted Sobolev spaces and based on the technique for the semi-linear system of parabolic equations which has been developed by M. A. Efendiev and S. V. Zelik, the global attractor in the corresponding phase space is constructed, the upper bound of its Kolmogorov’s ɛ-entropy is obtained, and the spatial chaos of the attractor for DCGL equation in ℝ² is detailed studied.      

Global existence and nonexistence of generalized coupled reaction‐diffusion systems

In this paper, we consider the initial boundary value problem for a class of reaction‐diffusion systems with generalized coupled source terms. The assumption on the coupled source terms refers to the single equations and includes many kinds of polynomial growth cases. Under this assumption, the reaction‐diffusion systems have a variational structure, which is the foundation of constructing the potential wells to classify the initial data. In subcritical energy level and critical energy level, which are divided from potential well theory, the global existence solution, blow‐up in finite time solution, and asymptotic behavior of solution are obtained, respectively. Furthermore, we show the sufficient conditions of global well posedness with supercritical energy level by combining with comparison principle and semigroup theory.     

Synchronization of delayed coupled reaction‐diffusion systems on networks

In this paper, the exponential synchronization problem of delayed coupled reaction‐diffusion systems on networks (DCRDSNs) is investigated. Based on graph theory, a systematic method is designed to achieve exponential synchronization between two DCRDSNs by constructing a global Lyapunov function for error system. Two different kinds of sufficient synchronization criteria are derived in the form of Lyapunov functions and coefficients of drive‐response systems, respectively. Finally, a numerical example is given to show the usefulness of the proposed criteria. Copyright © 2014 John Wiley & Sons, Ltd.     

Attractors for nonlinear reaction-diffusion systems in unbounded domains via the method of short trajectories

We consider a nonlinear reaction-diffusion equation on the whole space Rd. We prove the well-posedness of the corresponding Cauchy problem in a general functional setting, namely, when the initial datum is uniformly locally bounded in L² only. Then we adapt the short trajectory method to establish the existence of the global attractor and, if d?3, we find an upper bound of its Kolmogorov's ε-entropy.     

Local maxima of solutions to some nonsymmetric reaction‐diffusion systems

In this study, we examine the solution profile of some reaction‐diffusion systems. The systems are derived after approximating the Arrhenius term in our system which models the sintering process and consists of two coupled equations in terms of two unknowns. The unknowns describe the temperature of the solid and the concentration of the fuel. We show the evolution over time of local temperature hot spots. The key idea is to use ‘microscopic scaling’ around the temperature hot spot at the initial time along with asymptotic analysis. We also provide some numerical results that support the efficiency of our analysis. Copyright © 2013 John Wiley & Sons, Ltd.     

The effect of turbulence on mixing in prototype reaction‐diffusion systems

The effect of turbulence on mixing in prototype reaction‐diffusion systems is analyzed here in the special situation where the turbulence is modeled ideally with two separated scales consisting of a large‐scale mean flow plus a small‐scale spatiotemporal periodic flow. In the limit of fast reaction and slow diffusion, it is rigorously proved that the turbulence does not contribute to the location of the mixing zone in the limit and that this mixing zone location is determined solely by advection of the large‐scale velocity field. This surprising result contrasts strongly with earlier work of the authors that always yields a large‐scale propagation speed enhanced by small‐scale turbulence for propagating fronts. The mathematical reasons for these differences are pointed out. This main theorem rigorously justifies the limit equilibrium approximations utilized in non‐premixed turbulent diffusion flames and condensation‐evaporation modeling in cloud physics in the fast reaction limit. The subtle nature of this result is emphasized by explicit examples presented in the fast reaction and zero‐diffusion limit with a nontrivial effect of turbulence on mixing in the limit. The situation with slow reaction and slow diffusion is also studied in the present work. Here the strong stirring by turbulence before significant reaction occurs necessarily leads to a homogenized limit with the strong mixing effects of turbulence expressed by a rigorous turbulent diffusivity modifying the reaction‐diffusion equations. Physical examples from non‐premixed turbulent combustion and cloud microphysics modeling are utilized throughout the paper to motivate and interpret the mathematical results. © 2000 John Wiley & Sons, Inc.     

A fully adaptive approximation for quenching‐type reaction‐diffusion equations over circular domains

This article studies a fully adaptive finite difference method for solving quenching‐type nonlinear reaction‐diffusion equations over circular domains. Although an auxiliary condition at the origin and radial symmetry are imposed, adaptations are accomplished via arc‐length‐based monitoring functions in space and time, respectively. The monotonicity and positivity of the numerical solution are proved following a suitable grid constraint, and the numerical stability is ensured in the von Neumann sense. Theoretical bounds of the critical quenching radius are obtained and then refined through the computation. Computational examples are provided to illustrate the effectiveness and plausibility of the new adaptive computational procedure developed. © 2013 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 30: 472–489, 2014     

Nonlinear diffusion equations as asymptotic limits of Cahn‐Hilliard systems on unbounded domains via Cauchy's criterion

This paper develops an abstract theory for subdifferential operators to give existence and uniqueness of solutions to the initial‐boundary problem P for the nonlinear diffusion equation in an unbounded domain ( ), written as which represents the porous media, the fast diffusion equations, etc, where β is a single‐valued maximal monotone function on , and T>0. In Kurima and Yokota (J Differential Equations 2017; 263:2024‐2050 and Adv Math Sci Appl 2017; 26:221‐242) existence and uniqueness of solutions for P were directly proved under a growth condition for β even though the Stefan problem was excluded from examples of P . This paper completely removes the growth condition for β by confirming Cauchy's criterion for solutions of the following approximate problem _ε with approximate parameter ε>0: which is called the Cahn‐Hilliard system, even if ( ) is an unbounded domain. Moreover, it can be seen that the Stefan problem excluded from Kurima and Yokota (J Differential Equations 2017; 263:2024‐2050 and Adv Math Sci Appl 2017; 26:221‐242) is covered in the framework of this paper.     

Numerical simulation of reaction–diffusion equations on spherical domains

The two-variable reaction diffusion equations on the spherical domain is considered and simulated, using the semi-implicit Euler finite difference method. It is shown that the method keeps the kinetics from overshooting the stable branches when a large time step is used in the simulation.     

An efficient high‐order algorithm for solving systems of reaction‐diffusion equations

An efficient higher‐order finite difference algorithm is presented in this article for solving systems of two‐dimensional reaction‐diffusion equations with nonlinear reaction terms. The method is fourth‐order accurate in both the temporal and spatial dimensions. It requires only a regular five‐point difference stencil similar to that used in the standard second‐order algorithm, such as the Crank‐Nicolson algorithm. The Padé approximation and Richardson extrapolation are used to achieve high‐order accuracy in the spatial and temporal dimensions, respectively. Numerical examples are presented to demonstrate the efficiency and accuracy of the new algorithm. © 2002 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 18: 340–354, 2002; Published online in Wiley InterScience (www.interscience.wiley.com). DOI 10.1002/num.10012     

Two‐grid methods for mixed finite‐element solution of coupled reaction‐diffusion systems

We develop 2‐grid schemes for solving nonlinear reaction‐diffusion systems: where p = (p, q) is an unknown vector‐valued function. The schemes use discretizations based on a mixed finite‐element method. The 2‐grid approach yields iterative procedures for solving the nonlinear discrete equations. The idea is to relegate all the Newton‐like iterations to grids much coarser than the final one, with no loss in order of accuracy. The iterative algorithms examined here extend a method developed earlier for single reaction‐diffusion equations. An application to prepattern formation in mathematical biology illustrates the method's effectiveness. © 1999 John Wiley & Sons, Inc. Numer Methods Partial Differential Eq 15: 589–604, 1999     

Mathematical modeling of time fractional reaction–diffusion systems

We study a fractional reaction–diffusion system with two types of variables: activator and inhibitor. The interactions between components are modeled by cubical nonlinearity. Linearization of the system around the homogeneous state provides information about the stability of the solutions which is quite different from linear stability analysis of the regular system with integer derivatives. It is shown that by combining the fractional derivatives index with the ratio of characteristic times, it is possible to find the marginal value of the index where the oscillatory instability arises. The increase of the value of fractional derivative index leads to the time periodic solutions. The domains of existing periodic solutions for different parameters of the problem are obtained. A computer simulation of the corresponding nonlinear fractional ordinary differential equations is presented. For the fractional reaction–diffusion systems it is established that there exists a set of stable spatio-temporal structures of the one-dimensional system under the Neumann and periodic boundary conditions. The characteristic features of these solutions consist of the transformation of the steady-state dissipative structures to homogeneous oscillations or space temporary structures at a certain value of fractional index and the ratio of characteristic times of system.     

Wave simulations of Gray‐Scott reaction‐diffusion system

In this work, we study the numerical simulation of the one‐dimensional reaction‐diffusion system known as the Gray‐Scott model. This model is responsible for the spatial pattern formation, which we often meet in nature as the result of some chemical reactions. We have used the trigonometric quartic B‐spline (T4B) functions for space discretization with the Crank‐Nicolson method for time integration to integrate the nonlinear reaction‐diffusion equation into a system of algebraic equations. The solutions of the Gray‐Scott model are presented with different wave simulations. Test problems are chosen from the literature to illustrate the stationary waves, pulse‐splitting waves, and self‐replicating waves.     

Traveling wave solutions of n‐dimensional delayed reaction–diffusion systems and application to four‐dimensional predator–prey systems

This paper deals with the existence of traveling wave solutions for n‐dimensional delayed reaction–diffusion systems. By using Schauder's fixed point theorem, we establish the existence result of a traveling wave solution connecting two steady states by constructing a pair of upper–lower solutions that are easy to construct. As an application, we apply our main results to a four‐dimensional delayed predator–prey system and obtain the existence of traveling wave solutions. Copyright © 2016 John Wiley & Sons, Ltd.     

Parameter uniform approximations for time‐dependent reaction‐diffusion problems

Using discrete Green's functions techniques, we present a classification of fitted mesh methods for time‐dependent reaction diffusion problems, based on the analyses of Linß (Linß, Numer Algor 40 (2005), 23–32) for the analogous steady‐state problem and of Kopteva (Kopteva, Computing 66 (2001), 179–197) for time‐dependent convection‐diffusion problems. As examples of how to apply the analysis, we derive error estimates for the fitted meshes of Shishkin and Bakhvalov, and provide supporting numerical results. © 2007 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2007     

Nonstandard finite difference schemes for reaction‐diffusion equations

We construct finite difference schemes for a particular class of one‐space dimension, nonlinear reaction‐diffusion PDEs. The use of nonstandard finite difference methods and the imposition of a positivity condition constrain the schemes to be explicit and allow the determination of functional relations between the space and time step‐sizes. The general procedure is illustrated by applying it to several important model systems of PDEs © 1999 John Wiley & Sons, Inc. Numer Methods Partial Differential Eq 15: 201–214, 1999     

Viscous two‐fluid flows in perturbed unbounded domains

Two stationary plane free boundary value problems for the Navier‐Stokes equations are studied. The first problem models the viscous two‐fluid flow down a perturbed or slightly distorted inclined plane. The second one describes the viscous two‐fluid flow in a perturbed or slightly distorted channel. For sufficiently small data and under certain conditions on parameters the solvability and uniqueness results are proved for both problems. The asymptotic behaviour of the solutions is investigated. For the second problem an example of nonuniqueness is constructed. Computational results of flow problems that are very close to the above problems are presented. (© 2005 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)