Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Neutrals in plasma”

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 73 records · Page 4

A ground-based experimental test program to duplicate and study the spacecraft glow phenomenon

The use of a plasma device, the Advanced Concepts Torus-I, for producing atoms and molecules to study spacecraft glow mechanisms is discussed. A biased metal plate, located in the plasma edge, is used to accelerate and neutralize plasma ions, thus generating a neutral beam with a flux approx. 5 x 10 to the 14th power/sq cm/sec at the end of a drift tube. Our initial experiments are to produce a 10 eV molecular and atomic nitrogen beam directed onto material targets. Photon emission in the spectral range 2000 to 9000 A from excited species formed on the target surface will be investigated.

Langer, W. D.↗

Mean force kinetic theory applied to self-diffusion in supercritical Lennard-Jones fluids

A new kinetic-theory-based calculation of the self-diffusion coefficient for dense supercritical Lennard-Jones fluids is presented. The mean force kinetic theory, which was recently developed for transport in dense plasmas, is applied for the calculation of diffusion in dense neutral fluids. The calculation only requires the pair distribution function, a quantity that is readily calculable from equilibrium statistical mechanics for many systems, including the Lennard-Jones fluid. The self-diffusion coefficients are compared with calculations from molecular dynamics simulations, and good agreement at high density is demonstrated, even in the vicinity of the solid–fluid coexistence line. A comparison of different kinetic models with molecular dynamics simulations demonstrates that the transport coefficients have important contributions due to particle interaction via a potential of mean force and local correlations, which increase the collision rate. The new calculations compare well to those from free-volume theory and overcome a limitation of this theory that prevents its use in systems that interact via long range monotonic potentials. It is expected that this approach will also apply to other systems, including neutral-plasma and neutral-electrolyte mixtures.

79 ASTRONOMY AND ASTROPHYSICS↗

The warped neutral sheet and plasma sheet in the near-earth geomagnetic tail

An analysis of ISEE 2 plasma and magnetic field data indicates that the plasma sheet and neutral sheet in the near-earth magnetotail are warped in such a manner that in summer (winter) the neutral sheet rises above (dips below) the solar magnetospheric equatorial plane near the center of the tail, but dips below (rises above) the equatorial plane along the tail flanks. In the near tail, the neutral sheet crosses the equatorial plane at about 12 earth radii from the aberrated X axis, considerably closer to the center of the tail than has been inferred from data obtained farther downstream. This increase in the warp with decreasing distance from the earth is consistent with theoretical predictions. In the near tail, the warp is sufficiently strong when the dipole tilt angle is large that even in quiet times the upper or lower edge of the plasma sheet can be found close to the solar magnetospheric equatorial plane along the tail flanks. The seasonal dependence of the warp can produce certain dawn-dusk asymmetries in satellite data which are more apparent than real.

Gosling, J. T.↗

Eulerian simulations of electrostatic waves in plasmas with a single sign of charge

An Eulerian, numerical simulation is used to model the launching of plasma waves in a non-neutral plasma that is confined in a Penning–Malmberg trap. The waves are launched by applying an oscillating potential to an electrically isolated sector at one end of the conducting cylinder that bounds the confinement region and are received by another electrically isolated sector at the other end of the cylinder. The launching of both Trivelpiece–Gould waves and electron acoustic waves is investigated. Adopting a stratagem, the simulation captures essential features of the finite length plasma, while retaining the numerical advantages of a simulation employing periodic spatial boundary conditions. As a benchmark test of the simulation, the results for launched Trivelpiece–Gould waves of small amplitude are successfully compared to a linearized analytic solution for these fluctuations.

Physics↗

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↗

A buffer-gas trap for the NEPOMUC positron beam: optimization studies with electrons

