Search NASA⌕ Search

SEARCH · Search NASA

Results for “MHD simulation”

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.

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

At least 163 records · Page 9
All sourcesNASA NTRSDOE OSTI
NASA NTRS2013

Characterization of the Optical and X-ray Properties of the Northwestern Wisps in the Crab Nebula

We have studied the wisps to the northwest of the Crab pulsar as part of a multi-wavelength campaign in the visible and in X-rays. Optical observations were obtained using the Nordic Optical Telescope in La Palma and X-ray observations were made with the Chandra X-ray Observatory. The observing campaign took place from October 2010 until September 2012. About once per year we observe wisps forming and peeling off from (or near) the region commonly associated with the termination shock of the pulsar wind. We find that the exact locations of the northwestern wisps in the optical and in X-rays are similar but not coincident, with X-ray wisps preferentially located closer to the pulsar. This suggests that the optical and X-ray wisps are not produced by the same particle distribution. It is also interesting to note that the optical and radio wisps are also separated from each other (Bietenholz et al. 2004). Our measurements and their implications are interpreted in terms of a Doppler-boosted ring model that has its origin in MHD modeling. While the Doppler boosting factors inferred from the X-ray wisps are consistent with current MHD simulations of PWNe, the optical boosting factors are not, and typically exceed values from MHD simulations by about a factor of 4.

Weisskopf, M. C.↗
DOE OSTI2022

Long-term general relativistic magnetohydrodynamics simulations of magnetic field in isolated neutron stars

ABSTRACT Strong magnetic fields play an important role in powering the emission of neutron stars. Nevertheless, a full understanding of the interior configuration of the field remains elusive. In this work, we present general relativistic magnetohydrodynamics (MHD) simulations of the magnetic field evolution in neutron stars lasting ${\sim } {880}\,$ms (∼6.5 Alfvén crossing periods) and up to resolutions of $0.1155\,$km using Athena++. We explore two different initial conditions, one with purely poloidal magnetic field and the other with a dominant toroidal component, and study the poloidal and toroidal field energies, the growth times of the various instability-driven oscillation modes, and turbulence. We find that the purely poloidal setup generates a toroidal field, which later decays exponentially reaching $1{{\ \rm per\ cent}}$ of the total magnetic energy, showing no evidence of reaching equilibrium. The initially stronger toroidal field setup, on the other hand, loses up to 20 per cent of toroidal energy and maintains this state till the end of our simulation. We also explore the hypothesis, drawn from previous MHD simulations, that turbulence plays an important role in the quasi-equilibrium state. An analysis of the spectra in our higher resolution setups reveals, however, that in most cases we are not observing turbulence at small scales, but rather a noisy velocity field inside the star. We also observe that the majority of the magnetic energy gets dissipated as heat increasing the internal energy of the star, while a small fraction gets radiated away as electromagnetic radiation.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1991

Large-scale gas dynamical processes affecting the origin and evolution of gaseous galactic halos

Observations of galactic halo gas are consistent with an interpretation in terms of the galactic fountain model in which supernova heated gas in the galactic disk escapes into the halo, radiatively cools and forms clouds which fall back to the disk. The results of a new study of several large-scale gas dynamical effects which are expected to occur in such a model for the origin and evolution of galactic halo gas will be summarized, including the following: (1) nonequilibrium absorption line and emission spectrum diagnostics for radiatively cooling halo gas in our own galaxy, as well the implications of such absorption line diagnostics for the origin of quasar absorption lines in galactic halo clouds of high redshift galaxies; (2) numerical MHD simulations and analytical analysis of large-scale explosions ad superbubbles in the galactic disk and halo; (3) numerical MHD simulations of halo cloud formation by thermal instability, with and without magnetic field; and (4) the effect of the galactic fountain on the galactic dynamo.

Shapiro, Paul R.↗
NASA NTRS1994

Towards an MHD theory for the standoff distance of Earth's bow shock

