Search NASA⌕ Search

SEARCH · Search NASA

Results for “Numerical magnetohydrodynamics”

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
DOE OSTI2021

Formation of edge pressure pedestal and reversed magnetic shear due to toroidal rotation in a tokamak equilibrium

Toroidal rotation is well known to play significant roles in the edge transport and L-H transition dynamics of tokamaks. Our recent calculation finds that a sufficiently strong localized toroidal rotation can directly bring out the formation of edge pressure pedestal with reversed magnetic shear that is reminiscent of an H-mode plasma, purely through the effects of toroidal rotation on the tokamak MHD equilibrium itself. In particular, the enhanced edge toroidal rotation enables a substantial peaking of the parallel current profile near edge in higher β regimes, which leads to the flattening or reversal of the local q (safety factor) profile. Here the formation of pressure pedestal along with the reversed magnetic shear region is shown to be the natural outcome of the MHD tokamak equilibrium in a self-consistent response to the presence of a localized toroidal rotation typically observed in H-mode or QH-mode.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Benchmarking NIMROD continuum kinetic formulations through the steady-state poloidal flow

In this work, continuum kinetic formulations are employed as a mechanism to include closure physics in an extended magnetohydrodynamics model. Two continuum kinetic approaches have been implemented in the plasma fluid code NIMROD [Sovinec et al., “Nonlinear magnetohydrodynamics with high-order finite elements,” J. Comput. Phys. 195, 355 (2004)] including a Chapman–Enskog-like (CEL) formulation and a more conventional df approach. Ion kinetic closure schemes are employed to describe the neoclassical flow properties in axisymmetric toroidal geometry. In particular, predictions for steady-state values of poloidal flow profiles in tokamak geometry are provided using both the df formulation and two different solution techniques for the CEL approach. These results are benchmarked against analytic theory predictions as well as results from the drift kinetic code DK4D. The continuum kinetic formulations employed here show agreement with both the analytic theory and DK4D results, and offer a novel velocity space representation involving higher-order finite elements in pitch angle.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Augmenting machine learning of Grad–Shafranov equilibrium reconstruction with Green's functions

This work presents a method for predicting plasma equilibria in tokamak fusion experiments and reactors. The approach involves representing the plasma current as a linear combination of basis functions using principal component analysis of plasma toroidal current densities (J t ) from the EFIT-AI equilibrium database. Then utilizing EFIT's Green's function tables, basis functions are created for the poloidal flux (ψ) and diagnostics generated from the toroidal current (J t ). Similar to the idea of a physics-informed neural network (NN), this physically enforces consistency between ψ, J t , and the synthetic diagnostics. First, the predictive capability of a least squares technique to minimize the error on the synthetic diagnostics is employed. The results show that the method achieves high accuracy in predicting ψ and moderate accuracy in predicting J t with median R 2 = 0.9993 and R 2 = 0.978, respectively. A comprehensive NN using a network architecture search is also employed to predict the coefficients of the basis functions. The NN demonstrates significantly better performance compared to the least squares method with median R 2 = 0.9997 and 0.9916 for J t and ψ, respectively. The robustness of the method is evaluated by handling missing or incorrect data through the least squares filling of missing data, which shows that the NN prediction remains strong even with a reduced number of diagnostics. Additionally, the method is tested on plasmas outside of the training range showing reasonable results.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

On resistive interchange, double tearing, and resonant and non-resonant infernal modes in spherical tokamak negative central shear discharges

We examine the linear and nonlinear stability of a sequence of reversed shear NSTX equilibria with the same toroidal current density and pressure profile but with different toroidal field strengths. All equilibria have two q = 2 surfaces and have q > 1 everywhere. For this sequence, which all have β P ∼ 0.5, as the minimum value of q decreases below about 1.4, the localized resistive interchange criterion is strongly violated, and the unstable mode with toroidal mode number n=1 changes its dominant poloidal mode number m from a (m,n)=(2,1) double tearing mode to a non-resonant (1,1) infernal mode. This (1,1) mode nonlinearly flattens both the current density and the pressure near the magnetic axis. The higher-n unstable modes are all resistive infernal modes localized around the surface where q has a minimum and the magnetic shear vanishes. They generally saturate nonlinearly at low amplitude but can cause magnetic surface breakup near the axis.

Spitzer resistivity↗
DOE OSTI2024

Impact of various DIII-D diagnostics on the accuracy of neural network surrogates for kinetic EFIT reconstructions

