Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “plasma temperature”

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 577 records · Page 32
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

First tungsten radiation studies in DIII-D’s ITER baseline demonstration discharges

Abstract ITER Baseline Scenario plasmas were studied in DIII-D using krypton and xenon gases as a proxy for the tungsten that will be present in ITER. These impurities were chosen for having the same radiative loss rate Lz as tungsten would exhibit in the hotter ITER core. Results show that the scenario with these core radiators spans the range of impurity concentration and W radiated fraction expected for ITER, and up to 50% higher values, explored at zero injected torque, as well as 1 Nm and full co-torque injection with T ∼ 3 Nm. Stationary discharges with duration >2–4 τ R are achieved with f rad ⩾ 30% leading to a reduction in confinement of ∼10%, and a comparison with real metal radiators in the same range of f rad shows that the higher Lz at the lower temperatures in these plasmas yields too pessimistic results on the survivability and performance of this scenario in ITER. Simulations of ITER power balance including W radiation show that with concentration up to three times higher than in the DIII-D plasmas the scenario can be stationary, remaining at acceptable core radiated fraction values.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Numerical-heating effects in atmospheric pressure streamer discharges simulated with a PIC code

Artificial heating in plasma simulations is a well-known phenomenon which occurs when, among other things, the Debye length is poorly resolved by the simulation mesh. Here, in this work, the degree to which numerical-heating occurs during a simulation of a nanosecond atmospheric pressure streamer discharge is examined. The streamer is simulated using a two-dimensional finite-element, particle-in-cell code Empire, which uses direct simulation Monte Carlo for binary particle interactions. Initially, an estimate of the numerical-heating rate applied to Empire is performed using a simple plasma model. Second, a positive atmospheric pressure streamer discharge simulation is performed to study the effects of numerical heating on plasma density, electron temperature, and streamer velocity. The nominal Debye length is approximately 1 μm and the amount of numerical heating introduced in the simulation is varied by using mesh sizes ranging from 2 μm to 20 μm. A measurable numerical heating quantity is proposed that can be used to estimate the appropriate element size and quantify the numerical-heating that can be expected over the simulation time for an atmospheric pressure streamer. In conclusion while Δx/λ D violations can be an issue it is not likely to be an issue with streamer discharges that are temporally short and occur in environments where collision frequencies are high. This result validates the rationale of grid size choices for a large amount of previously published works where Δx/λ D violation was not clearly addressed. Primary finding of this work is that numerical heating is of minor concern for plasma simulations where electron–neutral collisions are numerous such that multiple collisions can occur within a single plasma period.

