Numerical boundary treatment for shock propagation in the fractional KdV-Burgers equation

In this paper, we propose a numerical boundary treatment for simulating shock propagation in the fractional KdV-Burgers equation, based on a machine learning strategy. Numerical boundary treatment is particularly challenging due to the nonlinear and global interaction among the fractional diffusion, convection and dispersion. We select a suitable number of boundary points and interior points, and correct the numerical value at the boundary ones by linear combination of those at the interior ones. The large set of parameters are trained from a numerical reference shock solution connecting the end states 1 and 0, by ridge regression. Numerical tests demonstrate that the boundary treatment is effective in reflection suppression, and preserves the correct shock speed even under perturbed initial profile. Furthermore, using the above parameters, we construct a modified boundary treatment to simulate a range of different end states, with effectiveness checked numerically.