Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “MHD”

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 469 records · Page 26
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

MHD-induced beta limits in the Large Helical Device

Using the extended-magnetohydrodynamics code, M3D-C1, we perform a systematic numerical study of the effect of externally applied heating on the achievable plasma beta in a ten field-period heliotron. Heat sources of varying intensity are applied to a vacuum magnetic field that is representative of the standard configuration of the Large Helical Device, with R 0 = 3.66 m, where R 0 is the radial position of the magnetic axis in vacuum. As the system is driven to a state that is unstable to low-n magnetohydrodynamic (MHD) modes, nonlinear mode interactions lead to the formation of chaotic magnetic fields. With sufficiently strong heating, a collapse of the electron temperature profile is observed. This demonstrates the necessity of simulating the self-consistent evolution of plasma profiles, without imposing assumptions on the structure of the magnetic field, to accurately determine transport properties in stellarator plasmas. It also highlights the value of these advanced simulation capabilities for accelerating the development of high-performance stellarator operating scenarios.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2012

UNDERSTANDING SDO/AIA OBSERVATIONS OF THE 2010 JUNE 13 EUV WAVE EVENT: DIRECT INSIGHT FROM A GLOBAL THERMODYNAMIC MHD SIMULATION

In this work, we present a comprehensive observation and modeling analysis of the 2010 June 13 extreme-ultraviolet (EUV) wave observed by the Atmospheric Imaging Assembly (AIA) aboard the Solar Dynamics Observatory (SDO). Due to extreme advances in cadence, resolution, and bandpass coverage in the EUV regime, the AIA instrument offers an unprecedented ability to observe the dynamics of large-scale coronal wave-like transients known as EUV waves. To provide a physical analysis and further complement observational insight, we conduct a three-dimensional, time-dependent thermodynamic MHD simulation of the eruption and associated EUV wave, and employ forward modeling of EUV observables to compare the results directly observations. We focus on two main aspects: (1) the interpretation of the stark thermodynamic signatures in the multi-filter AIA data within the propagating EUV wave front, and (2) an in-depth analysis of the simulation results and their implication with respect to EUV wave theories. Multiple aspects, including the relative phases of perturbed variables, suggest that the outer, propagating component of the EUV transient exhibits the behavior of a fast-mode wave. We also find that this component becomes decoupled from the evolving structures associated with the coronal mass ejection that are also visible, providing a clear distinction between wave and non-wave mechanisms at play.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2013

MAGNETOHYDRODYNAMIC WAVES AND CORONAL HEATING: UNIFYING EMPIRICAL AND MHD TURBULENCE MODELS

We present a new global model of the solar corona, including the low corona, the transition region, and the top of the chromosphere. The realistic three-dimensional magnetic field is simulated using the data from the photospheric magnetic field measurements. The distinctive feature of the new model is incorporating MHD Alfven wave turbulence. We assume this turbulence and its nonlinear dissipation to be the only momentum and energy source for heating the coronal plasma and driving the solar wind. The difference between the turbulence dissipation efficiency in coronal holes and that in closed field regions is because the nonlinear cascade rate degrades in strongly anisotropic (imbalanced) turbulence in coronal holes (no inward propagating wave), thus resulting in colder coronal holes, from which the fast solar wind originates. The detailed presentation of the theoretical model is illustrated with the synthetic images for multi-wavelength EUV emission compared with the observations from SDO AIA and STEREO EUVI instruments for the Carrington rotation 2107.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2013

FORCED FIELD EXTRAPOLATION: TESTING A MAGNETOHYDRODYNAMIC (MHD) RELAXATION METHOD WITH A FLUX-ROPE EMERGENCE MODEL

We undertake an attempt to reconstruct the Sun's non-force-free magnetic field. The solar corona is often considered to be magnetohydrostatic. We solve the full MHD equations with a semi-realistic atmosphere model to attain this stationary state. Our method is tested with a Sun-like model which simulates the emergence of a magnetic flux rope passing from below the photosphere into the corona. Detailed diagnostics shows that our method can model the forced field more successfully than the optimization and potential method, but it still needs to be applied to real data.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2013

