Analysis of Graphite Matrix Kinetics and Burnoff Products under Off-Normal High-Temperature Gas-Cooled Reactors Conditions
Explore the source record for details and available documents.
Brigham, Brian A.↗
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.
Catalytic hydrogenation of HMF to BHMF over copper catalysts
2,5-Bis(hydroxymethyl)furan (BHMF) is a bio-derived building block for polyester production, obtained via the hydrogenation of 5-hydroxymethylfurfural (HMF). First-principles thermodynamic equilibrium calculations indicate that this reaction is not thermodynamically limited under relevant conditions (e.g., 100 °C and high H 2 partial pressure). In this work, crude HMF was employed as the feedstock for BHMF synthesis. Initially, acidic impurities and humins were removed from unrefined HMF through filtration using a packed bed of γ-alumina. A comprehensive study of the filtration process is presented, including filtration kinetics, breakthrough curve analysis, and mathematical modeling. The purified HMF was subsequently hydrogenated over a 10 wt% CuZrO 2 catalyst, using ethanol as the reaction solvent. Batch reactions were first performed for collection of kinetic data to guide the transition to continuous flow operation. Kinetic data was collected in a fixed bed reactor at varying contact time, time on stream, temperature, and HMF concentration. This data was used to develop a kinetic model for HMF hydrogenation. Maximum BHMF production rates were achieved at 130 °C, accompanied by minor formation of byproducts from BHMF ring-opening reactions. The BHMF selectivity was 100 % at 100 °C although with lower reaction rates. Furthermore, catalyst stability tests revealed a loss of up to 50 % in catalytic activity within the first 24 h, likely due to the adsorption of HMF-derived oligomers that are not easily removed by filtration.
Crude HMF filtration↗
Unraveling the Ion Adsorption Kinetics in Microporous Carbon Electrodes: A Multiscale Quantum-Continuum Simulation and Experimental Approach
With the understanding that sorption in porous carbon electrodes is crucial to many environmental and energy technologies, such as capacitive deionization (CDI), supercapacitor energy storage, and activated carbon filters. In each of these examples, a practical model that can describe ion electrosorption kinetics is highly desirable for accelerating material design. Here, we proposed a multiscale model to study the ion electrosorption kinetics in porous carbon electrodes by combining quantum mechanical simulations with continuum approaches. Our model integrates the Butler-Volmer (BV) equation for sorption kinetics and a continuously stirred tank reactor (CSTR) formulation with atomistic calculations of ion hydration and ion-pore interactions based on density functional theory (DFT). We validated our model experimentally by using ion mixtures in a flow-through electrode CDI device and developed an in-line UV absorption system to provide unprecedented resolution of individual ions in the separation process. We showed that the multiscale model captures unexpected experimental phenomena that cannot be explained by the traditional ion electrosorption theory. The proposed multiscale framework provides a viable approach for modeling separation processes in systems where pore sizes and ion hydration effects strongly influence the sorption kinetics, which can be leveraged to explore possible strategies for improving carbon-based and, more broadly, pore-based technologies
36 MATERIALS SCIENCE↗
Importance of residence-time control of industrial screw-conveying reactors: Application to dilute-acid hydrolysis of biomass
Horizontal screw reactors are utilized in biorefineries for acid-catalyzed hydrolysis of xylan, which is a multi-step chemical reaction requiring accurate residence-time control. However, it is difficult to obtain online analytical measurement of reactant species. In this work, a residence-time distribution (RTD) is exhibited whose characteristics influence species yields. Sensitivity of product yield to RTD was investigated to understand the relative importance of operating control vs. inherent reactor dispersion. We find that reactor operation using a commonly used theoretical residence-time relationship can result in substantial yield losses. Instead, a model that accounts for the actual reactor RTD provides much improved results. The dispersion caused by reactor conditions only slightly hinders achieving theoretical optimal xylose yield (less than 3% yield loss for coefficient of variation less than 0.35), provided a validated RTD model is used to target the desired mean residence-time. In contrast, neglecting to account for the RTD by using the simplistic theoretical calculation results in xylose yields that are as much as 16% lower than the theoretical maximum.
09 BIOMASS FUELS↗
Predicting Morphological Evolution during Coprecipitation of MnCO 3 Battery Cathode Precursors Using Multiscale Simulations Aided by Targeted Synthesis
The performance of lithium-ion batteries is intimately linked to both the structure and the morphology of the cathode material, which in turn is critically linked to the synthesis conditions. However, few studies focus on understanding synthesis, especially during the coprecipitation of metal oxide precursors, a process that largely determines the final morphology of the material. In this paper, we go beyond the typical equilibrium particle shape analysis conducted in the literature and incorporate kinetic aspects of morphology evolution. We perform these studies using controlled synthesis on a well-defined metal salt system (MnCO 3 ) combined with multiscale simulations and high-resolution microscopy. Results show that with increasing metal concentration, the particles transition from rhombohedral to cubic to spherical shapes. Computational analysis using density functional theory (DFT) reveals that rhombohedral shaped particles evolve under equilibrium conditions. Phase field techniques indicate that at higher metal concentrations, fast growth kinetics of the precipitates result in the transition to cubic and, subsequently, spherical shapes, accompanied by a decrease in particle size. This study, while limited to the one metal salt system, provides an approach to shed light on the synthesis process of mixed transition metal salts, gradient materials, and other cathode materials of interest to the battery community.
25 ENERGY STORAGE↗
Initial exploration of a novel transient arrest system involving fuel heating
A preliminary analysis on a novel accident response system to diminish the severity of super- critical transients was conducted. The novel accident response system, called the instant shock arrest system, involves using electricity to heat the nuclear fuel at the onset of a large accidental reactivity insertion. This system is specifically designed for reactors with metallic fuel, such that the fuel is capable of conducting electricity, and being resistively heated. A reactor dynamics model of the advanced test reactor was created using the point kinetics equations and a linear reactivity feedback model to simulate how the system would effect the maximum fuel temperatures experienced during the transient. Transients with the instant shock arrest system were compared to those without it. It was found that the instant shock arrest system initially heated the fuel more than the unaffected transient but the negative reactivity inserted from such heating was enough to lower the maximum fuel temperature experienced during the transient. After simulating six different accident scenarios with reactivity insertions ranging from 0.5 to 1.3 dollar, it was found that an optimal system response could reduce peak fuel temperatures during the transient by 3.5% to 5%. Furthermore, discussion was given on how the optimal system response could be obtained using relatively simple numerical optimization algorithms due to the smoothness of the optimization problem. (authors)
22 GENERAL STUDIES OF NUCLEAR REACTORS↗
INITIAL EXPLORATION OF A NOVEL TRANSIENT ARREST SYSTEM INVOLVING FUEL HEATING
A preliminary analysis on a novel accident response system to diminish the severity of supercritical transients was conducted. The novel accident response system, called the instant shock arrest system, involves using electricity to heat the nuclear fuel at the onset of a large accidental reactivity insertion. This system is specifically designed for reactors with metallic fuel, such that the fuel is capable of conducting electricity, and being resistively heated. A reactor dynamics model of the advanced test reactor was created using the point kinetics equations and a linear reactivity feedback model to simulate how the system would effect the maximum fuel temperatures experienced during the transient. Transients with the instant shock arrest system were compare to those without it. It was found that the instant shock arrest system initially heated the fuel more than the unaffected transient but the negative reactivity inserted from such heating was enough to lower the maximum fuel temperature experienced during the transient. After simulating six different accident scenarios with reactivity insertions ranging from 0.5 \$ to 1.3 \$, it was found that the an optimal system response could reduce peak fuel temperatures during the transient by 3.5\% to 5\%. Furthermore, discussion was given on how the optimal system response could be obtained using relatively simple numerical optimization algorithms due to the smoothness of the optimization problem.
22 GENERAL STUDIES OF NUCLEAR REACTORS↗
INITIAL EXPLORATION OF A NOVEL TRANSIENT ARREST SYSTEM INVOLVING FUEL HEATING (Presentation)
A preliminary analysis on a novel accident response system to diminish the severity of supercritical transients was conducted. The novel accident response system, called the instant shock arrest system, involves using electricity to heat the nuclear fuel at the onset of a large accidental reactivity insertion. This system is specifically designed for reactors with metallic fuel, such that the fuel is capable of conducting electricity, and being resistively heated. A reactor dynamics model of the advanced test reactor was created using the point kinetics equations and a linear reactivity feedback model to simulate how the system would effect the maximum fuel temperatures experienced during the transient. Transients with the instant shock arrest system were compare to those without it. It was found that the instant shock arrest system initially heated the fuel more than the unaffected transient but the negative reactivity inserted from such heating was enough to lower the maximum fuel temperature experienced during the transient. After simulating six different accident scenarios with reactivity insertions ranging from 0.5 \$ to 1.3 \$, it was found that the an optimal system response could reduce peak fuel temperatures during the transient by 3.5% to 5%. Furthermore, discussion was given on how the optimal system response could be obtained using relatively simple numerical optimization algorithms due to the smoothness of the optimization problem.
22 GENERAL STUDIES OF NUCLEAR REACTORS↗
Kinetics of Radiation-induced Cr( II ) and Cr( III ) Redox Chemistry in Molten LiCl-KCl Eutectic
Chromium (Cr) is a frequent constituent of the metal alloys proposed for molten salt nuclear reactor (MSR) applications, and is typically the least noble metal ion present. Consequently, chromium is preferentially corroded into molten salt solutions. Here, the redox poise and redox cycling of chromium ions in the salt can greatly influence its corrosivity towards structural alloys, ultimately impacting the longevity of MSR systems. Radiation-induced chemistry is expected to play a significant role in determining the chromium oxidation state distribution during MSR operations. In the present research, electron pulse radiolysis techniques were employed to characterize the reactivity of Cr(II) and Cr(III) ions with primary radiolysis products in molten lithium chloride–potassium chloride (LiCl–KCl) eutectic over a temperature range of 400–600 °C. Both chromium oxidation states were found to rapidly react with the primary products of molten chloride salt radiolysis, i.e., the solvated electron (e S - ) and the dichlorine radical anion (Cl 2 ˙ - ). For reactions with the e S - , second-order rate coefficients (k) of k = (4.1 ± 0.2) and (6.1 ± 0.3) × 10 10 M -1 s -1 at 400 °C for Cr(II) and Cr(III), respectively, were determined. Temperature-dependent measurements allowed for the derivation of activation parameters for electron capture by Cr(II) and Cr(III). Both chromium ions also react with Cl 2 ˙ - , k = (7.2 ± 0.3) and (1.4 ± 0.1) × 10 9 M -1 s -1 at 400 °C for Cr(II) and Cr(III), respectively.
37 - INORGANIC, ORGANIC, PHYSICAL AND ANALYTICAL C↗
Transient kinetic insights into selective propene oxidation over industrial bismuth molybdate catalysts
Selective oxidation of propene to acrolein over industrial multicomponent bismuth molybdate (BMO) catalysts significantly depends on reaction conditions that include operating parameters and catalyst state. Here, this work investigates selective propene oxidation in the intrinsic kinetic paradigm of Temporal Analysis of Products (TAP) reactor by systematically varying catalyst redox state, temperature, and oxygen-to-propene feed ratio. A 93 % propene conversion with an acrolein yield of 80 % is achieved at elevated temperatures (450 °C) on oxidized catalysts under oxygen-rich conditions (O 2 :C 3 H 6 = 10). However, these conditions diminish acrolein-to-CO 2 selectivity due to enhanced total oxidation to CO 2 . In contrast, a reduced catalyst state, moderate temperature (350 °C), and lower oxygen feed (O 2 :C 3 H 6 = 1) nearly doubles the acrolein to-CO 2 selectivity, albeit at a lower acrolein yield (33 %). Transient kinetic studies together with a kinetic model that simplify the major product formation pathways in lumped non-elementary forms reveal that the availability of surface oxygen species plays a pivotal role in governing reaction pathways. Additionally, density functional theory (DFT) calculations on pure BMO catalysts inform the role of surface redox states on propene and oxygen activation barriers. Lattice oxygen at acrolein-selective sites drives both acrolein and CO 2 formation, while adsorbed oxygen at activation sites favors unselective CO 2 generation. This work establishes a critical relationship between transient product selectivity and surface oxygen availability, which is strongly influenced by catalyst redox state, feed ratio, and reaction temperature. These insights underscore the importance of dynamic reactor operation strategies and offer a foundation for designing next-generation propene oxidation processes with tunable acrolein selectivity.
42 - ENGINEERING↗
Application of an equation‐oriented framework to formulate and estimate parameters of chemical looping reaction models
Abstract Accurate, predictive reaction models are critical for the design and optimization of chemical looping combustion (CLC) reactors. The formulation and estimation of kinetic parameters for these reaction models using a first‐principles equation‐oriented (EO) approach is particularly beneficial as large amounts of experimental data spanning process‐relevant conditions can be used to estimate parameters in a computationally tractable way. This work demonstrates the application of a novel EO framework to develop reduction reaction kinetic models of an iron‐based CLC oxygen carrier (OC). An optimization problem is formulated to estimate kinetic parameters that provide the best fit to the experimental data. The model predicts the state of the OC with mean square error values of 2.5%–4.4% across the full range of validation data, including multiple reduction cycles.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Assessment of a Detailed Biomass Pyrolysis Kinetic Scheme in Multiscale Simulations of a Single-Particle Pyrolyzer and a Pilot-Scale Entrained Flow Pyrolyzer
A detailed biomass pyrolysis kinetic scheme was assessed in the multiscale simulations of a single-particle pyrolyzer with slow pyrolysis and a pilot-scale entrained flow pyrolyzer with fast pyrolysis. The detailed kinetic scheme of biomass pyrolysis developed by the CRECK group consists of 32 reactions and 58 species. A multiscale simulation model was developed, where the CRECK kinetics was employed to simulate biomass pyrolysis reactions, a one-dimensional particle model was utilized to simulate the intraparticle transport phenomena, and the particle-in-cell (PIC) model was employed to simulate the hydrodynamics. The multiscale model was first applied to simulate a single-particle pyrolysis experiment. The simulation with nonisothermal particles matched the experimental data better than the simulation with isothermal particles. Then the multiscale model was applied to simulate the pilot-scale entrained flow pyrolyzer. In this case, the simulation with isothermal particles matched the experimental data better than the simulation with nonisothermal particles. The reason for this difference might be that the kinetics itself already partially included the intraparticle transport effect as it was fitted using both TGA data (slow pyrolysis of small size biomass) and fluidized bed data (fast pyrolysis of relatively large size biomass). This study provides some insights into biomass pyrolysis kinetics development and pyrolyzer multiscale simulation for a future study.
09 BIOMASS FUELS↗
Comparison of Heterogeneous Hydroformylation of Ethylene and Propylene over RhCo 3 /MCM-41 Catalysts
The reaction mechanisms of heterogenous hydroformylation of ethylene and propylene were compared at 413-453 K using RhCo 3 /MCM-41 as catalysts. The reaction rates of propylene for both hydroformylation and the undesired side reaction of hydrogenation were found to be about one order of magnitude lower than those for ethylene in flow reactor studies. The difference in the kinetic behavior between ethylene and propylene was investigated by measuring the reaction orders and apparent activation energies, and these macro-kinetic observables were analyzed using the degree of rate control (DRC) method. In situ diffuse reflectance infrared Fourier transform spectroscopy (DRIFTS) experiments were performed to characterize the surface intermediates formed during the reactions. When the reactant was changed from ethylene to propylene, the IR peak corresponding to adsorbed CO exhibited a significant increase, while the IR peaks of alkyl group decreased in magnitude. Combined with the DRIFTS results, DRC analysis indicates that the first step of olefin hydroformylation, the formation of an alkyl group on the catalyst surface, plays a key role in the difference between ethylene and propylene. This step is kinetically non-relevant when ethylene is the reactant, but it is one of the rate-controlling steps for propylene. Furthermore, the low concentration of the adsorbed propyl group, which is a common intermediate shared by both hydroformylation and hydrogenation of propylene, decreases the rates of both reaction pathways as compared to ethylene.
10 SYNTHETIC FUELS↗
Flashback of Hydrogen-Methane Mixtures in a Fixed-Vane, Low-Swirl Burner at Atmospheric Pressure
Equivalence ratio at flashback, 𝜙FB, of low swirl injectors for 60% – 95% (H2, by mole) hydrogen-methane mixtures was determined experimentally in a < 100 kW atmospheric pressure burner. Pre-mixer velocities equaled 5, 7.5, and 10 m/s, with the reactants entering the burner at 294 K. The swirler design was selected following a literature review, and the test articles were additively manufactured. The test articles had two turning angles, 26° and 33°, with center-body hole diameters equaling 1.08 and 1.16 mm. A test article was built with each parameter combination, which yielded four test articles with measured swirl numbers varying from 0.43 to 0.49. The results showed agreement with prior work. 𝜙FB increased with increasing pre-mixer velocity and decreasing percent hydrogen. 𝜙FB was also observed to increase with swirl angle and decrease with increasing center-body hole diameter. The data was correlated using a semi-empirical quenching approach, with 𝐶quench=0.01, where the chemical time scale was calculated using a perfectly stirred reactor (PSR) model. Two different chemical kinetics mechanisms independently yielded the same results. The correlation makes the present small-scale laboratory results relevant to large, high-pressure gas turbines. The flashback data set is available for computational fluid dynamics (CFD) validation.
Searle, Matthew↗
The SAS4A/SASSYS-1 Version 5.8 Safety Analysis Code System
SAS4A/SASSYS-1 is a software simulation tool used to perform deterministic analysis of anticipated events as well as design basis and beyond design basis accidents for advanced nuclear reactors. Detailed, mechanistic models of steady-state and transient thermal, hydraulic, kinetic, and mechanical phenomena are employed to describe the response of the reactor core, the reactor primary and secondary coolant loops, the reactor control and protection systems, and the balance-of-plant to accidents caused by changes in coolant flow, loss of heat rejection, or reactivity insertion. The consequences of single and double-fault accidents can be modeled, including fuel and coolant heating, fuel and cladding mechanical behavior, core reactivity feedbacks, coolant loop performance including natural circulation, and decay heat removal. Analyses are typically terminated upon demonstration of reactor and plant shutdown to permanently coolable conditions, or upon violation of design basis margins. The objective of the analysis is to quantify accident consequences as measured by the transient behavior of system performance parameters, such as fuel and cladding temperatures, reactivity, and cladding strain. Originally developed for analysis of sodium cooled reactors with oxide fuel clad by stainless steel, the models were subsequently extended and specialized to metallic fuel clad with advanced alloys and to several other coolant options, including lead, LBE, and water.
22 GENERAL STUDIES OF NUCLEAR REACTORS↗
Computational Fluid Dynamics Combustion Modeling for Rotating Detonation Engines
This paper focuses on the development and validation of a combustion model for Computational Fluid Dynamics (CFD) modeling of Rotating Detonation Engines. A zero-dimensional Partially Stirred Reactor (PaSR) with a detailed chemical kinetic mechanism for hydrogen and air is used to model turbulent combustion. The model is computationally efficient and is based on the notion of partial mixing at the sub-grid level with turbulent exchange between mixed and unmixed regions. The ability of the PaSR model to accurately represent both detonative and deflagrative combustion is assessed by validating the results against experimental data. The effects of mesh resolution on the solution are also studied in order to determine if a mesh independent solution is obtainable with the Large Eddy Simulation (LES) approach to modeling turbulence. A comparison is made between the PaSR model and simply ignoring turbulence chemistry interactions which assumes that all species are perfectly mixed at the sub-grid level.
Strakey, Peter↗
Zeolite Membrane Reactor for Pre-Combustion Carbon Dioxide Capture (Final Scientific/Technical Report)
Water-gas shift (WGS) reaction followed by carbon dioxide (CO 2 ) separation is a critical step in the integrated gasification combined cycle (IGCC) process for fossil-fuel-fired electrical power generation with CO 2 capture. To intensify the IGCC process hydrogen-permselective zeolite membrane reactor offers promise to replace the conventional energy-intensive fixed-bed reactors and solvent-based CO 2 capture units. The objectives of this project were to develop a bench-scale zeolite membrane reactor (total membrane area: 932 cm 2 for a 21-tube membrane bundle) for the water-gas-shift reaction of raw syngas from an oxygen-blown coal-gasifier for H 2 production with simultaneous CO 2 separation at the capability of about 2 kilograms H 2 per day (equivalent to 2 kW IGCC power plant) and to demonstrate significant progress toward achieving overall performance goal of 90% CO 2 capture rate with 95% CO 2 purity at the cost of electricity 30% less than the baseline carbon capture approaches. This report summarizes results obtained in this project on scaling up the zeolite membrane reactor by a factor of 200 in membrane area, tests of the bench-scale zeolite membrane reactor for the water-gas-shift reaction at high temperature and high-pressure, and techno-economic analysis of the integration of the zeolite membrane reactor in IGCC power plant for the electrical generation with CO 2 capture. With effective pore modification by catalytic cracking deposition of MDES (methyldiethoxysilane), fabrication of MFI-type zeolite membranes was successfully scaled-up from a lab-scale disk type to bench-scale multiple-tube bundles on low-cost alumina supports. A Co-Mo based sour shift catalyst was evaluated and used in the zeolite membrane reactors for water-gas-shift reaction. The reaction kinetic and gas-permeation equations were developed and employed in mathematic models to guide/predict experiments/performance of zeolite membrane reactors for water-gas shift reaction. Multiple-tube zeolite membrane bundles and reactors were designed, fabricated and tested for gas separation and water-gas-shift reactions with real raw syngas from a coal-fired gasifier operated at high temperature and pressure. The zeolite membrane reactors demonstrated good long-term thermal and chemical stability and constant H 2 permeance (>300 GPU) together with a Knudsen selectivity (~4.7) for H 2 over CO 2 in the field test (cumulative time >28 hours) with high-sulfur coal-derived syngas. The zeolite-membrane-reactor integrated IGCC process was designed using the performance experimentally measured by the University of Cincinnati team with a single-tube zeolite membrane reactor that offers a CO conversion >98% with more than 90% of CO 2 and H 2 captured. With the above integrative approaches, a techno-economic analysis for a cost-benefit comparison was performed to uncover features that determine the power output, capital expenditure, operating expenditure, cost of electricity and cost of CO 2 capture in a 550-MW zeolite-membrane-reactor integrated IGCC process. The integration of the zeolite membrane reactor in IGCC could provide a significant reduction of 80% and 27% in the power consumption for Selexol™ Acid Gas Removal and CO 2 compression, respectively, which lowers the total auxiliary power consumption by 12.5%. However, the low pressure required at permeate stream for maintaining the driving force of hydrogen permeation through the zeolite membrane costs a huge power in permeate compressor, compensating the power consumption reduction achieved with the membrane reactor. Thus, for coal-fired IGCC for electricity generation with 90% CO 2 captured, the integration of the membrane reactor could provide a CO conversion ~99% and a significant drop in cost-of-electricity using zeolite membrane with H 2 permeance >600 GPU and the H 2 /CO 2 selectivity over 70.
01 COAL, LIGNITE, AND PEAT↗
Biomass pyrolysis devolatilization kinetics of herbaceous and woody feedstocks
Devolatilization kinetics were determined using a modified micropyrolyzer reactor for several biomass feedstocks: switchgrass, corn stover, red oak, and pine. Here, the micropyrolyzer was directly coupled to a flame ionization detector (FID) to track the release of volatiles from the biomass. Time series data from these experiments was analyzed to determine apparent devolatilization rates. Care was taken to assure the experiments were isothermal and kinetically limited calculating a Biot number less than 0.1 and pyrolysis numbers greater than 10, which simplifies the derivation of devolatilization rates. A single, first order reaction was able to model devolatilization rates at temperatures up to 500 °C. No correlation was found between the inorganic content of the biomass and its rate of devolatilization. Apparent activation energies were in the range of 54.9–88.4 kJ mol –1 . The rate coefficient at 500 °C was calculated as 1.90–5.14 s –1 for the four feedstocks.
09 BIOMASS FUELS↗