A Domain Decomposition Method for the Steady-State Navier-Stokes-Darcy Model with Beavers-Joseph Interface Condition

This paper proposes and analyzes a Robin-type multiphysics domain decomposition method (DDM) for the steady-state Navier-Stokes-Darcy model with three interface conditions. In addition to the two regular interface conditions for the mass conservation and the force balance, the Beavers-Joseph condition is used as the interface condition in the tangential direction. The major mathematical difficulty in adopting the Beavers-Joseph condition is that it creates an indefinite leading order contribution to the total energy budget of the system [Y. Cao et al., Comm. Math. Sci., 8 (2010), pp. 1-25; Y. Cao et al., SIAM J. Numer. Anal., 47 (2010), pp. 4239-4256]. In this paper, the well-posedness of the Navier-Stokes-Darcy model with Beavers-Joseph condition is analyzed by using a branch of nonsingular solutions. By following the idea in [Y. Cao et al., Numer. Math., 117 (2011), pp. 601-629], the three physical interface conditions are utilized together to construct the Robin-type boundary conditions on the interface and decouple the two physics which are described by Navier-Stokes and Darcy equations, respectively. Then the corresponding multiphysics DDM is proposed and analyzed. Three numerical experiments using finite elements are presented to illustrate the features of the proposed method and verify the results of the theoretical analysis

Modeling and a Domain Decomposition Method with Finite Element Discretization for Coupled Dual-Porosity Flow and Navier–Stokes Flow

In This Paper, We First Propose and Analyze a Steady State Dual-Porosity-Navier–Stokes Model, Which Describes Both Dual-Porosity Flow and Free Flow (Governed by Navier–Stokes Equation) Coupled through Four Interface Conditions, Including the Beavers–Joseph Interface Condition. Then We Propose a Domain Decomposition Method for Efficiently Solving Such a Large Complex System. Robin Boundary Conditions Are Used to Decouple the Dual-Porosity Equations from the Navier–Stokes Equations in the Coupled System. based on the Two Decoupled Sub-Problems, a Parallel Robin-Robin Domain Decomposition Method is Constructed and Then Discretized by Finite Elements. We Analyze the Convergence of the Domain Decomposition Method with the Finite Element Discretization and Investigate the Effect of Robin Parameters on the Convergence, Which Also Provide Instructions for How to Choose the Robin Parameters in Practice. Three Cases of Robin Parameters Are Studied, Including a Difficult Case Which Was Not Fully Addressed in the Literature, and the Optimal Geometric Convergence Rate is Obtained. Numerical Experiments Are Presented to Verify the Theoretical Conclusions, Illustrate How the Theory Can Provide Instructions on Choosing Robin Parameters, and Show the Features of the Proposed Model and Domain Decomposition Method

Decoupling the Stationary Navier-Stokes-Darcy System with the Beavers-Joseph-Saffman Interface Condition

This paper proposes a domain decomposition method for the coupled stationary Navier-Stokes and Darcy equations with the Beavers-Joseph-Saffman interface condition in order to improve the efficiency of the finite element method. The physical interface conditions are directly utilized to construct the boundary conditions on the interface and then decouple the Navier-Stokes and Darcy equations. Newton iteration will be used to deal with the nonlinear systems. Numerical results are presented to illustrate the features of the proposed method

A Parallel Robin-Robin Domain Decomposition Method for the Stokes-Darcy System

We propose a new parallel Robin-Robin domain decomposition method for the coupled Stokes-Darcy system with Beavers-Joseph-Saffman-Jones interface boundary condition. in particular, we prove that, with an appropriate choice of parameters, the scheme converges geometrically independent of the mesh size. © 2011 Society for Industrial and Applied Mathematics

On iterative subdomain methods for the Stokes--Darcy problem