VELOCITY-SHEAR-INDUCED MODE COUPLING IN THE SOLAR ATMOSPHERE AND SOLAR WIND: IMPLICATIONS FOR PLASMA HEATING AND MHD TURBULENCE

We analytically consider how velocity shear in the corona and solar wind can cause an initial Alfven wave to drive up other propagating signals. The process is similar to the familiar coupling into other modes induced by non-WKB refraction in an inhomogeneous plasma, except here the refraction is a consequence of velocity shear. We limit our discussion to a low-beta plasma, and ignore couplings into signals resembling the slow mode. If the initial Alfven wave is propagating nearly parallel to the background magnetic field, then the induced signals are mainly a forward-going (i.e., propagating in the same sense as the original Alfven wave) fast mode, and a driven signal propagating like a forward-going Alfven wave but polarized like the fast mode; both signals are compressive and subject to damping by the Landau resonance. For an initial Alfven wave propagating obliquely with respect to the magnetic field, the induced signals are mainly forward- and backward-going fast modes, and a driven signal propagating like a forward-going Alfven wave but polarized like the fast mode; these signals are all compressive and subject to damping by the Landau resonance. A backward-going Alfven wave, thought to be important in the development of MHD turbulence, is also produced, but it is very weak. However, we suggest that for oblique propagation of the initial Alfven wave the induced fast-polarized signal propagating like a forward-going Alfven wave may interact coherently with the initial Alfven wave and distort it at a strong-turbulence-like rate.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2014

Relativistic MHD simulations of poynting flux-driven jets

Relativistic, magnetized jets are observed to propagate to very large distances in many active galactic nuclei (AGNs). We use three-dimensional relativistic MHD simulations to study the propagation of Poynting flux-driven jets in AGNs. These jets are already assumed to be being launched from the vicinity (∼10{sup 3} gravitational radii) of supermassive black holes. Jet injections are characterized by a model described in Li et al., and we follow the propagation of these jets to ∼parsec scales. We find that these current-carrying jets are always collimated and mildly relativistic. When α, the ratio of toroidal-to-poloidal magnetic flux injection, is large the jet is subject to nonaxisymmetric current-driven instabilities (CDI) which lead to substantial dissipation and reduced jet speed. However, even with the presence of instabilities, the jet is not disrupted and will continue to propagate to large distances. We suggest that the relatively weak impact by the instability is due to the nature of the instability being convective and the fact that the jet magnetic fields are rapidly evolving on Alfvénic time scales. We present the detailed jet properties and show that far from the jet launching region, a substantial amount of magnetic energy has been transformed into kinetic energy and thermal energy, producing a jet magnetization number σ < 1. In addition, we have also studied the effects of a gas pressure supported 'disk' surrounding the injection region, and qualitatively similar global jet behaviors were observed. We stress that jet collimation, CDIs, and the subsequent energy transitions are intrinsic features of current-carrying jets.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2014

Basic Bell-MHD turbulence

Nonresonant current instability was identified by Bell as an important mechanism for magnetic field amplification in supernova remnants. In this paper we focus on studying the nonlinear stage of this instability using the incompressible MHD formulation. We demonstrate that the evolution of magnetic turbulence driven by the Bell instability resembles turbulence driven on large scales. More importantly, we demonstrate that the energy-containing scale for magnetic fields is proportional to the square root of the magnetic energy density. Given the observational constraints of the possible field amplification, this new relation allows us to directly estimate the maximum energy of particles scattered by such fields, and this estimate is normally below the average particle energy. This implies that, without taking into account the feedback to cosmic rays (CRs), the typical scales of Bell fields, in either the linear or nonlinear regime, will be too small to affect high-energy particle acceleration. We mention several scenarios of back reaction to CRs that could be important.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2014

Jet formation in GRBs: a semi-analytic model of MHD flow in Kerr geometry with realistic plasma injection

