Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Temperature reactivity coefficient”

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 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Methanol oxidation up to 100 atm in a supercritical pressure jet-stirred reactor

Methanol (CH 3 OH) has attracted considerable attention as a renewable fuel or fuel additive with low greenhouse gas emissions. Methanol oxidation was studied using a recently developed supercritical pressure jet-stirred reactor (SP-JSR) at pressures of 10 and 100 atm, at temperatures from 550 to 950 K, and at equivalence ratios of 0.1, 1.0, and 9.0 in experiments and simulations. Here, the experimental results show that the onset temperature of CH 3 OH oxidation at 100 atm is around 700 K, which is more than 100 K lower than the onset at 10 atm and this trend cannot be predicted by the existing kinetics models. Furthermore, a negative temperature coefficient (NTC) behavior was clearly observed at 100 atm at fuel rich conditions for methanol for the first time. To understand the observed temperature shift in the reactivity and the NTC effect, we updated some key elementary reaction rates of relevance to high pressure CH 3 OH oxidation from the literature and added some new low-temperature reaction pathways such as CH 2 O + HO 2 = HOCH 2 O 2 (RO 2 ), RO 2 + RO 2 = HOCH 2 O (RO) + HOCH 2 O (RO) + O 2 , and CH 3 OH + RO 2 = CH 2 OH + HOCH 2 O 2 H (ROOH). Although the model with these updates improves the prediction somewhat for the experimental data at 100 atm and reproduces well high-temperature ignition delay times and laminar flame speed data in the literature, discrepancies still exist for some aspects of the 100 atm low-temperature oxidation data. In addition, it was found that the pressure-dependent HO 2 chemistry shifts to lower temperature as the pressure increases such that the NTC effect at fuel-lean conditions is suppressed. Therefore, as shown in the experiments, the NTC phenomenon was only observed at the fuel-rich condition where fuel radicals are abundant and the HO 2 chemistry at high pressure is weakened by the lack of oxygen resulting in comparatively little HO 2 formation.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2025

Impacts of molecular architecture on the radiation-induced degradation and reaction kinetics of hydrophobic diglycolamides with the solvated electron and the dodecane radical cation

Given their proposed use as trivalent actinide–lanthanide separation ligands, the role of molecular architecture on the radiation robustness of diglycolamide (DGA) molecules has been investigated. This study examined three prototypical molecules with differences in their aliphatic chain architecture: N,N,N′,N′-tetra(n-octyl)diglycolamide (TODGA), N,N,N′,N′-tetra(2-ethylhexyl)diglycolamide (T2EHDGA), and N,N′-dimethyl-N,N′-dioctyldiglycolamide (DMDODGA). Rate coefficients and activation parameters are reported for the reactivity of each DGA with the solvated electron (e solv − ) and the corresponding dodecane radical cation (RH˙ + ) over the temperature range of 10.0 to 44.1 °C. These measurements indicate that DMDODGA is the most chemically reactive with both transient radicals, which may be attributed to this molecule's more accessible backbone. Complementary gamma dose accumulation studies (≤ 600 kGy) under envisioned process conditions—50 mM DGA in n-dodecane solvent—afforded dose constants for the loss of DGA of d = (3.41 ± 0.07) × 10 −3 , (4.19 ± 0.09) × 10 −3 , and (4.65 ± 0.09) × 10 −3 kGy −1 for T2EHDGA, DMDODGA, and TODGA, respectively. These dose constants indicate that varying DGA architecture affords subtle differences in chemical reactivity, leading to varying rates of radiolytic degradation under envisioned actinide–lanthanide separation conditions. However, more ambitious DGA frameworks, such as modifying the backbone, branching of the aliphatic chains, and/or changing the size of the chain may be required for larger gains in radiolytic longevity while optimizing actinide–lanthanide selectivity.

Arrhenius parameters↗
DOE OSTI2021

Quantum reactive scattering calculations for the cold and ultracold Li + LiNa → Li 2 + Na reaction

