Search NASA⌕ Search

SEARCH · Search NASA

Results for “PLASMA 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,171 records · Page 65
All sourcesNASA NTRSDOE OSTI
DOE OSTI2016

Extended calculations with spectroscopic accuracy: Energy levels and transition properties for the fluorine-like isoelectronic sequence with Z = 24–30

We have performed extensive multiconfiguration Dirac–Hartree–Fock calculations and second-order many-body perturbation calculations for F-like ions with Z = 24–30. Energy levels and transition rates for electric dipole (E1), electric-quadrupole (E2), electric-octupole (E3), magnetic dipole (M1), and magnetic-quadrupole (M2) transitions, as well as radiative lifetimes, are provided for the lowest 200 levels belonging to the 1s{sup 2}2s{sup 2}2p{sup 5}, 1s{sup 2}2s2p{sup 6}, 1s{sup 2}2s{sup 2}2p{sup 4}3l, 1s{sup 2}2s2p{sup 5}3l, 1s{sup 2}2p{sup 6}3l, and 1s{sup 2}2s{sup 2}2p{sup 4}4l configurations of each ion. The results from the two sets of calculations are in excellent agreement. Extensive comparisons are also made with other theoretical results and observed data from the CHIANTI and NIST databases. The present energies and wavelengths are believed to be accurate enough to aid line identifications involving the n = 3 and n = 4 configurations, for which observations are largely missing. The calculated wavelengths and transition data will be useful in the modeling and diagnostics of astrophysical and fusion plasmas.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2017

The Weakly Nonlinear Magnetorotational Instability in a Global, Cylindrical Taylor–Couette Flow

We conduct a global, weakly nonlinear analysis of the magnetorotational instability (MRI) in a Taylor–Couette flow. This is a multiscale, perturbative treatment of the nonideal, axisymmetric MRI near threshold, subject to realistic radial boundary conditions and cylindrical geometry. We analyze both the standard MRI, initialized by a constant vertical background magnetic field, and the helical MRI, with an azimuthal background field component. This is the first weakly nonlinear analysis of the MRI in a global Taylor–Couette geometry, as well as the first weakly nonlinear analysis of the helical MRI. We find that the evolution of the amplitude of the standard MRI is described by a real Ginzburg–Landau equation (GLE), whereas the amplitude of the helical MRI takes the form of a complex GLE. This suggests that the saturated state of the helical MRI may itself be unstable on long spatial and temporal scales.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

Modeling Emission of Heavy Energetic Neutral Atoms from the Heliosphere

Observations of energetic neutral atoms (ENAs) are a fruitful tool for remote diagnosis of the plasma in the heliosphere and its vicinity. So far, instruments detecting ENAs from the heliosphere were configured for observations of hydrogen atoms. Here, we estimate emissions of ENAs of the heavy chemical elements helium, oxygen, nitrogen, and neon. A large portion of the heliospheric ENAs is created in the inner heliosheath from neutralized interstellar pick-up ions (PUIs). We modeled this process and calculated full-sky intensities of ENAs for energies 0.2–130 keV/nuc. We found that the largest fluxes among considered species are expected for helium, smaller for oxygen and nitrogen, and smallest for neon. The obtained intensities are 50–10{sup 6} times smaller than the hydrogen ENA intensities observed by IBEX . The detection of heavy ENAs will be possible if a future ENA detector is equipped with the capability to measure the masses of observed atoms. Because of different reaction cross-sections among the different species, observations of heavy ENAs can allow for a better understanding of global structure of the heliosphere as well as the transport and energization of PUIs in the heliosphere.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

Fragmentation of a Filamentary Cloud Permeated by a Perpendicular Magnetic Field

We examine the linear stability of an isothermal filamentary cloud permeated by a perpendicular magnetic field. Our model cloud is assumed to be supported by gas pressure against self-gravity in the unperturbed state. For simplicity, the density distribution is assumed to be symmetric around the axis. Also for simplicity, the initial magnetic field is assumed to be uniform, and turbulence is not taken into account. The perturbation equation is formulated to be an eigenvalue problem. The growth rate is obtained as a function of the wavenumber for fragmentation along the axis and the magnetic field strength. The growth rate depends critically on the outer boundary. If the displacement vanishes in regions very far from the cloud axis (fixed boundary), cloud fragmentation is suppressed by a moderate magnetic field, which means the plasma beta is below 1.67 on the cloud axis. If the displacement is constant along the magnetic field in regions very far from the cloud, the cloud is unstable even when the magnetic field is infinitely strong. The cloud is deformed by circulation in the plane perpendicular to the magnetic field. The unstable mode is not likely to induce dynamical collapse, since it is excited even when the whole cloud is magnetically subcritical. For both boundary conditions, the magnetic field increases the wavelength of the most unstable mode. We find that the magnetic force suppresses compression perpendicular to the magnetic field especially in regions of low density.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Black Hole to Photosphere: 3D GRMHD Simulations of Collapsars Reveal Wobbling and Hybrid Composition Jets

