Search NASA⌕ Search

SEARCH · Search NASA

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 181 records · Page 10

Broken Symmetry and Coherent Structure in MHD Turbulence

Absolute equilibrium ensemble theory for ideal homogeneous magnetohydrodynamic (MHD) turbulence is fairly well developed. Theory and Simulation indicate that ideal MHD turbulence non-ergodic and contains coherent structure. The question of applicability real (i.e., dissipative) MHD turbulence is examined. Results from several very long time numerical simulations on a 64(exp 3) grid are presented. It is seen that coherent structure begins to form before decay dominates over nonlinearity. The connection with inverse spectral cascades and selective decay will also be discussed.

Shebalin, John V.↗

The Statistical Mechanics of Ideal MHD Turbulence

Turbulence is a universal, nonlinear phenomenon found in all energetic fluid and plasma motion. In particular. understanding magneto hydrodynamic (MHD) turbulence and incorporating its effects in the computation and prediction of the flow of ionized gases in space, for example, are great challenges that must be met if such computations and predictions are to be meaningful. Although a general solution to the "problem of turbulence" does not exist in closed form, numerical integrations allow us to explore the phase space of solutions for both ideal and dissipative flows. For homogeneous, incompressible turbulence, Fourier methods are appropriate, and phase space is defined by the Fourier coefficients of the physical fields. In the case of ideal MHD flows, a fairly robust statistical mechanics has been developed, in which the symmetry and ergodic properties of phase space is understood. A discussion of these properties will illuminate our principal discovery: Coherent structure and randomness co-exist in ideal MHD turbulence. For dissipative flows, as opposed to ideal flows, progress beyond the dimensional analysis of Kolmogorov has been difficult. Here, some possible future directions that draw on the ideal results will also be discussed. Our conclusion will be that while ideal turbulence is now well understood, real turbulence still presents great challenges.

Shebalin, John V.↗

The Effect of Magnetohydrodynamic (MHD) Energy Bypass on Specific Thrust for a Supersonic Turbojet Engine

This paper describes the preliminary results of a thermodynamic cycle analysis of a supersonic turbojet engine with a magnetohydrodynamic (MHD) energy bypass system that explores a wide range of MHD enthalpy extraction parameters. Through the analysis described here, it is shown that applying a magnetic field to a flow path in the Mach 2.0 to 3.5 range can increase the specific thrust of the turbojet engine up to as much as 420 N/(kg/s) provided that the magnitude of the magnetic field is in the range of 1 to 5 Tesla. The MHD energy bypass can also increase the operating Mach number range for a supersonic turbojet engine into the hypersonic flight regime. In this case, the Mach number range is shown to be extended to Mach 7.0.

Benyo, Theresa L.↗

Annular MHD Physics for Turbojet Energy Bypass

The use of annular Hall type MHD generator/accelerator ducts for turbojet energy bypass is evaluated assuming weakly ionized flows obtained from pulsed nanosecond discharges. The equations for a 1-D, axisymmetric MHD generator/accelerator are derived and numerically integrated to determine the generator/accelerator performance characteristics. The concept offers a shockless means of interacting with high speed inlet flows and potentially offers variable inlet geometry performance without the complexity of moving parts simply by varying the generator loading parameter. The cycle analysis conducted iteratively with a spike inlet and turbojet flying at M = 7 at 30 km altitude is estimated to have a positive thrust per unit mass flow of 185 N-s/kg. The turbojet allowable combustor temperature is set at an aggressive 2200 deg K. The annular MHD Hall generator/accelerator is L = 3 m in length with a B(sub r) = 5 Tesla magnetic field and a conductivity of sigma = 5 mho/m for the generator and sigma= 1.0 mho/m for the accelerator. The calculated isentropic efficiency for the generator is eta(sub sg) = 84 percent at an enthalpy extraction ratio, eta(sub Ng) = 0.63. The calculated isentropic efficiency for the accelerator is eta(sub sa) = 81 percent at an enthalpy addition ratio, eta(sub Na) = 0.62. An assessment of the ionization fraction necessary to achieve a conductivity of sigma = 1.0 mho/m is n(sub e)/n = 1.90 X 10(exp -6), and for sigma = 5.0 mho/m is n(sub e)/n = 9.52 X 10(exp -6).

