Search NASA⌕ Search

SEARCH · Search NASA

Results for “MHD modeling”

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 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
NASA NTRS1978

Dynamic MHD modeling of the solar wind disturbances during the August 1972 events

A time-dependent one-dimensional MHD theoretical model is tested by using plasma and magnetic field observations of Pioneer 9 and Pioneer 10 during the August 1972 events on the sun and in the interplanetary medium. These spacecraft were nearly aligned along a common heliocentric radius during these events, considered now to be the most spectacular and best-documented events during solar cycle 20. The observations of Pioneer 9 at 0.78 AU were used as input for the theoretical model. The plasma and magnetic field forcing functions were superimposed upon a preexisting ambient solar wind at this inner boundary, and the response was simulated as far as 8 AU. The simulated output at 2.2 AU is compared directly with the Pioneer 10 observations at 2.2 AU. Qualitative comparison is good, although several limitations of the one-dimensional theory are noted.

Dryer, M.↗
DOE OSTI2020

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↗
NASA NTRS2014

MHD Modeling of Coronal Loops: the Transition Region Throat

Context. The expansion of coronal loops in the transition region may considerably influence the diagnostics of the plasma emission measure. The cross-sectional area of the loops is expected to depend on the temperature and pressure, and might be sensitive to the heating rate. Aims. The approach here is to study the area response to slow changes in the coronal heating rate, and check the current interpretation in terms of steady heating models. Methods. We study the area response with a time-dependent 2D magnetohydrodynamic (MHD) loop model, including the description of the expanding magnetic field, coronal heating and losses by thermal conduction, and radiation from optically thin plasma. We run a simulation for a loop 50 Mm long and quasi-statically heated to about 4 millikelvin. Results. We find that the area can change substantially with the quasi-steady heating rate, e.g., by approx. 40% at 0.5 millikelvin as the loop temperature varies between 1 millikelvin and 4 millikelvin, and, therefore, affects the interpretation of the differential emission measure vs. temperature (DEM(T)) curves.

magnetohydrodynamics (MHD)↗
NASA NTRS1980

Dynamics of coronal transients - Two-dimensional non-plane MHD models

Numerical results are obtained for non-plane MHD responses to a sudden energy release in a stratified model atmosphere. In agreement with observations, it is shown that after the energy release, the magnetic field affected by the energy release relaxes toward the potential configuration while the outer atmospheric fields increase its shear. Additional results suggest a new way of interpreting the energy storage and release in repeated flares.

Nakagawa, Y.↗
DOE OSTI2023

MHD modeling of shattered pellet injection in JET

Abstract Nonlinear 3D MHD simulations of shattered-pellet injection (SPI) in JET show prototypical SPI-driven disruptions using the M3D-C1 and NIMROD extended-MHD codes. Initially, radiation-driven thermal quenches are accelerated by MHD activity as the pellet crosses rational surfaces, leading to a radiation spike, global stochasticization of the magnetic field, and a complete thermal quench. Eventually, current quenches, preceded by a current spike are seen as the Ohmic heating becomes equal to the radiative cooling. The results are qualitatively similar for both a single monolithic pellet, pencil-beam model, and a realistic shatter to represent the SPI plume. A scan in viscosity from 500 to 2000 m 2 s −1 for MHD simulations finds that reducing viscosity increases MHD activity and decreases thermal quench time slightly. A realistic cloud of fragments modeling shows that mixed-D–Ne pellet travels deeper into the plasma core before the thermal quench. At the slow pellet speeds, the pellet is found to be moving slowly enough inward that even the 5% neon in the mixed pellet is enough to effectively radiate the thermal energy available. Radiation toroidal peaking is predicted to be at levels consistent with experimental observations and reduced as the pellet travels deeper into the plasma. These simulations lay the ground work for more-sophisticated validative and predictive modeling of SPI in JET using both M3D-C1 and NIMROD.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS2003

Dynamical Evolution of the Inner Heliosphere Approaching Solar Activity Maximum: Interpreting Ulysses Observations Using a Global MHD Model

In this study we describe a series of MHD simulations covering the time period from 12 January 1999 to 19 September 2001 (Carrington Rotation 1945 to 1980). This interval coincided with: (1) the Sun s approach toward solar maximum; and (2) Ulysses second descent to the southern polar regions, rapid latitude scan, and arrival into the northern polar regions. We focus on the evolution of several key parameters during this time, including the photospheric magnetic field, the computed coronal hole boundaries, the computed velocity profile near the Sun, and the plasma and magnetic field parameters at the location of Ulysses. The model results provide a global context for interpreting the often complex in situ measurements. We also present a heuristic explanation of stream dynamics to describe the morphology of interaction regions at solar maximum and contrast it with the picture that resulted from Ulysses first orbit, which occurred during more quiescent solar conditions. The simulation results described here are available at: http://sun.saic.com.

