Search NASA⌕ Search

SEARCH · Search NASA

Results for “MAGNETOHYDRODYNAMICS”

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 577 records · Page 32
All sourcesNASA NTRSDOE OSTI
NASA NTRS2007

GRMHD/RMHD Simulations and Stability of Magnetized Spine-Sheath Relativistic Jets

A new general relativistic magnetohydrodynamics (GRMHD ) code "RAISHIN" used to simulate jet generation by rotating and non-rotating black holes with a geometrically thin Keplarian accretion disk finds that the jet develops a spine-sheath structure in the rotating black hole case. Spine-sheath structure and strong magnetic fields significantly modify the Kelvin-Helmholtz (KH) velocity shear driven instability. The RAISHIN code has been used in its relativistic magnetohydrodynamic (RMHD) configuration to study the effects of strong magnetic fields and weakly relativistic sheath motion, cl2, on the KH instability associated with a relativistic, Y = 2.5, jet spine-sheath interaction. In the simulations sound speeds up to ? c/3 and Alfven wave speeds up to ? 0.56 c are considered. Numerical simulation results are compared to theoretical predictions from a new normal mode analysis of the RMHD equations. Increased stability of a weakly magnetized system resulting from c/2 sheath speeds and stabilization of a strongly magnetized system resulting from d 2 sheath speeds is found.

Hardee, Philip↗
NASA NTRS2010

Coherent Eigenmodes in Homogeneous MHD Turbulence