Buffer-gas traps (BGTs) use inelastic interactions with nitrogen molecules to capture positrons from a continuous beam. These devices are invaluable for high-resolution studies of matter–antimatter interactions, antihydrogen research and positronium laser spectroscopy. We present a new project with the goal of producing a non-neutral plasma containing ${\sim }10^8$ low-energy positrons by installing a BGT on the NEPOMUC (NEutron induced POsitron source MUniCh) high-intensity positron beam. Details of the BGT are outlined and results are presented from experiments in which an electron beam, with a similar intensity and energy spread to the remoderated NEPOMUC beam, was used to create pulses of non-neutral electron plasma. The device is a vital component of the APEX (A Positron Electron eXperiment) project, which aims to create a low-temperature electron–positron pair plasma.

Physics↗

Modeling the Extended Neutral Atmosphere and Plasma Environment near Saturn

In the three years of this study we have published five papers in refereed journals. We have first examined satellite sources and their contribution to the observed neutral cloud. Based on the total calculated satellite sources and the spatial distribution of sputtered H20, we concluded that they cannot produce observed OH cloud. The E-ring contribution has been also studied in detail.In order to produce observed OH cloud we suggested that the E-ring might be the dominant source in inner Saturnian magnetosphere. We proposed a possible resupply mechanism which is needed to keep both E-ring and OH cloud in the present state: collisions between E-ring grains and remains of a disrupted satellite near Enceladus Lagrangian points. In this scenario a large amount of material, ranging from a few micrometers to hundred of meters, which is optically invisible at present, is likely to exist there. The fourth paper compares the magnetosheaths of the outer planets. A surprising result is that the hot proton component comprises about 40% of the total density, much larger than predicted by shock theory. Gas dynamic models of the boundaries show that the magnetospheres of Jupiter and Saturn are flattened at the poles. The last paper was published in GRL and is the first based of the model of neutrals developed as a main goal of this grant and which is now operational.This Monte Carlo collision code self- consistently determines the neutral distribution. from the rings and satellites until they are lost by ionization, by collisions with rings, moons, or Saturn, or by escape from Saturn. Our model is unique in that it includes the effects of plasma chemistry and both plasma-neutral and neutral neutral collisions to determine the dynamical evolution of the water group neutrals in Saturn's magnetosphere. The lifetimes of the neutrals against loss to photoionization, charge exchange, electron dissociation and electron impact dissociation are based on the model given by (Richardson et al. 1998) and vary with position in the magnetosphere. The dominant neutral dissociation channels H20->OH+H, H20-> O+H2, and OH->O+H2 are considered.

Richardson, John D.↗

SOLEDGE3X full vessel plasma boundary simulations of ITER non-active phase plasmas

The onset of detachment in the ITER machine is analyzed in this work through the help of 2D-axisymmetric boundary plasma simulations with the SOLEDGE3X-EIRENE code, which features a numerical domain for the plasma solver extending up to the first wall. The plasma boundary is computed in scenarios from the first non-active phase of ITER, in pure H and at 20 MW. This set of simulations is used in two aspects: first, to study the plasma detachment in the divertor, and second, the plasma conditions, fluxes, and beryllium erosion at the first wall. Here, the code results are also compared to those obtained with the well-established SOLPS-ITER code, which includes a plasma numerical domain only covering the main SOL. Results show an increase in the SOL width λ q with increasing density, and a detailed analysis is carried out, for the first time, on each of the different plasma-neutral interactions in the code’s physics model in EIRENE. The gross beryllium erosion rates of first wall panels are estimated from 2D simulations, with the aim of assessing their sensitivity to two parameters: the divertor density regime, and the presence of density shoulders in the far-SOL formed by enhanced perpendicular transport at this location. The erosion contributions from neutrals and ions are considered in each case, and the charge-exchange atoms fluxes and energy distributions are provided, highlighting the two atom populations (cold and charge-exchange).

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Transverse deflection and dissipation of small plasma beams and clouds in magnetized media

