Related papers: Average energy dissipation rates of additive impli…
This paper proposes an implicit family of sub-step integration algorithms grounded in the explicit singly diagonally implicit Runge-Kutta (ESDIRK) method. The proposed methods achieve third-order consistency per sub-step and thus the…
Multirate integration is an increasingly relevant tool that enables scientists to simulate multiphysics systems. Existing multirate methods are designed for equations whose fast and slow variables can be linearly separated using additive or…
Implicit-explicit Runge-Kutta (IMEX-RK) schemes are popular methods to treat multiscale equations that contain a stiff part and a non-stiff part, where the stiff part is characterized by a small parameter $\varepsilon$. In this work, we…
Isospectral Runge-Kutta methods are well-suited for the numerical solution of isospectral systems such as the rigid body and the Toda lattice. More recently, these integrators have been applied to geophysical fluid models, where their…
We introduce a class of high order accurate, semi-implicit Runge-Kutta schemes in the general setting of evolution equations that arise as gradient flow for a cost function, possibly with respect to an inner product that depends on the…
In this paper, we apply the Paired-Explicit Runge-Kutta (P-ERK) schemes by Vermeire et. al. (2019, 2022) to dynamically partitioned systems arising from adaptive mesh refinement. The P-ERK schemes enable multirate time-integration with no…
In this paper we discuss the use of implicit Runge-Kutta schemes for the time discretization of optimal control problems with evolution equations. The specialty of the considered discretizations is that the discretizations schemes for the…
New time integration methods are proposed for simulating incompressible multiphase flow in pipelines described by the one-dimensional two-fluid model. The methodology is based on 'half-explicit' Runge-Kutta methods, being explicit for the…
We present a dynamically load-balanced parallel $ p $-adaptive implicit high-order flux reconstruction method for under-resolved turbulence simulation. The high-order explicit first stage, singly diagonal implicit Runge-Kutta (ESDIRK)…
In this paper, we conduct an in-depth investigation of the structural intricacies inherent to the Invariant Energy Quadratization (IEQ) method as applied to gradient flows, and we dissect the mechanisms that enable this method to uphold…
We present MDIRK: a Multifluid second-order Diagonally-Implicit Runge-Kutta method to study momentum transfer between gas and an arbitrary number ($N$) of dust species. The method integrates the equations of hydrodynamics with an Implicit…
This work gives the asymptotic error distribution of the stochastic Runge--Kutta (SRK) method of strong order $1$ applied to Stratonovich-type stochastic differential equations. For dealing with the implicitness introduced in the diffusion…
In \cite{ZH2019}, we developed a boundary treatment method for implicit-explicit (IMEX) Runge-Kutta (RK) methods for solving hyperbolic systems with source terms. Since IMEX RK methods include explicit ones as special cases, this boundary…
A mixed accuracy framework for Runge--Kutta methods presented in [Grant, JSC 2022] has been shown to speed up the computation in diagonally implicit Runge--Kutta (DIRK) methods by using less expensive low accuracy approaches for the…
Many control, optimization, and learning algorithms rely on discretizations of continuous-time contracting systems, where preservation of contractivity under numerical integration is key for stability, robustness, and reliable fixed-point…
Compact Runge-Kutta (cRK) methods are a class of high order methods for solving hyperbolic conservation laws characterized by their compact stencil including only immediate neighboring finite elements. A Compact Runge-Kutta flux…
We present second-order optimally stable Implicit-Explicit (IMEX) Runge-Kutta (RK) schemes with application to a modified set of shallow water equations that can be used to model the dynamics of lava flows. The schemes are optimally stable…
The effects of kinetic-energy preservation errors due to Runge-Kutta (RK) temporal integrators have been analyzed for the case of large-eddy simulations of incompressible turbulent channel flow. Simulations have been run using the…
We present a paradigm for developing arbitrarily high order, linear, unconditionally energy stable numerical algorithms for gradient flow models. We apply the energy quadratization (EQ) technique to reformulate the general gradient flow…
An explicit numerical strategy that practically preserves invariants is derived for conservative systems by combining an explicit high-order Runge-Kutta (RK) scheme with a simple modification of the standard projection approach, which is…