We construct a semi-analytic model for magnetohydrodynamic (MHD) flows in Kerr geometry that incorporates energy loading via neutrino annihilation on magnetic field lines threading the horizon. We compute the structure of the double-flow established in the magnetisphere for a wide range of energy injection rates and identify the different operation regimes. At low injection rates, the outflow is powered by the spinning black hole via the Blandford-Znajek mechanism, whereas at high injection rates, it is driven by the pressure of the plasma deposited on magnetic field lines. In the intermediate regime, both processes contribute to the outflow formation. The parameter that quantifies the load is the ratio of the net power injected below the stagnation radius and the maximum power that can be extracted magnetically from the black hole.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2014

Modeling MHD accretion-ejection: episodic ejections of jets triggered by a mean-field disk dynamo

We present MHD simulations exploring the launching, acceleration, and collimation of jets and disk winds. The evolution of the disk structure is consistently taken into account. Extending our earlier studies, we now consider the self-generation of the magnetic field by an α{sup 2}Ω mean-field dynamo. The disk magnetization remains on a rather low level, which helps to evolve the simulations for T > 10, 000 dynamical time steps on a domain extending 1500 inner disk radii. We find the magnetic field of the inner disk to be similar to the commonly found open field structure, favoring magneto-centrifugal launching. The outer disk field is highly inclined and predominantly radial. Here, differential rotation induces a strong toroidal component, which plays a key role in outflow launching. These outflows from the outer disk are slower, denser, and less collimated. If the dynamo action is not quenched, magnetic flux is continuously generated, diffuses outward through the disk, and fills the entire disk. We have invented a toy model triggering a time-dependent mean-field dynamo. The duty cycles of this dynamo lead to episodic ejections on similar timescales. When the dynamo is suppressed as the magnetization falls below a critical value, the generation of the outflows and also accretion is inhibited. The general result is that we can steer episodic ejection and large-scale jet knots by a disk-intrinsic dynamo that is time-dependent and regenerates the jet-launching magnetic field.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

A steady-state picture of solar wind acceleration and charge state composition derived from a global wave-driven MHD model

The higher charge states found in slow (<400 km s{sup −1}) solar wind streams compared to fast streams have supported the hypothesis that the slow wind originates in closed coronal loops and is released intermittently through reconnection. Here we examine whether a highly ionized slow wind can also form along steady and open magnetic field lines. We model the steady-state solar atmosphere using the Alfvén Wave Solar Model (AWSoM), a global MHD model driven by Alfvén waves, and apply an ionization code to calculate the charge state evolution along modeled open field lines. This constitutes the first charge state calculation covering all latitudes in a realistic magnetic field. The ratios O{sup +7}/O{sup +6} and C{sup +6}/C{sup +5} are compared to in situ Ulysses observations and are found to be higher in the slow wind, as observed; however, they are underpredicted in both wind types. The modeled ion fractions of S, Si, and Fe are used to calculate line-of-sight intensities, which are compared to Extreme-ultraviolet Imaging Spectrometer (EIS) observations above a coronal hole. The agreement is partial and suggests that all ionization rates are underpredicted. Assuming the presence of suprathermal electrons improved the agreement with both EIS and Ulysses observations; importantly, the trend of higher ionization in the slow wind was maintained. The results suggest that there can be a sub-class of slow wind that is steady and highly ionized. Further analysis shows that it originates from coronal hole boundaries (CHBs), where the modeled electron density and temperature are higher than inside the hole, leading to faster ionization. This property of CHBs is global and observationally supported by EUV tomography.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

3D MHD SIMULATION OF FLARE SUPRA-ARCADE DOWNFLOWS IN A TURBULENT CURRENT SHEET MEDIUM

Supra-arcade downflows (SADs) are sunward, generally dark, plasma density depletions originated above posteruption flare arcades. In this paper, using 3D MHD simulations we investigate whether the SAD cavities can be produced by a direct combination of the tearing mode and Kelvin–Helmholtz instabilities leading to a turbulent current sheet (CS) medium or if the current sheet is merely the background where SADs are produced, triggered by an impulsive deposition of energy. We find that to give an account of the observational dark lane structures an addition of local energy, provided by a reconnection event, is required. We suggest that there may be a closed relation between characteristic SAD sizes and CS widths that must be satisfied to obtain an observable SAD.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

