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,603 records · Page 89
All sourcesNASA NTRSDOE OSTI
NASA NTRS1993

Plasma wave signatures in the magnetotail reconnection region - MHD simulation and ray tracing

An MHD simulation was performed to obtain a self-consistent model of magnetic field and plasma density near the X point reconnection region. The MHD model was used to perform extensive ray tracing calculations in order to clarify the propagation characteristics of the plasma waves near the X point reconnection region. The dynamic wave spectra possibly observed by the Geotail spacecraft during a typical cross-tail trajectory are reconstructed. By comparing the extensive ray tracing calculations with the plasma wave data from Geotail, it is possible to perform a kind of 'remote sensing' to identify the location and structure of potential X point reconnection regions.

Omura, Yoshiharu↗
DOE OSTI2021

Hole formation effect on surface morphological response of plasma-facing tungsten

We report simulation results on the effect of helium (He) bubble bursting-mediated surface hole formation on the surface morphological response of tungsten plasma-facing components (PFCs) in nuclear fusion devices. Our analysis is based on an atomistically informed, continuum-scale model, which is capable of accessing the spatiotemporal scales relevant to the fuzz nanostructure formation process on the surface of PFC tungsten. Our simulations account, in an empirical fashion, for two types of subsurface bubble dynamical phenomena in the nanobubble region of PFC tungsten during He plasma irradiation, involving bubble bursting and surface crater formation. We demonstrate that the hole formation effect on the PFC tungsten surface accelerates the growth rate of nanotendrils and the onset of fuzz formation. As a result, the predicted incubation time for surface nanotendril growth is shortened, in agreement with experimental data of incubation fluence at comparable plasma exposure conditions. We also explore systematically the dependence of the PFC surface morphological response on the areal density of holes introduced at regular time intervals onto the He-implanted tungsten surface, a parameter in our analysis that serves as a proxy for the rate of He bubble bursting. More importantly, our simulations capture fine surface features in the PFC tungsten surface morphology and predict that the average spacing between nanotendrils is on the order of 100 nm, consistent with the experimental findings.

Chen, Chao-Shou (ORCID:0000000307910042)↗
DOE OSTI2019

Real-time capable modeling of neutral beam injection on NSTX-U using neural networks

A new model of heating, current drive, torque and other effects of neutral beam injection on NSTX-U that uses neural networks has been developed. The model has been trained and tested on the results of the Monte Carlo code NUBEAM for the database of experimental discharges taken during the first operational campaign of NSTX-U. By projecting flux surface quantities onto empirically derived basis functions, the model is able to efficiently and accurately reproduce the behavior of both scalars, like the total neutron rate and shine through, and profiles, like beam current drive and heating. The model has been tested on the NSTX-U real-time computer, demonstrating a rapid execution time orders of magnitude faster than the Monte Carlo code that is well suited for the iterative calculations needed to interpret experimental results, optimization during scenario development activities, and real-time plasma control applications. Simulation results of a proposed design for a nonlinear observer that embeds the neural network calculations to estimate the poloidal flux profile evolution, as well as effective charge and fast ion diffusivity, are presented.

Neural networks↗
DOE OSTI2025

Comparison of laser-produced plasma spatio-temporal electron density evolution measured using interferometry with simulation results

Due to the difficulties associated with experimental measurements of laser-produced plasma (LPP) properties during the earliest stages of plasma evolution, radiation hydrodynamic codes are often used. However, although these codes have been extensively validated in the higher intensity regimes, validation at low to moderate intensities has been limited. In this study, the spatio-temporal electron density evolution in an LPP generated at moderate laser intensities and at various laser wavelengths was validated against the FLASH code for times up to 20 ns. The LPP was generated by focusing the fundamental and various harmonics radiation (1064, 532, and 266 nm) from a 6 ns full width half maximum Nd:YAG laser, at a laser intensity of 10 GW cm -2 , onto a copper target. The spatio-temporal density evolution of the expanding plasma was analyzed using Nomarski interferometry. Experimental measurements were found to be consistent with FLASH simulations, and the dependence of electron densities on wavelength was found to be in agreement with analytical models, varying as n e ∝ λ -0.7 . However, slight differences were noted in the widths and shape functions of the experimental and simulated electron density profiles.

2D fast imaging↗
NASA NTRS1983

Transverse acceleration of ions on auroral field lines