Long-duration γ-ray bursts (GRBs) accompany the collapse of massive stars and carry information about the central engine. However, no 3D models have been able to follow these jets from their birth via black hole (BH) to the photosphere. We present the first such 3D general-relativity magnetohydrodynamic simulations, which span over six orders of magnitude in space and time. The collapsing stellar envelope forms an accretion disk, which drags inwardly the magnetic flux that accumulates around the BH, becomes dynamically important, and launches bipolar jets. The jets reach the photosphere at ~10 12 cm with an opening angle θ j ~ 6° and a Lorentz factor Γ j ≲ 30, unbinding ≳90% of the star. We find that (i) the disk–jet system spontaneously develops misalignment relative to the BH rotational axis. As a result, the jet wobbles with an angle θ t ~ 12°, which can naturally explain quiescent times in GRB lightcurves. The effective opening angle for detection θ j + θ t suggests that the intrinsic GRB rate is lower by an order of magnitude than standard estimates. This suggests that successful GRBs are rarer than currently thought and emerge in only ~0.1% of supernovae Ib/c, implying that jets are either not launched or choked inside most supernova Ib/c progenitors. (ii) The magnetic energy in the jet decreases due to mixing with the star, resulting in jets with a hybrid composition of magnetic and thermal components at the photosphere, where ~10% of the gas maintains magnetization σ ≳ 0.1. This indicates that both a photospheric component and reconnection may play a role in the prompt emission.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Three-dimensional Dynamics of Strongly Twisted Magnetar Magnetospheres: Kinking Flux Tubes and Global Eruptions

Abstract The origins of the various outbursts of hard X-rays from magnetars (highly magnetized neutron stars) are still unknown. We identify instabilities in relativistic magnetospheres that can explain a range of X-ray flare luminosities. Crustal surface motions can twist the magnetar magnetosphere by shifting the frozen-in footpoints of magnetic field lines in current-carrying flux bundles. Axisymmetric (2D) magnetospheres exhibit strong eruptive dynamics, i.e., catastrophic lateral instabilities triggered by a critical footpoint displacement of ψ crit ≳ π . In contrast, our new three-dimensional (3D) twist models with finite surface extension capture important non-axisymmetric dynamics of twisted force-free flux bundles in dipolar magnetospheres. Besides the well-established global eruption resulting (as in 2D) from lateral instabilities, such 3D structures can develop helical, kink-like dynamics, and dissipate energy locally (confined eruptions). Up to 25% of the induced twist energy is dissipated and available to power X-ray flares in powerful global eruptions, with most of our models showing an energy release in the range of the most common X-ray outbursts, ≲10 43 erg. Such events occur when significant energy builds up while deeply buried in the dipole magnetosphere. Less energetic outbursts likely precede powerful flares, due to intermittent instabilities and confined eruptions of a continuously twisting flux tube. Upon reaching a critical state, global eruptions produce the necessary Poynting-flux-dominated outflows required by models prescribing the fast radio burst production in the magnetar wind—for example, via relativistic magnetic reconnection or shocks.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2019

Cluster, classify, regress: A general method for learning discontinuous functions

This paper presents a method for solving the supervised learning problem in which the output is highly nonlinear and discontinuous. Here, it is proposed to solve this problem in three stages: (i) cluster the pairs of input-output data points, resulting in a label for each point; (ii) classify the data, where the corresponding label is the output; and finally (iii) perform one separate regression for each class, where the training data corresponds to the subset of the original input-output pairs which have that label according to the classifier. It has not yet been proposed to combine these 3 fundamental building blocks of machine learning in this simple and powerful fashion. This can be viewed as a form of deep learning, where any of the intermediate layers can itself be deep. The utility and robustness of the methodology is illustrated on some toy problems, including one example problem arising from simulation of plasma fusion in a tokamak.

96 KNOWLEDGE MANAGEMENT AND PRESERVATION↗
DOE OSTI2022

Inputs and Outputs for paper "The DESC stellarator code suite. Part 1. Quick and accurate equilibria computations"

