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.

17 records

All sources NASA NTRS DOE OSTI

NASA NTRS2014

Reconstruction of a Broadband Spectrum of Alfvenic Fluctuations

Alfvenic fluctuations in the solar wind exhibit a high degree of velocities and magnetic field correlations consistent with Alfven waves propagating away and toward the Sun. Two remarkable properties of these fluctuations are the tendencies to have either positive or negative magnetic helicity (-1 less than or equal to sigma(sub m) less than or equal to +1) associated with either left- or right- topological handedness of the fluctuations and to have a constant magnetic field magnitude. This paper provides, for the first time, a theoretical framework for reconstructing both the magnetic and velocity field fluctuations with a divergence-free magnetic field, with any specified power spectral index and normalized magnetic- and cross-helicity spectrum field fluctuations for any plasma species. The spectrum is constructed in the Fourier domain by imposing two conditions-a divergence-free magnetic field and the preservation of the sense of magnetic helicity in both spaces-as well as using Parseval's theorem for the conservation of energy between configuration and Fourier spaces. Applications to the one-dimensional spatial Alfvenic propagation are presented. The theoretical construction is in agreement with typical time series and power spectra properties observed in the solar wind. The theoretical ideas presented in this spectral reconstruction provide a foundation for more realistic simulations of plasma waves, solar wind turbulence, and the propagation of energetic particles in such fluctuating fields.

waves↗

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.↗

NASA NTRS2008

The Adaptively-Refined MHD Solver (ARMS)

ARMS is a robust code for solving the single-fluid, fully compressible equations of magnetohydrodynamics. In this presentation I will describe how to access and run ARMS and how to use HelioSpace to analyze the results. The main features of the code are the following: 1. Equations are solved in either Cartesian or Spherical coordinates. For the latter, either a linear or a Log(r) coordinate can be used. 2. The latest flux-corrected-transport schemes developed at NRZ, are used to advance the convective equations. All waves and shocks are accurately captured. 3. A fully adaptive grid can be used. The code incorporates the PARAMESH technology for block-adaptation. Arbitrary adaptation criteria can be specified. 4. The code maintains a divergence-free magnetic field to machine accuracy throughout all refinement and de-refinement. 5. A uniform resistivity and thermal conduction can be included. Solution is either explicit or implicit. 6. Arbitrary boundary and initial conditions can be used in either Cartesian or spherical coordinates. 7. An accompanying visualization package, HelioSpace, is available for viewing and analyzing the native ARMS output, (without need for interpolating onto a structured grid).

Antichos, S. K.↗

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.↗

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.↗

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.↗

NASA NTRS2003

Efficient Low Dissipative High Order Schemes for Multiscale MHD Flows: Minimization of Delta * B Numerical Error - II

The generalization of a class of low-dissipative high order filter finite difference schemes for long time wave propagation of shock/turbulence/combustion compressible viscous gas dynamic flows to compressible MHD equations for structured curvilinear grids has been developed. The new scheme consists of a divergence free preserving high order spatial base scheme with a filter approach which can be divergence-free preserving depending on the type of filter operator being used, the method of applying the filter step, and the type of flow problem to be considered. Several variants of the filter approach that cater to different flow types are proposed. These filters provide a natural and efficient way for the minimization of the divergence of the magnetic field (Delta * B) numerical error in the sense that no standard divergence cleaning is required. Performance evaluation of these variants, and the key role that the proper treatment of their corresponding numerical boundary conditions can play will be illustrated. Many levels of grid refinement and detailed comparison with several commonly used compressible MHD shock-capturing schemes will be sought. For certain MHD 2-D test problems, divergence free preservation of the magnetic fields of these filter schemes has been achieved.

Sjoegreen, Bjoern↗