Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “extended MHD plasmas”

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 109 records · Page 6

Acceleration of the Fast Solar Wind by Solitary Waves in Coronal Holes

The purpose of this investigation is to develop a new model for the acceleration of the fast solar wind by nonlinear. time-dependent multidimensional MHD simulations of waves in solar coronal holes. Preliminary computational studies indicate that nonlinear waves are generated in coronal holes by torsional Alfv\'{e}n waves. These waves in addition to thermal conduction may contribute considerably to the accelerate the solar wind. Specific goals of this proposal are to investigate the generation of nonlinear solitary-like waves and their effect on solar wind acceleration by numerical 2.5D MHD simulation of coronal holes with a broad range of plasma and wave parameters; to study the effect of random disturbances at the base of a solar coronal hole on the fast solar wind acceleration with a more advanced 2.5D MHD model and to compare the results with the available observations; to extend the study to a full 3D MHD simulation of fast solar wind acceleration with a more realistic model of a coronal hole and solar boundary conditions. The ultimate goal of the three year study is to model the, fast solar wind in a coronal hole, based on realistic boundary conditions in a coronal hole near the Sun, and the coronal hole structure (i.e., density, temperature. and magnetic field geometry,) that will become available from the recently launched SOHO spacecraft.

Ofman, Leon↗

Impact of error fields and error field correction on heat fluxes in SPARC

Using a single toroidal array of coils to reduce the m,n = 2,1 resonant error field (EF) produced by the misalignment of the axisymmetric coils in SPARC can result in the enhancement of the local divertor heat fluxes. Managing high divertor heat fluxes (q ∥ $\simeq$ 10 GW m -2 ) poses a challenge for compact tokamak devices such as SPARC. The presence of non-axisymmetric magnetic field perturbations adds complexity to the problem by generating intricate 3D edge magnetic topologies that alter the heat flux distributions on the target plates. The aim of this work is to investigate the impact of the EF correction (EFC) on the heat fluxes at the divertor plates in SPARC. The MHD code M3DC1 has been used to simulate the 3D magnetic perturbations generated by the shift and tilt of several axisymmetric coils within specified tolerances, as well as from the array of EFC coils located at the midplane. Using a heuristic model that extends the concept of an axisymmetric heat flux layer to 3D plasmas, the resultant heat flux distributions is derived from magnetic footprints calculated with the MAFOT code. The results show that the EFC could either decrease or further enhance the local heat flux when used to correct the m,n = 2,1 resonant EF to enhance the core plasma performance.

3D fields↗

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↗

Indirect Magnetic Measurements on the Aerodynamic Heating Facility (AHF)

The arc jet complex at NASA Ames Research Center (ARC) investigates various thermal protection systems (TPS) and their performance in high enthalpy flow environments to improve heat shields used in atmospheric entry. The arc jet operates by inducing an arc between two electrodes and passing the heated gas through a supersonic nozzle. The arc plasma is typically a partially ionized air mixture at low temperature (< 2 eV) but operates for long periods of time (up to 30 min) at high power (< 60 MW). However, simulations using the in-house ARC Heater Simulator (ARCHeS) code predict that the arc is still unstable to MHD instabilities, primarily the kink instability, which limit the operational envelope of the facility [1]. Consequently, it is important to better characterize these instabilities to improve facility stability and extend operational ranges. Due to the high enthalpy and long timescales (i.e., minutes) arc diagnostics must be non-invasive (i.e., optical or magnetic). This work will cover the development of tri-axial magnetic hall sensors for arc jet instabilities and subsequent 3D characterization of arc motion. The latest results describing the kink instability amplitude, frequency and wavelength will be presented and represent the first measurements of this instability in a NASA arc jet facility. References: [1] Meurisse, J. (2018). 3D Unsteady Model of Arc Heater Plasma Flow Using the ARC Heater Simulator (ARCHeS), APS Gaseous Electronics Conference, 2018: APS.

Sebastian Vicente Colom↗

Energetic Electrons in Dipolarization Events: Spatial Properties and Anisotropy

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

anisotropy↗

FES 2023 Theory & Simulation Performance Target

