Long-time Behavior Research Articles

Long-time behavior of stochastic Hamilton-Jacobi equations

The long-time behavior of stochastic Hamilton-Jacobi equations is analyzed, including the stochastic mean curvature flow as a special case. In a variety of settings, new and sharpened results are obtained. Among them are (i) a regularization by noise phenomenon for the mean curvature flow with homogeneous noise which establishes that the inclusion of noise speeds up the decay of solutions, and (ii) the long-time convergence of solutions to spatially inhomogeneous stochastic Hamilton-Jacobi equations. A number of motivating examples about nonlinear stochastic partial differential equations are presented in the appendix.

Global well posedness for the semilinear edge-degenerate parabolic equations on singular manifolds

Abstract In this article, we study the long-time dynamical behavior of the solution for a class of semilinear edge-degenerate parabolic equations on manifolds with edge singularities. By introducing a family of potential well and compactness method, we reveal the dependence between the initial data and the long-time dynamical behavior of the solution. Specifically, we give the threshold condition for the initial data, which makes the solution exist globally or blowup in finite-time with subcritical, critical, and supercritical initial energy, respectively. Moreover, we also discussed the long-time behavior of the global solution, the estimate of blowup time, and blowup rate. Our results show that the relationship between the initial data and the long-time behavior of the solution can be revealed in the weighted Sobolev spaces for nonlinear parabolic equations on manifolds with edge singularities.

On long-time behavior for solutions of the Gear–Grimshaw system

On long-time behavior for solutions of the Gear–Grimshaw system

Improved uniform error estimates for the two-dimensional nonlinear space fractional Dirac equation with small potentials over long-time dynamics

We develop improved uniform error bounds on a second-order Strang splitting method for the long-time dynamics of the nonlinear space fractional Dirac equation (NSFDE) in two dimension (2D) with small electromagnetic potentials. First, a Strang splitting approach is implemented to discretize NSFDE in time. Afterwards the Fourier pseudospectral method is used to complete the discretization of NSFDE in space. With the aid of a second-order Strang splitting approach employed to the Dirac equation, the major local truncation error of the indicated numerical methods is established. Moreover, for the semi-discrete scheme and full-discretization, we rigorously demonstrate the improved, sharp uniform error estimates are O(ετ2) and O(h1m+h2m+ετ2) in virtue of the regularity compensation oscillation (RCO) technique. In the formulations, τ is the time step, hi(i=1,2) stands for spatial sizes in xi-directions, m is dependent on the regularity of solutions, and ε∈(0,1]. In order to verify our error bounds and to illustrate some fascinating long-time dynamical behaviors of the NSFDE with honeycomb lattice potentials for varied ε, numerical investigations are presented.

Dynamics of a reaction–diffusion waterborne pathogen model with free boundaries

In this paper, we consider a reaction–diffusion waterborne pathogen model with free boundary. We first prove that this problem has a unique global solution, and then we show that its longtime behavior is determined by a spreading–vanishing dichotomy. Moreover, we obtain sharp criteria for spreading and vanishing, which is very different from the results in Zhou et al. (2018) that the epidemic will spread for the basic reproduction number is larger than 1.

Spread of parasites affecting death and division rates in a cell population

We introduce a general class of branching Markov processes for the modelling of a parasite infection in a cell population. Each cell contains a quantity of parasites which evolves as a diffusion with positive jumps. The drift, diffusive function and positive jump rate of this quantity of parasites depend on its current value. The division rate of the cells also depends on the quantity of parasites they contain. At division, a cell gives birth to two daughter cells and shares its parasites between them. Cells may also die, at a rate which may depend on the quantity of parasites they contain. We study the long-time behaviour of the parasite infection.

Single-file diffusion in spatially inhomogeneous systems.

