Computational Science Laboratory
Permanent URI for this community
The mission of the Computational Science Laboratory (CSL) is to
develop innovative computational solutions for complex real-world problems, and to
foster a productive research and education environment emphasizing collaboration and innovation.
Browse
Browsing Computational Science Laboratory by Subject "0103 Numerical and Computational Mathematics"
Now showing 1 - 6 of 6
Results Per Page
Sort Options
- Alternating directions implicit integration in a general linear method frameworkSarshar, Arash; Roberts, Steven; Sandu, Adrian (Elsevier, 2021-05-15)Alternating Directions Implicit (ADI) integration is an operator splitting approach to solve parabolic and elliptic partial differential equations in multiple dimensions based on solving sequentially a set of related one-dimensional equations. Classical ADI methods have order at most two, due to the splitting errors. Moreover, when the time discretization of stiff one-dimensional problems is based on Runge–Kutta schemes, additional order reduction may occur. This work proposes a new ADI approach based on the partitioned General Linear Methods framework. This approach allows the construction of high order ADI methods. Due to their high stage order, the proposed methods can alleviate the order reduction phenomenon seen with other schemes. Numerical experiments are shown to provide further insight into the accuracy, stability, and applicability of these new methods.
- Linearly implicit GARK schemesSandu, Adrian; Guenther, Michael; Roberts, Steven (Elsevier, 2021-03-01)Systems driven by multiple physical processes are central to many areas of science and engineering. Time discretization of multiphysics systems is challenging, since different processes have different levels of stiffness and characteristic time scales. The multimethod approach discretizes each physical process with an appropriate numerical method; the methods are coupled appropriately such that the overall solution has the desired accuracy and stability properties. The authors developed the general-structure additive Runge–Kutta (GARK) framework, which constructs multimethods based on Runge–Kutta schemes. This paper constructs the new GARK-ROS/GARK-ROW families of multimethods based on linearly implicit Rosenbrock/Rosenbrock-W schemes. For ordinary differential equation models, we develop a general order condition theory for linearly implicit methods with any number of partitions, using exact or approximate Jacobians. We generalize the order condition theory to two-way partitioned index-1 differential-algebraic equations. Applications of the framework include decoupled linearly implicit, linearly implicit/explicit, and linearly implicit/implicit methods. Practical GARK-ROS and GARK-ROW schemes of order up to four are constructed.
- Multirate implicit Euler schemes for a class of differential-algebraic equations of index-1Hachtel, Christoph; Bartel, Andreas; Guenther, Michael; Sandu, Adrian (Elsevier, 2021-05-15)Systems of differential equations which consist of subsystems with widely different dynamical behaviour can be integrated by multirate time integration schemes to increase the efficiency. These schemes allow the usage of inherent step sizes according to the dynamical properties of the subsystem. In this paper, we extend the multirate implicit Euler method to semi-explicit differential–algebraic equations of index-1 where the algebraic constraints only occur in the slow changing subsystem. We discuss different coupling approaches and show that consistency and convergence order 1 can be reached. Numerical experiments validate the analytical results.
- Multirate linearly-implicit GARK schemesGuenther, Michael; Sandu, Adrian (Springer, 2021-12-28)Many complex applications require the solution of initial-value problems where some components change fast, while others vary slowly. Multirate schemes apply different step sizes to resolve different components of the system, according to their dynamics, in order to achieve increased computational efficiency. The stiff components of the system, fast or slow, are best discretized with implicit base methods in order to ensure numerical stability. To this end, linearly implicit methods are particularly attractive as they solve only linear systems of equations at each step. This paper develops the Multirate GARK-ROS/ROW (MR-GARK-ROS/ROW) framework for linearly-implicit multirate time integration. The order conditions theory considers both exact and approximative Jacobians. The effectiveness of implicit multirate methods depends on the coupling between the slow and fast computations; an array of efficient coupling strategies and the resulting numerical schemes are analyzed. Multirate infinitesimal step linearly-implicit methods, that allow arbitrarily small micro-steps and offer extreme computational flexibility, are constructed. The new unifying framework includes existing multirate Rosenbrock(-W) methods as particular cases, and opens the possibility to develop new classes of highly effective linearly implicit multirate integrators.
- Partitioned exponential methods for coupled multiphysics systemsNarayanamurthi, Mahesh; Sandu, Adrian (Elsevier, 2021-03-01)Multiphysics problems involving two or more coupled physical phenomena are ubiquitous in science and engineering. This work develops a new partitioned exponential approach for the time integration of multiphysics problems. After a possible semi-discretization in space, the class of problems under consideration is modeled by a system of ordinary differential equations where the right-hand side is a summation of two component functions, each corresponding to a given set of physical processes. The partitioned-exponential methods proposed herein evolve each component of the system via an exponential integrator, and information between partitions is exchanged via coupling terms. The traditional approach to constructing exponential methods, based on the variation-of-constants formula, is not directly applicable to partitioned systems. Rather, our approach to developing new partitioned-exponential families is based on a general-structure additive formulation of the schemes. Two method formulations are considered, one based on a linear-nonlinear splitting of the right hand component functions, and another based on approximate Jacobians. The paper develops classical (non-stiff) order conditions theory for partitioned exponential schemes based on particular families of T-trees and B-series theory. Several practical methods of third order are constructed that extend the Rosenbrock-type and EPIRK families of exponential integrators. Several implementation optimizations specific to the application of these methods to reaction-diffusion systems are also discussed. Numerical experiments reveal that the new partitioned-exponential methods can perform better than traditional unpartitioned exponential methods on some problems.
- Subspace adaptivity in Rosenbrock-Krylov methods for the time integration of initial value problemsTranquilli, Paul; Glandon, Ross; Sandu, Adrian (Elsevier, 2021-03-15)The Rosenbrock–Krylov family of time integration schemes is an extension of Rosenbrock-W methods that employs a specific Krylov based approximation of the linear system solutions arising within each stage of the integrator. This work proposes an extension of Rosenbrock–Krylov methods to address stability questions which arise for methods making use of inexact linear system solution strategies. Two approaches for improving the stability and efficiency of Rosenbrock–Krylov methods are proposed, one through direct control of linear system residuals and the second through a novel extension of the underlying Krylov space to include stage right hand side vectors. Rosenbrock–Krylov methods employing the new approaches show a substantial improvement in computational efficiency relative to prior implementations.