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 649 records · Page 36

Reacting CO 2 with Light Alkanes to Value-Added Products

Catalytic conversion of anthropogenic carbon dioxide (CO 2 ) into value-added products is a promising strategy to mitigate global carbon emissions. Concurrently, the shale gas revolution has provided an abundant supply of light alkanes (methane, ethane, propane, and butane), presenting a unique opportunity to employ these underutilized hydrocarbons as an effective, low-cost hydrogen source for CO 2 reduction. In this Perspective, we summarize past efforts, current state, and future opportunities for reacting CO 2 with light alkanes to generate a diverse range of value-added products. Compared with direct alkane conversion, the introduction of CO 2 fundamentally alters reaction thermodynamics and kinetics, enabling selective C–H and C–C bond activation while suppressing catalyst deactivation from coke formation. Building on decades of research in dry reforming and CO 2 -assisted dehydrogenation, recent advances in catalyst design have enabled CO 2 -assisted dehydrogenation processes that approach chemical equilibrium for the selective production of olefins and syngas. Importantly, advances in catalyst design and reactor engineering have further expanded the product scope beyond gas-phase (syngas and olefins) to include liquid-phase (oxygenates and aromatics), and solid-phase products (carbon nanomaterials). We highlight key catalyst design principles for controlling reaction pathways and discuss major challenges and opportunities in developing selective and versatile platforms for the simultaneous upgrading of CO 2 and light alkanes.

CO2↗

Development of an Advanced Computational Model for OMCVD of Indium Nitride

An advanced computational model is being developed to predict the formation of indium nitride (InN) film from the reaction of trimethylindium (In(CH3)3) with ammonia (NH3). The components are introduced into the reactor in the gas phase within a background of molecular nitrogen (N2). Organometallic chemical vapor deposition occurs on a heated sapphire surface. The model simulates heat and mass transport with gas and surface chemistry under steady state and pulsed conditions. The development and validation of an accurate model for the interactions between the diffusion of gas phase species and surface kinetics is essential to enable the regulation of the process in order to produce a low defect material. The validation of the model will be performed in concert with a NASA-North Carolina State University project.

Cardelino, Carlos A.↗

Insights into Constraining Rate Coefficients in Fuel Oxidation Mechanisms Using Genetic Algorithm Optimization

Accurate fuel oxidation mechanisms can enable predictive capabilities that aid in advancing combustion technologies. High-level computational kinetics can yield reasonable rate coefficients with uncertainties, in some cases, below a factor of 2. Computed rate coefficients can be constrained further by optimizing against experimental data. Here, we explore the application of genetic algorithm (GA) optimization to constrain computed rate coefficients in complex fuel oxidation mechanisms in conjunction with temperature-dependent species mole fractions from jet-stirred reactor (JSR) measurements. Cyclohexane is a model candidate for understanding the reactivity of cyclic fuels. In this work, we optimize the rate coefficients of the most recent literature cyclohexane mechanism, which incorporates theoretically computed rate coefficients for the reaction networks stemming from the first and second O 2 addition pathways, against the experimental results of two separate literature JSR studies. Optimization consistency is evaluated by carrying out three GA optimizations: fitting to the temperature-dependent species mole fractions in each JSR experiment separately and simultaneously fitting the species mole fractions in both experiments. Local sensitivity analyses are used to identify five influential low-temperature oxidation reactions for optimization. Although the three optimizations do not yield identical rate coefficients, the direction of change in all five rate coefficients is consistent among the three optimizations. Performance of the models from the three optimizations is assessed against literature ignition delay times with differences in the level of agreement observed among the different optimizations. In conclusion, comparisons are made with our recent optimization work of a cyclopentane oxidation master-equation model against time-resolved species concentrations, and insights and improvements of the strategy for constraining rate coefficients using GA optimization are discussed.

33 ADVANCED PROPULSION SYSTEMS↗

Determination of oxidation rates and volatile oxidation products for HTGR graphite matrix material exposed to steam atmospheres

