Engineering Papers⌕ Search

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.

At least 1,243 records · Page 69

MCTrans++: a 0-D model for centrifugal mirrors

The centrifugal mirror confinement scheme incorporates supersonic rotation of a plasma into a magnetic mirror device. This concept has been shown experimentally to drastically decrease parallel losses and increase plasma stability as compared with prior axisymmetric mirrors. MCTrans++ is a dimensionless (0-D) scoping tool which rapidly models experimental operating points in the Centrifugal Mirror Fusion Experiment (CMFX) at the University of Maryland. In the low-collisionality regime, parallel losses can be modelled analytically. A confining potential is set up that is partially ambipolar and partially centrifugal. Due to the stabilizing effects of flow shear, the perpendicular losses can be modelled as classical. Radiation losses such as bremsstrahlung and cyclotron emission are taken into account. A neutrals model is included, and, in some circumstances, charge-exchange losses are found to exceed all other loss mechanisms. We use the SUNDIALS ARKODE library to solve the underlying equations of this model; the resulting software is suitable for scanning large parameter spaces, and can also be used to model time-dependent phenomena such as a capacitive discharge. MCTrans++ has been used to verify results from prior centrifugal mirrors, create an experimental plan for CMFX and find configurations for future reactor-scale fusion devices.

Schwartz, Nick R. (ORCID:0000000296846435)↗

Transport physics dependence of Bohm speed in presheath–sheath transition

We report the ion exit flow speed at the sheath entrance is constrained by the Bohm criterion, which is used as a boundary condition for simulations that do not resolve the sheath region. Traditional Bohm criterion analysis invokes the equation of state and, thus, ignores transport physics in the sheath transition problem. An expression for the Bohm speed away from the asymptotic limit is derived from a set of anisotropic plasma transport equations. The thermal force, collisional temperature isotropization, and heat flux enter into the evaluation of the Bohm speed. By comparison with kinetic simulation results, this expression is shown to be accurate in the presheath–sheath transition region rather than a single point at the sheath entrance over a broad range of collisionality.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Velocity-space Origins of the Pressure–Strain Interaction in Multipopulation Distributions and Its Application to Magnetic Reconnection

A forefront research question is how energy evolves in weakly collisional plasmas for which departures from local thermodynamic equilibrium (LTE) are significant. The standard approach is studying the terms in the non-LTE energy evolution equation derived by taking the second moment of the Boltzmann equation, but the resultant fluid metrics do not retain information about which particles at which velocities drive energy evolution. A widely studied channel for internal energy density evolution is the pressure–strain interaction. Here, we employ the kinetic pressure–strain, a phase-space diagnostic whose velocity-space integral recovers the pressure–strain interaction to disambiguate the contributions to the pressure–strain interaction from disparate particle populations in composite phase-space densities. We develop phase-space analogs of the pressure–strain interaction decompositions to provide the phase-space origins of normal versus sheared flow. We introduce the “kinetic strain-rate” tensor, the phase-space analog of the strain-rate tensor, which we argue is needed to interpret the phase-space origins of the pressure–strain interaction. To demonstrate the utility of these quantities, we investigate them for composite electron distributions near the electron diffusion region in two-dimensional particle-in-cell simulations of antiparallel symmetric magnetic reconnection. We find that the phase-space-based diagnostics isolate the roles of distinct populations. These results contribute to a growing body of work providing new methods for quantifying phase-space energy evolution for a broad array of processes, from magnetic reconnection to collisionless shocks and turbulence, opening new pathways for answering longstanding problems of particle energization in weakly collisional plasmas.

79 ASTRONOMY AND ASTROPHYSICS↗

Ion-temperature- and density-gradient-driven instabilities and turbulence in Wendelstein 7-X close to the stability threshold

