Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “phase kinetics”

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,279 records · Page 71
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

CFD modeling of pre-spark heat release in a boosted direct-injection spark-ignition engine

Accurate predictions of low-temperature heat release (LTHR) are critical for modeling auto-ignition processes in internal combustion engines. While LTHR is typically obscured by deflagration, extremely late ignition phasing can lead to LTHR prior to the spark, a behavior known as pre-spark heat release (PSHR). In this research, PSHR in a boosted direct-injection spark-ignition engine was studied using 3-D computational fluid dynamics (CFD) and detailed chemical kinetics. The turbulent combustion was modeled via a hybrid approach that incorporates the G-equation model for tracking the turbulent flame front, and the well-stirred reactor model with detailed chemistry for assessing the low-temperature reactions in unburnt gas. Simulations were conducted using Co-Optima alkylate and E30 fuels at operating conditions characterized by different PSHR intensities. The predicted in-cylinder pressure and heat release rate were found to agree well with experiments. It was found the estimate of previous-cycle trapped residuals is of utmost importance for capturing PSHR correctly. A simulation best practice was developed which keeps the detailed chemistry solver active throughout the entire simulation, allowing to track the evolution of intermediate species from one cycle to the next. Following the validation, the dynamics of PSHR were analyzed in detail employing the pressure-temperature (P-T) trajectory framework. It was shown that PSHR correlated with the first-stage ignition delay of the fuel, hence showing close relation to the in-cylinder P-T trajectory and the chemical kinetics. Besides, it was indicated that LTHR is a self-limiting process that has the effect of attenuating the thermal stratification in the combustion chamber. Furthermore, it was observed the occurrence of PSHR caused the P-T trajectory of end-gas to overlap with the negative temperature coefficient region of the fuel’s ignition-delay maps. This effect was more significant in the fuel-rich regions where engine knock tendency would be generally higher, with potential implications on knock control and mitigation.

42 ENGINEERING↗
DOE OSTI2022

A physics-based model of cladding wastage layer formation rate

In this report, lower length scale simulations to inform an engineering-scale model of fuel-cladding chemical interaction (FCCI) conducted under the auspices of the Nuclear Energy Advanced Modeling and Simulation (NEAMS) program in FY22 are described. An overall strategy for the implementation of the BISON model is described, and the past and current lower length scale work on FCCI formation is put into the context of this overall strategy. Density functional theory and kinetic Monte Carlo simulations are used to determine the diffusion coefficient of neodymium in iron. Density functional theory calculations are also used to parameterize a parabolic approximation for the dependence of free energy on composition for the intermetallic compound Fe 17 Nd 2 . A phase-field model of the Fe-Nd system was developed that includes the random solid solution body-centered cubic phase and the intermetallic Fe 17 Nd 2 . The model was used to simulate growth of the intermetallic layer in a diffusion couple and the results were compared to experiment. The model’s prediction of the parabolic growth constant is within 40 % of the experimental value. The model is also used to simulate the formation of a grain boundary attack layer that is observed in experiment. The simulations show enhanced Nd concentration along the grain boundaries, but do not clearly show formation of a Fe 17 Nd 2 layer

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2021

Phase transitions in Schloegl's second model for autocatalysis on a Bethe lattice

Schloegl's second model (also known as the quadratic contact process) on a lattice involves spontaneous particle annihilation at rate p and autocatalytic particle creation at empty sites with n ≥ 2 occupied neighbors. The particle creation rate for exactly n occupied neighbors is selected here as n(n - 1)/[z(z - 1)] for lattice coordination number z. We analyze this model on a Bethe lattice. Precise behavior for stochastic models on regular periodic infinite lattices is usually surmised from kinetic Monte Carlo simulation on a finite lattice with periodic boundary conditions. However, the persistence of boundary effects for a Bethe lattice complicates this process, e.g., by inducing spatially heterogenous states. This motivates the exploration of various boundary conditions and unconventional simulation ensembles on the Bethe lattice to predict behavior for infinite size. Here, we focus on z = 3, and predict a discontinuous transition to the vacuum state on the infinite lattice when p exceeds a threshold value of around 0.053.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Product distribution, kinetics, and aerosol formation from the OH oxidation of dimethyl sulfide under different RO 2 regimes

