# American Institute of Mathematical Sciences

## Damping, stabilization, and numerical filtering for the modeling and the simulation of time dependent PDEs

 Laboratoire LAMFA (UMR CNRS 7352), Université de Picardie Jules Verne, 33 rue Saint Leu, 80039 Amiens Cédex, France

* Corresponding author: Jean-Paul Chehab

To the memory of Ezzeddine Zahrouni (1963-2018)

Received  June 2020 Revised  December 2020 Published  January 2021

We present here different situations in which the filtering of high or low modes is used either for stabilizing semi-implicit numerical schemes when solving nonlinear parabolic equations, or for building adapted damping operators in the case of dispersive equation. We consider numerical filtering provided by mutigrid-like techniques as well as the filtering resulting from operator with monotone symbols. Our approach applies to several discretization techniques and we focus on finite elements and finite differences. Numerical illustrations are given on Cahn-Hilliard, Korteweig-de Vries and Kuramoto-Sivashinsky equations.

Citation: Jean-Paul Chehab. Damping, stabilization, and numerical filtering for the modeling and the simulation of time dependent PDEs. Discrete & Continuous Dynamical Systems - S, doi: 10.3934/dcdss.2021002
##### References:

show all references

##### References:
Different filters
Different basis of orthogonal polynomials: Chebyshev polynomials (top left), Legendre polynomials (top right) and Fourier polynomials (bottom)
Hierarchy of the triangulation of $\Omega = ]0,1[^2$: solid line (coarse triangulation), dashed line (complementary triangulation)
3D output and iso-values for $u(x,y) = \cos \bigl(5(1-x^2-y^2)\bigr)$ on the unit disk. The function $u$ with $\mathbb{P}_2$ elements (top) and the $z_h$ components (bottom)
Function $u(x,y) = \cos \bigl(5(1-x^2-y^2)\bigr)$ on the unit disk. Eigenfunction (Fourier) coefficients on the fine mesh for the original function $u_h$ (red line) and the associated correction $z_h$ (blue line)
Decomposition of the signal $u(x) = \sin(2\pi x)+\sin(6\pi x)+\sin(12\pi x)+0.1\sin(20\pi x)+0.1\sin(30\pi x)+0.1\sin(120\pi x)$, $N = 100, m = 8$. Original signal (top), high frequencies fluctuent part (bottom left) and low frequencies mean part (bottom right)
Dirichlet BC (Top Left), Neumannn BC (Top Right), Periodic BC (Bottom). The original values are given at grid points $\times$, interpolated values are computed at grid points $o$. The boundary points are marked by a dot, in the Dirichlet case
Low-pass filters: (left) - Exponential of low pass filters = high-pass filter (right)
Coarse mesh (left) and fine mesh (right)
Initial solution (left) and Final Solution (right) - \hskip 1.cm $\mathbb{P}_2$ Elements. $\epsilon=0.08$, $\Delta t =1.e-4, \tau=4\epsilon^2, T=0.012$
Energy (left) and mass (right) vs time - \hskip 2.cm $\mathbb{P}_2$ Elements. $\epsilon=0.08$, $\Delta t =1.e-4, \tau=4\epsilon^2, \ T=0.012$
$\|u\|_{L^2}$ vs time with $u_0 = \cos(2\pi x/L)+\cos(6\pi x/L)+\cos(12\pi x/L)+\cos(20\pi x/L)$, $(\tau_0,\tau_1) = (10,100)$ (left) $(\tau_0,\tau_1) = (100,10)$ (right), $\Delta t = 1.e-1$, $T = 10$, $L = 100$
KdV $(\tau_0,\tau_1) = (0,0)$, $\mathbb{P}_1$ Elements. $\Delta t = 1.e-2$, $T = 40$. Mass $\int_0^Ludx$ (left) and $L^2$-norm $|u|_{L^2}$ (right) vs time
KdV $(\tau_0,\tau_1) = (0,100)$ Low pass damping, $\mathbb{P}_1$ Elements. $\Delta t = 1.e-2$, $T = 40$. Mass $\int_0^Ludx$ (left) and $L^2$-norm $|u|_{L^2}$ (right) vs time
KdV $(\tau_0,\tau_1) = (0,100)$ Low pass damping, $\mathbb{P}_1$ Elements. $\Delta t = 1.e-2$, $T = 4$. Mass $\int_0^Ludx$ (left) and $L^2$-norm $|u|_{L^2}$ (right) vs time
Heat Equation - $L^{\infty}$-norm of the error vs time - $n = 100, m = 4, \Delta t = 9.95 \, 10^{-5}, \tau = 1.6,10^{4}$
KSE Low and high frequency components of the solution at final time $T=140$ (left), time evolution of the mean value (right) - $n=128, \ m=8$, $\Delta t=0.01$, $L=10$ (line 1), $L=20$ (line 2)
KSE Low and high frequency components of the solution at final time $T=140$ (left), time evolution of the mean value (right) - $\Delta t=0.01$. Line 1: $L=50$, $n=128, \ m=8$; line 2: $L=100$, $n=256, \, m=10$
 [1] Takeshi Fukao, Shuji Yoshikawa, Saori Wada. Structure-preserving finite difference schemes for the Cahn-Hilliard equation with dynamic boundary conditions in the one-dimensional case. Communications on Pure & Applied Analysis, 2017, 16 (5) : 1915-1938. doi: 10.3934/cpaa.2017093 [2] Andreas C. Aristotelous, Ohannes Karakashian, Steven M. Wise. A mixed discontinuous Galerkin, convex splitting scheme for a modified Cahn-Hilliard equation and an efficient nonlinear multigrid solver. Discrete & Continuous Dynamical Systems - B, 2013, 18 (9) : 2211-2238. doi: 10.3934/dcdsb.2013.18.2211 [3] Laurence Cherfils, Madalina Petcu, Morgan Pierre. A numerical analysis of the Cahn-Hilliard equation with dynamic boundary conditions. Discrete & Continuous Dynamical Systems - A, 2010, 27 (4) : 1511-1533. doi: 10.3934/dcds.2010.27.1511 [4] Desheng Li, Xuewei Ju. On dynamical behavior of viscous Cahn-Hilliard equation. Discrete & Continuous Dynamical Systems - A, 2012, 32 (6) : 2207-2221. doi: 10.3934/dcds.2012.32.2207 [5] Laurence Cherfils, Alain Miranville, Sergey Zelik. On a generalized Cahn-Hilliard equation with biological applications. Discrete & Continuous Dynamical Systems - B, 2014, 19 (7) : 2013-2026. doi: 10.3934/dcdsb.2014.19.2013 [6] Álvaro Hernández, Michał Kowalczyk. Rotationally symmetric solutions to the Cahn-Hilliard equation. Discrete & Continuous Dynamical Systems - A, 2017, 37 (2) : 801-827. doi: 10.3934/dcds.2017033 [7] Pierre Garnier. Damping to prevent the blow-up of the korteweg-de vries equation. Communications on Pure & Applied Analysis, 2017, 16 (4) : 1455-1470. doi: 10.3934/cpaa.2017069 [8] Georgia Karali, Yuko Nagase. On the existence of solution for a Cahn-Hilliard/Allen-Cahn equation. Discrete & Continuous Dynamical Systems - S, 2014, 7 (1) : 127-137. doi: 10.3934/dcdss.2014.7.127 [9] Fausto Cavalli, Giovanni Naldi. A Wasserstein approach to the numerical solution of the one-dimensional Cahn-Hilliard equation. Kinetic & Related Models, 2010, 3 (1) : 123-142. doi: 10.3934/krm.2010.3.123 [10] Muhammad Usman, Bing-Yu Zhang. Forced oscillations of the Korteweg-de Vries equation on a bounded domain and their stability. Discrete & Continuous Dynamical Systems - A, 2010, 26 (4) : 1509-1523. doi: 10.3934/dcds.2010.26.1509 [11] Tomáš Roubíček. Cahn-Hilliard equation with capillarity in actual deforming configurations. Discrete & Continuous Dynamical Systems - S, 2021, 14 (1) : 41-55. doi: 10.3934/dcdss.2020303 [12] Dimitra Antonopoulou, Georgia Karali, Georgios T. Kossioris. Asymptotics for a generalized Cahn-Hilliard equation with forcing terms. Discrete & Continuous Dynamical Systems - A, 2011, 30 (4) : 1037-1054. doi: 10.3934/dcds.2011.30.1037 [13] Annalisa Iuorio, Stefano Melchionna. Long-time behavior of a nonlocal Cahn-Hilliard equation with reaction. Discrete & Continuous Dynamical Systems - A, 2018, 38 (8) : 3765-3788. doi: 10.3934/dcds.2018163 [14] Dimitra Antonopoulou, Georgia Karali. Existence of solution for a generalized stochastic Cahn-Hilliard equation on convex domains. Discrete & Continuous Dynamical Systems - B, 2011, 16 (1) : 31-55. doi: 10.3934/dcdsb.2011.16.31 [15] Alain Miranville, Sergey Zelik. The Cahn-Hilliard equation with singular potentials and dynamic boundary conditions. Discrete & Continuous Dynamical Systems - A, 2010, 28 (1) : 275-310. doi: 10.3934/dcds.2010.28.275 [16] S. Maier-Paape, Ulrich Miller. Connecting continua and curves of equilibria of the Cahn-Hilliard equation on the square. Discrete & Continuous Dynamical Systems - A, 2006, 15 (4) : 1137-1153. doi: 10.3934/dcds.2006.15.1137 [17] Hussein Fakih, Ragheb Mghames, Noura Nasreddine. On the Cahn-Hilliard equation with mass source for biological applications. Communications on Pure & Applied Analysis, 2021, 20 (2) : 495-510. doi: 10.3934/cpaa.2020277 [18] Amy Novick-Cohen, Andrey Shishkov. Upper bounds for coarsening for the degenerate Cahn-Hilliard equation. Discrete & Continuous Dynamical Systems - A, 2009, 25 (1) : 251-272. doi: 10.3934/dcds.2009.25.251 [19] Gianni Gilardi, A. Miranville, Giulio Schimperna. On the Cahn-Hilliard equation with irregular potentials and dynamic boundary conditions. Communications on Pure & Applied Analysis, 2009, 8 (3) : 881-912. doi: 10.3934/cpaa.2009.8.881 [20] Ivonne Rivas, Muhammad Usman, Bing-Yu Zhang. Global well-posedness and asymptotic behavior of a class of initial-boundary-value problem of the Korteweg-De Vries equation on a finite domain. Mathematical Control & Related Fields, 2011, 1 (1) : 61-81. doi: 10.3934/mcrf.2011.1.61

2019 Impact Factor: 1.233