Electrostatic gyrokinetic instabilities and turbulence in the Wendelstein 7-X stellarator are studied. Particular attention is paid to the ion-temperature-gradient (ITG) instability and its character close to marginal stability [Floquet-type turbulence (Zocco et al., Phys. Rev. E, vol. 106, 2022, p. L013202) with no electron temperature gradient]. The flux tube version of the δƒ code stella (Barnes et al., J. Comput. Phys., vol. 391, 2019, pp. 365–380) is used to run linear and nonlinear gyrokinetic simulations with kinetic electrons. The nature of the dominant instability depends on the wavelength perpendicular to the magnetic field, and the results are conveniently displayed in stability diagrams that take this dependence into account. This approach highlights the presence of universal instabilities, which are less unstable but have longer wavelengths than other modes. A quasi-linear estimate of the heat flux suggests they are relevant for transport. Close to the stability threshold, the linear eigenmodes and turbulence form highly extended structures along the computational domain if the magnetic shear is small. Numerical experiments and diagnostics are undertaken to assess the resulting radial localisation of the turbulence, which affects the interaction of the latter with zonal flows. Increasing the amplitude of the magnetic shear (e.g.through current drive) has a stabilising effect on the turbulence and, thus, reduces the nonlinear energy transport.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Generation of Subauroral Longitudinally Extended Emissions Following Intensifications of the Poleward Boundary of the Substorm Bulge and Streamer Production

Observations showing the development of Subauroral Longitudinally Extended Emigerssions on a global scale are presented. It is demonstrated that they occur as a separation of an east-west arc-like band of luminosity that detaches away from the equatorward edge of the auroral distribution following episodes of auroral streamer production. They persist for time intervals on the order of ~30 min and devolve into patchy disjointed segments before they fade. Emissions are seen in both 557.7 nm OI and 391.4 nm N + 2 1NG lines, and 630.0 nm emissions are often observed equatorward. The 391.4 nm emissions are typically weaker and fade away more quickly than the 557.7 nm emissions, suggesting that an auroral precipitation source is initially present but is rapidly depleted as the forms age. All cases are associated with enhanced Subauroral Polarization Stream (SAPS) flows, and one event shows clear association with large-scale Giant Undulations (GUs) and the formation of spur-like forms that fold equatorward under the main oval. A model is proposed for the production of subauroral longitudinally extended forms in which nonlinear growth of SAPS-induced surface waves on the plasmapause results in a disruption of the boundary. It is suggested that hot plasma of plasma sheet origin becomes entrained in the plasmasphere to produce transient precipitation-associated auroral emissions that may decay into STEVE emissions and that cold dense plasma from the plasmasphere becomes entrained onto open drift paths to feed long-lived drainage plumes. This process may occur quasi periodically during intense substorms and storm-time conditions.

79 ASTRONOMY AND ASTROPHYSICS↗

Hydrodynamic-dissipation relation for characterizing flow stagnation

Hydrodynamic stagnation converts flow energy into internal energy. Here we develop a technique to directly analyze this hydrodynamic-dissipation process, which also yields a lengthscale associated with the conversion of flow energy to internal energy. We demonstrate the usefulness of this analysis for finding and comparing the hydrodynamic-stagnation dynamics of implosions theoretically, and in a test application to Z-pinch implosion data.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Three-beam hybrid fs/ps coherent anti-Stokes Raman scattering of rotation–vibration non-equilibrium

In this work, we demonstrate time-resolved rotational and vibrational temperature measurements to probe plasma non-equilibrium through a simple three-beam hybrid femtosecond/picosecond coherent anti-Stokes Raman scattering (fs/ps CARS) system. A single pump/Stokes pair is employed to generate both the pure-rotational and ro-vibrational Raman coherence. A novel phase-matching scheme, to our knowledge, is employed to spatially overlap the two CARS signals, simplifying the optical design and allowing signal selection by tuning a single spectrometer grating. Measurements were performed near the electrodes in a N 2 DC glow discharge. Both the rotational temperature and the vibrational populations up to $v$ = 8 were calculated from the separately measured single-shot pure-rotational and ro-vibrational CARS spectra. Strong rotation–vibration non-equilibrium was observed at both electrodes, with the cathode showing higher vibrational and rotational temperatures. In addition, non-Boltzmann behaviors were observed at both electrodes. This simplified approach enables measurements of rotational and vibrational temperatures with high spatial resolution in non-equilibrium flows.

42 ENGINEERING↗

Initial transient stage of pin-to-pin nanosecond repetitively pulsed discharges in air

