Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “plasma properties”

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,135 records · Page 63

A flexible and differentiable coil proxy for stellarator equilibrium optimization

Balancing plasma performance and coil cost is a significant challenge when designing a stellarator power plant. Most current stellarator designs are produced through two-stage optimization: stage-1 for the equilibrium and stage-2 for a coil design that reproduces its magnetic configuration. Because few proxies connect both stages, two-stage optimization can produce plasmas that have high-quality physical properties but overly complex coils. In recent years, single-stage optimization has increasingly been used to optimize the plasma and coils simultaneously in order to improve the plasma–coil balance. However, all existing single-stage tools are specialized for filament coils, cannot model coil systems containing permanent magnets or dipole arrays, and continue to be challenged by numerical problems. The quasi-single-stage optimization finds a middle ground by integrating a coil optimization subproblem into stage-1 optimization. We present a flexible, differentiable coil complexity proxy based on the newly developed QUADCOIL coil optimization code. QUADCOIL is fast and can target realistic coil metrics and constraints that are unavailable to codes with comparable speed. We demonstrate the effectiveness and flexibility of the QUADCOIL proxy by presenting two QSS optimization studies. The first study produces a permanent magnet solution for the MUSE stellarator with 34% fewer magnets than previous solutions. The second study produces a coil solution for the ARIES-CS stellarator with 30% reductions in both peak and root-mean-square force.

