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 253 records · Page 14

First-principles study of L -shell iron and chromium opacity at stellar interior temperatures

Recently developed free-energy density functional theory (DFT)-based methodology for optical property calculations of warm dense matter has been applied for studying L-shell opacity of iron and chromium at T = 182 eV. We use Mermin–Kohn–Sham density functional theory with a ground-state and a fully temperature-dependent generalized gradient approximation exchange-correlation (XC) functionals. It is demonstrated that the role of XC at such a high-T is negligible due to the total free-energy of interacting system being dominated by the noninteracting free-energy term in agreement with estimations for the homogeneous electron gas. Furthermore, our DFT predictions are compared to the radiative emissivity and opacity of dense plasmas model, to the real-space Greens function method, and to experimental measurements. Good agreement is found between all three theoretical methods, and in the bound–continuum region for Cr when compared to the experiment, while the discrepancy between direct DFT calculations and the experiment for Fe remains essentially the same as for plasma-physics models.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Developing Antimatter Containment Technology: Modeling Charged Particle Oscillations in a Penning-Malmberg Trap

The NASA MSFC Propulsion Research Center (PRC) is conducting a research activity examining the storage of low energy antiprotons. The High Performance Antiproton Trap (HiPAT) is an electromagnetic system (Penning-Malmberg design) consisting of a 4 Tesla superconductor, a high voltage confinement electrode system, and an ultra high vacuum test section; designed with an ultimate goal of maintaining charged particles with a half-life of 18 days. Currently, this system is being experimentally evaluated using normal matter ions which are cheap to produce and relatively easy to handle and provide a good indication of overall trap behavior, with the exception of assessing annihilation losses. Computational particle-in-cell plasma modeling using the XOOPIC code is supplementing the experiments. Differing electrode voltage configurations are employed to contain charged particles, typically using flat, modified flat and harmonic potential wells. Ion cloud oscillation frequencies are obtained experimentally by amplification of signals induced on the electrodes by the particle motions. XOOPIC simulations show that for given electrode voltage configurations, the calculated charged particle oscillation frequencies are close to experimental measurements. As a two-dimensional axisymmetric code, XOOPIC cannot model azimuthal plasma variations, such as those induced by radio-frequency (RF) modulation of the central quadrupole electrode in experiments designed to enhance ion cloud containment. However, XOOPIC can model analytically varying electric potential boundary conditions and particle velocity initial conditions. Application of these conditions produces ion cloud axial and radial oscillation frequency modes of interest in achieving the goal of optimizing HiPAT for reliable containment of antiprotons.

Chakrabarti, S.↗

Modelling of edge plasma dynamics with active wall boundary conditions

A self-consistent 2D model is presented for transport in boundary plasma and plasma-facing material walls. Plasma dynamics in the domain is represented by a 2D collisional plasma fluid model in the edge-plasma code UEDGE, and transport of hydrogen and heat in the wall is represented by a system of reaction–diffusion equations in the 1D wall code FACE. To account for variation of parameters along the wall, in the coupled model multiple instances of the FACE code run in parallel. Here, the coupled model provides a tool for investigating a range of dynamic plasma–material interactions phenomena in 2D. For demonstration of its capability, one application of particular interest is the role of active wall in tokamak strike point sweeping proposed for mitigation of divertor heat loads. In the present study, the coupled calculations are applied to investigation of the impact of heat and hydrogen transport in the material wall on the divertor plasma and target heat load during sweeping of the target strike point for parameters of a high-power tokamak.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

A theoretical model for calculation of molecular stopping power

A modified local plasma model is established. The Gordon-Kim's molecular charged density model is employed to obtain a formula to evaluate the stopping power of many useful molecular systems. The stopping power of H2 and He gas was calculated for incident proton energy ranging from 100 keV to 2.5 MeV. The stopping power of O2, N2, and water vapor was also calculated for incident proton energy ranging from 40 keV. to 2.5 MeV. Good agreement with experimental data was obtained. A discussion of molecular effects leading to department from Bragg's rule is presented. The equipartition rule and the effect of nuclear momentum recoiling in stopping power are also discussed.

Xu, Y. J.↗

Plasma thermal transport with a generalized 8-moment distribution function

