A Numerical Algorithm for Computing α-Dissipative Solutions of the Hunter-Saxton Equation

A novel numerical algorithm is proposed for computing α-dissipative solutions of the Hunter-Saxton equation, where α belongs to W^{1,∞}(R, [0,1)). The method combines a projection operator with a generalized method of characteristics and an iteration scheme to handle the unknown energy loss at wave breaking occurrences.