Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “computational studies”

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 145 records · Page 8
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Experimental and computational study of phase space dynamics in strongly coupled plasmas with steep density gradients

Understanding how plasmas thermalize when density gradients are steep remains a fundamental challenge in plasma physics, with direct implications for fusion experiments and astrophysical phenomena. Standard hydrodynamic models break down in these regimes, and kinetic theories make predictions that have never been directly tested. Here, we present the first detailed phase-space measurements of a strongly coupled plasma as it evolves from sharp density gradients to thermal equilibrium. Using laser-induced fluorescence imaging of an ultracold calcium plasma, we track the complete ion distribution function f(x,v,t). We discover that commonly used kinetic models (Bhatnagar–Gross–Krook and Lenard–Bernstein) overpredict thermalization rates, even while correctly capturing the initial counterstreaming plasma formation. Our measurements reveal that the initial ion acceleration response scales linearly with electron temperature, and that the simulations underpredict the initial ion response. In our geometry we demonstrate the formation of well-controlled counterpropagating plasma beams. This experimental platform enables precision tests of kinetic theories and opens new possibilities for studying plasma stopping power and flow-induced instabilities in strongly coupled systems.

Bergeson, Scott (ORCID:0000000231249226)↗
DOE OSTI2025

The reaction NH 2 + CH 2 O: Kinetic measurements and computational studies

The reaction of amidogen with formaldehyde is relevant to astrochemistry and couples nitrogen and carbon chemistry in flames. The first measurements of the temperature dependence of the rate constant k were made, over 409–643 K, and are summarized as k = 7.1 × 10 −12 exp(−18.2 kJ mol −1 /RT) cm 3 molecule −1 s −1 with a confidence interval of ±20%. These results address a controversy over the presence of a significant barrier in the addition/elimination pathway leading to the production of formamide between theoretical models that include the zero-point energy of all modes at the transition state and a model that excludes transitional modes. The latter would lead to a negligible barrier, which is inconsistent with the experiments. Contrary to earlier claims, CCSDT(Q)-corrected energies combined with transition state theory yield quantitative accord with the measurements.

Astrochemistry↗
DOE OSTI2023

A combined experimental and computational study on the reaction dynamics of the 1-propynyl (CH 3 CC, X 2 A 1 ) – propylene (CH 3 CHCH 2 , X 1 A') system: formation of 1,3-dimethylvinylacetylene (CH 3 CCCHCHCH 3 , X 1 A') under single collision conditions

