Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Reactor 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 343 records · Page 19
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Measurements and kinetic modeling of O 2 vibrational kinetics in O 2 –Ar mixtures partially dissociated by a Ns pulse discharge

Vibrational kinetics of O 2 is studied during the O atom recombination in an O 2 –Ar mixture, partially dissociated by a burst of ns discharge pulses in a heated plasma flow reactor. The time-resolved temperature in the discharge afterglow is determined by Rayleigh scattering. Time-resolved O atom number density is measured by ps Two-Photon absorption Laser Induced Fluorescence, calibrated in xenon. Time-resolved vibrational level populations of molecular oxygen, O 2 (v= 8–20), are measured by ps Laser Induced Fluorescence (LIF), with the absolute calibration by NO LIF. Time-resolved ozone number density is monitored by broadband UV absorption. The results are compared with the predictions of a state-specific kinetic model. The experimental data indicate a rapid initial decay of O 2 (v) populations generated by electron impact in the discharge, due to the vibration-translation (V–T) relaxation by O atoms. This is followed by a slower population reduction, on the time scale much longer compared to that for V–T relaxation or vibration-vibration (V–V) exchange. Both O atoms and the O 2 (v) populations decay on the same time scale, indicating that chemical reactions initiated by the O atom recombination result in the generation of vibrationally excited O 2 molecules. These trends are reproduced by the kinetic model, which shows that the reaction of O atoms with ozone is the dominant pathway of O 2 (v) generation at the present conditions. The predicted relative O 2 (v) populations are close to the experimental results, but absolute number densities differ from the experimental data. This is likely due to uncertainties in the absolute calibration of LIF measurements and in the spectroscopic model used in the data reduction. The present work demonstrates the capability for the absolute, time-resolved measurements of vibrationally excited O 2 in recombining gas flows, to quantify the energy partition in the recombination reactions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Bridging Scales in Bioenergy and Catalysis: A Review of Mesoscale Modeling Applications, Methods, and Future Directions

Between the molecular and reactor scales, which are familiar to the chemical engineering community, lies an intermediate regime, here termed the “mesoscale,” where transport phenomena and reaction kinetics compete on similar time scales. Bioenergy and catalytic processes offer particularly important examples of mesoscale phenomena owing to their multiphase nature and the complex, highly variable porosity characteristic of biomass and many structured catalysts. In this review, we overview applications and methods central to mesoscale modeling as they apply to reaction engineering of biomass conversion and catalytic processing. A brief historical perspective is offered to put recent advances in context. Applications of mesoscale modeling are described, and several specific examples from biomass pyrolysis and catalytic upgrading of bioderived intermediates are highlighted. Methods including reduced order modeling, finite element and finite volume approaches, geometry construction and import, and visualization of simulation results are described; in each category, recent advances, current limitations, and areas for future development are presented. Owing to improved access to high-performance computational resources, advances in algorithm development, and sustained interest in reaction engineering to sustainably meet societal needs, we conclude that a significant upsurge in mesoscale modeling capabilities is on the horizon that will accelerate design, deployment, and optimization of new bioenergy and catalytic technologies.

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

Linear gyrokinetic stability of a high β non-inductive spherical tokamak

Spherical tokamaks (STs) have been shown to possess properties desirable for a fusion power plant such as achieving high plasma β and having increased vertical stability. To understand the confinement properties that might be expected in the conceptual design for a high β ST fusion reactor, a 1 GW ST plasma equilibrium was analyzed using local linear gyrokinetics to determine the type of micro-instabilities that arise. Kinetic ballooning modes and micro-tearing modes are found to be the dominant instabilities. The parametric dependence of these linear modes was determined and, from the insights gained, the equilibrium was tuned to find a regime marginally stable to all micro-instabilities at θ 0 = 0.0. This work identifies the most important micro-instabilities expected to generate turbulent transport in high β STs. The impact of such modes must be faithfully captured in first-principles-based reduced models of anomalous transport that are needed for predictive simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Insight into premixed diethoxymethane flames: Laminar burning velocities, temperatures, and emissions behaviour

