pith. sign in

arxiv: 2512.09813 · v2 · pith:HPP72HMLnew · submitted 2025-12-10 · ❄️ cond-mat.stat-mech · physics.flu-dyn

Burgers dynamics for Poisson point process initial conditions of the Weibull class

Patrick Valageas This is my paper

Pith reviewed 2026-05-16 23:11 UTC · model grok-4.3

classification ❄️ cond-mat.stat-mech physics.flu-dyn
keywords Burgers equationPoisson point processinviscid limitself-similar evolutionvelocity statisticsshock multiplicitycorrelation functionsstretched-exponential tails
0
0 comments X

The pith

Burgers equation with Poisson point process initial conditions yields exact expressions for velocity distributions, shock statistics and correlation functions.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper derives explicit analytical formulas for the one- and two-point velocity distributions, void and shock multiplicity functions, and the velocity and density correlation functions with their power spectra, all in the inviscid limit. It further shows that the entire hierarchy of n-point distributions factorizes into a sequence of two-point conditional probabilities. These results apply to initial velocity potentials generated by a Poisson point process whose intensity follows a power law with exponent alpha greater than minus one, producing self-similar evolution in which a characteristic length scale grows as a power law of time with exponent between zero and one half and probability distributions exhibit stretched-exponential tails whose exponents range from one to infinity.

Core claim

For initial conditions defined by a Poisson point process whose intensity follows a power law with exponent alpha greater than minus one, the geometrical construction of solutions in terms of first-contact parabolas supplies the exact inviscid solution to the one-dimensional Burgers equation. This construction directly yields closed-form expressions for the one- and two-point velocity probability distributions, the multiplicity functions of voids and shocks, and the two-point velocity and density correlation functions together with their power spectra. The full hierarchy of n-point distributions factorizes into a chain of two-point conditional probabilities, and the resulting statistics are,

What carries the argument

The first-contact parabolas geometrical construction, which identifies the solution by locating the lowest parabola that touches the initial potential at its points of contact.

If this is right

  • The velocity and density fields possess explicitly computable two-point correlation functions and power spectra at all times.
  • The multiplicity functions counting voids and shocks are given by closed-form expressions that depend only on the power-law exponent alpha.
  • All higher-order n-point distributions reduce exactly to a product of two-point conditional distributions.
  • The characteristic length scale of the velocity field grows as t to the power beta where beta lies between zero and one half.
  • The tails of all probability distributions are stretched exponentials whose decay exponents range continuously from one to infinity.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • The factorization into conditional two-point distributions may extend to other hyperbolic conservation laws that admit similar geometrical solutions.
  • Direct numerical simulations of the Burgers equation with the same Poisson initial conditions would provide a quantitative test of the predicted power-law growth of the length scale.
  • The stretched-exponential family obtained here could serve as a benchmark for approximate closures used in higher-dimensional or forced versions of the equation.
  • Because the initial conditions belong to the Weibull class, the same geometrical method may apply to related extreme-value problems in one-dimensional nonlinear wave equations.

Load-bearing premise

The first-contact parabola construction supplies the exact inviscid solution for initial conditions generated by a Poisson point process with power-law intensity.

What would settle it

A numerical integration of the inviscid Burgers equation starting from a Poisson point process with power-law intensity alpha greater than minus one, followed by direct comparison of the measured one-point velocity distribution against the predicted stretched-exponential form.

Figures

Figures reproduced from arXiv: 2512.09813 by Patrick Valageas.