In this work, evolution of parameters of nanosecond repetitively pulsed (NRP) discharges in pin-to-pin configuration in air was studied during the transient stage of initial 20 discharge pulses. Gas and plasma parameters in the discharge gap were measured using coherent microwave scattering, optical emission spectroscopy, and laser Rayleigh scattering for NRP discharges at repetition frequencies of 1, 10, and 100 kHz. Memory effects (when perturbations induced by the previous discharge pulse would not decay fully until the subsequent pulse) were detected for the repetition frequencies of 10 and 100 kHz. For 10 kHz NRP discharge, the discharge parameters experienced significant change after the first pulse and continued to substantially fluctuate between subsequent pulses due to rapid evolution of gas density and temperature during the 100 μs inter-pulse time caused by intense redistribution of the flow field in the gap on that time scale. For 100 kHz NRP discharge, the discharge pulse parameters reached a new steady-state at about five pulses after initiation. Further, this new steady-state was associated with well-reproducible parameters between the discharge pulses and substantial reduction in breakdown voltage, discharge pulse energy, and electron number density in comparison to the first discharge pulse. For repetition frequencies 1–100 kHz considered in this work, the memory effects can be likely attributed to the reduction in gas number density and increase in the gas temperature that cannot fully recover to ambient conditions before subsequent discharge pulses.

42 ENGINEERING↗

Construction of 3D MHD pressure drop correlation and flow characterization in the contraction region of a fusion blanket manifold

Inlet and outlet manifolds are typical components of liquid metal (LM) blanket designs of a fusion power reactor to be used to distribute the LM flow into breeding channels and collect it at the exit of the blanket. High pressure loss in the magnetohydrodynamic (MHD) flows featuring abrupt geometrical changes is one of the main feasibility issues of such designs. Recently, optimization studies were conducted to construct 3D MHD pressure drop correlations for a LM flow in an electrically insulating manifold with gradual expansion. Here, the 3D computational approach developed in that study is applied to the outlet manifold featuring gradual contraction. A systematic analysis was performed with a total number of 135 flow cases computed with COMSOL Multiphysics for Hartmann numbers 1000 < Ha < 10,000, Reynolds numbers 100 < Re < 12,000, and contraction angles 45° < θ < 75° for a fixed contraction ratio of 4. The effects of Ha, Re and θ on the flow recirculation, development length and the total pressure drop were carefully examined. A linear regression analysis was used to determine the power rule of pressure drop coefficient k related to Ha and Re, demonstrating a good match with the Ludford layer theory. Eventually, a correlation for the 3D MHD pressure drop coefficient was constructed as a function of Ha, Re and θ. Further, the results were compared against the inlet manifold. It was found that the flow in the inlet manifold exhibits larger recirculation zones. In the investigated range of Ha, Re and θ, the pressure drop coefficient k of the LM MHD flow in the gradual contraction is only slightly lower (< 8 %) than that in the gradual expansion.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Mutually magnetically insulated two-species Brillouin flow

In this work, we use the Brillouin flow analytic framework to examine the physics of Magnetically Insulated Transmission Lines (MITL). Here, we derive a model applicable to any particle species, including both positive and negative ions, in planar and cylindrical configurations. We then show how to self-consistently solve for two-species simultaneously, using magnetically insulated electrons and positive ions as an example. We require both layers to be spatially separated and magnetically insulated (mutually magnetically insulated); for a 7.5 cm gap with a 2 MV bias voltage, this condition requires magnetic fields in excess of 2.73 T. We see a close match between mutually insulated MITL performance and “superinsulated” (high degree of magnetic insulation) electron-only theory, as may be expected for these high magnetic fields. However, the presence of ions leads to several novel effects: (1) Opposite to electron-only theory, total electron currents increase rather than decrease as the degree of magnetic insulation becomes stronger. The common assumption of neglecting electrons for superinsulated MITL operation must be revisited when ions are present—we calculate up to 20× current enhancement. (2) The electron flow layer thickness increases up to double, due to ion space-charge enhancement. (3) The contributions from both ions and electrons to the MITL flow impedance are calculated. The flow impedance drops by over 50% when ions fill the gap, which can cause significant reflections at the load if not anticipated and degrade performance. Additional effects and results from the inclusion of the ion layer are discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Effects beyond ideal MHD on stability of wide and enhanced pedestal regimes in NSTX