Diethoxymethane ((CH 3 CH 2 O) 2 CH 2 , DEM) is a promising carbon-neutral fuel. DEM is a diether or acetal with a molecular structure similar to oxymethylene ethers (CH 3 O–(CH 2 O)n–CH 3 , OME n ). Thus, DEM can be expected to have a similar combustion behavior to OMEs, reducing harmful emissions such as NO x and particulate matter (PM) in internal combustion engines. From both experimental and kinetic modeling, fundamental studies on DEM are scarce in the literature. More studies are required to gain a detailed insight into the oxidation kinetics of DEM. Laminar burning velocity (LBV) is a critical property that allows a detailed assessment of the potential application of DEM in combustion devices. Unfortunately, the literature on the LBV of DEM is limited. Therefore, in this study we have investigated the LBV of DEM using two reactors for the first time, namely a heat flux burner and a combustion chamber. The experimental data is reported for equivalence ratio between 0.7 and 1.7, initial temperatures of 368–423 K, and initial pressure of 1–5 bar. In addition, we developed a detailed kinetic model extending our recent work of Shrestha et al. (Combust. Flame. 246 (2022) 112,426) to characterize the combustion behavior of DEM utilizing the new experimental data from this work and the literature data. Our model performs remarkably well in capturing the newly measured LBV experimental data over various experimental conditions. We found that DEM and dimethoxy methane (DMM) have similar values of LBVs (within ±1.5 cm/s) for a given condition, which indicates that intermediate chemistry governs the flame chemistry. Despite DEM being a larger molecule that is expected to have slightly lower LBVs than DMM, its effect on the measured values of LBVs is negligible. Finally, we experimentally measured NO x formation in DEM flame for the first time. The stochiometric flame has the highest NO x formation. The proposed model predicted the equivalence ratio dependence of NO x nicely. However, it overestimates the NO x formation for stoichiometric DEM/air mixtures by ~30 %. The model suggests that the thermal NO formation route is favored at lean and stochiometric conditions. In contrast, the prompt NO formation route is enhanced for rich mixtures.

10 SYNTHETIC FUELS↗
DOE OSTI2024

In-situ kinetic study of irradiation induced crystallization in amorphous Al 2 O 3

In the last ten years amorphous alumina coatings, deposited by Pulsed Laser Deposition, emerged as potential key enabling technology in the fields of heavy liquid metal fast reactors (lead and lead-bismuth) and fusion. In the former, as coating of the steel fuel cladding and in the latter as multifunctional coating providing a barrier against tritium permeation, steel corrosion and electrical insulation. Nevertheless, a detailed knowledge of the behavior of this thermodynamically metastable material at high temperatures and under neutron irradiation is still unknown. A knowledge gap that is mandatory to fill up for the deployment of this barrier technology. In the present work, we present a first step towards this goal, by the in-situ dynamic observation of the radiation induced crystallization processes of thin films of amorphous Al 2 O 3 , induced by ion-irradiation over an extensive range of temperatures (400-800 °C). The study was performed at the Intermediate Voltage Electron Microscope (IVEM)-Tandem Facility at Argonne National Laboratory. The experimental findings allow to elucidate the dependence of the grain growth on ion dose and temperature. A kinetic approach has been used to derive the process activation energies and other important parameters.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Transport and kinetics properties of LaF 3 in FLiNaK molten salt determined by electrochemical methods

It is found that lanthanum species can exist in the form of LaF 6 3- in FLiNaK (46.5 LiF-11.5 NaF-42 KF, mol%) molten salt in previous studies. The understanding of transport and reaction kinetics parameters of LaF 6 3- ions in FLiNaK molten salt is essential for the electrochemical separation of lanthanum from FLiNaK melt to achieve the coolant clean-up in fluoride salt cooled high-temperature reactor (FHR). Through conducting chronopotentiometry tests, the diffusion coefficients of LaF 6 3- ions in FLiNaK molten salt at the temperature range of 923 K–1023 K were determined. By conducting potentiodynamic polarization measurements at different concentrations and temperatures, the exchange current density, reaction rate constant, and charge transfer coefficient were obtained. The experiment data of potentiodynamic polarization were analyzed using a non-simplified electrode kinetics equation which incorporates both mass transfer and reaction kinetics.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
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

