SEARCH · Search NASA

Results for “Divergence-free magnetic fields”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 19 records

All sources NASA NTRS DOE OSTI

DOE OSTI2026

Asymptotic-preserving semi-implicit finite volume scheme for extended magnetohydrodynamics

A Finite Volume (FV) scheme is developed for solving the extended magnetohydrodynamic (XMHD) equations, yielding accurate results in the ideal, resistive, and Hall MHD limits. This is accomplished by first re-writing the XMHD equations such that it allows the algorithm to retain the use of ideal MHD Riemann solvers and the constrained transport method to preserve divergence-free magnetic fields. Incorporation of electron inertia and displacement current introduces additional numerical stiffness which motivates a semi-implicit FV scheme that re-formulates the XMHD model as a relaxation system. The equations are then advanced in time using an explicit 2nd-order Runge–Kutta scheme with operator splitting applied to the implicit source term updates at each sub-stage. For additional numerical stability, a density-dependent slope limiter is implemented to increase flux diffusivity at low density regions where non-ideal effects become significant. The algorithm is subsequently implemented in a scalable adaptive mesh refinement (AMR) framework. As the new algorithm retains many aspects of the ideal MHD formulations, it asymptotes naturally to the ideal MHD limit. Moreover, it shows promising results at the resistive and Hall MHD limits. This is verified against reference test problems for ideal, resistive and Hall MHD.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

NASA NTRS2006

Newly-Developed 3D GRMHD Code and its Application to Jet Formation

We have developed a new three-dimensional general relativistic magnetohydrodynamic code by using a conservative, high-resolution shock-capturing scheme. The numerical fluxes are calculated using the HLL approximate Riemann solver scheme. The flux-interpolated constrained transport scheme is used to maintain a divergence-free magnetic field. We have performed various 1-dimensional test problems in both special and general relativity by using several reconstruction methods and found that the new 3D GRMHD code shows substantial improvements over our previous model. The . preliminary results show the jet formations from a geometrically thin accretion disk near a non-rotating and a rotating black hole. We will discuss the jet properties depended on the rotation of a black hole and the magnetic field strength.

Mizuno, Y.↗

NASA NTRS2006

GRMHD Simulations of Jet Formation with a New Code

Mizuno, Y.↗

NASA NTRS2006

GRMHD Simulations of Jet Formation with a Newly-Developed GRMHD Code

We have developed a new three-dimensional general relativistic magnetohydrodynamic code by using a conservative, high-resolution shock-capturing scheme. The numerical fluxes are calculated using the HLL approximate Riemann solver scheme. The flux-CT scheme is used to maintain a divergence-free magnetic field. Various 1-dimensional test problems show significant improvements over our previous GRMHD code. We have performed simulations of jet formations from a geometrically thin accretion disk near a non-rotating and a rotating black hole. The new simulation results show that the jet is formed by the same manner as in previous works and propagates outward. As the magnetic field strength becomes weaker, larger amount of matter launches with the jet. On the other hand when the magnetic field strength becomes stronger, the jet has less-matter and becomes poynting flux dominated. We will also discuss how the jet properties depend on the rotation of a black hole.

Mizuno, Y.↗

NASA NTRS2006

RAISHIN: A High-Resolution Three-Dimensional General Relativistic Magnetohydrodynamics Code

We have developed a new three-dimensional general relativistic magnetohydrodynamic (GRMHD) code, RAISHIN, using a conservative, high resolution shock-capturing scheme. The numerical fluxes are calculated using the Harten, Lax, & van Leer (HLL) approximate Riemann solver scheme. The flux-interpolated, constrained transport scheme is used to maintain a divergence-free magnetic field. In order to examine the numerical accuracy and the numerical efficiency, the code uses four different reconstruction methods: piecewise linear methods with Minmod and MC slope-limiter function, convex essentially non-oscillatory (CENO) method, and piecewise parabolic method (PPM) using multistep TVD Runge-Kutta time advance methods with second and third-order time accuracy. We describe code performance on an extensive set of test problems in both special and general relativity. Our new GRMHD code has proven to be accurate in second order and has successfully passed with all tests performed, including highly relativistic and magnetized cases in both special and general relativity.

Mizuno, Yosuke↗

NASA NTRS2006

GRMHD Simulations of Jet Formation with RAISHIN

