Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “experimental evolution”

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 523 records · Page 29
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Beryllium-driven structural evolution at the divertor surface

Erosion of the beryllium first wall material in tokamak reactors has been shown to result in transport and deposition on the tungsten divertor. Experimental studies of beryllium implantation in tungsten indicate that mixed W–Be intermetallic deposits can form, which have lower melting temperatures than tungsten and can trap tritium at higher rates. To better understand the formation and growth rate of these intermetallics, we performed cumulative molecular dynamics (MD) simulations of both high and low energy beryllium deposition in tungsten. In both cases, a W–Be mixed material layer (MML) emerged at the surface within several nanoseconds, either through energetic implantation or a thermally-activated exchange mechanism, respectively. While some ordering of the material into intermetallics occurred, fully ordered structures did not emerge from the deposition simulations. Targeted MD simulations of the MML to further study the rate of Be diffusion and intermetallic growth rates indicate that for both cases, the gradual re-structuring of the material into an ordered intermetallic layer is beyond accessible MD time scales(≤1 μs). However, the rapid formation of the MML within nanoseconds indicates that beryllium deposition can influence other plasma species interactions at the surface and begin to alter the tungsten material properties. Therefore, beryllium deposition on the divertor surface, even in small amounts, is likely to cause significant changes in plasma-surface interactions and will need to be considered in future studies.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Heterogeneous microstructural evolution during hydrodynamic penetration of a high-velocity copper microparticle impacting copper

Microparticle hydrodynamic penetration (HDP) may be associated with the erosion regime in cold spray processing and other high-velocity impact events. Here, in an experimental approach where we can individually launch particles and study the impact sites, we explore copper microparticles impacted on copper substrates at velocities above 900 m/s where HDP begins. We lift cross-sectional lamellae from the impact sites with a focused-ion beam for further microstructural characterization using electron backscatter diffraction and scanning transmission electron microscopy. Due to the gradients of strain, strain rate, and temperature associated with HDP, heterogeneous microstructures result. The structural evolution processes observed include deformation twinning and multiple dislocation-mediated grain recrystallization mechanisms—geometric dynamic recrystallization (gDRX), discontinuous DRX (dDRX), and meta DRX (mDRX). The higher strains at the interface lead to the most significant structural changes and complex mechanisms. In contrast, there is a gradient to more conventional dislocation plasticity away from the interface (on either the particle or substrate side). Here, these microstructural observations are consistent with the deformation map for copper and extend the observations of impact-induced recrystallization across new regimes of behavior.

Cold spray process↗
DOE OSTI2025

Experimental characterization and atomistic simulation of grain boundary segregation in Mg-Y alloys

As a rare earth solute element in Mg alloys, Y has the beneficial effects of increasing both the strength and the ductility as well as weakening the crystallographic texture. To achieve a more fundamental understanding on how Y addition affects the microstructural evolution and mechanical properties, the Y segregation behavior at grain boundaries was investigated in Mg-1wt.%Y and Mg-7wt.%Y alloys at different conditions. The segregation intensity and its dependence on the grain boundary misorientation angle were experimentally characterized and computationally predicted. Strong segregation at grain boundaries was observed in both low and high Y-containing alloys. Y segregation was found to remain in alloy Mg-7Y after high-temperature annealing heat treatment at 540 °C. No direct correlation between the Y segregation intensity and the grain boundary misorientation angle could be established based on either the experimental characterization or the atomistic simulation with a spectral model. We thus conclude that grain boundary segregation of Y is independent of grain boundary misorientation angle.

Grain boundary↗
DOE OSTI2023

Predicting microstructurally sensitive fatigue‐crack path in WE43 magnesium using high‐fidelity numerical modeling and three‐dimensional experimental characterization