We study the effect of spatially varying potential and diffusivity on the dispersion of a tracer particle in single-file diffusion. Noninteracting particles in such a system exhibit normal diffusion at late times, which is characterized by an effective diffusion constant D_{eff}. Here we demonstrate the physically appealing result that the dispersion of single-file tracers in this system has the same long-time behavior as that for Brownian particles in a spatially homogeneous system with constant diffusivity D_{eff}. Our results are based on a late-time analysis of the Fokker-Planck equation, motivated by the mathematical theory of homogenization. The findings are confirmed by numerical simulations for both annealed and quenched initial conditions.

Coagulation equations with source leading to anomalous self-similarity

We study the long-time behaviour of the solutions to Smoluchowski coagulation equations with a source term of small clusters. The source drives the system out-of-equilibrium, leading to a rich range of different possible long-time behaviours, including anomalous self-similarity. The coagulation kernel is non-gelling, homogeneous, with homogeneity , and behaves like when with . Our analysis shows that the long-time behaviour of the solutions depends on the parameters γ and λ. More precisely, we argue that the long-time behaviour is self-similar, although the scaling of the self-similar solutions depends on the sign of and on whether or . In all these cases, the scaling differs from the usual one that has been previously obtained when or . In the last part of the paper, we present some conjectures supporting the self-similar ansatz also for the critical case .

Longtime behavior of the weakly coupled Euler-Bernoulli plate system with structural damping

In this paper, we study the longtime behavior of the weakly coupled Euler-Bernoulli plate system with one structural damping. The system consists of two plate equations coupled by zero order terms. First, we show that when the coupling coefficient function is a nonzero constant and the damping coefficient function is a positive constant, the optimal energy decay rate of the system is t−23. Then we consider the case that the coupling and damping coefficient functions have compact supports. By Carleman estimate and the frequency domain method, we show that the energy of the system decays logarithmically when the supports of the coupling and damping coefficient functions satisfy suitable assumptions.

Entropy dissipative higher order accurate positivity preserving time-implicit discretizations for nonlinear degenerate parabolic equations

We develop entropy dissipative higher order accurate local discontinuous Galerkin (LDG) discretizations coupled with Diagonally Implicit Runge–Kutta (DIRK) methods for nonlinear degenerate parabolic equations with a gradient flow structure. Using the simple alternating numerical flux, we construct DIRK-LDG discretizations that combine the advantages of higher order accuracy, entropy dissipation and proper long-time behavior. We theoretically prove the entropy dissipation of the implicit Euler-LDG discretization without any time-step restrictions when no positivity constraint is imposed. Next, in order to ensure the positivity of the numerical solution, we use the Karush–Kuhn–Tucker (KKT) limiter, which achieves a positive solution by coupling the positivity preserving KKT conditions with higher order accurate DIRK-LDG discretizations using Lagrange multipliers. In addition, mass conservation of the positivity-limited solution is ensured by imposing a mass conservation equality constraint to the KKT equations. Under a time step restriction, the unique solvability and entropy dissipation for implicit first order accurate in time, but higher order accurate in space, positivity-preserving LDG discretizations with periodic boundary conditions are proved, which provide a first theoretical analysis of the KKT limiter. Finally, numerical results demonstrate the higher order accuracy and entropy dissipation of the positivity-preserving DIRK-LDG discretizations for problems requiring a positivity limiter. In addition, we can observe from the numerical results that the implicit time-discrete methods alleviate the time-step restrictions needed for the stability of the numerical discretizations, which improves computational efficiency.

Efficient numerical description of the dynamics of interacting multispecies quantum gases

We present a highly efficient method for the numerical solution of coupled Gross–Pitaevskii equations describing the evolution dynamics of a multi-species mixture of Bose–Einstein condensates in time-dependent potentials. This method, based on a moving and expanding reference frame, compares favorably to a more standard but much more computationally expensive solution based on a frozen frame. It allows an accurate description of the long-time behavior of interacting, multi-species quantum mixtures including the challenging problem of long free expansions relevant to microgravity and space experiments. We demonstrate a successful comparison to experimental measurements of a binary Rb–K mixture recently performed with the payload of a sounding rocket experiment.

Large time behavior of signed fractional porous media equations on bounded domains

