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,225 records · Page 68
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Numerical analysis of soot emissions from gasoline-ethanol and gasoline-butanol blends under gasoline compression ignition conditions

In the present work, computational fluid dynamics (CFD) simulations of a single-cylinder gasoline compression ignition (GCI) engine were performed to investigate the impact of blending two biofuels, ethanol and n-butanol, with gasoline on the trade-off between combustion phasing and soot emissions under low load conditions. Here, in order to represent market gasoline (RD5-87), a four-component toluene primary reference fuel (TPRF) + ethanol (ETPRF) surrogate (with 20% ethanol by mole; E20) was formulated using a neural network based octane predictor such that the surrogate had the same ethanol content, Research Octane Number (RON) and Octane Sensitivity (S). In addition, a novel skeletal kinetic mechanism for ETPRF and TPRF + n-butanol (BTPRF) blends, incorporating polycyclic aromatic hydrocarbon (PAH) chemistry, was developed. A three-dimensional (3D) engine CFD formulation employing the skeletal mechanism, adaptive mesh refinement (AMR), finite-rate chemistry approach, and hybrid method of moments (HMOM) was adopted to capture the in-cylinder combustion phenomena and soot emissions. The engine CFD model was validated against RD5-87 experimental data for a broad range of start-of-injection (SOI) timings (-21/-27/-36/-45 crank angle degrees (CAD) after top-dead center (aTDC)), with respect to in-cylinder pressure, heat release rate, combustion phasing, and soot emissions. The closed-cycle simulation results were analyzed to elucidate the non-monotonic trend of soot emissions versus SOI timing: SOI-36 > SOI-45 > SOI-21 > SOI-27. Thereafter, the validated CFD model was employed to simulate the combustion of a gasoline-ethanol blend with 45% (by mole) ethanol (E45) and a gasoline-butanol blend with 45% (by mole) n-butanol (B45) under the same operating conditions to study the effects of fuel composition and SOI timing on combustion phasing and soot emissions. The sooting propensity followed the trend: B45 > E20 > E45 at all SOI timings. Overall, it was observed that the autoignition propensity was primarily related to fuel chemistry. On the other hand, sooting propensity showed strong coupling with both fuel chemistry and physical properties, with greater impact of fuel physical properties at advanced SOI timings.

30 DIRECT ENERGY CONVERSION↗
DOE OSTI2025

Orbit-averaging and deposition accuracy for runaway electron beams in hybrid kinetic-MHD simulations of the runaway plateau

We develop a new procedure that combines the kinetic orbit runaway electron code (KORC) and the NIMROD extended-magnetohydrodynamic code to simulate runaway electrons (REs) in the post-disruption plateau. KORC integrates guiding-center orbits, with a barycentric-based binary search strategy providing initial guesses for the Newton–Raphson logical-to-physical coordinate inversion, ensuring reliable particle-to-mesh mapping in NIMROD, whose fields remain static for the present study. Samples are drawn in accord with experimental parallel current profiles of RE beams during the plateau phase. Deposition in NIMROD is verified through comparison with a Python-based finite-element code that ensures periodicity in the poloidal direction and continuity at the magnetic axis. Accurate representation of near-axis fields requires finer mesh resolution to prevent under- and overshoots in current density from orbit inaccuracies. Yet, at a fixed particle count, increasing mesh resolution amplifies statistical noise in the deposited fields. An orbit-averaging method accumulates partial current deposits over multiple kinetic steps and reduces the statistical noise with little added computational cost. By coupling kinetic routines from KORC directly into the NIMROD codebase, these developments lay essential groundwork for future self-consistent KORC–NIMROD coupling.

Algorithms and data structure↗
DOE OSTI2022

Stability margin analysis of the Holos-quad microreactor design

