Related papers: Calliope: Pseudospectral shearing magnetohydrodyna…
We find and investigate via numerical simulations self-sustained two-dimensional turbulence in a magnetohydrodynamic flow with a maximally simple configuration: plane, noninflectional (with a constant shear of velocity) and threaded by a…
We study a turbulent helical dynamo in a periodic domain by solving the ideal magnetohydrodynamic (MHD) equations with the FLASH code using the divergence-cleaning eight-wave method and compare our results with direct numerical simulations…
The magnetorotational instability (MRI) is an important process in sufficiently ionized accretion disks, as it can create turbulence that acts as an effective viscosity, mediating angular momentum transport. Due to its local nature, it is…
Multigrid solvers are among the most efficient methods for solving the Poisson equation, which is ubiquitous in computational physics. For example, in the context of incompressible flows, it is typically the costliest operation. The present…
Especially in cold and high-density regions, the assumptions of ideal magnetohydrodynamics (MHD) can break down, making first order non-ideal terms such as Ohmic and ambipolar diffusion as well as the Hall effect important. In this study we…
In this work, a simple fourth-order accurate finite volume semi-discrete scheme is introduced to solve astrophysical magnetohydrodynamics (MHD) problems on Cartesian meshes. Hydrodynamic quantities like density, momentum and energy are…
We employ simulations of supersonic super-Alfvenic turbulence decay as a benchmark test problem to assess and compare the performance of nine astrophysical MHD methods actively used to model star formation. The set of nine codes includes:…
Turbulence in compressible plasma plays a key role in many areas of astrophysics and engineering. The extreme plasma parameters in these environments, e.g. high Reynolds numbers, supersonic and super-Alfvenic flows, however, make direct…
We propose a seamless multiscale method which approximates the macroscopic behavior of the passive advection-diffusion equations with steady incompressible velocity fields with multi-spatial scales. The method uses decompositions of the…
We introduce an innovative wavelet-based approach to dynamically adjust the local grid resolution to maintain a uniform specified error tolerance. Extending the work of Dubos and Kevlahan (2013), a wavelet multi-scale approximation is used…
Direct numerical simulations of decaying and forced magnetohydrodynamic (MHD) turbulence without and with mean magnetic field are analyzed by higher-order two-point statistics. The turbulence exhibits statistical anisotropy with respect to…
We present a hybrid machine learning framework that combines Physics-Informed Neural Operators (PINOs) with score-based generative diffusion models to simulate the full spatio-temporal evolution of two-dimensional, incompressible, resistive…
Many scientific and engineering problems involving multi-physics span a wide range of scales. Understanding the interactions across these scales is essential for fully comprehending such complex problems. However, visualizing multivariate,…
We present a new analysis of the anisotropic spectral energy distribution in incompressible magnetohydrodynamic (MHD) turbulence permeated by a strong mean magnetic field. The turbulent flow is generated by high-resolution pseudo-spectral…
We present a new numerical method of special relativistic resistive magnetohydrodynamics with scalar resistivity that can treat a range of phenomena, from nonrelativistic to relativistic (shock, contact discontinuity, and Alfv\'en wave).…
We extend recently-developed mesh-free Lagrangian methods for numerical magnetohydrodynamics (MHD) to arbitrary anisotropic diffusion equations, including: passive scalar diffusion, Spitzer-Braginskii conduction and viscosity, cosmic ray…
Fully realizing the potential of multigrid solvers often requires custom algorithms for a given application model, discretizations and even regimes of interest, despite considerable effort from the applied math community to develop fully…
Astrophysical turbulent flows display an intrinsically multi-scale nature, making their numerical simulation and the subsequent analyses of simulated data a complex problem. In particular, two fundamental steps in the study of turbulent…
We perform direct 3-dimensional numerical simulations for magnetohydrodynamic (MHD) turbulence in a periodic box of size $2\pi$ threaded by strong uniform magnetic fields. We use a pseudo-spectral code with hyperviscosity and…
Fluid discontinuities, such as shock fronts and vortex sheets, can reflect waves and become unstable to corrugation. Analytical calculations of these phenomena are tractable in the simplest cases only, while their numerical simulations are…