An MHD theory is developed for the stand-off distance a(sub s) of the bow shock and the thickness delta(sub ms) of the magnetosheath, using the empirical Spreiter et al. relation delta(sub ms) = kX and the MHD density ratio X across the shock. The theory includes as special cases the well-known gasdynamic theory and associated phenomenological MHD-like models for delta(sub ms) and a(sub s). In general, however, MHD effects produce major differences from previous models, especially at low Alfven (M(sub A)) and sonic (M(sub S)) Mach numbers. The magnetic field orientation, M(sub A), M(sub S) and the ratio of specific heats gamma are all important variables of the theory. Three principal conclusions are reached. First, the gasdynamic and phenomenological models miss important dependances on field orientation and M(sub S) and generally provide poor approximations to the MHD results. Second, changes in field orientation and M(sub S) are predicted to cause factor of approximately 4 changes in delta(sub ms) at low M(sub A). Third, using Spreiter et al.'s value for k in the MHD theory leads to maximum a(sub s) values at low M(sub A) and nominal M(sub S) that are much smaller than observations and MHD simulations require. Resolving this problem requires either the modified Spreiter-like relation and larger k found in recent MHD simulations and/or breakdown in the Spreiter-like relation at very low M(sub A).

Cairns, Iver H.↗
NASA NTRS1994

Towards an MHD Theory for the Standoff Distance of Earth's Bow Shock

A magnetohydrodynamic (MHD) theory is developed for the standoff distance a(s) of the bow shock and the thickness Delta(ms) of the magnetosheath, using the empirical Spreiter et al. relation Delta(ms) = kX and the MHD density ratio X across the shock. The theory includes as special cases the well-known gasdynamic theory and associated phenomenological MHD-like models for Delta(ms) and As. In general, however, MHD effects produce major differences from previous models, especially at low Alfev (Ma) and Sonic (Ms) Mach numbers. The magnetic field orientation Ma, Ms and the ratio of specific heats gamma are all important variables of the theory. In contrast, the fast mode Mach number need play no direct role. Three principle conclusions are reached. First the gasdynamic and phenomenological models miss important dependences of field orientation and Ms generally provide poor approximations to the MHD results. Second, changes in field orientation and Ms are predicted to cause factor of approximately 4 changes in Delta(ms) at low Ma. These effects should be important when predicting the shock's location or calculating gramma from observations. Third, using Spreiter et al.'s value for k in the MHD theory leads to maxima a(s) values at low Ma and nominal Ms that are much smaller than observations and MHD simulations require. Resolving this problem requires either the modified Spreiter-like relation and larger k found in recent MHD simulations and/or a breakdown in the Spreiter-like relation at very low Ma.

Carins, Iver H.↗
NASA NTRS2006

Plasma Circulation in the Magnetosphere

We investigate the global structure and dynamics of plasma circulation produced by prototypical solar wind disturbances of the interplanetary magnetic field and dynamic pressure. We track the global circulation and energization of solar wind, polar wind, and auroral wind plasmas throughout the magnetosphere, until they precipitate or escape into the downstream solar wind. We use the full equations of motion of the plasma ions within fields produced by a global MHD simulation of the dynamic solar wind interaction. We use the dynamic hot plasma density and Poynting energy flux specified at the inner boundary of the MHD simulation as drivers of conjugate ion outflow fluxes using local empirical relations obtained from the FAST and Polar missions. Birkeland currents computed by the MHD code are used to derive a field-parallel potential drop from a Knight-like relation [as modified by Lyons and Evans, 1980]. This potential drop is applied to each ion as an initial bulk energy, added to a thermal heating driven by the locally incident Poynting flux. The solar wind pressure increase case (B(sub Y) = 5; B(sub z) = 0 nT) produces an immediate substorm owing to compression of pre-existing plasmas. The SB(sub z), interval (embedded in NB(sub z)) produces a substorm after about one hour of development. Both disturbances enhance the auroral wind flux and heavy ion pressure of the magnetosphere substantially, with complex dynamic structuring by auroral acceleration vortexes and dynamic reconnection. Comparisons are made with observations during disturbed periods including the Halloween 2003 super-storm and other periods.

Moore, T. E.↗
NASA NTRS2014

Energetic Electrons in Dipolarization Events: Spatial Properties and Anisotropy