High-temperature gas-cooled reactors (HTGRs) in operation use tristructural isotropic (TRISO) particles embedded in graphite and carbonized resin matrix to form the fuel element. This graphite matrix material serves as a supportive structural element, heat transfer medium, and neutron moderator. In HTGR designs, fuel compacts are exposed to helium coolant, which facilitates high outlet temperatures (750°C 2 O, 800°C 2 , and H 2 , are quantified in varied oxidant atmospheres using a coupled thermogravimetric analyzer and mass spectrometer. Furthermore, oxidation rates reported here for varied steam (H 2 O [g]) atmospheres are predominantly linear and comparable with literature values in the range of tested temperatures (800–1200°C). Changes in dominant matrix oxidation products from primarily CO to a mixture of CO, CO 2 , and H 2 were observed at higher temperatures (≥1000°C) and steam atmospheres (≥5 kPa pH 2 O). Kinetic data indicates that there was no shift in oxidation regime with chemical oxidation occurring at all temperatures and H2O (g) atmospheres tested. These data provide insight into the oxidation behavior of graphite matrix material and will inform future testing conditions, notably mixed atmospheric conditions, of HTGR fuel elements.

36 MATERIALS SCIENCE↗

Modeling of thermal and kinetic processes in non-equilibrium plasma ignition applied to a lean combustion engine

In recent years novel ignition systems have been developed to enable stable and efficient engine operations with lean mixtures. Among them, radio-frequency corona ignition systems create discharges that involve a much wider region compared to traditional spark, and produce non equilibrium plasma with high levels of active radicals and excited species. These devices considerably increase the early flame growth speed and extend stable operating limits. With the aim of expanding the knowledge on high efficiency lean-burn SI engines, this paper investigates and compares the combustion development generated by spark and corona ignitions through computational fluid dynamics, within the Reynolds-Averaged Navier-Stokes framework for turbulence modeling. In order to simultaneously take thermal and chemical effects into account, the Perfectly Stirred Reactor combustion model is used. Experimental data are also collected for validation in an optical access engine, for different mixture levels, from stoichiometric to very lean. Furthermore, the faster burn rate generated by the corona system in the initial stage of the combustion is well predicted by the simulations, in all the relative air-fuel ratio conditions. Remarkably, as the mixture becomes lean, simulations are able to capture the non-linear transition from fast to slow kernel growth, before a self-sustainable flame propagation is established. This correlates very well with the measured engine cyclic variability and the corresponding steep change in the duration of the flame kernel formation. Ultimately, this study highlights the important role of the atomic oxygen, as active radical, in promoting and enhancing the combustion initiated by a corona discharge, in addition to the volumetric ignition effect. By contrast, the validated simulations allow to explain that the high-temperature thermal plasma generated in a traditional spark discharge is insensitive to kinetic aspects.

42 ENGINEERING↗

Defect formation and bending properties in graphite under He atom implantation investigated by molecular dynamics method

The formation and evolution of radiation defects and their effects on bending property of a single graphite under low energy helium atom implantation have been studied at atomic scale through molecular dynamics simulations. Various simple interstitial and vacancy defects are observed during the implantation process, which is a function of incident energy and cross section of helium-carbon interaction. Large defects are observed after continuous helium atom implantations. The effects of incident energy on the density, size and distribution of defects are analyzed in detail. Furthermore, the kinetic evolution process of radiation defects at different annealing temperatures have also been simulated, including the migration, aggregation, transformation, coalescence, and self-healing. Here, the large defect clusters, crossing several carbon layers, have been confirmed to be the primary reasons that affect the bending property of graphite. All these results provide a new understanding for further applications of graphite in nuclear reactors.

Atomic implantation↗

Structural Aspects of Neutron Survival Probabilities

The neutron survival probability (and related quantities including probabilities of extinction and initiation) is a central element of the broader stochastic theory of neutron populations and finds application in fields including reactor start-up, analysis of reactor power bursts and criticality accidents, and safeguards. In a full neutron transport formulation, the equation governing the single-neutron survival probability is a backward or adjoint-like integro-partial differential equation with the added complexity of being highly nonlinear. Analogous formulations of this equation exist in the context of many approximate theories of neutron transport, with the point kinetics formulation having received significant theoretical attention since the 1940s. This work continues this tradition by providing a novel analysis of the single-neutron survival probability equation using the tools of boundary layer theory. The analysis reveals that the “fully dynamic” solution of the single-neutron survival probability equation—and some key probability distributions derived from it—may be cast as a singular perturbation around the underlying quasi-static single-neutron probability of initiation. In this perturbation solution, the expansion parameter is the ratio of the neutron generation time to a macroscopic time scale characterizing the overall system evolution; this interpretation illuminates some of the fundamental structural aspects of neutron survival phenomena.

97 MATHEMATICS AND COMPUTING↗

Interactive computer modeling of combustion chemistry and coalescence-dispersion modeling of turbulent combustion

An interactive computer code for simulation of a high-intensity turbulent combustor as a single point inhomogeneous stirred reactor was developed from an existing batch processing computer code CDPSR. The interactive CDPSR code was used as a guide for interpretation and direction of DOE-sponsored companion experiments utilizing Xenon tracer with optical laser diagnostic techniques to experimentally determine the appropriate mixing frequency, and for validation of CDPSR as a mixing-chemistry model for a laboratory jet-stirred reactor. The coalescence-dispersion model for finite rate mixing was incorporated into an existing interactive code AVCO-MARK I, to enable simulation of a combustor as a modular array of stirred flow and plug flow elements, each having a prescribed finite mixing frequency, or axial distribution of mixing frequency, as appropriate. Further increase the speed and reliability of the batch kinetics integrator code CREKID was increased by rewriting in vectorized form for execution on a vector or parallel processor, and by incorporating numerical techniques which enhance execution speed by permitting specification of a very low accuracy tolerance.

Pratt, D. T.↗

Testing to Evaluate Processes Expected to Occur during MSR Salt Spill Accidents

Obtaining a license for a new nuclear reactor requires the identification and assessment of the potential consequences of specified accident scenarios, which are achieved using accident progression models. Those models need to be parameterized and validated using experimental data, but existing experimental data addressing processes relevant to molten salt reactor accidents are sparse. Specifically, experimental data that quantify the sensitives of important processes to the initial conditions of the spill, the ambient environment, and the containment features are needed to parameterize individual process models. Integrated experiments that simulate accident scenarios are also needed to provide data for model validation, but these experiments will require the use of proven methods to quantify the processes under evaluation. The overarching objectives of this work are to develop the methods for simulating the targeted processes, to determine the effectiveness of the methods in producing the data required for model development, to generate data that can be used to parameterize individual process models, and to provide key insights into the behavior of spilled molten salt that should be considered in models. Experimental methods were designed to quantify aspects of individual processes expected to occur during or after a molten salt spill accident that will affect the fate of spilled molten salt and the radionuclides within. These processes include 1) molten salt spreading and heat transfer, 2) molten salt flowing and freezing in tubing, 3) stainless steel corrosion kinetics in molten salt, and 4) molten salt splashing and aerosol generation. The initial tests described in this report were conducted using eutectic FLiNaK to demonstrate the test methods, the data that are generated, and the analyses of the data to derive values needed for modeling. The primary variables that were tested include initial salt temperature and the presence of volatile surrogate fission products (e.g., cesium and iodine). The developed methods are shown to be effective in quantifying the desired processes and can be applied to study more complex salt compositions of interest to molten salt reactor developers, a wide range of environmental conditions of interest to modelers, and additional variables relevant to salt spill accidents. The developed methods and insights gained from laboratory tests can also be incorporated in future large-scale integrated tests used to simulate molten salt spill accidents.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Mechanistic modeling of lifetime distribution of SiC/SiC composite claddings

