skip to main content

Title: An artificial compressibility ensemble algorithm for a stochastic Stokes‐Darcy model with random hydraulic conductivity and interface conditions

We propose and analyze an efficient ensemble algorithm with artificial compressibility (AC) for fast decoupled computation of multiple realizations of the stochastic Stokes‐Darcy model with random hydraulic conductivity (including the one in the interface conditions), source terms, and initial conditions. The solutions are found by solving three smaller decoupled subproblems with two common time‐independent coefficient matrices for all realizations, which significantly improves the efficiency for both assembling and solving the matrix systems. The fully coupled Stokes‐Darcy system can be first decoupled into two smaller subphysics problems by the idea of the partitioned time stepping, which reduces the size of the linear systems and allows parallel computing for each subphysics problem. The AC further decouples the velocity and pressure which further reduces storage requirements and improves computational efficiency. We prove the long time stability and the convergence for this new ensemble method. Three numerical examples are presented to support the theoretical results and illustrate the features of the algorithm, including the convergence, stability, efficiency, and applicability.

more » « less
Award ID(s):
1720001 1722647 2120413
Author(s) / Creator(s):
 ;  ;  
Publisher / Repository:
Wiley Blackwell (John Wiley & Sons)
Date Published:
Journal Name:
International Journal for Numerical Methods in Engineering
Page Range / eLocation ID:
p. 712-739
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Summary

    In this paper, we propose and analyze two stabilized mixed finite element methods for the dual‐porosity‐Stokes model, which couples the free flow region and microfracture‐matrix system through four interface conditions on an interface. The first stabilized mixed finite element method is a coupled method in the traditional format. Based on the idea of partitioned time stepping, the four interface conditions, and the mass exchange terms in the dual‐porosity model, the second stabilized mixed finite element method is decoupled in two levels and allows a noniterative splitting of the coupled problem into three subproblems. Due to their superior conservation properties and convenience of the computation of flux, mixed finite element methods have been widely developed for different types of subsurface flow problems in porous media. For the mixed finite element methods developed in this article, no Lagrange multiplier is used, but an interface stabilization term with a penalty parameter is added in the temporal discretization. This stabilization term ensures the numerical stability of both the coupled and decoupled schemes. The stability and the convergence analysis are carried out for both the coupled and decoupled schemes. Three numerical experiments are provided to demonstrate the accuracy, efficiency, and applicability of the proposed methods.

    more » « less
  2. Studying the propagation of uncertainties in a nonlinear dynamical system usually involves generating a set of samples in the stochastic parameter space and then repeated simulations with different sampled parameters. The main difficulty faced in the process is the excessive computational cost. In this paper, we present an efficient, partitioned ensemble algorithm to determine multiple realizations of a reduced Magnetohydrodynamics (MHD) system, which models MHD flows at low magnetic Reynolds number. The algorithm decouples the fully coupled problem into two smaller subphysics problems, which reduces the size of the linear systems that to be solved and allows the use of optimized codes for each subphysics problem. Moreover, the resulting coefficient matrices are the same for all realizations at each time step, which allows faster computation of all realizations and significant savings in computational cost. We prove this algorithm is first order accurate and long time stable under a time step condition. Numerical examples are provided to verify the theoretical results and demonstrate the efficiency of the algorithm.

    more » « less
  3. Abstract

    We study a second order ensemble method for fast computation of an ensemble of magnetohydrodynamics flows at small magnetic Reynolds number. Computing an ensemble of flow equations with different input parameters is a common procedure for uncertainty quantification in many engineering applications, for which the computational cost can be prohibitively expensive for nonlinear complex systems. We propose an ensemble algorithm that requires only solving one linear system with multiple right‐hands instead of solving multiple different linear systems, which significantly reduces the computational cost and simulation time. Comprehensive stability and error analyses are presented proving conditional stability and second order in time convergent. Numerical tests are provided to illustrate theoretical results and demonstrate the efficiency of the proposed algorithm.

    more » « less
  4. In this article, we consider a phase field model with different densities and viscosities for the coupled two-phase porous media flow and two-phase free flow, as well as the corresponding numerical simulation. This model consists of three parts: a Cahn–Hilliard–Darcy system with different densities/viscosities describing the porous media flow in matrix, a Cahn–Hilliard–Navier–Stokes system with different densities/viscosities describing the free fluid in conduit, and seven interface conditions coupling the flows in the matrix and the conduit. Based on the separate Cahn–Hilliard equations in the porous media region and the free flow region, a weak formulation is proposed to incorporate the two-phase systems of the two regions and the seven interface conditions between them, and the corresponding energy law is proved for the model. A fully decoupled numerical scheme, including the novel decoupling of the Cahn–Hilliard equations through the four phase interface conditions, is developed to solve this coupled nonlinear phase field model. An energy-law preservation is analyzed for the temporal semi-discretization scheme. Furthermore, a fully discretized Galerkin finite element method is proposed. Six numerical examples are provided to demonstrate the accuracy, discrete energy law, and applicability of the proposed fully decoupled scheme. 
    more » « less
  5. We consider a space-time domain decomposition method based on Schwarz waveform relaxation (SWR) for the time-dependent Stokes-Darcy system. The coupled system is formulated as a time-dependent interface problem based on Robin-Robin transmission conditions, for which the decoupling SWR algorithm is proposed and proved for the convergence. In this approach, the Stokes and Darcy problems are solved independently and globally in time, thus allowing the use of different time steps for the local problems. Numerical tests are presented for both non-physical and physical problems with various mesh sizes and time step sizes to illustrate the accuracy and efficiency of the proposed method. 
    more » « less