Stabilized Variational Formulations of Chorin-Type and Artificial Compressibility Methods for the Stochastic Stokes–Darcy Equations

In this paper, we consider two different types of numerical schemes for the nonstationary stochastic Stokes–Darcy equations with multiplicative noise. Firstly, we consider the Chorin-type time-splitting scheme for the Stokes equation in the free fluid region. The Darcy equation and an elliptic equation for the intermediate velocity of free fluid coupled with the interface conditions are solved, and then the velocity and pressure in free fluid region are updated by an elliptic system. Secondly, we further consider the artificial compressibility method (ACM) which separates the fully coupled Stokes–Darcy model into two smaller subphysics problems. The ACM reduces the storage and the computational time at each time step, and allows parallel computing for the decoupled problems. The pressure in free fluid region only needs to be updated at each time step without solving an elliptic system. We utilize the RT \(_1\) $_1^{}$ -P \(_1\) $_1^{}$ pair finite element space and the interior penalty discontinuous Galerkin (IPDG) scheme based on the BDM \(_1\) $_1^{}$ -P \(_0\) $_0^{}$ finite element space in the spatial discretizations. Under usual assumptions for the multiplicative noise, we prove that both of the Chorin-type scheme and the ACM are unconditionally stable. We present the error estimates for the time semi-discretization of the Chorin-type scheme. Numerical examples are provided to verify the stability estimates for both of schemes. Moreover, we test the convergence rate for the velocity in time for both of schemes, and the convergence rate for the pressure approximation in time average is also tested.