# American Institute of Mathematical Sciences

March  2021, 14(3): 1197-1212. doi: 10.3934/dcdss.2020234

## Stability estimates and a Lagrange-Galerkin scheme for a Navier-Stokes type model of flow in non-homogeneous porous media

 1 Division of Mathematical and Physics Science, Kanazawa University, Kakuma, Kanazawa 920-1192, Japan 2 Faculty of Mathematical and Physics, Kanazawa University, Kakuma, Kanazawa 920-1192, Japan 3 Japan Science and Technology Agency, PRESTO, Kawaguchi 332-0012, Japan

* Corresponding author: Imam Wijaya

Received  January 2019 Revised  May 2019 Published  December 2019

The purposes of this work are to study the $L^{2}$-stability of a Navier-Stokes type model for non-stationary flow in porous media proposed by Hsu and Cheng in 1989 and to develop a Lagrange-Galerkin scheme with the Adams-Bashforth method to solve that model numerically. The stability estimate is obtained thanks to the presence of a nonlinear drag force term in the model which corresponds to the Forchheimer term. We derive the Lagrange-Galerkin scheme by extending the idea of the method of characteristics to overcome the difficulty which comes from the non-homogeneous porosity. Numerical experiments are conducted to investigate the experimental order of convergence of the scheme. For both simple and complex designs of porosities, our numerical simulations exhibit natural flow profiles which well describe the flow in non-homogeneous porous media.

Citation: Imam Wijaya, Hirofumi Notsu. Stability estimates and a Lagrange-Galerkin scheme for a Navier-Stokes type model of flow in non-homogeneous porous media. Discrete & Continuous Dynamical Systems - S, 2021, 14 (3) : 1197-1212. doi: 10.3934/dcdss.2020234
##### References:

show all references

##### References:
Representative elementary volume (REV)
The order of convergence for scheme (22)
The boundary conditions and the finite element mesh
Time evolution of velocity magnitude
Computation domain and porosity value distribution
Time evolution of magnitude velocity
The unit of important symbols
 No Symbol Unit Name of the symbol [0.5ex] 1 $u$ ${\rm m} \cdot {\rm s}^{-1}$ Darcy velocity 2 $p$ ${\rm kg} \cdot {\rm m}^{-1}\cdot {\rm s}^{-2}$ Pressure 3 $\phi$ - porosity 4 $k_{D}$ ${\rm kg}^{-1} \cdot {\rm m}^3 \cdot {\rm s}$ Hydraulic conductivity 5 $K$ ${\rm m}^2$ Permeability 6 $\mu$ ${\rm kg} \cdot {\rm m}^{-1}\cdot {\rm s}^{-1}$ Dynamic viscosity 7 $\rho$ ${\rm kg} \cdot {\rm m}^{-3}$ Density 8 $d_{p}$ ${\rm m}$ Particle diameter 9 $F$ - Forchheimer constant 10 $B$ ${\rm kg} \cdot {\rm m}^{-2} \cdot {\rm s}^{-2}$ Drag force per unit volume
 No Symbol Unit Name of the symbol [0.5ex] 1 $u$ ${\rm m} \cdot {\rm s}^{-1}$ Darcy velocity 2 $p$ ${\rm kg} \cdot {\rm m}^{-1}\cdot {\rm s}^{-2}$ Pressure 3 $\phi$ - porosity 4 $k_{D}$ ${\rm kg}^{-1} \cdot {\rm m}^3 \cdot {\rm s}$ Hydraulic conductivity 5 $K$ ${\rm m}^2$ Permeability 6 $\mu$ ${\rm kg} \cdot {\rm m}^{-1}\cdot {\rm s}^{-1}$ Dynamic viscosity 7 $\rho$ ${\rm kg} \cdot {\rm m}^{-3}$ Density 8 $d_{p}$ ${\rm m}$ Particle diameter 9 $F$ - Forchheimer constant 10 $B$ ${\rm kg} \cdot {\rm m}^{-2} \cdot {\rm s}^{-2}$ Drag force per unit volume
