Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “plasma models”

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 649 records · Page 36
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Natural hot-ion modes in a rotating plasma

We report that in steady state, the fuel cycle of a fusion plasma requires inward particle fluxes of fuel ions. These particle flows are also accompanied by heating. In the case of classical transport in a rotating cylindrical plasma, this heating can proceed through several distinct channels depending on the physical mechanisms involved. Some channels directly heat the fuel ions themselves, whereas others heat electrons. Which channel dominates depends, in general, on the details of the temperature, density and rotation profiles of the plasma constituents. However, remarkably, under relatively few assumptions concerning these profiles, if the α particles, the by-products of the fusion reaction, can be removed directly by other means, then a hot-ion mode tends to emerge naturally.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Experimental observation and integrated modelling of proton-beryllium fusion in He and D plasmas at JET

Validated integrated modelling of JET ITER-like wall experiments in which fusion performance is driven by reactions between fast ions and intrinsically present metal wall impurities is presented. A steady-state L-mode plasma with dominant proton-beryllium fusion and neutron yields of up to ≈ 6·10 13 s -1 is developed in He and D, via radiofrequency heating of a H minority. The fusion drive is unambiguously confirmed by the neutral particle analyser, fast ion loss detector, and γ-ray diagnostics. Experiments are analysed via an integrated modelling framework, developed to model the two-stage proton beryllium-fusion chain and produce high-fidelity fusion product source terms. The modelling chain comprises TRANSP and JETTO for plasma core modelling, LOCUST for full orbit product tracking and collisional slowing-down, DRESS to resolve two- and three-body fusion kinematics, and MCNP for neutron transport calculations. Modelling shows that the primary 9 Be(p,n) 9 B reaction is the dominant neutron emitter at naturally present concentrations of beryllium in these experiments. The yield contribution of secondary reactions between fusion products and beryllium, 9 Be(d,n) 10 B and 9 Be(α,n) 12 C, is found to be negligible. The proton-deuteron knock-on effect in D plasmas is modelled, which is calculated to contribute ≈ 25% to the total neutron yield. For both He and D discharges the total computed neutron rates match fission chamber (FC) measurements within the combined experimental and computational uncertainty, with an average discrepancy of ≈ ± 20%. Realistic proton-beryllium neutron sources are propagated through JET’s MCNP neutron transport model which shows that 235 U FCs’ response is sensitive to p–Be source changes, with up to ≈ 10% variation compared to a D–D neutron source. We show that the high-energy tail of the fast proton minority can be studied with multi-foil neutron activation. The framework is also applied to the study of interactions between fast protons and boron impurities, of relevance to ITER. We calculate that in JET conditions a significant alpha source with DT-like energies could be generated through 11 B(p, α)2α fusion, and detected via γ-emission in secondary interactions between fast alphas and boron. The work represents an important step towards validating predictive integrated modelling capabilities for non-standard fusion reactions.

JET↗
DOE OSTI2019

A computational model for nanosecond pulse laser-plasma interactions

A multi-physics numerical model for laser-induced optical breakdown and laser-plasma interaction in a non-equilibrium gas is presented, accounting for: production of priming electrons via multi-photon ionization, energy absorption, cascade ionization, induced hydrodynamic response, and shock formation and propagation. The gas is governed by the Navier-Stokes equations, with non-equilibrium effects taken into account by means of a two-temperature model. The space-time dependence of the laser beam is modeled with a flux-tube formulation for the Radiative Transfer Equation. The flow governing equations are discretized in space using a second-order finite volume method. The semi-discrete equations are marched in time using an implicit-explicit (IMEX) dual time-stepping strategy, where diffusion and chemistry are solved implicitly, whereas convection is explicit. Application to a 20 ns long 50 mJ pulse laser-induced breakdown in quiescent O 2 shows the advantages of this temporal discretization during and just after the laser pulse, while a less-expensive symmetric Strang splitting (with implicit chemistry) is sufficient for the post-breakdown gas dynamics after ≃ 0.1 μs. The integrated model is shown to reproduce key features of corresponding experiments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Nonlinear and 3D plasmas. Final report

The project objectives are to characterize the nonlinear properties of edge instabilities in fusion plasmas from the standpoint of MHD and extended MHD models of plasmas and to develop, as much as possible, simplified models that may be useful for prediction and control.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1988

Modeling of the thermal plasma in the outer plasmasphere - A magnetospheric heat source