Catalytic co-pyrolysis of LDPE and PET with HZSM-5, H-beta, and HY: experiments and kinetic modelling

In this study, the catalytic pyrolysis of low-density polyethylene (LDPE), polyethylene terephthalate (PET) and their mixture (1 : 1 wt/wt) with three zeolite catalysts (HZSM-5, H-beta, HY) was investigated using a thermogravimetric analyzer (TGA) and a Pyroprobe® micro-reactor coupled to a gas chromatograph mass spectrometer (Py-GC/MS). The TGA results demonstrated that during pyrolysis at 10 °C min –1 , on average, zeolite catalysts decreased the maximum decomposition temperature by 149 °C for LDPE while only decreasing by 8 °C for PET. The derivative thermogravimetric (DTG) curve evidenced interactions when the two polymers were catalytically co-pyrolyzed for all the three catalysts. A lumped nth order reaction scheme was able to accurately model both non-catalytic and catalytic pyrolysis and co-pyrolysis by using the least squares fitting approach for determining the kinetic parameters. The kinetic model was able to model well the interaction effects observed during catalytic co-pyrolysis of LDPE and PET with HZSM-5, H-beta, and HY (Fit% Wt% > 96%, Fit% DTG > 94%). Py-GC/MS experiments for the catalytic fast pyrolysis of LDPE revealed that HZSM-5 resulted in the highest selectivity to aromatic hydrocarbons (31.6%) and HY resulted in the highest selectivity to gasoline range C 5 –C 10 paraffins and olefins (40.9%). Catalytic fast pyrolysis of PET showed high selectivity to benzene for all the catalysts (>43%) and that HZSM-5 resulted in the highest selectivity to polyaromatic hydrocarbons (24.7%). Finally, the catalytic fast co-pyrolysis of LDPE and PET revealed interaction effects for all the three catalysts evidenced by a positive synergy% for alkylated benzenes (3–142%) and polyaromatics (105–187%) with a concomitant negative synergy% for benzene (24–36%) and C 5 –C 10 paraffins and olefins (27–53%).

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

Neutronic Consideration of TREAT Facility Fuel SiC Recladding

The Transient Reactor Test (TREAT) Facility is an air-cooled reactor composed of graphite-urania fuel encapsulated in zirconium alloy canisters and surrounded by ~2 ft (~61 cm) of nuclear grade graphite. The primary role of TREAT is to provide in-pile nuclear-heated safety research via transient testing [1]. Benchmark efforts have demonstrated sensitivity of TREAT neutronics models to graphite cross section data and key neutron poisons found in the fuel assemblies [2]. During fabrication of TREAT fuel assemblies, a small concentration of boron was inadvertently introduced with variable concentration into the aggregate of fuel blocks prepared for TREAT use. The concentration of boron, and its respective uncertainty, represent the most sensitive component in a neutronics study; unfortunately, nothing can be done to reduce or otherwise mitigate its impact. Hydrogen in the fuel blocks reduces the mean free path of the core and its true quantity and uncertainty is approximated. The hydrogen remains from pitch binder used in the fuel fabrication process and was not completely released due to incomplete fuel graphitization. Destructive fuel testing is necessary to quantify the actual hydrogen content and to improve benchmark model quality. A third neutron poison of some significance is the hafnium impurity in the Zircaloy-3 (Zry-3) assembly cladding. The quantity of hafnium is expected to be within nominal manufacturing limits for nuclear grade zircaloy. Recladding the fuel assemblies with material that is less neutron absorbing is expected to benefit reactor neutron economy. Moreover, TREAT’s total transient energy is restricted based on the rate of Zry-3 oxidation in TREAT’s air coolant environment. Hence, cladding with more benign air reaction kinetics could dramatically uprate the facility’s experimental capability. Silicon carbide (SiC) has been proposed as a candidate material for recladding TREAT assemblies, as SiC/SiC demonstrate favorable properties within nuclear environments [3]. The current operational limits of TREAT are restricted by the oxidation of Zry-3. Change of the cladding to SiC could feasibly double the current temperature limitation. Recladding opportunities allow for previously damaged fuel assemblies to be returned to service, as Zry-3 is not commercially available. The initial effort required to consider proceeding with any recladding effort is to neutronically characterize key reactor parameters while replacing Zry-3 with SiC. A characterization summary, with respective results, is provided herein.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2023