Using the electromagnetic fields of an MHD simulation of magnetotail reconnection, flow bursts, and dipolarization, we further investigate the acceleration of electrons to suprathermal energies. Particular emphasis is on spatial properties and anisotropies as functions of energy and time. The simulation results are compared with Time History of Events and Macroscale Interactions during Substorms observations. The test particle approach successfully reproduces several observed injection features and puts them into a context of spatial maps of the injection region(s): a dominance of perpendicular anisotropies farther down the tail and closer to the equatorial plane, an increasing importance of parallel anisotropy closer to Earth and at higher latitudes, a drop in energy fluxes at energies below approximately 10 keV, coinciding with the plasma density drop, together with increases at higher energy, a triple peak structure of flux increases near 0 deg, 90 deg, and 180 deg, and a tendency of flux increases to extend to higher energy closer to Earth and at lower latitudes. We identified the plasma sheet boundary layers and adjacent lobes as a main source region for both increased and decreased energetic electron fluxes, related to the different effects of adiabatic acceleration at high and low energies. The simulated anisotropies tend to exceed the observed ones, particularly for perpendicular fluxes at high energies. The most plausible reason is that the MHD simulation lacks the effects of anisotropy-driven microinstabilities and waves, which would reduce anisotropies.

anisotropy↗
NASA NTRS2015

The Role of the Hall Effect in Global Structure and Dynamics of Planetary Magnetospheres: Ganymede as a Case Study

We present high-resolution Hall MHD simulations of Ganymede's magnetosphere demonstrating that Hall electric fields in ion-scale magnetic reconnection layers have significant global effects not captured in resistive MHD simulations. Consistent with local kinetic simulations of magnetic reconnection, our global simulations show the development of intense field-aligned currents along the magnetic separatrices. These currents extend all the way down to the moon's surface, where they may contribute to Ganymede's aurora. Within the magnetopause and magnetotail current sheets, Hall J x B forces accelerate ions to the local Alfven speed in the out-of-plane direction, producing a global system of ion drift belts that circulates Jovian magnetospheric plasma throughout Ganymede's magnetosphere. We discuss some observable consequences of these Hall-induced currents and ion drifts: the appearance of a sub-Jovian 'double magnetopause' structure, an Alfvenic ion jet extending across the upstream magnetopause, and an asymmetric pattern of magnetopause Kelvin-Helmholtz waves.

Magnetic Reconnection↗
NASA NTRS2001

Modeling Magnetospheric Sources

We have used global magnetohydrodynamic, simulations of the interaction between the solar wind and magnetosphere together with single particle trajectory calculations to investigate the sources of plasma entering the magnetosphere. In all of our calculations solar wind plasma primarily enters the magnetosphere when the field line on which it is convecting reconnects. When the interplanetary magnetic field has a northward component the reconnection is in the polar cusp region. In the simulations plasma in the low latitude boundary layer (LLBL) can be on either open or closed field lines. Open field lines occur when the high latitude reconnection occurs in only one cusp. In the MHD calculations the ionosphere does not contribute significantly to the LLBL for northward IMF. The particle trajectory calculations show that ions preferentially enter in the cusp region where they can be accelerated by non-adiabatic motion across the high latitude electric field. For southward IMF in the MHD simulations the plasma in the middle and inner magnetosphere comes from the inner (ionospheric) boundary of the simulation. Solar wind plasma on open field lines is confined to high latitudes and exits the tailward boundary of the simulation without reaching the plasma sheet. The LLBL is populated by both ionospheric and solar wind plasma. When the particle trajectories are included solar wind ions can enter the middle magnetosphere. We have used both the MHD simulations and the particle calculations to estimate source rates for the magnetosphere which are consistent with those inferred from observations.

Walker, Raymond J.↗
DOE OSTI2026

Nature of transonic sub-Alfvénic turbulence and density fluctuations in the near-Sun solar wind