Silicon carbide (SiC) fiber-reinforced SiC matrix (SiC/SiC) composites have emerged as a new material candidate for fuel claddings in light water reactors. Recent studies showed that the load capacity of SiC/SiC materials exhibits a considerable statistical variation. Therefore, reliability analysis plays a critical role in design of SiC/SiC composite claddings. Here, this paper presents a probabilistic model for the lifetime distribution of SiC/SiC composites. The model is anchored by a multiaxial stress-based failure criterion and subcritical damage accumulation mechanism. Based on the kinetics of subcritical damage growth, the lifetime distribution of a laboratory test specimen for any given loading history can be calculated. A finite weakest-link model is used to extrapolate the lifetime distribution of test specimens to full-length claddings. It is shown that the damage accumulation mechanism has a strong influence on the lifetime distribution of the cladding. This finding highlights the importance of understanding the static fatigue behavior of SiC/SiC composites. The present analysis also demonstrates an intricate length effect on the failure probability of the cladding, which is expected to play a crucial role in design extrapolation.

36 MATERIALS SCIENCE↗

Air oxidation of yttrium hydride as a high temperature moderator for thermal neutron spectrum fission reactors

Yttrium hydride (YH x ) is an attractive moderator material for thermal neutron spectrum fission reactors requiring a small reactor core volume and has been selected as the neutron moderator for the Transformational Challenge Reactor (TCR), an advanced gas-cooled microreactor. Before YH x can be used in this application, it is important to understand the material response to off-normal conditions. In the present study, 550–650 °C isothermal dry air oxidation was performed to simulate a depressurized loss of force circulation (DLOFC) event. The oxidation was performed using thermogravimetric analysis (TGA) on bulk crack-free YHx coupons. Oxidation studies were also performed on Y coupons to elucidate the impact of H on oxidation. Both the chemistry and distribution of processing impurities were found to strongly affect oxidization behavior on a batch-to-batch basis. Regardless of batch, YHx oxidized at a significantly lower rate than Y at all temperatures, and the lower rate was directly correlated with increased hydride content. Metallic Y exhibited complex exponential kinetics, whereas YH x also exhibited complex kinetics but gained considerably less mass. According to literature reports on protonic and native-ion conductivities of Y 2 O 3 and mass spectrometry analysis of gaseous reaction products formed during the oxidation of YH x , a mechanism for the reduced oxidation rate of yttrium hydride is suggested.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Mass transport in nanoarray monolithic catalysts: An experimental-theory study