Experimental and Computational Studies of the Production of 1,3-Butadiene from 2,3-Butanediol Using SiO2-Supported H3PO4 Derivatives

Silica-supported phosphoric acid and metal phosphate catalyzed 1,3-butadiene (BDE) production from 2,3-butanediol (2,3-BDO) was studied using experimental and computational techniques. The catalyst was initially tested in a continuous flow reactor using commercially available 2,3-BDO, leading to maximum BDE yields of 63C%. Quantum chemical mechanistic studies revealed 1,2-epoxybutane is a kinetically viable and thermodynamically stable intermediate, supported by experimental demonstration that this epoxide can be converted to BDE under standard reaction conditions. Newly proposed E2 and SN2' elementary steps were studied to rationalize the formation of BDE and all detected side-products. Additionally, using quantum mechanics/molecular mechanics (QM/MM) calculations, we modeled silica-supported phosphate catalysts to study the effect of the alkali metal center. Natural population analysis showed that phosphate oxygen atoms are more negatively charged in CsH2PO4/SiO2 than in H3PO4/SiO2. In combination with temperature-programmed desorption experiments using CO2, the results of this study suggest that the improved selectivity achieved when adding the metal center is related to an increase in the basicity of the catalyst.

butadiene↗
DOE OSTI2020

Accounting for erroneous model structures in biokinetic process models

In engineering practice, model-based design requires not only a good process-based model, but also a good description of stochastic disturbances and measurement errors to learn credible parameter values from observations. However, typical methods use Gaussian error models, which often cannot describe the complex temporal patterns of residuals. Consequently, this results in overconfidence in the identified parameters and, in turn, optimistic reactor designs. Here, we assess the strengths and weaknesses of a method to statistically describe these patterns with autocorrelated error models. This method produces increased widths of the credible prediction intervals following the inclusion of the bias term, in turn leading to more conservative design choices. However, we also show that the augmented error model is not a universal tool, as its application cannot guarantee the desired reliability of the resulting wastewater reactor design.

42 ENGINEERING↗
DOE OSTI2022

Development of a Gibbs Energy Minimiser for the MOOSE-based Corrosion Modelling App Yellowjacket and Validation of MSTDB

Nuclear materials are highly complex multiscale, multiphysics systems,and an effective prediction of nuclear reactor performance and safety requires simulation capabilities that tightly couple different physical phenomena. The Idaho National Laboratory’s Multiphysics Object Oriented Simulation Environment (MOOSE) provides the computational foundation for performing such simulations. With the move towards advanced reactors, such as the Molten Salt Reactor (MSR), that employ high temperature fluids compared to conventional reactors, corrosion has become a problem of great interest. A new application called Yellowjacket is currently under development to directly couple thermodynamic equilibrium and kinetics with phase field models in order to model corrosion in MSRs. As part of Yellowjacket, a Gibbs energy minimiser is being developed to perform thermochemical equilibrium calculations for a range of different materials, which is currently in its infancy. This report describes the further progress towards the development of Yellowjacket Gibbs energy minimiser. Ontario Tech University is developing a new Gibbs energy minimiser for Yellowjacket which is the primary contribution of this work. The aim to develop a thermochemistry solver for the MOOSE framework following the same development philosophy and using the same tools and libraries. A special focus is on performance, documentation and SQA. Furthermore through a scope extension partway through the fiscal year a thorough assessment of the MSTDB-TC v1.3 was performed at Ontario Tech University in the context of continuous improvement and quality assurance. The objective of the work was to have an arm’s length review of the database to give confidence that the database is performing as it was intended while assessing its current state to give recommendations to future developments. This assessment involved two parts: A) a quantitative assessment, and B) a qualitative assessment. Part A involved developing an automated test-suite that would compute values from the database using Thermochimica with comparisons to experimental measurements for validation purposes, which gives confidence to the database’s stakeholders that its functioning properly. Part B involved reviewing all binary systems in the database and making a qualitative assessment with two performance indicators: comprehensiveness and overall confidence. It is important to note that the models in the database are empirical, which is to say that the quality of any model is highly dependent on the experimental data used to inform its development.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2025

