Related papers: GiRaFFE: An Open-Source General Relativistic Force…
Modern astrophysical simulations aim to accurately model an ever-growing array of physical processes, including the interaction of fluids with magnetic fields, under increasingly stringent performance and scalability requirements driven by…
This paper describes the design and implementation of our new multi-group, multi-dimensional radiation hydrodynamics (RHD) code Fornax and provides a suite of code tests to validate its application in a wide range of physical regimes.…
Relativistic electron fluxes observed in geosynchronous orbit by GOES-8 in 1997 to 2000 were considered as a complex function of geomagnetic indices PC, Kp, and Dst as well as parameters of the magnetosphere size, subsolar Rs and terminator…
It is shown that force-free electrodynamics (FFE) breaks down in regions where $B^2 -E^2 <0$ (electric zones) even if $\pmb{E}\cdot\pmb{B} =0$. Spontaneous creation of such regions will inevitably lead to plasma oscillations that will…
Rotation measure (RM) synthesis is a widely used polarization processing algorithm for reconstructing polarized structures along the line of sight. Performing RM synthesis on large datasets produced by telescopes like LOFAR can be…
General relativistic Riemann solvers are typically complex, fragile and unwieldy, at least in comparison to their special relativistic counterparts. In this paper, we present a new high-resolution shock-capturing algorithm on curved…
The advent of high-fidelity imaging of supermassive black holes calls for efficient and robust data-analysis methods. In this work, we use $\texttt{Jipole}$, a differentiable, $\texttt{ipole}$-based radiative transfer code, to enable…
Electrical resistivity plays a fundamental role in many astrophysical systems, influencing the evolution of the magnetic field and energy dissipation processes. During the coalescence of two neutron stars, resistive effects can…
Force-Free Electrodynamics (FFE) is known to describe the highly magnetized plasma around pulsars and astrophysical black holes. The equations describing FFE are highly non-linear and the task of finding a well-defined analytical solution…
We implement general relativistic hydrodynamics in the moving-mesh code AREPO. We also couple a solver for the Einstein field equations employing the conformal flatness approximation. The implementation is validated by evolving isolated…
We extend the general-relativistic magnetohydrodynamics (GRMHD) capabilities of Athena++ to incorporate radiation. The intensity field in each finite-volume cell is discretized in angle, with explicit transport in both space and angle…
We perform two-dimensional (2D) implosion simulations using a Monte Carlo kinetic particle code. The paper is motivated by the importance of non-equilibrium effects in inertial confinement fusion (ICF) capsule implosions. These cannot be…
The auxiliary mode approach to time-dependent open quantum system calculations is implemented and refined to yield a feasible computational approach to simulate nanostructures far from equilibrium. It is done by a careful diagonalization of…
Soft robots, particularly magnetic soft robots, require specialized simulation tools to accurately model their deformation under external magnetic fields. However, existing platforms often lack dedicated support for magnetic materials,…
One of the main challenges in the numerical modeling of binary neutron-star (BNS) mergers is long-term simulations of the post-merger remnant over timescales of the order of seconds. When this modeling includes all the aspects of complex…
A Hamiltonian two-field gyrofluid model is used to investigate the dynamics of an electron-ion collisionless plasma subject to a strong ambient magnetic field, within a spectral range extending from the magnetohydrodynamic (MHD) scales to…
We present a new special relativistic hydrodynamics (SRHD) code capable of handling coexisting ultra-relativistically hot and non-relativistically cold gases. We achieve this by designing a new algorithm for conversion between primitive and…
We discuss the cosmological simulation code GADGET-2, a new massively parallel TreeSPH code, capable of following a collisionless fluid with the N-body method, and an ideal gas by means of smoothed particle hydrodynamics (SPH). Our…
Motivated by the desire for highly accurate numerical computations of compact binary spacetimes in the era of gravitational wave astronomy, we reexamine hyperbolicity and well-posedness of the initial value problem for popular models of…
Orbital-Free Density Functional Theory (OFDFT) has re-emerged as a viable alternative to Kohn-Sham DFT, driven by recent advances in kinetic energy density functionals (KEDFs). Nonlocal (NL) KEDFs have significantly extended OFDFT's…