No Arabic abstract
Although the spatially continuous version of the reaction-diffusion equation has been well studied, in some instances a spatially-discretized representation provides a more realistic approximation of biological processes. Indeed, mathematically the discretized and continuous systems can lead to different predictions of biological dynamics. It is well known in the continuous case that the incorporation of diffusion can cause diffusion-driven blow-up with respect to the $L^{infty}$ norm. However, this does not imply diffusion-driven blow-up will occur in the discretized version of the system. For example, in a continuous reaction-diffusion system with Dirichlet boundary conditions and nonnegative solutions, diffusion-driven blow up occurs even when the total species concentration is non-increasing. For systems that instead have homogeneous Neumann boundary conditions, it is currently unknown whether this deviation between the continuous and discretized system can occur. Therefore, it is worth examining the discretized system independently of the continuous system. Since no criteria exist for the boundedness of the discretized system, the focus of this paper is to determine sufficient conditions to guarantee the system with diffusion remains bounded for all time. We consider reaction-diffusion systems on a 1D domain with homogeneous Neumann boundary conditions and non-negative initial data and solutions. We define a Lyapunov-like function and show that its existence guarantees that the discretized reaction-diffusion system is bounded. These results are considered in the context of three example systems for which Lyapunov-like functions can and cannot be found.
Similarity solutions play an important role in many fields of science: we consider here similarity in stochastic dynamics. Important issues are not only the existence of stochastic similarity, but also whether a similarity solution is dynamically attractive, and if it is, to what particular solution does the system evolve. By recasting a class of stochastic PDEs in a form to which stochastic centre manifold theory may be applied we resolve these issues in this class. For definiteness, a first example of self-similarity of the Burgers equation driven by some stochastic forced is studied. Under suitable assumptions, a stationary solution is constructed which yields the existence of a stochastic self-similar solution for the stochastic Burgers equation. Furthermore, the asymptotic convergence to the self-similar solution is proved. Second, in more general stochastic reaction-diffusion systems stochastic centre manifold theory provides a framework to construct the similarity solution, confirm its relevance, and determines the correct solution for any compact initial condition. Third, we argue that dynamically moving the spatial origin and dynamically stretching time improves the description of the stochastic similarity. Lastly, an application to an extremely simple model of turbulent mixing shows how anomalous fluctuations may arise in eddy diffusivities. The techniques and results we discuss should be applicable to a wide range of stochastic similarity problems.
Motivated by the normal form of a fast-slow ordinary differential equation exhibiting a pitchfork singularity we consider the discrete-time dynamical system that is obtained by an application of the explicit Euler method. Tracking trajectories in the vicinity of the singularity we show, how the slow manifold extends beyond the singular point and give an estimate on the contraction rate of a transition mapping. The proof relies on the blow-up method suitably adapted to the discrete setting where a key technical contribution are precise estimates for a cubic map in the central rescaling chart.
We extend slow manifolds near a transcritical singularity in a fast-slow system given by the explicit Euler discretization of the corresponding continuous-time normal form. The analysis uses the blow-up method and direct trajectory-based estimates. We prove that the qualitative behaviour is preserved by a time-discretization with sufficiently small step size. This step size is fully quantified relative to the time scale separation. Our proof also yields the continuous-time results as a special case and provides more detailed calculations in the classical (or scaling) chart.
In this paper we are concerned with the regularity of solutions to a nonlinear elliptic system of $m$ equations in divergence form, satisfying $p$ growth from below and $q$ growth from above, with $p leq q$; this case is known as $p, q$-growth conditions. Well known counterexamples, even in the simpler case $p=q$, show that solutions to systems may be singular; so, it is necessary to add suitable structure conditions on the system that force solutions to be regular. Here we obtain local boundedness of solutions under a componentwise coercivity condition. Our result is obtained by proving that each component $u^alpha$ of the solution $u=(u^1,...,u^m)$ satisfies an improved Caccioppolis inequality and we get the boundedness of $u^{alpha}$ by applying De Giorgis iteration method, provided the two exponents $p$ and $q$ are not too far apart. Let us remark that, in dimension $n=3$ and when $p=q$, our result works for $frac{3}{2} < p < 3$, thus it complements the one of Bjorn whose technique allowed her to deal with $p leq 2$ only. In the final section, we provide applications of our result.
In this paper we prove the time-domain boundedness for noise-to-state exponentially stable systems, and further make an estimation of its lower bound function, which allows to answer the question that how long the solution of a stochastic noise-to-state exponentially stable system stays in the domain of attraction and what happens with it if it escapes from this region for a while. The results will complement the probability-domain boundedness of noise-to-state exponentially stable systems, and provide a new insight into noise-to-state exponential stability.