Context. Recent Parker Solar Probe (PSP) measurements have revealed that solar wind (SW) turbulence transits from a subsonic to a transonic regime near the Sun, while remaining sub-Alfvénic. These observations call for a revision of the existing SW models, where turbulence is considered to be both subsonic and sub-Alfvénic. Aims. In this work, we introduce a new magnetohydrodynamic (MHD) model of transonic sub-Alfvénic turbulence (TsAT). Methods. We used 3D MHD simulations initialized with parameters measured by PSP to investigate the properties of the new near-Sun SW transonic turbulent regime. We then derived a reduced set of MHD equations in the transonic sub-Alfvénic limit to interpret our numerical results. Results. Our TsAT model shows that turbulence is effectively nearly incompressible (NI) and has a 2D + slab (quasi-2D) geometry not only in the subsonic limit, but also in the transonic regime, as long as it remains sub-Alfvénic, a condition essentially enforced everywhere in the heliosphere by the strong local magnetic field. These predictions are consistent with 3D MHD simulations, showing that transonic turbulence is dominated by low-frequency quasi-2D incompressible structures, while compressible fluctuations are a minor component corresponding to low-frequency slow modes and high-frequency fast modes. Conclusions. Our new TsAT model extends existing NI theories of turbulence, and is potentially relevant for the theoretical and numerical modeling of space and astrophysical plasmas, including the near-Sun SW, the solar corona, and the interstellar medium.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1988

The effect of the solenoidal condition on the numerical magnetohydrodynamic simulation of coronal dynamics

Several MHD simulations of coronal dynamics performed more than 10 years ago violated the solenoidality condition for an 'open' magnetic field topology. Using an improved code, consideration is given to the effect on the physical validity of the numerical simulation (for a representative pulse disturbance) for the case when solenoidality is deliberately violated as compared to the case when it is preserved. It is found that the error incurred in this specific case ('open' topology) in the energy density and in the plasma density profiles is rather small and, hence, does not invalidate the earlier conclusions concerning mass and wave motion.

Wu, S. T.↗
DOE OSTI2025

UNR Grant Proposal Contribution

The work of Carrier et. al. (citation needed for Matt’s effect of surface roughness on phase transitions/ETI paper) demonstrates that 1D and 2D resistive magnetohydrodynamic (MHD) simulations can reliably model exploding aluminum rods driven by megaampere currents. Close agreement between simulation results and photonic Doppler velocimetry data from the Mykonos electrothermal instability (METI-II) campaign builds confidence in predictive modelling capabilities for pulsed-power HED experiments. Furthermore, 2D MHD simulations show how machined features and micro-scale surface roughness can seed the electrothermal instability (ETI). Surface roughness was observed to reduce the time of melt by 19%, with ETI growth driving the enhanced heating of the rod surface.

36 MATERIALS SCIENCE↗
NASA NTRS2011

Enhanced Spectral Anisotropies Near the Proton-Cyclotron Scale: Possible Two-Component Structure in Hall-FLR MHD Turbulence Simulations

Recent analysis of the magnetic correlation function of solar wind fluctuations at 1 AU suggests the existence of two-component structure near the proton-cyclotron scale. Here we use two-and-one-half dimensional and three-dimensional compressible MHD models to look for two-component structure adjacent the proton-cyclotron scale. Our MHD system incorporates both Hall and Finite Larmor Radius (FLR) terms. We find that strong spectral anisotropies appear adjacent the proton-cyclotron scales depending on selections of initial condition and plasma beta. These anisotropies are enhancements on top of related anisotropies that appear in standard MHD turbulence in the presence of a mean magnetic field and are suggestive of one turbulence component along the inertial scales and another component adjacent the dissipative scales. We compute the relative strengths of linear and nonlinear accelerations on the velocity and magnetic fields to gauge the relative influence of terms that drive the system with wave-like (linear) versus turbulent (nonlinear) dynamics.

Ghosh, Sanjoy↗
NASA NTRS2004

Fitting Flux Ropes to a Global MHD Solution: A Comparison of Techniques

Flux rope fitting (FRF) techniques are an invaluable tool for extracting information about the properties of a subclass of CMEs in the solar wind. However, it has proven difficult to assess their accuracy since the underlying global structure of the CME cannot be independently determined from the data. In contrast, large-scale MHD simulations of CME evolution can provide both a global view as well as localized time series at specific points in space. In this study we apply 5 different fitting techniques to 2 hypothetical time series derived from MHD simulation results. Independent teams performed the analysis of the events in "blind tests", for which no information, other than the time series, was provided. F rom the results, we infer the following: (1) Accuracy decreases markedly with increasingly glancing encounters; (2) Correct identification of the boundaries of the flux rope can be a significant limiter; and (3) Results from techniques that infer global morphology must be viewed with caution. In spite of these limitations, FRF techniques remain a useful tool for describing in situ observations of flux rope CMEs.