The Report for the 2023 Theory & Simulation Performance Target (TSPT) is presented herein. The 1st, 2nd and 3rd Quarter Reports are given in Sec. I, Sec. II and Sec. III. These are unchanged from those previously submitted, except that in this final document the equations, figures and citations have been automatically renumbered. The 4th Quarter Report is given in Sec. IV. The following material presents verification calculations of several widely used three-dimensional (3D) magnetohydrodynamic (MHD) equilibrium codes in strongly shaped configurations; for both fixed and free-boundary equilibria; for the so-called vacuum case and for the case of high plasma pressure; and, for magnetic fields that have nested magnetic surfaces and for fields that have a magnetic island. The configurations studied includes the recently discovered “precise-QA” case [Phys. Rev. Lett., 128, 035001 (2002)] and other highly optimized configurations. Because of the strong (extreme) shaping, the latter of these pose challenges to the 3D MHD codes. A new method is introduced for constructing coordinates in arbitrarily shaped toroidal domains, and this can improve the robustness of 3D equilibrium calculations. Beyond the equilibrium calculation, this report demonstrates the capability of performing non-linear extended-MHD modelling of optimized stellarator equilibria. This capability allows the assessment not only of linear stability thresholds but also of the nonlinear evolution of unstable equilibria to new configurations that might be far from the original configuration, with significantly different transport properties and magnetic topology. They also enable the self-consistent determination of pressure profiles in non-integrable magnetic topology through the use of highly anisotropic thermal-transport equations. The calculations presented here represent the first demonstration of extended-MHD calculations in strongly-shaped, optimized quasi-symmetric configurations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Solar Activity Modeling: From Subgranular Dynamical Scales to the Solar Cycles

The dynamical effects of solar magnetoconvection span a wide range spatial and temporal scales that extend from the interior to the corona and from fast turbulent motions to global magnetic activity. To study the solar activity on short temporal scales (from minutes to hours), we use 3D radiative MHD simulations that allow us to investigate complex turbulent interactions that drive various phenomena, such as plasma eruptions, spontaneous formation of magnetic structures, funnel-like structures and magnetic loops in the corona, and others. In particular, we focus on multi-scale processes of energy exchange across layers of the solar interior and atmosphere, which contribute to coronal heating and eruptive dynamics. For modeling global-scale activity, we use a data assimilation approach that has demonstrated great potential for building reliable long-term forecasts of solar activity. In particular, it has been shown that the Ensemble Kalman Filter (EnKF) method applied to the Parker-Kleeorin-Ruzmakin dynamo model is capable of predicting solar activity up to one sunspot cycle ahead in time, as well as estimating the properties of the next cycle a few years before it begins. In this presentation, using the available magnetogram data, we discuss development of the methodology and forecast quality criteria (including forecast uncertainties and sources of errors). We demonstrate the influence of observational limitations on prediction accuracy, and we present the EnKF predictions of the upcoming Solar Cycle (25) based on both the sunspot number series and observed magnetic fields and discuss the uncertainties and potential of the data assimilation approach for modeling and forecasting solar activity.

Kitiashvili, I. N.↗

Transmission of Alfven waves through the earth's bow shock - Theory and observation

From both theoretical and experimental bases, the transmission of Alfven waves through the bow shock is investigated. The theory of Alfven wave transmission through fast MHD shocks is extended to all cases of incident wave vectors. Particular consideration is given to Alfven waves propagating parallel to the ambient magnetic field with field perturbations polarized in the plane formed by the ambient magnetic field and the shock normal. An analysis is also made of magnetic field and plasma data from Explorer-35 in the vicinity of the bow shock. It is suggested that hydromagnetic waves are present in all of the 14 shock crossings studied, and that in upstream regions of at least 6 crossings, predominantly Alfvenic fluctuations exist. Average amplitudes of these fluctuations are measured on either side of the shock and the enhancement is measured by comparing their levels. Theoretical and experimental findings are compared and the apparent discrepancy in amplification factors may be explained by the strong damping of any transmitted magnetoacoustic modes downstream with relatively little damping of any transmitted Alfven waves.

Hassam, A. B.↗

Properties of the mantle-like magnetotail boundary layer: GEOTAIL data compared with a mantle model

Near the tail boundary beyond about 100 R(sub e), GEOTAIL often encounters a plasma mantle-like boundary layer in which the plasma flowing tailward transitions smoothly from magnetosheath values of speed and density to much smaller values, more characteristic of the tail lobe. This boundary layer had earlier been recognized on the basis of International Sun Earth Explorer 3 (ISEE 3) measurements. GEOTAIL confirms the existence of this layer and extends the documentation on its behavior. Boundary oscillations sweep the boundary layer over the spacecraft enabling GEOTAIL to 'sound' the layer's profile of plasma parameters. The density-versus-speed correlogram of the mantle-like part -- a useful diagnostic for comparisons -- is reasonably well simulated by a 1-D, MHD slow-mode expansion fan model of the plasma mantle. Flow directions are consistent with an open plasma mantle on the northern, duskside flank of the tail, as expected for the standard magnetic merging model of mantle formation.

Siscoe, G. L.↗

Stability and Control of Burning Tokamak Plasmas with Resistive Walls (Final Report)