Abstract Kinetic equilibrium reconstructions make use of profile information such as particle density and temperature measurements in addition to magnetics data to compute a self-consistent equilibrium. They are used in a multitude of physics-based modeling. This work develops a multi-layer perceptron (MLP) neural network (NN) model as a surrogate for kinetic Equilibrium Fitting (EFITs) and trains on the 2019 DIII-D discharge campaign database of kinetic equilibrium reconstructions. We investigate the impact of including various diagnostic data and machine actuator controls as input into the NN. When giving various categories of data as input into NN models that have been trained using those same categories of data, the predictions on multiple equilibrium reconstruction solutions (poloidal magnetic flux, global scalars, pressure profile, current profile) are highly accurate. When comparing different models with different diagnostics as input, the magnetics-only model outputs accurate kinetic profiles and the inclusion of additional data does not significantly impact the accuracy. When the NN is tasked with inferring only a single target such as the EFIT pressure profile or EFIT current profile, we see a large increase in the accuracy of the prediction of the kinetic profiles as more data is included. These results indicate that certain MLP NN configurations can be reasonably robust to different burning-plasma-relevant diagnostics depending on the accuracy requirements for equilibrium reconstruction tasks.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Modeling of Macroscopic Dynamics in Three-Dimensional Plasma Configurations (Final Technical Report)

