Numerical entropy production as smoothness indicator for shallow water equations

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

NUMERICAL ENTROPY PRODUCTION AS SMOOTHNESS INDICATOR FOR SHALLOW WATER EQUATIONS

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.

HPM-Based Dynamic Wavelet Transform and Its Application in Image Denoising

Wavelet-based multiscale interpolation operator is often employed to construct the adaptive numerical method for PDEs, in which the computational complexity of the wavelet transform is one of the main factors affecting the algorithm efficiency. As the wavelet transform just acts as the detector of the characteristic points in the interpolation operator, the multiscale wavelet interpolation operator can be viewed as a nonlinear problem. Based on this assumption, we construct an approximate dynamic interpolation operator with the homotopy perturbation method (HPM), which decreases the computational complexity of the wavelet transform appearing in the wavelet interpolation operator fromO((1/3)42J−1)toO(4J), whereJis the amount of the wavelet scales. Then an adaptive algorithm solving the Perona-Malik model on image denoising is constructed with the HPM-based interpolation operator. Last, the quasi-Shannon wavelet is employed to design the experiments on the medical image and some artificial images denoising. The experiment results show that the simplified wavelet interpolation operator based on HPM possesses the adaptability and nonsensitivity to the time step, which is helpful to improve the algorithm efficiency. This illustrates that the HPM-based wavelet interpolation operator is an effective tool to solve the problems in image processing.