A first-principles based quantum dynamics study of the Li + LiNa(ν = 0, j = 0) → Li 2 (ν', j') + Na reaction is reported for collision energies spanning the ultracold (1 nK) to cold (1 K) regimes. A full-dimensional ab initio potential energy surface for the ground electronic state of Li 2 Na is utilized that includes an accurate treatment of the long-range interactions. The Li + LiNa reaction is barrierless and exoergic and exhibits a deep attractive potential well that supports complex formation. Thus, significant reactivity occurs even for collision temperatures approaching absolute zero. The reactive scattering calculations are based on a numerically exact time-independent quantum dynamics methodology in hyperspherical coordinates. Total and rotationally resolved rate coefficients are reported at 56 collision energies and include all contributing partial waves. Several shape resonances are observed in many of the rotationally resolved rate coefficients and a small resonance feature is also reported in the total rate coefficient near 50 mK. Furthermore, of particular interest, the angular distributions or differential cross sections are reported as a function of both the collision energy and scattering angle. Unique quantum fingerprints (bumps, channels, and ripples) are observed in the angular distributions for each product rotational state due to quantum interference and shape resonance contributions. The Li + LiNa reaction is under active experimental investigation so that these intriguing features could be verified experimentally when sufficient product state resolution becomes feasible for collision energies below 1 K.

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

ARDP Natrium Neutronic Methodology: Argonne Neutronic Assessment of ABR-1000

The Natrium Sodium-cooled Fast Reactor (SFR) concept developed by TerraPower, LLC was selected as one of the advanced reactor designs for demonstration under the Advanced Reactor Demonstration Program Natrium Demo project. In collaboration with TerraPower, the Argonne National Laboratory (ANL) team provide independent reviews of the neutronic methodology being used for the Natrium design. As part of this activity, both the TerraPower and ANL teams agreed to perform independent neutronics analyses of the ABR-1000 reactor problem as a demonstration on the usage of the neutronics methodologies. Comparison of the results from TerraPower and ANL would expose any methodological differences in the modeling approaches. This report describes and presents the Argonne analysis work on the ABR-1000, and will be provided to TerraPower, along with data files, to complete the comparison study. The comparison will cover the methodology to evaluate the reactivity coefficients for a safety analysis, the shutdown margin, and the reconstructed pin power distribution. The methods and models used to calculate the kinetic parameters and the reactivity coefficients (density, Doppler, axial/radial expansion, etc.) are described in detail with recommendations of specific method options and modelling techniques. The control system reactivity worth and the excess reactivity due to temperature defect are calculated for the shutdown margin analysis. The method to evaluate the excess reactivity is discussed in detail. The pin power reconstruction methodology used in SE2ANL, SE2RCT, and DASSH are discussed in this report. The pin power distribution evaluated by different methods are compared and discussed in the report.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
NASA NTRS1994

Optical and Electrical Properties of Reactively Sputtered TiN, ZrN, and HfN Thin Films

Thin films of titanium, zirconium, and hafnium nitride are prepared by DC magnetron reactive sputtering at room temperature on fused silica, optical glass, and silicon substrates. Deposition parameters are investigated in order to obtain stoichiometric films. The optical and electrical properties of the films as a function of nitrogen partial pressure are determined. The results show that an inverse correlation exists between the optical reflectance and the electrical resistivity of the films. Ellipsometer data for all three films show that their refractive index in the visible spectrum is decreased by increasing the film thickness while the extinction coefficient is unchanged. Thin films of TiN have the lowest room temperature resistivity relative to ZrN and HfN thin films.

Reactive↗
DOE OSTI2022

Investigating the kinetic effect of prenol on iso -octane auto-ignition by means of an experimental and modeling study

A detailed experimental and kinetic modeling study was dedicated to understand the reported octane hyperboosting effect of prenol, by means of the measurement of the ignition delay times of its blends with iso -octane, and measurement of the mole fraction profiles of the fuels and intermediates inside the ULille rapid compression machine. Here, these results show that prenol addition leads to a reduction of the first-stage ignition phenomena and negative temperature coefficient behavior, which is only qualitatively captured by the model and is consistent with knock resistance improvement. It is suggested that this behavior is caused by two different factors. The first originates from gas-phase reactivity of prenol, and spans from the formation of unreactive unsaturated species through resonance-stabilized radicals, thereby constituting a competitive pathway for the radical pool generated by iso -octane. The second is of catalytic nature and cannot be captured by means of gas-phase kinetic modeling, but could also play an important role in the behavior of prenol in internal combustion engines.

