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.


Mathematics and Statistics

Research Center/Lab(s)

Center for High Performance Computing Research

Keywords and Phrases

Artificial Compressibility; Ensemble Method; Interface Conditions; Random Hydraulic Conductivity; Stokes-Darcy Model

International Standard Serial Number (ISSN)

0029-5981; 1097-0207

Document Type

Article - Journal

Document Version


File Type





© 2019 John Wiley & Sons Ltd, All rights reserved.

Publication Date

28 Feb 2020