SEARCH · Engineering Papers
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.
A recent publication presents the theoretical basis for extending the Newcomb criterion for the ideal magnetohydrodynamic stability of a straight, circular cylindrical plasma to an axisymmetric toroidal plasma. Here, this paper extends the method to a nonaxisymmetric plasma with stellarator symmetry and toroidal periodicity.
There are now more Particle-in-Cell (PIC) codes than ever before that researchers use to simulate intense laser-plasma interactions. To date, there have been relatively few direct comparisons of these codes in the literature, especially for relativistic intensity lasers interacting with thin overdense targets. To address this we perform a code comparison of three PIC codes: EPOCH, LSP, and WarpX for the problem of laser-driven ion acceleration in a 2D(3v) geometry for a 10 20 W cm -2 intensity laser. We examine the plasma density, ion energy spectra, and laser-plasma coupling of the three codes and find strong agreement. We also run the same simulation 20 times with different random seeds to explore statistical fluctuations of the outputs. We then compare the execution times and memory usage of the codes (without “tuning” to improve performance) using between 1 and 48 processors on one node. We provide input files to encourage larger and more frequent code comparisons in this field.
We performed simultaneous measurement of absorption, emission, and laser-induced fluorescence spectroscopic signatures for determining nanosecond and femtosecond laser-produced plasma’s physical properties throughout its lifecycle. Plasmas are produced by focusing either ~6 ns, 1064 nm pulses from an Nd:YAG or ~35 fs, ~800 nm pulses from a Ti:Sapphire laser on an Inconel target that contains Al as a minor alloying addition. A continuous-wave narrowband tunable laser was used for performing absorption and fluorescence spectroscopy while a fast-gated detection system was used for emission spectroscopy. The temporal evolution of emission, fluorescence, and absorbance of Al transitions are compared for both ns and fs laser-produced plasms. Time-resolved absorbance was also used for evaluating linewidth, lineshape, temperature, and column-averaged atomic number density at late times of ns and fs plasma evolution. Our results demonstrate that lower and excited-state populations of fs laser produced plasmas are short-lived in comparison to those in ns plasmas. The lower state population is observed to reach a maximum value earlier in time for the fs plasma versus the ns plasma, while the kinetic temperature for the ns plasma was higher than for the fs plasma at most times of the plasma evolution.
The properties and acceleration mechanisms of electrons (<200 keV) associated with a pair of tailward traveling flux ropes and accompanied reconnection X-lines in Earth's plasma sheet are investigated with MMS measurements. Energetic electrons are enhanced on both boundaries and core of the flux ropes. The power-law spectra of energetic electrons near the X-lines and in flux ropes are harder than those on flux rope boundaries. Theoretical calculations show that the highest energy of adiabatic electrons is a few keV around the X-lines, tens of keV immediately downstream of the X-lines, hundreds of keV on the flux rope boundaries, and a few MeV in the flux rope cores. The X-lines cause strong energy dissipation, which may generate the energetic electron beams around them. The enhanced electron parallel temperature can be caused by the curvature-driven Fermi acceleration and the parallel electric potential. Betatron acceleration due to the magnetic field compression is strong on flux rope boundaries, which enhances energetic electrons in the perpendicular direction. Electrons can be trapped between the flux rope pair due to mirror force and parallel electric potential. Electrostatic structures in the flux rope cores correspond to potential drops up to half of the electron temperature. The energetic electrons and the electron distribution functions in the flux rope cores are suggested to be transported from other dawn-dusk directions, which is a 3-dimensional effect. The acceleration and deceleration of the Betatron and Fermi processes appear alternately indicating that the magnetic field and plasma are turbulent around the flux ropes.
We investigate how growth of an oxide film will influence the deuterium recycling properties of Li on plasma facing surfaces. Lithium films on the walls or plasma-facing material surfaces of a fusion vacuum vessel improves plasma performance in part by removing residual impurity atoms from the plasma. Oxygen atoms, from residual water vapor or eroded oxide surfaces, are a prominent ubiquitous impurity that can create a layer of oxide on lithium films. The results of the work herein also indicate a possibility for determining the thickness of the oxide layer on such Li films by measuring the reflection probability of incident deuterium particles from the surface.
Forward modeling is used to interpret inversion patterns of the pedestal–Scrape of Layer (SOL) Electron Cyclotron Emission (ECE) in DIII-D H-mode experiments. The modeling not only significantly improves the ECE data interpretation quality but also leads to the potential measurements of (1) the magnetic field strength |B| at the separatrix, (2) the pedestal |B| evolution during an inter-Edge Localized Mode (ELM) period, and (3) the pedestal Magnetohydrodynamics (MHD) radial structure. The ECE shine-through effect leads to three types of pedestal–SOL radiation inversions that are discussed in this paper. The first type of inversion is the non-monotonic T e,rad profile with respect to the major radius. Using the ECE frequency at the minimum T e,rad , the inversion can be applied to measure the magnetic field |B| at the separatrix and calibrate the mapping of the ECE channels with respect to the separatrix. The second type of inversion refers to the opposite phase between the radiation fluctuations δT e,rad at the pedestal and SOL. This δT e,rad phase inversion is sensitive to density and temperature fluctuations at the pedestal foot and, thus, can be used to qualitatively measure the MHD radial structure. The third type of inversion appears when the pedestal and SOL T e,rad evolve in an opposite trend, which can be used to infer the pedestal |B| field change during an inter-ELM period. Finally, the bandwidth effect on measuring δT e,rad due to pedestal MHD is also investigated in the radiation modeling.
Experimental and modeling investigations of the high βp scenarios on the DIII-D and EAST tokamaks show advantages in high energy confinement, avoidance of n = 1 MHD, and core-edge integration with reduced heat flux, making this scenario an attractive option for CFETR steady-state operation. Experiments show that plasmas with high confinement and high density can be achieved with neutral beam injection on DIII-D (β p ~ 2.2, β N ~ 3.5, f BS ~ 50%, f Gw ~ 1.0, H 98y2 ~ 1.5) and pure RF power on EAST (β P ~ 2.0, β N ~ 1.6, f BS ~ 50%, f Gw ~ 0.8, H 98y2 > 1.3). By tailoring the current density profile, a q-profile with local (off-axis) negative shear is achieved that yields improved confinement and MHD stability. Transport analysis and simulation suggest that the combination of high density gradient and high Shafranov shift allows turbulence stabilization and higher confinement. Using on-axis ECH injection, tungsten accumulation is avoided on EAST, and this is reproduced in modeling. Reduced heat flux (by >40%) and maintenance of high core confinement is achieved with active feedback control of the radiated divertor, an important result for long pulse operation in tokamaks. In conclusion, the improved physics understanding and validated modeling tools are used to design a 1GW steady-state scenario for CFETR.
Zero-dimensional kinetic modeling of atmospheric pressure Ar–N 2 –H 2 nonthermal plasma was carried out to gain mechanistic insights into plasma-assisted catalytic synthesis of ammonia. Ar dilution is a common technique for tailoring plasma discharge properties and has been shown to enhance NH 3 formation when added to N 2 –H 2 plasma. The kinetic model was developed for a coaxial dielectric barrier discharge quartz wool-packed bed reactor operating at near room temperature using a kHz-frequency plasma source. With 30% Ar mixed in a 1:1 N 2 –H 2 plasma at 760 Torr, we find that NH 3 production is dominated by Eley–Rideal (E-R) surface reactions, which heavily involve surface NH x species derived from N and H radicals in the gas phase, while the influence of excited N 2 molecules is negligible. This is contrary to the commonly proposed mechanism that excited N 2 molecules created by Penning excitation of N 2 by Ar(4s) and Ar(4p) play a significant role in assisting NH 3 formation. Our model shows that the enhanced NH 3 formation upon Ar dilution is unlikely due to the interactions between Ar and H species, as excited Ar atoms have a weak effect on H radical formation through H 2 dissociation compared to electrons. We find that excited Ar atoms contribute to 28% of the N radical production in the gas phase via N 2 dissociation, while the rest are dominated by electron-impact dissociation. Furthermore, Ar species play a negligible role in the product NH 3 dissociation. N 2 conversion sensitivity analyses were carried out for electron number density (ne) and reduced electric field (E/N), and contributions from Ar to gas-phase N radical production were quantified. Finally, the model can provide guidance on potential reasons for observing enhanced NH 3 formation upon Ar dilution in N 2 –H 2 plasma beyond changes in the discharge characteristics.
In the presence of the substrate, various structural and dynamical properties of two-dimensional dusty plasma (2DDP) were investigated using Langevin dynamical simulations. This paper reviews a series of results of the structural and dynamical properties of 2DDP modified by one-dimensional periodic substrates (1DPSs) as follows. First, when the depth of the 1DPS increases gradually, it is found that the static structural order of 2DDP along each potential well of the 1DPS increases first, and then decreases gradually. When the width of the 1DPS increases gradually, the particle arrangement in each potential well of the 1DPS changes from one straight row gradually to two rows, including the stable zigzag structure. Second, when there are two rows in each potential well of the 1DPS, the phonon spectra are split into two branches, corresponding to the breathing and sloshing modes, respectively. In addition, due to the pure repulsion between dust particles, from the obtained phonon spectra, the sloshing wave propagates backward at small wave numbers. Third, the calculated mean-square displacement (MSD) shown that, at the intermediate timescale between the initial ballistic and final diffusive motions, the particle motion tends to be more subdiffusive while the depth of the 1DPS increases. While the width of the 1DPS increases gradually, the long-time diffusive motion first increases, then decreases, and finally increases again, exhibiting the oscillation-like diffusion, due to the stable zigzag structure. Finally, when an external direct-current driving force is applied on all particles of 2DDP, three different depinning phases are discovered, which are the pinned, disordered plastic flow, and moving ordered states, respectively, as the driving force increases from zero. In addition, the continuous/discontinuous property of the phase transition between these different depinning states is investigated, showing that the transition property is modulated by the depth of the applied 1DPS
The indirect drive approach to inertial confinement fusion has undergone important advances in the past few years. Improvements in temperature and density diagnostic methods are leading to more accurate measurements of the plasma conditions inside the Hohlraum and therefore to more efficient experimental designs. The implementation of dot spectroscopy has proven to be a versatile approach to extracting space- and time-dependent electron temperatures. In this method, a microdot of a mid-Z material is placed inside the Hohlraum and its K-shell emission spectrum is used to determine the plasma temperature. However, radiation transport of optically thick lines acting within the cylindrical dot geometry influences the outgoing spectral distribution in a manner that depends on the viewing angle. This angular dependence has recently been studied in the high energy density regime at the OMEGA laser facility, which allowed us to design and benchmark appropriate radiative transfer models that can replicate these geometric effects. By combining these models with the measurements from the dot spectroscopy experiments at the National Ignition Facility, we demonstrate here a novel technique that exploits the transport effects to obtain time-resolved measurements of the ion density of the tracer dots, without the need for additional diagnostics. Furthermore, we find excellent agreement between experiment and simulation, opening the possibility of using these geometric effects as a density diagnostic in future experiments.
Collisionality is one of the key parameters in determining turbulent transport in the plasma edge, regulating phenomena such as 'shoulder formation', separation of scale lengths in the scrape-off layer (SOL), turbulence damping and zonal flow dynamics. Understanding its role is therefore of primary importance for future reactors like ITER. Getting reliable predictions and a better characterization of plasma flow properties when varying collisionality remains, however, a critical challenge for the simulations. This paper focuses on the impact of varying collisionality in a non-isothermal three-dimensional fluid model of the plasma edge. A high field side limited configuration encompassing open and closed magnetic field lines with parameters typical of a medium-sized tokamak is considered. The present model can consistently account for the variations of collisionality and its impact on both the parallel resistivity η ∥ and the ion and electron parallel thermal conductivities χ ∥e,i . Details on mean flow and turbulence properties are given. Changing collisionality leads to significant changes in the flow properties both on the mean and fluctuating quantities. In particular, lowering collisionality decreases the size of coherent structures, the fluctuation levels of turbulence, and steepens the density and temperature equilibrium profiles around the separatrix leading to a global reduction of the turbulent transport. The SOL width is observed to increase with collisionality, eventually resulting in the disappearance of the scale lengths separation between near and far SOL, consistently with previous experimental observations. At low collisionality, where the presence of narrow feature is well-established, a contribution of heat conduction increases up to compete with heat convection.
Engineering features are known to cause jets of ablator materials to enter the fuel hot-spot in inertial confinement fusion implosions. The Biermann battery mechanism wraps them in a self-generated magnetic field. Here we show that higher-Z jets have an additional thermoelectric magnetic source term that is not present for hydrogen jets, verified here through a kinetic simulation. It has similar magnitude to the Biermann term. We then include this in an extended magneto-hydrodynamic approach to post process an xRAGE radiation-hydrodynamic implosion simulation. The simulation includes an accurate model for the capsule fill tube, producing a dense carbon jet that becomes wrapped in a 4000 T magnetic field. A simple spherical carbon mix model shows that this insulates the electron heat conduction enough to cause contraction of the jet to an optically thick equilibrium. The denser magnetized jet hydrodynamics could change its core penetration and therefore the final mix mass, which is known to be well correlated with fusion yield degradation. Fully exploring this will require self-consistent magneto-hydrodynamic simulations. Experimental signatures of this self-magnetization may emerge in the high energy neutron spectrum.
Motivated by the need for accurate, timely, and efficient calculations of plasma transport, predictions of plasma turbulence properties made using different TGLF saturation rules are benchmarked against corresponding predictions from linear and nonlinear gyrokinetic CGYRO simulations. This benchmarking is carried out using parameters taken from an inductive burning plasma scenario in a hypothetical compact high-field (R maj = 4 m, B T = 8 T) tokamak, lying in a much different regime of parameter space than either the TGLF calibration regime or current-day experiments. The core turbulent transport in this scenario is predicted to be dominated by ion temperature gradient (ITG) turbulence. In general, the ITG critical gradients predicted by various TGLF saturation rules are quite close to the CGYRO predictions. Both codes predict similar linear ITG growth rates and frequency spectra, as well as their scaling with R/L T i = −Rd ln(T i )/dr. However, TGLF systematically predicts unstable trapped-electron modes (TEMs) above k y ρ s ≃ 0.5 not seen by CGYRO for the same parameters, due to TGLF predicting a lower threshold in R/L T e than CGYRO for TEM onset. It is shown that for this scenario, nonlinear CGYRO simulations predict stiffer ITG turbulence than the TGLF SAT0 and SAT1 saturation rules, with energy fluxes close in magnitude and scaling with R/L T i to what is predicted by the SAT2 saturation rule. Self-consistent core profiles calculated using nonlinear CGYRO flux predictions and the PORTALS transport solver are shown to agree fairly well with corresponding predictions made using the TGLF SAT2 model, including a similar level of density peaking.
Abstract Similarity laws (SL) are essential for correlating the characteristics of plasmas at different dimensional scales, which have been validated for radio frequency (rf) discharges at low pressure but under limited conditions. In this work, we evaluate the effects of nonlinear collisions (e.g. stepwise ionization) on the similarity properties in capacitive rf discharges in argon across a wide range of pressure regimes via fully kinetic particle-in-cell simulations. The SL scalings of fundamental discharge parameters, e.g. the electron power absorption and electron energy probability function, are examined with and without nonlinear collisions, respectively. Without the nonlinear collisions, the similarity scalings are found to be rigorously valid. When the nonlinear collisions are considered, the similarity properties in rf discharges still exist approximately, which indicates that the violations caused by the nonlinear collisions are still minor in the studied cases. The reasons for the effectiveness of SL scalings with nonlinear collision processes are also discussed. The results from this study confirmed the validity of similarity transformations with more complex reaction kinetics, which may promote the applicability of the SL scalings for the design and fabrications of plasma devices.
We report here on the characterization of the conditions of an imploding cylindrical plasma by time-resolved x-ray emission spectroscopy. Knowledge about this implosion platform can be applied to studies of particle transport for inertial confinement fusion schemes or to astrophysical plasmas. A cylindrical Cl-doped CH foam within a tube of solid CH was irradiated by 36 beams (I total ~5×10 14 W/cm 2 , 1.5ns square pulse, and E total ~ 16.2kJ) of the OMEGA-60 laser to radially compress the CH toward the axis. The analysis of the time-resolved spectra showed that the compression can be described by four distinct phases, each presenting different plasma conditions. First the ablation of the cylinder is dominant; second, the foam is heated and induces a significant jump in emission intensities; third, the temperature and density of the foam reaches a maximum; and finally, the plasma expands. Ranges for the plasma temperature were inferred with the atomic physics code SCRAM (Spectroscopic Collisional-Radiative Atomic Model) and the experimental data have been compared to hydrodynamic simulations performed with the 2D code FLASH, which showed a similar implosion dynamic over time.
Ion Cyclotron Range of Frequencies (ICRF) heating is a critical system for the SPARC mission of demonstrating net positive fusion gain. In the first campaign, 20 MW of installed power at 120 MHz will be supplied to 10 four-strap antennas, arranged in poloidal pairs. The matching network is designed to accommodate for a range of loading conditions and a fast controller adjusts the matching during a pulse with frequency modulation of ±1 MHz. For the main L-mode scenario of the first experimental campaign, the power coupling and absorption is analysed in this work. Several feeding schemes are compared and limits on the amount of the total coupled power are evaluated. Coupling of the full 20 MW power is possible well within the limit of the peak electric field in the transmission line of 15 kV/cm. The 3 He minority and 2 nd harmonic T heating scheme is validated in full-wave modelling to be efficient for heating plasma in the SPARC L-mode Q>1 scenario, with dominant ion heating and efficient single-pass absorption. A range of ICRF wave parameters and minority concentration would be suitable for operation and allow tailoring the heating properties for plasma conditions and experimental objectives.
In this project, we have further developed GTC gyrokinetic particle simulations to study plasma confinement properties in the stellarators and tokamaks with 3D magnetic fields, where both neoclassical and turbulent transport are self-consistently treated on an equal footing.
Tungsten (W) and tungsten alloys are being considered as leading candidates for structural and functional materials in future fusion energy devices. The most attractive properties of tungsten for the design of magnetic and inertial fusion energy reactors are its high melting point, high thermal conductivity, low sputtering yield, and low long-term disposal radioactive footprint. Despite these relevant features, there is a lack of understanding of how the structural and mechanical properties of W-based alloys are affected by the temperature in fusion power plants. In this work, we present a study on the thermo-mechanical properties of five W-based plasma-facing materials. First-principles density functional theory (DFT) calculations are combined with the quasi-harmonic approximation (QHA) theory to investigate the electronic, structural, mechanical, and thermal properties of these W-based alloys as a function of temperature. The coefficient of thermal expansion, temperature-dependent elastic constants, and several elastic parameters, including bulk and Young’s modulus, are calculated. Our work advances the understanding of the structural and thermo-mechanical behavior of W-based materials, thus providing insights into the design and selection of candidate plasma-facing materials in fusion energy devices.