42 ENGINEERING↗
DOE OSTI2021

Inelastic, exchange, and reactive processes in rovibrationally excited collisions of HD with H

ABSTRACT The HD molecule is an important coolant in early universe chemistry models and a tracer of H2 in star-forming regions. Rate coefficients for collisional excitation and de-excitation of HD rotational and vibrational levels form important ingredients in astrophysical models. While collisions with He, H2, and H are the most important, available data for H + HD collisions are largely limited to temperatures less than 1000 K for the vibrational ground state, low-lying rotational levels of the v = 1 HD vibrational level, or computed without reactive contributions. Here, through explicit quantum scattering calculations, we report extensive data for rovibrational transitions in HD induced by H atoms for a range of rotational levels in v = 1 and some v = 0 levels for temperatures up to 1000 K. The significance of the computed results for astrophysical modeling is discussed.

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

Development and application of two-step uncertainty propagation and sensitivity analysis methodology for fast reactor safety analysis

Uncertainty quantification (UQ) in nuclear reactors for transients is directly linked with safety assessment through the cross-sections uncertainties, provided as a covariance matrix, which are propagated through the reactor system to output of interest pertaining to reactor safety, such as peak temperatures in fuel/clad/coolant. Using a two-step approach, uncertainties are first quantified and propagated from basic input variables (such as reaction cross-sections) to intermediate quantities (such as reactivity feedback coefficients) through lattice level calculations. Uncertainties of intermediate quantities (from the first step) are then propagated through the system transient calculations, in the second step, to obtain uncertainties on reactor safety output parameters of interest. The scope of this work consists of Uncertainty Quantification & Propagation of nuclear data uncertainties that are highly correlated through unprotected transient overpower and unprotected loss of flow to assess their impact on core safety parameters. This two-step approach in the presence of covariance renders the sensitivity analysis very challenging. In fact, usually the sensitivity analysis is restricted to each step, which limits its application since the sensitivities between the system output quantities and the basic input variables are difficult to obtain. Here, in this work, we address this issue by proposing a simple, general methodology to combine the sensitivity indices obtained in each step by assuming the model behavior being linear. For the first step Generalized Perturbation theory based indices are used while in the second step the recently studied Johnson indices. The uncertainty quantification and sensitivity methodologies discussed here are demonstrated on a generic LFR design which is based on the 500 MWth demonstration Lead-cooled fast reactor (DLFR) using oxide fuel, developed by Westinghouse Electric Company (WEC).

42 - ENGINEERING↗
NASA NTRS1993

Comprehensive investigation of HgCdTe metalorganic chemical vapor deposition

The principal objective of this experimental and theoretical research program was to explore the possibility of depositing high quality epitaxial CdTe and HgCdTe at very low pressures through metalorganic chemical vapor deposition (MOCVD). We explored two important aspects of this potential process: (1) the interaction of molecular flow transport and deposition in an MOCVD reactor with a commercial configuration, and (2) the kinetics of metal alkyl source gas adsorption, decomposition and desorption from the growing film surface using ultra high vacuum surface science reaction techniques. To explore the transport-reaction issue, we have developed a reaction engineering analysis of a multiple wafer-in-tube ultrahigh vacuum chemical vapor deposition (UHV/CVD) reactor which allows an estimate of wafer or substrate throughput for a reactor of fixed geometry and a given deposition chemistry with specified film thickness uniformity constraints. The model employs a description of ballistic transport and reaction based on the pseudo-steady approximation to the Boltzmann equation in the limit of pure molecular flow. The model representation takes the form of an integral equation for the flux of each reactant or intermediate species to the wafer surfaces. Expressions for the reactive sticking coefficients (RSC) for each species must be incorporated in the term which represents reemission from a wafer surface. The interactions of MOCVD precursors with Si and CdTe were investigated using temperature programmed desorption (TPD) in ultra high vacuum combined with Auger electron spectroscopy (AES). These studies revealed that diethyltellurium (DETe) and dimethylcadmium (DMCd) adsorb weakly on clean Si(100) and desorb upon heating without decomposing. These precursors adsorb both weakly and strongly on CdTe(111)A, with DMCd exhibiting the stronger interaction with the surface than DETe.