Gas-Cooled High-Temperature Pebble-Bed Reactor Reference Plant Model Updates

This work presents the latest improvements to, and investigations performed with, the pebble-bed high-temperature gas-cooled reactor (PB-HTGR) reference plant models for the United States Nuclear Regulatory Commission. These models serve as the foundation for the future development of detailed design evaluation models based on license applications. The reference plant models have been developed with the Comprehensive Reactor Analysis Bundle, or BlueCRAB, which is the code suite proposed for non-light-water reactor systems safety analysis. It incorporates various simulation tools developed by the Nuclear Energy Advanced Modeling and Simulation program, including the Griffin code for reactor physics, the Pronghorn and SAM codes for core thermal fluids, the BISON code for solid conduction and fuel performance, and the SAM code for system analysis. The primary objective of the work that was performed was to assess BlueCRAB’s level of readiness for modeling a PB-HTGR. To do so, we first developed numerical models in BlueCRAB that include the key physics for this technology to ensure an adequate level of fidelity for modeling PB-HTGR core performance and for performing multiphysics simulations for equilibrium core conditions and different accident scenarios. Then we simulated transient scenarios, including depressurized and pressurized loss of forced cooling accidents, over-cooling, and control rod withdrawal events with delayed and prompt supercritical reactivity insertions. The analysis in this report includes comparisons of the 2D thermal fluid porous media models in Pronghorn and SAM, and comparisons of coupled SAM/Griffin/SAM and coupled Pronghorn/Griffin for depressurized and pressurized loss of forced cooling, over-cooling, and control rod withdrawal events. In addition, we compare 3D, 2D, and 0D/PKE neutronic models for the two control rod withdrawal scenarios with coupled Pronghorn/Griffin. The comparisons show that the BlueCRAB models lead to physically intuitive solutions for the scenarios examined. The changes in the various scalar and vector fields, such as neutron flux, power, temperature, density, pressure, and velocity, are within the expected ranges, and their distributions can be explained by the system response of the transients and the geometric and material variations. Several comparisons suggest that the porous media models in Pronghorn and SAM can lead to similar solutions, even though they are based on different methodologies. This work further highlights the need for flexible tools with various levels of fidelity to cover the breadth and depth of needs that may arise in future technical evaluations of the PB-HTGR. We believe that the BlueCRAB capabilities will be a significant asset for confirmatory analyses in order to resolve important safety questions.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2023

Low-Temperature Plasma Assisted Kinetics Study of Ethanol

Abstract: Non- equilibrium plasma assisted kinetics of ethanol at 0.5 atm over a temperature range of 523-1203 K were studied at stoichiometric conditions by using a plasma flow reactor (PFR). The effects of plasma chemistry on fuel and its influence on neutral chemistry is decoupled with both experimental and modelling efforts. Enhanced reactivity was observed experimentally in both plasma-assisted oxidation and pyrolysis cases with fuel consumption beginning at temperatures as low as 523 K. Absolute fuel consumption was observed at temperatures nearly 200 K earlier than pure thermal reactions highlighting the effectiveness of non-equilibrium plasma to promote reactivity. Modelling results illustrated the effect of electron impact reactions towards accelerating fuel consumption and promoting low temperature chemistry by generating excited state N2 and reactive O and OH radicals.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2021

The role of gas flow distributions on CO 2 mineralization within monolithic cemented composites: coupled CFD-factorial design approach