The statistical mechanics of Fourier models of ideal, homogeneous, incompressible magnetohydrodynamic (MHD) turbulence is discussed, along with their relevance for dissipative magnetofluids. Although statistical theory predicts that Fourier coefficients of fluid velocity and magnetic field are zero-mean random variables, numerical simulations clearly show that certain coefficients have a non-zero mean value that can be very large compared to the associated standard deviation, i.e., we have coherent structure. We use eigenanalysis of the modal covariance matrices in the probability density function to explain this phenomena in terms of `broken ergodicity', which is defined to occur when dynamical behavior does not match ensemble predictions on very long time-scales. We provide examples from 2-D and 3-D magnetohydrodynamic simulations of homogeneous turbulence, and show new results from long-time simulations of MHD turbulence with and without a mean magnetic field

Shebalin, John V.↗
NASA NTRS2013

High-Speed Imaging of the First Kink Mode Instability in a Magnetoplasmadynamic Thruster

One of the biggest challenges to efficient magnetoplasmadynamic thruster (MPDT) operation is the onset of high-frequency voltage oscillations as the discharge current is increased above a threshold value. The onset regime is closely related to magnetohydrodynamic instabilities known as kink modes. This work documents direct observation of the formation and quasi-steady state behavior of an argon discharge plasma in a MPDT operating at discharge currents of 8 to 10 kA for a pulse length of approximately 4 ms. A high-speed camera images the quasi-steady-state operation of the thruster at 26,143 fps with a frame exposure time of 10 micro s. A 0.9 neutral density filter and 488-nm argon line filter with a 10-nm bandwidth are used on separate trials to capture the time evolution of the discharge plasma. Frame-by-frame analysis of the power flux incident on the CCD sensor shows both the initial discharge plasma formation process and the steady-state behavior of the discharge plasma. Light intensity levels on the order of 4-6 W/m2 indicate radial and azimuthal asymmetries in the concentration of argon plasma in the discharge channel. The plasma concentration exhibits characteristics that suggest the presence of a helical plasma column. This helical behavior has been observed in previous experiments that characterize plasma kink mode instabilities indirectly. Therefore, the direct imaging of these plasma kink modes further supports the link between MPDT onset behavior and the excitation of the magnetohydrodynamic instabilities.

Walker, Jonathan A.↗
NASA NTRS2018

Magnetosheath Propagation Time of Solar Wind Directional Discontinuities

Observed delays in the ground response to solar wind directional discontinuities have been explained as the result of larger than expected magnetosheath propagation times. Recently, Samsonov et al. (2017, https://doi.org/10.1002/2017GL075020) showed that the typical time for a southward interplanetary magnetic field (IMF) turning to propagate across the magnetosheath is 14 min. Here by using a combination of magnetohydrodynamic simulations, spacecraft observations, and analytic calculations, we study the dependence of the propagation time on solar wind parameters and near-magnetopause cutoff speed. Increases in the solar wind speed result in greater magnetosheath plasma flow velocities, decreases in the magnetosheath thickness and, as a result, decreases in the propagation time. Increases in the IMF strength result in increases in the magnetosheath thickness and increases in the propagation time. Both magnetohydrodynamic simulations and observations suggest that propagation times are slightly smaller for northward IMF turnings. Magnetosheath flow deceleration must be taken into account when predicting the arrival times of solar wind structures at the dayside magnetopause.

Samsonov, A. A.↗
NASA NTRS2018

Predicting Global Ground Geoelectric Field With Coupled Geospace and Three-Dimensional Geomagnetic Induction Models

We forecast the global effects of space weather on the geoelectric and geomagnetic fields using a novel combination of methods. We use a realistic three-dimensional (3-D) model of Earth's electrical conductivity and a realistic representation of magnetospheric and ionospheric current systems. Our scheme involves the following steps: (1) We run a global magnetohydrodynamic model of the magnetosphere coupled to an electrostatic model of the ionosphere. (2) We calculate a global time series of the ground magnetic field resulting from the ionospheric, field-aligned, and magnetospheric currents of the global magnetohydrodynamic model. (3) We approximate this external field by an equivalent source current flowing in a thin shell above Earth. (4) We calculate a global time series of geoelectric and geomagnetic fields from the equivalent current and a 3-D conductivity model of Earth that also takes into account the coast effect due to large horizontal conductivity gradient. We verify our implementation by comparing the results against known analytic and numeric solutions, and then apply our scheme to the geomagnetic storm of 14 and 15 December 2006. In particular, we show that accounting for 3-D structure of Earth's conductivity results in significantly enhanced geoelectric field at large lateral gradients of conductivity, especially in coastal regions, both at middle and high latitudes. In the studied geomagnetic storm the largest values of 3-D geoelectric field are detected at high latitudes reaching 2.5 volts per kilometer and the 3-D effect extends inland by a few hundred kilometers.

Honkonen, I.↗
NASA NTRS2019

Contribution of ULF Wave Activity to the Global Recovery of the Outer Radiation Belt During the Passage of a High‐Speed Solar Wind Stream Observed in September 2014

Energy coupling between the solar wind and the Earth's magnetosphere can affect the electron population in the outer radiation belt. However, the precise role of different internal and external mechanisms that leads to changes of the relativistic electron population is not entirely known. This paper describes how ultralow frequency (ULF) wave activity during the passage of Alfvenic solar wind streams contributes to the global recovery of the relativistic electron population in the outer radiation belt. To investigate the contribution of the ULF waves, we searched the Van Allen Probes data for a period in which we can clearly distinguish the enhancement of electron fluxes from the background. We found that the global recovery that started on 22 September 2014, which coincides with the corotating interaction region preceding a high‐speed stream and the occurrence of persistent substorm activity, provides an excellent scenario to explore the contribution of ULF waves. To support our analyses, we employed ground‐ and space‐based observational data and global magnetohydrodynamic simulations and calculated the ULF wave radial diffusion coefficients employing an empirical model. Observations show a gradual increase of electron fluxes in the outer radiation belt and a concomitant enhancement of ULF activity that spreads from higher to lower L‐shells. Magnetohydrodynamic simulation results agree with observed ULF wave activity in the magnetotail, which leads to both fast and Alfven modes in the magnetospheric nightside sector. The observations agree with the empirical model and are confirmed by phase space density calculations for this global recovery period.

Silva, L. A. Da↗
NASA NTRS2022

Quantifying Errors in 3D CME Parameters Derived from Synthetic Data Using White-Light Reconstruction Techniques

Current efforts in space weather forecasting of CMEs have been focused on predicting their arrival time and magnetic structure. To make these predictions, methods have been developed to derive the true CME speed, size, position, and mass, among others. Difficulties in determining the input parameters for CME forecasting models arise from the lack of direct measurements of the coronal magnetic fields and uncertainties in estimating the CME 3D geometric and kinematic parameters after eruption. White-light coronagraph images are usually employed by a variety of CME reconstruction techniques that assume more or less complex geometries. This is the first study from our International Space Science Institute (ISSI) team “Understanding Our Capabilities in Observing and Modeling Coronal Mass Ejections”, in which we explore how subjectivity affects the 3D CME parameters that are obtained from the Graduated Cylindrical Shell (GCS) reconstruction technique, which is widely used in CME research. To be able to quantify such uncertainties, the “true” values that are being fitted should be known, which are impossible to derive from observational data. We have designed two different synthetic scenarios where the “true” geometric parameters are known in order to quantify such uncertainties for the first time. We explore this by using two sets of synthetic data: 1) Using the ray-tracing option from the GCS model software itself, and 2) Using 3D magnetohydrodynamic (MHD) simulation data from the Magnetohydrodynamic Algorithm outside a Sphere code. Our experiment includes different viewing configurations using single and multiple viewpoints. CME reconstructions using a single viewpoint had the largest errors and error ranges overall for both synthetic GCS and simulated MHD white-light data. As the number of viewpoints increased from one to two, the errors decreased by approximately 4° in latitude, 22° in longitude, 14° in tilt, and 10° in half-angle. Our results quantitatively show the critical need for at least two viewpoints to be able to reduce the uncertainty in deriving CME parameters. We did not find a significant decrease in errors when going from two to three viewpoints for our specific hypothetical three spacecraft scenario using synthetic GCS white-light data. As we expected, considering all configurations and numbers of viewpoints, the mean absolute errors in the measured CME parameters are generally significantly higher in the case of the simulated MHD white-light data compared to those from the synthetic white-light images generated by the GCS model. We found the following CME parameter error bars as a starting point for quantifying the minimum error in CME parameters from white-light reconstructions: Δθ (latitude)=6° +2° -3° , Δϕ (longitude)=11° +18° -6° , Δγ (tilt)=25° +8° -7° , Δx (half-angle)=10° +12° -6° , Δh (height)=0.6 +1.2 -0.4 R ⨀ , and Δκ (ratio)=0.1 +0.03 -0.02 .

Coronal mass ejections↗
NASA NTRS2022

Quantifying Errors in 3D CME Parameters Derived From Synthetic Data Using White-Light Reconstruction Techniques

Current efforts in space weather forecasting of CMEs have been focused on predicting their arrival time and magnetic structure. To make these predictions, methods have been developed to derive the true CME speed, size, position, and mass, among others. Difficulties in determining the input parameters for CME forecasting models arise from the lack of direct measurements of the coronal magnetic fields and uncertainties in estimating the CME 3D geometric and kinematic parameters after eruption. White-light coronagraph images are usually employed by a variety of CME reconstruction techniques that assume more or less complex geometries. This is the first study from our International Space Science Institute (ISSI) team “Understanding Our Capabilities in Observing and Modeling Coronal Mass Ejections”, in which we explore how subjectivity affects the 3D CME parameters that are obtained from the Graduated Cylindrical Shell (GCS) reconstruction technique, which is widely used in CME research. To be able to quantify such uncertainties, the “true” values that are being fitted should be known, which are impossible to derive from observational data. We have designed two different synthetic scenarios where the “true” geometric parameters are known in order to quantify such uncertainties for the first time. We explore this by using two sets of synthetic data: 1) Using the ray-tracing option from the GCS model software itself, and 2) Using 3D magnetohydrodynamic (MHD) simulation data from the Magnetohydrodynamic Algorithm outside a Sphere code. Our experiment includes different viewing configurations using single and multiple viewpoints. CME reconstructions using a single viewpoint had the largest errors and error ranges overall for both synthetic GCS and simulated MHD white-light data. As the number of viewpoints increased from one to two, the errors decreased by approximately 4° in latitude, 22° in longitude, 14° in tilt, and 10° in half-angle. Our results quantitatively show the critical need for at least two viewpoints to be able to reduce the uncertainty in deriving CME parameters. We did not find a significant decrease in errors when going from two to three viewpoints for our specific hypothetical three spacecraft scenario using synthetic GCS white-light data. As we expected, considering all configurations and numbers of viewpoints, the mean absolute errors in the measured CME parameters are generally significantly higher in the case of the simulated MHD white-light data compared to those from the synthetic white-light images generated by the GCS model. We found the following CME parameter error bars as a starting point for quantifying the minimum error in CME parameters from white-light reconstructions: Δθ (latitude)=6° +2° -3° , Δϕ (longitude)=11° +18° -6° , Δγ (tilt)=25° +8° -7° , Δx (half-angle)=10° +12° -6° , Δh (height)=0.6 +1.2 -0.4 R ⨀ , and Δκ (ratio)=0.1 +0.03 -0.02 .

Coronal mass ejections↗
DOE OSTI2024

An MPMD approach coupling electromagnetic continuum mechanics approximations in ALEGRA

In this work, two complementary approximations for describing aspects of continuum electromagnetics in moving media are discussed: electroquasistatic and magnetoquasistatic. Each has been implemented in the finite element shock code ALEGRA for modeling dynamic electromechanical phenomena on typical engineering time scales, with fully integrated circuit coupling. The approximations can be obtained by consistent asymptotic balancing of Maxwell’s equations relative to timescales associated with magnetic diffusion, charge relaxation, and electromagnetic wave propagation. In ALEGRA, the electroquasistatic approximation is used for ferroelectric (FE) modeling, while the magnetoquasistatic approximation is used for magnetohydrodynamic (MHD) modeling. In this paper we introduce for the first time a detailed derivation of a useful quasi-steady “low-R m ” variant of the MHD approximation applicable for cases, such as with detonators, where the thermodynamic pressure arising from Joule heating dominates over magnetic forces. An additional purpose of this paper is to present a coupling mode using Multiple Program-Multiple Data (MPMD) message passing communication that allows the user to run 3D FE problems together with 2D and/or 3D MHD problems with the respective simulation domains coupled through a common circuit equation. The MPMD coupling capability is used here to model the dynamic coupling of a notional ferroelectric generator with an RP-87 exploding bridgewire detonator. The simulated bridgewire heats up and bursts under current generated by simulated depoling of the ferroelectric generator, as a demonstration of the MPMD capability.

42 ENGINEERING↗
DOE OSTI2025

Construction of 3D MHD pressure drop correlation and flow characterization in the contraction region of a fusion blanket manifold

Inlet and outlet manifolds are typical components of liquid metal (LM) blanket designs of a fusion power reactor to be used to distribute the LM flow into breeding channels and collect it at the exit of the blanket. High pressure loss in the magnetohydrodynamic (MHD) flows featuring abrupt geometrical changes is one of the main feasibility issues of such designs. Recently, optimization studies were conducted to construct 3D MHD pressure drop correlations for a LM flow in an electrically insulating manifold with gradual expansion. Here, the 3D computational approach developed in that study is applied to the outlet manifold featuring gradual contraction. A systematic analysis was performed with a total number of 135 flow cases computed with COMSOL Multiphysics for Hartmann numbers 1000 < Ha < 10,000, Reynolds numbers 100 < Re < 12,000, and contraction angles 45° < θ < 75° for a fixed contraction ratio of 4. The effects of Ha, Re and θ on the flow recirculation, development length and the total pressure drop were carefully examined. A linear regression analysis was used to determine the power rule of pressure drop coefficient k related to Ha and Re, demonstrating a good match with the Ludford layer theory. Eventually, a correlation for the 3D MHD pressure drop coefficient was constructed as a function of Ha, Re and θ. Further, the results were compared against the inlet manifold. It was found that the flow in the inlet manifold exhibits larger recirculation zones. In the investigated range of Ha, Re and θ, the pressure drop coefficient k of the LM MHD flow in the gradual contraction is only slightly lower (< 8 %) than that in the gradual expansion.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

3D hybrid fluid-particle jet simulations and the importance of synchrotron radiative losses

Context. Relativistic jets in active galactic nuclei are known for their exceptional energy output, and imaging the synthetic synchrotron emission of numerical jet simulations is essential for a comparison with observed jet polarization emission. Aims. Through the use of 3D hybrid fluid-particle jet simulations (with the PLUTO code), we overcome some of the commonly made assumptions in relativistic magnetohydrodynamic (RMHD) simulations by using non-thermal particle attributes to account for the resulting synchrotron radiation. Polarized radiative transfer and ray-tracing (via the RADMC-3D code) highlight the differences in total intensity maps when (i) the jet is simulated purely with the RMHD approach, (ii) a jet tracer is considered in the RMHD approach, and (iii) a hybrid fluid-particle approach is used. The resulting emission maps were compared to the example of the radio galaxy Centaurus A. Methods. We applied the Lagrangian particle module implemented in the latest version of the PLUTO code. This new module contains a state-of-the-art algorithm for modeling diffusive shock acceleration and for accounting for radiative losses in RMHD jet simulations. The module implements the physical postulates missing in RMHD jet simulations by accounting for a cooled ambient medium and strengthening the central jet emission. Results. We find a distinction between the innermost structure of the jet and the back-flowing material by mimicking the radio emission of the Seyfert II radio galaxy Centaurus A when considering an edge-brightened jet with an underlying purely toroidal magnetic field. We demonstrate the necessity of synchrotron cooling as well as the improvements gained when directly accounting for non-thermal synchrotron radiation via non-thermal particles.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Spot size measurement of a deuterium–tritium dense plasma focus using neutron radiography

Neutron radiography is a technique uniquely suited to applications in nuclear diagnostics, non-destructive testing, and subcritical experiments. The spatial resolution of neutron radiographs is degraded by optical blur in the imaging system and the neutron source size, where the ideal source is point-like to optimize the point-spread function. A potential neutron source for radiography is the dense plasma focus (DPF), a coaxial Z-pinch that produces thermonuclear and beam-target neutrons. To assess if the source size is suitable for radiography, a neutron imaging system was used to measure the source size of the 4 MA Sodium DPF at the Nevada National Security Site operating with deuterium–tritium gas-fill. The source size was measured using the edge-spread function of tungsten objects, each having a rolled (convex) edge. The spot size was found to be 7–12 mm full-width at half-max (FWHM) assuming a Gaussian source, though comparison is presented for Lorentzian and Bennett distributions. The average FWHM was found to be 8.6 ± 1.2 mm vertically and 10.8 ± 1.2 mm horizontally with respect to the image plane, averaging over varied edges and alignments. The results were sensitive to source alignment and edge metrology, which introduced notable uncertainties. These results are consistent with separate experimental measurements as well as magnetohydrodynamics simulations of this DPF, which suggest that neutron production can originate from pinches ∼5–7 mm off-axis. These results suggest that the DPF should be used for radiography at low magnification (M < 1) where spot size does not dominate spatial blur.

47 OTHER INSTRUMENTATION↗
DOE OSTI2024

Two-dimensional Thomson scattering measurements of misaligned electron density and temperature gradients and associated Biermann battery produced fields

We present optical Thomson scattering measurements of electron density and temperature in high Mach number laser-driven blast waves in homogeneous gases. Taylor–Sedov blast waves are launched in nitrogen (N2) or helium (He) at pressures between 0.4 mTorr and 10 Torr by ablating a solid plastic target with a high energy laser pulse (10 J, 1012 W/cm2). Experiments are performed at high repetition rate (1 Hz), which allows one-dimensional and two-dimensional Thomson scattering measurements over an area of several cm2 by automatically translating the scattering volume between shots. Electron temperature and density in the blast wave fronts were seen to increase with increasing background gas pressure. Measured electron density and temperature gradients were used to calculate ∂B/∂t∝∇Te×∇ne. The experimentally measured ∂B/∂t showed agreement with the magnetic field probe (B-dot) measurements, revealing that magnetic fields are generated in the observed blast waves via the Biermann battery effect. The results are compared to numerical three-dimensional collisional magnetohydrodynamic simulations performed with FLASH, and are discussed in the context of spontaneous magnetic field generation via the Biermann battery effect.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Radiatively cooled magnetic reconnection experiments driven by pulsed power

We present evidence for strong radiative cooling in a pulsed-power-driven magnetic reconnection experiment. Two aluminum exploding wire arrays, driven by a 20 MA peak current, 300 ns rise time pulse from the Z machine (Sandia National Laboratories), generate strongly driven plasma flows (MA≈7) with anti-parallel magnetic fields, which form a reconnection layer (SL≈120) at the mid-plane. The net cooling rate far exceeds the Alfvénic transit rate (τcool−1/τA−1≫1), leading to strong cooling of the reconnection layer. We determine the advected magnetic field and flow velocity using inductive probes positioned in the inflow to the layer, and inflow ion density and temperature from analysis of visible emission spectroscopy. A sharp decrease in x-ray emission from the reconnection layer, measured using filtered diodes and time-gated x-ray imaging, provides evidence for strong cooling of the reconnection layer after its initial formation. X-ray images also show localized hotspots, regions of strong x-ray emission, with velocities comparable to the expected outflow velocity from the reconnection layer. These hotspots are consistent with plasmoids observed in 3D radiative resistive magnetohydrodynamic simulations of the experiment. X-ray spectroscopy further indicates that the hotspots have a temperature (170 eV) much higher than the bulk layer (≤75 eV) and inflow temperatures (about 2 eV) and that these hotspots generate the majority of the high-energy (>1 keV) emission.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Verification of electromagnetic simulation capabilities in global gyrokinetic particle-in-cell code GTS

Recently, the numerical scheme presented by Mishchenko et al. enabled explicit gyrokinetic simulations of low-frequency electromagnetic instabilities in tokamaks at experimentally relevant values of plasma β⁠. This scheme resolved the long-standing cancellation problem that previously hindered gyrokinetic particle-in-cell code simulations of magnetohydrodynamic phenomena with inherently small parallel electric fields. Moreover, the scheme did not employ approximations that eliminate critical tearing-type instabilities. Here, we report on the implementation of this numerical scheme in the global gyrokinetic particle-in-cell code GTS. This implementation allows for a more complete and accurate picture of interaction between small scale turbulence and MHD modes in tokamaks. Additionally, we present a comprehensive set of verification simulations of numerous electromagnetic instabilities relevant to present-day tokamaks. These simulations encompass the kinetic ballooning mode, the internal kink mode, the tearing mode, the micro-tearing mode, and the toroidal Alfven eigenmode destabilized by energetic ions, which are all instrumental in understanding tokamak physics. We will also showcase the preliminary nonlinear simulations of kinetic ballooning instabilities and (2,1) island formation due to tearing mode instability. These simulations validate the accuracy of the scheme implementation and pave the way for studying how these instabilities affect plasma confinement and performance.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Transport coefficient sensitivities in a semi-analytic model for magnetized liner inertial fusion

Performance of magnetized liner inertial fusion (MagLIF) experiments is highly dependent on transport processes including magnetized heat flows and magnetic flux losses. Magnetohydrodynamic simulations used to model these experiments require a choice of model for the transport coefficients, which are the constants of proportionality relating driving terms, such as temperature gradients and currents, to the associated heat and magnetic field transport. The coefficients have been the subject of repeated recalculation using various methods throughout the years. Using a semi-analytic MagLIF model, we compare models for the transport coefficients. The choice of model modifies magnetic-flux losses caused by the Nernst thermoelectric effect and thermal conduction losses. We present simulated results from parameter scans conducted in order to compare the effects of the different models on parameters of interest in MagLIF. In some regions of parameter space, discrepancies of up to 38% are found in integrated quantities like the fusion yield. These results may serve as a guide for experimental validation of the various models, particularly as laser preheat energies and initial axial field strengths are increased on MagLIF experiments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Hall MHD waves: A fundamental departure from their MHD counterparts

It is demonstrated through a succinct derivation as to how the linear waves in Hall magnetohydrodynamics (HMHD) constitute a fundamental departure from the standard MHD waves. Apart from modifying the conventional MHD spectrum, the Hall current induces a distinct and new branch consisting of purely circularly polarized waves that may become the representative shear waves.

Astrophysics↗