Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Kinetic parameters”

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

High and Ultra-High Temperature Reaction Kinetics by Single Nanoparticle Mass Spectrometry

Methodology is presented for non-destructive, optically-detected single nanoparticle (NP) mass spectrometry, with the goal of extracting surface reaction kinetics for single NPs at high temperatures. Methods for determining the NP charge, mass, and temperature as a function of time are discussed, and the data are used to extract both the absolute kinetics for mass change, as well as the efficiencies of the surface processes that cause them. Factors that contribute to the uncertainties in absolute and relative mass determination, and in the resulting kinetic parameters, are discussed. The method allows the NP-to-NP variations in initial reactivity to be measured directly, along with the time evolution of reactivity resulting from NP structural/compositional changes that occur under reaction conditions. The strengths and limitations of single nanoparticle mass spectrometry as a high temperature surface kinetics tool are discussed in the context of sublimation and O2 oxidation kinetics experiments for single hafnium (Hf) NPs at temperatures ranging above 2400 K. The Hf oxidation kinetics are compared to analogous oxidation experiments for silicon, graphite, and carbon black NPs. In all four cases, the oxidation chemistry was dominated by processes that result in net mass loss, and the distinct mechanisms responsible are discussed. All four NPs also eventually passivated, i.e., the efficiencies for oxidative etching decreased by at least two orders of magnitude, relative to the initial efficiencies. Furthermore, the passivation mechanisms, which are quite different for carbon, compared to silicon or hafnium, are discussed. Carbon NP passivation is attributed to structural isomerization leading to fully coordinated, fullerene-like NP surfaces, while for silicon and hafnium, passivation results from delayed formation of an oxide layer, triggered by accumulation of oxygen in the NP sub-surface region.

36 MATERIALS SCIENCE↗
DOE OSTI2025

dynamics of organic-mineral interactions at the metal oxide-solution interface as studied via binding energetics (Final report)

This project focused on addressing longstanding fundamental and experimental uncertainties on how dissolved organic substances (DOS) interact with metal oxide surface under environmentally relevant conditions. By leveraging a custom-built real-time, in-tandem flow adsorption microcalorimetry-UV-Vis/fluorescence spectroscopy platform, we characterized the binding energetics, kinetics and mechanistic pathways driving DOS-metal oxide interactions at temporal resolution on the order of 1-5 seconds. We studied a diverse suite of model organic compounds/substances – including monocarboxylates (e.g. acetate and benzoate), di-carboxylates (oxalate and succinate), amino acids, amino-based nanparticles and natural organic matter – interacting at the mineral-water interface of structurally- and/or chemically distinct metal oxides (including SiO2, boehmite, ferrihydrite, and γ-Al2O3). Our results indicated that DOS-metal oxide interactions are governed by multi-step reaction pathways, often switching between distinct, resolvable enthalpy- and entropy-driven non-electrostatic or electrostatic configurations. To quantify these interactions, we developed and implemented an analytical workflow that integrates peak deconvolution and Monte-Carlo based error propagation to determine site-specific thermodynamic and kinetic parameters for individual binding/debinding events. In addition to resolving apparent first-order rate constants of each event, we were able to quantify associated apparent equilibrium constants as well as free energy, enthalpy and entropy contribution to the activation and subsequent progression of the binding/debinding process across compounds, compound class and metal oxide surfaces. The kinetic-thermodynamic data produced in this study captured how the interplay between oxide surface reactivity and DOS molecular structure jointly drives binding-debinding dynamics. Notably, that at pH below PZC of the oxide surface, neutral species were heavily involved in monocarboxylate binding, while anionic species drove dicarboxylate binding. Also, that among amino acids 1) positional isomers show distinctive binding characteristics to each other while enantiomers show no significant differences in binding characteristics, 2) molecules that bind via outer-sphere complexation show a larger entropic shift between binding and debinding with no impact on oxide surface while 3) inner-sphere interactions increased anion exchange capacity of the oxide surface. The new insights and data from this work has great potential for improving predictive modeling of carbon dynamics and specifically organic-mineral interactions in environmental and industrial systems.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Design considerations for a Digital Twin built to improve nitrification performance at a water resource recovery facility