Contains the DESC and VMEC input and output files used in the paper, as well as the plotting scripts used to create the figures in the paper (Current with the arXiv Mar 31 2022 version, https://arxiv.org/abs/2203.17173v1): 3D equilibrium codes are vital for stellarator design and operation, and high-accuracy equilibria are also necessary for stability studies. This paper details comparisons of two 3D equilibrium codes, VMEC, which uses a steepest-descent algorithm to reach a minimum-energy plasma state, and DESC, which minimizes the MHD force error in real space directly. Accuracy as measured by final plasma energy and satisfaction of MHD force balance, as well as other metrics, will be presented for each code, along with the computation time. It is shown that DESC is able to achieve more accurate solutions, especially near-axis. DESC's global Fourier-Zernike basis also yields the solution everywhere in the plasma volume, not just on discrete flux surfaces. Further, DESC can compute the same accuracy solution as VMEC in an order of magnitude less time.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1975

Ablative overlays for Space Shuttle leading edge ascent heat protection

Ablative overlays were evaluated via a plasma-arc simulation of the ascent pulse on the leading edge of the Space Shuttle Orbiter. Overlay concepts included corkboard, polyisocyanurate foam, low-density Teflon, epoxy, and subliming salts. Their densities ranged from 4.9 to 81 lb per cu ft, and the thicknesses varied from 0.107 to 0.330 in. Swept-leading-edge models were fabricated from 30-lb per cu ft silicone-based ablators. The overlays were bonded to maintain the surface temperature of the base ablator below 500 F during ascent. Foams provided minimum-weight overlays, and subliming salts provided minimum-thickness overlays. Teflon left the most uniform surface after ascent heating.

Strauss, E. L.↗
NASA NTRS1980

Waves in space plasmas - The mirror trapping of hot auroral electrons

A brief review is given of the problem of precipitation of auroral electrons by electrostatic Bernstein waves. Since the magnetospheric loss cone is small, only moderately small intense levels of wave turbulence are required to remove any large anisotropy sources of free energy and to maintain a weakly anisotropic electron distribution on strong diffusion precipitation. The electrostatic electron cyclotron harmonic waves are nonconvectively unstable for weak loss cone anisotropies and over a large range of parameters for both the hot and cold distributions. Since the instability is nonconvective, weak wave growth can be maintained independent of the flux level of the hot electrons, i.e., the instability does not have the stably trapped flux limit imposed by convective amplification. Recent plasma numerical simulations show that the nonlinear evolution of this instability involves both the pitch angle diffusion of the hot electrons and the heating of the cold electrons.

Ashour-Abdalla, M.↗
NASA NTRS1983

Performance of an ablator for Space Shuttle inorbit repair in an arc-plasma airstream

An ablator patch material performed well in an arc plasma environment simulating nominal Earth entry conditions for the Space Shuttle. Ablation tests using vacuum molded cones provided data to optimize the formulation of a two part polymer system for application under space conditions. The blunt cones were made using a Teflon mold and a state of the art caulking gun. Char stability of formulations with various amounts of catalyst and diluent were investigated. The char was found to be unstable in formulations with low amounts of catalyst and high amounts of diluent. The best polymer system determined by these tests was evaluated using a half tile patch in a multiple High Temperature Reusable surface Insulation tile model. It was demonstrated that this ablator could be applied in a space environment using a state of the art caulking gun, would maintain the outer mold line of the thermal protection system during entry, and would keep the bond line temperature at the aluminum tile interface below the design limit.

Stewart, D. A.↗
NASA NTRS1983

Acceleration of hydrogen ions and conic formation along auroral field lines

Electrostatic ion cyclotron turbulence and the formation of ion conics at low altitudes (about 1500 km) along auroral field lines have been investigated analytically and by plasma numerical simulations. Ion cyclotron waves are assumed to be driven unstable by the upgoing cold ionospheric electrons associated with the downward auroral current. When the electron drift speed is comparable to the electron thermal speed, it was found that the large-amplitude (the saturation level is approximately equal to unity) coherent (omega equals the ion gyrofrequency) ion cyclotron waves should exist along auroral field lines at low altitudes extending a few hundred kilometers. Ion conics are associated with ion cyclotron turbulence, and the ion bulk temperature is found to increase by a factor of 10 from the initial ionospheric temperature, while the temperature of the high-energy tail can be as much as 100 times the ionospheric temperature. Theory and simulations are in good agreement.

Okuda, H.↗