Propagation of a quasi-neutral plasma beam or cloud across a magnetic field is considered for the case where the transverse dimension of the beam or cloud is sufficiently small compared to ion gyroradii. This situation commonly arises for active experiments in near-earth space. Two mechanisms are presented for transverse deflection of a beam or cloud in the -v0 x B0 direction where v0 is the velocity relative to the ambient medium. In the first, asymmetric escape of ions from an electrically polarized beam or cloud causes transverse deflection by means of a rocket effect. The transverse deflection distance is estimated to be a few times the initial transverse dimension of the beam or cloud. Dissipation occurs within a few times the thermal ion transverse crossing time. In the second mechanism, asymmetric charging results from localized accumulation of incident ions from the ambient medium. This excess positive charge distorts electric equipotentials and drives electron Hall currents that maintain an asymmetric compressed magnetic field region. The asymmetry of the magnetic stress contributes to transverse deflection with the same sign as the rocket effect. The asymmetric magnetic field also focuses incident ions to yield the localized charge accumulation. These ideas are qualitatively consistent with observations of the Active Magnetospheric Particle Tracer Explorers artificial comet releases.

Cheng, Andrew F.↗

Stability of the Io torus

The inherent stability of the coupled neutral-plasma Io torus is investigated. Under the assumptions that the neutrals arise through sputtering by corotating torus ions, and that the ions move radially by centrifugally driven diffusion, it is found that the equilibrium state properties of the Io torus depend on two quantities: the sputtering yield, and Jupiter's ionospheric Pedersen conductance. For fixed values of these quantities, the torus equilibrium state is stable. When the torus neutral and ion densities deviate from the equilibrium state, the perturbations decay without oscillating, with a decay time near 1 day for the neutrals and near 60 days for the ions. The growths of the neutrals and ions from seed populations are computed. The torus evolves directly to its equilibrium state.

Huang, T. S.↗

Time integrator agnostic charge conserving finite element PIC

Developing particle-in-cell (PIC) methods using finite element basis sets, and without auxiliary divergence cleaning methods, was a longstanding problem until recently. It was shown that if consistent spatial basis functions are used, one can indeed create a methodology that was charge conserving, albeit using a leapfrog time stepping method. While this is a significant advance, leapfrog schemes are only conditionally stable and time step sizes are closely tied to the underlying mesh. Ideally, to take full advantage of advances in finite element methods (FEMs), one needs a charge conserving PIC methodology that is agnostic to the time stepping method. This is the principal contribution of this paper. In what follows, we shall develop this methodology, prove that both charge and Gauss’ laws are discretely satisfied at every time step, provide the necessary details to implement this methodology for both the wave equation FEM and Maxwell solver FEM, and finally demonstrate its efficacy on a suite of test problems. The method will be demonstrated by single particle evolution, non-neutral beams with space-charge, and adiabatic expansion of a neutral plasma, where the Debye length has been resolved, and real mass ratios are used.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Measurement of neutral gadolinium oscillator strengths using dual-comb absorption spectroscopy in laser-produced plasmas

The oscillator strengths of 43 transitions of neutral gadolinium from 530 to 535 nm were determined based on high-resolution absorption spectra measured by probing a laser-produced plasma (LPP) using time-resolved dual-comb spectroscopy (DCS). Absorption spectra were measured at various time delays, ranging from 33 to 252 µs after the onset of the plasma formation, over which the temperature and number density varied as the ablation plume evolved. A Boltzmann analysis was used to determine the excitation temperature and column density of Gd in the ablation plume, based on measured absorption from 20 Gd I spectral lines with known oscillator strengths. Oscillator strengths were then determined for additional dipole-allowed Gd I transitions identified in the absorption spectra, 19 of which were previously unreported. Oscillator strengths of all measured transitions were compared with literature values when available and show good agreement in most cases. Time-resolved measurement of atomic absorption in the cooling plasma provides access to a range of excitation temperatures and optical densities for the atoms in the ablation plume, allowing repeated measurements under different conditions, thereby improving confidence in the results. The high spectral bandwidth and resolution capabilities of DCS make it well-suited to measuring the dense optical spectrum of Gd, and the technique is applicable to other elements, especially lanthanides and actinides.

(280.4788) Optical Sensing, (280.5395) Plasma Diag↗