We have developed a new three dimensional general relativistic magnetohydrodynamic (GRMHD) code, RAISHIN, using a conservative, high-resolution shock capturing scheme. Numerical fluxes are calculated using the Harten, Lax, & van Leer (HLL) approximate Riemann solver scheme. The flux-interpolated, constrained transport scheme is used to maintain a divergence-free magnetic field. We describe code performance on some test problems in both special and general relativity. Our new GRMHD code has proven to be accurate to second order and has successfully passed several numerical test problems including highly relativistic and magnetized tests in both special and general relativity. We have performed several simulations of non-rotating and rotating black hole systems with a geometrically thin accretion disk. The simulations show the formation of jets driven by the Lorentz force and the gas pressure. It appears that the rotating black hole creates an additional faster, and more collimated outflow inside a broader, slower outflow that is also generated by the rotating accretion disk around a non-rotating black hole. The kinematic jet structure could thus be a sensitive function of black hole rotation.

Mizuno, Y.↗

NASA NTRS2007

GRMHD and GRPIC Simulations

We have developed a new three-dimensional general relativistic magnetohydrodynamic (GRMHD) code by using a conservative, high-resolution shock-capturing scheme. The numerical fluxes are calculated using the HLL approximate Riemann solver scheme. The flux-interpolated constrained transport scheme is used to maintain a divergence-free magnetic field. We have performed various 1-dimensional test problems in both special and general relativity by using several reconstruction methods and found that the new 3D GRMHD code shows substantial improvements over our previous code. The simulation results show the jet formations from a geometrically thin accretion disk near a nonrotating and a rotating black hole. We will discuss the jet properties depended on the rotation of a black hole and the magnetic field configuration including issues for future research. A General Relativistic Particle-in-Cell Code (GRPIC) has been developed using the Kerr-Schild metric. The code includes kinetic effects, and is in accordance with GRMHD code. Since the gravitational force acting on particles is extreme near black holes, there are some difficulties in numerically describing these processes. The preliminary code consists of an accretion disk and free-falling corona. Results indicate that particles are ejected from the black hole. These results are consistent with other GRMHD simulations. The GRPIC simulation results will be presented, along with some remarks and future improvements. The emission is calculated from relativistic flows in black hole systems using a fully general relativistic radiative transfer formulation, with flow structures obtained by GRMHD simulations considering thermal free-free emission and thermal synchrotron emission. Bright filament-like features protrude (visually) from the accretion disk surface, which are enhancements of synchrotron emission where the magnetic field roughly aligns with the line-of-sight in the co-moving frame. The features move back and forth as the accretion flow evolves, but their visibility and morphology are robust. We would like to extend this research using GRPIC simulations and examine a possible new mechanism for certain X-ray quasi-periodic oscillations (QPOs) observed in blackhole X-ray binaries.

Nishikawa, K.-I.↗

DOE OSTI2025

High order interpolation of magnetic fields with vector potential reconstruction for particle simulations

We propose a method for interpolating divergence-free continuous magnetic fields via vector potential reconstruction using Hermite interpolation, which ensures high-order continuity for applications requiring adaptive, high-order ordinary differential equation (ODE) integrators, such as the Dormand-Prince method. The method provides C(m) continuity and achieves high-order accuracy, making it particularly suited for particle trajectory integration and Poincaré section analysis under optimal integration order and timestep adjustments. Through numerical experiments, we demonstrate that the Hermite interpolation method preserves volume and continuity, which are critical for conserving toroidal canonical momentum and magnetic moment in guiding center simulations, especially over long-term trajectory integration. Furthermore, we analyze the impact of insufficient derivative continuity on Runge-Kutta schemes and show how it degrades accuracy at low error tolerances, introducing discontinuity-induced truncation errors. Lastly, we demonstrate performant Poincaré section analysis in two relevant settings of field data collocated from finite element meshes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

NASA NTRS2008

Resolving the Azimuthal Ambiguity in Vector Magnetogram Data with the Divergence-Free Condition: Theoretical Examination

We demonstrate that the azimuthal ambiguity that is present in solar vector magnetogram data can be resolved with line-of-sight and horizontal heliographic derivative information by using the divergence-free property of magnetic fields without additional assumptions. We discuss the specific derivative information that is sufficient to resolve the ambiguity away from disk center, with particular emphasis on the line-of-sight derivative of the various components of the magnetic field. Conversely, we also show cases where ambiguity resolution fails because sufficient line-of-sight derivative information is not available. For example, knowledge of only the line-of-sight derivative of the line-of-sight component of the field is not sufficient to resolve the ambiguity away from disk center.

Crouch, A.↗

DOE OSTI2025

Performance-portable Binary Neutron Star Mergers with AthenaK