The reaction of the 1-propynyl radical (CH 3 CC; X 2 A 1 ) with propylene (CH 3 CHCH 2 ; X 1 A') was studied in a crossed molecular beam machine at a collision energy of 37 ± 1 kJ mol –1 . Experimental data combined with high-level electronic structure (CCSD(T)-F12/cc-pVTZ-F12//ωB97X-D/6-311G(d,p)) and RRKM calculations reveal the reaction mechanism. Here, the overall barrierless and exoergic reaction involves indirect reaction dynamics and commences preferentially with addition of 1-propynyl with its radical centre to the carbon–carbon double bond at the terminal carbon atom of propylene. This work focuses on molecular mass growth process (hydrogen loss channels) although theory suggests methyl loss as a prevalent channel. In these processes, the C6H9 collision complexes either emit atomic hydrogen or undergo isomerisation followed by atomic hydrogen loss to preferentially yield the cis/trans isomers of 1,3-dimethylvinylacetylene (2-hexen-4-yne) as the primary product. Analysis of reaction dynamics of 1-propynyl and ethynyl radicals with propylene along with their fractional abundance in deep space suggests formation of methyl- and dimethyl derivatives of vinylacetylene in cold molecular clouds. Once formed they may engage in fundamental molecular mass growth processes via the barrierless Hydrogen Abstraction Vinylacetylene Addition mechanism that leads to the formation of methyl- and dimethylnaphthalenes thus providing a versatile route to methyl-substituted PAHs in interstellar medium.

1-propynyl↗
DOE OSTI2021

Computational study of Li 3 BO 3 and Li 3 BN 2 II: Stability analysis of pure phases and of model interfaces with Li anodes

Both Li 3 BO 3 and Li 3 BN 2 materials have promising properties for use in all solid-state batteries and other technologies dependent on electrolytes with significant ionic conductivity. As the second of a two-part study, the structural properties of Li 3 BO 3 and three reported phases of Li 3 BN 2 are investigated using first-principles modeling techniques. For alpha- Li 3 BN 2 , the tetragonal P4 2 / mnm structure reported in the literature is found to be unstable as evidenced by imaginary phonon modes near the M point of its Brillouin zone. Our simulations within the harmonic approximation suggest that the real α phase has the orthorhombic space group symmetry Pmmn formed with twice as many formula units and tiny adjustments of the equivalent lattice parameters and fractional coordinates. Extending the analysis of the Pmmn α- Li 3 BN 2 structure to the quasiharmonic approximation improves the agreement between the room-temperature x-ray pattern reported in the literature and the corresponding simulation results. In anticipation of the use of the monoclinic phases of Li 3 BO 3 and Li 3 BN 2 in Li ion conducting applications, chemical stability is investigated in terms of free-energy differences of possible decomposition and Li reaction processes, finding encouraging results. As further investigations of Li 3 BO 3 and β- Li 3 BN 2 as electrolyte or coating materials, particularly for use with Li metal anodes, idealized electrolyte/Li interfaces were investigated in terms of their geometric, energetic, and electronic properties. The results find the electrolyte/Li interfaces to be quite favorable, perhaps comparable to the pioneering LiPON/Li system.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Anti-Coulomb ion-ion interactions: A theoretical and computational study

The free energy of ion solvation can be decomposed into enthalpic and entropic contributions. This helps one to understand the connection between the dielectric properties and the underlying forces. We present a simple linear-response model of screened charge interactions that provides an alternative understanding of solvation barriers. Moreover, it explains the “anti-Coulomb” interactions (attraction between like-charged ions and repulsion between opposite-charged ions) observed in both simulations and experiments. We show that this is a universal behavior associated to the nonlocal response function of any dielectric or metallic systems. Published by the American Physical Society 2024

Wills, Alec (ORCID:0000000280184692)↗
DOE OSTI2021

A Computational Study of the Thermodynamic Conditions Leading to Autoignition in Nanosecond Pulsed Discharges

Nanosecond pulsed discharges have attracted the attention of engine manufacturers due to the possibility of attaining distributed ignition sites that accelerate burn rates while resulting in very little electrode erosion. Multidimensional modeling tools currently capture the electrical structure of such discharges accurately, but resolving the chemical structure remains a challenging problem owing to the disparity of time-scales in streamer propagation (nanoseconds) and ignition phenomena (microseconds). The purpose of this study is to extend multidimensional results toward resolving the chemical structure in the wake of streamers (or the afterglow) by using a batch reactor model (BRM). This can afford the use of very detailed chemical kinetic information. The full nonequilibrium nature of the electrons is taken into account, along with fast gas heating, shock wave propagation, and thermal diffusion. The results shed light on ignition phenomena brought about by such discharges.

42 ENGINEERING↗
DOE OSTI2020

A computational study of the thermodynamic conditions leading to autoignition in nanosecond pulsed discharges

Nanosecond pulsed discharges have attracted the attention of engine manufacturers due to the possibility of attaining distributed ignition sites that accelerate burn rates while resulting in very little electrode erosion. Multidimensional modeling tools currently capture the electrical structure of such discharges accurately, but resolving the chemical structure remains a challenging problem owing to the disparity of time-scales in streamer propagation (nanoseconds) and ignition phenomena (microseconds). The purpose of this study is to extend multidimensional results towards resolving the chemical structure in the wake of streamers (or the afterglow) by using a batch reactor model. This can afford the use of very detailed chemical kinetic information. The full non-equilibrium nature of the electrons is taken into account, along with fast gas heating, shock wave propagation, and thermal diffusion. The results shed light on ignition phenomena brought about by such discharges.

ignition, plasma, kinetics↗
DOE OSTI2023

Orientation-dependent etching of silicon by fluorine molecules: A quantum chemistry computational study

Anisotropic etching is a widely used process in semiconductor manufacturing, in particular, for micro- and nanoscale texturing of silicon surfaces for black silicon production. The typical process of plasma-assisted etching uses energetic ions to remove materials in the vertical direction, creating anisotropic etch profiles. Plasmaless anisotropic etching, considered here, is a less common process that does not use ions and plasma. The anisotropy is caused by the unequal etching rates of different crystal planes; the etching process, thus, proceeds in a preferred direction. In this paper, we have performed quantum chemistry modeling of gas-surface reactions involved in the etching of silicon surfaces by molecular fluorine. The results confirm that orientation-dependent etch rates are the reason for anisotropy. The modeling of F 2 dissociative chemisorption on F-terminated silicon surfaces shows that Si–Si bond breaking is slow for the Si(111) surface, while it is fast for Si(100) and Si(110) surfaces. Both Si(100) and Si(110) surfaces incorporate a larger number of fluorine atoms resulting in Si–Si bonds having a larger amount of positive charge, which lowers the reaction barrier of F 2 dissociative chemisorption, yielding a higher etch rate for Si(100) and Si(110) surfaces compared to Si(111) surfaces. Molecular dynamics modeling of the same reactions has shown that the chosen reactive bond order potential does not accurately reproduce the lower reaction barriers for F 2 dissociative chemisorption on Si(100) and Si(100) surfaces. Furthermore, reparameterization is necessary to model the anisotropic etching process that occurs at lower temperatures.

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

A computational study of turbulent separated flow over a wall-mounted cube at two different Reynolds numbers and incoming velocity profiles

Simulations of the canonical wall-mounted unit cube subjected to two distinct incident velocity profiles is performed at Reynolds numbers (Re) where the bulk flow characteristics are known to become relatively Re-insensitive using both direct numerical simulation (DNS) and wall-resolved large eddy simulation (WRLES). The aim of this work is to highlight the sensitivity of such bluff-body flows to mean shear and to provide a representative set of data for such flow scenarios where common turbulence modeling techniques often fail. Simple boundary conditions are selected to allow for easy comparisons for model development purposes and avoid confounding effects from freestream turbulence. In addition to mean velocity profiles, select turbulence statistics are presented in detail. We find that both mean velocity profiles and mean wake turbulence are Re-insensitive over the range examined here.

42 ENGINEERING↗
DOE OSTI2023

Formation of Zeolites Responsible for Waste Glass Rate Acceleration: An Experimental and Computational Study for Understanding Thermodynamic and Kinetic Processes

Before nuclear waste glasses can be disposed at geological disposal facilities, a performance assessment (PA) must be conducted to provide confidence to the regulators and public that the release rates of radioisotopes and hazardous materials into the environment are below the regulatory limits. Because the release rate of radionuclides is directly related to the dissolution rate of the radionuclide-containing glass, the corrosion behavior of the glass must be well understood for the disposal site lifetime. For simplicity, glass dissolution in static conditions has been categorized into three different stages, referred to as initial rate (Stage I), residual rate (Stage II) and, depending upon the secondary phases formed, a possible delayed resumption in the alteration rate (Stage III). The last of these stages is the least understood and has the possibility to drive radionuclide release rates above regulatory thresholds.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2023

Experimental and Computational Studies of Stress Corrosion Cracking of Alloys 308/309 and 82/182 Weldments in Corrosive and Radiation Environment

The goal of the project was to determine factors that influence SCC and IASCC in weldments found in LWR nuclear power plants. We focused on a SA508-304L SS weldment fabricated by EPRI using gas tungsten arc welding and used an aggressive BWR normal water chemistry (NWC) immersion environment. This weldment used 309L butter and 308L groove filler material. The microstructure of the 309L butter was non-uniform, exhibiting a 20–30μm thick martensitic layer closest to the SA508 interface, a 1–4 mm thick single γ austenite phase dilution zone, and a γ–δ duplex region extending to the 308L groove filler. The 308L groove filler had an entirely γ–δ duplex microstructure. Two approximately 1-inch thick 304L and SA508 plates approximately 12 by 6 square inches were joined using standard nuclear grade welding techniques. This included a post weld heat treatment of the 309L butter after application, a 0.32 cm fit-up root opening, and 57 bead lines of 308L groove filler applied in 18 layers. A 60 degree weld bevel angle was used and the 309L butter was 1.5 cm thick. Displacement cascade damage was induced using proton irradiation at the Michigan Ion Beam Laboratory. The incident proton energy was 2 MeV, the sample temperature was 360 ºC, and the calculated dpa value at 10 μm (60% of the Bragg peak depth of ~18 μm) was 5 dpa using the quick Kinchin-Pease model. Proton irradiation to this damage level required approximately 125 hours of beam time. Two types of samples were irradiated, tensile specimens and TEM bars. Samples were selected from all regions of the weldment, including the SA508-309L butter interface, the 309L-308L interface, and 308L-304L interface. The stainless steel alloys (304L, 308L, and 309L) within the heat affected zone are characterized by a duplex skeletal morphology of δ-ferrite and γ-austenite resulting from the recrystallization associated with weld fabrication. Approximately 7 to 8 mm of length along the specimen was irradiated. The gauge volume surfaces were mechanically polished and then electro-polished to remove mechanical damage from the mechanical polishing step prior to irradiation. Immersion tests were performed in a recirculating autoclave under BWR NWC conditions (2000 ppb wt. dissolved oxygen, neutral pH, 288 ºC, 10 MPa, and inlet water conductivity <100 nS/cm) to accelerate corrosion. Constant strain rate tests were performed either to failure or to approximately 5% strain. Strain rates of 10 -7 to 10 -6 mm/mm/s were used and typical immersion testing required four to six weeks to achieve failure or strains near 5%. Analysis primarily used advanced electron microscopy techniques of FIB lift out specimens.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2024

Universality of Collective Density Fluctuations in Liquids at and away from Equilibrium – An Integrated Neutron Scattering, Theoretical, and Computational Study (Final Technical Report)

Liquids, ubiquitous on earth, are prototypical disordered condensed matter. As the mediate phase between non-structured gases and ordered solids, liquids p ossess intriguing complexities that are absent on either extreme counter phase, especially when they are driven away from equilibrium. Despite decades of studies, the physics of liquids still keeps surprising us and challenges our understanding of condensed matter, emergent phenomena, and complex systems. Among the many interesting properties, the nature of collective dynamics in liquids remains elusive because of both the presence of strong interactions and the absence of translational invariance. This project studied the universal features of collective density fluctuations (or collective modes), which include both phonon-like excitations and collective relaxations, in the Q-dependent generalized hydrodynamic regime in liquids and liquid-like matter both at and away from equilibrium. We performed systematic studies of the fast collective dynamics (including both phonon-like excitations and collective relaxations) of four classes of liquids characterized by different fragilities and interactions using synergistically integrated coherent neutron scattering experiments (Inelastic Neutron Scattering (INS), Quasi-Elastic Neutron Scattering (QENS), and Neutron Spin Echo (NSE)), a ViscoElastic Hydrodynamics theory we have been developing, andMolecular Dynamics (MD) simulations and our analysis package LiquidLib. The transformative knowledge of the collective dynamics in liquids is important to understand and ultimately control the transport of mass, energy, and charge in liquids, both at and away from thermal equilibrium, which may will likely impact a wide range of disciplines.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Computational Study of Hydraulic Performance of South Carolina DOT Catch Basins Type CB 1, CB 16, and CB 17

Proper design of surface drainage of roadways is essential to minimize flooding and provide for traffic safety. Inlets collect the excess storm water from the drainage area of a roadway and discharge it to storm drains. Knowing the hydraulic efficiency of inlets, defined as the percentage of intercepted flow to the total street flow, is necessary in drainage design to determine inlet spacing such that the system can transport all or the majority of the road surface flow during rain events off of the road into the catch basins.

42 ENGINEERING↗