Simulations of ion heating due to ion-acoustic instabilities in presheaths

Particle-in-cell, direct simulation Monte Carlo simulations reveal that ion-acoustic instabilities excited in presheaths can cause significant ion heating. Ion-acoustic instabilities are excited by the ion flow toward a sheath when the neutral gas pressure is small enough and the electron temperature is large enough. A series of 1D simulations were conducted in which neutral plasma (electrons and ions) was uniformly sourced with an ion temperature of 0.026 eV and different electron temperatures (0.1 eV–50 eV). Ion heating was observed when the electron-to-ion temperature ratio exceeded the minimum value predicted by linear response theory to excite ion-acoustic instabilities at the sheath edge (T e /T i ≈ 28). When this threshold was exceeded, the temperature equilibration rate between ions and electrons rapidly increased near the sheath so that the local temperature ratio did not significantly exceed the threshold for instability. This resulted in significant ion heating near the sheath edge, which also extended back into the bulk plasma; presumably due to wave reflection from the sheath. Furthermore, this ion-acoustic wave heating mechanism was found to decrease for higher neutral pressures, where ion-neutral collisions damp the ion-acoustic waves and ion heating is instead dominated by inelastic collisions in the presheath.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Magnetohydrodynamics (MHD) Aerocapture System for Enabling Faster-Larger Planetary Science & Human Exploration Missions