Schneider, Steven J.↗

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↗

Time-discretization of a plasma-neutral MHD model with a semi-implicit leapfrog algorithm

The semi-implicit leapfrog time-discretization is a workhorse algorithm for initial-value MHD codes to bridge between vastly separated time scales. Inclusion of atomic interactions with neutrals breaks the functional structure of the MHD equations that exploited by the leapfrog. In this work, we address how to best integrate atomic physics into the semi-implicit leapfrog. Following the Crank-Nicolson method, one approach is to time-center the atomic interactions in the linear solver and use a Newton method to include the nonlinear contributions. Alternatively, another family of methods are based on operator-splitting the terms associated with the atomic interactions using a Strang-splitting technique. These methods naturally break equations into constituent ODE and PDE parts and preserve the structure exploited by the semi-implicit leapfrog. We study the accuracy and efficiency of these methods through a battery of 0D and 1D cases and show that a second-order-in-time Douglas-Rachford inspired coupling between the ODE and PDE advances is effective in reducing the time-discretization error to be comparable to that of Crank-Nicolson with Newton iteration of the nonlinear terms. Splitting ODE and PDE parts results in independent matrix solves for each field which reduces the computational cost considerably and provides parallelization over species relative to Crank-Nicolson.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Two-temperature effects in Hall-MHD simulations of the HIT-SI experiment

A two-temperature Hall-magnetohydrodynamic (Hall-MHD) model, which evolves the electron and ion temperatures separately, is implemented in the PSI-Tet 3D MHD code and used to model plasma dynamics in the Helicity Injected Torus–Steady Inductive (HIT-SI) experiment. The two-temperature model is utilized for HIT-SI simulations in both the PSI-Tet and NIMROD codes at a number of different injector frequencies in the 14.5–68.5 kHz range. At all frequencies, the NIMROD two-temperature model results in increased toroidal current, lower chord-averaged density, higher average temperatures, outward radial shift of the current centroid, and axial symmetrization of the current centroid, relative to the single-temperature NIMROD simulations. The two-temperature PSI-Tet model illustrates similar trends, but at high frequency operation, it exhibits lower electron temperature, smaller toroidal current, and decreased axial symmetrization with respect to the single-temperature PSI-Tet model. With all models, average temperatures and toroidal currents increase with the injector frequency. Power balance and heat fluxes to the wall are calculated for the two-temperature PSI-Tet model and illustrate considerable viscous and compressive heating, particularly at high injector frequency. Parameter scans are also presented for artificial diffusivity, wall temperature, and density. Both artificial diffusivity and the density boundary condition significantly modify the plasma density profiles, leading to larger average temperatures, toroidal current, and relative density fluctuations at low densities. As a result, a low density simulation achieves sufficiently high current gain (G > 5) to generate significant volumes of closed flux lasting 1–2 injector periods.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

MHD Analysis of Dual-Coolant Lead-Lithium Blanket for Spherical Tokamak Advanced Reactor

The tritium breeding blanket is vital for future fusion power plants, with the Spherical Tokamak Advanced Reactor (STAR) project highlighting the dual-coolant lead-lithium (DCLL) design. The DCLL blanket performs shielding, energy exhaust, and tritium breeding using a lead-lithium alloy , with lithium as the breeder and lead as the neutron multiplier. It also serves as the primary coolant, with helium providing supplemental cooling. Reduced-activation ferritic/martensitic steel is used for the blanket structure. Magnetohydrodynamic (MHD) phenomena influence the liquid metal flow in a magnetic field, affecting heat transfer in the breeder affected by energetic neutrons. Understanding key flow parameters in such conditions is critical for efficient DCLL design. This study uses three-dimensional thermofluid MHD analysis with ANSYS CFX software, modified at Princeton Plasma Physics Laboratory, to simulate high Hartmann flows. The neutronics code MCNP, coupled with plasma equilibrium, provides heat source distribution. In conclusion, we examine electromagnetic interactions in adjacent fluid domains and analyze the magnetic field’s impact on flow distribution in the inboard and outboard blanket layout, using detailed mesh generation for accurate results.