This research has focused on quantitative prediction of the stability, control, and equilibrium state evolution in toroidal burning plasmas. The stability of long pulse burning plasmas is highly sensitive to the physics of resonant layers in the plasma, sources of momentum and flow, kinetic effects of energetic particles, and boundary conditions at the wall, including feedback control and error fields. In ITER in particular, the low toroidal flow equilibrium state, sustained primarily by energetic alpha particles from fusion reactions, will require the consideration of all of these key elements to predict quantitatively the stability and evolution. The principal investigators on this proposal are leading experts in the relevant theoretical and computational areas, and aimed to perform computations guided by analytic modeling, to address this physics in realistic configurations. The overall goal is to understand the key physics mechanisms that describe resistive toroidal burning plasmas, surrounded by a resistive wall, under active feedback control. With the physics of the energetic ions, resonant layers, resistive wall, and toroidal momentum transport included, this study will extend from recent publications in theory and simulation of individual effects and move toward predictive modeling for burning plasmas.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Sawtooth suppression by flux pumping on HBT-EP

Abstract This study examines the mechanisms underlying sawtooth suppression in the High Beta Tokamak-Extended Pulse (HBT-EP) device. It is observed that strong-intensity sawtooth activities correlate with reduced-amplitude magnetohydrodynamics (MHD) edge modes which are identified as m / n = 3 / 1 external kink modes, while sawtooth suppression correlates with larger and saturated edge mode amplitudes. To further investigate these correlations, the plasma–wall coupling was manipulated by adjusting the positions of the conducting walls in HBT-EP. It was found that strong sawtooth events occur when the normalized wall radius b / a is within a critical value. This implies that the plasma–wall distance must be sufficiently small to ensure effective stabilization of the edge mode. Even slight differences in major radius result in significantly different discharge styles, categorized as ‘sawtoothing discharges’ and ‘sawtooth-suppressed discharges’ respectively. Through a series of mode structure analyses, we confirm the coexistence and coupling of the m / n = 1 / 1 helical core, m / n = 2 / 1 tearing mode, and m / n = 3 / 1 external kink mode during sawtooth-suppression, and that this coupling induces anomalous current broadening. Based on these findings, we conclude that sawtooth suppression in the HBT-EP tokamak is consistent with the process of magnetic flux pumping.

Physics↗

Prominence mass ejections and their effects on the corona. I - The eruptive prominence of 21 August 1973 and the surge of 4 December 1973

A previous treatment of prominence-induced coronal responses is extended to X-ray and H-alpha observations of an eruptive prominence and a surge by using a time-dependent two-dimensional single-fluid MHD computer code that neglects dissipation and radiation. The two events and their observation are described, and some physical parameters of the coronal plasma are derived. The observed coronal responses to ascending prominences are analyzed with the aid of the cited computer code. The results obtained indicate that: (1) the coronal response to an eruptive prominence may be simulated with a density-dominated pressure pulse of long duration; (2) the response to a surge may be simulated with a temperature pulse of short duration; and (3) the investigated eruptive-prominence disturbance injected a mass of approximately 10 to the 16th power g into the corona and was associated with a coronal transient having a calculated velocity of 275 km/s at a distance of 1.5 solar radii.

Smith, J. B., Jr.↗

Solar Activity Modeling: From Subgranular Dynamical Scales to the Solar Cycles

Dynamical effects of solar magnetoconvection span a wide range spatial and temporal scales that extends from the interior to the corona and from fast turbulent motions to the global-Sun magnetic activity. To study the solar activity on short temporal scales (from minutes to hours), we use 3D radiative MHD simulations that allow us to investigate complex turbulent interactions that drive various phenomena, such as plasma eruptions, spontaneous formation of magnetic structures, funnel-like structures and magnetic loops in the corona, and others. In particular, we focus on multi-scale processes of energy exchange across the different layers, which contribute to the corona heating and eruptive dynamics, as well as interlinks between different layers of the solar interior and atmosphere. For modeling the global-scale activity we use the data assimilation approach that has demonstrated great potential for building reliable long-term forecasts of solar activity. In particular, it has been shown that the Ensemble Kalman Filter (EnKF) method applied to the Parker-Kleeorin-Ruzmakin dynamo model is capable of predicting solar activity up to one sunspot cycle ahead in time, as well as estimating the properties of the next cycle a few years before it begins. In this presentation, using the available magnetogram data, we discuss development of the methodology and forecast quality criteria (including forecast uncertainties and sources of errors). We demonstrate the influence of observational limitation on the prediction accuracy. We present the EnKF predictions of the upcoming Solar Cycle 25 based on both the sunspot number series and observed magnetic fields, and discuss the uncertainties and potential of the data assimilation approach for modeling and forecasting the solar activity.