Following the methodology of Brasco (Adv Math 394:108029, 2022), we study the long-time behavior for the signed fractional porous medium equation in open bounded sets with smooth boundary. Homogeneous exterior Dirichlet boundary conditions are considered. We prove that if the initial datum has sufficiently small energy, then the solution, once suitably rescaled, converges to a nontrivial constant sign solution of a sublinear fractional Lane–Emden equation. Furthermore, we give a nonlocal sufficient energetic criterion on the initial datum, which is important to identify the exact limit profile, namely the positive solution or the negative one.

On the dynamics of a Leslie–Gower predator–prey ternary model with intraguild

AbstractIn this paper, a predator–prey model with intraguild predation describing the evolution between three interacting species—namely prey, mesopredator and top predator—is investigated, with the aim to model a complete food web. In particular, the longtime behaviour of the solutions is analysed, proving the existence of an absorbing set, and the linear and nonlinear stability analyses of the coexistence equilibrium are performed.

Existence and Internal Structure of the Deterministic Attracting Set for a Random Ant Colonies Model

This paper is concerned with the attracting set of an ant colonies model with bounded noisy perturbation. This perturbation is modeled by the well-known Ornstein–Uhlenbeck process and the arc tangent function. For the random model, we first verify the existence and uniqueness of the global positive solution, and then prove the existence of the deterministic attracting set. Furthermore, in order to reveal more detailed information about the long-time behavior of the ant colonies system, we analyze the internal structure of the attracting set and provide some conditions under which coexistence (or extinction) of the ant species exists in the ant colonies system.

Confined run-and-tumble model with boundary aggregation: Long-time behavior and convergence to the confined Fokker–Planck model

The motile micro-organisms such as Escherichia coli, sperm, or some seaweed are usually modeled by self-propelled particles that move with the run-and-tumble process. Individual-based stochastic models are usually employed to model the aggregation phenomenon at the boundary, which is an active research field that has attracted a lot of biologists and biophysicists. Self-propelled particles at the microscale have complex behaviors, while characteristics at the population level are more important for practical applications but rely on individual behaviors. Kinetic PDE models that describe the time evolution of the probability density distribution of the motile micro-organisms are widely used. However, how to impose the appropriate boundary conditions that take into account the boundary aggregation phenomena is rarely studied. In this paper, we propose the boundary conditions for a 2D confined run-and-tumble model (CRTM) for self-propelled particle populations moving between two parallel plates with a run-and-tumble process. The proposed model satisfies the relative entropy inequality and thus long-time convergence. We establish the relation between CRTM and the confined Fokker–Planck model (CFPM) studied in [J. Fu, B. Perthame and M. Tang, Fokker–Plank system for movement of micro-organism population in confined environment, J. Statist. Phys. 184 (2021) 1–25]. We prove theoretically that when the tumble is highly forward peaked and frequent enough, CRTM converges asymptotically to the CFPM. A numerical comparison of the CRTM with aggregation and CFPM is given. The time evolution of both the deterministic PDE model and individual-based stochastic simulations are displayed, which match each other well.

A multiple time scale approach for anisotropic inertial wave turbulence

Wave turbulence is the study of the long-time statistical behaviour of equations describing a set of weakly nonlinear interacting waves. Such a theory, which has a natural asymptotic closure, allows us to probe the nature of turbulence more deeply than the exact Kolmogorov laws by rigorously proving the direction of the cascade and the existence of an inertial range, predicting stationary spectra for conserved quantities, or evaluating the Kolmogorov constant. An emblematic example is given by fast rotating fluids for which a wave turbulence theory has been derived by Galtier (Phys. Rev. E, vol. 68, issue 1, 2003, p. 015301). This work involves non-trivial analytical developments for a problem that is anisotropic by nature. We propose here a new path for the derivation of the kinetic equation by using the anisotropy at the beginning of the analysis. We show that the helicity basis is not necessary to obtain the wave amplitude equation for the canonical variables that involve a combination of poloidal and toroidal fields. The multiple time scale method adapted to this anisotropic problem is then used to derive the kinetic equation that is the same as the original work when anisotropy is eventually taken into account. This result proves the commutativity between asymptotic closure and anisotropy. In addition, the multiple time scale method informs us that the kinetic equation can be derived without imposing restrictions on the probability distribution of the wave amplitude such as quasi-Gaussianity, or on the phase such as random phase approximation that naturally occurs dynamically.