Moment equations that model plasma transport require an ansatz distribution function to close the system of equations. The resulting transport is sensitive to the specific closure used, and several options have been proposed in the literature. Two different 8-moment distribution functions can be generalized to form a single-parameter family of distribution functions. The transport coefficients resulting from this generalized distribution function can be expressed in terms of this free parameter. This provides the flexibility of matching the 8-moment model to some validating result at a given magnetization value, such as Braginskii’s transport, or the more recent results of Davies et al. [Physics of Plasma, 28, 012305 (2021)]. Here, this process can be thought of as solving for the 8-moment distribution function that matches the value of a transport coefficient given by a Chapman-Enskog expansion, while retaining the improved physical properties, such as finite propagation speeds and time dependence which belong to the hyperbolic moment models. Since the presented generalized distribution function only has a single free parameter, only a single transport coefficient can be matched at a time. However, this generalization process may be extended to provide multiple free parameters. The focus of this Brief Communication is on the dramatically improved thermal conductivity of the proposed model compared to the two base moment models.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Near field interaction of microwave signals with a bounded plasma plume

The objective was to study the effect of the arcjet thruster plume on the performance of an onboard satellite reflector antenna. A project summary is presented along with sections on plasma and electromagnetic modeling. The plasma modeling section includes the following topics: wave propagation; plasma analysis; plume electron density model; and the proposed experimental program. The section on electromagnetic modeling includes new developments in ray modeling and the validation of three dimensional ray results.

Ling, Hao↗

HiPACE++: A portable, 3D quasi-static particle-in-cell code

Modeling plasma accelerators is a computationally challenging task and the quasi-static particle-in-cell algorithm is a method of choice in a wide range of situations. In this work, we present the first performance-portable, quasi-static, three-dimensional particle-in-cell code HiPACE++. By decomposing all the computation of a 3D domain in successive 2D transverse operations and choosing appropriate memory management, HiPACE++ demonstrates orders-of-magnitude speedups on modern scientific GPUs over CPU-only implementations. The 2D transverse operations are performed on a single GPU, avoiding time-consuming communications. The longitudinal parallelization is done through temporal domain decomposition, enabling near-optimal strong scaling from 1 to 512 GPUs. HiPACE++ is a modular, open-source code enabling efficient modeling of plasma accelerators from laptops to state-of-the-art supercomputers.

97 MATHEMATICS AND COMPUTING↗

Ion Transport Modeling in Non-Collisional RF Sheaths

A reactor scale ICP plasma model does not adequately resolve the submillimeter scale rf sheaths due to severe disparities in length and times scales. These non-collisional sheaths determine the spread in the ion bombardment energy distribution on the wafer. Analytical models based on "damped potential" has been used by various authors to simplify sheath modeling. However, inherent in these models are some assumptions that may be invalid. In this work we show that simplifying assumptions such as a uniform ion flux may be in error by as much as 200% at 13.56MHz. This work involves solving ion transport in non-collisional rf sheaths. A comprehensive study of the dependence of ion flux and energy on uncertainties in the sheath-presheath (i.e. sheath thickness) boundary will also be done.

Meyyappan, M.↗

Three-dimensional ray tracing of the Jovian magnetosphere in the low-frequency range

Three-dimensional ray tracing of the Jovian DAM emission has been performed utilizing the O-4 magnetic field model (Acuna and Ness, 1979) and a realistic plasma model. Minimal assumptions about the emission mechanism have been made that include radiation in the right-hand extraordinary mode, propagating nearly perpendicular to the field line at source points located just above the RX cutoff frequency along Io flux tubes. Ray tracing has been performed in the frequency range from 2-35 MHz from successive Io flux tubes separated by ten degrees of central meridian longitude for a full circumference of northern hemisphere sources. The results show unusual complexity in model arc spectra that is displayed in a constant Io phase format with many similarities to the Voyager PRA data. The results suggest much of the variation in observed DAM spectral features is a result of propagation effects rather than emission process differences.

Menietti, J. D.↗

A Coupled Plasma and Sheath Model for High Density Reactors

