A Convergent Lagrangian Discretization for a Nonlinear Fourth-Order Equation

A fully discrete Lagrangian scheme for numerical solution of the nonlinear fourth-order DLSS equation in one space dimension is analyzed. The discretization is based on the equation’s gradient flow structure in the L2\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$L^2$$\end{document}-Wasserstein metric. By construction, the discrete solutions are strictly positive and mass conserving. A further key property is that they dissipate both the Fisher information and the logarithmic entropy. Our main result is a proof of convergence of fully discrete to weak solutions in the limit of vanishing mesh size. Convergence is obtained for arbitrary nonnegative, possibly discontinuous initial data with finite entropy, without any CFL-type condition. The key estimates in the proof are derived from the dissipations of the two Lyapunov functionals. Numerical experiments illustrate the practicability of the scheme.