A numerical simulation-based study to investigate the stability properties of inherently three-dimensional plasma magnetic configurations is presented. These configurations, called stellarators and heliotrons, use external helical coils to provide the twisting magnetic field that is required to confine ionized plasma particles. In contrast to the tokamak configuration, which nominally has the property of axisymmetry, electrical current does not need to be induced in stellarator and heliotron plasma. Experiments have shown that these three-dimensional configurations are remarkably robust to macroscopic instability, even when driven to relatively high levels of plasma pressure. Our understanding of how the plasma is able to withstand macroscopic instability will be enhanced by the work proposed here. The effort will generalize existing numerical modeling capability, and it will apply the new capability to model macroscopic dynamics in these three-dimensional configurations. The development work has been applied to the NIMROD simulation code (https://nimrodteam.org) that was originally based on axisymmetric problem domains. It has been generalized to accommodate the three-dimensional shaping of the plasma volume in stellarators and heliotrons. Its use of equilibrium data has also been generalized, and the implementation allows different numerical representation of the magnetic field. Analysis has been performed to inform the selection of numerical methods, and a method to accelerate the iterative solution of the necessary algebraic systems has been developed and implemented. Verification of the new NIMSTELL code on benchmark problems has also been an important part of the effort.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1986

Formation of active region and quiescent prominence magnetic field configurations

To investigate the formation of prominences, researchers studied chromospheric mass injection into an overlying coronal dipole magnetic field using a 2-D ideal magnetohydrodynamic (MHD) numerical model. Researchers propose that active region prominences are formed by chromospheric plasmas injected directly into the overlying coronal magnetic field and that quiescent prominences are formed by plasmas evaporated at the interface between spicules and corona. Hence, for the simulation of an active region prominence magnetic field we inject the mass from one side, but use a symmetric mass injection to form a quiescent prominence field configuration. Researchers try to find optimum conditions for the formation of Kippenhahn-Schuluter(K-S)type field configuration for stable support of the injection plasmas. They find that the formation of K-S type field configuration by mass injection requires a delicate balance between injection velocity, density, and overlying magnetic fields. These results may explain why a prominence does not form on every neutral line.

An, C.-H.↗
NASA NTRS1985

Shock drift acceleration

This is a review of the fundamental physics of the interactions of charged particles treated individually while they interact with fast mode magnetohydrodynamic shocks. Numerical simulation and analytical theory are used to develop predictions of the expected characteristics of this process strong upstream anisotropies directed along the magnetic field and downstream anisotropies tending to be peaked more perpendicular to the field; relatively more of the enhancement of higher-energy particles occurring upstream; sensitive dependence on shock normal to magnetic field angle of the efficiency of energization. Observations which display all of the above characteristics are reviewed. Also discussed is the relationship of shock drift acceleration to the models for stochastic transport of charged particles in the vicinity of shocks. Extensions of this work in both the observational and theoretical approaches are discussed.

Armstrong, Thomas P.↗
NASA NTRS1987

Comparison of MHD simulation for the February 1986 events with interplanetary observations by the spacecraft Sakigake

During the period 3-10 February 1986 a series of major solar flares occurred on the Sun and several intense geomagnetic storms took place on the Earth. To examine the causality between the solar activity and the geomagnetic activity in this period, a magnetohydrodynamic (MHD) numerical simulation was performed using a 2 1/2 -D numerical code. In that period of February 1986, the Japanese spacecraft Sakigake was at 0.84 AU, 57 deg west of the Earth. Besides the in-situ measurements of the interplanetary plasma, Sakigake also provided Doppler scintillation observations. Comparisons between the results of the MHD simulation and the measurements made by the spacecraft Sakigake are presented.

Smith, Z.↗
NASA NTRS1994

Global, quasi-steady dynamics of the distant solar wind 2: Deformation of the heliospheric current sheet

A three-dimensional (3-D) magnetohydrodynamic (MHD) numerical model is used to trace global deformations of the heliospheric current sheet (HCS) caused by large-scale dynamical interactions associated with corotating solar wind flows. Configurations incorporating the tilted-dipole geometry are investigated out past 30 AU for a variety of dipole tilt angles, alpha. Inclusion of the full, 3-D interplanetary dynamics allows north-south displacements and the east-west warping of the HCS by advective corotational effects to be accurately assessed for the first time. It is found that large-scale spatial correlations between velocity and density imposed at the coronal source (i.e., the geometric arrangement whereby regions of high-speed low-density material lying adjacent to areas of slow, dense flow interact obliquely under the influence of solar rotation) result in a distinctive pattern of deformation of the HCS. For an alpha = 30 deg tilted-dipole example, it is shown that typical zonal variations in radial velocity lead to significant folding of the HCS within about 5 AU of the Sun. By 10 AU, additional sharp bends appear near the latitudinal extremes of the HCS surface, where it is overtaken by shock fronts driven by 3-D corotating interaction regions (CIRs). Moreover, the model suggests that inside about 20 AU, major plasma structures are systematically organized, about the HCS, such that the greatest concentrations of material and magnetic field (the centroids of the 3-D CIR structures) are coincident with the folded crests of the HCS (near heliographic latitudes lambda = +/- alpha). Thus, in these circumstances many of the more interesting dynamical features inconveniently lie well away from the heliographic equator. At larger heliocentric distances, where neighboring CIRs begin to interact strongly, the warping of the HCS abates dramatically and the association between folds in the HCS and major field and density concentrations is weakened and ultimately breaks down. On the basis of the model calculations, quasi-steady deformations and related phenomena should figure prominently in the interplanetary medium, except during those periods when the magnetic dipole and solar spin axes are nearly aligned (alpha is approximately less than 10 deg). Even then, they may be of significance if there are any substantial local warps or kinks in the streamer belt.

Pizzo, V. J.↗
NASA NTRS1994

3-D simulation of high-latitude interaction regions: Comparison with Ulysses results

A three-dimensional (3-D) magnetohydrodynamic (MHD) numerical model is used to simulate the global evolution of a steady, tilted-dipole solar wind flow configuration similar to that prevalent in interplanetary space in 1993. Systematic latitudinal changes in the structure of a corotating interaction region (CIR) near 5 AU is shown to agree well with recent Ulysses observations. The abrupt disappearance of forward shocks and continued persistence of reverse shocks poleward of the latitude where Ulysses crossed the southern edge of the coronal streamer belt is explained as a natural consequence of the 3-D flow geometry.

Pizzo, V. J.↗
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↗
NASA NTRS2000

Two-Species, 3D, MHD Simulation of Europa's Interaction with Jupiter's Magnetosphere

The interaction of Europa with the Jovian a magnetosphere has been studied by using a two species in ideal magnetohydrodynamic (MHD) numerical model. This model considers the upstream plasma in the Jovian magnetosphere and the molecular oxygen ions in the ionosphere of Europa, separately. We present results a from simulation studies, which take into account impact ionization, recombination, and the effect of a possible induced dipole magnetic field of Europa. The total mass loading of the magnetospheric flow and the ionization frequency used in the model are consistent with the estimates of Europa's ionosphere and atmosphere. The multi-species MHD equations are solved by using a finite volume, high-order, Godunov-type method on an adoptively refined unstructured grid, which allows a detailed modeling of the region near Europa's surface, while still resolving both the upstream region and the satellite's wake. We have paid special attention to the wake of Europa, in order to be able to make comparisons with the Galileo's E4 flyby observations, as well as other model calculations. The calculated escape flux of a O2+ down the tail was found to be about 5.6 x 10(exp 25) s(sup -1).

Liu, Yifan↗
NASA NTRS2016

A Model for Straight and Helical Solar Jets: II. Parametric Study of the Plasma Beta

Context. Jets are dynamic, impulsive, well-collimated plasma events that develop at many different scales and in different layers of the solar atmosphere. Aims. Jets are believed to be induced by magnetic reconnection, a process central to many astrophysical phenomena. Within the solar atmosphere, jet-like events develop in many different environments, e.g. in the vicinity of active regions as well as in coronal holes, and at various scales, from small photospheric spicules to large coronal jets. In all these events, signatures of helical structure and/or twisting/rotating motions are regularly observed. The present study aims to establish that a single model can generally reproduce the observed properties of these jet-like events. Methods. In this study, using our state-of-the-art numerical solver ARMS, we present a parametric study of a numerical tridimensional magnetohydrodynamic (MHD) model of solar jet-like events. Within the MHD paradigm, we study the impact of varying the atmospheric plasma beta on the generation and properties of solar-like jets. Results. The parametric study validates our model of jets for plasma beta ranging from 10(sup 3) to 1, typical of the different layers and magnetic environments of the solar atmosphere. Our model of jets can robustly explain the generation of helical solar jet-like events at various beta less than or equal to 1. We show that the plasma beta modifies the morphology of the helical jet, explaining the different observed shapes of jets at different scales and in different layers of the solar atmosphere. Conclusions. Our results allow us to understand the energisation, triggering, and driving processes of jet-like events. Our model allows us to make predictions of the impulsiveness and energetics of jets as determined by the surrounding environment, as well as the morphological properties of the resulting jets.

Magnetic Field↗
NASA NTRS2016

Do We Know the Actual Magnetopause Position for Typical Solar Wind Conditions?

We compare predicted magnetopause positions at the subsolar point and four reference points in the terminator plane obtained from several empirical and numerical MHD (magnetohydrodynamics) models. Empirical models using various sets of magnetopause crossings and making different assumptions about the magnetopause shape predict significantly different magnetopause positions (with a scatter greater than 1 Earth radius (R (sub E)) even at the subsolar point. Axisymmetric magnetopause models cannot reproduce the cusp indentations or the changes related to the dipole tilt effect, and most of them predict the magnetopause closer to the Earth than non axisymmetric models for typical solar wind conditions and zero tilt angle. Predictions of two global non axisymmetric models do not match each other, and the models need additional verification. MHD models often predict the magnetopause closer to the Earth than the non axisymmetric empirical models, but the predictions of MHD simulations may need corrections for the ring current effect and decreases of the solar wind pressure that occur in the foreshock. Comparing MHD models in which the ring current magnetic field is taken into account with the empirical Lin et al. model, we find that the differences in the reference point positions predicted by these models are relatively small for B (sub z) equals 0 (note: B (sub z) is when the Earth's magnetic field points north versus Sun's magnetic field pointing south). Therefore, we assume that these predictions indicate the actual magnetopause position, but future investigations are still needed.

Samsonov, A. A.↗
NASA NTRS1981

Magnetohydrodynamic interaction of high-speed streams

Numerical solutions of a magnetohydrodynamic model are carried out to describe the nonlinear interaction of corotating high-speed streams near the solar equatorial plane. Two problems are studied. The first problem is to simulate the evolution of an idealized high-speed stream. Numerical solutions are obtained to represent the variations of flow velocity, magnetic field, plasma density, temperature, and conduction heat flux in the interaction region. They demonstrate that the dynamical interaction and heat conduction process are responsible for the thermal structure of a high-speed stream. The second problem deals with the formation of corotating shock waves near the leading edge of a broad stream resulting from the merging of characteristic curves. Corotating shocks do not necessarily occur in pairs; a reverse shock can be formed without a forward shock nearby.

Whang, Y. C.↗
NASA NTRS1972

Numerical calculation of the lunar wake in a magnetohydrodynamic model.

Following previous interpretation, a model is developed which has the solar wind flowing directly into the lunar surface and the magnetic field passing unimpeded through the moon. The collapse of the plasma void behind the moon is calculated numerically using the equations of magnetohydrodynamics. The model is three-dimensional, allowing complete calculations to be made without special restrictions on the magnetic field direction. A cursory quantitative comparison with data shows agreement to be satisfactory. The calculations indicate that the rate of collapse of the void does not depend strongly on the solar wind-magnetic field angle. The penumbral decrease depends very strongly on both the solar wind-magnetic field angle and the angle of passage through the wake.

Beers, B. L.↗
NASA NTRS1994

Oblique magnetohydrodynamic cosmic-ray-modified shocks: Two-fluid numerical simulations

We present the first results of time-dependent two-fluid cosmic-ray (CR) modified magnetohydrodynamic (MHD) shock simulations. The calculations were carried out with a new numerical code for one-dimensional ideal MHD. By coupling this code with the CR energy transport equation we can simulate the time-dependent evolution of MHD shocks, including the acceleration of the CR and their feedback on the shock structures. We report tests of the combined numerical method including comparisons with analytical steady state results published earlier by Webb, as well as internal consistency checks for more general MHD CR shock structures after they apppear to have converged to dynamical steady states. We also present results from an initial time-dependent simulation which extends the parameter space domain of previous analytical models. These new results support Webb's suggestion that equilibrium oblique shocks are less effective than parallel shocks in the acceleration of CR. However, for realistic models of anisotropic CR diffusion, oblique shocks may achieve dynamical equilibrium on shorter timescales than parallel shocks.

Frank, Adam↗