Related papers: Implementation of the Shearing Box Approximation i…
For more than a decade, the so-called shearing box model has been used to study the fundamental local dynamics of accretion discs. This approach has proved to be very useful because it allows high resolution and long term studies to be…
The set of equations for magnetohydrodynamic (MHD) waves in a shear flow is consecutively derived. The proposed scenario involves the presence of a self-sustained turbulence and magnetic field. In the framework of Langevin--Burgers approach…
We present a new algorithm to solve the equations of radiation hydrodynamics (RHD) in a frequency-integrated, two-moment formulation. Novel features of the algorithm include i) the adoption of a non-local Variable Eddington Tensor (VET)…
We present a formulation and numerical algorithm to extend the scheme for grey radiation magneto-hydrodynamics (MHD) developed by Jiang (2021) to include the frequency dependence via the multi-group approach. The entire frequency space can…
Momentum Stochastic Gradient Descent (MSGD) algorithm has been widely applied to many nonconvex optimization problems in machine learning, e.g., training deep neural networks, variational Bayesian inference, and etc. Despite its empirical…
We propose a general-purpose, self-adaptive approach to construct variational wavefunction ans\"atze for highly accurate quantum dynamics simulations based on McLachlan's variational principle. The key idea is to dynamically expand the…
Rotationally supported, cold, gaseous disks are ubiquitous in astrophysics and appear in a diverse set of systems, such as protoplanetary disks, accretion disks around black holes, or large spiral galaxies. Capturing the gas dynamics…
Energy-preserving numerical methods for solving the Hodge wave equation is developed in this paper. Based on the de Rham complex, the Hodge wave equation can be formulated as a first-order system and mixed finite element methods using…
The trapezoidal rule, which is a special case of the Newmark family of algorithms, is one of the most widely used methods for transient hyperbolic problems. In this work, we show that this rule conserves linear and angular momenta and…
In certain astrophysical systems the commonly employed ideal magnetohydrodynamics (MHD) approximation breaks down. Here, we introduce novel explicit and implicit numerical schemes of ohmic resistivity terms in the moving-mesh code AREPO. We…
We present a new numerical algorithm to solve the recently derived equations of two-moment cosmic ray hydrodynamics (CRHD). The algorithm is implemented as a module in the moving mesh Arepo code. Therein, the anisotropic transport of cosmic…
We present the first three-dimensional magnetohydrodynamic (MHD) simulations of a circumbinary disk surrounding an equal mass binary. The binary maintains a fixed circular orbit of separation $a$. As in previous hydrodynamical simulations,…
We develop structure-preserving time integration schemes for Gaussian wave packet dynamics associated with the magnetic Schr\"odinger equation. The variational Dirac--Frenkel formulation yields a finite-dimensional Hamiltonian system for…
Non-smooth optimization models play a fundamental role in various disciplines, including engineering, science, management, and finance. However, classical algorithms for solving such models often struggle with convergence speed,…
We introduce a second-order numerical scheme for compressible atmospheric motions at small to planetary scales. The collocated finite volume method treats the advection of mass, momentum, and mass-weighted potential temperature in…
Vortices, turbulence, and unsteady non-laminar flows are likely both prominent and dynamically important features of astrophysical disks. Such strongly nonlinear phenomena are often difficult, however, to simulate accurately, and are…
We propose Quantum Riemannian Hamiltonian Descent (QRHD), a quantum algorithm for continuous optimization on Riemannian manifolds that extends Quantum Hamiltonian Descent (QHD) by incorporating geometric structure of the parameter space via…
We present a new magnetohydrodynamic (MHD) code for the simulation of wave propagation in the solar atmosphere, under the effects of electrical resistivity, but not dominant, and heat transference in a uniform 3D grid. The code is based on…
We present a new computational method for smoothly matching general relativistic ideal magnetohydrodynamics (MHD) to its force-free limit. The method is based on a flux-conservative formalism for MHD and its force-free limit, and a vector…
We propose a variational alternative to the Trotter-Suzuki decomposition that provides greater control over errors while preserving the unitary structure of time evolution. The variational parameters in our ansatz are derived from a global…