Riley, Pete↗
DOE OSTI2022

Non-linear MHD modelling of edge localized modes suppression by resonant magnetic perturbations in ITER

Edge localized modes (ELMs) suppression by resonant magnetic perturbations (RMPs) was studied with the non-linear magneto-hydro-dynamic (MHD) code JOREK for the ITER H-mode scenarios at 15 MA, 12.5 MA, 10 MA/5.3 T. In this study, the main aim was to demonstrate that ELMs can be suppressed by RMPs while the divertor 3D footprints of heat and particle fluxes remain within divertor material limits. The unstable peeling–ballooning modes responsible for ELMs without RMPs were modelled first for each scenario using numerically accessible parameters for ITER. Then the stabilization of ELMs by RMPs was modelled with the same parameters. RMP spectra, optimized by the linear MHD MARS-F code, with main toroidal harmonics N = 2, N = 3, N = 4 have been used as boundary conditions of the computational domain of JOREK, including realistic RMP coils, main plasma, scrape off layer (SOL) divertor and realistic first wall. The model includes all relevant plasma flows: toroidal rotation, two fluid diamagnetic effects and neoclassical poloidal friction. With RMPs, the main toroidal harmonic and the non-linearly coupled harmonics remain dominant at the plasma edge, producing saturated modes and a continuous MHD turbulent transport thereby avoiding ELM crashes in all scenarios considered here. The threshold for ELM suppression was found at a maximum RMP coils current of 45 kAt–60 kAt compared to the coils maximum capability of 90 kAt. In the high beta poloidal steady-state 10 MA/5.3 T scenario, a rotating QH-mode without ELMs was observed even without RMPs. In this scenario with RMPs N = 3, N = 4 at 20 kAt maximum current in RMP coils, similar QH-mode behaviour was observed however with dominant edge harmonic corresponding to the main toroidal number of RMPs. The present MHD modelling was limited in time by few tens of ms after RMPs were switched on until the magnetic energy of the modes saturates. As a consequence the thermal energy was still evolving on this time scale, far from the ITER confinement time scale and hence only the form of 3D footprints on the divertor targets can be indicated within this set-up. Also note, that the divertor physics was missing in this model, so realistic values of fluxes are out of reach in this modelling. However the stationary 3D divertor and particle fluxes could be simply extrapolated from these results to the stationary situation considering that a large power fraction should be radiated in the core and SOL and only about 50 MW power is going to the divertor, which is an arbitrary, but reasonable number used here. The 3D footprints with RMPs show the characteristic splitting with the main RMP toroidal symmetry. The maximum radial extension of the footprints typically was ~20 cm in inner divertor and ~40 cm in outer divertor with stationary heat fluxes decreasing further out from the initial strike point from ~5 MW m ₋2 to ~1 MW m ₋2 assuming a total power in the divertor and walls is 50 MW. The heat fluxes remain within the divertor target and baffle areas, however with rather small margin in the outer divertor which could be an issue for the first wall especially in transient regimes when part of the plasma thermal energy is released due to switching on the RMP coils. This fact should be considered when RMPs are applied with a more favorable application before or soon after the L–H transition, although optimization is required to avoid increasing the L–H power threshold with RMPs.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1981

A three dimensional MHD model of the earth's magnetosphere

The results of a global MHD calculation of the steady state solar wind interaction with a dipole magnetic field are presented. The computer code used, being much faster than previous codes, makes it possible to increase the number of grid points in the system by an order of magnitude. The resulting model qualitatively reproduces many of the observed features of the quiet time magnetosphere including the bow shock, magnetopause, and plasma sheet.

Wu, C. C.↗
NASA NTRS1992

MHD modeling of magnetotail instability for anisotropic pressure

Results of three-dimensional MHD simulations of magnetotail dynamics with anisotropic pressure are presented. The pressure tensor is assumed to be gyrotropic, satisfying a modified double adiabatic approximation including Ohmic heating. It is found that these constraints tend to stabilize the tail with respect to the resistive tearing instability. Including an increasing level of additional isotropization results in a gradual transition to the fast instability found in an isotropic model, which is discussed for comparison. Possible isotropization mechanisms yielding magnetotail instability are discussed also. The stabilizing anisotropies are found primarily in the plasma sheet boundary region. This fact indicates a possibly important role of isotropization mechanisms operating in this region in destabilizing the magnetotail and initiating a tearing instability without the necessity of changes in fluctuation levels in the neutral sheet.