DCLL↗

Low- n global ideal MHD instabilities in the CFETR baseline scenario

This article reports an evaluation on the linear ideal magnetohydrodynamic (MHD) stability of the China Fusion Engineering Test Reactor (CFETR) baseline scenario for various first-wall locations. The initial-value code NIMROD and eigen-value code AEGIS are employed in this analysis. Despite the distinctly different approaches in modeling the scrape off layer region, the dominant growth in each of the low-n (n=1–10) modes are consistent between the two codes. We report the higher-n modes are dominated by ballooning modes and localized in the pedestal region, while the lower-n modes have more prominent external kink components and broader mode profiles. The influence of the plasma-vacuum profile and wall shape are also examined using NIMROD. In the presence of a resistive wall, the low-n ideal MHD instabilities are further studied using AEGIS. For the designed first-wall location, the n=1 resistive wall mode is found to be unstable, which could be fully stabilized by uniform toroidal rotation above the 2.9% core Alfven speed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Interpretive MHD modeling of dispersive shell pellet injection for rapid shutdown in tokamaks

Dispersive shell pellet (DSP) injection is modeled with the extended-MHD code NIMROD for interpretive insight into the results of recent DIII-D DSP experiments and to explore the dynamics of an inside-out thermal quench for disruption mitigation in tokamaks. Simulations of the pre-thermal quench (TQ) phase indicate that the upper bound for the quantity of ablated carbon shell material that will not perturb the flux surfaces is in the ballpark of, but somewhat below the experimental quantity. Even below this quantity, sufficient electrons are added to the plasma by the shell material to produce significant dilution cooling before the TQ is triggered. Simulations carried through the end of the TQ have very large amplitude MHD fluctuations (δB/B>10 -2 ) at the time of the plasma current spike associated with current profile redistribution. Finally, after the plasma current spike, which is of comparable amplitude to that measured in DIII-D experiments, none of the runaway electron test-particles whose orbits are tracked throughout the simulation remain confined.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

MHD stability of negative triangularity DIII-D plasmas

Abstract Negative triangularity (NT) experiments in DIII-D point to an emergent reactor scenario free of sawteeth, endowed with benign, nondisruptive n = 2 tearing modes, which experience q min ⩾ 1 similar to the positive triangularity hybrid scenario. Plasmas exhibiting this behavior attain β N > 3 , high enough to reconsider long held views of NT stability. Ideal MHD and tearing stability analysis of well-diagnosed equilibrium reconstructions of experimental hybrid-like plasmas predict that among shape parameters, MHD stability limits are only sensitive to average triangularity. Operation is predicted to be possible at β N relevant to upcoming tokamaks and commercial NT reactors.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Validation of 3D MHD simulations of Ne/D 2 mixed shattered pellet injection in JET

Nonlinear 3D magnetohydrodynamic (MHD) modeling of shattered pellet injection (SPI) in JET tokamak plasmas is performed with the JOREK code. The study focuses on the validation of simulation result with respect to experimental observations, addressing in particular figures of merit for the efficiency of the SPI, as a technique to mitigate thermal loads by injecting radiative impurities like neon during the thermal quench (TQ) phase of disruptions. A set of JET pulses with neon/deuterium atomic mixture ratio in the shattered pellet varying from 10% to 100% provides the experimental data to be compared with. Simulations using different models for the ablation of solid fragments and the radiation of neon impurities are considered. Synthetic diagnostics are employed for a direct quantitative comparison with key JET diagnostic systems providing in particular radiation, electron temperature and density, and magnetic measurements. 3D radiation structures from MHD simulations are analyzed to estimate the toroidal asymmetry of radiation, which is difficult to measure in JET due to the presence of only two toroidally displaced multi-channel bolometry systems. JOREK simulations show that, before the TQ, the bulk of the radiation is concentrated in a toroidal region enclosing the radiating fragments that is missed by both multi-channel bolometers, suggesting a possible underestimation of the total radiation by such diagnostic systems in the pre-TQ phase. On the other hand, JOREK predicts that during the TQ the radiation becomes more toroidally symmetric, with positive implications for the heat load to plasma facing components.