This paper examines transverse ion heating on auroral field lines associated with current-driven electrostatic ion cyclotron waves theoretically and by numerical simulations. The auroral plasma is assumed to consist of drifting electrons and stationary hydrogen and oxygen ions. Depending on the ratio of the electron drift speed to the thermal speed and the ratio of hydrogen to oxygen concentrations, preferential heating of either hydrogen or oxygen ions can take place. It is found that unless the oxygen ions are a minority species, oxygen transverse heating generally exceeds that of hydrogen ions. Theory and numerical simulations are in good agreement.

Ashour-Abdalla, M.↗
DOE OSTI2016

Comparisons of time explicit hybrid kinetic-fluid code Architect for Plasma Wakefield Acceleration with a full PIC code

Architect, a time explicit hybrid code designed to perform quick simulations for electron driven plasma wakefield acceleration, is described. In order to obtain beam quality acceptable for applications, control of the beam-plasma-dynamics is necessary. Particle in Cell (PIC) codes represent the state-of-the-art technique to investigate the underlying physics and possible experimental scenarios; however PIC codes demand the necessity of heavy computational resources. Architect code substantially reduces the need for computational resources by using a hybrid approach: relativistic electron bunches are treated kinetically as in a PIC code and the background plasma as a fluid. Cylindrical symmetry is assumed for the solution of the electromagnetic fields and fluid equations. In this paper both the underlying algorithms as well as a comparison with a fully three dimensional particle in cell code are reported. The comparison highlights the good agreement between the two models up to the weakly non-linear regimes. In highly non-linear regimes the two models only disagree in a localized region, where the plasma electrons expelled by the bunch close up at the end of the first plasma oscillation.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS2011

Global Variation of Meteor Trail Plasma Turbulence

We present the first global simulations on the occurrence of meteor trail plasma irregularities. These results seek to answer the following questions: when a meteoroid disintegrates in the atmosphere will the resulting trail become plasma turbulent, what are the factors influencing the development of turbulence, and how do they vary on a global scale. Understanding meteor trail plasma turbulence is important because turbulent meteor trails are visible as non-specular trails to coherent radars, and turbulence influences the evolution of specular radar meteor trails, particularly regarding the inference of mesospheric temperatures from trail diffusion rates, and their usage for meteor burst communication. We provide evidence of the significant effect that neutral atmospheric winds and density, and ionospheric plasma density have on the variability of meteor trail evolution and the observation of nonspecular meteor trails, and demonstrate that trails are far less likely to become and remain turbulent in daylight, explaining several observational trends using non-specular and specular meteor trails.

Dyrud, L. P.↗
NASA NTRS2005

Processing and Synthesis of Pre-Biotic Chemicals in Hypervelocity Impacts

Hypervelocity impacts (HVIs) may have played a significant role in establishing the initial organic inventory for pre-biotic chemistry on the Earth and other planetary bodies. In addition to the delivery of organic compounds intact to planetary surfaces, generally at velocities below approx.20 km/s, HVIs also enable synthesis of new molecules. The cooling post-impact plasma plumes of HVIs in the interstellar medium (ISM), the protosolar nebula (PSN), and the early solar system comprise pervasive conditions for organic synthesis. Such plasma synthesis (PS) can operate over many length scales (from nm-scale dust to planets) and energy scales (from molecular rearrangement to atomization and recondensation). HVI experiments with the flexibility to probe the highest velocities and distinguish synthetic routes are a high priority to understand the relevance of PS to exobiology. We describe here recent studies of PS at small spatial scales and extremely high velocities with pulsed laser ablation (PLA). PLA can simulate the extreme plasma conditions generated in impacts of dust particles at speeds of up to 100 km/s or more. When applied to carbonaceous solids, new and pre-biotically relevant molecular species are formed with high efficiency [1,2].

Brickerhoff, W. B.↗
DOE OSTI2015

Influence of wall plasma on microwave frequency and power in relativistic backward wave oscillator

The RF breakdown of the slow wave structure (SWS), which will lead to the generation of the wall plasma, is an important cause for pulse shortening in relativistic backward wave oscillators. Although many researchers have performed profitable studies about this issue, the influence mechanism of this factor on the microwave generation still remains not-so-clear. This paper simplifies the wall plasma with an “effective” permittivity and researches its influence on the microwave frequency and power. The dispersion relation of the SWS demonstrates that the introduction of the wall plasma will move the dispersion curves upward to some extent, which is confirmed by particle-in-cell (PIC) simulations and experiments. The plasma density and volume mainly affect the dispersion relation at the upper and lower frequency limits of each mode, respectively. Meanwhile, PIC simulations show that even though no direct power absorption exists since the wall plasma is assumed to be static, the introduction of the wall plasma may also lead to the decrease in microwave power by changing the electrodynamic property of the SWS.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2015