A simplified point model of the Holos-Quad microreactor is introduced. The model is based on the point kinetics equations coupled to three heat balance equations, representing the mean temperatures of the fuel particles, the graphite moderator and the helium coolant. The differential equations are converted to the frequency domain, enabling the construction of the closed-loop reactor transfer function. Using this function the stability margins of the core design is analyzed for various power levels. It is shown that the gain margins approaches infinity, demonstrating the stability of the reactor for all power levels. The phase margin at nominal power is about 60 degrees, however it shows a non-monotonous dependence on power, with minimal value obtained for about 10% of nominal power. As power level is further increased, the phase margin also increases, demonstrating the reactor becomes more stable. This behavior may be of high importance in load-follow scenarios, where the power level of the reactor changes with time. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2020

Direct Observation of Alternating Octahedral and Prismatic Sodium Layers in O3-Type Transition Metal Oxides

The oxygen stacking of O3-type layered sodium transition metal oxides (O3-NaTMO 2 ) changes dynamically upon topotactic Na extraction and reinsertion. While the phase transition from octahedral to prismatic Na coordination that occurs at intermediate desodiation by transition metal slab gliding is well understood, the structural evolution at high desodiation, crucial to achieve high reversible capacity, remains mostly uncharted. In this work, the phase transitions of O3-type layered NaTMO 2 at high voltage are investigated by combining experimental and computational approaches. An OP2-type phase that consists of alternating octahedral and prismatic Na layers is directly observed by in situ X-ray diffraction and high-resolution scanning transmission electron microscopy. The origin of this peculiar phase is explained by atomic interactions involving Jahn–Teller active Fe 4+ and distortion tolerant Ti 4+ that stabilize the local Na environment. Finally, the path-dependent desodiation and resodiation pathways are also rationalized in this material through the different kinetics of the prismatic and octahedral layers, presenting a comprehensive picture about the structural stability of the layered materials upon Na intercalation.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Sr(Ag 1− x Li x ) 2 Se 2 and [Sr 3 Se 2 ][(Ag 1− x Li x ) 2 Se 2 ] Tunable Direct Band Gap Semiconductors

Abstract Synthesizing solids in molten fluxes enables the rapid diffusion of soluble species at temperatures lower than in solid‐state reactions, leading to crystal formation of kinetically stable compounds. In this study, we demonstrate the effectiveness of mixed hydroxide and halide fluxes in synthesizing complex Sr/Ag/Se in mixed LiOH/LiCl. We have accessed a series of two‐dimensional Sr(Ag 1− x Li x ) 2 Se 2 layered phases. With increased LiOH/LiCl ratio or reaction temperature, Li partially substituted Ag to form solid solutions of Sr(Ag 1− x Li x ) 2 Se 2 with x up to 0.45. In addition, a new type of intergrowth compound [Sr 3 Se 2 ][(Ag 1− x Li x ) 2 Se 2 ] was synthesized upon further reaction of Sr(Ag 1− x Li x ) 2 Se 2 with SrSe. Both Sr(Ag 1− x Li x ) 2 Se 2 and [Sr 3 Se 2 ][(Ag 1− x Li x ) 2 Se 2 ] exhibit a direct band gap, which increases with increasing Li substitution ( x ). Therefore, the band gap of Sr(Ag 1− x Li x ) 2 Se 2 can be precisely tuned via fine‐tuning x that is controlled by only the flux ratio and temperature.

Zhou, Xiuquan↗
DOE OSTI2022

Assessing entropy for catalytic processes at complex reactive interfaces

When chemical reactions are accelerated by a catalyst, entropy differences between reactants and their transient intermediates can be the driving force behind the promotion or inhibition of desired and parasitic chemical pathways. Understanding and controlling catalytic processes therefore requires both a fundamental and practicable understanding of entropy in addition to enthalpy. In unstructured media such as the vapor phase equilibrated with sparsely covered surfaces, entropy can be adequately accounted for by well-established approaches based on translational, rotational, and harmonic vibrational partition functions. However, these approximations become inadequate in more complex condensed phase environments, e.g., solid liquid interfaces of confined reaction spaces. In this chapter, we provide an overview of the state-of-art in the computational quantification of entropy and its known ramifications on catalysis. The fundamental roles of thermodynamics and kinetics in catalysis are covered in enough detail to appreciate and contextualize the computational methods employed to compute chemically accurate estimates of entropy. These methods are discussed in appropriate detail and range from the ubiquitous harmonic oscillator approximation where entropy unrelated to high frequency oscillations is typically underestimated, to enhanced free energy sampling with molecular dynamics where the desired accuracy must be weighed against the associated computational cost of obtaining it. The rising importance of machine learning and artificial intelligence in accelerating methodological progress in this field is touched upon, as well. Finally, applications, successes, and pitfalls of using these methods are provided to showcase past and present accomplishments while clarifying where improvements in both understanding and methodology are still needed.