Regularity theory for a new class of fractional parabolic stochastic evolution equations

A new class of fractional-order parabolic stochastic evolution equations of the form (∂t+A)γX(t)=W˙Q(t)\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$(\\partial _t + A)^\\gamma X(t) = {\\dot{W}}^Q(t)$$\\end{document}, t∈[0,T]\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$t\\in [0,T]$$\\end{document}, γ∈(0,∞)\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\gamma \\in (0,\\infty )$$\\end{document}, is introduced, where -A\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$-A$$\\end{document} generates a C0\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$C_0$$\\end{document}-semigroup on a separable Hilbert space H and the spatiotemporal driving noise W˙Q\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$${\\dot{W}}^Q$$\\end{document} is the formal time derivative of an H-valued cylindrical Q-Wiener process. Mild and weak solutions are defined; these concepts are shown to be equivalent and to lead to well-posed problems. Temporal and spatial regularity of the solution process X are investigated, the former being measured by mean-square or pathwise smoothness and the latter by using domains of fractional powers of A. In addition, the covariance of X and its long-time behavior are analyzed. These abstract results are applied to the cases when A:=Lβ\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$A:= L^\\beta $$\\end{document} and Q:=L~-α\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$Q:={\\widetilde{L}}^{-\\alpha }$$\\end{document} are fractional powers of symmetric, strongly elliptic second-order differential operators defined on (i) bounded Euclidean domains or (ii) smooth, compact surfaces. In these cases, the Gaussian solution processes can be seen as generalizations of merely spatial (Whittle–)Matérn fields to space–time.

Krylov complexity and spectral form factor for noisy random matrix models

We study the spectral properties of two classes of random matrix models: non-Gaussian RMT with quartic and sextic potentials, and RMT with Gaussian noise. We compute and analyze the quantum Krylov complexity and the spectral form factor for both of these models. We find that both models show suppression of the spectral form factor at short times due to decoherence effects, but they differ in their long-time behavior. In particular, we show that the Krylov complexity for the non-Gaussian RMT and RMT with noise deviates from that of a Gaussian RMT. We discuss the implications and limitations of our results for quantum chaos and quantum information in open quantum systems. Our study reveals the distinct sensitivities of the spectral form factor and complexity to non-Gaussianity and noise, which contribute to the observed differences in the different time domains.

Pyramidal traveling waves around an obstacle

This paper is considered with a three-dimensional pyramidal traveling wave of the homogeneous reaction–diffusion equation with a compact obstacle in R3. The long-time behavior of the entire solution derived from a pyramidal traveling wave after passing an obstacle is explored. Mainly relying on auxiliary functions, some new forms of super- and sub-solutions with good characteristics are proposed, and some properties of the pyramidal traveling wave are obtained. The results show that when the propagation of the entire solution is complete, the solution will recover to the same pyramidal traveling wave as before.

The Navier–Stokes equations with body forces decaying coherently in time

The long-time behavior of solutions of the three-dimensional Navier–Stokes equations in a periodic domain is studied. The time-dependent body force decays, as time t tends to infinity, in a coherent manner. In fact, it is assumed to have a general and complicated asymptotic expansion which involves complex powers of et, t, ln⁡t, or other iterated logarithmic functions of t. We prove that all Leray–Hopf weak solutions admit an asymptotic expansion which is independent of the solutions and is uniquely determined by the asymptotic expansion of the body force. The proof makes use of the complexifications of the Gevrey–Sobolev spaces together with those of the Stokes operator and the bilinear form of the Navier–Stokes equations.