Reducing the mass transfer resistance globally of a catalyst is a key to enhancing the catalytic reaction kinetics and fully utilizing the catalyst activity. Despite the success in tailoring the external mass transfer in the widely studied washcoat monoliths, the internal mass transfer resistance is difficult to be reduced due to the requirement of increasing macroporosity while maintaining high specific surface area and mechanical stability. Therefore, nanostructured array-based monolithic catalysts (nanoarray catalysts) have been developed in the past decade as a promising class of structured catalysts that may complement or substitute washcoat catalysts. This work fundamentally elucidates the enhanced mass transport properties of the nanoarray monolithic catalysts by a combination of experimental measurements and theoretical modeling. Using a low-dimensional model, the relative contributions of resistances were quantified in terms of chemical kinetics, internal and external mass transfers based on a probe model of C2H4 oxidation over the TiO2 supported Pt-based monolithic catalysts. The nanoarray catalysts displayed a lower internal mass transfer resistance than the washcoat counterparts as a result of the high macroporosity and small thickness of nanoarray layers. Finally, the nanoarray configuration provides a new pathway towards designing high-performance monolithic reactors and catalysts with low internal diffusion limitations for various gas phase reactions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

PAH formation from jet stirred reactor pyrolysis of gasoline surrogates

Soot particles and their precursor polycyclic aromatic hydrocarbon (PAH) species, formed during combustion, are responsible for particulate emissions in gasoline direct injection (GDI) engines. To better understand the effects of fuel composition on formation of soot in GDI engines, the pyrolysis of several gasoline surrogates was studied in a jet-stirred reactor across a broad temperature range at atmospheric pressure and 1 s residence time. Fuel and intermediate species, including aromatics up to naphthalene, were measured using gas chromatography (GC). PAH concentrations from pyrolysis of surrogate fuels were compared to gain insight into the effects of fuel composition on PAH formation. In addition, synergistic effects were observed in pyrolysis experiments of binary blends. A detailed kinetic model, recently developed at Lawrence Livermore National Laboratory (LLNL), successfully captured the effects of blending and the concentration of major PAHs. Major reaction pathways are discussed, as well as the role of important intermediate species, such as acetylene, and resonantly stabilized radicals such as allyl, propargyl, cyclopentadienyl, and benzyl in the formation of PAH.