Fu, Lanke [Princeton Plasma Physics Laboratory (PP↗

Invariant criteria for bound states, degree of ionization, and plasma phase transition

Basis invariant characterizations of bound states and bound fraction of a partially ionized hydrogen plasma are given in terms of properties of the spectrum of eigenvalues and eigenfunctions of the equilibrium quantum statistical one-proton-one-electron reduced density matrix. It is suggested that these can be used to place theories of a proposed plasma-ionization phase transition on a firm foundation. This general approach may be relevant to cosmological questions such as the quark deconfinement-confinement transition.

Girardeau, M. D.↗

Linear gyrokinetic stability of a high β non-inductive spherical tokamak

Spherical tokamaks (STs) have been shown to possess properties desirable for a fusion power plant such as achieving high plasma β and having increased vertical stability. To understand the confinement properties that might be expected in the conceptual design for a high β ST fusion reactor, a 1 GW ST plasma equilibrium was analyzed using local linear gyrokinetics to determine the type of micro-instabilities that arise. Kinetic ballooning modes and micro-tearing modes are found to be the dominant instabilities. The parametric dependence of these linear modes was determined and, from the insights gained, the equilibrium was tuned to find a regime marginally stable to all micro-instabilities at θ 0 = 0.0. This work identifies the most important micro-instabilities expected to generate turbulent transport in high β STs. The impact of such modes must be faithfully captured in first-principles-based reduced models of anomalous transport that are needed for predictive simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Fabrication and thermophysical properties of UO 2 -UB 2 and UO 2 -UB 4 composites sintered via spark plasma sintering

Uranium dioxide (UO 2 ) composites with uranium diboride (UB 2 ) and uranium tetraboride (UB 4 ) have been proposed as advanced fuel candidates due to their high thermal conductivity, high melting point, high fissile density and their ability to incorporate a built-in burnable poison by tailoring the targeted 10 B/ 11 B ratio. As such, it is important to assess the fabrication, and thermal and micromechanical properties of such composites. In this work, UO 2 -UB 2 and UO 2 -UB 4 samples with boride phase fractions of 5, 15 and 30 wt% were fabricated to high densities (above 95 % theoretical density) via spark plasma sintering (SPS). This enabled sintering at relatively low temperatures and short timescales. SPS also aided in maintaining the target phase fractions of the samples as reactions between the constituent phases were suppressed due to the short timescales and reducing environment during sintering. Here, thermal diffusivity measurements from 299 to 1273 K were conducted through laser flash analysis (LFA). The diffusivity increased as a function of boride weight fraction, and UB 2 additions increased the thermal diffusivity of the composites more than UB 4 additions. Assessment of the LFA results indicated in-situ reactions between the UO 2 and boride phases that suppress the thermal diffusivity occur above 800 K for all samples. Oxidation of the boride phase was proposed as the underlying reaction. This was supported by thermodynamic assessments from the literature, as well as microstructural, crystallographic, and nanoindentation characterization performed on these samples.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Properties of AKR from an interferometer analysis of the ISEE-1 and -2 plasma wave data

Both ISEE 1 and 2 are equipped with identical wide-band analog plasma wave receivers which makes long baseline interferometry possible. A cross-correlation of the received signals from each spacecraft during periods of intense Auroral Kilometric Radiation (AKR) provides information on the source size of AKR. From the magnitude of the cross-correlation coefficient, measured over baselines from 40 Km to 8000 Km, an upper limit to the statistical source size can be determined. Thirty-five AKR events at 125 kHz are being analyzed. Preliminary results indicate that the apparent source sizes are smaller than .1 Re.

Baumback, M. M.↗

On the mechanism of ionization oscillations in Hall thrusters

Low-frequency ionization oscillations involving plasma and neutral density (breathing modes) are the most violent perturbations in Hall thrusters for electric propulsion. Because of its simplicity, the zero-dimensional (0D) predator–prey model of two nonlinearly coupled ordinary differential equations for plasma and neutral density has often been used for the characterization of such oscillations and scaling estimates. In this work, we investigate the properties of its continuum analog, the one-dimensional (1D) system of two nonlinearly coupled equations in partial derivatives (PDEs) for plasma and neutral density. This is a more general model, of which the standard 0D predator–prey model is a special limit case. We show that the 1D model is stable and does not show any oscillations for the boundary conditions relevant to Hall thrusters and the uniform ion velocity. We then propose a reduced 1D model based on two coupled PDEs for plasma and neutral densities that is unstable and exhibit oscillations if the ion velocity profile with the near-the-anode back-flow (toward the anode) region is used. Comparisons of the reduced model with the predictions of the full model that takes into account the self-consistent plasma response show that the main properties of the breathing mode are well captured. In particular, it is shown that the frequency of the breathing mode oscillations is weakly dependent on the final ion velocity but shows a strong correlation with the width of the ion back-flow region.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Relativistic, perpendicular shocks in electron-positron plasmas

One-dimensional particle-in-cell plasma simulations are used to examine the mechanical structure and thermalization properties of collisionless relativistic shock waves in electron-positron plasmas. Shocks propagating perpendicularly to the magnetic field direction are considered. It is shown that these shock waves exist, and that they are completely parameterized by the ratio of the upstream Poynting flux to the upstream kinetic energy flux. The way in which the Rankine-Hugoniot shock jump conditions are modified by the presence of wave fluctuations is shown, and they are used to provide a macroscopic description of these collisionless shock flows. The results of a 2D simulation that demonstrates the generality of these results beyond the assumption of the 1D case are discussed. It is suggested that the thermalization mechanism is the formation of a synchrotron maser by the coherently reflected particles in the shock front. Because the downstream medium is thermalized, it is argued that perpendicular shocks in pure electron-positron plasmas are not candidates as nonthermal particle accelerators.

Gallant, Yves A.↗

Voltage waveform tailoring for high aspect ratio plasma etching of SiO 2 using Ar/CF 4 /O 2 mixtures: Consequences of low fundamental frequency biases

The use of non-sinusoidal waveforms in low pressure capacitively coupled plasmas intended for microelectronics fabrication has the goal of customizing ion and electron energy and angular distributions to the wafer. One such non-sinusoidal waveform uses the sum of consecutive harmonics of a fundamental sinusoidal frequency, f 0 , having a variable phase offset between the fundamental and even harmonics. In this paper, we discuss results from a computational investigation of the relation between ion energy and DC self-bias when varying the fundamental frequency f 0 for capacitively coupled plasmas sustained in Ar/CF 4 /O 2 and how those trends translate to a high aspect ratio etching of trenches in SiO 2 . The fundamental frequency, f 0 , was varied from 1 to 10 MHz and the relative phase from 0° to 180°. Two distinct regimes were identified. Average ion energy onto the wafer is strongly correlated with the DC self-bias at high f 0 , with there being a maximum at φ = 0° and minimum at φ = 180°. In the low frequency regime, this correlation is weak. Average ion energy onto the wafer is instead dominated by dynamic transients in the applied voltage waveforms, with a maximum at φ = 180° and minimum at φ = 0°. The trends in ion energy translate to etch properties. In both, the high and low frequency regimes, higher ion energies translate to higher etch rates and generally preferable final features, though behaving differently with phase angle.

36 MATERIALS SCIENCE↗

Near-Surface Plasma Characterization of the 12.5-kW NASA TDU1 Hall Thruster

To advance the state-of-the-art in Hall thruster technology, NASA is developing a 12.5-kW, high-specific-impulse, high-throughput thruster for the Solar Electric Propulsion Technology Demonstration Mission. In order to meet the demanding lifetime requirements of potential missions such as the Asteroid Redirect Robotic Mission, magnetic shielding was incorporated into the thruster design. Two units of the resulting thruster, called the Hall Effect Rocket with Magnetic Shielding (HERMeS), were fabricated and are presently being characterized. The first of these units, designated the Technology Development Unit 1 (TDU1), has undergone extensive performance and thermal characterization at NASA Glenn Research Center. A preliminary lifetime assessment was conducted by characterizing the degree of magnetic shielding within the thruster. This characterization was accomplished by placing eight flush-mounted Langmuir probes within each discharge channel wall and measuring the local plasma potential and electron temperature at various axial locations. Measured properties indicate a high degree of magnetic shielding across the throttle table, with plasma potential variations along each channel wall being less than or equal to 5 eV and electron temperatures being maintained at less than or equal to 5 eV, even at 800 V discharge voltage near the thruster exit plane. These properties indicate that ion impact energies within the HERMeS will not exceed 26 eV, which is below the expected sputtering threshold energy for boron nitride. Parametric studies that varied the facility backpressure and magnetic field strength at 300 V, 9.4 kW, illustrate that the plasma potential and electron temperature are insensitive to these parameters, with shielding being maintained at facility pressures 3X higher and magnetic field strengths 2.5X higher than nominal conditions. Overall, the preliminary lifetime assessment indicates a high degree of shielding within the HERMeS TDU1, effectively mitigating discharge channel erosion as a life-limiting mechanism.

Lifetime↗

Accelerated Synthesis of Ordered Mesoporous Carbons Using Plasma

Conventional ordered mesoporous carbon (OMC) production usually requires long processing times in the carbonization step to achieve desired temperatures through controlled ramps. To enable expedited materials discovery, developing advanced manufacturing capability with significantly improved throughput is highly desired. Current approaches for accelerating the synthesis of OMCs include using microwave and Joule heating. However, both methods rely on the introduction of additional components, such as microwave absorbers and electrically conductive agents, within the bulk materials to impart the ability to reach high carbonization temperatures. This work demonstrates accelerated synthesis and functionalization of OMCs through the use of a dielectric barrier discharge plasma, where carbonization can be accomplished within 15 min using 30 W plasma sources, representing more than an order of magnitude increase in polymer-to-carbon conversion kinetics compared to that of a traditionally pyrolyzed analogue. Particularly, the ability of performing rapid carbonization without the use of additional substrates within the OMC precursor systems is advantageous. A systematic investigation of how plasma power, time, and gas atmosphere impact the resulting OMC pore textures and properties is performed, demonstrating the broad applicability of plasma-enabled carbonization methods. Furthermore, we demonstrate that the plasma treatment strategy can be extended to incorporate heteroatoms into the carbon framework by introducing ammonia gas, resulting in OMCs with a nitrogen content up to 4.7 at %, as well as non-Pluronic templating systems for synthesizing OMC with pore sizes larger than 10 nm. As employing a plasma source for materials pyrolysis is an industrially relevant approach, our system can be extended toward scaled synthesis of OMCs with much faster production rates.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Jet quenching and medium response in high-energy heavy-ion collisions: a review

Jet quenching has been used successfully as a hard probe to study properties of the quark-gluon plasma (QGP) in high-energy heavy-ion collisions at both the relativistic heavy-ion collider and the large hadron collider. Furthermore, we will review recent progresses in theoretical and phenomenological studies of jet quenching with jet transport models. Special emphasis is given to effects of jet-induced medium response on a wide variety of experimental observables and their implications on extracting transport properties of the QGP in heavy-ion collisions.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

dc electrical conductivity in strongly magnetized plasmas

A generalized Ohm's law is derived to treat strongly magnetized plasmas in which the electron gyrofrequency significantly exceeds the electron plasma frequency. Here, strong magnetization of electrons causes the frictional drag between electrons and ions due to Coulomb collisions to shift, producing an additional transverse resistivity term in the generalized Ohm's law that is perpendicular to both the current ($J$) and the Hall ($JxB$) direction. In the limit of very strong magnetization the parallel resistivity is found to increase by a factor of 3/2, and the perpendicular resistivity by a factor of $\frac{3}{4} \ln (m_i/m_e)$, where $m_i$ and $m_e$ are the ion and electron masses. These results suggest that strong magnetization significantly changes the magnetohydrodynamic evolution of a plasma.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Validation of a comprehensive first-principles-based framework for predicting the performance of future stellarators

This paper presents the validation of the GENE-KNOSOS-Tango framework for recovering both the steady-state plasma profiles in the considered radial domain and selected turbulence trends in a stellarator. This framework couples the gyrokinetic turbulence code GENE, the neoclassical transport code KNOSOS, and the transport solver Tango in a multi-timescale simulation feedback loop. Ion-scale kinetic-electron and electron-scale adiabatic-ion flux-tube simulations were performed to evolve the density and temperature profiles for four OP1.2b W7-X scenarios. The simulated density and temperature profiles showed good agreement with the experimental data using a reasonable set of boundary conditions. Equally important was the reproduction of observed trends for several turbulence properties, such as density fluctuations and turbulent heat diffusivities. Key effects were also touched upon, such as electron-scale turbulence and the neoclassical radial electric field shear. The validation of the GENE-KNOSOS-Tango framework enables credible predictions of physical phenomena in stellarators and reactor performance based on a given set of edge parameters.

Gyrokinetic simulations↗

Near-Net-Shape HIP Manufacturing for sCO2 Turbomachinery Cost Reduction

sCO2 turbomachinery that operates above 650°C requires the use of γ' strengthened Ni-based superalloys, leading to high cost and barrier of market adoption. Near-net-shape (NNS) HIP manufacturing with 282® alloy powder is being developed for sCO2 turbine components, with a significant estimated cost reduction. Tensile, creep, low cycle fatigue properties of argon gas atomized and plasma atomized powders were evaluated and compared to sand cast HAYNES® 282®. While tensile strength and fatigue life outperformed sand cast material, a 10~25% debit in the creep stress capability was observed due to the fine grain size and presence of prior particle boundaries (PPBs). Finite element model calibrated by powder rheological properties accurately predicted the nonuniform shrinkage during HIP, providing HIP tooling design to achieve the target dimension. A 20 lbs. turbine nozzle ring was successfully demonstrated within 0.01 inch dimensional tolerance at the vanes. A 1700lbs. turbine casing with complex internal struts and manifolds was also demonstrated being close to the target dimension.

14 SOLAR ENERGY↗

ETG turbulence in a tokamak pedestal

This paper explores the fundamental characteristics of electron-temperature-gradient (ETG)-driven turbulence in the tokamak pedestal. The extreme gradients in the pedestal produce linear instabilities and nonlinear turbulence that are distinct from the corresponding ETG phenomenology in the core plasma. The linear system exhibits multiple (greater than ten) unstable eigenmodes at each perpendicular wave vector, representing different toroidal and slab branches of the ETG instability. Proper orthogonal decomposition of the nonlinear fluctuations reveals no clear one-to-one correspondence between the linear and nonlinear modes for most wave vectors. Moreover, nonlinear frequencies deviate strongly from those of the linear instabilities, with spectra peaking at positive frequencies, which is opposite the sign of the ETG instability. The picture that emerges is one in which the linear properties are preserved only in a narrow range of k-space. Outside this range, nonlinear processes produce strong deviations from both the linear frequencies and eigenmode structures. This is interpreted in the context of critical balance, which enforces alignment between the parallel scales and fluctuation frequencies. We also investigate the nonlinear saturation processes. We observe a direct energy cascade from the injection scale to smaller scales in both perpendicular directions. However, in the bi-normal direction, there is also nonlocal inverse energy transfer to larger scales. Neither streamers nor zonal flows dominate the saturation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Coatings for high-temperature bearings and seals

Criteria are discussed for predicting the probable lubricating ability of candidate solid materials from a consideration of their basic chemical and physical properties. The properties considered to be important in the model are thermochemical potential, adhesion, low hardness, plasticity, yield strength in shear, and brittle-to-ductile transition characteristics. A review of the selection and tribological testing of materials, which were selected for use in self-lubricating composite coatings by employing this model, is given. Two series of plasma-sprayed coatings with good tribological properties over a wide temperature spectrum are described. The PS 100 series of coatings contain oxidatively stable solid lubricants in a nichrome matrix. The PS 200 series contains the same solid lubricants in a very wear resistant metal-bonded chromium carbide matrix. Examples are given of applications of these coatings in high speed shaft seals, sliding contact bearings, and Stirling engine cylinder liners.

Sliney, Harold E.↗

A Particle-in-Cell Method for Plasmas with a Generalized Momentum Formulation, Part I: Model Formulation

Here, this paper formulates a new particle-in-cell method for the Vlasov–Maxwell system. Under the Lorenz gauge condition, Maxwell’s equations for the electromagnetic fields can be written as a collection of scalar and vector wave equations. The use of potentials for the fields motivates the adoption of a Hamiltonian formulation for particles that employs the generalized (conjugate) momentum. A notable advantage offered by the Hamiltonian formulation is the elimination of time derivatives in the Lorenz gauge formulation that are required by the standard Newton–Lorentz treatment of the particles. This allows the fields to retain the full time-accuracy guaranteed by the field solver. The resulting updates for particles require only knowledge of the fields and their spatial derivatives. An analytical method for constructing these spatial derivatives is presented that exploits the underlying integral solution used in the field solver for the wave equations. Moreover, these derivatives are demonstrated to converge at the same rate as the fields in both time and space. The Method of Lines Transpose field solver we consider in this work is globally first-order accurate in time and high-order accurate in space (e.g., fourth- and fifth-order) and belongs to a larger class of methods which are unconditionally stable, can address geometry, and leverage $\mathcal {O}(N)$ fast summation methods for efficiency. We demonstrate the method on several well-established benchmark problems on bounded domains, including a plasma sheath as well as a relativistic particle beam. The efficacy of the proposed formulation is established by comparing with a second-order accurate finite-difference time-domain method that employs a leapfrog time advance for particles and a charge conserving map suitable for bounded domains. The new method shows mesh-independent numerical heating properties even in cases where the plasma Debye length is smaller than the grid spacing. This is an important feature of the new method for problems defined on bounded domains, because it permits the use of coarser grids in space in the representation of the fields. Such a capability has significant implications for the simulation of plasmas in bounded domains with complex geometry, where the ratio between the largest and smallest cells can vary significantly. The use of high-order spatial approximations in the new method also means that fewer grid points are required in order to achieve a fixed accuracy. Our results also suggest that the new method can be used with fewer simulation particles per cell compared to the benchmark explicit method, which permits further computational savings.

97 MATHEMATICS AND COMPUTING↗