We present a coupled plasma and collisionless; sheath model for the simulation of high density plasma processing reactors. Due to inefficiencies in numerical schemes and the resulting computational burden, a coupled multidimensional plasma and sheath simulation has not been possible model for gas mixtures and high density reactors of practical interest. In this work we demonstrate that with a fully implicit algorithm and a refined computational mesh, a self-consistent plasma and sheath simulation is feasible. We discuss the details of the model equations, the importance of ion inertia, and the resulting sheath profiles for argon and chlorine plasmas. We find that at low operating pressures (10-30 mTorr), the charge separation occurs only within a 0.5 mm layer near the surface in a 300 mm inductively coupled plasma etch reactor. A unified model eliminates the use of off-line or loosely coupled sheath models with simplifying assumptions which generally lead to uncertainties in ion flux and sheath electrical properties.

Deepak, Bose↗

CCD soft X-ray observations of the Puppis a supernova remnant

We present the first images and spectra of an astronomical object, other than the Sun, acquired with a charge coupled device (CCD) imaging X-ray spectrometer. During a 230 s sounding rocket observation, we have acquired moderate-resolution spectra and moderate-resolution images of a portion of the Pup A supernova remnant (SNR). Based on these data, we conclude that the X-ray spectrum of Pup A is inconsistent with any single-temperature equilibrium or nonequilibrium plasma model. We find evidence for variations in the emitting plasma on scales as small as 5.0 min and as large as 30.0 min. The spatial structure of the spectral variations in the remnant is found to be inconsistent with the standard Sedov model for the evolution of a SRN into a homogeneous interstellar medium (ISM). We suggest that the remnant is expanding into a region of the ISM having a density of approximately 1 cm(exp -3) with inhomogeneities on the order of 50%. We have found evidence for the presence of a knot of plasma enriched in neon, but require more data to be conclusive.

Berthiaume, G. D.↗

Plasma transport in the equatorial ionosphere during the great magnetic storm of March 1989

We have modeled plasma transport in the low-latitude and equatorial ionosphere during the great magnetic storm of March 1989. Our goal was to provide a consistent explanation for the DMSP (Defense Meteorological Satellite Program) observations of dramatic decreases in ion density and rapid ion drifts in the low latitude ionosphere over South America during the storm. The modeling effort supports the hypothesis that abnormally large upward drifts lifted F region plasma above the satellite's altitude and created the density depletions observed by DMSP. Modeled O(+) densities at the satellite's altitude have a strong qualitative resemblance to DMSP observations. Both the model and the observations indicate a deep density trough with extremely sharp boundaries surrounding the equator. The widths of both the modeled and the observed equatorial troughs increase with time. Vertical ion drifts predicted by the model also have been compared with DMSP measurements. Like the observed vertical drifts, the modeled drifts reversed sign near the trough boundaries. The modeled vertical drifts are of the same order and direction as the vertical component of E x B convection near the equator, but of opposite direction (downward) near the trough boundaries and outside of the trough.

Rasmussen, C. E.↗

Atmospheric pressure plasma functionalization of polystyrene

Atmospheric pressure plasma jets (APPJs) are used to improve the adhesive and hydrophilic properties of commodity hydrocarbon polymers such as polypropylene, polyethylene, and polystyrene (PS). These improvements largely result from adding oxygen functional groups to the surface. PS functionalization is of interest to produce high value biocompatible well-plates and dishes, which require precise control over surface properties. In this paper, we discuss results from a computational investigation of APPJ functionalization of PS surfaces using He/O 2 /H 2 O gas mixtures. A newly developed surface reaction mechanism for functionalization of PS upon exposure to these plasmas is discussed. A global plasma model operated in plug-flow mode was used to predict plasma-produced species fluxes onto the PS surface. A surface site balance model was used to predict oxygen-functionalization of the PS following exposure to the plasma and ambient air. We found that O-occupancy on the surface strongly correlates with the O-atom flux to the PS, with alcohol groups and cross-linked products making the largest contributors to total oxygen fraction. Free radical sites, such as alkoxy and peroxy, are quickly consumed in the post-plasma exposure to air through passivation and cross-linking. O-atom fluences approaching 10 17 cm -2 saturate the O-occupancy on the PS surface, creating functionality that is not particularly sensitive to moderate changes in operating conditions.

36 MATERIALS SCIENCE↗

Multi-Scale Modeling for Plasma-Enhanced Ammonia Decomposition over Carbides and Nitrides

