Numerical entropy production as smoothness indicator for shallow water equations

2020 ◽  
Vol 61 ◽  
pp. 398-415
Author(s):  
Sudi Mungkasi ◽  
Stephen Gwyn Roberts

The numerical entropy production (NEP) for shallow water equations (SWE) is discussed and implemented as a smoothness indicator. We consider SWE in three different dimensions, namely, one-dimensional, one-and-a-half-dimensional, and two-dimensional SWE. An existing numerical entropy scheme is reviewed and an alternative scheme is provided. We prove the properties of these two numerical entropy schemes relating to the entropy steady state and consistency with the entropy equality on smooth regions. Simulation results show that both schemes produce NEP with the same behaviour for detecting discontinuities of solutions and perform similarly as smoothness indicators. An implementation of the NEP for an adaptive numerical method is also demonstrated. doi:10.1017/S1446181119000154

2019 ◽  
Vol 61 (4) ◽  
pp. 398-415
Author(s):  
SUDI MUNGKASI ◽  
STEPHEN GWYN ROBERTS

The numerical entropy production (NEP) for shallow water equations (SWE) is discussed and implemented as a smoothness indicator. We consider SWE in three different dimensions, namely, one-dimensional, one-and-a-half-dimensional, and two-dimensional SWE. An existing numerical entropy scheme is reviewed and an alternative scheme is provided. We prove the properties of these two numerical entropy schemes relating to the entropy steady state and consistency with the entropy equality on smooth regions. Simulation results show that both schemes produce NEP with the same behaviour for detecting discontinuities of solutions and perform similarly as smoothness indicators. An implementation of the NEP for an adaptive numerical method is also demonstrated.


2015 ◽  
Vol 771 ◽  
pp. 157-160 ◽  
Author(s):  
Sudi Mungkasi

We extend the application of numerical entropy production, as a smoothness indicator, from conservation laws to balance laws. We aim to indicate the smoothness of solutions to the shallow water equations involving varying width, which are a system of balance laws. The numerical entropy production appears to be accurate to detect discontinuities. As a numerical test, a radial dam break is considered. We assume that there is a higher level of water inside a radial dam than water outside the dam wall. If the radial dam is totally broken, then water flows from inside to outside. The flow results in a solution having shock discontinuities. Finding the positions of the discontinuities is our interest. They are the positions where numerical solutions, such as those generated by a finite volume method, decrease their accuracy. Detecting the position of the discontinuity can help in the improvement of the numerical solution in terms of its accuracy. We obtain that the numerical entropy production is simple to implement but give an accurate detection. The discontinuity of the stage (free water surface) is clearly detected by large values of the numerical entropy production as the smoothness indicator.


2016 ◽  
Vol 2016 ◽  
pp. 1-7 ◽  
Author(s):  
Sudi Mungkasi

This paper presents a numerical entropy production (NEP) scheme for two-dimensional shallow water equations on unstructured triangular grids. We implement NEP as the error indicator for adaptive mesh refinement or coarsening in solving the shallow water equations using a finite volume method. Numerical simulations show that NEP is successful to be a refinement/coarsening indicator in the adaptive mesh finite volume method, as the method refines the mesh or grids around nonsmooth regions and coarsens them around smooth regions.


2014 ◽  
Vol 178 (3) ◽  
pp. 278-298 ◽  
Author(s):  
Yu. A. Chirkunov ◽  
S. Yu. Dobrokhotov ◽  
S. B. Medvedev ◽  
D. S. Minenkov

Sign in / Sign up

Export Citation Format

Share Document