Related papers: Adaptive mesh refinement with spectral accuracy fo…
We recently developed a new general relativistic magnetohydrodynamic code with adaptive mesh refinement that evolves the electromagnetic (EM) vector potential (A) instead of the magnetic fields directly. Evolving A enables one to use any…
We present numerical experiments for geophysics electromagnetic (EM) modeling based upon high-order edge elements and supervised $h+p$ refinement approaches on massively parallel computers. Our high-order $h+p$ refinement strategy is based…
The present paper introduces a class of finite volume schemes of increasing order of accuracy in space and time for hyperbolic systems that are in conservation form. This paper specifically focuses on Euler system that is used for modeling…
We present Variable Eddington Tensor-closed Transport on Adaptive Meshes (\texttt{VETTAM}), a new algorithm to solve the equations of radiation hydrodynamics (RHD) with support for adaptive mesh refinement (AMR) in a frequency-integrated,…
A new implementation for the time evolution of the magnetic vector potential is obtained for smoothed particle magnetohydrodynamics by considering the induction equation in integral form. Galilean invariance is achieved through proper gauge…
There is great interest in numerical relativity simulations involving matter due to the likelihood that binary compact objects involving neutron stars will be detected by gravitational wave observatories in the coming years, as well as to…
The plasma response to Resonant Magnetic Perturbations (RMPs) in ASDEX Upgrade is modeled with the non-linear resistive MHD code JOREK, using input profiles that match those of the experiments as closely as possible. The RMP configuration…
In this paper two theoretical approaches for the calculation of the rate of quasi-stationary, two-dimensional magnetic reconnection with nonuniform anomalous resistivity are considered in the framework of incompressible magnetohydrodynamics…
Adaptive mesh refinement is central to the efficient solution of partial differential equations (PDEs) via the finite element method (FEM). Classical $r$-adaptivity optimizes vertex positions but requires solving expensive auxiliary PDEs…
We present a fourth-order projection method with adaptive mesh refinement (AMR) for numerically solving the incompressible Navier-Stokes equations (INSE) with subcycling in time. Our method features (i) a reformulation of INSE so that the…
We introduce a hybrid "Modified Genetic Algorithm-Multilevel Stochastic Gradient Descent" (MGA-MSGD) training algorithm that considerably improves accuracy and efficiency of solving 3D mechanical problems described, in strong-form, by PDEs…
We solve the general relativistic magnetohydrodynamics equations using distributed parallel adaptive mesh refinement. We discuss strong scaling tests of the code, and present evolutions of Michel accretion and a TOV star.
Algebraic Multigrid (AMG) methods are state-of-the-art algebraic solvers for partial differential equations. Still, their efficiency depends heavily on the choice of suitable parameters and/or ingredients. Paradigmatic examples include the…
Worldwide very long baseline radio interferometry arrays are expected to obtain horizon-scale images of supermassive black hole candidates as well as of relativistic jets in several nearby active galactic nuclei. This motivates the…
Adaptive mesh refinement (AMR) is often used when solving time-dependent partial differential equations using numerical methods. It enables time-varying regions of much higher resolution, which can be used to track discontinuities in the…
Solving the Euler equations of ideal hydrodynamics as accurately and efficiently as possible is a key requirement in many astrophysical simulations. It is therefore important to continuously advance the numerical methods implemented in…
We develop an improved lattice-Boltzmann numerical scheme to solve magnetohydrodynamic (MHD) equations in the regime of low magnetic Reynolds numbers, grounded on a manifestly Galilean covariant modeling of the Navier-Stokes equations. The…
We modify an existing magnetohydrodynamics algorithm to make it more compatible with a dimensionally-split (DS) framework. It is based on the standard reconstruct-solve-average strategy (using a Riemann solver), and relies on constrained…
We study dendritic microstructure evolution using an adaptive grid, finite element method applied to a phase-field model. The computational complexity of our algorithm, per unit time, scales linearly with system size, rather than the…
Higher order finite volume schemes for magnetohydrodynamics (MHD) and relativistic magnetohydrodynamics (RMHD) are very valuable because they allow us to carry out astrophysical simulations with very high accuracy. However, astrophysical…