Abstract. The atmospheric oxidation of dimethyl sulfide (DMS) represents a major natural source of atmospheric sulfate aerosols. However, there remain large uncertainties in our understanding of the underlying chemistry that governs the product distribution and sulfate yield from DMS oxidation. Here, chamber experiments were conducted to simulate gas-phase OH-initiated oxidation of DMS under a range of reaction conditions. Most importantly, the bimolecular lifetime (τbi) of the peroxy radical CH3SCH2OO was varied over several orders of magnitude, enabling the examination of the role of peroxy radical isomerization reactions on product formation. An array of analytical instruments was used to measure nearly all sulfur-containing species in the reaction mixture, and results were compared with a near-explicit chemical mechanism. When relative humidity was low, “sulfur closure” was achieved under both high-NO (τbi<0.1 s) and low-NO (τbi>10 s) conditions, though product distributions were substantially different in the two cases. Under high-NO conditions, approximately half the product sulfur was in the particle phase, as methane sulfonic acid (MSA) and sulfate, with most of the remainder as SO2 (which in the atmosphere would eventually oxidize to sulfate or be lost to deposition). Under low-NO conditions, hydroperoxymethyl thioformate (HPMTF, HOOCH2SCHO), formed from CH3SCH2OO isomerization, dominates the sulfur budget over the course of the experiment, suppressing or delaying the formation of SO2 and particulate matter. The isomerization rate constant of CH3SCH2OO at 295 K is found to be 0.13±0.03 s−1, in broad agreement with other recent laboratory measurements. The rate constants for the OH oxidation of key first-generation oxidation products (HPMTF and methyl thioformate, MTF) were also determined (kOH+HPMTF=2.1×10-11 cm3 molec.−1 s−1, kOH+MTF=1.35×10-11 cm3 molec.−1 s−1). Product measurements agree reasonably well with mechanistic predictions in terms of total sulfur distribution and concentrations of most individual species, though the mechanism overpredicts sulfate and underpredicts MSA under high-NO conditions. Lastly, results from high-relative-humidity conditions suggest efficient heterogenous loss of at least some gas-phase products.

Ye, Qing (ORCID:0000000337978988)↗
DOE OSTI2020

Precipitation-site competition in duplex stainless steels: Cu clusters vs spinodal decomposition interfaces as nucleation sites during thermal aging

Competing microstructural evolution mechanisms can exist simultaneously when duplex stainless steels are operating for several decades in a high temperature service environment. Such competition between different microstructural evolution pathways can be difficult to ascertain using simple model alloy systems necessitating detailed structural and compositional analysis of phase transformation mechanisms in complex alloys Thus, in this study duplex stainless steels with complex but well understood chemistries were used to investigate the relative importance of different heterogeneous nucleation sites – specifically spinodal decomposition and Cu clustering – on precipitation of Ni-Si-Mn particles during long term thermal aging. Precipitation of Ni-Si-Mn particles in duplex stainless steels and ferritic steels during thermal aging and irradiation is known to have profound effect on mechanical properties. Using duplex stainless steels with custom modified compositions along with atom probe tomography characterization and first-passage kinetic Monte Carlo (FPKMC) simulations it is revealed that while the interface between Cr and Fe formed during spinodal decomposition can be a pathway for diffusion of Ni, Si, Mn, and Cu, it is not a site that encourages precipitation of Ni-Si-Mn rich particles. Instead, the presence of a higher concentration of Cu in these steels leads to formation of small Cu clusters with high energy interfaces that act as nucleation sites for Ni-Si-Mn particles because of its strong interaction particularly with Ni increasing the chemical driving force for precipitation. These results will help inform predictive models for the use of duplex stainless steels and other precipitation-hardened alloys for extended operation at high temperatures.