Molecular dynamics simulations of Si etching in Cl- and Br-based plasmas: Cl{sup +} and Br{sup +} ion incidence in the presence of Cl and Br neutrals

Classical molecular dynamics (MD) simulations have been performed for Cl{sup +} and Br{sup +} ions incident on Si(100) surfaces with Cl and Br neutrals, respectively, to gain a better understanding of the ion-enhanced surface reaction kinetics during Si etching in Cl- and Br-based plasmas. The ions were incident normally on surfaces with translational energies in the range E{sub i} = 20–500 eV, and low-energy neutrals of E{sub n} = 0.01 eV were also incident normally thereon with the neutral-to-ion flux ratio in the range Γ{sub n}{sup 0}/Γ{sub i}{sup 0} = 0–100, where an improved Stillinger--Weber potential form was employed for the interatomic potential concerned. The etch yields and thresholds presently simulated were in agreement with the experimental results previously reported for Si etching in Cl{sub 2} and Br{sub 2} plasmas as well as in Cl{sup +}, Cl{sub 2}{sup +}, and Br{sup +} beams, and the product stoichiometry simulated was consistent with that observed during Ar{sup +} beam incidence on Si in Cl{sub 2}. Moreover, the surface coverage of halogen atoms, halogenated layer thickness, surface stoichiometry, and depth profile of surface products simulated for Γ{sub n}{sup 0}/Γ{sub i}{sup 0} = 100 were in excellent agreement with the observations depending on E{sub i} reported for Si etching in Cl{sub 2} plasmas. The MD also indicated that the yield, coverage, and surface layer thickness are smaller in Si/Br than in Si/Cl system, while the percentage of higher halogenated species in product and surface stoichiometries is larger in Si/Br. The MD further indicated that in both systems, the translational energy distributions of products and halogen adsorbates desorbed from surfaces are approximated by two Maxwellians of temperature T{sub 1} ≈ 2500 K and T{sub 2} ≈ 7000–40 000 K. These energy distributions are discussed in terms of the desorption or evaporation from hot spots formed through chemically enhanced physical sputtering and physically enhanced chemical sputtering, which have so far been speculated to both occur in the ion-enhanced surface reaction kinetics of plasma etching.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1991

The formation of quasi-parallel shocks

In a collisionless plasma, the coupling between a piston and the plasma must take place through either laminar or turbulent electromagnetic fields. Of the three types of coupling (laminar, Larmor and turbulent), shock formation in the parallel regime is dominated by the latter and in the quasi-parallel regime by a combination of all three, depending on the piston. In the quasi-perpendicular regime, there is usually a good separation between piston and shock. This is not true in the quasi-parallel and parallel regime. Hybrid numerical simulations for hot plasma pistons indicate that when the electrons are hot, a shock forms, but does not cleanly decouple from the piston. For hot ion pistons, no shock forms in the parallel limit: in the quasi-parallel case, a shock forms, but there is severe contamination from hot piston ions. These results suggest that the properties of solar and astrophysical shocks, such as particle acceleration, cannot be readily separated from their driving mechanism.

Cargill, Peter J.↗
NASA NTRS2019

Plasma Science in Planetary Entry

Spacecraft entering a planetary atmosphere dissipate a great deal of energy into the surrounding gas. In the frame of reference of the vehicle, the atmospheric gas suddenly decelerates from hypersonic (Mach ~5-50) to subsonic velocities. The kinetic energy of the gas is rapidly converted to thermal and chemical energy, forming a bow shock behind which a plasma with energies on the order of one electron volt (eV) is produced. The resulting shock layer relaxes from strong thermal non-equilibrium that is translationally hot but internally cold and un-ionized toward a thermochemically equilibrated plasma over a distance of a few centimeters. Composition is dependent upon the planetary atmosphere – Air for Earth, CO2/N2 for Mars and Venus, N2/CH4 for Titan and H2/He/CH4 for Saturn, Neptune and Jupiter. Typical velocities of entry may range from 3-7 km/s (4-25 MJ/kg) for Titan/Mars, 8-14 km/s (30-100 MJ/kg) for Earth/Venus, and 25-40 km/s (300-800 MJ/kg) for outer planets. The equilibrium plasmas produced from these conditions are highly dissociated (up to and above 99%) and ionized (0.1- 15%), with temperatures from 7,000-15,000K and pressures from 0.1-1.0 bar. Understanding the behavior of these plasmas – the way in which they approach equilibrium, how they radiate, and how they interact with materials – is an active area of research necessitated by requirements to predict and test the performance of thermal protection systems (TPS) that enable spacecraft to deliver scientific instruments, and people, to foreign worlds and back to Earth. The endeavor is a multi-physics problem, with key processes highlighted in Fig. 1. This white paper describes the current state of the art in simulating shock layer plasmas both computationally and in ground test facilities. Gaps requiring further research and development are identified.

