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 415 records · Page 23
All sourcesNASA NTRSDOE OSTI
NASA NTRS2011

Multi-MW Closed Cycle MHD Nuclear Space Power Via Nonequilibrium He/Xe Working Plasma

Prospects for a low specific mass multi-megawatt nuclear space power plant were examined assuming closed cycle coupling of a high-temperature fission reactor with magnetohydrodynamic (MHD) energy conversion and utilization of a nonequilibrium helium/xenon frozen inert plasma (FIP). Critical evaluation of performance attributes and specific mass characteristics was based on a comprehensive systems analysis assuming a reactor operating temperature of 1800 K for a range of subsystem mass properties. Total plant efficiency was expected to be 55.2% including plasma pre-ionization power, and the effects of compressor stage number, regenerator efficiency and radiation cooler temperature on plant efficiency were assessed. Optimal specific mass characteristics were found to be dependent on overall power plant scale with 3 kg/kWe being potentially achievable at a net electrical power output of 1-MWe. This figure drops to less than 2 kg/kWe when power output exceeds 3 MWe. Key technical issues include identification of effective methods for non-equilibrium pre-ionization and achievement of frozen inert plasma conditions within the MHD generator channel. A three-phase research and development strategy is proposed encompassing Phase-I Proof of Principle Experiments, a Phase-II Subscale Power Generation Experiment, and a Phase-III Closed-Loop Prototypical Laboratory Demonstration Test.

Litchford, Ron J.↗
NASA NTRS2012

Test of Shi et al. Method to Infer the Magnetic Reconnection Geometry from Spacecraft Data: MHD Simulation with Guide Field and Antiparallel Kinetic Simulation

When analyzing data from an array of spacecraft (such as Cluster or MMS) crossing a site of magnetic reconnection, it is desirable to be able to accurately determine the orientation of the reconnection site. If the reconnection is quasi-two dimensional, there are three key directions, the direction of maximum inhomogeneity (the direction across the reconnection site), the direction of the reconnecting component of the magnetic field, and the direction of rough invariance (the "out of plane" direction). Using simulated spacecraft observations of magnetic reconnection in the geomagnetic tail, we extend our previous tests of the direction-finding method developed by Shi et al. (2005) and the method to determine the structure velocity relative to the spacecraft Vstr. These methods require data from four proximate spacecraft. We add artificial noise and calibration errors to the simulation fields, and then use the perturbed gradient of the magnetic field B and perturbed time derivative dB/dt, as described by Denton et al. (2010). Three new simulations are examined: a weakly three-dimensional, i.e., quasi-two-dimensional, MHD simulation without a guide field, a quasi-two-dimensional MHD simulation with a guide field, and a two-dimensional full dynamics kinetic simulation with inherent noise so that the apparent minimum gradient was not exactly zero, even without added artificial errors. We also examined variations of the spacecraft trajectory for the kinetic simulation. The accuracy of the directions found varied depending on the simulation and spacecraft trajectory, but all the directions could be found within about 10 for all cases. Various aspects of the method were examined, including how to choose averaging intervals and the best intervals for determining the directions and velocity. For the kinetic simulation, we also investigated in detail how the errors in the inferred gradient directions from the unmodified Shi et al. method (using the unperturbed gradient) depended on the amplitude of the calibration errors. For an accuracy of 3 for the maximum gradient direction, the calibration errors could be as large as 3% of reconnection magnetic field, while for the same accuracy for the minimum gradient direction, the calibration errors could only be as large as 0.03% of the reconnection magnetic field. These results suggest that the maximum gradient direction can normally be determined by the unmodified Shi et al. method, while the modified method or some other method must be used to accurately determine the minimum gradient direction. The structure velocity was found with magnitude accurate to 2% and direction accurate to within 5%.

Denton, R.↗
NASA NTRS2017

The Substorm Cycle as Reproduced by Global MHD Models

Recently, Gordeev et al. (2015) suggested a method to test global MHD models against statistical empirical data. They showed that four community-available global MHD models supported by the Community Coordinated Modeling Center (CCMC) produce a reasonable agreement with reality for those key parameters (the magnetospheric size, magnetic field, and pressure) that are directly related to the large-scale equilibria in the outer magnetosphere. Based on the same set of simulation runs, here we investigate how the models reproduce the global loading-unloading cycle. We found that in terms of global magnetic flux transport, three examined CCMC models display systematically different response to idealized2 h north then 2 h south interplanetary magnetic field (IMF) Bz variation. The LFM model shows a depressed return convection and high loading rate during the growth phase as well as enhanced return convection and high unloading rate during the expansion phase, with the amount of loaded unloaded magnetotail flux and the growth phase duration being the closest to their observed empirical values during isolated substorms. Two other models exhibit drastically different behavior. In the BATS-R-US model the plasma sheet convection shows a smooth transition to the steady convection regime after the IMF southward turning. In the Open GGCM a weak plasma sheet convection has comparable intensities during both the growth phase and the following slow unloading phase. We also demonstrate potential technical problem in the publicly available simulations which is related to post processing interpolation and could affect the accuracy of magnetic field tracing and of other related procedures.