36 MATERIALS SCIENCE↗

Structured light approaches in laser-based plasma diagnostics

There is a growing demand for plasma diagnostics suitable for industrial plasma reactors employed in semiconductor nanofabrication, especially relevant to microelectronics and quantum information systems. Such reactors typically have limited optical access and pose considerable diagnostic challenges, including intense background emission, significant thermal loads, and contamination of optical viewports. In this study, we outline research into structured light techniques (laser beams with tailored spatial, temporal, or phase characteristics) that effectively overcome these issues using laser-induced fluorescence (LIF) as an example. The focus of presented diagnostics is on ion kinetics analysis within an industrial plasma source, although this approach is broadly applicable to other plasma systems and diagnostic contexts. We present a confocal LIF implementation using an axicon-generated Bessel annular beam, achieving spatial resolutions of approximately 5 mm at a focal distance of 300 mm, with potential improvements to about 1 mm. This approach matches conventional orthogonal LIF performance but requires only one optical port. Wavelength-modulation LIF employs nonlinear laser wavelength tuning to measure spectral line derivatives, suppressing background emission and enhancing details of spectral line shape. Additionally, we present new results on applying vortex beams (laser beams carrying orbital angular momentum, OAM) for LIF measurements in an industrial plasma device. These measurements enable simultaneous axial and tangential velocity determination using a single laser beam and have been tested with xenon ion transition. Initial quantification of results was performed. Together, these structured-light approaches provide robust, background-resilient, multi-dimensional diagnostics for complex plasma environments.