Abstract Microstructurally small fatigue‐crack growth in polycrystalline materials is highly three‐dimensional due to sensitivity to local microstructural features (e.g., grains). One requirement for modeling microstructurally sensitive crack propagation is establishing the criteria that govern crack evolution, including crack deflection. Here, a high‐fidelity finite‐element modeling framework is used to assess the performance and validity of various crack‐growth criteria, including slip‐based metrics (e.g., fatigue‐indicator parameters), as potential criteria for predicting three‐dimensional crack paths in polycrystalline materials. The modeling framework represents cracks as geometrically explicit discontinuities and involves voxel‐based remeshing, mesh‐gradation control, and a crystal‐plasticity constitutive model. The predictions are compared to experimental measurements of WE43 magnesium samples subject to fatigue loading, for which three‐dimensional grain structures and fatigue‐crack surfaces were measured post‐mortem using near‐field high‐energy x‐ray diffraction microscopy and x‐ray computed tomography. Findings from this work are expected to improve the predictive capabilities of simulations involving microstructurally small fatigue‐crack growth in polycrystalline materials.

Engineering↗
DOE OSTI2024

Dynamic data-driven multiscale modeling for predicting the degradation of a 316L stainless steel nuclear cladding material

Here, we have developed a long short-term memory stacked ensemble (LSTM-SE) surrogate modeling approach that can provide rapid predictions of microstructural evolution and the resultant mechanical properties of American Iron and Steel Institute (AISI) 316L series stainless steel (316LSS) fuel cladding under conditions of varying temperature and radiation dose rate. To acquire training data, we developed and implemented a kinetic Monte Carlo (KMC) model to simulate precipitation kinetics of M 23 C 6 , γ', and G phases within SS316L cladding. Experimentally reported precipitation kinetics of SS316L in literature were linked to the kinetic parameters of the simulated precipitation in our KMC model. The model was then used to simulate microstructure evolution under synthetically generated treatments of varying temperature and radiation dose rate, for periods of up to 3000 hours. Changes in volume fraction, number density, and particle size of precipitates were recorded, and particle area fractions were correlated using statistical methods to develop the surrogate model. Simultaneously, the mechanical properties of the simulated microstructures were evaluated using microstructure-based finite element method (FEM) analysis to determine the elastic modulus, yield stress, ultimate tensile strength, and elongation to failure of the aged microstructures. Using this approach, our surrogate model can predict precipitation behavior within 0.25% volume fraction and mechanical properties within 6% relative error from the values predicted by the KMC and FEM models using 50 training simulations as input. The trained recurrent neural network-based model can return estimations of precipitation kinetics and mechanical properties ~1000 times faster than the physics-based codes. This work demonstrates, as a proof of concept, that reactor material service lifetimes under variable service conditions can be predicted for a statistics-based model from a practicably obtainable dataset.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Dynamic Data-Driven Multiscale Modeling for Predicting Microstructure and Mechanical Properties of 316L Stainless Steel

We have developed a long short-term memory stacked ensemble (LSTM-SE) surrogate modeling approach that can provide rapid predictions of microstructural evolution and the resultant mechanical properties of American Iron and Steel Institute (AISI) 316L series stainless steel (316LSS) fuel cladding under conditions of varying temperature and radiation dose rate. To acquire training data, we developed and implemented a kinetic Monte Carlo (KMC) model to simulate precipitation kinetics of M 23 C 6 , γ', and G phases within SS316L cladding. Experimentally reported precipitation kinetics of SS316L in literature were linked to the kinetic parameters of the simulated precipitation in our KMC model. The model was then used to simulate microstructure evolution under synthetically generated treatments of varying temperature and radiation dose rate, for periods of up to 3000 hours. Changes in volume fraction, number density, and particle size of precipitates were recorded, and particle area fractions were correlated using statistical methods to develop the surrogate model. Simultaneously, the mechanical properties of the simulated microstructures were evaluated using microstructure-based finite element method (FEM) analysis to determine the elastic modulus, yield stress, ultimate tensile strength, and elongation to failure of the aged microstructures. Using this approach, our surrogate model can predict precipitation behavior within 0.25% volume fraction and mechanical properties within 6% relative error from the values predicted by the KMC and FEM models using 50 training simulations as input. The trained recurrent neural network-based model can return estimations of precipitation kinetics and mechanical properties ~1000 times faster than the physics-based codes. This work demonstrates, as a proof of concept, that reactor material service lifetimes under variable service conditions can be predicted for a statistics-based model from a practicably obtainable dataset.

36 MATERIALS SCIENCE↗
DOE OSTI2024

High-temperature 205 Tl decay clarifies 205 Pb dating in early Solar System