Gordeev, E.↗
NASA NTRS2017

Electron Drift Resonance in the MHD-Coupled Comprehensive Inner Magnetosphere-Ionosphere Model

Relativistic electrons in the outer radiation belt are highly dynamic and respond to interplanetary solar wind structures interacting with the Earth's magnetic field. A known mechanism dictating electron dynamics is the drift-resonant interaction with ultralow frequency (ULF) waves. The present work simulates the ring current and radiation belt electron populations in the bounce-averaged, kinetic Comprehensive Inner Magnetosphere-Ionosphere model coupled with the Block Adaptive Tree Solar Wind Roe-type Upwind Scheme global magnetospheric magnetohydrodynamic (MHD) code using an idealized ULF wave solar wind density driver. ULF waves generated with 10 min periods (at 1.67 mHz frequencies) in the MHD model are characterized and the corresponding energization of electrons and radial transport of electron phase space density is presented. The drift-resonant electron energy is determined in the simulation and is consistent with the electron resonance conditions in dipolar magnetic fields. The present results will be an important component of understanding inner magnetospheric dynamics and how these inner magnetospheric populations interact with ULF waves resulting from interplanetary solar wind structures.

Komar, C. M.↗
NASA NTRS2017

A Method to Predict Magnetopause Expansion in Radial IMF Events by MHD Simulations

This paper presents a method for taking into account changes of solar wind parameters in the foreshock using global MHD simulations. We simulate four events with very distant subsolar magnetopause crossings that occurred during quasi-radial interplanetary magnetic field (IMF) intervals lasting from one to several hours. Using previous statistical results, we suggest that the density and velocity in the foreshock cavity decrease to approx. 60% and approx. 94% of the ambient solar wind values when the IMF cone angle falls below 50 deg. This diminishes the solar wind dynamic pressure to 53% and causes a corresponding magnetospheric expansion. We change the upstream solar wind parameters in a global MHD model to take these foreshock effects into account. We demonstrate that the modified model predicts magnetopause distances during radial IMF intervals close to those observed by THEMIS. The strong total pressure decrease in the data seems to be a local, rather than a global, phenomenon. Although the simulations with decreased solar wind pressure generally reproduce the observed total pressure in the magnetosheath well, the total pressure in the magnetosphere often agrees better with results for nonmodified boundary conditions. The last result reveals a limitation of our method: we changed the boundary conditions along the whole inflow boundary, although a more correct approach would be to vary parameters only in the foreshock. A model with the suggested global modification of the boundary conditions better predicts the location of part of the magnetopause behind the foreshock but may fail in predicting the rest of the magnetopause.

Samsonov, A. A.↗
NASA NTRSDate not supplied

Observations and 3D MHD Modeling of Small-scale Solar Prominence Oscillations

Small scaleoscillations were recently observed in prominences with high resolution Hinode/SOT and IRIS SJI in unprecedented detailand identified as fast magnetosonic waves. The high-resolution observations provide evidence of flows of cool material, fluid instabilities such as Kelvin-Helmholtz and Rayleigh-Taylor, and the spectroscopic information provides evidence of Doppler shifts, heating, and cooling of the prominence material. Recently nonlinear fast magnetosonic waves were observed withperiods in the range of 5–11 minutes and wavelengths in the plane of the sky (POS) of ~2000kmand flows, and the flows in narrow threads withPOS speed ∼16–46 km/s. Thenonlinear fast magnetosonic waves and flows were modeled using 2.5D MHD (Ofman & Kucera 2020). Here, we extend the studyof the nonlinear fast magnetosonic waves using more realistic 3D MHD model. The model includes the 3D structure of the density and the magnetic fieldas well as mode coupling anddemonstrate the propagating density compressions associated with the nonlinear fast magnetosonicwaves, as well as the guided modesin the pillar. The 3DMHD modelingimprove understanding ofthe POSeffects in limb prominence oscillations observations.The results of the modeling are useful for coronal seismology of the prominence structure.

Leon Ofman↗
NASA NTRSDate not supplied

Effect of Plasma Sheaths on Earth Re-entry MHD Processes