The carbonation kinetics of monolithic cementing composites are strongly affected by gas transport which is, in turn, influenced by microstructural resistances and the presence of liquid water within pore networks. The non-uniform gas flow distribution within the CO 2 mineralization reactor can impart mass transfer resistance in the monolith microstructure, which affects the uptake of CO 2 (“carbonation”) of the cementing composites. This paper demonstrates how the gas spatial distribution (velocity and flow rate; quantified by CFD analysis) and processing conditions (temperature, relative humidity, and flow rate; quantified by factorial design) affect drying and carbonation, and in turn, the engineering properties of a representative ‘monolithic’ carbonate-cemented concrete component (i.e., herein concrete masonry unit: CMUs, also known as concrete block). It is shown that the gas flow distribution affects drying front penetration and results in moisture and carbonation gradients within the monolith. Particularly, variations in drying kinetics caused by non-uniformity of the contacting gas velocity impose gradients in moisture saturation, which results in increasing microstructural resistance to CO 2 transport. The resultant non-uniform carbonate-mineral formation (i.e., carbonate cementation), if not controlled, can produce gradients in mechanical properties and may alter failure patterns upon loading. Finally, these insights inform the optimal design of gas flow distribution systems and processing conditions within CO 2 mineralization reactors for the manufacturing of low-CO 2 concrete components using CO 2 -dilute industrial flue gas streams.

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

Experimental and computational investigations of ethane and ethylene kinetics with copper oxide particles for Chemical Looping Combustion

In this work, reaction pathways for the oxidation of methane, ethane, and ethylene with CuO was obtained by ReaxFF Molecular Dynamics (MD) simulations between temperatures of 1000 K and 2000 K. Experiments in a fixed-bed flow reactor were preformed with methane, ethane, and ethylene at temperatures ranging from 500 K to 1000 K with time-dependent species measurements from an Electron-Ionization Molecular Beam Mass Spectrometer (MBMS), and species validation with Gas Chromatography (GC) for detection of complete and intermediate combustion products. The MBMS and GC allow for the detection of oxygenated species and larger species produced from radical reformation. The simulation and experiment agree on the production of such species as CH 3 CHO, CH 2 O, CO, and H 2 O, which allow for the creation of simple C1 and C2 reaction pathways, which can be used in kinetic models of C2 species and larger fuels such as biofuels, which inherently depend on C1 and C2 kinetics and reaction pathway. The simulation and experiment disagree on the formation of C 2 H 2 , CH 3 OH, and CO 2 with large amounts of C 2 H 2 being measured in the ethylene oxidation simulations and CH 3 OH being formed in methane oxidation simulations, while neither species were experimentally found. In the case of CO 2 large amounts of CO 2 are rapidly produced in experiments with C2 fuels at 800 K, while little-to-no CO 2 was observed in simulations. This is believed to be resulting from the extremely short timescale of the simulations, preventing total oxidation of the fuel. Here, the differences in products produced between simulation and experiment allow for the potential to modify the ReaxFF potential functions to more accurately model the experimental products of Cu–H–O–C reaction kinetics.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

SAS4A/SASSYS-1 Commercial Grade Dedication Example Report for a Generic Sodium Pool-Type Fast Reactor Application