We introduce an extension to the AthenaK code for general-relativistic magnetohydrodynamics (GRMHD) in dynamical spacetimes using a 3+1 conservative Eulerian formulation. Like the fixed-spacetime GRMHD solver, we use standard finite-volume methods to evolve the fluid and a constrained-transport scheme to preserve the divergence-free constraint for the magnetic field. We also utilize a first-order flux correction (FOFC) scheme to reduce the need for an artificial atmosphere and optionally enforce a maximum principle to improve robustness. We demonstrate the accuracy of AthenaK using a set of standard tests in flat and curved spacetimes. Using a SANE accretion disk around a Kerr black hole, we compare the new solver to the existing solver for stationary spacetimes using the so-called "HARM-like" formulation. We find that both formulations converge to similar results. We also include the first published binary neutron star (BNS) mergers performed on graphical processing units (GPUs). Thanks to the FOFC scheme, our BNS mergers maintain a relative error of $\mathcal{O}$(10 –11 ) or better in baryon mass conservation up to collapse. Finally, we perform scaling tests of AthenaK on OLCF Frontier, where we show excellent weak scaling of ≥80% efficiency up to 32,768 GPUs and 74% up to 65,536 GPUs for a GRMHD problem in dynamical spacetimes with six levels of mesh refinement. AthenaK achieves an order-of-magnitude speedup using GPUs compared to CPUs, demonstrating that it is suitable for performing numerical relativity problems on modern exascale resources.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

NASA NTRS1993

Reconnection of vorticity lines and magnetic lines

Magnetic field and fluid vorticity share many features. First, as divergence-free vector fields they are conveniently visualized in terms of their field lines, curves that are everywhere tangent to the field. The lines indicate direction and their density indicates field strength. The question arises of the extent to which the evolution of the fields can be treated in terms of the evolution of their field lines. Newcomb (1958) derived the general conditions on the evolution of vector fields that permit the identification of field lines from one instant to the next. The equations of evolution of the vorticity field and the magnetic field fall within Newcomb's analysis. The dynamics of the flows differ between these two systems, so that geometrically similar phenomena happen in different ways in the two systems. In this paper the geometrical similarities are emphasized. Reconnection will be defined here as evolution in which it is not possible to preserve the global identification of some field lines. There is a close relation between reconnection and the topology of the vector field lines. Nontrivial topology occurs where the field has null points or there are field lines that are closed loops.

Greene, John M.↗

NASA NTRS1993

A time domain, weighted residual formulation of Maxwell's equations

A finite element model is developed and used to simulate two-dimensional electromagnetic wave propagation and scattering. The spatial discretization of the time-domain electrodynamic equations is accomplished by a Galerkin approach. The semi-discrete equations are solved explicitly using a second-order Runge-Kutta scheme. Both the electric and magnetic fields are discretized using a single grid, with the divergence-free conditions satisfied through a correction approach. Examples depicting the scattering of plane waves in 2D geometries are given to demonstrate the validity of the methodology.

Young, Jeffrey L.↗

DOE OSTI2025

GR-Athena++: General-relativistic Magnetohydrodynamics Simulations of Neutron Star Spacetimes

We present the extension of GR-Athena++ to general-relativistic magnetohydrodynamics (GRMHD) for applications to neutron star spacetimes. The new solver couples the constrained transport implementation of Athena++ to the Z4c formulation of the Einstein equations to simulate dynamical spacetimes with GRMHD using oct-tree adaptive mesh refinement. We consider benchmark problems for isolated and binary neutron star spacetimes demonstrating stable and convergent results at relatively low resolutions and without grid symmetries imposed. The code correctly captures magnetic field instabilities in nonrotating stars with total relative violation of the divergence-free constraint of 10 –16 . It handles evolutions with a microphysical equation of state and black hole formation in the gravitational collapse of a rapidly rotating star. For binaries, we demonstrate correctness of the evolution under the gravitational radiation reaction and show convergence of gravitational waveforms. We showcase the use of adaptive mesh refinement to resolve the Kelvin–Helmholtz instability at the collisional interface in a merger of magnetised binary neutron stars. GR-Athena++ shows strong scaling efficiencies above 80% in excess of 10 5 CPU cores and excellent weak scaling is shown up to ~5 × 10 5 CPU cores in a realistic production setup. GR-Athena++ allows for the robust simulation of GRMHD flows in strong and dynamical gravity with exa-scale computers.

79 ASTRONOMY AND ASTROPHYSICS↗

NASA NTRS1991

Flux-corrected transport techniques for multidimensional compressible magnetohydrodynamics

A prescription is given for conservatively integrating generalized hydromagnetic equations using flux-corrected transport (FCT) techniques. By placing the magnetic-field components at the interface locations of the finite-difference grid, the field is kept divergence-free to within machine roundoff error. The use of FCT techniques allows an integration scheme of high accuracy to be employed, while the numerical ripples associated with large dispersion errors are avoided. The method is particularly well suited for problems involving magnetohydrodynamic shocks and other discontinuities.