The first study of the full coupling between the aerothermodynamics, the magnetohydrodynamics, and the plasma sheaths within earth entry flows is here per-formed. The problem tackled herein is representative of a force-generating MHD patch located between the plenum and the frustum of a capsule entering the earth’s atmosphere at Mach 34. The reactions are obtained from a modified 11-species Park chemical solver. The physical model fully couples the drift-diffusion model for the sheaths to the multispecies Navier-Stokes equations for the plasma flow. The Hall and ion slip effects are taken into consideration within the plasma flow and within the sheaths. The effect of the electrode material on the MHD process is studied. Using thoriated tungsten instead of graphite leads to a four-fold increase in the Lorentz forces and also leads to significantly reduced heat fluxes on the cathode. This is attributed to the much higher electrical conductivity of the thoriated tungsten sheath reducing by orders of magnitude the plasma electrical resistance near the surfaces.

Bernard Parent↗
NASA NTRSDate not supplied

3D Radiative MHD Modeling of the Solar Atmospheric Dynamics and Structure

Dramatic dynamical phenomena accompanied by strong thermodynamic and magnetic structuring are the primary drivers of great interest in studying the solar atmosphere with high spatial and temporal resolutions. Using current computational capabilities, it has become possible to model the magnetized solar plasma in different regimes with a high degree of realism. To study the fine structuring of the solar atmosphere and dynamics, we use 3D MHD radiative models covering all layers from the upper convection zone to the corona. Realistic 3D radiative MHD modeling of the solar magnetoconvection and atmosphere allows us to generate synthetic observables that directly link the physical properties of the solar plasma to spectroscopic observables. In the presentation, we discuss qualitative and quantitative changes of the atmospheric structure and dynamics at different layers of the solar atmosphere, properties of acoustic and surface gravity waves, sources of local heating in the chromosphere-corona transition region, formation of shocks, and high-frequency oscillations in the corona, as well as manifestation of these phenomena in the modeled observables.

SMD↗
NASA NTRS2022

Strategies for Determining the Cascade Rate in MHD Turbulence: Isotropy, Anisotropy, and Spacecraft Sampling

Exact laws for evaluating cascade rates, tracing back to the Kolmogorov "4/5" law, have been extended to many systems of interest including magnetohydrodynamics (MHD), and compressible flows of the magnetofluid and ordinary fluid types. It is understood that implementations may be limited by the quantity of available data and by the lack of turbulence symmetry. Assessment of the accuracy and feasibility of such third-order (or Yaglom) relations is most effectively accomplished by examining the von Kármán–Howarth equation in increment form, a framework from which the third-order laws are derived as asymptotic approximations. Using this approach, we examine the context of third-order laws for incompressible MHD in some detail. The simplest versions rely on the assumption of isotropy and the presence of a well-defined inertial range, while related procedures generalize the same idea to arbitrary rotational symmetries. Conditions for obtaining correct and accurate values of the dissipation rate from these laws based on several sampling and fitting strategies are investigated using results from simulations. The questions we address are of particular relevance to sampling of solar wind turbulence by one or more spacecraft.

Yanwen Wang↗
NASA NTRS2023

Comparing Magnetopause Predictions From Two MHD Models During A Geomagnetic Storm and A Quiet Period

Magnetopause location is an important prediction of numerical simulations of the magnetosphere, yet the models can err, either under-predicting or over-predicting the motion of the boundary. This study compares results from two of the most widely used magnetohydrodynamic (MHD) models, the Lyon–Fedder–Mobarry (LFM) model and the Space Weather Modeling Framework (SWMF), to data from the GOES 13 and 15 satellites during the geomagnetic storm on 22 June 2015, and to THEMIS A, D, and E during a quiet period on 31 January 2013. The models not only reproduce the magnetopause crossings of the spacecraft during the storm, but they also predict spurious magnetopause motion after the crossings seen in the GOES data. We investigate the possible causes of the over-predictions during the storm and find the following. First, using different ionospheric conductance models does not significantly alter predictions of the magnetopause location. Second, coupling the Rice Convection Model (RCM) to the MHD codes improves the SWMF magnetopause predictions more than it does for the LFM predictions. Third, the SWMF produces a stronger ring current than LFM, both with and without the RCM and regardless of the LFM spatial resolution. During the non-storm event, LFM predicts the THEMIS magnetopause crossings due to the southward interplanetary magnetic field better than the SWMF. Additionally, increasing the LFM spatial grid resolution improves the THEMIS predictions, while increasing the SWMF grid resolutions does not.

magnetohydrodynamics↗
DOE OSTI2017