Values of $Er1$ and $Er2$, their slopes, and CPU times for the problem in Subsection 6.1 by scheme (22)
 $N$ $Er1$ $Er2$ Slope of $Er1$ Slope of $Er2$ CPU time [s] 4 $3.4\times10^{-1}$ $1.6\times10^{-1}$ $-$ $-$ 1.9 8 $7.1\times10^{-2}$ $5.8\times10^{-3}$ 2.26 4.76 16.4 16 $1.4\times10^{-2}$ $1.2\times10^{-3}$ 2.34 2.30 174.8 32 $3.5\times10^{-3}$ $2.9\times10^{-4}$ 2.00 2.05 577.4 64 $1.0\times10^{-3}$ $6.3\times10^{-5}$ 1.81 2.20 5,953.9 128 $2.8\times10^{-4}$ $1.5\times10^{-5}$ 1.84 2.07 58,150.9
 $N$ $Er1$ $Er2$ Slope of $Er1$ Slope of $Er2$ CPU time [s] 4 $3.4\times10^{-1}$ $1.6\times10^{-1}$ $-$ $-$ 1.9 8 $7.1\times10^{-2}$ $5.8\times10^{-3}$ 2.26 4.76 16.4 16 $1.4\times10^{-2}$ $1.2\times10^{-3}$ 2.34 2.30 174.8 32 $3.5\times10^{-3}$ $2.9\times10^{-4}$ 2.00 2.05 577.4 64 $1.0\times10^{-3}$ $6.3\times10^{-5}$ 1.81 2.20 5,953.9 128 $2.8\times10^{-4}$ $1.5\times10^{-5}$ 1.84 2.07 58,150.9
 [1] Yueqiang Shang, Qihui Zhang. A subgrid stabilizing postprocessed mixed finite element method for the time-dependent Navier-Stokes equations. Discrete & Continuous Dynamical Systems - B, 2021, 26 (6) : 3119-3142. doi: 10.3934/dcdsb.2020222 [2] Ling-Bing He, Li Xu. On the compressible Navier-Stokes equations in the whole space: From non-isentropic flow to isentropic flow. Discrete & Continuous Dynamical Systems, 2021, 41 (7) : 3489-3530. doi: 10.3934/dcds.2021005 [3] Daoyuan Fang, Ting Zhang. Compressible Navier-Stokes equations with vacuum state in one dimension. Communications on Pure & Applied Analysis, 2004, 3 (4) : 675-694. doi: 10.3934/cpaa.2004.3.675 [4] Thomas Y. Hou, Ruo Li. Nonexistence of locally self-similar blow-up for the 3D incompressible Navier-Stokes equations. Discrete & Continuous Dynamical Systems, 2007, 18 (4) : 637-642. doi: 10.3934/dcds.2007.18.637 [5] Cheng Wang. Convergence analysis of Fourier pseudo-spectral schemes for three-dimensional incompressible Navier-Stokes equations. Electronic Research Archive, , () : -. doi: 10.3934/era.2021019 [6] Xin-Guang Yang, Rong-Nian Wang, Xingjie Yan, Alain Miranville. Dynamics of the 2D Navier-Stokes equations with sublinear operators in Lipschitz-like domains. Discrete & Continuous Dynamical Systems, 2021, 41 (7) : 3343-3366. doi: 10.3934/dcds.2020408 [7] Xin Zhong. Global strong solution and exponential decay for nonhomogeneous Navier-Stokes and magnetohydrodynamic equations. Discrete & Continuous Dynamical Systems - B, 2021, 26 (7) : 3563-3578. doi: 10.3934/dcdsb.2020246 [8] Huancheng Yao, Haiyan Yin, Changjiang Zhu. Stability of rarefaction wave for the compressible non-isentropic Navier-Stokes-Maxwell equations. Communications on Pure & Applied Analysis, 2021, 20 (3) : 1297-1317. doi: 10.3934/cpaa.2021021 [9] Mirela Kohr, Sergey E. Mikhailov, Wolfgang L. Wendland. Dirichlet and transmission problems for anisotropic stokes and Navier-Stokes systems with L∞ tensor coefficient under relaxed ellipticity condition. Discrete & Continuous Dynamical Systems, 2021  doi: 10.3934/dcds.2021042 [10] Francis Hounkpe, Gregory Seregin. An approximation of forward self-similar solutions to the 3D Navier-Stokes system. Discrete & Continuous Dynamical Systems, 2021  doi: 10.3934/dcds.2021059 [11] Mikhail Dokuchaev, Guanglu Zhou, Song Wang. A modification of Galerkin's method for option pricing. Journal of Industrial & Management Optimization, 2021  doi: 10.3934/jimo.2021077 [12] Derrick Jones, Xu Zhang. A conforming-nonconforming mixed immersed finite element method for unsteady Stokes equations with moving interfaces. Electronic Research Archive, , () : -. doi: 10.3934/era.2021032 [13] Matthias Erbar, Jan Maas. Gradient flow structures for discrete porous medium equations. Discrete & Continuous Dynamical Systems, 2014, 34 (4) : 1355-1374. doi: 10.3934/dcds.2014.34.1355 [14] Armin Lechleiter, Tobias Rienmüller. Factorization method for the inverse Stokes problem. Inverse Problems & Imaging, 2013, 7 (4) : 1271-1293. doi: 10.3934/ipi.2013.7.1271 [15] Tomáš Roubíček. An energy-conserving time-discretisation scheme for poroelastic media with phase-field fracture emitting waves and heat. Discrete & Continuous Dynamical Systems - S, 2017, 10 (4) : 867-893. doi: 10.3934/dcdss.2017044 [16] Nouressadat Touafek, Durhasan Turgut Tollu, Youssouf Akrour. On a general homogeneous three-dimensional system of difference equations. Electronic Research Archive, , () : -. doi: 10.3934/era.2021017 [17] Zhiming Guo, Zhi-Chun Yang, Xingfu Zou. Existence and uniqueness of positive solution to a non-local differential equation with homogeneous Dirichlet boundary condition---A non-monotone case. Communications on Pure & Applied Analysis, 2012, 11 (5) : 1825-1838. doi: 10.3934/cpaa.2012.11.1825 [18] Lunji Song, Wenya Qi, Kaifang Liu, Qingxian Gu. A new over-penalized weak galerkin finite element method. Part Ⅱ: Elliptic interface problems. Discrete & Continuous Dynamical Systems - B, 2021, 26 (5) : 2581-2598. doi: 10.3934/dcdsb.2020196 [19] Kaifang Liu, Lunji Song, Shan Zhao. A new over-penalized weak galerkin method. Part Ⅰ: Second-order elliptic problems. Discrete & Continuous Dynamical Systems - B, 2021, 26 (5) : 2411-2428. doi: 10.3934/dcdsb.2020184 [20] Abderrazak Chrifi, Mostafa Abounouh, Hassan Al Moatassime. Galerkin method of weakly damped cubic nonlinear Schrödinger with Dirac impurity, and artificial boundary condition in a half-line. Discrete & Continuous Dynamical Systems - S, 2021  doi: 10.3934/dcdss.2021030

2019 Impact Factor: 1.233