JOREK↗

Diagnosis of 3D magnetic field and mode composition in MHD turbulence with Y -parameter

ABSTRACT Magnetic fields are crucial in numerous astrophysical processes within the interstellar medium (ISM). However, the detailed determination of magnetic field geometry is notoriously challenging. Based on the modern magnetohydrodynamic (MHD) turbulence theory, we introduce a novel statistical technique, the ‘Y-parameter’, to decipher the magnetic field inclination in the ISM and identify dominant turbulence modes. The Y-parameter, calculated as the ratio of anisotropies of different Stokes parameter combinations, displays contrasting trends with the mean-field inclination angle in Alfvénic and compressible turbulence modes. A Y-parameter value around 1.5 ± 0.5 provide a statistical boundary to determine the dominant MHD turbulence modes. We have discovered specific correlations between the Y-parameter value and the inclination angle that unveil the dominant turbulence mode. This methodology, when applied to future radio polarization surveys such as LOFAR and SKA, promises to significantly enhance our knowledge of 3D magnetic field in the ISM and improve our understanding of interstellar turbulence.

79 ASTRONOMY AND ASTROPHYSICS↗

Parametric dependencies of resonant layer responses across linear, two-fluid, drift-MHD regimes

Non-axisymmetric magnetic fields arising in a tokamak either by external or internal perturbations can induce complex non-ideal MHD responses in their resonant surfaces while remaining ideally evolved elsewhere. This layer response can be characterized in a linear regime by a single parameter called the inner-layer Delta, which enables outer-layer matching and the prediction of torque balance to non-linear island regimes. Here, we follow strictly one of the most comprehensive analytic treatments including two-fluid and drift MHD effects and keep the fidelity of the formulation by incorporating the numerical method based on the Riccati transformation when quantifying the inner-layer Delta. The proposed scheme reproduces not only the predicted responses in essentially all asymptotic regimes but also with continuous transitions as well as improved accuracies. In particular, the Delta variations across the inertial regimes with viscous or semi-collisional effects have been further resolved, in comparison with additional analytic solutions. The results imply greater shielding of the electromagnetic torque at the layer than what would be expected by earlier work when the viscous or semi-collisional effects can compete against the inertial effects, and also due to the intermediate regulation by kinetic Alfven wave resonances as rotation slows down. These are important features that can alter the nonaxisymmetric plasma responses including the field penetration by external fields or island seeding process in rotating tokamak plasmas.

MHD instability↗

Soft X-ray measurements of transport and MHD activity in the core and edge NSTX plasma (Final Technical Report)

A comprehensive suite of innovative diagnostic instrumentation and techniques has been developed by the Johns Hopkins collaboration for measurements of transport and MHD activity in the core and edge of the National Spherical Torus Experiment (NSTX) at the Princeton Plasma Physics Laboratory (PPPL). These include ultrafast multi-energy poloidal soft X-ray arrays, high spatial resolution multi-energy tangential soft X-ray arrays, high spectral resolution transmission grating imaging VUV spectrometers and ultrafast bolometric arrays, which measure over a broad range of spatial, spectral, and temporal scales the emission from ohmic, RF, and beam driven NSTX plasmas. The instrumentation was designed to fulfill the specific measurement needs of the Spherical Torus (ST) magnetic fusion concept and is accompanied by a comprehensive set of numerical tools for data analysis and modeling, aimed at determining as many as possible plasma parameters and at extracting as much as possible physics information from NSTX experiments. The Johns Hopkins diagnostics have been and continue to be instrumental for essentially all the NSTX experiments. In addition, our own research and experiments using these diagnostics has led over the years to fundamental conclusions for the ST fusion physics, such as the establishment of near neoclassical impurity transport in the NSTX core and edge plasma, and the critical role of the Global Alfven Eigenmode (GAE) MHD activity for the ST core electron thermal transport. The Johns Hopkins Diagnostics continue to be operational on the upgraded NSTX-U and will provide key plasma information in its upcoming high field, current and power experimental campaigns. The present report summarizes the main technical characteristics and measurement capabilities of the diagnostic instrumentation and methods developed within our research and describes a key physics result achieved during this research. Success in our research was made possible also by the excellent support and cooperation of the NSTX research team and of PPPL researchers, who are warmly acknowledged here.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Realistic Modeling of Multi-Scale MHD Dynamics of the Solar Atmosphere