Figure 1
Figure 1. Figure 1: A realization of the system for the cases [PITH_FULL_IMAGE:figures/full_fig_p004_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: One-point probability distribution P0(q) = P0(v) of the Lagrangian coordinate q, or of the velocity v, from Eq.(16). We display our results on a linear scale (left panel) and a logarithmic scale (right panel), for the cases α = −0.5, 0, and 2. The dotted lines in the right panel are the asymptotic results (18). F. Limit α → ∞ In the limit α → ∞, the Poisson intensity (12) implies that very few points lie b… view at source ↗
Figure 3
Figure 3. Figure 3: Void probability Pvoid(x) from Eq.(33), for the cases α = −1/2, 0, and 2, as in [PITH_FULL_IMAGE:figures/full_fig_p008_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: Left panel: cumulative void multiplicity function nvoid(> x) from Eq.(35). Middle panel: void multiplicity function nvoid(x) from Eq.(35). Right panel: number density of voids Nvoid as a function of α (solid blue line). In the left and middle panels, the dotted lines are the stretched exponentials associated with Eq.(36). In the right panel the dotted and dashed lines are the asymptotic regimes (38). 2. Mu… view at source ↗
Figure 5
Figure 5. Figure 5: Velocity correlation Bv(x) for the cases α = −0.5, 0, and 2. In the right panel the dotted lines are the asymptotic stretched exponentials (44). Its limiting behaviors are α → −1 : Nvoid ∼ 1 α + 1 , and for α → ∞ : Nvoid ∼ p α/π, (38) so that the number of voids per unit length diverges in both limits α → −1 and α → ∞. The mean void size is given by ⟨x⟩void = R ∞ 0 dx nvoid(x) x R ∞ 0 dx nvoid(x) = 1 Nvoid… view at source ↗
Figure 6
Figure 6. Figure 6: Velocity power spectrum E(k) on linear and logarithmic scales. where Aν is defined in Eq.(28). Performing the integration by parts over q ′ ⋆ yields x ≥ 0 : Bv(x) = 1 α + 1 d dx [xRα+1(x)] , (43) with Rν defined in Eq.(28). By parity, Bv(−x) = Bv(x) = Bv(|x|). This gives the small- and large-scale asymptotic behaviors |x| ≪ 1 : Bv(x) = Rα+1(0) α + 1 + 2R′ α+1(0) α + 1 |x| + . . . , |x| ≫ 1 : Bv(x) ∼ −x −2α… view at source ↗
Figure 7
Figure 7. Figure 7: Left and middle panels: density correlation function ξ ̸= δ (x) on linear and logarithmic scales for x > 0. Right panel: density power spectrum Pδ(k). D. Density correlation and power spectrum The conservation of matter, encoded in the continuity equation, implies that the density field ρ(x) (normalized by the mean density of the system) and the density contrast δ = ρ − 1 are given by ρ(x) = dq dx, δ(x) = … view at source ↗
Figure 8
Figure 8. Figure 8: Left panel: probability distribution P ̸= x (q) for the case α = 0, for the three scales x = 0.2, 1, and 3. The black dotted lines are the small-q linear asymptote (58) and the black dashed line for x = 3 is the large-separation asymptote (59). Right panel: probability distribution P ̸= x (ρ) for the same cases. E. Lagrangian increment We now consider the distribution of the Lagrangian increment q = q2 − q… view at source ↗
Figure 9
Figure 9. Figure 9: Left panel: probability distribution P ̸= q (x) of the Eulerian increment x, for the case α = 0 and the three scales q = 0.2, 1, and 3, as in [PITH_FULL_IMAGE:figures/full_fig_p015_9.png] view at source ↗
read the original abstract

We derive the statistical properties of one-dimensional Burgers dynamics with stochastic initial conditions for the velocity potential defined by a Poisson point process whose intensity follows a power law with exponent $\alpha > -1$. Working in the inviscid limit and exploiting the geometrical construction of solutions in terms of first-contact parabolas, we derive explicit analytical expressions for a broad set of statistical quantities. These include the one- and two-point probability distributions of the velocity, the multiplicity functions of voids and shocks, and the velocity and density correlation functions together with their associated power spectra. We also show that the full hierarchy of $n$-point distributions factorizes into a sequence of two-point conditional probabilities. This class of initial conditions leads to self-similar evolution and produces probability distributions characterized by stretched-exponential tails, with tail exponents spanning the full range from unity to infinity. The associated characteristic length scale grows as a power law of time, with an exponent lying between zero and one half.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

2 major / 3 minor

Summary. The manuscript analyzes the one-dimensional inviscid Burgers equation with initial velocity potential generated by a Poisson point process whose intensity scales as a power law |x|^α with α > -1. Exploiting the standard geometrical construction via first-contact parabolas, the authors obtain explicit analytical expressions for the one- and two-point velocity PDFs, the multiplicity functions of voids and shocks, the velocity and density correlation functions together with their power spectra, and demonstrate that the full n-point hierarchy factorizes into a chain of two-point conditional probabilities. The resulting statistics exhibit self-similar evolution, stretched-exponential tails whose exponents range from 1 to ∞, and a characteristic length that grows as a power law in time with exponent in (0, 1/2).

Significance. If the derivations hold, the paper supplies exact, closed-form results for a wide family of random initial data in Burgers turbulence, a canonical model for shock formation and nonlinear wave steepening. The explicit formulas for PDFs, multiplicity functions, spectra, and the factorization property furnish precise benchmarks for numerical codes and analytic insight into intermittency and tail behavior. The work extends the classical geometrical approach to the Weibull-class Poisson processes while preserving parameter-free derivations once α is fixed, which is a notable strength.

major comments (2)
  1. [§3.1, Eq. (8)] §3.1, Eq. (8): the claimed factorization of the n-point velocity distribution into successive two-point conditionals is derived under the independent-increments property of the Poisson process, but the proof sketch does not explicitly verify that the first-contact parabola ordering preserves this factorization when multiple parabolas compete at the same location; an expanded step showing the conditional independence would remove any ambiguity.
  2. [§4.3, Eq. (22)] §4.3, Eq. (22): the power spectrum of the density field is stated to decay as k^{-(2+α)} at large k, yet the derivation appears to interchange the order of averaging over the Poisson process and the Fourier transform without justifying the dominated-convergence step for α near -1; a brief estimate of the remainder term would confirm the result holds uniformly in the allowed range.
minor comments (3)
  1. [Figure 3] Figure 3: the curves for different α are plotted on the same panel without a clear legend or line-style key; adding explicit labels for α = 0, 1, 2 would improve readability.
  2. [§2.2] §2.2: the definition of the stretched-exponential tail exponent β(α) is introduced only in the text; placing the explicit formula β(α) = (α+1)/2 next to the first appearance would help readers track the dependence on the single free parameter.
  3. [References] The reference list omits the classic work of Sinai (1992) on the geometrical solution of Burgers with random initial data; adding this citation would place the present results in clearer context.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the careful reading of our manuscript and the constructive comments. We address each of the major comments below and have made the corresponding revisions.

read point-by-point responses

  1. Referee: §3.1, Eq. (8): the claimed factorization of the n-point velocity distribution into successive two-point conditionals is derived under the independent-increments property of the Poisson process, but the proof sketch does not explicitly verify that the first-contact parabola ordering preserves this factorization when multiple parabolas compete at the same location; an expanded step showing the conditional independence would remove any ambiguity.

    Authors: We appreciate this observation. The independent-increments property of the Poisson point process ensures that increments over disjoint intervals are independent, and the first-contact parabola construction selects the global infimum of the parabolic potentials. Because the selection at each step depends only on the current position and the future increments (which are independent of the past by the Poisson property), the conditional distribution of the next contact point given the previous one remains independent of the earlier history. We have expanded the derivation in the revised §3.1 with an explicit paragraph that walks through this conditional-independence argument step by step, confirming that the n-point factorization is preserved. revision: yes

  2. Referee: §4.3, Eq. (22): the power spectrum of the density field is stated to decay as k^{-(2+α)} at large k, yet the derivation appears to interchange the order of averaging over the Poisson process and the Fourier transform without justifying the dominated-convergence step for α near -1; a brief estimate of the remainder term would confirm the result holds uniformly in the allowed range.

    Authors: We thank the referee for this technical remark. The density field is a sum of positive contributions from the Poisson points, so the interchange between expectation and Fourier transform is justified by the monotone convergence theorem. For the boundary behavior as α → -1^+, we have added a short estimate in the revised §4.3: the remainder after truncation at a large but finite number of points is bounded by an integrable term whose Fourier transform decays at least as k^{-(2+α)} times a slowly varying factor, which remains o(k^{-(2+α)}) uniformly for α ≥ -1 + δ with any fixed δ > 0. This confirms the stated large-k asymptotics throughout the allowed range. revision: yes

Circularity Check

0 steps flagged

No significant circularity; derivation is self-contained

full rationale

The paper starts from the standard first-contact parabola construction, which is known to yield the exact weak solution of the inviscid Burgers equation for any suitable initial potential, and applies it to a Poisson point process whose intensity is a power law with fixed input exponent alpha > -1. All explicit expressions for PDFs, multiplicity functions, correlation functions, and spectra are obtained by direct integration over the geometry of contact parabolas and the independent-increments property of the Poisson process. No parameters are fitted to data, no predictions are defined in terms of themselves, and no load-bearing step reduces to a self-citation whose validity is presupposed by the present work. The restriction alpha > -1 is a regularity condition on the input measure, not a derived result.

Axiom & Free-Parameter Ledger

1 free parameters · 1 axioms · 0 invented entities

The central claim rests on the standard geometrical solution method for the inviscid 1D Burgers equation applied to a Poisson point process with power-law intensity; alpha > -1 is the defining parameter of the initial-condition class.

free parameters (1)
  • alpha
    Power-law exponent of the Poisson intensity, with alpha > -1 defining the Weibull class; treated as an input parameter of the model family.
axioms (1)
  • domain assumption Geometrical construction via first-contact parabolas yields the exact solution in the inviscid limit
    Invoked throughout to obtain the velocity field from the initial potential and to derive all statistical quantities.

pith-pipeline@v0.9.0 · 5460 in / 1407 out tokens · 37131 ms · 2026-05-16T23:11:40.319296+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

What do these tags mean?
matches
The paper's claim is directly supported by a theorem in the formal canon.
supports
The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
extends
The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
uses
The paper appears to rely on the theorem as machinery.
contradicts
The paper's claim conflicts with a theorem or certificate in the canon.
unclear
Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.

Reference graph

Works this paper leans on

34 extracted references · 34 canonical work pages · 2 internal anchors

  1. [1]

    For the first term in Eq.(24), the factorλ(ψ⋆)actually impliesψ ⋆ >0, since this configuration requires a Poisson point at the intersection(q⋆, ψ⋆)

    = Z ∞ −∞ dq′ ⋆ Z ∞ ψmin(q′⋆) dψ⋆ λ(ψ⋆)δD(q′ 1 −q ′ ⋆)δD(q′ 2 −q ′ ⋆) +xθ(q ′ 1 < q ′ ⋆)θ(q′ 2 > q ′ ⋆)ψ−(q′ 1)αψ+(q′ 2)α e−I,(24) 7 with ψ−(q′) = max 0, ψ⋆ + 1 2 q′ ⋆ + x 2 2 − 1 2 q′ + x 2 2 , ψ +(q′) = max 0, ψ⋆ + 1 2 q′ ⋆ − x 2 2 − 1 2 q′ − x 2 2 , I(ψ ⋆, q′ ⋆) = 1 α+ 1 "Z q′ ⋆ −∞ dq′ψ−(q′)α+1 + Z ∞ q′⋆ dq′ψ+(q′)α+1 # ,(25) and |q′ ⋆| ≥ x 2 :ψ min(q′ ⋆...

    work page

  2. [2]

    If the two parabolas share the same contact point,q1 =q 2 =q ⋆, then the density vanishes,ρ= 0, and the interval [x1, x2]contains no matter

    Probability of an empty interval The overdensity within the Eulerian interval[x1, x2]is defined by ρx1,x2 = q2 −q 1 x2 −x 1 ≥0,(32) where the density has been rescaled by the mean densityρ0, though we retain the simpler notationρin the following. If the two parabolas share the same contact point,q1 =q 2 =q ⋆, then the density vanishes,ρ= 0, and the interv...

    work page

  3. [3]

    It is related to the void probabilityPvoid(x)by Pvoid(x) = Z ∞ x dx′nvoid(x′) (x′ −x),whencen void(> x) =− dPvoid dx =−R ′ α(x), n void(x) = d2Pvoid dx2 =R ′′ α(x)

    Multiplicity function of voids and distance between shocks Letn void(x)dxdenote the number of voids per unit length with sizes in the interval[x, x+dx]. It is related to the void probabilityPvoid(x)by Pvoid(x) = Z ∞ x dx′nvoid(x′) (x′ −x),whencen void(> x) =− dPvoid dx =−R ′ α(x), n void(x) = d2Pvoid dx2 =R ′′ α(x). (35) Using the asymptotic expression (3...

    work page

  4. [4]

    Configurations withψ ⋆ ≤0do not contribute, as in this case the parabolic arcsP 1 andP 2 in the upper half-plane are disjoints. Then, the integral over either symmetric arc yields a vanishing mean velocity: the contributions from the two sides ofx1 alongP 1 have identical weights but opposite velocitiesv 1 =x 1 −q 1. As for the Fréchet-type case studied i...

    work page

  5. [5]

    J. M. Burgers,The Nonlinear Diffusion Equation(Springer Netherlands, 1974)

    work page 1974

  6. [6]

    Hopf, The partial differential equation ut + uux =µxx, Communications on Pure and Applied Mathematics3, 201 (1950)

    E. Hopf, The partial differential equation ut + uux =µxx, Communications on Pure and Applied Mathematics3, 201 (1950)

    work page 1950

  7. [7]

    J. D. Cole, On a quasi-linear parabolic equation occurring in aerodynamics, Quarterly of Applied Mathematics9, 225 (1951)

    work page 1951

  8. [8]

    S. N. Gurbatov, A. I. Saichev, S. N. Gurbatov, and A. I. Saichev, Degeneracy of one-dimensional acoustic turbulence under large Reynolds numbers, ZhETF80, 689 (1981)

    work page 1981

  9. [9]

    G. B. Whitham,Linear and Nonlinear Waves(wiley, 1999)

    work page 1999

  10. [10]

    Frisch and J

    U. Frisch and J. Bec,New trends in turbulence Turbulence: nouveaux aspects(Springer, Berlin, Heidelberg, 2001) pp. 341–383, arXiv:0012033 [nlin]

    work page 2001

  11. [11]

    Bec and K

    J. Bec and K. Khanin, Burgers turbulence, Physics Reports447, 1 (2007)

    work page 2007

  12. [12]

    R. H. Kraichnan, Lagrangian-History Statistical Theory for Burgers’ Equation, The Physics of Fluids11, 265 (1968). 17

    work page 1968

  13. [13]

    Kida, Asymptotic properties of Burgers turbulence, Journal of Fluid Mechanics93, 337 (1979)

    S. Kida, Asymptotic properties of Burgers turbulence, Journal of Fluid Mechanics93, 337 (1979)

    work page 1979

  14. [14]

    Frachebourg, P

    L. Frachebourg, P. A. Martin, and J. Piasecki, Ballistic aggregation: a solvable model of irreversible many particles dynamics, Physica A: Statistical Mechanics and its Applications279, 69 (2000)

    work page 2000

  15. [15]

    Valageas, Ballistic aggregation for one-sided Brownian initial velocity, Physica A: Statistical Mechanics and its Appli- cations388, 1031 (2009)

    P. Valageas, Ballistic aggregation for one-sided Brownian initial velocity, Physica A: Statistical Mechanics and its Appli- cations388, 1031 (2009)

    work page 2009

  16. [16]

    S. N. Gurbatov, A. I. Saichev, S. F. Shandarin, S. N. Gurbatov, A. I. Saichev, and S. F. Shandarin, The large-scale structure of the universe in the frame of the model equation of non-linear diffusion, MNRAS236, 385 (1989)

    work page 1989

  17. [17]

    Gurbatov, A

    S. Gurbatov, A. Malakhov, and A. I. Saichev,Nonlinear random waves and turbulence in nondispersive media : waves, rays, particles(Manchester University Press ; Distributed exclusively in the US and Canada by St. Martin’s Press, 1991)

    work page 1991

  18. [18]

    Vergassola, B

    M. Vergassola, B. Dubrulle, U. Frisch, A. Noullez, M. Vergassola, B. Dubrulle, U. Frisch, and A. Noullez, Burgers’ equation, Devil’s staircases and the mass distribution for large-scale structures., A&A289, 325 (1994)

    work page 1994

  19. [19]

    Density fields and halo mass functions in the Geometrical Adhesion toy Model

    P. Valageas and F. Bernardeau, Density fields and halo mass functions in the geometrical adhesion toy model, Physical Review D83, 043508 (2011), arXiv:1009.1974

    work page internal anchor Pith review Pith/arXiv arXiv 2011

  20. [20]

    Y. B. Zel’dovich, Gravitational instability: An approximate theory for large density perturbations., A&A5, 84 (1970)

    work page 1970

  21. [21]

    S. F. Shandarin and Y. B. Zeldovich, The large-scale structure of the universe: Turbulence, intermittency, structures in a self-gravitating medium, Reviews of Modern Physics61, 185 (1989)

    work page 1989

  22. [22]

    Z. S. She, E. Aurell, and U. Frisch, The inviscid Burgers equation with initial data of Brownian type, Communications in Mathematical Physics148, 623 (1992)

    work page 1992

  23. [23]

    S. N. Gurbatov, S. I. Simdyankin, E. Aurell, U. Frisch, and G. Tóth, On the decay of Burgers turbulence, Journal of Fluid Mechanics344, 339 (1997), arXiv:9709002 [physics]

    work page 1997

  24. [24]

    Bertoin, The inviscid burgers equation with Brownian initial velocity, Communications in Mathematical Physics193, 397 (1998)

    J. Bertoin, The inviscid burgers equation with Brownian initial velocity, Communications in Mathematical Physics193, 397 (1998)

    work page 1998

  25. [25]

    Valageas, Statistical properties of the burgers equation with brownian initial velocity, Journal of Statistical Physics134, 589 (2009)

    P. Valageas, Statistical properties of the burgers equation with brownian initial velocity, Journal of Statistical Physics134, 589 (2009)

    work page 2009

  26. [26]

    Frachebourg and P

    L. Frachebourg and P. A. Martin, Exact statistical properties of the Burgers equation, Journal of Fluid Mechanics417, 323 (2000), arXiv:9905056 [cond-mat]

    work page 2000

  27. [27]

    Some statistical properties of the Burgers equation with white-noise initial velocity

    P. Valageas, Some statistical properties of the Burgers equation with white-noise initial velocity, Journal of Statistical Physics137, 729 (2009), arXiv:0903.0956

    work page internal anchor Pith review Pith/arXiv arXiv 2009

  28. [28]

    S. A. Molchanov, D. Surgailis, and W. A. Woyczynski, The large-scale structure of the universe and quasi-Voronoi tessel- lation of shock fronts in forced Burgers turbulence in Rd, Annals of Applied Probability7, 200 (1997)

    work page 1997

  29. [29]

    Gueudré and P

    T. Gueudré and P. Le Doussal, Statistics of shocks in a toy model with heavy tails, Physical Review E - Statistical, Nonlinear, and Soft Matter Physics89, 42111 (2014)

    work page 2014

  30. [30]

    Valageas, Burgers dynamics for poisson point process initial conditions, (2025), arXiv:2511.03647

    P. Valageas, Burgers dynamics for poisson point process initial conditions, (2025), arXiv:2511.03647

    work page arXiv 2025

  31. [31]

    Bernard and K

    D. Bernard and K. G. Gawe¸dzki, Scaling and exotic regimes in decaying Burgers turbulence, J. Phys. A: Math. Gen31, 8735 (1998)

    work page 1998

  32. [32]

    Bauer and D

    M. Bauer and D. Bernard, Sailing the deep blue sea of decaying burgers turbulence, Journal of Physics A: Mathematical and General32, 5179 (1999)

    work page 1999

  33. [33]

    S. N. Gurbatov, O. V. Rudenko, and A. I. Saichev,Waves and Structures in Nonlinear Nondispersive Media, Nonlinear Physical Science (Springer Berlin Heidelberg, Berlin, Heidelberg, 2011)

    work page 2011

  34. [34]

    S. A. Molchanov, D. Surgailis, and W. A. Woyczynski, Hyperbolic asymptotics in Burgers’ turbulence and extremal processes, Communications in Mathematical Physics168, 209 (1995)

    work page 1995