Radioactive nuclei with lifetimes on the order of millions of years can reveal the formation history of the Sun and active nucleosynthesis occurring at the time and place of its birth. Among such nuclei whose decay signatures are found in the oldest meteorites, 205 Pb is a powerful example, as it is produced exclusively by slow neutron captures (the s process), with most being synthesized in asymptotic giant branch (AGB) stars. However, making accurate abundance predictions for 205 Pb has so far been impossible because the weak decay rates of 205 Pb and 205 Tl are very uncertain at stellar temperatures. To constrain these decay rates, we measured for the first time the bound-state β - decay of fully ionized 205 Tl 81+ , an exotic decay mode that only occurs in highly charged ions. The measured half-life is 4.7 times longer than the previous theoretical estimate and our 10% experimental uncertainty has eliminated the main nuclear-physics limitation. With new, experimentally backed decay rates, we used AGB stellar models to calculate 205 Pb yields. Propagating those yields with basic galactic chemical evolution (GCE) and comparing with the 205 Pb/ 204 Pb ratio from meteorites, we determined the isolation time of solar material inside its parent molecular cloud. We find positive isolation times that are consistent with the other s-process short-lived radioactive nuclei found in the early Solar System. Our results reaffirm the site of the Sun’s birth as a long-lived, giant molecular cloud and support the use of the 205 Pb– 205 Tl decay system as a chronometer in the early Solar System.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

In Situ Scattering Studies of Crystallization Kinetics in a Phase-Separated Zr–Cu–Fe–Al Bulk Metallic Glass

In this work, thermal stability and the crystallization kinetics of a phase-separated Zr-Cu-Fe-Al bulk metallic glass were investigated using in situ high-energy synchrotron X-ray and neutron diffraction, as well as small-angle synchrotron X-ray scattering. It was revealed that this glass with excellent glass-forming ability possesses a two-step crystallization behavior. The crystalline products and their evolution sequence are more complicated than a homogeneous Zr-Cu-Al glass with average glass-forming ability. The experimental results indicate that a finely distributed nanometer-sized cubic Zr 2 Cu phase forms first and then transforms to a tetragonal Zr 2 Cu phase, while the matrix transforms to an orthorhombic Zr3Fe phase. The strength of the Zr-Cu-Fe-Al composite containing cubic Zr 2 Cu phase and glass matrix increases, and the plasticity also improves compared to the as-cast Zr-Cu-Fe-Al bulk metallic glass. Our results suggest that the formation of multiple and complex crystalline products would be the characteristics of the Zr-Cu-Fe-Al glass with better glass-forming ability. Our study may shed light on the synthesis of bulk-sized glass-nanocrystals composites of high strength and good plasticity.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Multi-objective surrogate-assisted calibration of CPFEM models using macroscopic response and in situ EBSD measurements of grain reorientation trajectories

Crystal plasticity finite element method (CPFEM) models are widely used to simulate the deformation behaviour of polycrystalline materials, but their calibration is often limited by their high computational cost and the non-convexity of the optimisation landscape. Here, this study develops a multi-objective surrogate-assisted calibration workflow that couples a multi-objective genetic algorithm (MOGA) with an adaptively trained deep neural network (DNN) surrogate model to efficiently identify CPFEM parameters from experimental data. The workflow is demonstrated on three crystal plasticity (CP) formulations of increasing complexity — Voce hardening (VH), two-coefficient latent hardening (LH2), and six-coefficient latent hardening (LH6) — using in situ electron backscatter diffraction (EBSD) measurements of Alloy 617 under uniaxial tensile loading. The CPFEM models are calibrated against the experimentally observed stress–strain response and reorientation trajectories of eight grains, then validated against eight additional trajectories and overall texture evolution. Across the CP formulations, the macroscopic response was reproduced reliably, while differences emerged in the robustness and accuracy of the grain-scale predictions. Including grain reorientation trajectories in the multi-objective calibration improved texture evolution predictions and filtered out physically inconsistent parameter sets that can arise from calibrating against only the stress–strain data. The workflow also demonstrates good transferability of calibrated parameters from a low- to a high-fidelity microstructural model. These results provide practical guidance for integrating in situ microstructural data into CPFEM through efficient, repeatable, and physically meaningful multi-objective calibration.