Realistic 3D radiative MHD simulations open new perspectives for understanding the turbulent dynamics of the solar surface, its coupling to the atmosphere, and the physical mechanisms of generation and transport of non-thermal energy. Traditionally, plasma eruptions and wave phenomena in the solar atmosphere are modeled by prescribing artificial driving mechanisms using magnetic or gas pressure forces that might arise from magnetic field emergence or reconnection instabilities. In contrast, our 'ab initio' simulations provide a realistic description of solar dynamics naturally driven by solar energy flow. By simulating the upper convection zone and the solar atmosphere, we can investigate in detail the physical processes of turbulent magnetoconvection, generation and amplification of magnetic fields, excitation of MHD waves, and plasma eruptions. We present recent simulation results of the multi-scale dynamics of quiet-Sun regions, and energetic effects in the atmosphere and compare with observations. For the comparisons we calculate synthetic spectro-polarimetric data to model observational data of SDO, Hinode, and New Solar Telescope.

Realistic Modeling↗

Toward Hall MHD Power Generation and Drag Augmentation Using a Coaxial-Electrode Configuration During Hypersonic Entry

Planetary EDL requires that a spacecraft must reduce its relative velocity from multiple km/s to zero for safe landing. The entry phase requires slowing down through an intense hypersonic flight regime (Mach 25 to Mach 5). Magnetohydrodynamics (MHD) is a next-generation concept with a wide untapped set of technologies that can expand planetary exploration capabilities and create more robust entry vehicles by utilizing in-situ resources, increasing control authority, and reducing heat flux. There is therefore a need for an experimentally validated simulation tool for MHD technologies in hypersonic entry.

MHD drag↗

Is the Global MHD Modeling of the Magnetosphere Adequate for GIC Prediction: the May 27–28, 2017 Storm

Practical steps taken by the international community to reduce the damage to technological systems from space weather include the development of numerical models capable of real-time predictions of electromagnetic disturbances at the Earth’s surface. Here we examine the feasibility of a version of the Space Weather Modeling Framework (SWMF) global MHD simulation code similar to that used by the NOAA Space Weather Prediction Center to predict the level of geomagnetic field variability, and consequently geomagnetically induced currents (GICs). We consider the contribution of geomagnetic disturbances to the bursts of GIC in the electric power line of the Kola Peninsula during the May 27–28, 2017 storm and compare the observations with results of the global MHD model. During the maximal disturbance magnetic field variations at East Scandinavian stations become more irregular, as intense Pi3 pulsations are superposed on the magnetic bay. These pulsations are not quasi-sinusoidal waves like typical Pc5 pulsations, but they are rather a quasi-periodic sequence of magnetic impulses with time scales ~5–15 min. During this period with elevated Pi3 activity very high values of GIC were recorded (variations >100 A) in the electric power transmission line. The SWMF modeling reasonably well reproduces the global magnetospheric parameters, such as SYM-H index or cross-polar potential. However, the magnetic field variability dB/dt in the East Scandinavia predicted by the modeling has turned out to be more than order of magnitude less than that observed. Thus, the version of SWMF with the grid used by NOAA SWPC still cannot adequately predict for the May 27–28 event the fine structure of the storm/substorm—Pi3 geomagnetic disturbances, and consequently the magnitude of the GIC that they drive.

geomagnetically induced currents↗