A posteriori error estimation for the Stokes–Darcy coupled problem on anisotropic discretization

This paper presents an a posteriori error analysis for the stationary Stokes–Darcy coupled problem approximated by finite element methods on anisotropic meshes in or 3. Korn's inequality for piecewise linear vector fields on anisotropic meshes is established and is applied to non‐conforming finite element method. Then the existence and uniqueness of the approximation solution are deduced for non‐conforming case. With the obtained finite element solutions, the error estimators are constructed and based on the residual of model equations plus the stabilization terms. The lower error bound is proved by means of bubble functions and the corresponding anisotropic inverse inequalities. In order to prove the upper error bound, it is vital that an anisotropic mesh corresponds to the anisotropic function under consideration. To measure this correspondence, a so‐called matching function is defined, and its discussion shows it to be useful tool. With its help, the upper error bound is shown by means of the corresponding anisotropic interpolation estimates and a special Helmholtz decomposition in both media. Copyright © 2016 John Wiley & Sons, Ltd.     

Modified two-grid method for solving coupled Navier-Stokes/Darcy model based on Newton iteration

A new decoupled two-gird algorithm with the Newton iteration is proposed for solving the coupled Navier-Stokes/Darcy model which describes a fluid flow filtrating through porous media. Moreover the err...     

Integral representation of a solution to the Stokes–Darcy problem

With methods of potential theory, we develop a representation of a solution of the coupled Stokes–Darcy model in a Lipschitz domain for boundary data in H^?1/2. Copyright © 2014 John Wiley & Sons, Ltd.     

Superconvergence analysis of FEMs for the Stokes–Darcy system

We consider a superconvergence analysis for quadratic finite element approximations of the Stokes–Darcy system. The superclose property of an extra half order is proven for uniform triangular meshes. Based on the result of the superclose property, global superconvergence is derived by applying a postprocessing technique. In addition, some numerical examples are presented to demonstrate our theoretical results. Copyright © 2010 John Wiley & Sons, Ltd.     

Finite element analysis of thermally coupled nonlinear Darcy flows

We consider a coupled system describing nonlinear Darcy flows with temperature dependent viscosity and with viscous heating. We first establish existence, uniqueness, and regularity of the weak solution of the system of equations. Next, we decouple the coupled system by a fixed point algorithm and propose its finite element approximation. Finally, we present convergence analysis with an error estimate between continuous solution and its iterative finite element approximation.© 2009 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2010     

Coupling nonlinear Stokes and Darcy flow using mortar finite elements

We study a system composed of a nonlinear Stokes flow in one subdomain coupled with a nonlinear porous medium flow in another subdomain. Special attention is paid to the mathematical consequence of the shear-dependent fluid viscosity for the Stokes flow and the velocity-dependent effective viscosity for the Darcy flow. Motivated by the physical setting, we consider the case where only flow rates are specified on the inflow and outflow boundaries in both subdomains. We recast the coupled Stokes–Darcy system as a reduced matching problem on the interface using a mortar space approach. We prove a number of properties of the nonlinear interface operator associated with the reduced problem, which directly yield the existence, uniqueness and regularity of a variational solution to the system. We further propose and analyze a numerical algorithm based on mortar finite elements for the interface problem and conforming finite elements for the subdomain problems. Optimal a priori error estimates are established for the interface and subdomain problems, and a number of compatibility conditions for the finite element spaces used are discussed. Numerical simulations are presented to illustrate the algorithm and to compare two treatments of the defective boundary conditions.     

Domain decomposition for an asymptotic geological fault modeling

We present in this paper a domain decomposition method to treat faults in geological basin modeling. The particularity of this model is that the faults whose widths are very small in comparison with the basin size, are not characterized as subdomains any more but as interfaces between sedimentary blocks. The originality of this work lies in the formulation of this new fault model and in the definition and the computation of the interface conditions between the subdomains. To cite this article: E. Flauraud et al., C. R. Mecanique 331 (2003).     

Singular nature of nonlinear macroscale effects in high-rate flow through porous mediaNature singulière des effets macroscopiques non linéaires dans un écoulement à hautes vitesses en milieux poreux

The quadratic law of laminar flow through porous media at high Reynolds numbers, which is well confirmed by the multiple experimental data, is shown to give rise to three fundamental paradoxes. All them can be resolved by assuming the singular structure of flow. The singularity is produced by the formation of jet brunches which invade the stagnant zones and sharply loss their kinetic energy. The numerical simulation confirms this effect. To cite this article: M. Panfilov et al., C. R. Mecanique 331 (2003).     

沉降分布孔隙率多孔介质模型及数值研究

研究沉降分布孔隙率多孔介质流动和传热,根据"O"形圈理论和现场测定确定孔隙率系数,建立坐标方向孔隙率分布函数;考虑流体密度变化,并引入Brinkman-Forchheimer的扩展Darcy模型,能量方程采用界面连续条件,建立沉降分布孔隙率多孔介质流动和传热求解模型.采用差分法对模型进行离散化,应用高斯-赛德尔方法迭代求解.数值分析表明:沉降分布孔隙率条件下多孔介质内流体流动速度在壁面附近较大,中心部位较小,壁面附近孔隙率的增大使得低流速区域减小,较高流速区域增大;当孔隙率小值时,温度按线性减小;当孔隙率大值时,温度在高低温壁面附近迅速减小,在中部减小较缓,热量按导热和对流共同传递;孔隙率增大能使平均怒谢尔数增大,对流换热作用增强.     

Poroacoustic solitary waves under the unidirectional Darcy–Jordan model

Working in the context of poroacoustics, we present new, physically relevant, explicit solutions to the Cauchy problem for the model we term the (1D) damped Riemann equation. The solitary waveforms that evolve from both Lorentzian ($C^{\infty }$-smooth) and symmetric-exponential ($C^0$-smooth) initial conditions are analyzed, the focus being on wave overturning and the evolution/structure of the shocks which develop thereafter. In addition to those for both the multi- and single-valued forms of each solution, expressions for the shock amplitude, velocity, and critical values of the physical parameters are derived/compared. Lastly, links to other areas of continuum physics, and possible follow-on investigations, are noted.