A structure preserving discretization for the Derrida-Lebowitz-Speer-Spohn equation based on diffusive transport

Daniel Matthes, Eva-Maria Rott, Giuseppe Savaré and André Schlichting

Numer. Math. 157(4), 1347–1395 (2025)

Abstract

We propose a spatial discretization of the fourth-order nonlinear DLSS equation on the circle. Our choice of discretization is motivated by a novel gradient flow formulation with respect to a metric that generalizes martingale transport. The discrete dynamics inherits this gradient flow structure, and in addition further properties, such as an alternative gradient flow formulation in the Wasserstein distance, contractivity in the Hellinger distance, and monotonicity of several Lypunov functionals. Our main result is the convergence in the limit of vanishing mesh size. The proof relies an a discrete version of a nonlinear functional inequality between integral expressions involving second order derivatives.

Logarithmic plot of numerical solution the DLSS equation started from a discretization of Z_m,ε^-1[ε^1/2+ (1+cos(2π x)/2)^m]^2 with parameters ε=0.001, m=8,16.
Logarithmic plot of numerical solution the DLSS equation started from a discretization of Zm,ε1[ε1/2+(1+cos(2πx)2) ⁣m] ⁣2Z_{m,\varepsilon}^{-1} \bigl[\varepsilon^{1/2}+ \bigl(\frac{1+\cos(2\pi x)}{2}\bigr)^{\!m}\bigr]^{\!2} with parameters ε=0.001\varepsilon=0.001 , m=8,16m=8,16 .

Publication history

Preprint
2023-12-20
Received
2023-12-27
Revised
2025-02-14
Accepted
2025-06-08
Published online
2025-07-09
Published
2025-07

Topics: structure-preserving discretization, optimal transport, gradient flows, numerical analysis