skip to main content

Title: Global solutions of aggregation equations and other flows with random diffusion

Aggregation equations, such as the parabolic-elliptic Patlak–Keller–Segel model, are known to have an optimal threshold for global existence versus finite-time blow-up. In particular, if the diffusion is absent, then all smooth solutions with finite second moment can exist only locally in time. Nevertheless, one can ask whether global existence can be restored by adding a suitable noise to the equation, so that the dynamics are now stochastic. Inspired by the work of Buckmaster et al. (Int Math Res Not IMRN 23:9370–9385, 2020) showing that, with high probability, the inviscid SQG equation with random diffusion has global classical solutions, we investigate whether suitable random diffusion can restore global existence for a large class of active scalar equations in arbitrary dimension with possibly singular velocity fields. This class includes Hamiltonian flows, such as the SQG equation and its generalizations, and gradient flows, such as those arising in aggregation models. For this class, we show global existence of solutions in Gevrey-type Fourier–Lebesgue spaces with quantifiable high probability.

more » « less
Author(s) / Creator(s):
Publisher / Repository:
Springer Science + Business Media
Date Published:
Journal Name:
Probability Theory and Related Fields
Page Range / eLocation ID:
p. 1219-1262
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract

    We consider the focusing energy-critical quintic nonlinear wave equation in 3D Euclidean space. It is known that this equation admits a one-parameter family of radial stationary solutions, called solitons, which can be viewed as a curve in $ \dot H^s_x({{\mathbb{R}}}^3) \times H^{s-1}_x({{\mathbb{R}}}^3)$, for any $s> 1/2$. By randomizing radial initial data in $ \dot H^s_x({{\mathbb{R}}}^3) \times H^{s-1}_x({{\mathbb{R}}}^3)$ for $s> 5/6$, which also satisfy a certain weighted Sobolev condition, we produce with high probability a family of radial perturbations of the soliton that give rise to global forward-in-time solutions of the focusing nonlinear wave equation that scatter after subtracting a dynamically modulated soliton. Our proof relies on a new randomization procedure using distorted Fourier projections associated to the linearized operator around a fixed soliton. To our knowledge, this is the 1st long-time random data existence result for a focusing wave or dispersive equation on Euclidean space outside the small data regime.

    more » « less
  2. Abstract

    We give sharp conditions for the large time asymptotic simplification of aggregation-diffusion equations with linear diffusion. As soon as the interaction potential is bounded and its first and second derivatives decay fast enough at infinity, then the linear diffusion overcomes its effect, either attractive or repulsive, for large times independently of the initial data, and solutions behave like the fundamental solution of the heat equation with some rate. The potential$$W(x) \sim \log |x|$$W(x)log|x|for$$|x| \gg 1$$|x|1appears as the natural limiting case when the intermediate asymptotics change. In order to obtain such a result, we produce uniform-in-time estimates in a suitable rescaled change of variables for the entropy, the second moment, Sobolev norms and the$$C^\alpha $$Cαregularity with a novel approach for this family of equations using modulus of continuity techniques.

    more » « less
  3. Simulation of flow and transport in petroleum reservoirs involves solving coupled systems of advection-diffusion-reaction equations with nonlinear flux functions, diffusion coefficients, and reactions/wells. It is important to develop numerical schemes that can approximate all three processes at once, and to high order, so that the physics can be well resolved. In this paper, we propose an approach based on high order, finite volume, implicit, Weighted Essentially NonOscillatory (iWENO) schemes. The resulting schemes are locally mass conservative and, being implicit, suited to systems of advection-diffusion-reaction equations. Moreover, our approach gives unconditionally L-stable schemes for smooth solutions to the linear advection-diffusion-reaction equation in the sense of a von Neumann stability analysis. To illustrate our approach, we develop a third order iWENO scheme for the saturation equation of two-phase flow in porous media in two space dimensions. The keys to high order accuracy are to use WENO reconstruction in space (which handles shocks and steep fronts) combined with a two-stage Radau-IIA Runge-Kutta time integrator. The saturation is approximated by its averages over the mesh elements at the current time level and at two future time levels; therefore, the scheme uses two unknowns per grid block per variable, independent of the spatial dimension. This makes the scheme fairly computationally efficient, both because reconstructions make use of local information that can fit in cache memory, and because the global system has about as small a number of degrees of freedom as possible. The scheme is relatively simple to implement, high order accurate, maintains local mass conservation, applies to general computational meshes, and appears to be robust. Preliminary computational tests show the potential of the scheme to handle advection-diffusion-reaction processes on meshes of quadrilateral gridblocks, and to do so to high order accuracy using relatively long time steps. The new scheme can be viewed as a generalization of standard cell-centered finite volume (or finite difference) methods. It achieves high order in both space and time, and it incorporates WENO slope limiting. 
    more » « less
  4. Abstract

    A physically based subgrid variability model for snow process using the Fokker‐Planck equation (FPE) approach was proposed. This FPE can express the evolution of the probability density function (PDF) of snow depth within a finite area, possibly a grid cell of distributed models or a small basin, whose shape can be irregular. The main advantage of this approach is that it does not rely on a given PDF but dynamically computes the PDF through an advection‐diffusion‐type equation, the FPE, which was derived from point‐scale process‐based governing equations. Snow depth was treated as a random variable, while the snow redistribution and snowmelt rate were treated as the sources of stochasticity. The main challenge in solving this FPE is evaluating the time‐space covariances appearing in the diffusion coefficient. In this study, approximations to evaluate the covariance terms, accounting for snowmelt and snow redistribution, were proposed. The simulated results of the FPE model were validated by the measured time series of snow depth at one site and the spatial distributions of snow depth measured by ground penetrating radar and airborne light detection and ranging (Lidar). It was shown that the point‐observed snow depth fell within the simulated range during most of the 2‐year study period. The simulated PDFs of snow depth within the study area were similar to the observed PDFs of snow depth by ground penetrating radar and Lidar. In summary, these results demonstrate the efficacy of the proposed FPE model representing the subgrid variability of snow depth.

    more » « less
  5. Abstract

    We consider a higher‐dimensional version of the Benjamin‐Ono (HBO) equation in the 2D setting:, which is‐critical, and investigate properties of solutions both analytically and numerically. For a generalized equation (fractional 2D gKdV) after deriving the Pohozaev identities, we obtain nonexistence conditions for solitary wave solutions, then prove uniform bounds in the energy space or conditional global existence, and investigate the radiation region, a specific wedge in the negative‐direction. We then introduce our numerical approach in a general context, and apply it to obtain the ground state solution in the 2D critical HBO equation, then show that its mass is a threshold for global versus finite time existing solutions, which is typical in the focusing (mass‐)critical dispersive equations. We also observe that globally existing solutions tend to disperse completely into the radiation in this nonlocal equation. The blow‐up solutions travel in the positive‐direction with the rescaled ground state profile while also radiating dispersive oscillations into the radiative wedge. We conclude with examples of different interactions of two solitary wave solutions, including weak and strong interactions.

    more » « less