Crystal plasticity finite element method↗
DOE OSTI2022

Critical heat flux on zircaloy and accident tolerant fuel cladding under prototypical conditions of pressurized and boiling water reactors

First of a kind high spatial and temporal resolution temperature measurements of Pressurized Water Reactor and Boiling Water Reactor simulated fuel pins (high heat flux cosine profile heaters) in annular rod type geometry during nucleate boiling, departure from nucleate boiling and re-wetting has been performed. Distributed fiber optic measurements of the cladding axial temperature with spatial resolutions down to 2.5 mm axially at up to three equally spaced azimuthal locations per rod have demonstrated the similarities and differences of the different accident tolerant fuel cladding materials. These experiments show the performance of accident tolerant cladding with a precise measurement of the location and evolution of dryout followed by the re-wetting phenomena (see supplementary time history of a critical heat flux event). Experimental parameters varied coolant mass fluxes from 1695 to 2712 kg/m 2 s, pressures from 10 to 20 MPa, and inlet subcooling from 10 to 55 °C. It was found that the coatings and different material have little effect on the value of the critical heat flux. Any discrepancies between the local critical heat flux verified among the cladding material are attributed to different departure from nucleate boiling phenomena captured during the experiment with the high-resolution optical fiber temperature data. Substantial differences in the performance of the different cladding material to survive a critical heat flux event were however observed by post-test X-ray and optical inspection. External evident damage and internal cracks to the simulated fuel pin having the bare zircaloy cladding were observed in regions associated with the occurrence of the critical heat flux for longer periods. Such damages were associated with zirconium oxidation, and not noticed on the Cr coated zircaloy or FeCrAl claddings. These results indicate that the use of FeCrAl or a Cr coated zircaloy as cladding in Pressurized Water Reactor and Boiling Water Reactor may increase the survivability of the cladding material for short transients.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2025

Collaborative or competitive phase transformation processes in an additively manufactured-maraging steel M300

Although retained austenite is present in minor amounts in as-built additively manufactured-maraging steels, its evolution during ageing is key to controlling the final microstructure. This study investigates the relationship between austenite reversion and precipitation evolution during ageing, particularly whether these processes compete or cooperate. Using Scanning/Transmission Electron Microscopy, Atom Probe Tomography and Thermo-Calc® PRISMA simulations, the presence and evolution of retained austenite, Ni 3 Ti and Fe 7 Mo 6 intermetallic phases was characterised across different ageing conditions. Experimental results revealed the presence of non-enriched austenite in the as-built condition. Nevertheless, with increasing ageing temperature, Ni and Mo segregation became evident. High-Resolution Transmission Electron Microscopy revealed that Ni 3 Ti precipitation primarily occurred within the martensitic matrix, while Fe 7 Mo 6 nucleated preferentially at the austenite-martensite interface. PRISMA simulations indicated early and rapid precipitation in solute-enriched areas (i.e., intercellular regions), with Ni 3 Ti forming prior to Fe 7 Mo 6 . Both experimental and simulation results suggested Ni diffusion controls retained austenite growth, while Ti and Mo drive Ni 3 Ti and Fe 7 Mo 6 precipitation, respectively. The evidence reported in this work supports a both collaborative and independent phase transformation mechanism, where austenite reversion and precipitation co-evolve.

42 ENGINEERING↗
DOE OSTI2022

A level set approach for the computational study of a yield stress fluid filling a thin mold

Many important engineering and scientific applications such as cement slurries, foams, crude oil, and granular avalanches involve the concept of yield stress. Therefore, modeling yield stress fluids in different flow configurations, including the accurate prediction of the yield surface, is important. In this paper, we present a computational model based on the finite element method to study the flow of yield stress fluids in a thin mold and compare the results with data from flow visualization experiments. Here, we use the level set method to describe the interface between the filling fluid and air. We use polypropylene glycol as a model Newtonian fluid and Carbopol for the model yield stress fluid, as the Carbopol solution demonstrates yielding without thixotropy. To describe the yielding and shear-thinning behavior, we use a generalized Newtonian constitutive equation with a Bingham–Carreau–Yasuda form. We compare the results obtained from the mold filling experiments with the results from the three-dimensional (3D) model and from a reduced-order Hele-Shaw (HS) model that is two-dimensional, including the effect of shear-thinning along the thin direction only approximately. We show that both the 3D and the HS model can capture the experimental meniscus shape reasonably well for all the fluids considered at three different flow rates. This indicates that the shape evolution is insensitive to the dimensionality of the model. However, the viscosity and yield surfaces predicted by the 3D and HS models are different. The HS model underestimates the high viscosity and unyielded regions compared to the estimation by the 3D model.