Cruden, Brett A.↗
DOE OSTI2017

Simulation of perturbation produced by an absorbing spherical body in collisionless plasma

A steady plasma state reached in the course of charging of an absorbing spherical body is found using computational methods. Numerical simulations provide complete information on this process, thereby allowing one to find the spatiotemporal dependences of the physical quantities and observe the kinetic phenomena accompanying the formation of stable electron and ion distributions in phase space. The distribution function of trapped ions is obtained, and their contribution to the screening of the charged sphere is determined. The sphere charge and the charge of the trapped-ion cloud are determined as functions of the unperturbed plasma parameters.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Full-wave modeling of RF waves in fusion plasmas with finite element method: Progress in past decades and its future role

This paper reviews the progress in computing radio frequency (RF) wave fields in fusion plasmas, specifically focusing on simulations utilizing the finite element method (FEM) over the past few decades. Computing RF wave fields in fusion plasmas presents unique challenges due to large simulation domains, complex antenna structures, non-local dielectric properties, and wide ranges of spatial scales. It highlights key developments and outlines future directions, primarily addressing waves in the ion cyclotron (IC) to lower hybrid (LH) frequency range. We begin with briefly revisiting earlier developments before the widespread availability of modern computer-aided engineering (CAE) software based on FEM. This historical perspective illuminates early progress and the physics difficulties that motivated ongoing work within the community. Modern wave simulations for RF antennas based on FEM are characterized by the use of detailed 3D antenna model geometry generated from engineering CAD software and localized wave dielectric model. Significant advancements have also been made in improving physics models to include phenomena such as RF sheath rectification and wave scattering. FEM-based RF simulations have also been applied to compute wave propagation in the core region, where the inclusion of non-local dielectric response is crucial. This is a challenging goal, and several promising approaches have been proposed in this area. Additionally, RF simulation development initiatives based on open-source libraries have gained popularity, demonstrating scalability and flexibility in extending physics models. This paper will discuss the advantages and disadvantages of using such a publicly available FEM library.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1995

Numerical simulation of filling a magnetic flux tube with a cold plasma: Anomalous plasma effects

Large-scale models of plasmaspheric refilling have revealed that during the early stage of the refilling counterstreaming ion beams are a common feature. However, the instability of such ion beams and its effect on refilling remain unexplored. In order to learn the basic effects of ion beam instabilities on refilling, we have performed numerical simulations of the refilling of an artificial magnetic flux tube. (The shape and size of the tube are assumed so that the essential features of the refilling problem are kept in the simulation and at the same time the small scale processes driven by the ion beams are sufficiently resolved.) We have also studied the effect of commonly found equatorially trapped warm and/or hot plasma on the filling of a flux tube with a cold plasma. Three types of simulation runs have been performed.

Singh, Nagendra↗
NASA NTRS1974

Thermal relaxation of a two-dimensional plasma in a d.c. magnetic field. I - Theory. II - Numerical simulation

Theoretical considerations relevant to the rate of thermal relaxation of a two-dimensional plasma in a strong uniform dc magnetic field are developed. The Vahala-Montgomery (1971) kinetic description is completed by providing a cut-off time for the time of interaction of two particles contributing to the collision term. The kinetic equation is shown to predict that thermal relaxation varies as a function of defined dimensionless time.

Hsu, J.-Y.↗
NASA NTRS1982

Structure of field rotating disturbances in warm plasma

A model in which thermal effects are simulated through use of a multibeam plasma distribution function is developed and investigated to see if solutions which take an initially uniform magnetized plasma to a new uniform state with different field orientation are possible. The momentum conservation integrals are found to admit two classes of such solutions, but only one class exhibits appropriate asymptotic behavior. Extensive numerical integrations have failed to demonstrate the existence of the desired solutions.

Wolfson, R.↗
DOE OSTI2015

Beyond Orbital-Motion-Limited theory effects for dust transport in tokamaks

Dust transport in tokamaks is very important for ITER. Can many kilograms of dust really accumulate in the device? Can the dust survive? The conventional dust transport model is based on Orbital-Motion-Limited theory (OML). But OML can break in the limit where the dust grain becomes positively charged due to electron emission processes because it overestimates the dust collected power. An OML + approximation of the emitted electrons trapped/passing boundary is shown to be in good agreement with PIC simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