Related papers: A Fully Explicit Integrator for Modeling Astrophys…
We analyze the behaviour of an ensemble of time integrators applied to the semi-discrete problem resulting from the spectral discretization of the equations describing Boussinesq convection in a cylindrical annulus. The equations are cast…
We study the conservation properties of the Hermite-discontinuous Galerkin (Hermite-DG) approximation of the Vlasov-Maxwell equations. In this semi-discrete formulation, the total mass is preserved independently for every plasma species.…
A wide range of physical phenomena exhibit auxiliary admissibility criteria, such as conservation of entropy or various energies, which arise implicitly under the exact solution of their governing PDEs. However, standard temporal schemes,…
Computational chemical combustion problems are known to be stiff, and are typically solved with implicit time integration methods. A novel exponential time integrator, EPI3V, is introduced and applied to a spatially homogeneous isobaric…
We propose a novel approach for modeling chemical reactions within the particle-based Fokker-Planck framework for gas flow simulations which conserves mass, momentum, and energy while retaining the performance advantages of the…
(Abridged) In this paper we present a compilation of results from our most advanced neutron star merger simulations, including a description of the employed numerical procedures and a more complete overview over a large number of computed…
The structural flexibility of the exponential propagation iterative methods of Runge-Kutta type (EPIRK) enables construction of particularly efficient exponential time integrators. While the EPIRK methods have been shown to perform well on…
The paper explores the numerical stability and the computational efficiency of a direct method for unfolding the resolution function from the measurements of the neutron induced reactions. A detailed resolution function formalism is laid…
The present work is devoted to developing the integrated hydrokinetic model (iHKM) for relativistic nucleus-nucleus collisions. While the previous cycle of works on this topic focused on ultra-relativistic collisions at the top RHIC and…
Gradient flows of unit vector fields arise in a wide range of physical models such as harmonic map heat flows, nematic liquid crystals, and magnetization dynamics. Designing numerical schemes that simultaneously preserve the unit length…
Accurately modelling cold and ultracold reactive collisions occuring over deep potential wells, such as \ce{D+ + H2 -> H+ + HD}, requires the development of new theoretical and computational methodologies. One potentially useful framework…
We present an improved method for computing incompressible viscous flow around suspended rigid particles using a fixed and uniform computational grid. The main idea is to incorporate Peskin's regularized delta function approach [Acta…
This paper develops a new class of linearly implicit time integration schemes called Linearly-Implicit Runge-Kutta-W (LIRK-W) methods. These schemes are based on an implicit-explicit approach which does not require a splitting of the right…
This study explores in depth rotating detonation engines (RDEs) fueled by premixed stoichiometric hydrogen/air mixtures through two-dimensional numerical simulations including a detailed chemical kinetic mechanism. To model the spatial…
Physics based emulators offer a fast and reliable replacement for an exact solution of the scattering problem in nuclear physics. Previous work developed a reduced-basis emulator for single-channel elastic scattering using an optical…
In two preceding papers we have shown that, when reaction networks are well-removed from equilibrium, explicit asymptotic and quasi-steady-state approximations can give algebraically-stabilized integration schemes that rival standard…
In many natural and industrial applications, turbulent flows encompass some form of dispersed particles. Although this type of multiphase turbulent flow is omnipresent, its numerical modeling has proven to be a remarkably challenging…
We present an implementation of a damped response framework for calculating resonant inelastic X-ray scattering (RIXS) at the equation-of-motion coupled cluster singles and doubles (CCSD) and second-order approximate coupled cluster singles…
Dense plasmas occur in stars, giant planets and in inertial fusion experiments. Accurate modeling of the electronic structure of these plasmas allows for prediction of material properties that can in turn be used to simulate these…
This paper presents a rank-adaptive implicit-explicit integrator for the tensor approximation of three-dimensional convection-diffusion equations. In particular, the recently developed Reduced Augmentation Implicit Low-rank (RAIL)…