Iterative subdomain methods for the Stokes--Darcy problem that use Robin boundary conditions on the interface are reviewed. Their common underlying structure and their main differences are identified. In particular, it is clarified that there are different updating strategies for the interface conditions. For small values of fluid viscosity and hydraulic permeability, which are relevant in applications from geosciences, it is shown in numerical studies that only one of these updating strategies leads to an efficient numerical method, if this strategy is used in combination with appropriate parameters in the Robin boundary conditions. In particular, it is observed that the values of appropriate parameters are larger than those proposed so far. Not only the size but also the ratio of appropriate Robin parameters depends on the coefficients of the problem

Coupling Biot and Navier-Stokes equations for modelling fluid-poroelastic media interaction

The interaction between a fluid and a poroelastic structure is a complex problem that couples the Navier–Stokes equations with the Biot system. The finite element approximation of this problem is involved due to the fact that both subproblems are indefinite. In this work, we first design residual-based stabilization techniques for the Biot system, motivated by the variational multiscale approach. Then, we state the monolithic Navier–Stokes/Biot system with the appropriate transmission conditions at the interface. For the solution of the coupled system, we adopt both monolithic solvers and heterogeneous domain decomposition strategies. Different domain decomposition methods are considered and their convergence is analyzed for a simplified problem. We compare the efficiency of all the methods on a test problem that exhibits a large added-mass effect, as it happens in hemodynamics applications

A Multigrid Multilevel Monte Carlo Method for Stokes–Darcy Model with Random Hydraulic Conductivity and Beavers–Joseph Condition

A multigrid multilevel Monte Carlo (MGMLMC) method is developed for the stochastic Stokes–Darcy interface model with random hydraulic conductivity both in the porous media domain and on the interface. Three interface conditions with randomness are considered on the interface between Stokes and Darcy equations, especially the Beavers–Joesph interface condition with random hydraulic conductivity. Because the randomness through the interface affects the flow in the Stokes domain, we investigate the coupled stochastic Stokes–Darcy model to improve the fidelity. Under suitable assumptions on the random coefficient, we prove the existence and uniqueness of the weak solution of the variational form. To construct the numerical method, we first adopt the Monte Carlo (MC) method and finite element method, for the discretization in the probability space and physical space, respectively. In order to improve the efficiency of the classical single-level Monte Carlo (SLMC) method, we adopt the multilevel Monte Carlo (MLMC) method to dramatically reduce the computational cost in the probability space. A strategy is developed to calculate the number of samples needed in MLMC method for the stochastic Stokes–Darcy model. In order to accomplish the strategy for MLMC method, we also present a practical method to determine the variance convergence rate for the stochastic Stokes–Darcy model with Beavers–Joseph interface condition. Furthermore, MLMC method naturally provides the hierarchical grids and sufficient information on these grids for multigrid (MG) method, which can in turn improve the efficiency of MLMC method. In order to fully make use of the dynamical interaction between these two methods, we propose a multigrid multilevel Monte Carlo (MGMLMC) method with finite element discretization for more efficiently solving the stochastic model, while additional attention is paid to the interface and the random Beavers–Joesph interface condition. The computational cost of the proposed MGMLMC method is rigorously analyzed and compared with the SLMC method. Numerical examples are provided to verify and illustrate the proposed method and the theoretical conclusions

modelling fluid-poroelastic media interaction

The interaction between a fluid and a poroelastic structure is a complex problem that couples the Navier-Stokes equations with the Biot system. The finite element approximation of this problem is involved due to the fact that both subproblems are indefinite. In this work, we first design residual-based stabilization techniques for the Biot system, motivated by the variational multiscale approach. Then, we state the monolithic Navier-Stokes/Biot system with the appropriate transmission conditions at the interface. For the solution of the coupled system, we adopt both monolithic solvers and heterogeneous domain decomposition strategies. Different domain decomposition methods are considered and their convergence 1 is analyzed for a simplified problem. We compare the efficiency of all the methods on a test problem that exhibits a large added-mass effect, as it happens in hemodynamics applications.