Duplex stainless steel, Thermal aging, Precipitati↗
DOE OSTI2022

Electrical sensing, tracking, and actuation of droplets

Devices, techniques, and processes are disclosed that use electrical impedance to detect of the presence and contents of droplets including cells, nucleic acids, proteins, or solute concentrations in an array of retrievable, trackable, trapped droplets in a fluidic system. Electrodes may be positioned underneath individual droplet traps in a microchannel to assay droplet contents and/or actuating droplets for the release of the droplets from corresponding traps. The disclosed technology may be used for detection of the results of solvent extraction processes including time-dependent quantification of metal ion concentration in the aqueous and organic phases, for wastewater treatment, heavy metal detection, pharmaceutical industry, and/or biotechnology, or for environmental monitoring of wastewater for toxic metal, monitoring of biological cell viability and proliferation, monitoring of extraction processes used in heavy metal mining, monitoring of extraction processes used in nuclear fuel processing, monitoring kinetics of enzyme processes, and/or assessing pharmacodynamics and drug efficacy.

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

Stretched exponential magnetic relaxation dynamics in artificial square ice revealed through x-ray photon correlation spectroscopy

X-ray photon correlation scattering measurements are undertaken on a thermally active artificial spin ice based on the square lattice, referred to as artificial square ice, to probe the fluctuation timescales as a function of temperature as the system passes through the paramagnetic-antiferromagnetic phase transition, which belongs to the two-dimensional Ising universality class. In the paramagnetic regime, a single exponential timescale is seen, whereas at and below the critical temperature, a stretched exponential decorrelation is observed, with the stretching exponent decreasing from unity down to below one-half as the temperature reduces. This trend is confirmed by kinetic Monte Carlo simulations of a simplified point-dipolar square ice system, and is in agreement with past theoretical work on the kinetic Ising model where stretched exponential relaxation due to equilibrium domain wall dynamics below the critical temperature was found.

36 MATERIALS SCIENCE↗
DOE OSTI2022

An electrochemical mesoscale tool for modeling the corrosion of structural alloys by molten salt

Understanding the impact of microstructure on corrosion rates can aid the development of corrosion-resistant alloys for molten salt reactors. Here in this work, we develop an electrochemical phase-field model for capturing the microstructure-dependent corrosion of structural alloys by molten salts. As a demonstration problem, we apply this model to capture the selective depletion of Cr from Ni-Cr grain boundaries during corrosion in molten FLiBe salt. We perform sensitivity analysis and model verification on 1D simulations to confirm that the model predicts diffusion-limited kinetics. The model is validated using 1D, 2D, and 3D simulations against experimental data for Ni-5Cr and Ni-20Cr corrosion in molten FLiBe. The 1D simulations predict the corrosion behavior with reasonable accuracy when using an effective diffusion coefficient that accurately represents the grain boundary diffusion. 2D simulations that represent the grain structure underpredict the corrosion. 3D simulations that represent the grain structure predict the corrosion with reasonable accuracy. The corrosion rate predicted by the 3D simulations is proportional to the average grain size at the alloy/salt interface.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2025

Single-Crystallization of O3-Type Layered Oxide Cathode for Na-Ion Battery

