Numerical simulation of conservation laws with moving grid nodes: Application to tsunami wave modelling

Abstract

In the present article we describe a few simple and efficient finite volume type schemes on moving grids in one spatial dimension combined with appropriate predictor-corrector method to achieve higher resolution. The underlying finite volume scheme is conservative and it is accurate up to the second order in space. The main novelty consists in the motion of the grid. This new dynamic aspect can be used to resolve better the areas with large solution gradients or any other special features. No interpolation procedure is employed, thus unnecessary solution smearing is avoided, and therefore, our method enjoys excellent conservation properties. The resulting grid is completely redistributed according the choice of the so-called monitor function. Several more or less universal choices of the monitor function are provided. Finally, the performance of the proposed algorithm is illustrated on several examples stemming from the simple linear advection to the simulation of complex shallow water waves. The exact well-balanced property is proven. We believe that the techniques described in our paper can be beneficially used to model tsunami wave propagation and run-up.

Keywords

Bibliographic record

BibTeX Citation

@article{Khakimzyanov2019numericalsimulation,
  author = {Khakimzyanov, G. and Dutykh, D. and Mitsotakis, D. and Shokina, N.},
  title = {Numerical simulation of conservation laws with moving grid nodes: Application to tsunami wave modelling},
  journal = {Geosciences},
  year = {2019},
  volume = {9},
  number = {5},
  pages = {197},
  doi = {10.3390/geosciences9050197},
  abstract = {In the present article we describe a few simple and efficient finite volume type schemes on moving grids in one spatial dimension combined with appropriate predictor-corrector method to achieve higher resolution. The underlying finite volume scheme is conservative and it is accurate up to the second order in space. The main novelty consists in the motion of the grid. This new dynamic aspect can be used to resolve better the areas with large solution gradients or any other special features. No interpolation procedure is employed, thus unnecessary solution smearing is avoided, and therefore, our method enjoys excellent conservation properties. The resulting grid is completely redistributed according the choice of the so-called monitor function. Several more or less universal choices of the monitor function are provided. Finally, the performance of the proposed algorithm is illustrated on several examples stemming from the simple linear advection to the simulation of complex shallow water waves. The exact well-balanced property is proven. We believe that the techniques described in our paper can be beneficially used to model tsunami wave propagation and run-up.},
  keywords = {Moving grid nodes, advection, wave run-up, moving grids, Shallow water waves, shallow water equations, finite volumes, Wave propagation, Predictor-corrector method, Conservation laws, adaptivity, Finite volume method, conservative finite differences, Tsunami wave modelling, Numerical simulation, Run-up, Monitor function, Well-balanced property}
}