Hesse, Michael↗
NASA NTRS1999

Review of Recent Results in Global MHD Modeling: ISTP Project Scientist for Theory and Ground-Based Observations

Global MHD (magnetohydrodynamic) simulations have shown a remarkable ability to describe the global dynamics of geospace. The limitations of the physical approximations underlying MHD would seem to limit the effectiveness of these codes, since kinetic and hybrid effects should manifest themselves by cross-scale coupling from microscales to mesoscales to global scales. However three effects appear to allow the codes to operate much more successfully than one would at first believe. They are:(l) the globally self-consistent nature of the codes with very well defined exterior boundary conditions (the solar wind) which allows the proper intercommunication between magnetospheric regions on MHD scales, (2) the control by global dynamics of the boundary layer locations where micro and meso scale processes operate, and (3) the critical role of numerical diffusion and with a sufficiently high resolution grid, the use of an empirical resistivity term, which if set at a level where the major magnetosphere boundaries properly calibrate against their observed locations, appear to well represent the effects of kinetic and hybrid processes on the global dynamics. The effectiveness of the global MHD codes, which have been developed under the ISTP mission, in describing Wind, Polar and Geotail observations, as well as ground-based observations are described. Particular emphasis is placed upon the Polar imaging data which when combined with ground-based data and global MHD-based synthetic aurora and convection patterns provide a powerful tool in understanding the final link in the solar-terrestrial chain: coupling into the atmosphere and ionosphere.

Curtis, Steven↗
NASA NTRS1994

MHD modeling of magnetotail instability for localized resistivity

We present results of a three-dimensional magnetohydrodynamic (MHD) simulation of magnetotail evolution initiated by a sudden occurrence or increase of spatially localized resistivity as the major expected concequence of some localized microinstability. Because of the absence of a quantitative model, possible variations of resistivity levels with current density, or the reduction thereof, are not incorporated in the present investigation. The emphasis of the study is on an investigation of the changes to the overall evolution brought about by this localization, in particular, on the disruption and diversion of the cross-tail current and the nonlinear evolution of the magnetotail instability. The immediate consequences of the occurrence of the localized resistance and the resulting electric field are a reduction and diversion of the electric current around the region of high resistivity, associated with an increase of B(sub z) ('dipolarization') at the earthward edge and a decrease of B(sub z) at the tailward edge of this region. These effects, however, are localized and do not involve a reduction of the total cross-tail current and hence do not lead to the global development of a 'substorm current wedge,' which includes not only the reduction of the cross-tail current but also the buildup of a global field-aligned current system of 'regional 1' type (toward the Earth on the dawnside and away on the duskside of the tail). Such signatures develop at a later time, as consequences of a three-dimensional tearing instability, which is triggered by the occurrence of the resistivity. These features are found in combination with plasmoid formation and ejection, quite similar to results of earlier simulations with uniform resistivity. Differences are found in the timescale of the evolution, which tends to be shorter for localized resistivity, and in the propagation of the dipolarization effects in the equatorial plane. Whereas for uniform resistivity the temporal increase in B(sub z) tends to propagate tailward, apparently due to a pileup effect in the near-Earth region, an earthward propagation is found for the localized resistivity. This propagation results from an earthward convection of the gradient in B(sub z), which is set up by the reconnection process further tailward.

Hesse, Michael↗
NASA NTRS1998

Multidimensional MHD Model Studies of the Ionospheres of Venus and Mars

Continuing efforts have been made towards an increased understanding of the solar wind interaction and ionospheric processes at Venus and Mars. This work centered on a systematic development of a new generation of three dimensional magnetohydrodynamic (MHD) numerical code, which models the interaction processes of the solar wind with non-magnetic planets, such as Venus and Mars. We have also worked on a number of different, more specific and discrete studies, as various opportunities arose. We have developed new numerical codes for magnetospheric and cometary studies. As a first step in this process we built an axisymmetric model in which the solar wind interacts with a hard, perfectly conducting sphere. Even that model provided, in certain respects, significant improvements over previous ones.

Nagy, Andrew↗
DOE OSTI2024

Non-linear MHD modelling of transients in tokamaks: a review of recent advances with the JOREK code