Boride Based Electrode Materials with Enhanced Stability Under Extreme Conditions for MHD Direct Power Extraction. Final Scientific Report

Coal power currently makes up almost 40% of global electricity production. This method of extracting energy from coal can be added to any coal plant, improving its energy output. Not only would this prepare us to keep up with the ever-increasing energy demands but it could also allow us to significantly cut the greenhouse gas emissions for the coal power industry. Hafnium diboride (HfB 2 ), Zirconium diboride (ZrB 2 ), and their solid solution mixtures were studied, primarily for electrode coating applications in coal-fired magnetohydrodynamic (MHD) direct power extraction plants. MHD works by directly extracting current from superheated plasma flame (2725 °C) seeded with potassium salts. Obviously, the electrodes near this flame are subject to extreme temperatures and corrosive environments. The coating material needs to protect the electrodes from salt corrosion, oxidation, and melting, while still providing electrical conductivity. Transition metal diborides have many of the desired properties but need improvement for oxidation resistance above 1200 °C. Two main methods of oxidation resistance improvements were investigated during this research. The first was through electrochemical anodization, a process where a thin protective oxide layer is grown on the surface to prevent further detrimental oxidation at high temperatures. Several solutions, voltages, and methods were used to find an oxide layer that would protect at high temperatures. The second method was to use various additives and mixture ratios to improve densification and oxidation resistance. Ratios of HfB 2 and ZrB 2 tested were 1:4, 1:1, 4:1 and 1:1 with additives of Ta, Hf, and Zr metals, bringing the boron to metal ratio down to 1.86. Additionally, 1.8 mol% Gd 2 O 3 and LaB 6 were added to the 1:1 ratio HfB 2 -ZrB 2 . Hydrogen evolution reaction kinetics were also investigated to shed some light on the ability of these materials to resist corrosion. Aqueous corrosion tests were also performed in several solutions of different pH. The 1:1 ratio solution added with LaB 6 gave the best overall results both in high temperature and room temperature aqueous experiments. Further work will be needed to evaluate the efficacy of these materials for such applications.

01 COAL, LIGNITE, AND PEAT↗
DOE OSTI2017

On performance of Krylov smoothing for fully-coupled AMG preconditioners for VMS resistive MHD

Here, this study explores the performance and scaling of a GMRES Krylov method employed as a smoother for an algebraic multigrid (AMG) preconditioned Newton- Krylov solution approach applied to a fully-implicit variational multiscale (VMS) nite element (FE) resistive magnetohydrodynamics (MHD) formulation. In this context a Newton iteration is used for the nonlinear system and a Krylov (GMRES) method is employed for the linear subsystems. The efficiency of this approach is critically dependent on the scalability and performance of the AMG preconditioner for the linear solutions and the performance of the smoothers play a critical role. Krylov smoothers are considered in an attempt to reduce the time and memory requirements of existing robust smoothers based on additive Schwarz domain decomposition (DD) with incomplete LU factorization solves on each subdomain. Three time dependent resistive MHD test cases are considered to evaluate the method. The results demonstrate that the GMRES smoother can be faster due to a decrease in the preconditioner setup time and a reduction in outer GMRESR solver iterations, and requires less memory (typically 35% less memory for global GMRES smoother) than the DD ILU smoother.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2017

MHD Collimation Mechanism in Arched Flux Ropes Characterized Using Volumetric, Time-Dependent B -Vector Measurements

Laboratory measurements of B(x,t) in a volume enclosing portions of two arched flux ropes show flux rope collimation driven by gradients in axial current density. These measurements verify the three predictions of a proposed MHD collimation mechanism: (1) axial magnetic forces exist in current channels with spatially varying minor radius, (2) these forces can drive counterpropagating axial flows, and (3) this process collimates the flux rope. As a result, this mechanism may explain the axial uniformity of solar loops and is relevant to other systems with current channels of varying minor radius such as solar prominences and astrophysical jets.

36 MATERIALS SCIENCE↗
DOE OSTI2023

On a fully-implicit VMS-stabilized FE formulation for low Mach number compressible resistive MHD with application to MCF

