Stability and convergence of the mark and cell finite difference scheme for Darcy‐Stokes‐Brinkman equations on non‐uniform grids

In this paper, we consider the mark and cell (MAC) method for Darcy‐Stokes‐Brinkman equations and analyze the stability and convergence of the method on nonuniform grids. Firstly, to obtain the stability for both velocity and pressure, we establish the discrete inf‐sup condition. Then we introduce an auxiliary function depending on the velocity and discretizing parameters to analyze the super‐convergence. Finally, we obtain the second‐order convergence in L2 norm for both velocity and pressure for the MAC scheme, when the perturbation parameter ? is not approaching 0. We also obtain the second‐order convergence for some terms of ∥·∥_? norm of the velocity, and the other terms of ∥·∥_? norm are second‐order convergence on uniform grid. Numerical experiments are carried out to verify the theoretical results.     

Partitioned penalty methods for the transport equation in the evolutionary Stokes–Darcy‐transport problem

There has been a surge of work on models for coupling surface‐water with groundwater flows which is at its core the Stokes–Darcy problem, as well as methods for uncoupling the problem into subdomain, subphysics solves. The resulting (Stokes–Darcy) fluid velocity is important because the flow transports contaminants. The numerical analysis and algorithm development for the evolutionary transport problem has, however, focused on a quasi‐static Stokes–Darcy model and a single domain (fully coupled) formulation of the transport equation. This report presents a numerical analysis of a partitioned method for contaminant transport for the fully evolutionary system. The algorithm studied is unconditionally stable with one subdomain solve per step. Numerical experiments are given using the proposed algorithm that investigates the effects of the penalty parameters on the convergence of the approximations.     

A conforming mixed finite-element method for the coupling of fluid flow with porous media flow

Salim Meddahi ^{We consider a porous medium entirely enclosed within a fluidregion and present a well-posed conforming mixed finite-elementmethod for the corresponding coupled problem. The interfaceconditions refer to mass conservation, balance of normal forcesand the Beavers–Joseph–Saffman law, which yieldsthe introduction of the trace of the porous medium pressureas a suitable Lagrange multiplier. The finite-element subspacesdefining the discrete formulation employ Bernardi–Raugeland Raviart–Thomas elements for the velocities, piecewiseconstants for the pressures and continuous piecewise-linearelements for the Lagrange multiplier. We show stability, convergenceand a priori error estimates for the associated Galerkin scheme.Finally, we provide several numerical results illustrating thegood performance of the method and confirming the theoreticalrates of convergence.}     

A second‐order partitioned method with different subdomain time steps for the evolutionary Stokes‐Darcy system

A second‐order decoupled algorithm for the nonstationary Stokes‐Darcy system, which allows different time steps in two subregions, is proposed and analyzed in this paper. The algorithm, which is a combination of the second‐order backward differentiation formula and second‐order extrapolation method, uncouples the problem into two decoupled problems per time step. We prove the unconditional stability and long‐time stability of the decoupled scheme with different time steps and derive error estimates of this decoupled algorithm using finite element spatial discretization. Numerical experiments are provided to illustrate the accuracy, effectiveness, and stability of the decoupled algorithm and show its advantages of increasing accuracy and efficiency.     

Darcy–Forchheimer 3D flow of carbon nanotubes with homogeneous and heterogeneous reactions

This article deals with Darcy–Forchheimer three dimensional (3D) flow of water-based carbon nanotubes (CNTs) with heterogeneous–homogeneous reactions. A bidirectional nonlinear extendable surface has been employed to create the flow. Flow in porous space is represented by Darcy–Forchheimer expression. Heat transfer mechanism is explored through convective heating. Equal diffusion coefficients are considered for both auto catalyst and reactants. Results for single-wall (SWCNT) and multi-wall (MWCNT) carbon nanotubes have been presented and compared. The diminishment of partial differential framework into nonlinear ordinary differential framework is made through suitable transformations. Optimal homotopy scheme is used for arrangements development of governing flow problem. Optimal homotopic solution expressions for velocities and temperature are studied through plots by considering various estimations of physical variables. Moreover the surface drag coefficients and heat transfer rate are analyzed through plots.     

Non-Darcy free convection of Fe3O4-water nanoliquid in a complex shaped enclosure under impact of uniform Lorentz force

Effect of Lorentz forces on natural convection in a complex shaped cavity filled with nanoliquid immersed in porous medium is investigated by means of Control volume based finite element method (CVFEM). Non Darcy model is taken into account for porous media. The working fluid is Fe₃O₄ –water and its viscosity considered as function of magnetic field. Figures are illustrated for different values of Darcy number (Da), Fe₃O₄ -water volume fraction (?), Rayleigh (Ra) and Hartmann (Ha) numbers. Results depict that enhancing in Lorentz forces results in reduce in nanofluid motion and increase the thickness of thermal boundary. Convective heat transfer enhances with rise of Darcy number.     

Numerical analysis of the coupling of free fluid with a poroelastic material

This paper presents a numerical solution of the coupled system of the time-dependent Stokes and fully dynamic Biot equations. The numerical scheme is based on standard inf-sup stable finite elements in space and the Backward Euler scheme in time. We establish stability of the scheme and derive error estimates for the fully discrete coupled scheme. To handle realistic parameters which may cause nonphysical oscillations in the pore fluid pressure, a heuristic stabilization technique is considered. Numerical errors and convergence rates for smooth problems as well as tests on realistic material parameters are presented.     

Nonlinear systems arising from nonisothermal, non-Newtonian Hele-Shaw flows in the presence of body forces and sources

In this paper, we give a formal derivation of several systems of equations for injection moulding. This is done starting from the basic equations for nonisothermal, non-Newtonian flows in a three-dimensional domain. We derive systems for both (T⁰, p⁰) and (T¹, p¹) in the presence of body forces and sources. We find that body forces and sources have a nonlinear effect on the systems. We also derive a nonlinear “Darcy law”. Our formulation includes not only the pressure gradient, but also body forces and sources, which play the role of a nonlinearity. Later, we prove the existence of weak solutions to certain boundary value problems and initial-boundary value problems associated with the resulting equations for (T⁰,p⁰) but in a more general mathematical setting.     

Convergence order estimates of the local discontinuous Galerkin method for instationary Darcy flow

We present a new a priori stability and convergence analysis for the local discontinuous Galerkin method applied to the instationary Darcy problem formulated on a d‐dimensional polytope with nonhomogeneous Neumann and Dirichlet boundary conditions. In addition to including a spatially and temporally varying permeability tensor into all estimates, the utilized analysis technique produces a convergence order result for the primary and the flux variables. The only stabilization in the proposed scheme is represented by a penalty term in the primary unknown, and our analysis provides some insights into the role played by this particular choice of stabilization. © 2016 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 33: 1374–1394, 2017     

Spectral discretization of Darcy equations coupled with Stokes equations by vorticity–velocity–pressure formulations

We propose to make the numerical analysis of a model coupling the Darcy equations in a porous medium with the Stokes equations in the cracks. The coupling is provided by a pressure continuity on the interface. We describe a discretization by spectral element methods. We derive a priori optimal error estimates and we present some numerical experiments which confirm the results of the analysis.© 2017 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 33: 1628–1651, 2017