Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “MHD simulation”

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 1,045 records · Page 58
All sourcesNASA NTRSDOE OSTI
NASA NTRS1983

The hydrodynamics of magnetic nonequilibrium

Because the dynamical nonequilibrium of closely packed twisted flux tubes is essentially the problem of two-dimensional MHD turbulence, the extensive literature on two-dimensional turbulence represents the hydrodynamics of magnetic nonequilibrium. Many features of the turbulence can conversely be understood as a consequence of the dynamical nonequilibrium. Attention is presently given to the dynamics of the strong fluid jets, issuing from reconnection points, that are exhibited by numerical simulations of two-dimensional turbulence. Simple examples are given which show that the uniformity of jet width observed is a consequence of the form of the pressure of the magnetic fields between which the jets are confined. The general hydrodynamics of a layer of fluid confined by steady pressures is reducible to a quadrature, so that a variety of effects may be considered, including gravitation. The dissipation of the solar corona is discussed.

Parker, E. N.↗
DOE OSTI2021

On Ohm’s law in reduced plasma fluid models

Drift-reduced MHD models are widely used to study magnetised plasma phenomena, in particular for magnetically confined fusion applications, as well as in solar and astrophysical research. This letter discusses the choice of Ohm's law in these models, the resulting dispersion relations for the dynamics parallel to the magnetic field, and the implications for numerical simulations. We find that if electron pressure is included in Ohm's law, then both electromagnetic and finite electron mass effects must also be included in order to obtain physical dispersion relations. A simple modification to the plasma vorticity is also found which improves handling of low density regions, of particular relevance to the simulation of the boundary region of magnetised plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Modeling of carbon pellets disruption mitigation in an NSTX-U plasma