The development of high-energy-density Na-ion batteries places significant demands on single-crystal layered oxide cathodes, especially for further high-voltage, solid-state battery scenarios. In the O3-type structure, due to the original sluggish Na ion diffusion kinetics (approximately 1 order of magnitude lower than that of Li-ion), and further hindrance against diffusion kinetics caused by single-crystal architecture, these inherent defects lead to the decline in the electrochemical performance. Herein, we demonstrated that the single crystallization of O3-type NaNi 1/3 Fe 1/3 Mn 1/3 O 2 cathode (d 50 = 5.04 μm) aggravates surface-to-bulk phase inhomogeneity distribution, which is attributed to the uneven Na ions extraction. Moreover, the Na-depletion of the surface/shell region not only aggravates Na ion diffusion resistance but also leads to a higher valence state of transition-metal elements (e.g., Ni/Fe) near the surface of the single-crystal particle, which further compromises the cathode-electrolyte interface stability. Furthermore, not limited to revealing the challenges, tuning the particle size and moderating quasi-single-crystal strategies have been proven to effectively mitigate the negative uneven distributions of Na ions, phases, and valence/oxidative states, resulting in efficient modification for single crystallization of Na-layered oxide cathodes.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Lattice distortion inducing exciton splitting and coherent quantum beating in CsPbI 3 perovskite quantum dots

It is well known that shape or crystal anisotropy in nanocrystal quantum dots, which are tiny semiconductor nanoparticles, results in energy splitting of their optically bright excitons (bound electron-hole pairs), which is known as fine structure splitting (FSS). Studying FSS usually requires single or just a few quantum dots because of its sensitivity to nanocrystal size and shape. In this work, by using ensemble-level femtosecond polarized transient absorption, the authors observed clear bright-exciton FSS in solution-processed CsPbI 3 perovskite QDs, which is manifested as exciton quantum beats (periodic oscillations of kinetic traces). Remarkably, the beat frequency of a given sample, as determined by the FSS energy, can be continuously controlled by changing the temperature. The FSS is linked experimentally for the first time to temperature dependent lattice distortion in the orthorhombic crystal structure. Because these orthorhombic-phase nanocrystal quantum dots are bounded by the pseudocubic family of crystal planes, the lattice distortion results in a misalignment between the lattice vectors and the nanocrystal edge directions which breaks the symmetry, leading to previously unrecognized avoided crossing “fine-structure gap” between bright exciton. Furthermore, this fine structure gap is responsible for the ability to observe quantum beats despite the shape inhomogeneity within the measured ensemble. Lattice-distortion dependent FSS in CsPbI 3 QDs has important consequences for their application as a platform for coherence and entanglement control in quantum information science.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Investigating Gas Dynamics and Reactions in Supersonic Jets

UF6 is highly reactive and readily hydrolyzes in contact with moisture to form uranyl fluoride (UO2 F2 ) and hydrofluoric acid (HF), however there is disagreement between experimental and computation investigations into the intermediates of the reaction. • As the hydrolysis reaction occurs at very fast time scales, it is difficult to investigate the intermediates and reaction kinetics experimentally. • Previous experiments at SRNL have had some success at halting and studying the UF6 hydrolysis reaction using cryogenic techniques, however isotopic substitution studies limit possible identity of intermediate compounds. • Primary goal of this project is to capture and study reaction intermediates produced during the gas phase hydrolysis of UF6 via supersonic expansion.

Waldron, Abigail M.↗
DOE OSTI2020

Coupling of Pronghorn and RELAP-7 for a Pebble Bed Reactor

High temperature gas cooled reactors (HTGR) are a candidate for timely Gen-IV reactor technology deployment because of high technology readiness and walk-away safety. Among HTGRs, pebble bed reactors (PBRs) have attractive features such as low excess reactivity and online refueling. Pebble bed reactors pose unique challenges to analysts and reactor designers such as continuous burnup distribution depending on pebble motion and recirculation, radiative heat transfer across a variety of gas-filled gaps, and long design basis transients such as pressurized and depressurized loss of forced circulation. Modeling and simulation is essential for both the PBR’s safety case and design process. In order to verify and validate the new generation codes the Nuclear Energy Agency (NEA) Data bank provide a set of benchmarks data together with solutions calculated by the participants using the state of the art codes of that time. An important milestone to test the new PBR simulation codes is the OECD NEA PBMR-400 benchmark which includes thermal hydraulic and neutron kinetic standalone exercises as well as coupled exercises and transients scenarios. In this work, the reactor multiphysics code MAMMOTH and the thermal hydraulics code Pronghorn, both developed by the Idaho National Laboratory (INL) within the multiphysics object-oriented simulation environment (MOOSE), have been used to solve Phase 1 exercises 1 and 2 of the PBMR-400 benchmark. The steady state results are in agreement with the other participants’ solutions demonstrating the adequacy of MAMMOTH and Pronghorn for simulating PBRs.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2023