Nikic, Dejan [University of New Mexico, Albuquerqu↗
DOE OSTI2023

Jet thermalization in QCD kinetic theory

We perform numerical studies in the framework of QCD kinetic theory to investigate the energy and angular profiles of a high energy parton — as a proxy for a jet produced in heavy ion collisions — passing through a Quark-Gluon Plasma (QGP). We find that the fast parton loses energy to the plasma mainly via a radiative turbulent quark and gluon cascade that transports energy locally from the jet down to the temperature scale where dissipation takes place. In this first stage of the system time evolution, the angular structure of the turbulent cascade is found to be relatively collimated. However, when the lost energy reaches the plasma temperature it is rapidly transported to large angles w.r.t. the jet axis and thermalizes. We investigate the contribution of the soft jet constituents to the total jet energy. We show that for jet opening angles of about 0.3 rad or smaller, the effect is negligible. Conversely, larger opening angles become more and more sensitive to the thermal component of the jet and thus to medium response. Our result showcases the importance of the jet cone size in mitigating or enhancing the details of dissipation in jet quenching observables.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Plasma catalysis: separating plasma and surface contributions for an Ar/N 2 /O 2 atmospheric discharge interacting with a Pt catalyst

Atmospheric pressure non-equilibrium plasmas can form nitrogen oxide (NO x ) compounds directly from nitrogen and oxygen without a catalyst, and at lower catalyst temperatures than would be possible without plasma. In this work, the oxidation of plasma-produced NO from an Ar/N 2 /O 2 non-equilibrium atmospheric-pressure plasma-jet (APPJ) over a platinum-on-alumina powder catalyst was investigated with in-situ infrared spectroscopy. Products downstream of the catalyst bed were analyzed along with catalyst surface species. The catalyst was exposed to plasma at both constant temperature and a cyclic temperature ramp in order to study long-lasting and transient surface changes. Primary incident reactive species to the catalyst were assessed to be NO and O 3 . Pt-Al 2 O 3 at 350 °C increased oxidation of NO relative to Al 2 O 3 or an empty chamber. The surface state of Pt-Al 2 O 3 evolves during plasma-effluent exposure and requires upwards of 20 min exposure for stabilization compared to Al 2 O 3 . Once stable surface conditions are achieved, thermal cycling reveals a repeatable hysteresis pattern in downstream products. At low temperature, oxygen and NO x accumulate on the catalyst surface and react at elevated temperatures to form NO 2 . Increasing plasma power and O 2 :N 2 ratio increases the hysteresis of the heating relative to the cooling curves in the pattern of NO 2 formation. The limitation on NO oxidation at high temperatures was assessed to be Pt-O which is depleted as the catalyst is heated. Once stored species have been depleted, NO oxidation rates are determined by incoming reactants. Two overlapping NO oxidation patterns are identified, one determined by surface reactants formed at low temperature, and the other by reactants arriving at the surface at high temperature. As a result, the plasma is responsible for providing the reactants to the catalyst surface, while the catalyst enables reaction at high temperature or storage at low temperature for subsequent reaction.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Non-thermal atmospheric pressure plasma–liquid synthesis of organic acids in aqueous solution from carbon monoxide

This work aims at understanding the conversion of CO to organic acids, namely oxalic acid and formic acid, using non-thermal atmospheric pressure plasma over aqueous solutions. CO exhibited significantly higher conversion to organic acids (more than 15×) compared to CO 2 under the same reaction conditions. The result bolsters a proposed two-step process for CO 2 fixation, whereby CO 2 is first converted to CO, and then CO is converted to organic acids. The organic acids produced from CO are intermediates in the water–gas shift (WGS) reaction of CO in the presence of an aqueous solution to dissolved CO 2 and hydrogen gas. Based on a simple thermodynamic analysis, the organic acid yield was increased by lowering the plasma–liquid reaction temperature using an ice bath to cool the reaction flask. The composition of the organic acids could be varied by changing the pH of the solution. Oxalate was formed in higher concentrations with increasing solution pH above the pK a of the radical species (CO 2 )˙ − . Below the pK a value, formate was the exclusive organic acid formed. The production of formate has a rather weak pH dependence but is enhanced slightly at a basic pH above 10. Furthermore, at basic pH, the effect of electrolyte concentration comes into play. Higher electrolyte concentrations, leading to shorter electrolyte Debye lengths, resulted in lowered organic acid yields. The highest yields of organic acids obtained in our system were 122 mg L −1 for oxalate and 77 mg L −1 for formate at an optimum 1 mM NaOH concentration in the starting solution. This work is a successful pioneering example of CO to organic acids conversion using non-thermal plasmas, which opens the pathway for a promising two-step conversion process of CO 2 to organic acids.

Sudagar, Alcina Johnson [Washington University in ↗
DOE OSTI2022

Multi-species temperature and number density analysis of a laser-produced plasma using dual-comb spectroscopy

Dual-comb spectroscopy (DCS) represents a novel method of using absorption spectroscopy as a diagnostic tool for multispecies analysis of excitation temperatures and column densities in laser-produced plasmas (LPPs). DCS was performed on a LPP generated by ablating a multielement alloy containing Nd, Gd, and Fe. Transitions from all three elements were observed in absorption spectra measured from 530.08 to 535.19 nm at seven time-delays from 31 to 250 μs after ablation. The spectra were fit using a nonlinear regression algorithm to determine peak areas, and excitation temperatures and column densities were determined for the three atomic species separately using Boltzmann plots. The measured excitation temperatures of Nd I and Gd I showed good agreement at all time-delays, whereas the Fe I temperature was found to be higher, and the ratios between the column densities varied with delay. The observations are understood via effects of LPP spatial averaging, elemental fractionation, and molecular formation and are compared and contextualized with previous work studying LPPs using other spectroscopic techniques. A brief discussion of the precision and accuracy of the determined excitation temperatures and column densities is also presented.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Vapor-Phase Aggregation of Cerium Oxide Nanoparticles in a Rapidly Cooling Plasma

Local conditions, such as temperature and oxygen availability, have a pronounced effect on the formation and evolution of fallout following a nuclear explosion. While the behavior of nuclear-relevant materials such as uranium has begun to be explored under a wider range of environments, little is known about the behavior of plutonium. Here, using cerium as a surrogate, we track the vapor-phase aggregation of cerium oxide nanoparticles created in a plasma flow reactor under conditions of controlled temperature at two different oxygen fugacities. In situ optical emission spectroscopy is used to measure the variations in the spectral intensity of atomic and molecular species with temperature and oxygen content. We find that the relative rate of gas-phase oxidation of cerium is highly dependent on both temperature and local redox conditions within the flow reactor, to the extent that doubling the oxygen availability effectively doubles the amount of vapor-phase cerium monoxide at high temperatures (>2000 K). Condensed cerium oxide nanoparticles are also collected and analyzed ex situ via transmission electron microscopy and grazing-incidence small-angle X-ray scattering to determine their elemental composition, crystal structure, and size distribution. The size and morphology of the condensed nanoparticles are independent of local redox conditions, forming the same crystal type with the same size distribution regardless of oxygen availability. Postcondensation particle evolution, however, is found to be predominantly driven by temperature, with the average particle size increasing as particles cool and subsequently aggregate. These results expand our understanding of the chemical and physical behavior of refractory oxides that form during the early stages of fallout formation.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2026

Simulations of grain growth in tungsten armor materials under ARC plasma edge operation conditions using an integrated plasma-edge/materials model

An integrated model of grain growth deuterium-exposed tungsten polycrystals, consisting of a two-dimensional vertex dynamics model fitted to atomistic data, has been developed to assess the grain growth kinetics of deuterium-exposed polycrystalline tungsten (W). The model tracks the motion of grain boundaries under the effect of driving forces stemming from grain boundary curvature and differential deuterium concentration accumulation. Here, we apply the model to experimentally synthesized tungsten polycrystals under deuterium-saturated conditions relevant to the ARC concept design. The results indicate rapid grain growth kinetics in the near-surface region adjacent to the plasma, where the temperature reaches 1400 K, whereas the microstructure remains stable deeper in the material with the lower temperature of 1000 K. The combined modeling and analysis further reveal that monolithic tungsten produced via conventional fabrication routes is highly susceptible to grain coarsening at temperatures exceeding 1000 K, largely independent of the magnitude of the applied driving force. Moreover, the accumulation of deuterium near grain boundaries has a pronounced inhibitory effect on grain boundary migration. High-angle grain boundaries ( > 50°) contribute more significantly to the overall grain growth process.

36 MATERIALS SCIENCE↗
DOE OSTI2021

An Eulerian Vlasov-Fokker–Planck algorithm for spherical implosion simulations of inertial confinement fusion capsules

Here, we present a numerical algorithm that enables a phase-space adaptive Eulerian Vlasov–Fokker–Planck (VFP) simulation of inertial confinement fusion (ICF) capsule implosions. The approach relies on extending a recent mass, momentum, and energy conserving phase-space moving-mesh adaptivity strategy to spherical geometry. In configuration space, we employ a mesh motion partial differential equation (MMPDE) strategy while, in velocity space, the mesh is expanded/contracted and shifted with the plasma’s evolving temperature and drift velocity. The mesh motion is dealt with by transforming the underlying VFP equations into a computational (logical) coordinate, with the resulting inertial terms carefully discretized to ensure conservation. To deal with the spatial and temporally varying dynamics in a spherically imploding system, we have developed a novel nonlinear stabilization strategy for MMPDE in the configuration space. The strategy relies on a nonlinear optimization procedure that optimizes between mesh quality and the volumetric rate change of the mesh to ensure both accuracy and stability of the solution. Implosions of ICF capsules are driven by several boundary conditions: (1) an elastic moving wall boundary; (2) a time-dependent Maxwellian Dirichlet boundary; and (3) a pressure-driven Lagrangian boundary. Of these, the pressure-driven Lagrangian boundary driver is new to our knowledge. The implementation of our strategy is verified through a set of test problems, including the Guderley and Van-Dyke implosion problems — the first-ever reported using a Vlasov–Fokker–Planck model.

1D2V↗
DOE OSTI2020

Progress in DIII-D towards validating divertor power exhaust predictions

In this study, UEDGE simulations highlight the role of cross-field drifts on the onset of detached conditions, and new calibrated divertor vacuum ultra violet (VUV) spectroscopy is used to challenge the predictions of radiative constituents in these simulations. UEDGE simulations for DIII-D H-mode plasmas with the open divertor with the ion ∇B-drift towards the X-point show a bifurcated onset of the low field side (LFS) divertor detachment, consistent with experimentally observed step-like detachment onset (Jaervinen A.E. et al 2018 Phys. Rev. Lett. 121 075001). The divertor plasma in the simulations exhibits hysteresis in upstream separatrix density between attached and detached solution branches. Reducing the drift magnitude by a factor of 3 eliminates the step-like detachment onset in the simulations, confirming the strong role of drifts in the bifurcated detachment onset. When measured local plasma densities and temperatures are within proximity of predicted values in the simulations, there is no shortfall of the local emission of the dominant resonant radiating lines. However, the simulations systematically predict a factor of two lower total integrated radiated power than measured by the bolometer with the difference lost through radial heat flow out of the computational domain. Even though there is no shortfall in the emission of the dominant lines, a shortfall of total radiated power can be caused by underpredicted spatial extent of the radiation front, indicating a potential upstream or divertor transport physics origin for the radiation shortfall, or shortfall of radiated power in the spectrum between the dominant lines. In addition to the underpredicted spatial extent, in detached conditions, the simulations overpredict the peak radiation and dominant carbon lines near the X-point, which can be alleviated by manually increasing divertor diffusivity in the simulations, highlighting the ad hoc cross-field transport as one of the key limitations of the predictive capability of these divertor fluid codes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Autonomous Aerosol and Plasma Co‐Jet Printing of Metallic Devices at Ambient Temperature

Abstract Additive manufacturing of metallic materials holds the potential to revolutionize the fabrication of functional devices unattainable via traditional methods. Despite recent advancements, printing metallic materials typically requires thermal processing at elevated temperatures to form dense structures with desired properties, which presents a major challenge for direct printing and integration with temperature‐sensitive materials. Herein, a unique co‐jet printing (CJP) method is reported integrating an aerosol jet and a non‐thermal, atmospheric pressure plasma jet to enable concurrent aerosol deposition of metal nanoparticle inks and in situ sintering at ambient temperature. A machine learning algorithm is integrated with the CJP to perform real‐time defect detection and autonomous correction, enhancing the yield of printed films with high electrical conductivity from 44% to 94%. Concurrent printing and sintering eliminate the need for post‐printing processing, reducing the overall manufacturing time by multiple folds depending on product size. CJP enables direct printing of functional devices on a variety of temperature‐sensitive materials including biological materials. Direct printing of hydration sensors on living plant leaves is demonstrated for long‐duration monitoring of hydration level in the plant. The versatile CJP method opens tremendous opportunities to harmoniously integrate abiotic and biotic materials for emerging applications in wearable/implantable devices and biohybrid systems.

Du, Yipu [Department of Aerospace and Mechanical E↗
DOE OSTI2025

Particle balance of deuterium during deuterium shattered pellet injection shutdown in DIII-D

A particle balance analysis was conducted during a deuterium (D 2 ) shattered pellet injection-induced plasma shutdown on the DIII-D tokamak to determine why less than 20% of the pellet material is assimilated into the core plasma by the mid-current quench (CQ). Initially, most of the D 2 is injected as frozen shards and ionized upon entering the vessel. During the thermal quench, ionized particles move to the divertors and subsequently to the center post (CP) walls, where they rapidly recycle and partially accumulate as neutrals without assimilating into the core plasma. In contrast, the particle flux to the outer midplane walls is negligible, despite being accompanied by hot plasma with electron temperatures exceeding 100 eV. During mid-CQ, volume recombination effects, although not large enough to impact overall particle balance, were significant enough to require accounting for accurate interpretation of fast-framing camera D-alpha signals and the estimation of the CP wall particle flux. In addition, toroidal asymmetries, observed in measurements of toroidal electron density perturbations and the phase of magnetohydrodynamic modes, are present throughout the shutdown and can account for a discrepancy in the assimilation rate for up to 50% of the observed D 2 particle inventory. These sources and sinks of particles and fluxes were identified using absolutely calibrated D-alpha brightness and Langmuir probes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Physics-based adaptivity of a spectral method for the Vlasov–Poisson equations based on the asymmetrically-weighted Hermite expansion in velocity space

We propose a spectral method for the 1D-1V Vlasov–Poisson system where the discretization in velocity space is based on asymmetrically-weighted Hermite functions, dynamically adapted via a scaling α and shifting u of the velocity variable. Specifically, at each time instant an adaptivity criterion selects new values of α and u based on the numerical solution of the discrete Vlasov–Poisson system obtained at that time step. Once the new values of the Hermite parameters α and u are fixed, the Hermite expansion is updated and the discrete system is further evolved for the next time step. The procedure is applied iteratively over the desired temporal interval. The key aspects of the adaptive algorithm are: the map between approximation spaces associated with different values of the Hermite parameters that preserves total mass, momentum and energy; and the adaptivity criterion to update α and u based on physics considerations relating the Hermite parameters to the average velocity and temperature of each plasma species. For the discretization of the spatial coordinate, we rely on Fourier functions and use the implicit midpoint rule for time stepping. The resulting numerical method possesses intrinsically the property of fluid-kinetic coupling, where the low-order terms of the expansion are akin to the fluid moments of a macroscopic description of the plasma, while kinetic physics is retained by adding more spectral terms. Moreover, the scheme features conservation of total mass, momentum and energy associated in the discrete, for periodic boundary conditions. A set of numerical experiments confirms that the adaptive method outperforms the non-adaptive one in terms of accuracy and stability of the numerical solution.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Uniform high-k amorphous native oxide synthesized by oxygen plasma for top-gated transistors

The integration of high-k gate dielectrics with two-dimensional (2D) semiconducting channel materials is essential for high-performance and low-power electronics. However, the conformal deposition of a uniform high-k dielectric with sub-1-nm equivalent oxide thickness (EOT) and high interface quality on high-mobility 2D semiconductors is still challenging. Here we report a facile approach to synthesize a uniform high-k (ε r ~ 22) amorphous native oxide Bi 2 SeO x on the 2D semiconducting Bi2O2Se using O 2 plasma at room temperature. The conformal native oxide can directly serve as gate dielectrics with EOT of ~0.9 nm, while the original properties of underlying 2D Bi 2 O 2 Se is preserved. Furthermore, high-resolution area-selective oxidation of Bi 2 O 2 Se is achieved to fabricate discrete electronic components. This facile integration of a 2D semiconductor and its high-k native oxide holds high promise for next-generation nanoelectronics.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Material migration in W and Mo during bubble growth and fuzz formation

Growth of helium (He) induced bubbles and fuzz in tungsten (W) and molybdenum (Mo) is investigated using samples of W films on Mo substrates and Mo films on W substrates exposed to He-containing plasma in the temperature range of 340 to 1075 K, fluence range of 1.0–14 × 10 25 He·m -2 , and incident ion energy of <50 eV. No fuzz (only up to 2 nm diameter bubbles) and no material transport occur in W films at ≤750 K, while precursors-of or fully-developed fuzz and material mixing occur in W and Mo films at ≥800 K. This suggests that fuzz forms in multi-material systems as long as one material meets the conditions for fuzz formation, namely T s / T m ~ 0.27–0.5 where T s and T m are the sample exposure and material melting temperatures, respectively. Larger He bubbles, more material mixing, and further-developed fuzz occur at higher temperature due to increased mobility of He atoms and small He clusters. Accumulation of substrate material at the surface of fuzzy W and Mo thin-film (<80 nm) samples suggests fuzz growth by material transport from the bubble layer in the bulk up to the fiber tip, likely by a two-step process: (i) diffusion of punched dislocation loops in the bulk toward the fuzz base and (ii) diffusion of adatoms along the fuzz base and fiber surface (with effective transport of adatoms upwards due to trapping of adatoms at curved surfaces of fiber tips and/or due to the continuous generation of adatoms at the fuzz base). While the bubble size and fuzz thickness increase with reduced W concentration in Mo thin-film samples at 838 K likely due to an increase in trap mutation and dislocation loop punching in Mo compared to W, the fuzz thickness decreases with reduced W concentration at 1075 K despite an increase in the bubble size likely due to slower diffusion of interstitial loops in Mo.

Physics↗
DOE OSTI2023

Bayesian calibration of viscous anisotropic hydrodynamic simulations of heavy-ion collisions

Owing to large pressure gradients at early times, standard hydrodynamic model simulations of relativistic heavy-ion collisions do not become reliable until O(1) fm/c after the collision. To address this one often introduces a prehydrodynamic stage that models the early evolution microscopically, typically as a conformal, weakly interacting gas. In such an approach the transition from the prehydrodynamic to the hydrodynamic stage is discontinuous, introducing considerable theoretical model ambiguity. Alternatively, fluids with large anisotropic pressure gradients can be handled macroscopically using the recently developed viscous anisotropic hydrodynamics (VAH). In high-energy heavy-ion collisions VAH is applicable already at very early times, and at later times transitions smoothly into conventional second-order viscous hydrodynamics. Here we present a Bayesian calibration of the VAH model with experimental data for Pb-Pb collisions at the LHC at $\sqrt{^{S}NN}$ = 2.76 A TeV. We find that the VAH model has the unique capability of constraining the specific viscosities of the quark-gluon plasma at higher temperatures than other previously used models.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Study of Shock Formation Parameters With Drive Conditions in Magnetically Accelerated Plasma Flows

We present experimental data regarding the formation of high-energy-density shocks in magnetically accelerated plasma flows using pulsed power drivers. We quantify the flow velocity and temperature of the ablated plasma using optical Thomson scattering and gated emission imaging across two different generators. We show that, regardless of the drive parameters, the plasma flows show continuous acceleration over centimeter spatial scales, in line with trends in published simulation work. When stationary targets are placed in these supersonic flows, bow-shock formation is observed at all drive parameters in a range of materials. In the higher density flow generated on the 1-MA COBRA generator at Cornell University, heating of the upstream flow ahead of the shock is observed and quantified, which is not observed at the lower density flow on the 0.2-MA Bertha driver at UC San Diego. Here, when combined with previous work on the XP generator at Cornell, we can show that these three experimental setups allow control of the effect of radiation loss and upstream absorption on the formation of the bow shock.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Spectroscopic Characterization of Plasmoid Properties During Pellet Fueling in W7-X

This study utilizes a spectroscopic approach to investigate the properties of plasmoids that are formed during the process of cryogenic hydrogen pellet fueling in the Wendelstein 7-X (W7-X) stellarator. An analysis of the Balmer series emissions was conducted using a diagnostic that was installed during the 2024 operational campaign. Electron temperature, density, and plasma beta ( β ) values can be inferred from the emissions of radiation from the ablation plasmoid. These values are essential for validating pellet ablation models and, in the future, optimizing fueling strategies in steady-state fusion devices.

Cryogenic hydrogen pellets↗