A case study has been carried out using data from the Dynamics Explorer 1 and 2 spacecraft to study the effect of Coulomb interactions between ring current and suprathermal O(+) and thermal protons on the plasmasphere. Results from a one-dimensional plasmaspheric model suggest that heating due to Coulomb collisions may be sufficient to raise the ion and electron temperatures to observed values. The resultant high temperature produced enhancements in the model O(+) and O(++) densities in agreement with observations.

Chandler, M. O.↗
NASA NTRS1994

Steepened structures in equatorial spread F. 2: Theory

This paper investigates the properties of a one-dimensional fluid model of plasma convection in the equatorial F region ionosphere. The model equations are similar in form to Burgers equation except for additional higher-order spatial derivatives. Like Burgers equation, solution to the model have the form of propagating, shocklike structures. Numerical simulations of the model closely resemble the steepened structures observed by sounding rocket plasma density probes within equatorial spread F. Simulated denstiy power spectra, like the spectra computed from in situ data, seem to possess power law forms with a break at wavelengths of about 100 m. The precise wavenumber of the spectral break is determined by the ambipolar diffusion coefficient. The model predicts that electric field fluctuations perpendicular to the direction of plasma steepening should be proportional to the plasma density fluctuations. Electric field fluctuations parallel to the steepening will be due primarily to the ambipolar field and have a Boltzmann relationship with density (square of the absolute value of delta E) approximately equal to (K(exp 2))(square of the absolute value of (delta n/n)). At wavelengths less than about 300 m, the ambipolar field should be the dominant component of the total field intensity.

Hysell, D. L↗
DOE OSTI2022

Laser-driven ionization mechanisms of aluminum for single particle aerosol mass spectrometry

Single particle aerosol mass spectrometry (SPAMS), an analytical technique for measuring the size and composition of individual micron-scale particles, is capable of analyzing atmospheric pollutants and bioaerosols much more efficiently and with more detail than conventional methods which require the collection of particles onto filters for analysis in the laboratory. Despite SPAMS’ demonstrated capabilities, the primary mechanisms of ionization are not fully understood, which creates challenges in optimizing and interpreting SPAMS signals. In this paper, we present a well-stirred reactor model for the reactions involved with the laser-induced vaporization and ionization of an individual particle. The SPAMS conditions modeled in this paper include a 248 nm laser which is pulsed for 8 ns to vaporize and ionize each particle in vacuum. The ionization of 1 μm, spherical Al particles was studied by approximating them with a 0-dimensional plasma chemistry model. The primary mechanism of absorption of the 248 nm photons was pressure-broadened direct photoexcitation to Al(y 2 D). Atoms in this highly excited state then undergo superelastic collisions with electrons, heating the electrons and populating the lower energy excited states. We found that the primary ionization mechanism is electron impact ionization of various excited state Al atoms, especially Al(y 2 D). Because the gas expands rapidly into vacuum, its temperature decreases rapidly. The rate of three-body recombination (e - + e - + Al + → Al + e - ) increases at low temperature, and most of the electrons and ions produced recombine within several μs of the laser pulse. The importance of the direct photoexcitation indicates that the relative peak heights of different elements in SPAMS mass spectra may be sensitive to the available photoexcitation transitions. We also discuss the effects of laser intensity, particle diameter, and expansion dynamics.

47 OTHER INSTRUMENTATION↗
NASA NTRS1980

A computer model of solar panel-plasma interactions

High power solar arrays for satellite power systems are presently being planned with dimensions of kilometers, and with tens of kilovolts distributed over their surface. Such systems face many plasma interaction problems, such as power leakage to the plasma, particle focusing, and anomalous arcing. These effects cannot be adequately modeled without detailed knowledge of the plasma sheath structure and space charge effects. Laboratory studies of 1 by 10 meter solar array in a simulated low Earth orbit plasma are discussed. The plasma screening process is discussed, program theory is outlined, and a series of calibration models is presented. These models are designed to demonstrate that PANEL is capable of accurate self consistant space charge calculations. Such models include PANEL predictions for the Child-Langmuir diode problem.

Cooke, D. L.↗
DOE OSTI2021

Fully Convolutional Spatio-Temporal Models for Representation Learning in Plasma Science