Observational diagnostics of self-gravitating MHD turbulence in giant molecular clouds

We study the observable signatures of self-gravitating magnetohydrodynamics (MHD) turbulence by applying the probability density functions (PDFs) and the spatial density power spectrum to synthetic column density maps. We find that there exists three characterizable stages of the evolution of the collapsing cloud which we term “early,” “intermediate,” and “advanced.” At early times, i.e., t<0.15t{sub ff}, the column density has a power spectral slope similar to nongravitating supersonic turbulence and a lognormal distribution. At an intermediate stage, i.e., 0.15t{sub ff} 0.35t{sub ff}, the power spectral slope is positive valued, and a dramatic increase is observed in the PDF moments and the Tsallis incremental PDF parameters, which gives rise to deviations between PDF-sonic Mach number relations. Finally, we show that the imprint of gravity on the density power spectrum can be replicated in non-gravitating turbulence by introducing a delta-function with amplitude equivalent to the maximum valued point in a given self-gravitating map. We find that the turbulence power spectrum restored through spatial filtering of the high density material.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

MHD MODELING OF THE OUTER HELIOSPHERIC STRUCTURES AROUND THE HELIOPAUSE

Voyager 1 (V1) observed the first step increase of the galactic cosmic ray (GCR) flux for energies >70 MeV and >211 MeV near the heliopause at 2012.31, just before the magnetic field polarity changed from “away” to “toward” at 2012.47. This event was not accompanied by an anomalous cosmic ray (ACR) flux change. The second GCR step increase, together with an ACR decrease, was observed at 2012.65. The magnetic field change was observed with its polarity changing from a “toward” to an apparent “away” direction at 2012.56, at ∼122 AU from the Sun, with an eventual orientation of 22°. The later magnetic field has been identified as the interstellar magnetic field and the current sheet at the polarity change is identified as the boundary of the interstellar magnetic field. We present a three-dimensional stationary MHD simulation with a heliopause located at ∼122 AU and a magnetic field polarity change across the heliopause. The apparent “away” polarity of the magnetic field is confirmed to be that of the interstellar magnetic field in this simulation. We present a detailed analysis of the magnetic field line topology in the vicinity of the heliopause and argue for a favorable condition of the first GCR step increase with no concurrent ACR flux change before the time corresponding to the magnetic polarity change at 2012.47.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

Synergy between fast-ion transport by core MHD and test blanket module fields in DIII-D experiments

We measured fast-ion transport caused by the combination of MHD and a mock-up test-blanket module (TBM) coil in the DIII-D tokamak. The primary diagnostic is an infrared camera that measures the heat flux on the tiles surrounding the coil. The combined effects of the TBM and four other potential sources of transport are studied: neoclassical tearing modes, Alfvén eigenmodes, sawteeth, and applied resonant magnetic perturbation fields for the control of edge localized modes. A definitive synergistic effect is observed at sawtooth crashes where, in the presence of the TBM, the localized heat flux at a burst increases from 0.36 ± 0.27 to 2.6 ± 0.5 MW/m -2 .

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2017

Analysis of ELM stability with extended MHD models in JET, JT-60U and future JT-60SA tokamak plasmas

