Ensemble Domain Decomposition Algorithm for the Fully Mixed Random Stokes–Darcy Model with the Beavers–Joseph Interface Conditions
Feng Shi,Yizhong Sun,Haibiao Zheng
DOI: https://doi.org/10.1137/22m1482846
IF: 3.039
2023-06-14
SIAM Journal on Numerical Analysis
Abstract:SIAM Journal on Numerical Analysis, Volume 61, Issue 3, Page 1482-1512, June 2023. In this paper, an efficient ensemble domain decomposition algorithm is proposed for fast solving the fully mixed random Stokes–Darcy model with the physically realistic Beavers–Joseph interface conditions. We utilize the Monte Carlo method for the coupled model with random inputs to derive some deterministic Stokes–Darcy numerical models and use the idea of the ensemble to realize the fast computation of multiple problems. One remarkable feature of the algorithm is that multiple linear systems share a common coefficient matrix in each deterministic numerical model, which significantly reduces the computational cost and achieves comparable accuracy with the traditional methods. Moreover, by domain decomposition, we can decouple the Stokes–Darcy system into two smaller subphysics problems naturally. Both mesh-dependent and mesh-independent convergence rates of the algorithm are rigorously derived by choosing suitable Robin parameters. Optimized Robin parameters are derived and analyzed to accelerate the convergence of the proposed algorithm. Especially, for small hydraulic conductivity in practice, the almost optimal geometric convergence can be obtained by finite element discretization. Finally, two groups of numerical experiments are conducted to validate the exclusive features of the proposed algorithm.
mathematics, applied