A Convergent Numerical Algorithm for $\alpha$-Dissipative Solutions of the Hunter-Saxton Equation

Abstract

A convergent numerical method for $\alpha$-dissipative solutions of the Hunter-Saxton equation is derived. The method is based on applying a tailor-made projection operator to the initial data, and then solving exactly using the generalized method of characteristics. The projection step is the only step that introduces any approximation error. It is therefore crucial that its design ensures not only a good approximation of the initial data, but also that errors due to the energy dissipation at later times remain small. Furthermore, it is shown that the main quantity of interest, the wave profile, converges in $L^{\infty }$ for all $t\ge 0$, while a subsequence of the energy density converges weakly for almost every time.

arXiv comment51 pages, 9 figures. Corrected some typos in the appendix

References in corpus (5)

• Global conservative solutions of the Camassa-Holm equation for initial data nonvanishing asymptotics
• Existence and Lipschitz stability for $\alpha$-dissipative solutions of the two-component Hunter-Saxton system
• Numerical conservative solutions of the Hunter--Saxton equation
• Lipschitz Stability for the Hunter-Saxton Equation
• On the equivalence of Eulerian and Lagrangian variables for the two-component Camassa-Holm system

Cited by in corpus (2)

• A numerical view on α-dissipative solutions of the Hunter-Saxton equation
• Rate of convergence for numerical $\alpha$-dissipative solutions of the Hunter-Saxton equation