Raupp, Gregory B.↗
NASA NTRS1999

Thermophysical Property Measurements of Molten Semiconductors in 1-g and Reduced-g Condition

Understanding and controlling the formation kinetics of varieties of crystal imperfections such as point defects, non uniform distribution of doping atoms, and impurity atoms in growing crystals are very important. Theoretical (numerical) modeling of the crystal growth process is an essential step to achieving these objectives. In order to obtain reliable modeling results, input parameters, i.e. various thermophysical parameters, must be accurate. The importance of accurate thermophysical properties of semiconductors in crystal growth cannot be overly emphasized. The total hemispherical emissivity, for instance, has a dramatic impact on the thermal environment. It determines the radiative emission from the surface of the melt which determines to a large extent the profile of the solidified crystal. In order to understand the convection and the turbulence in a melt, viscosity becomes an important parameter. The liquid surface tension determines the shape of the liquid-atmosphere interface near the solid-liquid-atmosphere triple point. Currently used values for these parameters are rather inaccurate, and this program intends to provide more reliable measurements of these thermophysical properties. Thus, the objective of this program is in the accurate measurements of various thermophysical properties which can be reliably used in the modeling of various crystal growth processes. In this program, thermophysical properties of molten semiconductors, such as Si, Ge, Si-Ge, and InSb will be measured as a function of temperature using the High Temperature Electrostatic Levitator at JPL. Each material will be doped by different kinds of impurities at various doping levels. Thermophysical properties which will be measured include: density, thermal expansion coefficient, surface tension, viscosity, specific heat, hemispherical total emissivity, and perhaps electrical and thermal conductivities. Many molten semiconductors are chemically reactive with crucibles. As a result, these dispersed impurities in the melts tend to substantially modify the properties of pure semiconductors. Sample levitation done in a vacuum clearly helps maintain the sample purity. However, in the 1-g environment, all gravity caused effects such as convection, sedimentation and buoyancy are still present in the sample. In addition, large forces needed to levitate a sample in the presence of the gravity can cause additional flows in the melt. The use of the High Temperature Electrostatic Levitator (HTESL) for the present research is a recent development and little is known about the flows induced by the electrostatic forces. In this ground base program, we will define the limits of HTESL technology as various thermophysical properties of molten semiconductors are measured.

Rhim, Won-Kyu↗
DOE OSTI2019

Modeling reactivity insertion experiments of TRISO particles in NSRR using BISON

This research presents BISON simulations of stress, temperature, and failure probability of tristructural isotropic (TRISO) fuel particles subjected to transient power pulse conditions in the Nuclear Safety Research Reactor (NSRR). We report that by modifying the default elastic properties of the PyC and the UO 2 coefficient of thermal expansion correlation, BISON was found to produce suitable agreement with the observed and independently simulated results when appropriate shape and scale Weibull parameters for SiC failure were chosen. The Weibull parameter effects were also explored within the range of SiC failure and were applied using experimental data from hemispherical crush testing.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2022

Transient Radiation-Induced Berkelium(III) and Californium(III) Redox Chemistry in Aqueous Solution

Despite the significant impact of radiation-induced redox reactions on the accessibility and lifetimes of actinide oxidation states, fundamental knowledge of aqueous actinide metal ion radiation chemistry is limited, especially for the late actinides. A quantitative understanding of these intrinsic radiation-induced processes is essential for investigating the fundamental properties of these actinides. We present here a picosecond electron pulse reaction kinetics study into the radiation-induced redox chemistry of trivalent berkelium (Bk(III)) and californium (Cf(III)) ions in acidic aqueous solutions at ambient temperature. New, and first of a kind, second-order rate coefficients are reported for the transient radical-induced reduction of Bk(III) and Cf(III) by the hydrated electron (eaq−) and hydrogen atom (H•), demonstrating significant reactivity (up to 1011M-1s-1) indicative of a preference of these metals to adopt divalent states. Additionally, we reportthe first-ever second-order rate coefficients for the transient radical-induced oxidation of these elements by reaction with hydroxyl (•OH) and nitrate (NO3•) radicals, which also exhibited fast reactivity (ca. 108 M-1s-1). Transient Cf(II), Cf(IV), and Bk(IV) absorption spectra are also reported. Overall, the presented data highlight the existence of rich, complex, intrinsic late actinide radiation-induced redox chemistry that has the potential to influence the findings of other areas of actinide science.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS2001