Since our completing the NIAC Phase I NIAC on this Advanced Aerocapture System, NASA Langley Research Center has funded or supported a number of studies and code enhancements through its Center Innovation Fund (CIF) and NASA’s NSTGRO and Internship Programs to mature the analysis capabilities and quantify the merits of the MHD Aerocapture System technology. These efforts have resulted in a plug and play analysis capability for assessing MHD aerocapture system performance for arrival at many planetary bodies of interest. Our efforts have especially focused on the potential mass savings for improving the capacity for science observations at Neptune and Triton. A re-cent Forbes article published “‘Orbital mechanics is probably going to decide for us whether we go to Uranus or Neptune because we need to flyby Jupiter,’ said Kunio Sayanagi at Hampton University, Virginia, who also worked on the Neptune Odyssey proposal…. Exactly when a mission can be sent to Uranus, or Neptune, depends on the relative position of Jupiter, which can help give a spacecraft a gravitational slingshot. That drastically shortens the cruise phase.” [1] Since shortening the cruise phase is important for these science missions, any mass savings enabled by the MHD Aerocapture System could be reallocated to increasing Thermal Protection System mass to allow faster arrival speeds and/or for onboarding additional payloads such as science instruments, batteries, or propellant for conducting more science for longer durations in the desired orbits. The analysis steps and codes for conducting trades and sizing vehicles for aerocapture are as follows: Step 1 is to conduct aeroheating analysis using LAURA of the selected entry vehicle shape to identify locations on the forebody where ionization and flow velocity are sufficient for producing Lo-rentz forces. LAURA is a multiblock structured grid finite-volume CFD solver developed at the NASA Langley Research Center. [2] LAURA has been used for aerothermal analysis support of the entry, de-scent and landing (EDL) phase of interplanetary missions over the last three decades [3-7]. Step 2 is to port the LAURA results into CFDWARP to calcu-late electrical and thermal conductivities of ionized flow for sizing MHD patch system and calculating Lorentz forces needed for controls analysis. CFDWARP is a CFD code that uses advanced nu-merical methods that enable the simulation of the full coupling between the aerodynamics, the magne-tohydrodynamics, and the non-neutral plasma sheaths. CFDWARP has the unique capability to simulate efficiently the non-neutral sheaths (near the electrodes) in coupled form with the quasi-neutral bulk MHD flow [8-11]. Step 3 is to link re-sults from LAURA and CFDWARP into POST2 for calculating entry trajectories and comparing MHD control results with other aerodynamic control strategies. The Program to Optimize Simulated Tra-jectories II (POST2) is a generalized point mass, discrete parameter targeting and optimization pro-gram. POST2 provides the capability to target and optimize point mass trajectories for multiple pow-ered or un-powered vehicles near an arbitrary rotat-ing, oblate planet [12]. Step 4: TPS sizing was per-formed using the Fully Implicit Ablation and Ther-mal-response code (FIAT) tool which computes the transient one-dimensional thermal response and surface thermochemistry of a multilayer stackup of thermal protection, bonding, and structural materi-als subject to aeroheating on one surface [13]. The sizing and margining methodology used was based on the approach documented by Mahzari and Milos [14] for the dual-layered heatshield for extreme entry environment technology (DL-HEEET) TPS concept. TPS analysis utilizes trajectory information from POST2. Using this step-wise plug and play MHD Aerocapture performance assessment process, our analysis targets a Neptune aerocapture trajectory that will place the spacecraft in an observation orbit for Triton. [15]. Magnetohydrodynamic (MHD) control of a 4.5-meter diameter MSL-style capsule resulted in TPS mass savings of nearly 2000 kg when using an MHD system mass of under 200 kg. The flight path for a vehicle using the MHD control strategy has a much lower heat rate and heat load compared to the conventional aerodynamic aerocapture strategies known as bank angle con-trolled (BAC) and direct force controlled (DFC). Both BAC and DFC have heat rates significantly greater than 1500 W/cm2 typically used as an upper limit for PICA. Thus, DL-HEEET TPS concept was required for the BAC and DFC control strategies. However, considering the more benign environ-ments for the MHD case, additional TPS concepts with improved mass efficiency were also assessed. PICA was considered for the MHD controlled strat-egy since the maximum heat rate was well within the limits (<1500 W/cm2) of PICA. TPS sizing re-sulted in a significant mass reduction. The PICA layer for this sizing case was about 7.8 cm. As a point of reference, the Mars 2020 mission, which used this same PICA concept, had a PICA thickness of 3.18 cm [16]. The trajectories used for the TPS sizing originat-ed from the POST2 simulations. The current, I, to an electromagnet configuration can be manipulated to allow for active control of the vehicle. Manipula-tion of the current, I, changes the magnetic field, B, which affects the Lorentz force and therefore the MHD drag force on the vehicle. Our analysis in-cluded both open-loop and close-loop control. Closed-loop control will enable improved overall performance when taking into account mission level uncertainties, such as interplanetary delivery errors and atmospheric modeling uncertainties. The open-loop and closed-loop MHD control cases do not dip as deep into the atmosphere as the aerodynamic cases. Three types of aerodynamic-only approach-es are investigated: bank angle modulation (BAM), director force control (DFC), and Drag Modulated. BAM and DFC make use of vehicle aerodynamic angles to steer the vehicle. Thus, changing the aer-odynamic forces acting on the vehicle for control, aerodynamic drag modulated case requires a vary-ing drag area to modulate the drag force. The MHD drag modulated case modulates MHD generated drag force that adds to the aerodynamic drag. This higher atmospheric activation of drag forces by the MHD patch results in significantly less heat flux on the vehicle. The MHD technology will enable shorter cruise times and deceleration of larger payloads for increasing the capacity for science at the Ice Giants or for returning astronauts to Earth from cislunar space or from Mars. The purpose of this presentation is to provide more details about this work and to highlight plans for further research and development including a flight demonstration.

R. W. Moses↗

Rocket/Radar Sporadic-E Experiment Conducted during the El Coqui 2 Campaign

