Singularities and the distribution of density in the Burgers/adhesion model

Abstract

Physica D in pressWe are interested in the tail behavior of the pdf of mass density within the one and d-dimensional Burgers/adhesion model used, e.g., to model the formation of large-scale structures in the Universe after baryon-photon decoupling. We show that large densities are localized near "kurtoparabolic" singularities residing on space-time manifolds of codimension two (d ≤ 2) or higher (d ≥ 3). For smooth initial conditions, such singularities are obtained from the convex hull of the Lagrangian potential (the initial velocity potential minus a parabolic term). The singularities contribute universal power-law tails to the density pdf when the initial conditions are random. In one dimension the singularities are preshocks (nascent shocks), whereas in two and three dimensions they persist in time and correspond to boundaries of shocks; in all cases the corresponding density pdf has the exponent −7/2, originally proposed by E, Sinai (1997 Phys. Rev. Lett. 78, 1904) for the pdf of velocity gradients in one-dimensional forced Burgers turbulence. We also briefly consider models permitting particle crossings and thus multi-stream solutions, such as the Zel'dovich approximation and the (Jeans)-Vlasov-Poisson equation with single-stream initial data: they have singularities of codimension one, yielding power-law tails with exponent −3.su quell' immenso baratro di stelle sopra quei gruppi, sopra quegli ammassi, quel seminío, quel balenío di stelle Giovanni Pascoli, from La Vertigine [1]

New sequential exact Euclidean distance transform algorithms based on convex analysis

“…Historically, its computation was first motivated by the study of Hamilton-Jacobi equations, and gave birth to several numerical algorithms [28][29][30][31][32], and later to a linear-time algorithm [33], with applications in numerous fields: numerical simulation of Burger's equation (see for example [34][35][36][37]), robotics [38], network communication [39], pattern recognition [40], numerical simulation of multiphasic flows [41], and analysis of the distribution of chemical compounds in the atmosphere [42]. The algorithms presented here rely on convex analysis properties to achieve a worst-case linear computation time.…”

What Shape Is Your Conjugate? A Survey of Computational Convex Analysis and Its Applications

“…The tools used are the Fenchel conjugate, the convex envelope, and other Convex Analysis arguments. The same algorithm is key to numerous numerical simulations for the Burgers' equation [15,28,87,88,103,104,205].…”

On the Hamilton–Jacobi Equation for a Harmonic Oscillator

“…For the author, a source of motivation is the system of the Hamilton-Jacobi equation S t + 1 2 |∇S| 2 + U (x) = 0 and the continuity equation ρ t + div(ρ∇S) = 0 appearing, e.g., in the high frequency approximation of Schrödinger's equation [8,10]. A special feature of this system is that mass accumulates on the singular set of S. Thus, the mass density ρ(t) is a measure even if initially the density is a regular function [2,3,5,7,14,16]. The one-dimensional case has been investigated in some detail in [5].…”