Kollias, Loukas↗
DOE OSTI2023

Realistic large eddy and dispersion simulation experiments during project sagebrush phase 1

Typical dispersion simulations use mesoscale models at grid spacing greater than 1 km to simulate pollutant transport. Advances in computing power have allowed Large Eddy Simulations (LES), (grid spacing less than 100 m) to be run with boundary conditions in the grayscale regime (grid size of 1 km–100 m). We used the WRF mesoscale/LES model and the HYSPLIT particle dispersion model to simulate SF 6 transport from Project Sagebrush Phase 1, using a grayscale-aware scheme with 4 nests of 3 km, 1 km, 200 m, 67 m. The WRF/LES simulation was in satisfactory agreement with meteorological observations but the HYSPLIT SF6 simulations indicated insufficient sub-grid turbulence kinetic energy, especially in the lowest 100 m. This resulted in over predictions of the surface SF6 concentration and insufficient vertical transport. Best agreement with the Sagebrush observations was obtained with the WRF/LES large-scale and sub-grid turbulence, enhanced with a HYSPLIT parameterized turbulence profile.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Constraining the release of Sn to the ambient melting point following shock loading using time-resolved x-ray diffraction

The formation of liquid following release from a shocked state governs the transition from spall to cavitation and the formation of ejecta in metals. In order to build physics-based models of these processes, it is necessary to critically evaluate the relative importance of kinetics and entropy generation during the release along with the accuracy of multiphase equations of state. Tin (Sn) has served as a testbed for a variety of experiments examining strength and ejecta due to its accessible melt boundary and solid–solid phase transitions. This work presents experiments examining the phase evolution of high purity Sn following the shock and release to ambient stress near the melting point. Sn is found to release to states between its ambient solidus and liquidus from approximately 19 to 33 GPa under uniaxial loading, with the two-phase region being characterized by a reduction in the intensity of the (220), (211) β -Sn doublet. Jetting experiments performed at 27–28 GPa exhibit comparable diffraction patterns with what is observed following the uniaxial release. The solid fractions of β-Sn in the ambient mixed phase region are found to decrease linearly with increasing shock stress as increasing liquid Sn is formed. Here, the results provide much needed information for interpreting measurements of dynamic strength at a high strain rate and experiments examining cavitation and shallow bubble collapse in Sn.

, spall↗
DOE OSTI2020

Dendrite Growth Morphology Modeling in Liquid and Solid Electrolytes

The main goal of this project is to develop a multi-scale modeling approach that connects micron-scale phase-field models and atomic-scale density functional theory (DFT)-based simulations via parameter- and relationship-passing in order to predict Li-metal dendrite morphology evolution, in both liquid and solid electrolytes. The key hypothesis of the DFT-informed phase-field multiscale modeling approach is that it can capture the electrochemical-mechanical driving forces and incorporate the roles of nano-meter-thin solid electrolyte interphase (SEI) in liquid electrolytes as well as of the microstructures of micro-meter-thick solid electrolytes (SEs) for all-solid-state batteries. In this project, we have formulated and implemented phase-field models to incorporate the electrochemical driving forces in liquid electrolytes and then incorporate mechanical driving forces to simulate dendrite growth in solid electrolytes with resolved microstructures. We have implemented two treatments for the SEI: an explicit model to include the microstructure of the SE or SEI in the phase field model and an implicit model to simulate the impact of nano-meter thick SEI in liquid electrolytes by varying the electrode/electrolyte interfacial properties. The key interfacial properties, including the electronic and ionic transport properties, the charge transfer reaction kinetics, and mechanical properties, were computed by DFT-based calculations. At the DFT-based model, one key advancement is to directly predict the charge transfer reaction kinetics at a complex Li/SEI/electrolyte interface by linking DFT with density functional tight binding (DFTB) calculations. As the main accomplishments, we have demonstrated two successful predictions in both solid electrolyte and liquid electrolyte based on this multiscale approach. The predicted intergranular Li dendrite growth in LLZO revealed the importance of trapped electrons at internal interfaces in the microstructure of LLZO. The predicted electroplating morphology of mossy Li and faceted Mg agreed well with experiments. The insights provided by the multiscale model and the model enabled electrolyte and SEI design will accelerate the development of Li-metal electrode for high energy density batteries, that meet DOE’s target on cell density (>350 Wh/kg) and cost below $100/kWhuse for EV applications.