Kitiashvili, I. N.↗

Io's Plasma Environment During the Galileo Flyby: Global Three-Dimensional MHD Modeling with Adaptive Mesh Refinement

The first results for applying a three-dimensional multimedia ideal MHD model for the mass-loaded flow of Jupiter's corotating magnetospheric plasma past Io are presented. The model is able to consider simultaneously physically realistic conditions for ion mass loading, ion-neutral drag, and intrinsic magnetic field in a full global calculation without imposing artificial dissipation. Io is modeled with an extended neutral atmosphere which loads the corotating plasma torus flow with mass, momentum, and energy. The governing equations are solved using adaptive mesh refinement on an unstructured Cartesian grid using an upwind scheme for AHMED. For the work described in this paper we explored a range of models without an intrinsic magnetic field for Io. We compare our results with particle and field measurements made during the December 7, 1995, flyby of to, as published by the Galileo Orbiter experiment teams. For two extreme cases of lower boundary conditions at Io, our model can quantitatively explain the variation of density along the spacecraft trajectory and can reproduce the general appearance of the variations of magnetic field and ion pressure and temperature. The net fresh ion mass-loading rates are in the range of approximately 300-650 kg/s, and equivalent charge exchange mass-loading rates are in the range approximately 540-1150 kg/s in the vicinity of Io.

Combi, M. R.↗

Final Report for Center for Tokamak Transient Simulations at USU

Providing plasma fluid codes like NIMROD with continuum drift kinetic (CDK) physics that is quantitatively valid and computationally feasible throughout the spatial domain is difficult. Work at Utah State University (USU), in collaboration with the Center for Tokamak Transient Simulations (CTTS), focused on applying CDK closures in disruption-related calculations. Three examples where kinetic physics is paramount are (1) the electron stress tensor closure in Ohms law for accurately describing neoclassical tearing mode (NTM) evolution, (2) runaway electron (RE) density (nRE) and current (jRE) moments in NIMROD’s extended MHD model for self-consistent evolution of RE populations during disruptions and, (3) energetic ion effects on a myriad of MHD instabilities. While NTM simulations and continuum and PIC approaches to energetic ions in NIMROD have been a major goals of USU’s closure work for several years, the development of self-consistent CDK RE capability in NIMROD was started and extended considerably during the CTTS effort. Some goals of CDK RE in NIMROD are to explore the effects of the 2D relativistic phase space in 4D simulations and compare with NIMROD’s fluid RE model. Four publications and two PhD theses came out of the USU CTTS effort.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Experimental studies of plasma-antenna coupling with the JET Alfven Eigenmode Active Diagnostic

This paper presents a dedicated study of plasma-antenna (PA) coupling with the Alfven Eigenmode Active Diagnostic (AEAD) in JET. Stable AEs and their resonant frequencies f, damping rates gamma < 0, and toroidal mode numbers n are measured for various PA separations and limiter versus X-point magnetic configurations. Two stable AEs are observed to be resonantly excited at distinct low and high frequencies in limiter plasmas. The values of f and n do not vary with PA separation. However, |gamma| increases with PA separation for the low-f, but not high-f, mode, yet this may be due to slightly different edge conditions. The high-f AE is detected throughout the transition from limiter to X-point configuration, though its damping rate increases; the low-f mode, on the other hand, becomes unidentifiable. The linear resistive MHD code CASTOR is used to simulate the frequency scan of an AEAD-like external antenna. For the limiter pulses, the high-f mode is determined to be an n = 0 GAE, while the low-f mode is likely an n = 2 TAE. During the transition from limiter to X-point configuration, CASTOR indicates that n = 1 and 2 EAEs are excited in the edge gap. These results extend previous experimental studies in JET and Alcator C-Mod; validate the computational work performed by Dvornova et al 2020 Phys. Plasmas 27 012507; and provide guidance for the optimization of PA coupling in upcoming JET energetic particle experiments, for which the AEAD will aim to identify the contribution of alpha particles to AE drive during the DT campaign.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Hot relativistic winds and the Crab Nebula

Efforts to formulate a self-consistent model of pulsar magnetospheres which links the particle source near the pulsar to the outflowing relativistic wind and couples the wind to the surrounding nebula are reviewed. The use of a relativistic MHD wind is recommended to account for global photon emission and the invisibility of the method of plasma transport. Consideration of a magnetic monopole relativistic wind due to an axially symmetric aligned rotator is combined with calculations of the initial velocity of the wind to show that the flow velocity in such a model will never exceed Mach 1. Extending the solution to the case of a hot relativistic wind at supersonic speeds is noted to yield results consistent with observations of the Crab Nebula

Fujimura, F. S.↗