An Efficient and Long-Time Accurate Third-Order Algorithm for the Stokes–Darcy System

A third order in time numerical IMEX-type algorithm for the Stokes–Darcy system for flows in fluid saturated karst aquifers is proposed and analyzed. a novel third-order Adams–Moulton scheme is used for the discretization of the dissipative term whereas a third-order explicit Adams–Bashforth scheme is used for the time discretization of the interface term that couples the Stokes and Darcy components. the scheme is efficient in the sense that one needs to solve, at each time step, decoupled Stokes and Darcy problems. Therefore, legacy Stokes and Darcy solvers can be applied in parallel. the scheme is also unconditionally stable and, with a mild time-step restriction, long-time accurate in the sense that the error is bounded uniformly in time. Numerical experiments are used to illustrate the theoretical results. to the authors\u27 knowledge, the novel algorithm is the first third-order accurate numerical scheme for the Stokes–Darcy system possessing its favorable efficiency, stability, and accuracy properties

Weak Galerkin finite element methods for elasticity and coupled flow problems

Includes bibliographical references.2020 Summer.We present novel stabilizer-free weak Galerkin finite element methods for linear elasticity and coupled Stokes-Darcy flow with a comprehensive treatment of theoretical results and the numerical methods for each. Weak Galerkin finite element methods take a discontinuous approximation space and bind degrees of freedom together through the discrete weak gradient, which involves solving a small symmetric positive-definite linear system on every element of the mesh. We introduce notation and analysis using a general framework that highlights properties that unify many existing weak Galerkin methods. This framework makes analysis for the methods much more straightforward. The method for linear elasticity on quadrilateral and hexahedral meshes uses piecewise constant vectors to approximate the displacement on each cell, and it uses the Raviart-Thomas space for the discrete weak gradient. We use the Schur complement to simplify the solution of the global linear system and increase computational efficiency further. We prove first-order convergence in the L2 norm, verify our analysis with numerical experiments, and compare to another weak Galerkin approach for this problem. The method for coupled Stokes-Darcy flow uses an extensible multinumerics approach on quadrilateral meshes. The Darcy flow discretization uses a weak Galerkin finite element method with piecewise constants approximating pressure and the Arbogast-Correa space for the weak gradient. The Stokes domain discretization uses the classical Bernardi-Raugel pair. We prove first-order convergence in the energy norm and verify our analysis with numerical experiments. All algorithms implemented in this dissertation are publicly available as part of James Liu's DarcyLite and Darcy+ packages and as part of the deal.II library

Turbulence: Numerical Analysis, Modelling and Simulation

The problem of accurate and reliable simulation of turbulent flows is a central and intractable challenge that crosses disciplinary boundaries. As the needs for accuracy increase and the applications expand beyond flows where extensive data is available for calibration, the importance of a sound mathematical foundation that addresses the needs of practical computing increases. This Special Issue is directed at this crossroads of rigorous numerical analysis, the physics of turbulence and the practical needs of turbulent flow simulations. It seeks papers providing a broad understanding of the status of the problem considered and open problems that comprise further steps