25 ENERGY STORAGE↗
DOE OSTI2024

1D modeling of plasma streamers at ammonia-air flame conditions

Abstract Self-consistent 1D modeling of streamers in ammonia-oxygen-nitrogen-water mixtures has been performed in this work. A fluid model that includes species transport, electrostatic potential, and detailed chemistry was developed and verified. This model is then used to simulate the avalanche, streamer formation and propagation phases, driven by a nanosecond voltage pulse, at different thermochemical conditions derived from a 1D laminar premixed ammonia-air flame. The applicability of the Meek’s criterion in predicting the streamer inception location was successfully confirmed. Streamer formation and propagation duration were found to vary significantly with different thermochemical conditions, due to the difference in ionization rates. The thermochemical state also affected the breakdown characteristics which was tested by maintaining the background reduced electric field constant. Detailed kinetic analyses revealed the importance of O ( 1 D ) in the production of key radicals, such as O, OH, and NH 2 . Furthermore, the contributions of the dissociative electronic excitation of NH 3 towards the production of H and NH 2 radicals have also been reported. Spatial and temporal evolution of the electron energy loss fractions for various inelastic collision processes at different thermochemical states uncovered the input plasma energy spent of fuel dissociation and the large variability in the dominant processes during the avalanche and streamer propagation phases. The methodology and analyses reported in this work are key towards developing effective strategies for controlled nanosecond-pulsed non-equilibrium plasma sources used for ammonia ignition and flame stabilization.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Bulk and surface diffusion of neodymium in alpha-uranium: Ab initio calculations and kinetic Monte Carlo simulations

A fundamental understanding of lanthanide transport in metallic fuels is critical for high fidelity modeling of the fuel-cladding chemical interaction (FCCI) phenomenon, which can lead to the formation of brittle intermetallic compounds and premature failure of the cladding. Here we report a combined ab initio density functional theory (DFT) and kinetic Monte Carlo (KMC) study of the bulk diffusivity of Nd in α-U, fully taking into account the effect of radiation enhanced diffusion. In this work, the vacancy mechanism is considered to be the dominant mechanism for the bulk diffusion of Nd since a Nd interstitial is found to be intrinsically unstable in α-U. The surface diffusivity of a Nd adatom on α-U (001) surface has been further predicted using KMC simulations parameterized by DFT calculations. The present study suggests that Nd transport via the surface diffusion mechanism can be many orders of magnitude faster than bulk diffusion. Furthermore, the results from the present lower length scale study can be used to inform mesoscale phase-field simulations to determine the effective diffusion coefficient of Nd through α-U with a porous microstructure.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Precise Control of Phase Separation Enables 12% Efficiency in All Small Molecule Solar Cells