We have trained a fully convolutional spatio-temporal model for fast and accurate representation learning in the challenging exemplar application area of fusion energy plasma science. The onset of major disruptions is a critically important fusion energy science issue that must be resolved for advanced tokamak plasmas such as the $25B burning plasma international thermonuclear experimental reactor (ITER) experiment. While a variety of statistical methods have been used to address the problem of tokamak disruption prediction and control, recent approaches based on deep learning have proven particularly compelling. In the present paper, we introduce further improvements to the fusion recurrent neural network (FRNN) software suite, which delivered cross-machine disruption predictions with unprecedented accuracy using a large database of experimental signals from two major tokamaks. Up to now, FRNN was based on the long short-term memory (LSTM) variant of recurrent neural networks to leverage the temporal information in the data. Here, we implement and apply the "temporal convolutional neural network (TCN)" architecture to the time-dependent input signals. Furthermore, this allows highly optimized convolution operations to carry the majority of the computational load of training, thus enabling a reduction in training time, and the effective use of high-performance computing resources for hyperparameter tuning. At the same time, the TCN-based architecture achieves better predictive performance when compared with the LSTM architecture for various tasks for a representative fusion database.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Closure models for the feedback of energetic particles on plasma turbulence

Energetic particles interact with the plasma surrounding them, resonating with certain types of plasma waves to stabilize them while destabilizing others, and changing the character of the background turbulence in ways that have not been fully quantified or understood. Interaction with the turbulent background plasma is key to the acceleration of many types of energetic particles including high-energy cosmic rays, solar energetic particles, and pick-up ions. The acceleration of particles is a process that would ideally be described by a kinetic model, a type of model that follows a probability distribution function (PDF) for all particles in 7-dimensional (x, y, z, v x , v y , v z , t) space. Because of the high dimensionality of a kinetic model, simulations that solve kinetic equations use the largest computational resources currently available, and are yet unable to simulate a realistic number of particles, reach the large scales necessary for astrophysical problems, and use high-precision numerical methods. Two available alternatives to kinetic plasma models have been explored for this problem, with limited success. One is a multi-fluid model produced by a cumulant discarding closure, which evolves coupled equations for the velocity, magnetic field, and internal energy for both the background plasma and the fluid of energetic particles. However, simulations that solve multi-fluid magnetohydrodynamic (MHD) equations are able to include the interaction with energetic particles only in crude ways, typically as an add-on pressure term. The second alternative is to use a hybrid method to couple a fluid description of the background plasma to a kinetic model or a Fokker–Planck model for the energetic particles. These methods are hampered by the physical modeling of the coupling. In this work, we develop a new model, which follows the PDF for all particles; this can be viewed as a step toward physical realism above a multi-fluid MHD model, while also being more computationally efficient than a kinetic model. The equations we develop model both the background plasma and the energetic particles self-consistently. Over the last decade, similar PDF methods have been developed to a high level of sophistication to model reactive flows and turbulent combustion for engineering applications. For treatment of the feedback of the energetic particles on a background plasma, a PDF closure approach should evaluate the mean characteristics, including the density, with better statistical quality than will particle-sampling procedures.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1995

Chromospheric Evaporation and Decimetric Radio Emission in Solar Flares

We have discovered decimetric signatures of the chromospheric evaporation process. Evidence for the radio detection of chromospheric evaporation is based on the radio-inferred values of (1) the electron density, (2) the propagation speed, and (3) the timing, which are found to be in good agreement with statistical values inferred from the blueshifted Ca xix soft X-ray line. The physical basis of our model is that free-free absorption of plasma emission is strongly modified by the steep density gradient and the large temperature increase in the upflowing flare plasma. The steplike density increase at the chromospheric evaporation front causes a local discontinuity in the plasma frequency, manifested as almost infinite drift rate in decimetric type III bursts. The large temperature increase of the upflowing plasma considerably reduces the local free-free opacity (due to the T-(exp -3/2) dependence) and thus enhances the brightness of radio bursts emitted at the local plasma frequency near the chromospheric evaporation front, while a high-frequency cutoff is expected in the high-density regions behind the front, which can be used to infer the velocity of the upflowing plasma. From model calculations we find strong evidence that decimetric bursts with a slowly drifting high-frequency cutoff are produced by fundamental plasma emission, contrary to the widespread belief that decimetric bursts are preferentially emitted at the harmonic plasma level. We analyzed 21 flare episodes from 1991-1993 for which broadband (100-3000 MHz) radio dynamic spectra from Phoenix, hard X-ray data from BATSE/CGRO, and soft X-ray data from GOES were available. We detected slowly drifting high-frequency cutoffs between 1.1 and 3.0 GHz, with drift rates of -41 +/- 32 MHz/s, extending over time intervals of 24 +/- 23 s. Developing a density model for type III-emitting flare loops based on the statistically observed drift rate of type III bursts by Alvarez & Haddock, we infer velocities of up to 360 km/s for the upflowing plasma, with an average of v(sub CE) = 236 +/- 130 km /s for episodes with 5-15 s duration. The mean electron density of the upflowing plasma is n(sub e) = 5.2(+/-3.1) x 10(exp 10) /cu cm when it is first detected in radio, at coronal altitudes of h(sub 0) = 9.2 +/- 2.3 Mm.

