Many practical problems in science and engineering are modeled by large systems of ordinary differential equations (ODEs) which arise from discretization in space of partial differential equations (PDEs) by finite difference methods, finite elements or finite volume methods, or pseudospectral methods. For such systems there are often natural splittings of the right hand sides of the differential systems into two parts, one of which is non-stiff or mildly stiff, and suitable for explicit time integration, and the other part is stiff, and suitable for implicit time integration. The efficient solution can be provided by implicit-explicit (IMEX) schemes. In present research we consider the class of general linear methods (GLMs) for ordinary differential equations. We construct IMEX GLMs of order p = 1, 2, …,4 with desired stability properties. We assume that the explicit and implicit parts of the IMEX scheme have the same abscissa vector c and coefficient matrices U and B. We look for implicit methods of order p and stage order q = p, which have the property of so-called inherent Runge–Kutta stability (IRKS), and which are A-stable and if possible, also L-stable. We also require that the vector of external approximation is of Nordsieck form. Next, we attempt to maximize the combined region of absolute stability. Finally, we apply constructed methods to a series of test problems. This is a joint work with A. Cardone, Z. Jackiewicz and P. Pierzchała.
Search for implicit-explicit general linear methods with inherent Runge-Kutta stability
Angelamaria Cardone;Zdzislaw Jackiewicz;
2019-01-01
Abstract
Many practical problems in science and engineering are modeled by large systems of ordinary differential equations (ODEs) which arise from discretization in space of partial differential equations (PDEs) by finite difference methods, finite elements or finite volume methods, or pseudospectral methods. For such systems there are often natural splittings of the right hand sides of the differential systems into two parts, one of which is non-stiff or mildly stiff, and suitable for explicit time integration, and the other part is stiff, and suitable for implicit time integration. The efficient solution can be provided by implicit-explicit (IMEX) schemes. In present research we consider the class of general linear methods (GLMs) for ordinary differential equations. We construct IMEX GLMs of order p = 1, 2, …,4 with desired stability properties. We assume that the explicit and implicit parts of the IMEX scheme have the same abscissa vector c and coefficient matrices U and B. We look for implicit methods of order p and stage order q = p, which have the property of so-called inherent Runge–Kutta stability (IRKS), and which are A-stable and if possible, also L-stable. We also require that the vector of external approximation is of Nordsieck form. Next, we attempt to maximize the combined region of absolute stability. Finally, we apply constructed methods to a series of test problems. This is a joint work with A. Cardone, Z. Jackiewicz and P. Pierzchała.I documenti in IRIS sono protetti da copyright e tutti i diritti sono riservati, salvo diversa indicazione.