Inter-Area Oscillation Damping with Type-5 Wind Power Plant: Preprint

This paper investigates the potential of using brushless excitation (BLE) for not only riding through the fault but also to damp inter-area power oscillation with the help of a synchronous generator (SG) used in Type 5 Wind Power Plant (WPP). In BLE, an auxiliary synchronous generator (ASG) behaving like an exciter is coupled and driven by the rotor of the main SG. The BLE's field current's ASG is fed by two separate loops of the automatic voltage regulator (AVR) and Power System Stabilizer (PSS). The AVR system implements a control loop to regulate the generator terminal voltage, VT. For the PSS, the kinetic energy of the wind turbine is utilized according to the estimated rotational speed of the synchronous generator shaft of the SG. It mitigates the necessity of any curtailment of active power for damping. The effectiveness of the proposed control scheme is verified with a three-phase short circuit fault in a two-area power system.

brushless↗
DOE OSTI2021

Manipulating the Crystallization Kinetics by Additive Engineering toward High-Efficient Photovoltaic Performance

Additive processing is proven to be an effective method to improve the efficiency and stability of perovskite solar cells; however, its intrinsic role in directing the crystallization pathway and thus morphology formation remains unknown. In situ grazing-incidence wide-angle x-ray scattering (GIWAXS) is applied to study the function of a 1,8-diiodooctane (DIO) additive in manipulating the crystallization behavior of perovskite thin films. It is seen that the DIO additive could induce multi-stage intermediate crystallization phases and increases the activation energy for nucleation and growth, which postpones the perovskite phase transformation time and broadens the transition zone. The elongated crystallization process affords improved perovskite thin film crystallinity and reduces defect density, which enables a longer carrier diffusion length. As a result, improved device efficiency, moisture, and thermal stability can be achieved. The following study provides a new prospective in understanding the additive function in perovskite thin film morphology control from fundamental parameters, indicating the importance of minor processing conditions in global property management toward high device performance.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Hybrid particle-spectral method for kinetic plasma simulations

A hybrid model for numerical solutions of the Vlasov–Poisson equations is presented, which blends spectral and particle approaches. The model splits the distribution function for plasma species into both spectral and particle representations in the velocity space to combine the advantages of each approach. The spectral representation leverages asymmetrically weighted Hermite basis, whereas the particle representation leverages the particle-in-cell method. Configuration phase space is decomposed with the Fourier method, which is well suited for periodic problems. We derive conservation equations for mass, momentum, and energy for the proposed combined method. It is shown that the coupling error between the two methods is absent in the semi-discrete setting (not taking into account time discretization). Finally, numerical test cases are presented simulating a weak electron beam interaction with plasma, leading to beam–plasma instability. The initially localized electron beam evolved into a highly non-equilibrium distribution function in the velocity space. A small growth rate and the resonance nature of instability make it difficult to obtain accurate solutions for purely particle methods due to noise, which falls as ∼1/Np with a number of particles. At the same time, purely spectral methods may require a large number of modes to capture the highly non-equilibrium state of the evolved beam. We show that the hybrid method is well suited for such problems: it reproduces the linear stage as well as nonlinear dynamics with sufficient accuracy using a highly non-equilibrium distribution function.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Hybrid particle-spectral method for kinetic plasma simulations