In the U.S., a key component of the commercialization of advanced reactors is completion of a license application, which must ultimately be approved by the Nuclear Regulatory Commission (NRC). The approval of the license application by the NRC is contingent on satisfactory demonstration of the design basis and the response of the advanced reactor design to transient and accident scenarios using accepted codes and methods. This report describes the qualification and dedication requirements that the advanced reactor safety analysis system software SAS4A/SASSYS-1 are expected to need to fulfill to be used for sodium-cooled pool-type fast reactor licensing. The qualification and dedication requirements are identified through performance critical characteristics and evaluation model acceptance criteria representative of the advanced reactor design considered for licensing. This document captures, additionally, the verification process developed to demonstrate that the software fulfills the qualification and dedication requirements for a generic sodium-cooled pool-type fast reactor as part of the commercial grade dedication process. Like most software that has primarily existed in the research and development space, the most significant challenge facing SAS4A/SASSYS-1 for use in a licensing framework is the availability of a documentation basis describing the code pedigree. SAS4A/SASSYS-1 has been used for licensing of the fast flux test facility (FFTF) and the JOYO sodium-cooled fast reactor in Japan, as well as the design of the CRBR Plant. However, the historical verification and validation (V&V) activities supporting SAS4A/SASSYS-1 development do not align with modern software quality assurance (SQA) and V&V requirements. Two approaches to use of SAS4A/SASSYS-1 in a commercial licensing framework have been identified: commercial-grade dedication (CGD) and software qualification. The methods and requirements prescribed in the ASME NQA-1-2008/2009 Standard and Regulatory Guide 1.203 on the evaluation model development and assessment process (EMDAP) have been used as guidance to define the CGD and qualification processes, respectively. A qualification and dedication requirements matrix has been developed which utilizes fundamental software verification. In this process, software verification is defined as a software quality process aimed at defining software requirement specifications, developing software design documentation, and performing and documenting acceptance testing of the code against requirements. A key element of software qualification and dedication includes determination of software acceptance with respect to critical characteristics relevant to the functional requirements of the software. To assist with identification of cross-cutting transient phenomena and functional requirements, domestic SFR vendor designs have been reviewed to identify a reference SFR design. For this report, the reference design is defined as a pool-type reactor with metal alloy fuel, a liquid-metal intermediate heat transport system, and passive decay heat rejection systems. Given this reference, a series of high-level cross-cutting phenomena was identified for a general class of single-fault undercooling or reactivity insertion transients that scopes the design basis space, with the goal of assisting with prioritization of documentation development efforts for key transient models in SAS: 1) Reactivity feedback response prior to scram; 2) System-wide thermal inertia; 3) Transition in natural circulation flow regime in heat removal systems; 4) Decay heat generation; 5) Steady-state fuel characterization; 5) Clad/fuel behavior at elevated temperatures; 6) Point kinetics and decay heat; 7) Pump coastdown behavior; 8) Core flow redistribution in loss of forced convection; 9) Pool stratification. As a demonstration of CGD of SAS4A/SASSYS-1 for a sodium pool reactor, a software qualification and dedication gap analysis as it relates to code documentation has been performed. This effort leverages the framework established as part of the SAS4A/SASSYS-1 SQA Program. This CGD demonstration provides a framework that vendors can build upon to demonstrate the applicability of the SAS4A/SASSYS-1 software for licensing a sodium-cooled pool-type fast reactor.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2020

Probing the low-temperature chemistry of methyl hexanoate: Insights from oxygenate intermediates

Understanding the combustion of methyl esters is crucial to elucidate kinetic pathways and predict combustion parameters, soot yields, and fuel performance of biodiesel, however most kinetic studies of methyl esters have focused on smaller, surrogate model esters. Methyl hexanoate is a larger methyl ester approaching the chain length of methyl esters found in biodiesel and has not received as much research attention as other smaller esters. The purpose of this work is to present the first atmospheric pressure combustion data of methyl hexanoate, CH 3 CH 2 CH 2 CH 2 CH 2 COOCH 3 . Mixtures of 2% methyl hexanoate in O 2 and N 2 are studied using a plug flow reactor at atmospheric pressure, wall temperatures from 573 to 973 K, residence times from roughly 1-2 s., and fuel equivalence ratios of 1, 1.5, and 2. Exhaust gases are analyzed by a gas chromatograph-mass spectrometer system and species mole fractions are presented. Here, the literature model shows satisfactory agreement with the experimental species profiles and improvements for future mechanistic studies are suggested. In particular, this work proposes new unimolecular decomposition pathways of methyl hexanoate to form methanol or methyl acetate. Furthermore, the experiment detected three unsaturated esters that are direct products of the low temperature oxidation chemistry and it provides more insight into branching ratios for the formation of methyl hexanoate radicals and for the decomposition of hydroperoxyalkyl radicals.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