Related papers: Scalable explicit implementation of anisotropic di…
We study symplectic numerical integration of mechanical systems with a Hamiltonian specified in non-canonical coordinates and its application to guiding-center motion of charged plasma particles in magnetic confinement devices. The…
A novel 3-D higher-order finite-difference time-domain framework with complex frequency-shifted perfectly matched layer for the modeling of wave propagation in cold plasma is presented. Second- and fourth-order spatial approximations are…
Rapidly rotating Rayleigh-B\'enard convection admits a class of exact steady single-mode solutions describing high-amplitude convection cells. Using a matched asymptotic analysis in the high-Rayleigh-number limit, we obtain a rigorous…
Simulating complex astrophysical reacting flows is computationally expensive -- reactions are stiff and typically require implicit integration methods. The reaction update is often the most expensive part of a simulation, which motivates…
In this master thesis we have compared different second order stabilized explicit Runge-Kutta methods when applied to the incompressible Navier-Stokes equations by means of a projection method and a differential algebraic approach. We…
In this work, we introduce two novel reformulations of a recent weakly hyperbolic model for two-phase flow with surface tension. In the model, the tracking of phase boundaries is achieved by using a vector interface field, rather than a…
Thermal radiative transfer (TRT) presents significant computational challenges due to the stiff, nonlinear coupling between radiation and material energy, particularly in multigroup, high-fidelity transport models. In this work, we develop…
This work introduces a new class of Runge-Kutta methods for solving nonlinearly partitioned initial value problems. These new methods, named nonlinearly partitioned Runge-Kutta (NPRK), generalize existing additive and component-partitioned…
We propose an experimental study of adaptive time-stepping methods for efficient modeling of the aggregation-fragmentation kinetics. Precise modeling of this phenomena usually requires utilization of the large systems of nonlinear ordinary…
When dealing with shallow water simulations, the velocity profile is often assumed to be constant along the vertical axis. However, since in many applications this is not the case, modeling errors can be significant. Hence, in this work, we…
Singly-TASE operators for the numerical solution of stiff differential equations were proposed by Calvo et al. in J.Sci. Comput. 2023 to reduce the computational cost of Runge-Kutta-TASE (RKTASE) methods when the involved linear systems are…
Based on the recent development of Jacobian-free Lax-Wendroff (LW) approaches for solving hyperbolic conservation laws [Zorio, Baeza and Mulet, Journal of Scientific Computing 71:246-273, 2017], [Carrillo and Par\'es, Journal of Scientific…
Implicit schemes have been extensively used in building physics to compute the solution of moisture diffusion problems in porous materials for improving stability conditions. Nevertheless, these schemes require important sub-iterations when…
Explicit Runge-Kutta schemes become impractical when a stiff linear operator is present in the dynamics. This failure mode is quite common in numerical simulations of fluids and plasmas. Lawson proposed Generalized Runge-Kutta Processes for…
The recently-introduced relaxation approach for Runge-Kutta methods can be used to enforce conservation of energy in the integration of Hamiltonian systems. We study the behavior of implicit and explicit relaxation Runge-Kutta methods in…
It is difficult to design high order numerical schemes which could preserve both the maximum bound property (MBP) and energy dissipation law for certain phase field equations. Strong stability preserving (SSP) Runge-Kutta methods have been…
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 develop entropy dissipative higher order accurate local discontinuous Galerkin (LDG) discretizations coupled with Diagonally Implicit Runge-Kutta (DIRK) methods for nonlinear degenerate parabolic equations with a gradient flow structure.…
We develop numerical methods for reaction-diffusion systems based on the equations of fluctuating hydrodynamics (FHD). While the FHD formulation is formally described by stochastic partial differential equations (SPDEs), it becomes similar…
In this paper, we develop a general framework for constructing higher-order, unconditionally energy-stable exponential time differencing Runge-Kutta methods applicable to a range of gradient flows. Specifically, we identify conditions…