Devore, C. Richard↗

NASA NTRS1973

Quantitative magnetospheric models derived from spacecraft magnetometer data

Quantitative models of the external magnetospheric field were derived by making least-squares fits to magnetic field measurements from four IMP satellites. The data were fit to a power series expansion in the solar magnetic coordinates and the solar wind-dipole tilt angle, and thus the models contain the effects of seasonal north-south asymmetries. The expansions are divergence-free, but unlike the usual scalar potential expansions, the models contain a nonzero curl representing currents distributed within the magnetosphere. Characteristics of four models are presented, representing different degrees of magnetic disturbance as determined by the range of Kp values. The latitude at the earth separating open polar cap field lines from field lines closing on the dayside is about 5 deg lower than that determined by previous theoretically-derived models. At times of high Kp, additional high latitude field lines are drawn back into the tail.

Mead, G. D.↗

NASA NTRS1975

A quantitative magnetospheric model derived from spacecraft magnetometer data

The model is derived by making least squares fits to magnetic field measurements from four Imp satellites. It includes four sets of coefficients, representing different degrees of magnetic disturbance as determined by the range of Kp values. The data are fit to a power series expansion in the solar magnetic coordinates and the solar wind-dipole tilt angle, and thus the effects of seasonal north-south asymmetries are contained. The expansion is divergence-free, but unlike the usual scalar potential expansion, the model contains a nonzero curl representing currents distributed within the magnetosphere. The latitude at the earth separating open polar cap field lines from field lines closing on the day side is about 5 deg lower than that determined by previous theoretically derived models. At times of high Kp, additional high-latitude field lines extend back into the tail. Near solstice, the separation latitude can be as low as 75 deg in the winter hemisphere. The average northward component of the external field is much smaller than that predicted by theoretical models; this finding indicates the important effects of distributed currents in the magnetosphere.

Mead, G. D.↗

NASA NTRS2005

Adaptive Numerical Dissipation Control in High Order Schemes for Multi-D Non-Ideal MHD

The required type and amount of numerical dissipation/filter to accurately resolve all relevant multiscales of complex MHD unsteady high-speed shock/shear/turbulence/combustion problems are not only physical problem dependent, but also vary from one flow region to another. In addition, proper and efficient control of the divergence of the magnetic field (Div(B)) numerical error for high order shock-capturing methods poses extra requirements for the considered type of CPU intensive computations. The goal is to extend our adaptive numerical dissipation control in high order filter schemes and our new divergence-free methods for ideal MHD to non-ideal MHD that include viscosity and resistivity. The key idea consists of automatic detection of different flow features as distinct sensors to signal the appropriate type and amount of numerical dissipation/filter where needed and leave the rest of the region free from numerical dissipation contamination. These scheme-independent detectors are capable of distinguishing shocks/shears, flame sheets, turbulent fluctuations and spurious high-frequency oscillations. The detection algorithm is based on an artificial compression method (ACM) (for shocks/shears), and redundant multiresolution wavelets (WAV) (for the above types of flow feature). These filters also provide a natural and efficient way for the minimization of Div(B) numerical error.

Yee, H. C.↗

NASA NTRS1999

Constructing the Coronal Magnetic Field by Correlating Parameterized Magnetic Field Lines with Observed Coronal Plasma Structures

A method is presented for constructing the coronal magnetic field from photospheric magnetograms and observed coronal loops. A set of magnetic field lines generated from magnetogram data is parameterized and then deformed by varying the parameterized values. The coronal flux tubes associated with this field are adjusted until the correlation between the field lines and the observed coronal loops is maximized. A mathematical formulation is described which ensures that (1) the normal component of the photospheric field remains unchanged, (2) the field is given in the entire corona over an active region, (3) the field remains divergence-free, and 4electric currents are introduced into the field. It is demonstrated that a parameterization of a potential field, comprising a radial stretching of the field, can provide a match for a simple bipolar active region, AR 7999, which crossed the central meridian on 1996 November 26. The result is a non-force-free magnetic field with the Lorentz force being of the order of 10(exp -5.5) g per s(exp 2) resulting from an electric current density of 0.79 micro A per m(exp 2). Calculations show that the plasma beta becomes larger than unity at a strong non-radial currents requires low height of about 0.25 solar radii supporting the non-force-free conclusion. The presence of such strong non-radial currents requires large transverse pressure gradients fo maintain a magnetostatic atmosphere, required by the relatively persistent nature of the coronal structures observed in AR 7999. This scheme is an important tool in generating a magnetic field solution consistent with the coronal flux tube observations and the observed photospheric magnetic field.

Allen, Gary G.↗