We propose and study two second order in time implicit-explicit methods for the coupled Stokes-Darcy system that governs flows in karst aquifers and other subsurface flow systems. the first method is a combination of a second-order backward differentiation formula and the second order Gear's extrapolation approach. the second is a combination of the second-order Adams-Moulton and second-order Adams-Bashforth methods. Both algorithms only require the solution of decoupled Stokes and Darcy problems at each time-step. Hence, these schemes are very efficient and can be easily implemented using legacy codes. We establish the unconditional and uniform in time stability for both schemes. the uniform in time stability leads to uniform in time control of the error which is highly desirable for modeling physical processes, e.g., contaminant sequestration and release, that occur over very long-time scales. Error estimates for fully discretized schemes using finite element spatial discretization's are derived. Numerical examples are provided that illustrate the accuracy, efficiency, and long-time stability of the two schemes. © 2013 Society for Industrial and Applied Mathematics.


Mathematics and Statistics


National Science Foundation, Grant DMS10008852

Keywords and Phrases

Adams-Moulton and Adams-Bashforth Methods; Backward Differentiation Formulas; Finite Element Methods; Gear's Extrapolation; Karst Aquifers; Long-Time Stability; Stokes-Darcy Systems; Unconditional Stability; Uniform in Time Error Estimates

International Standard Serial Number (ISSN)


Document Type

Article - Journal

Document Version

Final Version

File Type





© 2023 Society for Industrial and Applied Mathematics, All rights reserved.

Publication Date

25 Dec 2013