Romadanov, Ivan [Princeton Plasma Physics Laborato↗

Solar Decarbonization of Paraffin Dehydrogenation Through Particle Heat Carriers (Final Technical Report)

This project focuses on solutions to decarbonize high-temperature catalytic processes using solar thermal heat. The primary project goal is to show the validity of a moving packed bed reactor for propane dehydrogenation using catalyst particles as the heat carrier for the reaction, which can be heated by concentrated solar energy in a particle receiver. This concept, if further developed, may provide a cost-effective pathway for converting lower value gases to important chemical precursors for industrial materials using only renewable energy. The project was divided into six tasks. In Task 1, DFT calculations were performed to understand the role of Pt and Sn in the catalytic dehydrogenation reaction. In Task 2 chemical kinetics measurements were made for several catalyst formulations at high temperatures. In Task 3, the solar absorptance of catalyst particles was compared to the absorptance of commonly used materials in particle receivers. In Task 4, numerical models were developed which could predict performance of the complete system and predict specific temperatures in the system. In Task 5, a prototype system was designed, fabricated, and tested to show the validity of the concept. Task 6 concerned project management activities. Experiments with the prototype showed repeatable thermal performance at temperatures targeted for the reaction. A limited set of tests were done with active catalyst and propane dehydrogenation, showing conversion of propane to propylene with a range of conversions and selectivities. The results are promising, and the prototype designed was reliable during testing, and the team expects that further development of the prototype would yield improved results. A numerical model framework based on coupled fluid and particle mechanics was developed with high computational efficiency using GPU calculations. The model may prove highly useful for evaluating other high-temperature particle systems. However, it was determined that simpler porous media models were good fits for the needs of the current moving packed bed concept. Data showing strong solar absorption of the particles validates the plan of using existing solar particle receivers with only a change in the particle type. Catalyst investigation showed that Pt 1 Sn 3 is the most viable candidate for developing PtSn catalysts for high temperature propane dehydrogenation, considering the balance of activity, selectivity, and deactivation. This project completed an initial study of various factors needed to incorporate a moving bed catalytic reactor for propane dehydrogenation into a concentrated solar thermal particle system. Future developments may allow this technology to be scaled up and help to use solar thermal energy to decarbonize not only the propane dehydrogenation reaction, but other gas-solid catalytic reactions at similar temperatures.

14 SOLAR ENERGY↗

Discovery and Characterization of a Metastable Cubic Interstitial Nickel–Carbon System with an Expanded Lattice

Metastable, i.e., kinetically favored but thermodynamically not stable, interstitial solid solutions of carbon in iron are well-understood. Carbon can occupy the interstitial atoms of the host metal, altering its properties. Alloying of the host metal results in the stabilization of the FeC x phases, widening its application. Pure nickel finds niche applications, mainly focusing on catalysis, while nickel alloys are widely applied, e.g., in gas turbines, reactors, and seawater piping. Nickel carbide (Ni 3 C) is the well-known stable Ni–C system displaying a rhombohedral (R3̅c) crystal structure. Some reports describe an elusive cubic Ni–C system, observed during certain catalytic reactions occurring on nickel and formed by the occupation of the interstitials of the metal with carbon: to date, the stabilization and characterization of this phase have not been accomplished. Hereby, we report on the synthesis of a cubic metastable NiC x phase using chemical vapor deposition of methane on supported nickel nanoparticles. The structure was predicted by DFT/ReaxFF, synthesized and monitored with in situ time-resolved synchrotron XRD, and experimentally confirmed by Rietveld refinement and (S)TEM-EELS under ambient conditions. The results show an Fm3̅m phase with a lattice parameter of a = 3.749 ± 0.037 Å at room temperature, with the highest ever reported atomic percentage of carbon occupying the octahedral interstices of 23.1%, resulting in a NiC 0.3 phase. The degree of occupation of the interstitial voids by carbon can be controlled, enabling the tuning of the host metal’s d-spacing and composition, highlighting the applicability of this synthesis route for catalytic nanoparticle preparation.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Modelling ignition characteristics of rich H2/O2 mixture in a monolithic catalytic reactor

Ignition of rich hydrogen-oxygen mixture in a monolithic catalytic reactor is studied using a transient combustion model. The model assumes a quasi-steady gas phase and a thermally-thin substrate. The ignition time lag is due to the thermal inertia of the substrate. One-step global chemical reaction is assumed both on the surface and in the gas phase. An effectiveness factor is introduced to account for the transition from a kinetically-limited catalytic surface reaction to a diffusion-limit one during ignition. Results presented include a catalytic ignition boundary, ignition delay time and the transient response in the catalytic bed.

Tien, Ta-Ching↗

Updated Chemical Kinetics and Sensitivity Analysis Code

An updated version of the General Chemical Kinetics and Sensitivity Analysis (LSENS) computer code has become available. A prior version of LSENS was described in "Program Helps to Determine Chemical-Reaction Mechanisms" (LEW-15758), NASA Tech Briefs, Vol. 19, No. 5 (May 1995), page 66. To recapitulate: LSENS solves complex, homogeneous, gas-phase, chemical-kinetics problems (e.g., combustion of fuels) that are represented by sets of many coupled, nonlinear, first-order ordinary differential equations. LSENS has been designed for flexibility, convenience, and computational efficiency. The present version of LSENS incorporates mathematical models for (1) a static system; (2) steady, one-dimensional inviscid flow; (3) reaction behind an incident shock wave, including boundary layer correction; (4) a perfectly stirred reactor; and (5) a perfectly stirred reactor followed by a plug-flow reactor. In addition, LSENS can compute equilibrium properties for the following assigned states: enthalpy and pressure, temperature and pressure, internal energy and volume, and temperature and volume. For static and one-dimensional-flow problems, including those behind an incident shock wave and following a perfectly stirred reactor calculation, LSENS can compute sensitivity coefficients of dependent variables and their derivatives, with respect to the initial values of dependent variables and/or the rate-coefficient parameters of the chemical reactions.

Radhakrishnan, Krishnan↗