Riley, Pete↗
NASA NTRS2004

Fitting Flux Ropes to a Global MHD Solution: A Comparison of Techniques

Flux rope fitting (FRF) techniques are an invaluable tool for extracting information about the properties of a sub-class of coronal mass ejections (CMEs) in the solar wind. However, it has proven difficult to assess their accuracy since the underlying global structure of the CME cannot be independently determined from the data. In contrast, large-scale MHD simulations of CME evolution can provide both a global view as well as localized time series at specific points in space. In this study we apply five different fitting techniques to two hypothetical time series derived from MHD simulation results. Independent teams performed the analysis of the events in "blind tests", for which no information, other than time series, was provided. From the results, we infer the following: 1) Accuracy decreases markedly with increasingly glancing encounters; 2) Correct identification of the boundaries of the flux rope can be a significant limiter; 3) Results from techniques that infer global morphology must be viewed with caution. In spite of these limitations, FRF techniques remain a useful tool for describing in situ observations of flux rope CMEs.

Riley, Pete↗
NASA NTRS1989

Slow shock characteristics as a function of distance from the X-line in the magnetotail

Both particle and MHD simulations are performed to study the characteristics of slow shocks in the magnetotail. The particle simulations indicate that switch-off shocks exhibit large amplitude rotational wave trains, while magnetotail slow shocks with an intermediate Mach number M(An) less than M(c) of about 0.98 do not display such rotational wave trains. The MHD simulations show that the spontaneous reconnection process in the near-earth plasma sheet leads to the formation of a pair of slow shocks tailward of the reconnection line (X-line). The properties of slow shocks are found to vary as a function of the distance from X-line due to the formation of plasmoid. Slow shocks in most regions of magnetotail are found to be nonswitch-off shocks with M(An) of less than 0.98. The present results are used to discuss the lack of large amplitude rotational wave trains at slow shocks in the deep magnetotail.

Lee, L. C.↗
DOE OSTI2025

Burn propagation in magnetized high-yield inertial fusion

Recent experiments at the National Ignition Facility (NIF) have demonstrated ignition for the first time in an inertial confinement fusion (ICF) experiment, a major milestone allowing the possibility of high energy gain through burn propagation. Use of external magnetic fields, applied primarily to reduce thermal losses, could increase hotspot temperature and ease requirements for ignition, opening up the capsule design space for high energy gain. However, this same restriction of thermal transport has the potential to inhibit burn propagation, which is vital in the attainment of high gain. In this work, radiation-magnetohydrodynamics (MHD) simulations carried out using the code Chimera are used to investigate the effect of a pre-imposed magnetic field on ignition and burn propagation. This paper studies the propagation of burn using both an idealized planar model and in fully integrated 2D MHD simulations of an igniting NIF capsule. A study of magnetized burn propagation in the idealized planar model identifies three regimes of magnetized burn propagation: (1) thermal conduction driven; (2) alpha transport driven; and (3) fully suppressed burn. Simulations of NIF shot N210808 with an applied 40 T axial field show clear indication of burn suppression perpendicular to field lines, with rapid burn observed along field lines. Implosion shape is altered by the field, and anisotropic conduction causes significant modification to the rate of ablation during stagnation. These results highlight the fundamental changes to implosion dynamics in high-yield magnetized ICF and motivate further study to better optimize future magnetized target designs for high gain.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS2023

The Complexity of the Day-Side X-Line During Southward Interplanetary Magnetic Field

High-resolution global magnetohydrodynamics (MHD) simulations include both meso- and global-scale processes occurring at the magnetopause, which interact to determine the time-dependent orientation of the day-side x-line (DXL). This study demonstrates that the global orientation of the DXL in GAMERA global MHD simulations varies on a time scale of minutes during steady southward interplanetary magnetic field conditions. This behavior manifests in observational data when reconnection outflows indicate that the direction to the x-line is opposite to the prediction from a steady-state model of the reconnection location. Because steady-state models of the DXL do not capture dynamics that are independent of solar wind variations, particularly surface waves and flux transfer events, they represent a time-averaged state of the system.

Brandon Burkholder↗