SEARCH · Engineering Papers
Results for “plasma flow”
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.
In-situ coating of silicon-rich films on tokamak plasma-facing components with real-time Si material injection
Abstract Experiments have been conducted in the DIII-D tokamak to explore the in-situ growth of silicon-rich layers as a potential technique for real-time replenishment of surface coatings on plasma-facing components (PFCs) during steady-state long-pulse reactor operation. Silicon (Si) pellets of 1 mm diameter were injected into low- and high-confinement (L-mode and H-mode) plasma discharges with densities ranging from 3.9– 7.5 × 10 19 m −3 and input powers ranging from 5.5 to 9 MW. The small Si pellets were delivered with the impurity granule injector at frequencies ranging from 4 to 16 Hz corresponding to mass flow rates of 5–19 mg s −1 (1– 4.2 × 10 20 Si s −1 ) at cumulative amounts of up to 34 mg of Si per five-second discharge. Graphite samples were exposed to the scrape-off layer and private flux region plasmas through the divertor material evaluation system to evaluate the Si deposition on the divertor targets. The Si II emission at the sample correlates with silicon injection and suggests net surface Si-deposition in measurable amounts. Post-mortem analysis showed Si-rich coatings containing silicon oxides, of which SiO 2 is the dominant component. No evidence of SiC was found, which is attributed to low divertor surface temperatures. The in-situ and ex-situ analysis found that Si-rich coatings of at least 0.4–1.2 nm thickness have been deposited at 0.4–0.7 nm s −1 . The technique is estimated to coat a surface area of at least 0.94 m 2 on the outer divertor. These results demonstrate the potential of using real-time material injection to form Si-enriched layers on divertor PFCs during reactor operation.
Dissipative magnetohydrodynamics for nonresistive relativistic plasmas: An implicit second-order flux-conservative formulation with stiff relaxation
Based on a 14-moment closure for nonresistive (general-) relativistic viscous plasmas, here we describe a new numerical scheme that is able to handle all first-order dissipative effects (heat conduction, bulk and shear viscosities), as well the anisotropies induced by the presence of magnetic fields. The latter is parametrized in terms of a thermal gyrofrequency or, equivalently, a thermal Larmor radius and allows to correctly capture the thermal Hall effect. By solving an extended Israel-Stewart-like system for the dissipative quantities that enforces algebraic constraints via stiff-relaxation, we are able to cast all first-order dissipative terms in flux-divergence form. This allows us to apply traditional high-resolution shock capturing methods to the equations, making the system suitable for the numerical study of highly turbulent flows. We present several numerical tests to assess the robustness of our numerical scheme in flat spacetime. The 14-moment closure can seamlessly interpolate between the highly collisional limit found in neutron star mergers, and the highly anisotropic limit of relativistic Braginskii magnetohydrodynamics appropriate for weakly collisional plasmas in black-hole accretion problems. We believe that this new formulation and numerical scheme will be useful for a broad class of relativistic magnetized flows.
Novel features of attractors and transseries in nonconformal Bjorken flows
In this work we investigate the impact of conformal symmetry breaking on hydrodynamization of a far-from-equilibrium fluid. We find a new kind of transseries solutions for the nonconformal hydrodynamic equations of a longitudinal boost invariant expanding plasma. The new transseries solutions unveil a rich physical structure which arises due to the interplay of different physical scales. In the perfect fluid case the nonconformal speed of sound slows down the cooling of the temperature due to the emergence of logarithmic corrections that depends on the mass of the particle. Further, these terms propagate into the perturbative and nonperturbative sectors of the transseries once viscous corrections are included. The logarithmic mass contributions increase the asymptotic value of the Knudsen number while decreasing the damping rate of the transient nonhydrodynamic modes and thus, yielding to an extremely slow hydrodynamization process where flow lines merge to their forward attractor at extremely late times. The early time free streaming expansion is modified and receives logarithmic mass corrections induced by the shear-bulk couplings. The global flow structure and numerical analyses carried out in our work demonstrate the existence of local attraction at early and late times for the shear viscous tensor and bulk viscous pressure.
Mitigating Hydrogen-Induced Degradation of Iridium Anodes in Proton Exchange Membrane Water Electrolyzers
One promising method for reducing precious metal usage in proton exchange membrane water electrolysis is lowering iridium (Ir) loading at the anode. However, low-loading Ir catalysts often suffer from poor stability under high current densities. In this study, hydrogen (H2) crossover from the cathode to the anode is identified as a key degradation pathway. Temperature-programmed reduction confirms the reduction of IrO2 at 80 °C in a hydrogen environment, highlighting the vulnerability of IrO2-based catalysts to H2 exposure. To mitigate this effect, palladium (Pd) is introduced as an anode additive, acting as an H2 oxidation catalyst and mitigating IrO2 reduction. This protective role is verified by inductively coupled plasma optical emission spectroscopy and in situ X-ray absorption spectroscopy, showing significantly suppressed Ir dissolution at 80 °C under H2 flow when an O-covered Pd surface is present at oxygen evolution reaction potentials. Results from the current study identify a new strategy in improving activity and durability of catalysts in electrolyzers.
Theory and simulations of linear and nonlinear two-dimensional Rayleigh–Taylor dynamics with variable acceleration
Interfacial Rayleigh–Taylor mixing is crucial to describing important natural and engineering processes, such as exploding supernovae, laser micromachining, hot spots in inertial confinement fusion, and optical telecommunications. These require the characterization of the time dependence of the driving acceleration. We compare our theoretical formulation based on group theory foundations with interface-capturing numerical simulations for linear and nonlinear two-dimensional Rayleigh–Taylor instabilities in a finite-sized domain with time-varying acceleration over broad ranges of Atwood numbers and acceleration exponents. Detailed corroboration between theory and simulations is provided for this foundational case. Both demonstrate the strong interfacial nature of Rayleigh–Taylor instabilities, which suggests that practical flow fields can be reconstructed from the derived fluid potential using the proposed theory. A robust agreement is also obtained for the early and late-time evolution of the amplitudes of the bubble and spike, which demonstrate that the Rayleigh–Taylor flow can transition to the mixing regime even for a single-mode initial perturbation. Corroboration with experiments of high energy density plasmas motivated by studies of supernovae is also achieved. In addition, a long-standing puzzle in Rayleigh–Taylor dynamics on the interplay between the acceleration, the shear, and the interface morphology in the theory and simulations is resolved by accounting for finite viscosity of the fluids. The characterization of Rayleigh–Taylor instabilities as a highly interfacial phenomenon provides valuable insight into its multiscale nature, which enhances the design and understanding of numerous processes of practical interest.
Validating methods for modeling composition gradients in planar shock experiments
An interface is Rayleigh–Taylor (RT) unstable when acceleration pushes a less dense material into a more dense one, and the growth of the instability is governed partly by the Atwood number gradient. Double-shell inertial confinement fusion capsules have a foam spacer layer pushing on an inner capsule composed of a beryllium tamper and high-Z inner shell, and so have RT unstable interfaces that require benchmarking. To this end, the results of a planar shock experiment with beryllium/tungsten targets are presented. One target had the normal bilayer construction of beryllium and tungsten in two distinct layers; the second target had the beryllium grading into tungsten with a quasi-exponential profile, motivated by the potential for reduced RT growth with the gradient profile. Simulations mimic the shock profiles for both targets and match the shock velocity to within 5%. These results validate the ability of our simulations to model double-shell capsules with bilayer or graded layer Be/W inner shells, which are needed to design future experiments at the National Ignition Facility.
How does negative triangularity mitigate ITG turbulence and transport?
Improved confinement in negative triangularity (NT) experiments is attributed to reduced fluxes driven by micro-turbulence. The physical mechanism of why thermal confinement improves in NT relative to PT is unknown. This study employs gyrokinetic flux tube simulations using the GENE code with local Miller equilibrium to elucidate the physical mechanisms behind the beneficial effects of NT flux surface shapes. The focus is on collisionless ion temperature gradient (ITG) driven turbulence with adiabatic electrons. The kinetic profiles are held fixed across a scan of triangularity values, thus enabling comparisons on a level playing field. The reduced linear growth rates for NT is shown to be due to a reduced eigenmode averaged magnetic drift frequency and a wider, stronger negative local magnetic shear region about the outboard mid-plane. The nonlinear heat flux is lower for NT than that for PT, due to reduced radial correlation length and increased correlation time (τ c ) of fluctuations. These, in turn, are due to a comparatively higher level of self-generated zero-frequency E × B zonal shearing rate ω E in NT as compared to PT. Though the linear zonal potential residual is lower for NT, the nonlinearly generated E × B zonal shearing rate is higher for NT than for PT. This outcome is linked to the distinctive features of the radial wavenumber spectra of the zonal potential and the zonal shearing rate. The dimensionless parameter ω E τ c is suggested as a figure of merit. This is higher for NT than for PT. Thus, the reduced heat diffusivity for NT is linked to increased ω E τ c . Self-generated temperature corrugations (i.e. zonal temperature gradients) are much weaker than the background mean temperature gradient. Nevertheless, temperature corrugations are more pronounced in NT than in PT.
Probing nuclear structure of heavy ions at energies available at the CERN Large Hadron Collider
We perform high-statistics simulations to study the impacts of nuclear structure on the ratios of anisotropic flow observables in 208 Pb + 208 Pb and 129 Xe + 129 Xe collisions at the Large Hadron Collider. Even with 40% difference in atomic numbers between 208 Pb and 129 Xe nuclei, the ratios of anisotropic flow in the same centrality class between the two collision systems are strongly affected by the nuclear structure inputs in the initial state. The ratios of v 2 {4}/v 2 {2} in these collisions are sensitive to the nuclear skin thickness of the colliding nuclei, providing indirect constraints on the nuclei's neutron skin. In conclusion, our model predictions serve as a benchmark to compare with experimental measurements.
Observable signatures of initial state momentum anisotropies in nuclear collisions
In this work, we show that the correlation between the elliptic momentum anisotropy v 2 and the average transverse momentum [ p T ] at fixed multiplicity in small system nuclear collisions carries information on the origin of the observed momentum anisotropy. A calculation using a hybrid IP-Glasma+ music +UrQMD model that includes contributions from final state response to the initial geometry as well as initial state momentum anisotropies of the color glass condensate predicts a characteristic sign change of the correlator ρ ^ ( v 2 2 , [ p T ] ) as a function of charged particle multiplicity in p + Au and d + Au collisions at s = 200 GeV , and p + Pb collisions at s = 5.02 TeV . This sign change is absent in calculations without initial state momentum anisotropies. The model further predicts a qualitative difference between the centrality dependence of ρ ^ ( v 2 2 , [ p T ] ) in Au + Au collisions at s = 200 GeV and Pb + Pb collisions at s = 5.02 TeV , with only the latter showing a sign change in peripheral events. Predictions for O + O collisions at different collision energy show a similar behavior. Experimental observation of these distinct qualitative features of ρ ^ ( v 2 2 , [ p T ] ) in small and large systems would constitute strong evidence for the presence and importance of initial state momentum anisotropies predicted by the color glass condensate effective theory.
Transverse Momentum Differential Global Analysis of Heavy-Ion Collisions
Not Available
Quark flavor equilibration of the quark-gluon plasma
The early stage of a heavy-ion collision is marked by rapid entropy production and the transition from a gluon saturated initial condition to a plasma of quarks and gluons that evolves hydrodynamically. However, during the early times of the hydrodynamic evolution, the chemical composition of the QCD medium is still largely unknown. We present a study of quark chemical equilibration in the (Q)GP using a novel model of viscous hydrodynamic evolution in partial chemical equilibrium. Motivated by the success of gluon saturated initial condition models, we initialize the QCD medium as a completely gluon dominated state. Local quark production during the hydrodynamic phase is then simulated through the evolution of time-dependent fugacities for each independent quark flavor, with the timescales set as free parameters to compare different rates of equilibration. We present the results of complete heavy-ion collision simulations using this partial chemical equilibrium model, and show the effects on hadronic and electromagnetic observables. In particular, we show that the development of flow is sensitive to the equilibration timescale, providing an empirical way to probe the chemical equilibration of the QCD medium.
Electromagnetic Inhomogeneous-Energy-Density-Driven Instability in the UCLA Large Plasma Device (EM-IEDDI in LAPD) (Final Technical Report)
Using the Large Plasma Device (LAPD), the electromagnetic signatures associated with the inhomogeneous-energy-density-driven (IEDD) instability were targeted for documentation in a lab experiment. Exploitable features of LAPD are the extreme-length, high-beta, low-collisionality, and excellent reproducibility. The signatures of interest were mode amplitude, mode structure, frequency spectrum, phase velocity, and excitation threshold. The abrupt inhomogeneity at the interface between two coaxial plasma columns, that have disparate space potential values, were expected to generate a strong, localized, transverse, dc electric field and form a strong ExB-flow inhomogeneity that was expected to trigger broadband, intermediate-frequency waves that propagate many gyroradii beyond the interfacial source region. WVU was unable to destabilize the IEDDI mode in LAPD. Experimental limitation of three different parameters conspire to prevent IEDDI-favorable conditions even though Alfven-favorable conditions are met. Electron-neutral collision frequency is too high, Structuring of the radial profile of the radial electric field is negligible and the control of that structure is minimal, and the Magnetic field is too large for some purposes and too small for other purposes. WVU designed a new antenna and relied on externally-driving the IEDDI (instead of spontaneous destabilization). The antenna operation failed for unforeseen reasons and lack of electrician time and availability.
Turbulence characteristics and flow dynamics impacts on the H-mode transition in favourable magnetic geometry with lower power threshold
We report the L-H transition power threshold (P LH ) in favourable magnetic geometry (ion ∇B drift pointing towards X-point) is much lower than in the unfavourable magnetic geometry (ion ∇B drift pointing away from X-point) on multiple tokamaks. In a systematic experiment on DIII-D, the ion ∇B drift direction was changed continuously from the unfavourable to favourable configuration during plasma discharges. During such process, the input neutral beam power was kept constant at a value that was above P LH for favourable configuration, but lower than P LH for unfavourable configuration. Toroidal field and plasma current were also kept constant and there was little change in the edge electron density ne and electron temperature T e profiles. The density fluctuation amplitude was reduced approaching the transition, while a large increase of turbulence Reynolds stress and flow shear were simultaneously observed. The turbulence decorrelation rate was found to increase as the ion ∇B drift direction was moving towards the favourable configuration, but the flow shear also increased and exceeded the turbulence decorrelation rate. These measurements demonstrate an important correlation between turbulence and turbulence-driven flow and a lowering of PLH, provide insights into the underlying physics behind the hidden parameters and inform a more complete physics-based model of the L-H transition power threshold.
Active Aerodynamic Load Control for Wind Turbines
The goal of this project was to develop and demonstrate an advanced dielectric barrier discharge (DBD) plasma actuator technology. We set out to demonstrate the efficacy and impact of the new actuator technology as a key component of an active load control system for wind turbines. A DBD plasma actuator consists of a thin layer of dielectric material separating a pair of offset electrodes. One of those electrodes is embedded between the dielectric layer and a non-conductive substrate, while the other is exposed to air. When driven by an appropriate high voltage waveform, the device ionizes the air near the surface of the dielectric and adjacent to the exposed electrode. Collisions between the ions in the plasma and neutral air molecules result in a wall-jet – a region of induced air velocity that can be used to modify the flow around a lifting body. Use of the device near the trailing edge of a wind turbine blade, designed in such a way as to amplify the effects of the flow-modifying device, can result in large changes to the global forces experienced by the blade. Because of the fast response time of the device, it can allow the turbine to react in real time to changes in the wind associated with turbulence, wind shear, gusts, and the like, when paired with appropriate sensors and control algorithms. The goal of this project was to increase the capacity of the device to induce velocity on its surface at the levels required by large, utility-scale wind turbines. The first technical goal was to increase the induced velocity from the current industry-best of about 4 m/s to 10 m/s by modifying the electrical waveform used to drive the device and by introducing a semi-conductive surface coating to control electrostatic charge build-up on the surface. The second goal was to use the device to modify the lift on a representative airfoil in a wind tunnel, demonstrating a reduction in lift coefficient of 0.2 or better. The third goal was to produce an actuator-induced velocity of 20 m/s. Alongside our partners at the University of Texas at Dallas, we also applied modern, advanced design methods to optimize the impact of the technology on the design of wind turbines. We also had planned to install a segmented, active load control system on a test turbine to demonstrate the ability to reduce unsteady aerodynamic forces on the turbine associated with changes in the wind. The highest induced velocity achieved during this research was 11 m/s. However, practical design constraints limited the change in lift coefficient to about 0.12 for a representative airfoil in the wind tunnel at a Reynolds number of 400,000. The primary conclusion is that the plasma actuator control authority remains insufficient for practical purposes when extrapolated to Reynolds numbers over 1 million. The potential impact of the active lift control concept was evaluated through detailed simulations for three wind turbine sizes: a 3.4MW onshore turbine, a 10MW offshore turbine, and a 15MW offshore turbine. A feedback control system was designed for each turbine within two scenarios: one where the active lift control is used as a retrofit capability on the baseline design, and the other where the designers were permitted to “upscale” the turbines in order increase annual energy production. The levelized cost of energy was then evaluated for the range of turbine sizes and design configuration. It was found that the LCOE reduction associated with active lift control fell in the range of 0.7% to 7.2%, with the highest reduction associated with upscaling the 3.4MW turbine.
Electrostatic turbulence in EAST plasmas with internal transport barrier
Based on first-principles nonlinear gyrokinetic simulations, the electrostatic turbulence properties in the internal transport barrier (ITB) region of an Experimental Advanced Superconducting Tokamak discharge (#93890) are investigated. Specifically, ITBs with steep density and temperature gradients are located in the weakly negative magnetic shear region at the plasma center. In the linear stage, the growth rate and frequency of the ion temperature gradient (ITG) mode increase significantly due to resonant excitation by trapped electrons. That is, the resonance between trapped electrons and the ITG becomes strong due to the precession drift reversal of trapped electrons by the negative magnetic shear and Shafranov shift. Meanwhile, the trapped electron mode is stable in the ITB region due to only a very small fraction of electrons precessing in the direction of the electron diamagnetic drift. Nonlinear simulations show that, after considering the non-adiabatic effect of trapped electrons, the heat conductivity of ions and the turbulence intensity increase by at least a factor of 7 compared with the results only considering the adiabatic effect of electrons. The zonal charge density of trapped electrons can partially cancel that of ions, which weakens the intensity of the zonal flow, and consequently reduces the zonal flow regulation and enhances the turbulent transport.
Simulating Charm Quarks in IP-Glasma Initial Stage and Quark-Gluon Plasma: A Hybrid Approach for charm quark phenomenology
We present phenomenological findings on charm quark transport while including its energy loss in both pre-equilibrium and hydrodynamic stages of the evolution. We employed the MARTINI event generator for the production and evolution of heavy quarks in the relativistic heavy-ion collisions. The sensitivity of the heavy meson nuclear modification factor and flow coeffcient to the early stage of heavy-ion collisions and bulk medium evolution is analyzed for Pb+Pb collisions at 5.02 TeV. Our study provides insights into the interaction strength of charm quarks during the early phase and within the quark-gluon plasma.
High-Fidelity Arc-Discharge Model for Hydrogen-Plasma-Smelting-Reduction of Iron Ore
Electrification and use of renewable hydrogen is currently a necessity for decarbonizing the iron-and-steel industry. In this regard, hydrogen plasma smelting reduction (HPSR) is a novel pathway that is being explored for reduction of iron ore. HPSR provides several decarbonization merits compared to conventional blast furnaces. Firstly, the use of renewable hydrogen drastically reduces the CO2 emissions compared to the use of coke. Secondly, renewable electricity in the form of a thermal plasma for making reactive hydrogen species (radicals, ions) are more efficient at reducing iron ore compared to neutral H2. Thirdly, a molten product compatible with downstream processes is obtained from the intense heat transfer from the plasma. However, the scale-up of this technology requires fundamental exploration of hydrogen plasma dynamics and its interaction with complex solid material that include phase changing iron-ore and slag. In this work, we present a first principles continuum scale model for thermal plasmas in Ar/H2 gas mixtures typically used for HPSR. The thermal plasma governing equations for mass, momentum and energy with Lorentz force and Joule heating source terms are solved along with electromagnetic equations for electrostatic and magnetic vector potential. Our solver will be based on Pele, a suite of reacting flow solvers designed for advanced scientific computing architectures (Henry De Frahan et al., Proceedings of SIAM Parallel Processing, 13-25, 2024), and will utilize adaptive mesh generation for enhanced resolutions at locations of intense physicochemical interactions. This study will present the impact of Ar to H2 ratios on excited/dissociated hydrogen species concentrations, plasma temperature and conductivity along with the impact of outgassed species (water, metal vapor, O, OH radicals) from ore surface on gas phase chemistry. Furthermore, the heat and species flux to the surface will be quantified as a function of applied voltages in a transferred arc configuration.