A hybrid model for numerical solutions of the Vlasov–Poisson equations is presented, which blends spectral and particle approaches. The model splits the distribution function for plasma species into both spectral and particle representations in the velocity space to combine the advantages of each approach. The spectral representation leverages asymmetrically weighted Hermite basis, whereas the particle representation leverages the particle-in-cell method. Configuration phase space is decomposed with the Fourier method, which is well suited for periodic problems. We derive conservation equations for mass, momentum, and energy for the proposed combined method. It is shown that the coupling error between the two methods is absent in the semi-discrete setting (not taking into account time discretization). Finally, numerical test cases are presented simulating a weak electron beam interaction with plasma, leading to beam–plasma instability. The initially localized electron beam evolved into a highly non-equilibrium distribution function in the velocity space. A small growth rate and the resonance nature of instability make it difficult to obtain accurate solutions for purely particle methods due to noise, which falls as ~1/$\sqrt{N_p}$ with a number of particles. At the same time, purely spectral methods may require a large number of modes to capture the highly nonequilibrium state of the evolved beam. We show that the hybrid method is well suited for such problems: it reproduces the linear stage as well as nonlinear dynamics with sufficient accuracy using a highly non-equilibrium distribution function.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Stabilizing High–Nickel Cathodes with High–Voltage Electrolytes

Electrolytes connect the two electrodes in a lithium battery by providing Li + transport channels between them. Advanced electrolytes are being explored with high-nickel cathodes and the lithium-metal anode to meet the high energy density and cycle life goals, but the origin of the performance differences with different electrolytes is not fully understood. Here, the mechanisms involved in protecting the high-capacity, cobalt-free cathode LiNiO 2 with a model high-voltage electrolyte (HVE) are delineated. The kinetic barrier posed by a thick surface degradation layer with poor Li + -ion transport is found to be the major contributor to the fast capacity fade of LiNiO 2 with the conventional carbonate electrolyte. In contrast, HVE reduces the side reactions between the electrolyte and the electrodes, leading to a thinner nano-interphase layer comprised of more beneficial species. Crucially, the HVE leads to a different surface reorganization pathway involving the formation of a thinner nanoscale LiNi 2 O 4 spinel phase on the LiNiO 2 surface. Finally, with a high 3D Li + -ion and electronic conductivity, the spinel LiNi 2 O 4 reorganization nanolayer preserves fast Li + transport across the cathode–electrolyte interface, reduces reaction heterogeneity in the electrode and alleviates intergranular cracking within secondary particles, resulting in superior long-term cycle life.

25 ENERGY STORAGE↗
DOE OSTI2025

Directed energy deposition of functionally graded V-4Cr-4Ti to Fe-9Cr transition for fusion power systems

This study proposes a graded structure via additive manufacturing for divertor and first wall blanket applications in fusion reactors. Materials were selected based on thermodynamic calculations to operate from 1100 °C at the plasma-facing level to 550 °C at the structural steel level. Conventional joining methods often lead to failures due to discrete reaction layers with significant mechanical property differences. Using laser beam-directed energy deposition (LB-DED), this study demonstrates the fabrication of a VCrTi-Gr91 steel functionally graded component through a novel process parameter optimization framework. A systematic approach included powder characterization, single-track depositions, and construction of printability maps. Near full-density specimens of each interlayer were additively manufactured, and a transition from V-based alloys to reduced activation ferritic martensitic steels was achieved. Computational material selection of interlayer alloys and thermodynamic/diffusion kinetics simulations prevented most interface incompatibilities. A brittle intermetallic formed at one interface, causing cracking, which was not predicted by current thermodynamic models. Transition alloy design approach was updated with a more recent database and a mitigation strategy has been proposed to eliminate the formation of deleterious intermetallic phases. Ultimately, LB-DED has proven effective for producing multi-material graded systems for fusion applications, with the demonstrated process parameter optimization framework applicable to various materials.

Additive manufacturing↗