Aschwanden, Markus J.↗
DOE OSTI2021

Volumetric recombination in EMC3-EIRENE: Implementation and first application to the pre-fusion power operation phase in ITER

An upgrade to the 3D plasma boundary model EMC3-EIRENE for detached divertor plasmas is introduced and evaluated. Stabilization of the iterative approximation of a self-consistent plasma–neutral gas solution is required at low divertor plasma temperatures, and this has been achieved by linearization of the energy loss term associated with the electron–neutral gas interaction. An implicit method for volumetric recombination has been implemented, which improves agreement with SOLPS-ITER results for the pre-fusion power operation phase in ITER. Resonant magnetic perturbations for control of edge localized modes are found to result in a partially detached divertor plasma with a significantly broader reduction of particle loads in the traditional strike zone, which is attributed to volumetric recombination. The additional non-axisymmetric far scrape-off layer strike points, on the other hand, remain attached at relatively high temperature so that volumetric recombination is not effective there.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Plasma breakdown in bubbles passing between two pin electrodes

The ignition of plasmas in liquids has applications from medical instrumentation to manipulation of liquid chemistry. Formation of plasmas directly in a liquid often requires prohibitively large voltages to initiate breakdown. Producing plasma streamers in bubbles submerged in a liquid with higher permittivity can significantly lower the voltage needed to initiate a discharge by reducing the electric field required to produce breakdown. The proximity of the bubble to the electrodes and the shape of the bubbles play critical roles in the manner in which the plasma is produced in, and propagates through, the bubble. In this paper, we discuss results from a three-dimensional direct numerical simulation (DNS) used to investigate the shapes of bubbles formed by injection of air into water. Comparisons are made to results from a companion experiment. A two-dimensional plasma hydrodynamics model was then used to capture the plasma streamer propagation in the bubble using a static bubble geometry generated by the DNS The simulations showed two different modes for streamer formation depending on the bubble shape. In an elliptical bubble, a short electron avalanche triggered a surface ionization wave (SIWs) resulting in plasma propagating along the surface of the bubble. In a circular bubble, an electron avalanche first traveled through the middle of the bubble before two SIWs began to propagate from the point closest to the grounded electrode where a volumetric streamer intersected the surface. In an elliptical bubble approaching a powered electrode in a pin-to-pin configuration, we experimentally observed streamer behavior that qualitatively corresponds with computational results. Optical emission captured over the lifetime of the streamer curve along the path of deformed bubbles, suggesting propagation of the streamer along the liquid/gas boundary interface. Plasma generation supported by the local field enhancement of the deformed bubble surface boundaries is a mechanism that is likely responsible for initiating streamer formation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Ion temperature clamping in Wendelstein 7-X electron cyclotron heated plasmas