Stability of edge-localized modes (ELMs) in spherical tokamaks is explored using the extended MHD model. Linear NIMROD simulations have been performed for three NSTX discharges 132543, 132588, and 141133, to investigate the role of resistivity, diffusivity, and shear flows on the onset of ELMs. The first discharge represents the wide pedestal regime and the later two discharges represent ELM-free enhanced pedestal H-mode. We first present the effect of toroidal rotation shear and find a flow shear destabilizing effect in these NSTX discharges. Simulations are also extended to include the two-fluid and ion gyroviscosity effects. Simulations show that the flow shear can shift the mode spectrum and alter the critical condition of ELM onset. We also uncover that ELM onset prediction in spherical tokamaks requires effects beyond MHD, in particular gyroviscosity and diamagnetic terms could stabilize Peeling-Ballooning modes consistent with the experimental observation of ELM-free regimes in NSTX. The findings give new insight into the nature of the interplay between resistivity, flows, and diamagnetic stabilization in ELM suppression and have potential applications to ELM control schemes in NSTX-U and next-generation spherical tokamaks. This study identifies the essential physical effects that must be included in future predictive and validation simulations.

NSTX↗

UO 2 + 5 vol% ZrB 2 nano composite nuclear fuels with full boron retention and enhanced oxidation resistance

The boron isotope ( 10 B) can be used as a neutron absorber in UO 2 to control the reactivity of nuclear fuel pellets, however, the boron source can react with oxygen source in UO 2 to form B 2 O 3 that vaporize readily at temperatures above 1200 °C. Unfortunately, the sintering of UO 2 fuel requires hours holding at high temperature (>1700 °C), resulting in an inevitable B loss during sintering and unpredictable B concentration in final product. It is challenging to incorporate boron through a conventional sintering method. In this work, we demonstrated that spark plasma sintering (SPS), a field assisted sintering technology, can effectively densify UO 2 + 5 vol% ZrB 2 composite fuel pellets by rapid consolidation at 1600 °C for a short duration of 5 min under an applied pressure of 40 MPa. Further, thermogravimetric analysis (TGA) measurements confirm that ZrB 2 is fully retained inside the composite fuel pellets. Inside the composite fuel pellets, nano sized ZrB 2 particles are uniformly distributed along the grain boundaries of the UO 2 matrix. The ZrB 2 particle transforms to a glassy B 2 O 3 phase covering the sample surface and grain boundaries of UO 2 matrix after a simple post-sintering annealing at 1000 °C in flowing Argon gas for 4 h. The formed glassy B 2 O 3 slows down the diffusion of oxygen ions and postpones the onset temperature for oxidation of UO 2 from 400 °C to 550 °C. This study demonstrates the capability of SPS, an advanced fuel manufacturing technique, to achieve a full retention of ZrB 2 in UO 2 oxide fuel and increase oxidation resistance through a simple post-sintering annealing. The reported work holds great engineering potential for development of advanced oxide fuel for nuclear application.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Acceleration mechanisms of energetic ion debris in laser-driven tin plasma EUV sources

Laser-driven tin plasmas are driving new-generation nanolithography as sources of extreme ultraviolet (EUV) radiation centered at 13.5 nm. A major challenge facing industrial EUV source development is predicting energetic ion debris produced during the plasma expansion that may damage the sensitive EUV channeling multilayer optics. Gaining a detailed understanding of the plasma dynamics and ion acceleration mechanisms in these sources could provide critical insights for designing debris mitigation strategies in future high-power EUV sources. Here, we develop a fully kinetic model of tin-EUV sources using one-dimensional particle-in-cell simulations to study ion debris acceleration, which will be valuable for cross-validation of radiation-hydrodynamic simulations. An inverse-bremsstrahlung heating operator is used to model the interaction of a tin target with an Nd:YAG laser, and thermal conduction is included through a Monte Carlo Coulomb collision operator. While the large-scale evolution is in reasonable agreement with analogous hydrodynamic simulations, the significant timescale for collisional equilibration between electrons and ions allows for the development of prominent two-temperature features. A collimated flow of energetic ions is produced with a spectrum that is significantly enhanced at high energies compared to fluid simulations. The dominant acceleration mechanism is found to be a large-scale electric field supported mainly by the electron pressure gradient, which is enhanced in the kinetic simulations due to the increased electron temperature. We discuss the implications of these results for future modeling of tin-EUV sources and the development of debris mitigation schemes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Energy Flow in Thin Shell Implosions and Explosions