42 ENGINEERING↗
DOE OSTI2021

Dislocation Loops in Proton Irradiated Uranium-Nitrogen-Oxygen System

Here in this study, we investigated the dislocation loop types formed in the proton-irradiated uranium-nitrogen-oxygen (U-N-O) system, which involves uranium mononitride (UN), uranium sesquinitride (a-U2N3), and uranium dioxide (UO2) phases. The dislocation loop formation is examined using specimens irradiated at 400°C and 710°C. Based on the detailed transmission-based electron microscopy characterization with i) the morphology-based on-zone and ii) the invisibility-criterion based two-beam condition imaging techniques, only a single type of dislocation loop in each phase is found: a/2?110?, a/2?111?, or a/3?111? dislocation loops in UN, a-U2N3, and UO2 phases, respectively. Molecular statics calculations for the formation energy of perfect and faulted dislocation loops in UN phases indicate a critical loop size of ~ 6 nm, above which perfect loops are energetically favorable. This could explain the absence of faulted loops in the experimental observation of the irradiated UN phase at two temperatures. This work will enhance the understanding of irradiation induced microstructural evolution for uranium nitrides as advanced nuclear fuels for the next-generation nuclear reactors.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2025

Tracking Dendritic Growth in Hydrogen-Based Hematite Reduction via Computer Vision

The reduction of hematite to metallic iron using hydrogen (H2) as a reducing agent presents a promising pathway for decarbonizing steel production. In this study, we employ a combination of in situ confocal scanning laser microscopy (CSLM) and advanced computer vision techniques to quantitatively analyze dendritic growth of ferrite during H2-based reduction of iron oxide at high temperatures. A workflow integrating Watershed Image Segmentation (WIS) and Lucas-Kanade Optical Flow (LKOF) is developed to extract both global and local kinetic information from time-resolved micrograph sequences. H2 reduction experiments conducted at 1400 degrees C and 1500 degrees C demonstrate a clear correlation between temperature and reduction rate, as evidenced by accuracy of fitted Johnson-Mehl-Avrami-Kolmogorov (JMAK) parameters. Optical flow analysis further elucidates the anisotropic and branched nature of dendritic growth, providing spatially resolved velocity fields that correlate well with global transformation kinetics. The proposed methodology demonstrates strong agreement with experimental measurements and literature values, offering a robust framework for automated image-based analysis to study kinetics through microstructural evolution in the reduction of iron ore, and likely other reaction-diffusion phenomena.

08 HYDROGEN↗
DOE OSTI2024

Femtosecond Core-Level Spectroscopy Reveals Involvement of Triplet States in the Gas-Phase Photodissociation of Fe(CO) 5

Excitation of iron pentacarbonyl [Fe(CO) 5 ], a prototypical photocatalyst, at 266 nm causes the sequential loss of two CO ligands in the gas phase, creating catalytically active, unsaturated iron carbonyls. Despite numerous studies, major aspects of its ultrafast photochemistry remain unresolved because the early excited-state dynamics have so far eluded spectroscopic observation. This has led to the long-held assumption that ultrafast dissociation of gas-phase Fe(CO) 5 proceeds exclusively on the singlet manifold. Herein, we present a combined experimental–theoretical study employing ultrafast extreme ultraviolet transient absorption spectroscopy near the Fe M 2,3 -edge, which features spectral evolution on 100 fs and 3 ps time scales, alongside high-level electronic structure theory, which enables characterization of the molecular geometries and electronic states involved in the ultrafast photodissociation of Fe(CO) 5 . We assign the 100 fs evolution to spectroscopic signatures associated with intertwined structural and electronic dynamics on the singlet metal-centered states during the first CO loss and the 3 ps evolution to the competing dissociation of Fe(CO) 4 along the lowest singlet and triplet surfaces to form Fe(CO) 3 . Calculations of transient spectra in both singlet and triplet states as well as spin–orbit coupling constants along key structural pathways provide evidence for intersystem crossing to the triplet ground state of Fe(CO) 4 . Finally, our work presents the first spectroscopic detection of transient excited states during ultrafast photodissociation of gas-phase Fe(CO) 5 and challenges the long-standing assumption that triplet states do not play a role in the ultrafast dynamics.

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