Abstract Compared to conjugated polymers, small‐molecule organic semiconductors present negligible batch‐to‐batch variations, but presently provide comparatively low power conversion efficiencies (PCEs) in small‐molecular organic solar cells (SM‐OSCs), mainly due to suboptimal nanomorphology. Achieving precise control of the nanomorphology remains challenging. Here, two new small‐molecular donors H13 and H14 , created by fluorine and chlorine substitution of the original donor molecule H11 , are presented that exhibit a similar or higher degree of crystallinity/aggregation and improved open‐circuit voltage with IDIC‐4F as acceptor. Due to kinetic and thermodynamic reasons, H13 ‐based blend films possess relatively unfavorable molecular packing and morphology. In contrast, annealed H14 ‐based blends exhibit favorable characteristics, i.e., the highest degree of aggregation with the smallest paracrystalline π–π distortions and a nanomorphology with relatively pure domains, all of which enable generating and collecting charges more efficiently. As a result, blends with H13 give a similar PCE (10.3%) as those made with H11 (10.4%), while annealed H14 ‐based SM‐OSCs have a significantly higher PCE (12.1%). Presently this represents the highest efficiency for SM‐OSCs using IDIC‐4F as acceptor. The results demonstrate that precise control of phase separation can be achieved by fine‐tuning the molecular structure and film formation conditions, improving PCE and providing guidance for morphology design.

Bin, Haijun↗
DOE OSTI2023

Separating Xylene Isomers with a Calcium Metal‐Organic Framework

Abstract The purification of p‐xylene (pX) from its xylene isomers represents a challenging but important industrial process. Herein, we report the efficient separation of pX from its ortho‐ and meta‐ isomers by a microporous calcium‐based metal–organic framework material (HIAM‐203) with a flexible skeleton. At 30 °C, all three isomers are accommodated but the adsorption kinetics of o‐xylene (oX) and m‐xylene (mX) are substantially slower than that of pX, and at an elevated temperature of 120 °C, oX and mX are fully excluded while pX can be adsorbed. Multicomponent column breakthrough measurements and vapor‐phase/liquid‐phase adsorption experiments have demonstrated the capability of HIAM‐203 for efficient separation of xylene isomers. Ab initio calculations have provided useful information for understanding the adsorption mechanism.

Yu, Liang↗
DOE OSTI2023

Fast Charge Induced Phase Evolution and Element Contribution of Nickel-Rich Layered Cathode for Lithium-Ion Batteries

In nickel-rich layered cathode materials, three transition metals (TM = Ni, Co and Mn) play critical roles in Li storage performance. However, rate-dependent phase evolution and elemental contribution of nickel-rich cathode materials are not well understood, but very important for further design and development of these cathodes in high-power applications. Here, the rate-dependent phase evolution and elemental contribution of LiNi 0.5 Mn 0.3 Co 0.2 O 2 (NMC532) cathode are investigated thoroughly by using time-resolved synchrotron-based in situ X-ray diffraction and absorption techniques. The increase of Ni content in NMC532 has been found to be the main cause of the complex structural changes, resulting in distortion of TM-O 6 octahedron and strong static vibration between TM and O coordination. It is revealed that the fast charge (10C and 30C) of NMC532 leads to intermediate Li-poor phases, shrinking of H2 phase region, and prolonged O1 phase. During high-rate charging, Co is oxidized easily in low voltage region, while Ni mainly dominates charge compensation in high voltage region. It is found that delithiation-induced local structure distortion transfers from Ni to Co sites in fast kinetic process. In conclusion, these findings provide in-depth understanding for the fast charge behavior of Ni-rich layered cathode materials and help to guide further development of advanced high-power lithium-ion batteries.

25 ENERGY STORAGE↗
DOE OSTI2023

Unraveling two distinct polymorph transition mechanisms in one n-type single crystal for dynamic electronics

Cooperativity is used by living systems to circumvent energetic and entropic barriers to yield highly efficient molecular processes. Cooperative structural transitions involve the concerted displacement of molecules in a crystalline material, as opposed to typical molecule-by-molecule nucleation and growth mechanisms which often break single crystallinity. Cooperative transitions have acquired much attention for low transition barriers, ultrafast kinetics, and structural reversibility. However, cooperative transitions are rare in molecular crystals and their origin is poorly understood. Crystals of 2-dimensional quinoidal terthiophene (2DQTT-o-B), a high-performance n-type organic semiconductor, demonstrate two distinct thermally activated phase transitions following these mechanisms. Here we show reorientation of the alkyl side chains triggers cooperative behavior, tilting the molecules like dominos. Whereas, nucleation and growth transition is coincident with increasing alkyl chain disorder and driven by forming a biradical state. We establish alkyl chain engineering as integral to rationally controlling these polymorphic behaviors for novel electronic applications.

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