This study presents the development and evaluation of a fully-implicit variational multiscale (VMS) stabilized unstructured finite element (FE) formulation for compressible magnetohydrodynamics (MHD) model, at low Mach number regime. The model describes the dynamics of a compressible conducting fluid in the low Mach number limit in the presence of electromagnetic fields and can be used to study aspects of astrophysical phenomena, important science and technology applications, and basic plasma physics phenomena. The specific applications that motivate this study are macroscopic simulations of the longer time-scale stability and disruptions of magnetic confinement fusion (MCF) devices, specifically the ITER tokamak. The discussion considers the development of the VMS FE representation, the structure of the stabilizing terms that deal with significant convective flows, the stabilization of the nearly incompressible response of the fluid flow, and the stabilization of the constraint that enforces the solenoidal involution on the magnetic field. The nonlinear discretized system is solved with scalable preconditioned Newton–Krylov iterative methods, which employs a multiphysics block preconditioning method based on approximate block factorizations and Schur complements. The study presents an evaluation of the VMS method on a 2D cartesian tearing mode instability, and illustrates the scalability of the solvers on MCF relevant problems. A set of results are also presented for longer time-scale stability and disruptions for the ITER tokamak. These include a vertical displacement event (VDE), and a (1,1) internal kink mode. Here, the formulation is demonstrated to be scalable and also reasonably robust with respect to the Lundquist number scaling.

42 ENGINEERING↗
DOE OSTI2025

Influence of initial conditions on data-driven model identification and information entropy for ideal mhd problems

Data-driven methods of model identification are able to discern governing dynamics of a system from data. Such methods are well suited to help us learn about systems with unpredictable evolution or systems with ambiguous governing dynamics given our current understanding. Many plasma problems of interest fall into these categories as there are a wide range of models that exist, however each model is only useful in a certain regime and often limited by computational complexity. To ensure data-driven methods align with theory, they must be consistent and predictable when acting on data whose governing dynamics are known. Weak Sparse Identification of Nonlinear Dynamics (WSINDy) is a recently developed data-driven method that has shown promise in learning governing dynamics from data with high noise levels [1]. This work examines how WSINDy acts on ideal MHD test problems as the initial conditions are varied and specifies limiting requirements for successful equation identification. Furthermore, it is hard to recover the governing dynamics from data that emphasize a single dominant behavior. In these low information cases, Shannon information entropy is able to pick up on the redundancies in the data that affect recoverability.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2019

An Extended MHD Study of the 16 October 2015 MMS Diffusion Region Crossing

The Magnetospheric Multiscale (MMS) mission has given us unprecedented access to high cadence particle and field data of magnetic reconnection at Earth's magnetopause. MMS first passed very near an X-line on 16 October 2015, the Burch event, and has since observed multiple X-line crossings. Subsequent 3-D particle-in-cell (PIC) modeling efforts of and comparison with the Burch event have revealed a host of novel physical insights concerning magnetic reconnection, turbulence-induced particle mixing, and secondary instabilities. In this study, we employ the Gkeyll simulation framework to study the Burch event with different classes of extended, multifluid magnetohydrodynamics (MHD), including models that incorporate important kinetic effects, such as the electron pressure tensor, with physics-based closure relations designed to capture linear Landau damping. Such fluid modeling approaches are able to capture different levels of kinetic physics in global simulations and are generally less costly than fully kinetic PIC. Herein we focus on the additional physics one can capture with increasing levels of fluid closure refinement via comparison with MMS data and existing PIC simulations. In particular, we find that the ten-moment model well captures the agyrotropic structure of the pressure tensor in the vicinity of the X-line and the magnitude of anisotropic electron heating observed in MMS and PIC simulations. However, the ten-moment model is found to have difficulty resolving the lower hybrid drift instability, which plays a fundamental role in heating and mixing electrons in the current layer.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

An ideal MHD δW stability analysis that bypasses the Newcomb equation

In the following work, we demonstrate the efficacy of a Riccati δW ideal MHD stability analysis that bypasses the numerically intractable integration of the Newcomb equation. By transforming the linear Newcomb equation into a quadratic Riccati equation, an accurate and equivalent δW analysis is performed that is shown to enjoy some numerical advantages. Additionally, we demonstrate that the Riccati approach is better conditioned than its Newcomb counterpart at the magnetic axis, and we apply dynamical systems insights to examine its behavior at singular points. We further discuss the constraints involved in solving for the extremal admissible perturbations of a δW analysis.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

MHD mode tracking using high-speed cameras and deep learning

Abstract We present a new algorithm to track the amplitude and phase of rotating magnetohydrodynamic (MHD) modes in tokamak plasmas using high speed imaging cameras and deep learning. This algorithm uses a convolutional neural network (CNN) to predict the amplitudes of the n = 1 sine and cosine mode components using solely optical measurements from one or more cameras. The model was trained and tested on an experimental dataset consisting of camera frame images and magnetic-based mode measurements from the High Beta Tokamak - Extended Pulse (HBT-EP) device, and it outperformed other, more conventional, algorithms using identical image inputs. The effect of different input data streams on the accuracy of the model’s predictions is also explored, including using a temporal frame stack or images from two cameras viewing different toroidal regions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