Implosion dynamics of triple-nozzle gas-puff z pinches on COBRA

Experiments on the 1-MA, 220-ns COBRA generator at Cornell University were conducted to provide detailed measurements of structured cylindrical gas-puff z pinches. In the experiments, a 7 cm diameter triple-nozzle gas valve assembly with concentric outer and inner annular nozzles and a central gas jet initialize the z-pinch load with various working gases, radial density profiles, and externally applied axial magnetic fields. Planar laser-induced fluorescence provides a measure of the initial neutral gas density of the load, while three-frame laser shearing interferometry and multi-frame extreme ultraviolet (XUV) cameras reveal the formation and propagation of a magneto-Rayleigh–Taylor (MRT) unstable shock layer. Implosion trajectories are compared to simple, experimentally informed models and found to be in good agreement. Differences in the structure of the accelerating plasma sheath and evolution of the MRT instability are observed for different gas species and axial magnetic field strengths, correlating with differences in pinch uniformity and x-ray emission. Here, the average instability growth is compared to linear MRT theory predictions using the instantaneous acceleration of the best-fit implosion models and characteristic instability wavelength, with the effective Atwood number and seed perturbation size as fit parameters. For high density argon center jets, ionization prior to the arrival of the imploding plasma sheath suggests a heating mechanism consistent with photoionization by XUV self-emission.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Prediction of the energetic particle redistribution by an improved critical gradient model and analysis of the transport threshold

Based on the theory of critical gradient model (CGM) and following the simulation method proposed by Waltz et al., a combination of TGLFEP and EPtran code is employed to predict the energetic particle (EP) transport induced by Alfvén eigenmodes (AEs). To be consistent with the experiment, recent improvements to the simulation method include consideration of threshold evolution and orbit loss due to finite orbit width. The revised CGM is applied to simulate two DIII-D experimental discharges (#142111 and #153071). It well reproduces the experimental profiles with multiple unstable AEs and large-scale EP transport. Discharge #142111 had previously been simulated using a nonlinear MHD-kinetic code MEGA with a transport mechanism based on stochasticity induced by overlapping AE. By comparing the simulated EP profiles, we find that the AE transport threshold is approximated by both the MEGA nonlinear stability threshold and the proposed CGM threshold (error <5% for single n and <17% for multiple n simulation). Both of them are larger than the linear stability threshold of the most unstable AE mode by a quantity of the order of the flux needed to sustain EP transport by the background turbulence. Additionally, we have also applied the improved CGM to simulate the α particle redistribution for a China Fusion Engineering Test Reactor steady state scenario. Because of the clear separation between the AE unstable region and the loss cone, only a moderate α particle loss of ~9.6% is predicted.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Disorder-induced heating as a mechanism for fast neutral gas heating in atmospheric pressure plasmas

Abstract Recent findings suggest that ions are strongly correlated in atmospheric pressure plasmas if the ionization fraction is sufficiently high ( ≳ 10 − 5 ). A consequence is that ionization causes disorder-induced heating (DIH), which triggers a significant rise in ion temperature on a picosecond timescale. This is followed by a rise in the neutral gas temperature on a longer timescale of up to nanoseconds due to ion–neutral temperature relaxation. The sequence of DIH and ion–neutral temperature relaxation suggests a new mechanism for ultrafast neutral gas heating. Previous work considered only the case of an instantaneous ionization pulse, whereas the ionization pulse extends over nanoseconds in many experiments. Here, molecular dynamics simulations are used to analyze the evolution of ion and neutral gas temperatures for a gradual ionization over several nanoseconds. The results are compared with published experimental results from a nanosecond pulsed discharge, showing good agreement with a measurement of fast neutral gas heating.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