Decoupling methods for the time-dependent Navier-Stokes-Darcy interface model

In this research, several decoupling methods are developed and analyzed for approximating the solution of time-dependent Navier-Stokes-Darcy (NS-Darcy) interface problems. This research on decoupling methods is motivated to efficiently solve the complex Stokes-Darcy or NS-Darcy type models, which arise from many interesting real world problems involved with or even dominated by the coupled porous media flow and free flow. We first discuss a semi-implicit, multi-step non-iterative domain decomposition (NIDDM) to solve a coupled unsteady NS-Darcy system with Beavers-Joseph-Saffman-Jones (BJSJ) interface condition and obtain optimal error estimates. Second, a parallel NIDDM is developed to solve unsteady NS-Darcy model with Beavers-Joseph (BJ) interface condition, which is much more complicated than BJSJ interface condition. We overcome the major difficulties in the analysis which arise from nonlinear terms and BJ interface condition. Furthermore, a Lagrange multiplier method is proposed under the framework of the domain decomposition method to overcome the difficulty of non-unique solutions arising from the defective boundary condition. Meanwhile, we propose and analyze an efficient ensemble algorithm, which can significantly improve the computational efficiency, for fast computation of multiple realizations of the stochastic Stokes-Darcy model with a random hydraulic conductivity tensor. Furthermore, we utilize the idea of artificial compressibility, which decouples the velocity and pressure, to construct the decoupled ensemble algorithm to improve computational efficiency further. We prove that the proposed ensemble methods offer long time stability and optimal error estimates under a time-step condition and two parameter conditions --Abstract, page iii

The LifeV library: engineering mathematics beyond the proof of concept

LifeV is a library for the finite element (FE) solution of partial differential equations in one, two, and three dimensions. It is written in C++ and designed to run on diverse parallel architectures, including cloud and high performance computing facilities. In spite of its academic research nature, meaning a library for the development and testing of new methods, one distinguishing feature of LifeV is its use on real world problems and it is intended to provide a tool for many engineering applications. It has been actually used in computational hemodynamics, including cardiac mechanics and fluid-structure interaction problems, in porous media, ice sheets dynamics for both forward and inverse problems. In this paper we give a short overview of the features of LifeV and its coding paradigms on simple problems. The main focus is on the parallel environment which is mainly driven by domain decomposition methods and based on external libraries such as MPI, the Trilinos project, HDF5 and ParMetis. Dedicated to the memory of Fausto Saleri.Comment: Review of the LifeV Finite Element librar

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 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

Optimized Schwarz Methods in the Stokes-Darcy Coupling

This article studies optimized Schwarz methods for the Stokes–Darcy problem. Robin transmission conditions are introduced, and the coupled problem is reduced to a suitable interface system that can be solved using Krylov methods. Practical strategies to compute optimal Robin coefficients are proposed, which take into account both the physical parameters of the problem and the mesh size. Numerical results show the effectiveness of our approach.European Union Seventh Framework Programme (FP7/2007-2013; grant 294229) to M. Discacciat

Reactive Flows in Deformable, Complex Media

Many processes of highest actuality in the real life are described through systems of equations posed in complex domains. Of particular interest is the situation when the domain is variable, undergoing deformations that depend on the unknown quantities of the model. Such kind of problems are encountered as mathematical models in the subsurface, or biological systems. Such models include various processes at different scales, and the key issue is to integrate the domain deformation in the multi-scale context. Having this as the background theme, this workshop focused on novel techniques and ideas in the analysis, the numerical discretization and the upscaling of such problems, as well as on applications of major societal relevance today

Optimized Schwarz methods for the Stokes-Darcy coupling

This paper studies Optimized Schwarz methods for the Stokes-Darcy problem. Robin transmission conditions are introduced and the coupled problem is reduced to a suitable interface system that can be solved using Krylov methods. Practical strategies to compute optimal Robin coefficients are proposed which take into account both the physical parameters of the problem and the mesh size. Numerical results show the effectiveness of our approach

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