Methanol Uptake by Low Temperature Aqueous Sulfuric Acid Solutions

The global methanol budget is currently unbalanced, with source terms significantly larger than the sinks terms. To evaluate possible losses of gaseous methanol to sulfate aerosols, the solubility and reactivity of methanol in aqueous sulfuric acid solutions representative of upper tropospheric and lower stratospheric aerosols is under investigation. Methanol will partition into sulfate aerosols according to its Henry's law solubility. Using standard uptake techniques in a Knudsen cell reactor, we have measured the effective Henry's law coefficient, H*, for cold (196 - 220 K) solutions ranging between 45 and 70 wt % H2SO4. We have found that methanol solubility ranges from approx. 10(exp 5) - 10(exp 7) M/atm for UT/LS conditions. Solubility increases with decreasing temperature and with increasing sulfuric acid content. Although methanol is slightly more soluble than are acetone and formaldehyde, current data indicate that uptake by clean aqueous sulfuric acid particles will not be a significant sink for methanol in the UT/LS. These solubility measurements include uptake due to physical solvation and any rapid equilibria which are established in solution. Reaction between primary alcohols and sulfuric acid does occur, leading to the production of alkyl sulfates. Literature values for the rate of this reaction suggest that formation of CH3OSO3H is not significant over our experimental time scale for solutions below 80 wt % H2SO4. To confirm this directly, results obtained using a complementary equilibrium measurement technique will also be presented.

Iraci, L. T.↗
DOE OSTI2022

Replicating HCCI-like autoignition behavior: What gasoline surrogate fidelity is needed?

This work seeks to characterize the fidelity needed in a gasoline surrogate with the intent to replicate the complex autoignition behavior exhibited within advanced combustion engines, and specifically Homogeneous Charge Compression Ignition (HCCI). A low-temperature gasoline combustion (LGTC) engine operating in HCCI mode and a rapid compression machine (RCM) are utilized to experimentally quantify fuel reactivity, through autoignition and preliminary heat release characteristics. Fuels considered include a research grade E10 U.S. gasoline (RD5-87), three multi-component surrogates (PACE-1, PACE-8, PACE-20), and a binary surrogate (PRF88.4). Each fuel was studied at lean/HCCI-like conditions covering a wide range of temperatures and pressures that are representative of naturally aspirated to high boost engine operation. Detailed chemical kinetic modeling is also undertaken using a recently updated gasoline surrogate kinetic model to simulate the RCM experiments and to provide chemical insight into surrogate-to-surrogate differences. The LGTC engine experiments demonstrate nearly identical reactivity between PACE-20 and RD5-87 across studied conditions, while faster phasing is seen for both PACE-1 and PACE-8 due to their stronger intermediate- and low-temperature heat release (ITHR/LTHR) at naturally aspirated and boosted conditions, respectively. The RCM experiments reveal typical low-temperature, negative temperature coefficient (NTC) and intermediate-temperature autoignition behaviors at all pressure conditions for RD5-87, which are qualitatively reproduced by all surrogates. Quantitative discrepancies in both autoignition and preliminary heat release are observed for all surrogates, while their ability to replicate RD5-87 autoignition behavior follows the order of PACE-20 > PACE-1 > PACE-8 > PRF88.4. Excellent mapping is obtained between the LGTC engine and the RCM, where the engine pressure-time trajectories can be characterized by the regimes represented by the RCM autoignition isopleths. The kinetic model performs commendably when simulating both autoignition and preliminary heat release of PACE-20, while typically overpredicting ignition delay times for PACE-1, PACE-8 and PRF88.4 at high-pressure and low-temperature/NTC conditions. Sensitivity and rate of production (ROP) analyses highlight surrogate-to-surrogate differences in the governing chemical kinetics where n-pentane initiates rapid OH branching at a faster rate and an earlier timing for PACE-20 than iso-pentane does for PACE-1 and PACE-8, making it computationally more reactive than the other surrogates. The current study highlights the need to include non-standardized properties, such as the lean/HCCI-like autoignition characteristics, in addition to ASTM properties (e.g., RON, MON) as metrics of fuel reactivity and targets to be matched when formulating high-fidelity surrogates that fully capture gasoline advanced combustion behavior such as HCCI-like autoignition.

