Rate of convergence for numerical $\alpha$-dissipative solutions of the Hunter-Saxton equation

We prove that $\alpha$-dissipative solutions to the Cauchy problem of the Hunter-Saxton equation, where $\alpha \in W^{1, \infty}(\mathbb{R}, [0, 1))$, can be computed numerically with order $\mathcal{O}(\Delta x^{{1}/{8}}+\Delta x^{{\beta}/{4}})$ in $L^{\infty}(\mathbb{R})$, provided there exist constants $C > 0$ and $\beta \in (0, 1]$ such that the initial spatial derivative $\bar{u}_{x}$ satisfies $\|\bar{u}_x(\cdot + h) - \bar{u}_x(\cdot)\|_2 \leq Ch^{\beta}$ for all $h \in (0, 2]$. The derived convergence rate is exemplified by a number of numerical experiments.
Comment: 45 pages, 10 figures; Corrected a crucial typo in Definition 3.5 and some other typos in the text