The neoclassical transport optimization of the Wendelstein 7-X stellarator has not resulted in the predicted high energy confinement of gas fueled electron-cyclotron-resonance-heated (ECRH) plasmas as modelled in (Turkin et al 2011 Phys. Plasmas 18 022505) due to high levels of turbulent heat transport observed in the experiments. The electron-turbulent-heat transport appears non-stiff and is of the electron temperature gradient (ETG)/ion temperature gradient (ITG) type (Weir et al 2021 Nucl. Fusion 61 056001). As a result, the electron temperature T e can be varied freely from 1 keV–10 keV within the range of P ECRH = 1–7 MW, with electron density n e values from 0.1–1.5 × 10 20 m –3 . By contrast, in combination with the broad electron-to-ion energy-exchange heating profile in ECRH plasmas, ion-turbulent-heat transport leads to clamping of the central ion temperature at T i ~ 1.5 keV ± 0.2 keV. In a dedicated ECRH power scan at a constant density of $\langle n_{e} \rangle$ = 7 × 10 19 m –3 , an apparent 'negative ion temperature profile stiffness' was found in the central plasma for (r/a < 0.5), in which the normalized gradient ∇T i /T i decreases with increasing ion heat flux. The experiment was conducted in helium, which has a higher radiative density limit compared to hydrogen, allowing a broader power scan. This 'negative stiffness' is due to a strong exacerbation of turbulent transport with an increasing ratio of T e /T i in this electron-heated plasma. This finding is consistent with electrostatic microinstabilities, such as ITG-driven turbulence. Theoretical calculations made by both linear and nonlinear gyro-kinetic simulations performed by the GENE code in the W7-X three-dimensional geometry show a strong enhancement of turbulence with an increasing ratio of T e /T i . The exacerbation of turbulence with increasing T e /T i is also found in tokamaks and inherently enhances ion heat transport in electron-heated plasmas. This finding strongly affects the prospects of future high-performance gas-fueled ECRH scenarios in W7-X and imposes a requirement for turbulence-suppression techniques.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1988

A model for the spatial and energy distributions of energetic neutral atoms produced within the Saturn/Titan plasma system

A model for the production of energetic neutral atoms (ENAs) in the Saturn/Titan system, that takes into account existing data on the spatial distribution of the background neutral density, the energy and spatial distributions of the energetic ions, and the energy dependence of the charge-exchange cross sections, is presented. Hydrogen production via H(+), H is presented in detail. The extended Saturn/Titan system is 'imaged' in energetic neutral H at 66 R(s) and 25-deg latitude. From this vantage point, the energetic H flux is 6 + or - 1 H/sq cm per sec between 10 and 100 keV. Comparing the predicted H spectrum with that observed by Kirsch et al. (1981), the production of neutral O is estimated, and it is concluded that O(+), H charge exchange in the inner magnetosphere is sufficient to account for all the ENA O production.

Hsieh, K. C.↗
NASA NTRS1990

Plasmasphere-ionosphere coupling. II - Ion composition measurements at plasmaspheric and ionospheric altitudes and comparison with modeling results

This paper presents Dynamic Explorer data on the plasma coupling between the plasmasphere and ionosphere. DE 1 measurements of ion composition and temperatures at 1.4-3.5 R(F) in the plasmasphere were combined with DE-2 measurements of ion composition and electron and ion temperatures in the upper F region/topside ionosphere, closely spaced in universal and local time for cases in the November 6-11, 1981 period. The observations are compared directly with the field-line interhemispheric plasma (FLIP) model calculations of altitudinal ion density and temperature profiles. It was found that, when the FLIP model permitted fractional trapping of ionospheric photoelectrons and consequent plasmaspheric heating, good agreement with the observations was obtained.

Horwitz, J. L.↗
NASA NTRS1985

Numerical modeling of the hybrid plume plasma rocket

This paper reports on ongoing research in the area of numerical simulations of hybrid plume dynamics, and describes the basic mechanisms for energy and particle transport. The solution approach assumes cylindrical geometry and includes a multiplicity of atomic reactions, as well as the presence of a strong magnetic field. The principal reactions are electron impact ionization, and resonant charge exchange between hot ions and cold neutrals. Radial particle and energy transport is mainly by diffusion. The present description includes a modified Bohm diffusion model for plasma in the core of the plume, and classical neutral particle diffusion in the cooler regions of the flow. Neutrals are allowed to 'free stream' in the low density regions, where the collision mean-free-path becomes comparable, or even larger than the characteristic dimensions of the system. A parameteric scan of the various regimes of interest has been completed and the most significant results are presented.

Chang, F. R.↗
NASA NTRS1989

A numerical model of electrodynamics of plasma within the contaminant gas cloud of the Space Shuttle Orbiter at low earth orbit

A two-dimensional cloud was used to study the plasma dynamics within the outgas cloud associated with the Orbiter. It is shown that the polarization field is not symmetric about the direction of motion of the outgas cloud. It rotates in a way that can be predicted in simple cases by the ratio of the Hall and Pederson currents within the outgas cloud. The polarization field magnitude produced in the model was not large.

Eccles, J. Vincent↗