Single carbon pellet disruption mitigation simulations using M3D-C1 were conducted for an NSTX-U-like plasma to support the electromagnetic pellet injection concept (EPI). A carbon ablation model has been implemented in M3D-C1 and tested with available data. 2D simulations were conducted in order to estimate the amount of carbon needed to quench the plasma, finding that the content in a 1 mm radius vitreous carbon pellet (~3.2x10E20 atoms) would be enough if it is entirely ablated. 3D simulations were performed, scanning over pellet velocity and parallel thermal conductivity, as well as different injection directions and pellet concepts (solid pellets and shell pellets). The sensitivity of the thermal quench and other related quantities to these parameters has been evaluated. A 1 mm radius solid pellet only partially ablates at velocities of 300 m/s or higher, thus being unable to fully quench the plasma. To further enhance the ablation, approximations to an array of pellets and the shell pellet concept were also explored. 3D field line stochastization plays an important role in both quenching the center of the plasma and in heat flux losses, thus lowering the amount of carbon needed to mitigate the plasma when compared to the 2D case. This study constitutes an important step forward in `predict-first' simulations for disruption mitigation in NSTX-U and other devices, such as ITER.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1991

Magnetosphere-ionosphere coupling during plasmoid evolution - First results

The influence of magnetosphere-ionosphere coupling on the dynamic evolution of the magnetotail is investigated by a three-dimensional resistive MHD code that includes the effects of the closure of field-aligned currents in a simple resistive model ionosphere. Particular emphasis is on the role of this coupling during substorm evolution and the modification of the latter by the convection driven by the ionospheric electric fields. For comparison, results are presented from a simulation which uses an infinitely conducting ionosphere but is otherwise identical. Comparison of the two simulations shows that the major impact of magnetosphere-ionospheric communication is an acceleration of magnetotail evolution. Otherwise, phenomena in the two models are qualitatively similar. It is concluded that ionospheric effects do not significantly affect substorm associated magnetotail dynamics.

Hesse, Michael↗
NASA NTRS2003

Large-Eddy Simulation of Conductive Flows at Low Magnetic Reynolds Number

In this paper we study the LES method with dynamic procedure in the context of conductive flows subject to an applied external magnetic field at low magnetic Reynolds number R(sub m). These kind of flows are encountered in many industrial applications. For example, in the steel industry, applied magnetic fields can be used to damp turbulence in the casting process. In nuclear fusion devices (Tokamaks), liquid-lithium flows are used as coolant blankets and interact with the surrounding magnetic field that drives and confines the fusion plasma. Also, in experimental facilities investigating the dynamo effect, the flow consists of liquid-sodium for which the Prandtl number and, as a consequence, the magnetic Reynolds number is low. Our attention is focused here on the case of homogeneous (initially isotropic) decaying turbulence. The numerical simulations performed mimic the thought experiment described in Moffatt in which an initially homogeneous isotropic conductive flow is suddenly subjected to an applied magnetic field and freely decays without any forcing. Note that this flow was first studied numerically by Schumann. It is well known that in that case, extra damping of turbulence occurs due to the Joule effect and that the flow tends to become progressively independent of the coordinate along the direction of the magnetic field. Our comparison of filtered direct numerical simulation (DNS) predictions and LES predictions show that the dynamic Smagorinsky model enables one to capture successfully the flow with LES, and that it automatically incorporates the effect of the magnetic field on the turbulence. Our paper is organized as follows. In the next section we summarize the LES approach in the case of MHD turbulence at low R(sub m) and recall the definition of the dynamic Smagorinsky model. In Sec. 3 we describe the parameters of the numerical experiments performed and the code used. Section 4 is devoted to the comparison of filtered DNS results and LES results. Conclusions are presented in Sec. 5.

Knaepen, B.↗
NASA NTRS1999

A Two-Fluid, MHD Coronal Model

We describe first results from a numerical two-fluid MHD model of the global structure of the solar Corona. The model is two-fluid in the sense that it accounts for the collisional energy exchange between protons and electrons. As in our single-fluid model, volumetric heat and Momentum sources are required to produce high speed wind from Corona] holes, low speed wind above streamers, and mass fluxes similar to the empirical solar wind. By specifying different proton and electron heating functions we obtain a high proton temperature in the coronal hole and a relatively low proton temperature above the streamer (in comparison with the electron temperature). This is consistent with inferences from SOHO/UltraViolet Coronagraph Spectrometer instrument (UVCS), and with the Ulysses/Solar Wind Observations Over the Poles of the Sun instrument (SWOOPS) proton and electron temperature measurements which we show from the fast latitude scan. The density in the coronal hole between 2 and 5 solar radii (2 and 5 R(sub S)) is similar to the density reported from SPARTAN 201.-01 measurements by Fisher and Guhathakurta [19941. The proton mass flux scaled to 1 AU is 2.4 x 10(exp 8)/sq cm s, which is consistent with Ulysses observations. Inside the closed field region, the density is sufficiently high so that the simulation gives equal proton and electron temperatures due to the high collision rate. In open field regions (in the coronal hole and above the streamer) the proton and electron temperatures differ by varying amounts. In the streamer the temperature and density are similar to those reported empirically by Li et al. [1998], and the plasma beta is larger than unity everywhere above approx. 1.5 R(sub S), as it is in all other MHD coronal streamer models [e.g., Steinolfson et al., 1982; also G. A. Gary and D. Alexander, Constructing the coronal magnetic field, submitted to Solar Physics, 1998].

Suess, S. T.↗
NASA NTRS1987

Magnetohydrodynamic (MHD) analyses of various forms of activity and their propagation through helio spheric space

Theoretical and numerical modeling of solar activity and its effects on the solar atmosphere within the context of magnetohydrodynamics were examined. Specifically, the scientific objectives were concerned with the physical mechanisms for the flare energy build-up and subsequent release. In addition, transport of this energy to the corona and solar wind was also investigated. Well-posed, physically self-consistent, numerical simulation models that are based upon magnetohydrodynamics were sought. A systematic investigation of the basic processes that determine the macroscopic dynamic behavior of solar and heliospheric phenomena was conducted. A total of twenty-three articles were accepted and published in major journals. The major achievements are summarized.

Wu, S. T.↗
DOE OSTI2022

Potassium carbonate decomposition modeling within rotating detonation engines for direct power extraction applications

We report that recently, there has been a significant interest in detonation-based combustion systems, such as rotating detonation engines (RDE), due to potential performance advantages in propulsion and energy applications. This includes improvements in thermodynamic efficiency, inherently high gas velocities, and the possibility of achieving an increase in total pressure (i.e. pressure gain) through the combustion process. Due to approximating a constant volume combustion process, these devices also tend to exhibit extremely high local gas temperatures relative to a comparable constant pressure combustor. These advantages overlap well with desirable performance characteristics of direct power extraction (DPE) technologies, such as a magnetohydrodynamic (MHD) generator. Typically, in DPE systems hot combustion products are seeded with an easily ionizable material such as potassium carbonate (K 2 CO 3 ) in order to boost the electrical conductivity. However, due to the short gas residence times within an RDE, it was unclear whether forming an electrically conductive combustion plasma would be feasible for integration with a downstream MHD generator. A model is presented which describes the heating, decomposition, and ionization of solid K 2 CO 3 particles and aqueous solutions of K 2 CO 3 in water, for a given initial particle diameter. This model was combined with available computational fluid dynamics (CFD) data for an oxygen-methane RDE in a one-way coupled Eulerian-Lagrangian framework to predict particle trajectories and the corresponding heating, decomposition, and ionization histories. Electrical conductivities were computed using a previously developed model, and a method was proposed to determine an equivalent average electrical conductivity. Results show that particle sizes below ~30 µm are able to fully decompose before reaching the exit of the RDE. While the one-way coupled nature of the simulations precluded rigorous evaluation of the effects of seed material on detonability, a substantial temperature reduction is expected at the detonation wave front due to heating and decomposition. A preliminary comparison is presented between the RDE and an equivalent constant pressure adiabatic combustor, showing a potential performance advantage for the RDE.

42 ENGINEERING↗
DOE OSTI2022

Hot-tail electrons’ impact on assimilation and injection penetration of D2 shattered pellet injections

The assimilation and the penetration of injected materials play significant roles in the mitigation efficiency of shattered pellet injection (SPI) as a disruption mitigation system (DMS), hence an accurate description of the fragment ablation rate is desirable. Current mainstream 3D MHD codes modelling SPIs mostly assume instantaneous thermalization between the previously hot ambient electrons and the newly released cold electrons. This could result in underestimation of the ablation rate if the hot electron thermalization time is comparable or even longer than the fragment flying time across flux tubes, putting the validity of these 3D SPI simulations into question. To resolve this doubt, we, in this study, investigate the thermalization dynamics and the overall hot-electron impact. The finite-time collisional thermalization of hot-tail electrons in a rapidly cooling plasma, as well as the so-called ‘self-limiting’ effect are considered. The former effect tends to deplete the colder population within a hot-tail species, while the latter is found to preferentially deplete the higher energy population. The combined result is found to cause an almost self-similar decay of the hot electron distribution function, while its shape does not deviate much from that of Maxwellian distribution and the mean energy does not change much during the thermalization process. Based on this observation, axisymmetric JOREK D2 SPI simulations were carried out with additional hot-tail contribution to evaluate their overall impact onto the injection assimilation and penetration. It is found that the hot-tail effect indeed causes enhanced assimilation and shallower penetration, although the overall effect depends on the exact injection configuration, with the slow injection showing negligible hot-tail effect while the fast single non-shattered pellet case shows drastic hot-tail ablation enhancement. For ITER-like SPI parameters, there is no significant deviation in the total assimilation, but some deviation in the injection penetration is observed for the fast injection velocity cases.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1991

Thermoelectric effects in decaying homogeneous magneto-gas turbulence

In the formulation of compressible MHD (i.e., magneto-gas dynamics), a 'generalized Ohm's law is required. In particular, an electron pressure term and a 'Hall effect' term may appear as non-negligible additions to the Ohm's law that is conventionally used for incompressible MHD. In 'high-beta' (i.e., relatively low magnetic energy) situations, the Hall term may be neglected (at least initially) but, as it turns out, the electron pressure term cannot be neglected. Here, three-dimensional, high-beta, homogeneous, decaying, magneto-gas turbulence is examined with regard to this additional term. Through numerical simulation, it is found that 'thermoelectric effects' are produced that significantly alter the evolution of the magnetic field and electric current strengths.

Shebalin, John V.↗
DOE OSTI2021

Photo Doppler Velocimetry and Gamma/Neutron Yield Measurements of Godiva-IV Critical Assembly

The purpose of this experiment was to collect data for a dynamic benchmark of the Godiva-IV burst assembly. Data were taken using two different measurement systems. The gamma & neutron flux was measured using a calibrated, high-dynamic-range scintillator detector, the MHD-240. The measurement will be used to determine the total fission yield and alpha throughout the burst and shutdown sequence of Godiva-IV. The Photon Doppler Velocimetry (PDV) measurement detected the surface motion of the burst assembly as a function of time. The combined data will be used in a multi-physics simulation of the Godiva-IV burst assembly.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2023

NIMROD Development AND Applications for Advanced Simulations of Tokamak Plasmas

Magnetic fusion experiments have the largest gradients in the world that are in steady state. These gradients lead to instabilities, and thus, most of plasma theory over the past four decades have been devoted to the development of instability theory and their use in interpreting plasma results. Instabilities in tokamaks lead to 3 different phenomenological outcomes: coherent saturation, turbulent saturation, and sudden relaxation (which may include disruptions). For long-wavelength instabilities, the NIMROD extended MHD code has emerged as an important tool for understanding tokamak instabilities. Because of the long history of NIMROD, a new version of it has been started at Tech-X to be able to address multi-species capabilities as well as exploit modern GPU systems. The new version builds on the previous version and improves the workflow by enabling new equilibria from experiments to be resolved. The fundamental formulation of the equations underlying NIMROD has also been explored in multiples ways. First, understanding instabilities in tokamaks through numerical methods is aided by being able to understand the rich history of analytic studies. A paper was published which aid in understanding the literature by simplifying some of the analytic machinery inherent in these studies. Extended MHD, like gyrokinetics, are quasineutral models. From a theoretical point of view, plasma quasineutrality approximation is best understood as the same as the magneto-quasistatic approximation of the Maxwell equations. Finally, the standard model of tokamak theory is that of instability theory. A simple dynamical systems model has been developed to better illustrate the strengths and weaknesses of this model.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Predicting operational windows of ELMs suppression by resonant magnetic perturbations in the DIII-D and KSTAR tokamaks

A newly developed plasma response model, combining the nonlinear two-fluid MHD code TM1 and toroidal MHD code GPEC run in ideal mode, quantitatively predicts the narrow isolated q 95 windows (Δq 95 ~ 0.1) of ELM suppression by n = 1, 2 and 3 resonant magnetic perturbations (RMPs) in both DIII-D and KSTAR tokamaks across a wide range of plasma parameters. The key physics that unites both experimental observations and our simulations is the close alignment of essential resonant q-surfaces and the location of the top of the pedestal prior to an ELM. This alignment permits an applied RMP to produce field penetration due to the lower E × B rotation at the pedestal top rather than being screened. The model successfully predicts that narrow magnetic islands form when resonant field penetration occurs at the top of pedestal, and these islands are easily screened when q 95 moves off resonance, leading to very narrow windows of ELM suppression (typically Δq 95 ~ 0.1). Furthermore, the observed reduction in the pedestal height is also well captured by the calculated classical collisional transport across the island. We recover observed q 95 , β N and plasma shape dependence of ELM suppression due to the effect of magnetic islands on pedestal transport and Peeling-Ballooning- Mode (PBM) stability. Importantly, experiments do occasionally observe wide windows of ELM suppression (Δq 95 > 0.5). Our model reveals that at low pedestal-top density multiple islands open, leading to wide operational windows of ELM suppression consistent with experiment. The model indicates that wide q95 windows of ELM suppression can be achieved at substantially higher pedestal pressure with less confinement degradation in DIII-D by operating at higher toroidal mode number (n = 4) RMPs. As a result, this can have significant implications for the operation of the ITER ELM control coils for maintaining high confinement together with ELM suppression.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Ideal internal kink stability in presence of plasma flow and neoclassical toroidal viscosity due to energetic particles

Abstract The influence of energetic particles (EPs) on the ideal internal kink mode, in rotating tokamak plasmas, is numerically investigated by simultaneously solving MHD-kinetic hybrid equations together with a toroidal momentum balance equation utilizing the MARS-Q code (Liu et al 2013 Phys. Plasmas 20 042503). The neoclassical toroidal viscous (NTV) torque, induced by precessional drift resonances of trapped energetic particles, acts as the momentum sink term to damp the plasma flow. Quasi-linear initial value simulations show local reduction of the flow amplitude and enhancement of the flow shear near the q = 1 rational surface ( q is the safety factor) due to EP induced NTV. Both effects in turn destabilize the internal kink mode. These numerical findings are robust against the initial linear stability of internal kink, the initial plasma flow profile, as well as the equilibrium distribution model for EPs.

Physics↗
NASA NTRS1984

The collisionless macroscopic Kelvin-Helmholtz instabiity. I Transverse electrostatic mode

The electrostatic limit of the Kelvin-Helmholtz instability in a collisionless plasma is investigated for a magnetic field perpendicular to the velocity flow. A two-fluid eigenvalue analysis of the transverse electrostatic Kelvin-Helmholtz instability is performed using the Chew-Goldberger-Low hydromagnetic equations with first-order finite Larmer radius corrections, establishing that the transverse electrostatic mode reduces to the incompressible transverse MHD mode in the limit of zero ion gyroradius. The linear growth rates are determined and a strong stabilization of the instability for (k-perpendicular) (rho sub i) of 0.2 or greater is observed. The nonlinear stage of the instability displays large vortices whose size is many times the initial velocity shear length. At saturation the simulation is dominated by the longest-wavelength mode permitted in the system. The simulation results can be understood in terms of a simple model which assumes that saturation occurs when the angular frequency of the vortex motion reaches the linear growth rate of the instability.

Prichett, P. L.↗
DOE OSTI2020

Modeling of the ECCD injection effect on the Heliotron J and LHD plasma stability

The aim of the study is to analyze the stability of the energetic particle modes (EPM) and Alfven Eigenmodes (AE) in Helitron J and LHD plasma if the electron cyclotron current drive (ECCD) is applied. Additionally, the analysis is performed using the code FAR3d that solves the reduced MHD equations describing the linear evolution of the poloidal flux and the toroidal component of the vorticity in a full 3D system, coupled with equations of density and parallel velocity moments for the energetic particle (EP) species, including the effect of the acoustic modes. The Landau damping and resonant destabilization effects are added via the closure relation. The simulation results show that the n = 1 EPM and n = 2 global AE (GAE) in Heliotron J plasma can be stabilized if the magnetic shear is enhanced at the plasma periphery by an increase (co-ECCD injection) or decrease (ctr-ECCD injection) of the rotational transform at the magnetic axis ($\rlap{-} \iota_{0}$). In the ctr-ECCD simulations, the EPM/AE growth rate decreases only below a given $\rlap{-} \iota_{0}$, similar to the ECCD intensity threshold observed in the experiments. In addition, ctr-ECCD simulations show an enhancement of the continuum damping. The simulations of the LHD discharges with ctr-ECCD injection indicate the stabilization of the n = 1 EPM, n = 2 toroidal AE (TAE) and n = 3 TAE, caused by an enhancement of the continuum damping in the inner plasma leading to a higher EP β threshold with respect to the co- and no-ECCD simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1977

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.↗
DOE OSTI2023

On the relationship between the multi-region relaxed variational principle and resistive inner-layer theory

We show that the variational energy principle of the multi-region relaxed magnetohydrodynamic (MRxMHD) model can be used to predict finite-pressure linear tearing instabilities. In this model, the plasma volume is sliced into sub-volumes separated by 'ideal interfaces', and in each volume the magnetic field relaxes to a Taylor state, where the pressure gradient $\nabla p = 0$. The MRxMHD model is implemented in the Stepped-Pressure Equilibrium Code (SPEC) so that the equilibrium solution in each region is computed while preserving the force balance across the interfaces. As SPEC computes the Hessian matrix (a discretized stability matrix), the stability of an MRxMHD equilibrium can also be computed with SPEC. In this article, using SPEC, we investigate the effect of local pressure gradients and the $\nabla p = 0$ in the vicinity of the resonant surface of a tearing mode. For low-beta plasma, we have been able to illustrate a relationship between the resistive singular-layer theory (Coppi et al 1966 Nucl. Fusion6 101; Glasser et al 1975 Phys. Fluids18 875–88) and the MRxMHD model. Within the singular layer, the volume-averaged magnetic helicity and the flux-averaged toroidal flux are shown to be the invariants for the linear tearing modes in SPEC simulations. Our technique to compute MRxMHD stability is first tested numerically in a cylindrical tokamak and its application in toroidal geometry is demonstrated. We demonstrate an agreement between the stability boundary obtained with SPEC simulation and the resistive inner-layer theories.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