skip to main content

Title: Dispersive Hydrodynamics of Soliton Condensates for the Korteweg–de Vries Equation

We consider large-scale dynamics of non-equilibrium dense soliton gas for the Korteweg–de Vries (KdV) equation in the special “condensate” limit. We prove that in this limit the integro-differential kinetic equation for the spectral density of states reduces to theN-phase KdV–Whitham modulation equations derived by Flaschka et al. (Commun Pure Appl Math 33(6):739–784, 1980) and Lax and Levermore (Commun Pure Appl Math 36(5):571–593, 1983). We consider Riemann problems for soliton condensates and construct explicit solutions of the kinetic equation describing generalized rarefaction and dispersive shock waves. We then present numerical results for “diluted” soliton condensates exhibiting rich incoherent behaviors associated with integrable turbulence.

more » « less
Author(s) / Creator(s):
; ; ;
Publisher / Repository:
Springer Science + Business Media
Date Published:
Journal Name:
Journal of Nonlinear Science
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Commuting integral and differential operators connect the topics of signal processing, random matrix theory, and integrable systems. Previously, the construction of such pairs was based on direct calculation and concerned concrete special cases, leaving behind important families such as the operators associated to the rational solutions of the Korteweg–de Vries (KdV) equation. We prove a general theorem that the integral operator associated to every wave function in the infinite-dimensional adelic Grassmannian G r a d of Wilson always reflects a differential operator (in the sense of Definition 1 below). This intrinsic property is shown to follow from the symmetries of Grassmannians of Kadomtsev–Petviashvili (KP) wave functions, where the direct commutativity property holds for operators associated to wave functions fixed by Wilson’s sign involution but is violated in general. Based on this result, we prove a second main theorem that the integral operators in the computation of the singular values of the truncated generalized Laplace transforms associated to all bispectral wave functions of rank 1 reflect a differential operator. A 9 0 ○ rotation argument is used to prove a third main theorem that the integral operators in the computation of the singular values of the truncated generalized Fourier transforms associated to all such KP wave functions commute with a differential operator. These methods produce vast collections of integral operators with prolate-spheroidal properties, including as special cases the integral operators associated to all rational solutions of the KdV and KP hierarchies considered by [Airault, McKean, and Moser, Commun. Pure Appl. Math. 30, 95–148 (1977)] and [Krichever, Funkcional. Anal. i Priložen. 12, 76–78 (1978)], respectively, in the late 1970s. Many examples are presented. 
    more » « less
  2. Abstract

    We study thin films with residual strain by analyzing the$$\Gamma -$$Γ-limit of non-Euclidean elastic energy functionals as the material’s thickness tends to 0. We begin by extending prior results (Bhattacharya et al. in Arch Ration Mech Anal 228: 143–181, 2016); (Agostiniani et al. in ESAIM Control Opt Calculus Var 25: 24, 2019); (Lewicka and Lucic in Commun Pure Appl Math 73: 1880–1932, 2018); (Schmidt in J de Mathématiques Pures et Appliquées 88: 107–122, 2007) , to a wider class of films, whose prestrain depends on both the midplate and the transversal variables. The ansatz for our$$\Gamma -$$Γ-convergence result uses a specific type of wrinkling, which is built on exotic solutions to the Monge-Ampere equation, constructed via convex integration (Lewicka and Pakzad in Anal PDE 10: 695–727, 2017). We show that the expression for our$$\Gamma -$$Γ-limit has a natural interpretation in terms of the orthogonal projection of the residual strain onto a suitable subspace. We also show that some type of wrinkling phenomenon is necessary to match the lower bound of the$$\Gamma -$$Γ-limit in certain circumstances. These results all assume a prestrain of the same order as the thickness; we also discuss why it is natural to focus on that regime by considering what can happen when the prestrain is larger.

    more » « less
  3. Abstract

    We present probabilistic interpretations of solutions to semi-linear parabolic equations with polynomial nonlinearities in terms of the voting models on the genealogical trees of branching Brownian motion (BBM). These extend McKean’s connection between the Fisher–KPP equation and BBM (McKean 1975Commun. Pure Appl. Math.28323–31). In particular, we present ‘random outcome’ and ‘random threshold’ voting models that yield any polynomial nonlinearityfsatisfyingf(0)=f(1)=0and a ‘recursive up the tree’ model that allows to go beyond this restriction onf. We compute several examples of particular interest; for example, we obtain a curious interpretation of the heat equation in terms of a nontrivial voting model and a ‘group-based’ voting rule that leads to a probabilistic view of the pushed-pulled transition for a class of nonlinearities introduced by Ebert and van Saarloos.

    more » « less
  4. null (Ed.)
    Abstract We formulate a class of stochastic partial differential equations based on Kelvin’s circulation theorem for ideal fluids. In these models, the velocity field is randomly transported by white-noise vector fields, as well as by its own average over realizations of this noise. We call these systems the Lagrangian averaged stochastic advection by Lie transport (LA SALT) equations. These equations are nonlinear and non-local, in both physical and probability space. Before taking this average, the equations recover the Stochastic Advection by Lie Transport (SALT) fluid equations introduced by Holm (Proc R Soc A 471(2176):20140963, 2015). Remarkably, the introduction of the non-locality in probability space in the form of momentum transported by its own mean velocity gives rise to a closed equation for the expectation field which comprises Navier–Stokes equations with Lie–Laplacian ‘dissipation’. As such, this form of non-locality provides a regularization mechanism. The formalism we develop is closely connected to the stochastic Weber velocity framework of Constantin and Iyer (Commun Pure Appl Math 61(3):330–345, 2008) in the case when the noise correlates are taken to be the constant basis vectors in $$\mathbb {R}^3$$ R 3 and, thus, the Lie–Laplacian reduces to the usual Laplacian. We extend this class of equations to allow for advected quantities to be present and affect the flow through exchange of kinetic and potential energies. The statistics of the solutions for the LA SALT fluid equations are found to be changing dynamically due to an array of intricate correlations among the physical variables. The statistical properties of the LA SALT physical variables propagate as local evolutionary equations which when spatially integrated become dynamical equations for the variances of the fluctuations. Essentially, the LA SALT theory is a non-equilibrium stochastic linear response theory for fluctuations in SALT fluids with advected quantities. 
    more » « less
  5. Abstract

    Motivated by the challenge of sampling Gibbs measures with nonconvex potentials, we study a continuum birth–death dynamics. We improve results in previous works (Liuet al2023Appl. Math. Optim.8748; Luet al2019 arXiv:1905.09863) and provide weaker hypotheses under which the probability density of the birth–death governed by Kullback–Leibler divergence or byχ2divergence converge exponentially fast to the Gibbs equilibrium measure, with a universal rate that is independent of the potential barrier. To build a practical numerical sampler based on the pure birth–death dynamics, we consider an interacting particle system, which is inspired by the gradient flow structure and the classical Fokker–Planck equation and relies on kernel-based approximations of the measure. Using the technique of Γ-convergence of gradient flows, we show that on the torus, smooth and bounded positive solutions of the kernelised dynamics converge on finite time intervals, to the pure birth–death dynamics as the kernel bandwidth shrinks to zero. Moreover we provide quantitative estimates on the bias of minimisers of the energy corresponding to the kernelised dynamics. Finally we prove the long-time asymptotic results on the convergence of the asymptotic states of the kernelised dynamics towards the Gibbs measure.

    more » « less