Transient magneto-hydrodynamic (MHD) events like edge localized modes (ELMs) or disruptions are a concern for magnetic confinement fusion power plants. Research with the MHD code JOREK towards understanding control of such instabilities is reviewed here in a concise way to provide a complete overview, while we refer to the original publications for details. Experimental validation for unmitigated vertical displacement events progressed. The mechanism of vertical force mitigation by impurity injection was identified. Two-way eddy current coupling to CARIDDI was completed. Shattered pellet injection was simulated in JET, KSTAR, ASDEX Upgrade (AUG) and ITER. Benign runaway electron beam termination in JET and ITER was studied. Coupling of kinetic REs to the MHD is ongoing and a virtual RE synchrotron radiation diagnostic was developed. Regarding pedestal physics, regimes devoid of large ELMs in AUG were simulated and predictive JT60-SA simulations are ongoing. For ELM suppression by resonant magnetic perturbations (RMPs), AUG, ITER and EAST simulations were performed. A free boundary RMP model was validated against experiments. Evidence for penetrated magnetic islands at the pedestal top based on AUG experiments and simulations was found. Simulations of the naturally ELM-free quiescent H-mode in AUG and HL-3 show external kink mode formation prevents pedestal build-up towards an ELM within windows of the edge safety factor. With kinetic neutral particles, high field side high density formation in ITER was simulated and with kinetic impurities, tungsten transport in AUG RMP plasmas was studied. To capture turbulent transport, electro-static full-f particle in cell models for ion temperature gradient and trapped electron modes were established and benchmarked. Application to RMP plasmas shows enhanced turbulence in comparison to unperturbed states. Energetic particle interactions with MHD were studied. Flux pumping that prevents the safety factor on axis from dropping below unity was simulated. First non-linear stellarator applications include current relaxation in $l$ = 2 stellarators, while verification for advanced stellarators progresses.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1978

Dynamic MHD modeling of solar wind corotating stream interaction regions observed by Pioneer 10 and 11

The use of the Pioneer 10 and 11 projects to test an MHD one-dimensional time-dependent model of corotating solar wind streams during the period from Sept. 30 to Nov. 25, 1973 is described. During this period, five or six corotating interaction regions streamed past the two spacecraft, and, as a result of multiple-spacecraft radial alignment and temporally varying conditions at the solar wind source, the pattern predicted by the Steinolfson et al. (1975) model could be compared with observations. The results, in general, support the validity of the model, although the neglect of thermal energy exchange leads to incorrect values for the proton temperature. A detailed analysis of a stream is discussed.

Dryer, M.↗
DOE OSTI2022

Nonlinear MHD modeling of n = 1 RMP -induced pedestal transport and mode coupling effects on ELM suppression in KSTAR

We report fully suppressing edge-localized modes (ELMs), e.g., with resonant magnetic perturbations (RMPs), is essential to reach and sustain high-performance steady-state H-mode plasmas because large ELMs can significantly reduce the lifetime of divertor components in future tokamak reactors. RMP-driven ELM suppression in KSTAR has been modeled by coupling the neoclassical transport code PENTRC to the nonlinear 3D MHD code JOREK. We have found that the radial transport from the combined effects of the kink-peeling, tearing response, and neoclassical toroidal viscosity can explain the pedestal degradation observed in experiments. In addition, it has been found that the RMP response can increase the inter-ELM heat flux on the lower outer divertor by redistributing the heat transport between the divertor plates. In addition to the degraded pedestal, ELM suppression is also attributable to the RMP-induced mode interactions. While the linear stability of peeling-ballooning mode (PBMs) improves owing to the degraded pedestal, the PBM and RMP interaction increases the spectral transfer between edge harmonics, preventing catastrophic growth and the crash of unstable modes. Here, it turns out that the magnetic islands near the pedestal top can play a vital role in mediating the mode interactions.

mode coupling↗
NASA NTRS1974

The ratio of specific heats for postshock plasmas of a detached bow shock - An MHD model

The empirical relationship between the standoff distance of a detached bow shock (generated by the flow of a supersonic gas past an impenetrable obstacle), the size of the obstacle, the Mach number of the gas, and the ratio of specific heats has been generalized to include the magnetic field. The value of the ratio of specific heats (gamma-prime) in the postshock plasma has been calculated in terms of the preshock Alfvenic and sonic Mach numbers and orientation of the magnetic field. The empirical relationship is further generalized by taking into consideration the normal momentum and energy flux due to waves and/or turbulence and/or heat flow in association with high Mach number shocks. The computed value of gamma prime is substantially modified in comparison with that given by the MHD or the gas dynamic model. For this generalized model the computed gamma prime can be considered to be a more precise thermodynamic quantity, since the macroscopic parameters of the plasma have been separated out. Application of this empirical relationship to the earth's bow shock has been given.

Chao, J. K.↗