A Digital Twin built around Activated Sludge Model No. 1 was deployed at a full-scale water resource recovery facility. Its design included a waste rate recommender system based on automatic scenario analyses, where influent loads and waste rates are varied to determine their impact on nitrification. At the same frequency as these scenario analyses, scheduled auto-calibrations allow for nitrifier maximum specific growth rate (μmax-NITO) soft sensing, the only kinetic parameter shown to require adjustment if the objective is aeration tank effluent ammonia forecasting accuracy. By integrating temperature forecasting over the next three sludge ages, this Digital Twin approach creates opportunities for advancing waste rate decisions in anticipation of seasonal temperature changes, optimizing ammonia control authority under varying influent loads, and furnishing valuable insights for future capital projects requiring nitrifier kinetic understanding and modelling.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Pitfalls in parameters: practical process development in chemical vapor processing of SiC

Vapor processing of high-value materials, such as silicon carbide (SiC), is of interest for many industries, including aerospace and energy production. Chemical vapor infiltration (CVI) of additively manufactured components is an especially attractive manufacturing process currently in development. Here, a novel development workflow is demonstrated with the H 2 –CH 3 SiCl 3 gas system for SiC to accelerate the process optimization of CVI SiC. A combination of calculated thermodynamics and high-throughput experimental chemical vapor deposition (CVD) coatings substantially reduced the experiments required with slow CVI processes. The computational results accurately predicted changes in the thermodynamic conditions tested, while CVD coatings – characterized by Raman spectroscopy – addressed changes in kinetic parameters. This workflow is also applicable to other vapor-processing systems, such as pyrolytic carbon, ZrC, or Si 3 N 4 .