In order to investigate the complex electrodynamics and neutral-plasma coupling inherent to sporadic-E layers in the earth's mid-latitude ionosphere, a series of rocket/radar experiments were planned as part of the NASA El Coqui H Campaign from Tortuguero Launch Range, Puerto Rico, in March-April, 1998. The rocket experiments consisted of two pairs of "mother-daughter" payloads with limited apogees so that the payloads "hovered" in the sporadic-E region (95-125 km). Each payload pair included vector DC and AC electric field detectors, a highly accurate flux-gate DC magnetometer, an ion mass spectrometer, an ionization gauge, and spaced-electric field receivers to measure the wavelength and phase velocity of the unstable plasma waves. Separate rockets were included to simultaneously carry aloft TMA trails to measure the neutral wind and its velocity shear, believed responsible for the sporadic-E layer formation. In addition to the rocket experiments, incoherent scatter radar measurements of plasma density and drift velocity were gathered almost every night during the 3 week campaign. Continuous VHF backscatter radar operations were carried out from a site near Salinas, Puerto Rico, where 3-m backscatter echoes were observed associated with sporadic-E and other types of low altitude ionospheric layers. Other radars that operated during the campaign included an HF backscatter system near Ponce, Puerto Rico, and a second VHF backscatter radar set up near Aguadila Puerto Rico. On 24 March 1998, one of the instrumented rockets was launched, attaining an apogee of 129 km. The payloads successfully pierced an intense sporadic-E layer observed by both the Arecibo radar and the in-situ density and ion mass spectrometer probes. In-situ DC electric fields revealed very low (about 1-2 mV/m) ambient fields with small amplitude structures of the same order. No high frequency (short scale) waves were observed, consistent with the VHF backscatter observations at the time of the launch. An overview of the observations will be presented.

Pfaff, R.F.↗

Two-Dimensional Transport Studies for the Composition and Structure of the Io Plasma Torus

Research efforts in the second quarterly period have been focused primarily upon reviewing and assessing the neutral-plasma reactions and the plasma-plasma reactions that are important in determining the production and loss rates for the primary heavy ion species S(+), S(++), S(+++), O(+), and O(++) in the plasma torus and in acquiring new and updating old cross sections for the important processes.

Smyth, William H.↗

A reduced-model (nSOLT) simulation of neutral recycling effects on plasma turbulence in the divertor region of MAST-U

The 2D scrape-off-layer turbulence code (nSOLT), which includes 1D kinetic neutral–plasma interactions, is applied to study effects of neutral recycling on plasma turbulence for parameters illustrative of the MAST-U divertor region. Neutral recycling is modeled by injecting a fraction of the parallel plasma flux to the divertor back into the simulation domain as a source of Franck–Condon-distributed neutrals. Stationary sources, concentrated at the magnetic separatrix, model plasma streaming into the divertor region from the upstream scrape-off-layer and sustain plasma turbulence absent neutral recycling. Starting from one such no-neutrals equilibrium, we initiate recycling in a numerical experiment designed to diagnose and identify the effects of various neutral–plasma interactions on the divertor plasma, divertor turbulence, and plasma exhaust. The onset of recycling triggers an initial burst of enhanced cross field plasma transport that is quelled by ionization cooling and charge–exchange (CX) friction, with growing neutral pressure, leading to a quiescent, turbulence-free state. Diagnosis of this transient burst reveals that (1) the sudden increase in plasma density due to ionization dominates the onset of the burst, (2) electron cooling due to ionization increases collisionality and disconnects blob filaments from the sheath, and (3) CX friction drives tripole polarization of a blob that can dominate the curvature-driven dipole polarization, leading to the stagnation of blob propagation and reduced radial turbulent transport. It is shown that CX friction is negligible compared to sheath physics in determining equilibrium mean flow shearing rates, for parameters considered herein (specifically a short connection length to the divertor target), while it can significantly reduce interchange-instability growth rates.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Enhanced tuneable rotatory power in a rotating plasma

We report that the gyrotropic properties of a rotating magnetized plasma are derived analytically. Mechanical rotation leads to a new cutoff for wave propagation along the magnetic field, and polarization rotation above this cutoff is the sum of the classical magneto-optical Faraday effect and the mechanico-optical polarization drag. Exploiting the very large effective group index near the cutoff, we expose here that polarization drag can be 10 4 larger than Faraday rotation at GHz frequency. The rotation leads to weak absorption while allowing direct frequency control, demonstrating the unique potential of rotating plasmas for nonreciprocal elements. The very large rotation frequency of a dense non-neutral plasma could enable unprecedented gyrotropy in the THz regime.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