Related papers: An Intercomparison Between Divergence-Cleaning and…
Invariant discretization schemes are derived for the one- and two-dimensional shallow-water equations with periodic boundary conditions. While originally designed for constructing invariant finite difference schemes, we extend the usage of…
Meshing is a critical, but user-intensive process necessary for stable and accurate simulations in computational fluid dynamics (CFD). Mesh generation is often a bottleneck in CFD pipelines. Adaptive meshing techniques allow the mesh to be…
The ideal Magnetohydrodynamics (MHD) equations are challenging because one needs to maintain the divergence free condition, $\nabla \cdot \Bv = 0$. Many numerical methods have been developed to enforce this condition. In this work, we…
We consider the application of the magnetic flux leakage (MFL) method to the detection of defects in ferromagnetic (steel) tubulars. The problem setup corresponds to the cases where the distance from the casing and the point where the…
Magnetic reconnection occurring in collisionless environments is a multi-scale process involving both ion and electron kinetic processes. Because of their small mass, the electron scales are difficult to resolve in numerical and satellite…
The advent of robust, reliable and accurate higher order Godunov schemes for many of the systems of equations of interest in computational astrophysics has made it important to understand how to solve them in multi-scale fashion. This is so…
Computations have helped elucidate the dynamics of Earth's mantle for several decades already. The numerical methods that underlie these simulations have greatly evolved within this time span, and today include dynamically changing and…
In this paper we present a new high order semi-implicit DG scheme on two-dimensional staggered triangular meshes applied to different nonlinear systems of hyperbolic conservation laws such as advection-diffusion models, incompressible…
A high-fidelity finite volume scheme based on the BVD (boundary variation diminishing) concept is proposed in this study to solve the ideal magnetohydrodynamics (MHD) equations. A hybrid spatial reconstruction profile, consisting of a…
We continue the investigation of how to use the divergence-free condition to resolve the azimuthal ambiguity present in vector magnetogram data. In previous articles, by Crouch, Barnes, and Leka (Solar Physics, 260, 271, 2009) and Crouch…
We investigate how the azimuthal ambiguity in solar vector magnetogram data can be resolved by using the divergence-free property of magnetic fields. In a previous article, by Crouch, Barnes, and Leka (Solar Phys. 260, 271, 2009),…
Discretized numerical simulations are a powerful tool for investigation of nonlinear MHD turbulence in accretion disks. However, confidence in their quantitative predictions requires a demonstration that further refinement of the spatial…
A description is given of the algorithms implemented in the AstroBEAR adaptive mesh refinement code for ideal magnetohydrodynamics. The code provides several high resolution, shock capturing schemes which are constructed to maintain…
Scale invariance is a hallmark of many natural systems, including solar flares, where energy release spans a vast range of scales. Recent computational advances, at the level of both algorithmics and hardware, have enabled high-resolution…
We present a new algorithm for numerical magnetohydrodynamics on staggered meshes preserving $\nabla \cdot B = 0$. Our algorithm is based on the constrained transport method and supports both cell-based adaptive mesh refinement and temporal…
Numerical simulations of the propagation of charged particles through magnetic fields solving the equation of motion often leads to the usage of an interpolation in case of discretely defined magnetic fields, typically given on a…
Preserving scalar boundedness is important for numerical schemes used in turbulent compressible multi-component flow simulations to prevent unphysical results and unstable simulations. However, ensuring scalar boundedness for high-order,…
The generation of large-scale magnetic field in the kinematic regime in the absence of an alpha-effect is investigated by following two different approaches, namely the test-field method and multiscale stability theory relying on the…
Numerical methods for solving the ideal magnetohydrodynamic (MHD) equations in more than one space dimension must confront the challenge of controlling errors in the discrete divergence of the magnetic field. One approach that has been…
The tearing mode instability is a key process for magnetic energy conversion in magnetohydrodynamics, once anti-parallel components are allowed to reconnect, leading to the formation of magnetic islands. It has been employed to explain…