Lamm, Benjamin W. [Oak Ridge National Laboratory (↗
DOE OSTI2024

Mechanistic and Kinetic Analysis of Complete Methane Oxidation on a Practical PtPd/Al 2 O 3 Catalyst

A PtPd/Al 2 O 3 catalyst developed for the complete oxidation of methane from the ventilation air of underground coal mines is compared against a model PdO/Al 2 O 3 catalyst. Although the PtPd/Al 2 O 3 catalyst is substantially more active and stable than the model catalyst, the nature of active sites between the two catalysts is deemed to be fundamentally the same based on their response to different feed gas compositions and the evolution of surface CO adsorption complexes during time-resolved CO adsorption DRIFTS experiment. For both catalysts, coordinatively unsaturated Pd sites are considered the active centers for methane activation and the subsequent oxidation reaction. H 2 O competes with CH 4 for the same active sites, resulting in severe inhibition. Additionally, the CH 4 oxidation reaction also causes self-inhibition. Taking both inhibition effects into consideration, a relatively simple kinetic model is developed. The model provides a good fit of the 72 sets of kinetic data collected on the PtPd/Al 2 O 3 catalyst under practically relevant reaction conditions with CH 4 concentration in the range of 0.05–0.4%, H 2 O concentration of 1.0–5.0%, and reaction temperatures of 450–700 °C. Kinetic parameters based on the model suggest that the CH 4 activation energy on the PtPd/Al 2 O 3 catalyst is 96.7 kJ/mol, and the H 2 O adsorption energy is –31.0 kJ/mol. Both values are consistent with the parameters reported in the literature. The model can be used to develop catalyst sizing guidelines and be incorporated into the control algorithm of the catalytic system.

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

Kinetics Modeling and Reactor Design Study of Glucose-to-Terpenes Cell-Free Conversion

Cell-free systems offer many advantages over traditional biological conversion by eliminating biological growth constraints. It also offers easy manipulation and finetuning of the reaction conditions for each individual enzyme. The conversion of cellulosic glucose to Limonene, a terpene, is a promising pathway for producing fuels and chemicals. Recent advances in developing cell-free systems focuses on bench scale optimization of terpene yield and to demonstrate its feasibility towards commercialization [1,2]. There is significant knowledge gap regarding reaction kinetics of these cell-free systems to further study how it will perform at larger scale. We present here, our studies on reaction kinetics and reactor design implications of cell-free glucose to Limonene conversion to facilitate the further development and commercialization of this process. We developed a novel kinetic model based on the metabolic-network structure of the cell-free system with multi-substrate reversible Michaelis-Menten rate law. To estimate kinetic parameters for this system of rate equations, we employed Bayesian optimization to perform global search with the assistance of gaussian processes to balance exploration and exploitation. The model parameters estimated showed good results compared with experimental data. The estimated parameters were used to perform sensitivity analysis. We found that Hexokinase is one of the most critical enzymes that affect the conversion of the glucose. We also observed that abundance of co-factors is also critical to the conversion of glucose to limonene. We investigated packed bed reactors with enzymes immobilized on the surface of particles to convert glucose stream into Limonene for larger scale production. The reactor design such as particle size, enzyme loading, and flow rate are found to be critical for improving yields. [1] Dudley, Q.M., Nash, C.J. and Jewett, M.C., 2019. Synthetic Biology, 4(1), p.ysz003. [2] Korman, T.P., Opgenorth, P.H. and Bowie, J.U., 2017. Nature communications, 8(1), p.15526.

09 BIOMASS FUELS↗
DOE OSTI2025

Kinetics of Pyrolysis and Thermal Evolution of Negev Desert Lithologies

The Negev desert in Israel is home to large quantities of organic-rich, shallow marine sedimentary lithologies that could potentially accommodate the disposal of spent nuclear fuel. Previous thermal analyses of Negev carbonates have focused on industrially relevant considerations such as natural gas and oil extraction or pyrolysis for recovering hydrocarbon fuels. Here, this study addresses thermal evolution of the Negev organic-rich carbonate, siliceous, and phosphorite rocks and associated chemical, mineralogical, and microstructural changes that may occur under prolonged thermal loading in the vicinity of spent nuclear fuel disposal systems. Our employed methods include high-temperature X-ray diffraction, high-temperature infrared spectroscopy, and thermal analysis integrating thermogravimetry, differential scanning calorimetry, and mass spectrometry. Further, we apply iterative iso-conversional model-free methods to derive kinetic parameters for thermal decomposition of the Negev organic-rich carbonate rocks from 200 to 550 °C. Our results have provided mechanistic insights into the thermal evolution encompassing water desorption, decomposition of organic matter, and decarbonation of carbonate phases.

58 GEOSCIENCES↗
DOE OSTI2025

SPARTA: A flux adjustment methodology to interpret complex experiments

For the accurate determination of reactivity from a detector count rate, correction of spatial effects is of prime importance. This spatial correction is often provided using simulation methodologies, but this may introduce a bias if the result of the experiment is also used as input data for the simulation. Here, this work presents a flux adjustment methodology able to infer experimental reactivity and correction of spatial effects without the need for a simulation. It can process the signal from a complex experiment such as a heat balance measurement in the TREAT reactor, where control rods are continuously adjusted to maintain a constant power. In the present work, this methodology successfully computed the reactivity and the local spatial variation of the flux of a generated signal. It also proved to be robust against noise and errors on kinetic parameters and provides a credible interpretation of a heat balance experiment in TREAT. Efficiency of flux adjustment methods for complex experiment enable a better experiment interpretation less reliant on nuclear data evaluation.

73 - NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Enzyme property prediction using artificial intelligence

Artificial intelligence (AI)-driven enzyme property prediction enables rapid discovery and engineering of enzymes for a wide range of biotechnological and therapeutic applications. Here, we first introduce the key components in AI model development, including enzyme datasets, protein representation methods, and model architectures. We then highlight a variety of AI tools developed for the prediction of enzyme properties and functional annotations, including enzyme structure, kinetic parameters, substrate specificity, thermostability, solubility, Enzyme Commission number, and Gene Ontology term. Moreover, we describe representative downstream applications enabled by these AI tools. Finally, we discuss some challenges and opportunities as well as future prospects.

Yuan, Le [University of Illinois at Urbana-Champai↗
DOE OSTI2025

Best of both worlds: Enforcing detailed balance in machine learning models of transition rates

The slow microstructural evolution of materials often plays a key role in determining material properties. When the unit steps of the evolution process are slow, direct simulation approaches such as molecular dynamics become prohibitive and Kinetic Monte-Carlo (kMC) algorithms, where the state-to-state evolution of the system is represented in terms of a continuous-time Markov chain, are instead frequently relied upon to efficiently predict long-time evolution. The accuracy of kMC simulations however relies on the complete and accurate knowledge of reaction pathways and corresponding kinetics. This requirement becomes extremely stringent in complex systems such as concentrated alloys where the astronomical number of local atomic configurations makes the a priori tabulation of all possible transitions impractical. Machine learning models of transition kinetics have been used to mitigate this problem by enabling the efficient on-the-fly prediction of kinetic parameters. While conventional KMC methods based on transition state theory naturally yield reversible dynamics that exactly obey the detailed balance criterion, providing strong guarantees on the properties of the stationary distribution, many recently-proposed ML-based approaches to barrier predictions provide no such guarantees. In this study, we derive conditions under which physics-informed ML architectures exactly enforce the detailed balance condition by construction, even when relying on non-extensive descriptions of states in terms of local environments around mobile defects. In conclusion, using the diffusion of a vacancy in a concentrated alloy as an example, we show that such ML architectures also exhibit superior performance in terms of prediction accuracy, demonstrating that the imposition of physical constraints can facilitate the accurate learning of barriers at no increase in computational cost.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Mutation of an active site-adjacent residue in VIM indirectly dictates interactions with and blunts inhibition by D-captopril

Activity assays and X-ray crystallographic studies were undertaken to elucidate the inhibitory mechanism of captopril stereoisomers on Verona integron-encoded metallo-β-lactamases, specifically VIM-20, VIM-31, and VIM-15. All three VIM-2-like variants (VIM-20, VIM-31, and VIM-15) and VIM-2 expressed in Escherichia coli exhibited catalytic activity with comparable steady-state kinetic parameters. Among the tested thiol drugs (L- and D-captopril, D,L-thiorphan, and 2,3-dimercaprol), IC50 analyses indicated that D-captopril and 2,3-dimercaprol were more potent inhibitors against the VIM enzymes examined in this study. Notably, the IC50 value of D-captopril against VIM-31 was an exception, closely resembling that of L-captopril. To elucidate this exceptional inhibitory potency of D-captopril and its binding mode in the active site of VIM-31, high-resolution crystal structures of VIM-20, VIM-31, and VIM-15 in complex with both L- and D-captopril are reported. These findings will help evaluate whether the identified potent inhibitor D-captopril could be further developed as a pan inhibitor targeting the VIM-family enzymes.

Silwal, Surendra Bikram↗
DOE OSTI2026

Hygrothermal aging effects on polyimide and acrylate-based adhesive materials in high-performance flat-flex cable assemblies

Flat-flex assemblies are widely used across a diverse range of technologies, including consumer electronics, automotive systems, aerospace and defense applications, and medical devices, due to their compact form factor and flexibility. However, environmental stressors like temperature and humidity are known to degrade performance over time, though the main mechanisms for degradation remain unknown. Here, in this work, we examined the aging behavior of copper/polyimide/adhesive laminate systems under varying temperature and humidity conditions, focusing on material degradation and its effects on mechanical and dielectric properties. Results show that peel strength and breakdown voltage decrease with exposure time, temperature, and humidity; spectroscopic characterization revealed that the adhesive, and not the polyimide, is the weak link in the material stack-up. The adhesive, identified as a butyl acrylate-acrylonitrile (BA-AN) copolymer, exhibited age-related spectral changes that correlated with exposure severity and duration. Hydrolysis at BA ester and AN nitrile groups was identified as the primary degradation mechanism, producing amides, acids, and alcohols. We developed a mechanistic model that connects kinetic parameters for BA-AN spectral band decay to peel strength and breakdown voltage, and implicates hydrolysis at the AN moiety, not BA, in performance decrements. Modeling predictions at ambient conditions indicate faster decline in peel strength than dielectric strength at early times, followed by plateauing of both properties at longer times. These findings provide critical insights into laminate aging mechanisms and highlight the importance of addressing BA-AN degradation to improve the long-term reliability of these systems in high-humidity environments.

organic↗
DOE OSTI2026

Impact of Rock Wettability and Mineralization on CO 2 Storage Efficiency in Basalt Reservoirs

Basalt formations have emerged as highly promising targets for geological CO 2 storage due to their abundance of reactive silicate minerals that can rapidly convert dissolved CO 2 into stable carbonate minerals. The efficiency of CO 2 trapping in basalts, however, depends not only on their geochemical reactivity but also on reservoir wettability, which governs how CO 2 partitions among structural, residual, solubility, and mineral trapping mechanisms. However, wettability and mineral kinetics have largely been examined separately, wettability for plume migration and kinetics in terms of mineralization, leaving their coupled impact unquantified. Here, we systematically integrate th ese processes in a reservoir-scale reactive transport model to quantify the impact of three distinct wettability conditions (water-wet, mixed-wet, and CO 2 -wet) on plume migration, mineral precipitation, pore-scale blockage, and long-term immobilization. Results show that under ideal conditions (no mineralization), water-wet systems exhibited plume aspect ratios up to 28% higher than CO 2 -wet systems and promoted at least 60% greater residual trapping. When mineralization was incorporated, mineral trapping accounted for 57%, 53%, and 46% of the total immobilized CO 2 in water-wet, mixed-wet, and CO 2 -wet systems, respectively. Dissolved CO 2 contributed 43-46%, while residual trapping remained below 5%, and mobile CO 2 was substantially reduced due to mineralization, demonstrating that as the system becomes more water-wet, CO 2 is more effectively converted into stable carbonate minerals, although extensive mineral precipitation may lead to pore clogging. A global Sobol-Morris sensitivity analysis further reveals that the kinetic parameters of reactive minerals are the most influential factor, followed by wettability, while salinity has only a minor influence. Furthermore, these findings highlight that wettability is not merely a flow property but a fundamental control on CO 2 mineralization in basalts, with direct implications for site selection for storage projects.

Basalts↗
DOE OSTI2025

Degrees of Rate Control in Interconnected Reaction Networks

Overall reactions in interconnected networks exhibit net, forward, and reverse rates that are governed by both constitutive elementary steps in the pathway of interest and branching elementary steps that lead to alternative products. Accordingly, steps in branching pathways exhibit negative net, forward, and reverse degrees of rate control, as they reduce reaction flux to the desired product. We here contextualize the forward and reverse degrees of rate control in terms of kinetic resistances (inverse of rates) and leverage the additive nature of kinetic resistance to decouple kinetic driving forces contributed by constitutive elementary steps and branching points (nodal species) in interconnected networks. Regardless of the network connectivity, forward and reverse degrees of rate control are shown to converge at equilibrium. Away from equilibrium, we identify two critical features of interconnected networks: stoichiometric regularity─condition where all stoichiometric numbers are unity─and pathway symmetry around nodal species─condition where branching pathways share the same rate constants, stoichiometry, and species concentrations/activities─that result in (i) equal forward, reverse, and consequently net degrees of rate control and (ii) forward and reverse degrees of rate control that exhibit constant offsets, respectively, across all extents of reaction. Furthermore, our discourse further provides a mathematical description for the influence of stoichiometric irregularity and pathway asymmetry on forward and reverse degrees of rate control. Altogether, the presented work details the effects of network (inter)connectivity and stoichiometry on reaction kinetics and, in doing so, establishes general protocols for capturing these effects as additive terms in the formulation of forward and reverse degrees of rate control.

10 SYNTHETIC FUELS↗
DOE OSTI2025

Catalytic Resonance Theory: Experimental and Kinetic Interpretation of Programmable Catalysis

The complexity of programmable catalysts that change in the physical or electronic state on the time scale of a catalytic turnover results in uncertainty in the kinetic outcome of common measurements with dynamic catalyst experiments. To provide an interpretation of utility to experimentalists, a general programmable catalytic reaction exhibiting high turnover efficiency was simulated to understand the kinetic behavior of catalysts with varying frequency, binding energy, oscillation amplitude, and minimum oscillation binding energy at different temperatures and reactant gas pressures. The time-averaged catalytic rate for varying temperature in the form of an Arrhenius plot exhibited three distinct kinetic regimes, with the intermediate temperatures or applied frequencies exhibiting zero slope, indicative of barrier-free kinetic control. Transitions in experimentally measurable kinetic regimes with both temperature and applied catalyst oscillation frequency were associated with transitions in the sensitivity of reaction rate constants and degrees of rate control. Furthermore, these distinct kinetic macroscopic features specific to programmable catalysts were identified for observation in experimental dynamic kinetic measurements.

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

A High-Fidelity Molecular Model of the Cu(111) Repeating Unit

Dynamic processes at surfaces are central to heterogeneous catalysis, but their atomistic mechanism(s) can prove difficult to elucidate due to variations in material structure and the corresponding impact on reactivity. Moreover, disparities between reaction conditions and those employed for spectroscopic characterization at surfaces can inhibit detailed understanding of catalysis-relevant chemistries. Herein, we substantiate the so-called “cluster-surface” analogy by leveraging a low-valent tricopper architecture ( 1 ) as a model system for small molecule activation at Cu(111). Two reaction classes are explored: the adsorption of carbon monoxide (CO) and the dissociative adsorption of dihydrogen (H 2 ). These processes serve as an ideal testbed to compare the reactivity of a molecular cluster ( 1 ) to that of a heterogeneous surface, as both reactions have empirical data from measurements performed on crystalline Cu(111). Cluster 1 reversibly binds CO. Variable temperature NMR analysis with 13 CO reveals a favorable enthalpy but large negative entropy (−5.1 kcal × mol –1 and −22.9 cal × mol –1 × K –1 , respectively) for CO binding, affording a process that is marginally endergonic at room temperature (ΔG ads (298.15 K) = 1.7 ± 0.5 kcal × mol –1 ). Similarly, analogous to a Cu(111) surface, 1 is shown to oxidatively add (chemisorb) H 2 . Kinetic parameters were determined for this process and the activation enthalpy (8.4 ± 0.5 kcal × mol –1 ) closely mirrors that established for H 2 binding at the Cu(111) facet (6.0 to 12.4 kcal × mol –1 ). Together, these results showcase that a trinuclear cluster can reproduce the small molecule binding and activation energetics of a bulk crystalline surface, setting the stage for studying less-defined surface processes in an atomically precise molecular setting.

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

Processing behavior evolution of recycled polypropylene: An integrated experimental and Computer-Aided engineering simulation study

Polypropylene (PP) comprises 21% of global plastics production and 18% of plastics waste, yet less than 1% of solid-waste PP is recycled in the United States (U.S.), representing significant environmental and economic challenges. Mechanical recycling, the most prevalent recycling method, subject's materials to thermomechanical stresses, which typically degrade polymer properties, affecting the quality of polymer products. This study replicates the impact of mechanical recycling through multiple extrusion cycles to examine the effects on PP's processing behavior. Dynamic scanning calorimetry (DSC) measurements showed stable melting behavior across all processing conditions, while crystallization analysis exhibited consistent shifts in kinetic parameters. Rheological characterization demonstrated progressive viscosity reductions through successive cycles, particularly pronounced at elevated reprocessing temperatures. Here, the integration of this experimental data into injection molding simulations showed that recycled PP maintains viable processing characteristics. Our findings establish quantitative correlations between processing history and material behavior, enabling optimization of processing parameters directly rather than relying on trial-and-error approaches. While these results reflect idealized recycling conditions with minimal contamination, they provide a framework for understanding fundamental property evolution during mechanical recycling.

42 ENGINEERING↗
DOE OSTI2025

Transient Modeling and Simulation of a Generic Stable Salt Reactor

A SAM system-level model of a generic stable salt reactor has been developed to investigate thermal-hydraulic behavior and safety performance under steady and transient conditions. The model integrates information generated from a reactor physics analysis using PROTEUS and PERSENT, and a computation fluid dynamics (CFD) analysis using STAR-CCM+. A loose, iterative coupling scheme between PROTEUS and SAM is implemented to calculate the equilibrium power and temperature distributions in the steady-state critical core condition. The converged steady-state model is then used in PERSENT to calculate the four reactivity feedback temperature coefficients (Doppler, fuel density, coolant density, and core radial expansion) and kinetic parameters that are needed in SAM to model the temperature feedback effects in transient simulations. Within the fully enclosed liquid fuel pins, natural convection is the dominant heat transfer mechanism. The STAR-CCM+ model of the fuel pin considers conjugate heat transfer from the liquid fuel salt to the pin cladding and external reactor coolant. The CFD results of the axial and radial temperature profiles are used to empirically determine an effective fuel salt thermal conductivity in the SAM fuel pin model so that the temperatures predicted by the SAM model match as closely as possible the CFD results. In the central region of the fuel pin, the effective thermal conductivity is as high as similar to 60 times the physical fuel salt thermal conductivity. The whole-plant SAM model is then used to simulate an unprotected station blackout transient. The results of this simulation showed that the large negative fuel axial expansion reactivity feedback reduces fission power to similar to 2.4% nominal power. The core is cooled by natural circulation, which removes heat in the core to the emergency heat removal system, and ultimately, to the ambient. However, peak fuel salt and cladding temperatures can potentially reach as high as 1500 K, albeit briefly, if the shutdown mechanism fails to operate.

stable salt reactor; transient simulations; system↗