Energy flow and balance in convergent systems beyond petapascal energy densities controls the fate of late-stage stars and the potential for controlling thermonuclear inertial fusion ignition. Timeresolved x-ray self-emission imaging combined with a Bayesian inference analysis is used to describe the energy flow and the potential information stored in the rebounding spherical shock at 0.22 petaPascal (2.2 Gbar or billions of atmospheres pressure). This analysis, together with a simple mechanical model, describes the trajectory of the shell and the time history of the pressure at the fuel-shell interface, ablation pressure, and energy partitioning including kinetic energy of the shell and internal energy of the fuel. The techniques used here provide a fully self-consistent uncertainty analysis of integrated implosion data, a thermodynamic-path independent measurement of pressure in the petaPascal range, and can be used to deduce the energy flow in a wide variety of implosion systems to petapascal energy densities.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Low-Temperature Plasma Treatment for Enhanced Recovery of Highly Valued Critical REEs from Coal (Final Technical Report)

The principle objective of this project was to develop a novel process flowsheet using low-temperature plasma treatment (LTP) integrated with conventional hydrometallurgical processes to recover rare earth elements (REEs), especially highly valued REEs (i.e. scandium and critical REEs), from coal and coal byproducts. LTP treatment was used to pretreat the solid feed material prior to leaching in an effort to improve the leaching characteristics of the REEs and potentially reduce acid consumption. Initially, a detailed experimental program was performed on a laboratory LTP furnace to identify a set of process parameter values to be used for assessing the enhanced leaching performance achieved on four specific gravity fractions of West Kentucky No. 13 and Fire Clay coarse coal refuse materials. High-temperature oxidation (HTO)) was also evaluated as a means for improving REE leachability on the same fractions of both coal sources. The lixiviants used to assess leachability improvements included i) untreated water, ii) 0.1 M ammonium sulfate solution and iii) 1.2 M sulfuric acid solution at a temperature of 75°C. The laboratory data was used to develop multiple process flow sheets using HTO only to produce a rare earth oxide mix concentrate or utilizing both LTP and HTO treatment processes to generate a high purity heavy REO mix concentrate and a light REO mix concentrate. A techno-economic analysis was performed on each circuit.

01 COAL, LIGNITE, AND PEAT↗

Visualization of Fast Ion Phase-Space Flow Driven by Alfvén Instabilities

Fast ion phase-space flow, driven by Alfven eigenmodes (AEs), is measured by an imaging neutral particle analyzer in the DIII-D tokamak. The flow firstly appears near the minimum safety factor at the injection energy of neutral beams, and then moves radially inward and outward by gaining and losing energy, respectively. The flow trajectories in phase space align well with the intersection lines of the constant magnetic moment surfaces and constant E – (ω/n)P ζ surfaces, where E, P ζ are energy and toroidal canonical momentum of ions; ω and n are angular frequencies and toroidal mode numbers of AEs. It is found that the flow is so destructive that the thermalization of fast ions is no longer observed in regions of strong interaction. Here, the measured phase-space flow is consistent with nonlinear hybrid kinetic-magnetohydrodynamics simulation. Calculations of the relatively narrow phase-space islands reveal that fast ions must transition between different flow trajectories to experience large-scale phase-space transport.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Quark Mass Dependence of Heavy Quark Diffusion Coefficient from Lattice QCD

We present the first study of the quark mass dependence of the heavy quark momentum and spatial diffusion coefficients using lattice QCD with light dynamical quarks corresponding to a pion mass of 320 MeV. We find that, for the temperature range 195 MeV < 𝑇 < 293 MeV, the spatial diffusion coefficients of the charm and bottom quarks are smaller than those obtained in phenomenological models that describe the 𝑝 𝑇 spectra and elliptic flow of open heavy flavor hadrons.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