Ammonia is a carbon-free hydrogen carrier, but its decomposition typically requires high temperatures over costly Ru-based catalysts due to the large barrier for N≡N bond formation. We develop a multiscale framework combining density functional theory, zero-dimensional plasma kinetics, and microkinetic modeling to elucidate how non-thermal plasma (NTP) enables low-temperature NH 3 decomposition over Co-based carbides and nitrides, benchmarked against Ru and Co. Under thermal conditions, all catalysts are limited by N≡N bond formation, with Co 3 C(001) most active owing to its negatively charged surface, strong N* binding, and low activation barriers of N≡N bond formation. Plasma-induced vibrational excitation of NH 3 and its reactive radicals promotes a radical-driven •NH 2 –N* coupling pathway that dominates on Co 3 C(001) and Co 3 N(001), shifting the rate-limiting step to NH 3 (v1) dissociation, increasing turnover frequencies by up to 6 orders of magnitude, and reducing the temperature needed to reach a turnover frequency of 5 s –1 from >680 °C (Ru and Co under thermal condition) to 267 °C (Co 3 C) and 415 °C (Co 3 N). These results identify Co-based carbides and nitrides as promising plasma-active catalysts for energy-efficient hydrogen production from ammonia.

ammonia decomposition↗

X-Ray Spectra of VY Scl Stars Are Not Blackbodies

Using ASCA data, we find, contrary to other researchers using ROSAT data, that the X, ray spectra of the VY Scl stars TT Ari and KR Aur are poorly fit by an absorbed blackbody model but are well fit by an absorbed thermal plasma model. The different conclusions about the nature of the X-ray spectrum of KR Aur may be due to differences in the accretion rate, since this star was in a high optical state during the ROSAT observation, but in an intermediate optical state during the ASCA observation. TT Ari, on the other hand, was in a high optical state during both observations, so directly contradicts the hypothesis that the X-ray spectra of VY Sol stars in their high optical states are blackbodies. Instead, based on theoretical expectations and the ASCA, Chandra, and XMM spectra of other nonmagnetic cataclysmic variables, we believe that the X-ray spectra of VY Sol stars in their low and high optical states are due to hot thermal plasma in the boundary layer between the accretion disk and the surface of the white dwarf, and appeal to the acquisition of Chandra and XMM grating spectra to test this prediction.

Mauche, C. W.↗

Multi-harmonic electron cyclotron instabilities

The reported investigation constitutes an extension of studies conducted by Ashour-Abdalla and Kennel (1975, 1976, 1978) with respect to a basic plasma model of Young et al. (1973). The model involves a combination of a cold Maxwellian background plasma, a hot plasma, and a 'loss cone' type of free energy source. Previous results on the first cyclotron harmonic bands are extended to multiharmonics. The significance of the obtained relations is discussed and tentative conclusions are presented. Given that the spatial growth rates of the convective modes are comparable, and that simultaneous nonconvective instability (NCI) is possible, it is concluded that multiharmonic emissions ought to be a common feature of the magnetospheric electrostatic wave observations. Since the volume of parameter space for which the first harmonic is NCI, and the volume for which the convective first harmonic mode has significant spatial growth rates, exceed those for the higher harmonics, first harmonic waves should be the most commonly observed and the higher harmonics should usually be accompanied by the first harmonic.

Ashour-Abdalla, M.↗

Detection of an extended soft X-ray source H 2326-79 in the southern sky

A new extended soft X-ray source H 2326-79 has been detected with the low energy detectors of the HEAO1 A-2 experiment. The source is found to have a size of approximately 2.5 arcsec and is centered at RA(1950)=23 h 26 m, Decl.(1950)= -79 deg 24 min. The observed energy flux, seen mainly in the 0.18-0.4 keV band, corrected for finite source size, is found to be about 2.6 x 10 to the -11th erg/sq cm-sec. It is determined that the source spectrum fits well with either an exponential model (gaunt factor included) or with a Raymond Smith plasma model. The observed characteristics of the source are studied for two physical models: a very old low surface brightness supernova remnant undetected in the radio and optical regions, and an enhanced X-ray emission region surrounding an H I cloud embedded in the hot cavity of Loop I. The observed source is also interpreted as a group of two or more unresolved stellar X-ray sources.

Agrawal, P. C.↗