42 ENGINEERING↗
DOE OSTI2024

Xenon compatibility in magmatic processes: Hadean to current contexts

Xenon (Xe) behaviour in petrological processes, albeit essential to constrain mantle ingassing and degassing models, is elusive due to its volatile nature, and lack of direct investigation at the pressures (P) and temperatures (T) relevant to magma formation and crystallisation at depth. Xenon stands out amongst noble gases due to its unique reactivity with silicates of the lower crust and upper mantle, which could at least partially explain that published mineral/melt partitioning coefficients span up to six orders of magnitude. We report partition coefficients of Xe using in situ X-ray fluorescence at high P and T, and mass spectrometry analyses. Xenon is found to be moderately incompatible in anorthite-clinopyroxene mix in equilibrium with basalt (partition coefficient value of 0.16 ± 0.06), and compatible in olivine in equilibrium with basalt (partition coefficients in the range 88 ± 22 to 302 ± 46). While Xe is, thus, concentrated in basaltic melts coexisting with crystallising pyroxenes and feldspars, it is strongly retained in olivine at depth. Consequently, Xe originally contained in solid Earth has been preferentially retained at depth throughout Earth’s history, from the magma ocean stages to present day partial mantle melting processes.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1981

Review of literature surface tension data for molten silicon

Measurements of the surface tension of molten silicon are reported. For marangoni flow, the important parameter is the variation of surface tension with temperature, not the absolute value of the surface tension. It is not possible to calculate temperature coefficients using surface tension measurements from different experiments because the systematic errors are usually larger than the changes in surface tension because of temperature variations. The lack of good surface tension data for liquid silicon is probably due to its extreme chemical reactivity. A material which resists attack by molten silicon is not found. It is suggested that all of the sessile drip surface tension measurements are probably for silicon which is contaminated by the substrate materials.

Hardy, S.↗
DOE OSTI2023

Sensitivity study on the fin effect of thermocouple mounted on the heated surface under film boiling condition

In the reactivity-initiated accident (RIA) scenario, right after the power burst with the occurrence of departure from nucleate boiling (DNB), the cladding temperature quickly reaches its maximum and is maintained in a steady film boiling regime. Then, under limited heat transfer performance of film boiling, the local cooling effect by the thermocouple wire may not be neglected. A computational model is developed to estimate the local cooling effect of the thermocouple wire mounted on the cladding surface under film boiling condition. The model predicts the temperature profile along the cladding surface as physically expected and is verified against the benchmark problem. Sensitivity studies reveal the significant effects of the thickness of vapor film layer and the diameter of thermocouple on the temperature drop near the thermocouple, for they increase the heat transfer area to the liquid water region. The impact of film boiling heat transfer coefficient, however, is not significant, as the primary heat transfer mechanism near the thermocouple is the heat conduction to the wire. Future study will include the contribution of radiative heat transfer between the cladding surface and the vapor film layer, as well as a transient power history which represents an actual RIA scenario.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2026

Ion correlations explain kinetic selectivity in diffusion-limited solid-state synthesis reactions

Establishing viable solid-state synthesis pathways for novel inorganic materials remains a major challenge in materials science. Previous pathway design methods using pairwise reaction approaches have navigated the thermodynamic landscape with first-principles data but lack kinetic information, limiting their effectiveness. This gap leads to suboptimal precursor selection and predictions, especially for reactions forming competing phases with similar formation energies, where ion diffusion is a critical influence. Here we demonstrate an inorganic synthesis framework by incorporating machine learning-derived transport properties through ‘liquid-like’ product layers into a thermodynamic cellular reaction model. In the Ba–Ti–O system, known for its competitive polymorphism, we obtain accurate predictions of phase formation with varying BaO:TiO2 ratios as a function of time and temperature. We find that diffusion–thermodynamics interplay governs phase compositions, with cross-ion transport coefficients critical for predicting diffusion-limited selectivity. This work bridges length scales and timescales by integrating solid-state reaction kinetics with first-principles thermodynamics and spatial reactivity.

Atomistic models↗