Abstract:We propose a new numerical method for $\alpha$-dissipative solutions of the Hunter-Saxton equation, where $\alpha$ belongs to $W^{1, \infty}(\mathbb{R}, [0, 1))$. The method combines a projection operator with a generalized method of characteristics and an iteration scheme, which is based on enforcing minimal time steps whenever breaking times cluster. Numerical examples illustrate that these minimal time steps increase the efficiency of the algorithm substantially. Moreover, convergence of the wave profile is shown in $C([0, T], L^{\infty}(\mathbb{R}))$ for any finite $T \geq 0$.
What problem does this paper attempt to address?
The problem that this paper attempts to solve is to develop a new numerical method to deal with the α - dissipative solutions of the Hunter - Saxton (HS) equation, where α belongs to the space \( W^{1,\infty}(\mathbb{R},[0,1)) \). Specifically, the objectives of the paper include:
1. **Solve the wave - breaking problem**: In the case of weak solutions, the HS equation usually develops singularities (called wave - breaking) in a finite time, resulting in the loss of uniqueness of the solution. Therefore, a method needs to be found to continue to describe the behavior of these weak solutions after wave - breaking.
2. **Introduce an energy - dissipation mechanism**: For each wave - breaking event, a certain proportion of energy is removed according to its spatial location. This proportion is determined by the function α, that is, at the wave - breaking point \((t_c, x_c)\), the removed energy is the \( \alpha(x_c) \)-proportion of the concentrated energy.
3. **Ensure the effectiveness and accuracy of the numerical method**: The proposed method combines projection operators, the generalized characteristic line method and an iterative scheme, and improves the algorithm efficiency by forcing a minimum time step. In addition, the paper also proves the convergence of the waveform in the space \( C([0, T], L^\infty(\mathbb{R})) \).
4. **Handle the cumulative wave - breaking time**: Since the wave - breaking time may accumulate, the paper introduces an iterative scheme to calculate the energy that should be removed at each wave - breaking, and avoids repeated calculations by extracting a finite sequence of wave - breaking times.
### Summary of Mathematical Formulas
- HS equation:
\[
u_t(t, x) + uu_x(t, x) = \frac{1}{4} \int_{-\infty}^x u_x^2(t, z) \, dz - \frac{1}{4} \int_x^\infty u_x^2(t, z) \, dz, \quad u|_{t = 0} = u_0.
\]
- Wave - breaking time function τ:
\[
\tau(\xi) =
\begin{cases}
0, & \text{if } y_{0,\xi}(\xi) = U_{0,\xi}(\xi) = 0, \\
-\frac{2y_{0,\xi}(\xi)}{U_{0,\xi}(\xi)}, & \text{if } U_{0,\xi}(\xi) < 0, \\
\infty, & \text{otherwise}.
\end{cases}
\]
- Energy - loss function V:
\[
V(t, \xi) = \int_{-\infty}^\xi \left(1 - \alpha(y(\tau(\eta), \eta)) \right) \chi_{\{\omega : t \geq \tau(\omega) > 0\}}(\eta) V_{0,\xi}(\eta) \, d\eta.
\]
Through these methods, the paper aims to provide an effective numerical framework for dealing with the α - dissipative solutions of the HS equation and ensure the accuracy and stability of the numerical results.