The stability with respect to a peeling–ballooning mode (PBM) was investigated numerically with extended MHD simulation codes in JET, JT-60U and future JT-60SA plasmas. The MINERVA-DI code was used to analyze the linear stability, including the effects of rotation and ion diamagnetic drift (${\omega }_{* {\rm{i}}}$), in JET-ILW and JT-60SA plasmas, and the JOREK code was used to simulate nonlinear dynamics with rotation, viscosity and resistivity in JT-60U plasmas. It was validated quantitatively that the ELM trigger condition in JET-ILW plasmas can be reasonably explained by taking into account both the rotation and ${\omega }_{* {\rm{i}}}$ effects in the numerical analysis. When deuterium poloidal rotation is evaluated based on neoclassical theory, an increase in the effective charge of plasma destabilizes the PBM because of an acceleration of rotation and a decrease in ${\omega }_{* {\rm{i}}}$. The difference in the amount of ELM energy loss in JT-60U plasmas rotating in opposite directions was reproduced qualitatively with JOREK. By comparing the ELM affected areas with linear eigenfunctions, it was confirmed that the difference in the linear stability property, due not to the rotation direction but to the plasma density profile, is thought to be responsible for changing the ELM energy loss just after the ELM crash. A predictive study to determine the pedestal profiles in JT-60SA was performed by updating the EPED1 model to include the rotation and ${\omega }_{* {\rm{i}}}$ effects in the PBM stability analysis. It was shown that the plasma rotation predicted with the neoclassical toroidal viscosity degrades the pedestal performance by about 10% by destabilizing the PBM, but the pressure pedestal height will be high enough to achieve the target parameters required for the ITER-like shape inductive scenario in JT-60SA.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Elimination of MHD current sheets by modifications to the plasma wall in a fixed boundary model

Models of magnetohydrodynamic (MHD) equilibia that for computational convenience assume the existence of a system of nested magnetic flux surfaces tend to exhibit singular current sheets. These sheets are located on resonant flux surfaces that are associated with rational values of the rotational transform. We study the possibility of eliminating these singularities by suitable modifications of the plasma boundary, which we prescribe in a fixed boundary setting. We find that relatively straightforward iterative procedures can be used to eliminate weak current sheets that are generated at resonant flux surfaces by the nonlinear interactions of resonating wall harmonics. These types of procedures may prove useful in the design of fusion devices with configurations that enjoy improved stability and transport properties.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Influences of δB contribution and parallel inertial term of energetic particles on MHD-kinetic hybrid simulations: a case study of the 1/1 internal kink mode

The magnetohydrodynamic-kinetic (MHD-kinetic) hybrid model (Park et al 1992 Phys. Fluids B 4 2033–7) has been widely applied in studying energetic particles (EPs) problems in fusion plasmas for past decades. The pressure-coupling scheme or the current-coupling scheme is adopted in this model. However, two noteworthy issues arise in the model application: firstly, the coupled term introduced in the pressure-coupling scheme, (∇•P h ) ⟂ , is often simplified by ∇•P h , which is equivalent to neglecting the parallel inertial term of EPs; secondly, besides the $δf$ contribution caused by changing in the EP distribution function, the magnetic field perturbation (the $δB$ contribution) generated during development of the instabilities should also be considered, but it is often ignored in existing hybrid simulations. In this paper, we derive the analytical formulations under these two coupling schemes and then numerically study the representative case of the linear stability of the $m/n$ = $1/1$ internal kink mode (IKM) (Fu et al 2006 Phys. Plasmas 13 052517) by using the CLT-K code. Further, it is found that the approximated models can still yield reasonable results when EPs are isotopically distributed. But it fails completely in cases with anisotropic EP distributions. In addition, we further investigate the influence of EP's orbit width on the stability of IKM and verify the equivalence between pressure-coupling scheme and the current-coupling scheme.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2017

Resistive MHD modelling of the quasi-single helicity state in the KTX regimes

Here, the potential formation of a quasi-single-helicity (QSH) state in the Keda Torus eXperiment (KTX) is investigated in resistive MHD simulations using the NIMROD code. We focus on the effects of finite resistivity on the mode structure and characteristics of the dominant linear and nonlinear resistive tearing-mode in a finite β, cylindrical configuration of a reversed field pinch model for KTX. In the typical resistive regimes of KTX where the Lundquist number $S=5 \times 10^4$ , the plasma transitions to a steady QSH state after evolving through an initial transient phase with multiple helicities. The dominant mode of the QSH state develops from the dominant linear tearing mode instability. In the lower β regime, the QSH state is intermittent and short in duration; in the higher β regime, the QSH state persists for a longer time and should be more easily observed in experiments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