Tunable angle-dependent electrochemistry at twisted bilayer graphene with moiré flat bands

Tailoring electron transfer dynamics across solid–liquid interfaces is fundamental to the interconversion of electrical and chemical energy. Stacking atomically thin layers with a small azimuthal misorientation to produce moiré superlattices enables the controlled engineering of electronic band structures and the formation of extremely flat electronic bands. Here, we report a strong twist-angle dependence of heterogeneous charge transfer kinetics at twisted bilayer graphene electrodes with the greatest enhancement observed near the ‘magic angle’ (~1.1°). Furthermore, this effect is driven by the angle-dependent tuning of moiré-derived flat bands that modulate electron transfer processes with the solution-phase redox couple. Combined experimental and computational analysis reveals that the variation in electrochemical activity with moiré angle is controlled by a structural relaxation of the moiré superlattice at twist angles of <2°, and ‘topological defect’ AA stacking regions, where flat bands are localized, produce a large anomalous local electrochemical enhancement that cannot be accounted for by the elevated local density of states alone.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Continuum shock mixture models for Ni+Al multilayers: Individual layers and bulk equations of state

Continuum shock mixture models are reviewed and applied to determine the equations of state for five different compositions of Ni x Al y ⁠, as well as bulk Ni+Al reactive multilayers, by combining the fundamental property data for elemental nickel and aluminum. From the literature, we down-select and evaluate two analytical models for the mixture Hugoniot, i.e., the well-known method of kinetic energy averaging (KEA) and a recent model proposed by Jordan and Baer [J. Appl. Phys. 111, 083516 (2012)]. Fundamentally, the former method assumes pressure equilibrium, whereas the latter assumes a common particle velocity and mixture sound speed from compressible two-phase cavitating flows. Additionally, we construct thermodynamically complete equations of state by fitting Einstein oscillator series models for the specific heat at constant volume. Finally, the solid solution approximation is invoked for intermetallic compositions, which are not strictly physical mixtures. Overall, the KEA model provides a better fit to the available Ni x Al y and Ni+Al multilayer shock compression data; however, there are combinations of material properties where the performance of these two models is thought to be reversed. Moreover, the results of this work include the first analytical solution of Jordan–Baer that does not require numerical root finding, as well as proposed modifications to the Einstein oscillator series to incorporate some effects of local pressure–temperature equilibrium and reaction–diffusion. Future work is planned that will use these equations of state in mesoscale simulations to study shock-induced reaction in Ni+Al multilayers, and the intended application is illustrated with a brief 2D hydrocode example.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Small to Full-Height Scale Comparisons of Cesium Ion Exchange Performance with Crystalline Silicotitanate

Crystalline silicotitanate (CST) ion exchanger is planned to be used to remove cesium (137Cs) from the aqueous phases of Hanford tank wastes in preparation for vitrification at the Waste Treatment and Immobilization Plant (WTP) Pretreatment Facility. Column scale up testing was conducted to evaluate performance of Cs exchange onto the CST. Batch contact testing was conducted to assess exchange kinetics of Cs exchange at four different CST PSDs. Column testing was conducted at three column sizes, small (2.5% full height), medium (12% full height), and full height, to assess Cs load performance behavior. Testing at the small scale was compared to actual Hanford tank waste testing to verify the validity of simulant tests to accurately represent full height column performance. A change in CST particle size was essential in scaling the small column dynamics up to full scale. Based on these results, a determination of intraparticle and film diffusion impacts on overall mass transfer coefficients will allow future modeling of the breakthrough performance at a range of process conditions.

Westesen, Amy M.↗