Search NASA⌕ Search

SEARCH · Search NASA

Results for “Kinetic Modelling”

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 145 records · Page 8

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↗

Probabilistic and maximum entropy modeling of chemical reaction systems: Characteristics and comparisons to mass action kinetic models

We demonstrate and characterize a first-principles approach to modeling the mass action dynamics of metabolism. Starting from a basic definition of entropy expressed as a multinomial probability density using Boltzmann probabilities with standard chemical potentials, we derive and compare the free energy dissipation and the entropy production rates. We express the relation between entropy production and the chemical master equation for modeling metabolism, which unifies chemical kinetics and chemical thermodynamics. Because prediction uncertainty with respect to parameter variability is frequently a concern with mass action models utilizing rate constants, we compare and contrast the maximum entropy model, which has its own set of rate parameters, to a population of standard mass action models in which the rate constants are randomly chosen. We show that a maximum entropy model is characterized by a high probability of free energy dissipation rate and likewise entropy production rate, relative to other models. We then characterize the variability of the maximum entropy model predictions with respect to uncertainties in parameters (standard free energies of formation) and with respect to ionic strengths typically found in a cell.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Identifying factors controlling the selective ethane dehydrogenation on Pt-based catalysts from DFT based micro-kinetic modeling

Herein, a detailed EDH reaction mechanism on Pt(111), Pt(211), Pt 3 Sn(111) and Pt 3 Sn(211) surfaces has been systematically computed with DFT. A mean-field micro-kinetic model has been established to simulate the product distribution and selectivity on each surface at experimentally relevant conditions. We have uncovered a number of different factors controlling the catalytic performance of Pt-based EDH catalysts, thus leading to a better understanding of experimentally detected phenomenon. In this study, the (111) and (211) model surfaces of Pt and Pt 3 Sn catalysts were chosen to represent the terrace and step-edge sites when simulating the EDH reaction mechanisms.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Fermi-Based Kinetic Model for the Sabatier Reaction: Sabatier Principle and Beyond It

The CO 2 methanation reaction (aka the Sabatier reaction) is one of the important carbon utilization technologies to reduce atmospheric greenhouse gases. In this study, kinetic modeling of the Sabatier reaction proceeding via the CO route has been presented by considering weak and strong chemisorption concepts rooted in electronic theory. Based on the single-route mechanism of the Sabatier reaction, the dependence of the reaction rate on the position of the catalyst Fermi level with respect to the antibonding states of the reagent (CO 2 ) and the product (CH 4 ), is demonstrated with the classic Sabatier volcano-shape. Furthermore, the model was applied to explain experimental results of reaction rates for Ni- and Ru-based catalytic systems and to understand origin of enhancement of turnover frequencies established for dynamic catalysis.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Phase transformation kinetics model for metals

We develop a new model for phase transformation kinetics in metals by generalizing the Levitas–Preston (LP) phase field model of martensite phase transformations (see Levitas and Preston (2002a,b) and Levitas et al. (2003)) to arbitrary pressure. Furthermore, we account for and track: the interface speed of the pressure-driven phase transformation, properties of critical nuclei, as well as nucleation at grain sites and on dislocations and homogeneous nucleation. The volume fraction evolution of each phase is described by employing KJMA (Kolmogorov, 1937; Johnson and Mehl, 1939; Avrami, 1939, 1940, 1941) kinetic theory. We then test our new model for iron under ramp loading conditions and compare our predictions for the α → ϵ iron phase transition to experimental data of Smith et al. (2013). In conclusion, more than one combination of material and model parameters (such as dislocation density and interface speed) led to good agreement of our simulations to the experimental data, thus highlighting the importance of having accurate microstructure data for the sample under consideration.

36 MATERIALS SCIENCE↗

Encapsulation of highly viscous CO 2 capture solvents for enhanced capture kinetics: Modeling investigation of mass transfer mechanisms

The encapsulation of highly viscous liquid-like Nanoparticle Organic Hybrid Materials (NOHMs) inside a gas permeable polymer to form SIPs (Solvent Impregnated Polymers) significantly enhanced the CO 2 capture kinetics of NOHMs, leading to a remarkable 50-fold increase in CO 2 flux compared to the neat NOHMs. To understand the mechanism for enhanced CO 2 mass transfer within these hybrid materials, kinetic modeling of CO 2 uptake into SIPs containing polyethylenimine functionalized NOHMs, denoted NPEI-SIP, was conducted. CO 2 mass transfer into NPEI-SIP films was found to conform, both qualitatively and quantitatively, with a diffusion-controlled moving front model. The diffusion-controlled model was also used to simulate CO 2 uptake within a fixed bed containing polydisperse NPEI-SIP particles, and this model accurately predicted experimentally measured breakthrough curves at 25 °C and 50 °C. The 50-fold increase in gas flux was shown to be a consequence of the very large CO 2 permeability within the polymer–solvent composite. Furthermore, the increase in gas flux is also dependent on the diffusion–reaction regime in which the chemical solvent operates, and the largest improvement will occur when immobilizing solvents, such as NOHMs, which operate in the instantaneous-reaction regime. The CO 2 capacity (~3 mol CO 2 /kg) and saturation time (~5 min) of 430 μm SIP particles were comparable to popular CO 2 chemisorption materials such as amine grafted silicates, in spite of the slow kinetics of NOHM-I-PEI in liquid form (CO 2 saturation time ~24 h for a 1 mm thin film).

36 MATERIALS SCIENCE↗

A Kinetic Model of the Long-Term Corrosion of Glass-Ceramic Materials

Multiphase waste forms show promise for increased waste loading and for the ability to dispose of contaminated solid and particulate waste through direct densification. However, achieving predictive capability for long-term durability of multiphase waste forms, and thus assessing their possible deployment, requires expanding the current, limited knowledge base. Here, we describe the development of a corrosion model of a two-phase waste form consisting of crystals of known volume fraction embedded in a glass matrix. This model accounts for the dissolution of both the crystalline and glass phases as well as the hydration of the glass phase through an ion exchange reaction. Because of the large difference in solubility between the two phases, the reactive surface of the crystalline phase is a function of the extent of dissolution of the glass phase in this model. Model parameterization was performed using corrosion data, such as from single-pass flow-through tests, for the individual phases. The parameterized corrosion model was evaluated against static dissolution test data for a glass-ceramic multiphase waste form. This evaluation demonstrated the model’s ability to reproduce the time-dependent release of key tracers of glass and crystalline phase dissolution. Hence, the development of a kinetic model provides a pathway for long-term durability predictions and thus the use of multiphase waste forms in nuclear cleanup missions.

Kerisit, Sebastien N.↗

Capability in Theory, Modeling, and Validation for a Range of Innovative Fusion Concepts using High-Fidelity Moment-Kinetic Models

A computational modeling capability is created and available to the fusion community to understand and design lower-cost and innovative fusion concepts. The approach uses high- fidelity kinetic, moment-kinetic, and moment models and includes sophisticated plasma- boundary interactions. A majority of fusion-relevant simulations are performed with magnetohydrodynamic models and hybrid particle-in-cell codes, with limited-fidelity electron and kinetic physics. However, in fusion configurations like Z-pinches, field-reversed- configurations, plasma jet magneto-inertial fusion, spinning mirrors, and others, kinetic effects (both electron and ions) are critical to understand the physics and design scaling into the highly kinetic regime of a burning fusion plasma. Furthermore, as present fusion machines move towards a burning plasma regime, liquid-metal blankets are needed to handle first-wall heat- flux, reduce erosion, and eventually for energy conversion and fuel breeding. The work performed under this ARPA-E BETHE Capability Team advances the state-of-the-art in modeling and understanding plasma dynamics in fusion devices and its coupling with liquid-metal dynamics. These are critical areas of research for fusion energy to become realizable. To address these complex problems, we have leveraged and extended computational capabilities through the code, Gkeyll (developed jointly with Princeton Plasma Physics Laboratory and academic partners), for kinetic and moment modeling of fusion plasmas. The Concept Teams supported by this Capability Team include the Wisconsin High-field Axisymmetric Mirror (WHAM), Centrifugal Mirror Experiment (CFME), Plasma-Jet Magneto- Inertial Fusion (PJMIF), and solid and liquid wall plasma-material interaction studies relevant to a number of fusion concepts including Zap Energy’s Z-pinch. This software is open-source and available to the fusion community as a high-fidelity tool for the design of lower-cost fusion experiments. 3D gyrokinetic simulations of WHAM are now possible for long enough time scales to understand the evolution of interchange instabilities. 3D multi-fluid simulations of CMFE at higher Mach numbers are now possible for detailed design iterations with the goal of stability. The state-of-the-art in understanding shock formation and shock mitigation regimes in merging liners for PJMIF have been furthered by our kinetic simulations. Our novel models and frameworks studying plasma-material interaction by incorporating wall emission for various solid wall materials of relevance to pulsed and steady fusion concepts have advanced the state-of-the-art in our understanding of particle fluxes, heat fluxes, and other quantities at cathodes and anodes. The results from this work may explain discrepancies between experimental and theoretical predictions of achieved current densities in pulsed concepts such as Z-pinches. Another significant contribution of this Capability Team is the development and deployment of a novel experimental platform, LEX (Liquid Electrode eXperiment), at Virginia Tech to understand liquid metal free-surface response to electromagnetic pulses. The novel experiments along with model validation quantified the effect of different materials and sizes of liquid metal droplets on the radiative power balance of fusion plasmas for pulsed concepts. Furthermore, these experiments provided mitigation strategies for violent liquid metal response for high current pulses as would be expected in fusion regimes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Dynamic Model of Basic Oxygen Steelmaking Process Based on Multizone Reaction Kinetics: Modeling of Manganese Removal

In the earlier study, a dynamic model for the BOF process based on the multizone reaction kinetics has been developed. In the preceding part, the mechanism of manganese transfer in three reactive zones of the converter has been analyzed. The study predicts that temperature at the slag–metal reaction interface plays a major role in the Mn reaction kinetics. Further, mathematical treatments to simulate the transient rate parameters associated with each reaction interface have been developed. The model calculations of Mn removal rate obtained from different zones of the converter predicts that the first stage of the blow is dominated by the oxidation of Mn at the jet impact zone, albeit some additional Mn refining has been observed as a result of the oxidation of metal droplets in emulsion phase. The simulation result shows that the reversion of Mn from slag to metal primarily takes place at the metal droplet in the emulsion due to an increase in slag–metal interface temperature during the middle stage of blowing. In the final stage of the blow, the competition between simultaneous reactions in jet impact and emulsion zone controls the direction of mass flow of manganese. Further, the model prediction shows that the Mn refining in the emulsion is a strong function of droplet diameter and the residence time. Smaller sized droplets approach equilibrium quickly and thus contribute to a significant Mn conversion between slag and metal compared to the larger sized ones. The overall model prediction for Mn in the hot metal has been found to be in good agreement with two data sets obtained from different size converters reported in the literature.

36 MATERIALS SCIENCE↗

Experimental Validation of the Plasma Kinetic Model

Electrostatic discharge (ESD) events pose serious safety concerns for nuclear weapon handling and maintenance. Current safety analyses operate under overly conservative assumptions resulting in costly safety procedures to inadvertent discharges. Colorado School of Mines has developed a plasma kinetic model coupled with a Maxwell’s equations solver in COMSOL to address these assumptions. A spherical capacitor system was designed to measure current throughout a threshold breakdown. The results will be compared to the simulations run by CSM. This paper details the design process and preliminary results for this validation effort.

42 ENGINEERING↗

Interfacial mixing in high-energy-density matter with a multiphysics kinetic model

We have extended here a recently developed multispecies, multitemperature Bhatnagar-Gross-Krook model [Haack et al., J. Stat. Phys. 168, 822 (2017)], to include multiphysics capabilities that enable modeling of a wider range of physical conditions. In terms of geometry, we have extended from the spatially homogeneous setting to one spatial dimension. In terms of the physics, we have included an atomic ionization model, accurate collision physics across coupling regimes, self-consistent electric fields, and degeneracy in the electronic screening. We apply the model to a warm dense matter scenario in which the ablator-fuel interface of an inertial confinement fusion target is heated, but for larger length and time scales and for much higher temperatures than can be simulated using molecular dynamics. Relative to molecular dynamics, the kinetic model greatly extends the temperature regime and the spatiotemporal scales over which we are able to model. In our numerical results we observe hydrogen from the ablator material jetting into the fuel during the early stages of the implosion and compare the relative size of various diffusion components (Fickean diffusion, electrodiffusion, and barodiffusion) that drive this process. We also examine kinetic effects, such as anisotropic distributions and velocity separation, in order to determine when this problem can be described with a hydrodynamic model.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Dynamic Model of Basic Oxygen Steelmaking Process Based on Multizone Reaction Kinetics: Modeling of Decarburization

In a previous study by the authors (Rout et al. in Metall Mater Trans B 49:537–557, 2018), a dynamic model for the BOF, employing the concept of multizone kinetics was developed. In the current study, the kinetics of decarburization reaction is investigated. The jet impact and slag–metal emulsion zones were identified to be primary zones for carbon oxidation. The dynamic parameters in the rate equation of decarburization such as residence time of metal drops in the emulsion, interfacial area evolution, initial size, and the effects of surface-active oxides have been included in the kinetic rate equation of the metal droplet. A modified mass-transfer coefficient based on the ideal Langmuir adsorption equilibrium has been proposed to take into account the surface blockage effects of SiO{sub 2} and P{sub 2}O{sub 5} in slag on the decarburization kinetics of a metal droplet in the emulsion. Further, a size distribution function has been included in the rate equation to evaluate the effect of droplet size on reaction kinetics. The mathematical simulation indicates that decarburization of the droplet in the emulsion is a strong function of the initial size and residence time. A modified droplet generation rate proposed previously by the authors has been used to estimate the total decarburization rate by slag–metal emulsion. The model’s prediction shows that about 76 pct of total carbon is removed by reactions in the emulsion, and the remaining is removed by reactions at the jet impact zone. The predicted bath carbon by the model has been found to be in good agreement with the industrially measured data.

36 MATERIALS SCIENCE↗

A Kinetic Model for Modification of MgAl{sub 2}O{sub 4} Spinel Inclusions During Calcium Treatment in the Ladle Furnace

A common approach in the processing of Al-killed steels is the modification of alumina inclusions by calcium treatment. A model for this process has already been developed by the current authors and validated by means of plant data. Often, magnesium aluminate spinels are also present and transformed to magnesium calcium aluminate inclusions during calcium treatment. In the current study, a fundamental kinetic model is proposed for modification of spinel inclusions considering mass transfer of solute through the boundary layer to the inclusion and diffusion within the product layer. Sensitivity analysis has been carried out to determine rate-controlling step for spinel modification. The rate of transformation of spinel inclusions is compared with alumina inclusions reported in the current authors’ prior model.

36 MATERIALS SCIENCE↗

Development of Tailor-Made Superabsorbent Polymers: Review of Key Aspects from Raw Material to Kinetic Model

The ability of superabsorbent polymers (SAP) to absorb and retain a large amount of water and/or aqueous solution enables their applications in areas such as agriculture, medicine and water treatment. This versatility has attracted the attention of researchers interested in developing new materials, especially bio-based polymers. As a result, there are numerous studies reporting different types of monomers and crosslinks to produce SAP. Conversely, in the literature there is a lack of discussion of the modeling approaches to describe SAP production. This paper makes a comprehensive and systematic review of the main issues for the synthesis of new SAP materials: monomer and crosslink types, as well as the kinetic model. A database on the raw materials to produce SAP for the main applications should be useful for those interested in producing new SAP materials. The mathematical model correlates structural properties with the synthesis conditions, thus supporting the development of tailor-made products.

36 MATERIALS SCIENCE↗

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

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

09 BIOMASS FUELS↗

Probing the Effects of NO x and SO x Impurities on Oxy-Fuel Combustion in Supercritical CO 2 : Shock Tube Experiments and Chemical Kinetic Modeling

The direct-fired supercritical carbon dioxide cycles are one of the most promising power generation methods in terms of their efficiency and environmental friendliness. Two important challenges in implementing these cycles are the high pressure (300 bar) and high CO 2 dilution (>80%) in the combustor. The design and development of supercritical oxy-combustors for natural gas require accurate reaction kinetic models to predict the combustion outcomes. The presence of a small amount of impurities in natural gas and other feed streams to oxy-combustors makes these predictions even more complex. During oxy-combustion, trace amounts of nitrogen present in the oxidizer is converted to NO x and gets into the combustion chamber along with the recirculated CO 2 . Similarly, natural gas can contain a trace amount of ammonia and sulfurous impurities that get converted to NO x and SO x and get back into the combustion chamber with recirculated CO 2 . In this work, a reaction model is developed for predicting the effect of impurities such as NO x and SO x on supercritical methane combustion. The base mechanism used in this work is GRI Mech 3.0. H 2 S combustion chemistry is obtained from Bongartz et al. while NOx chemistry is from Konnov. The reaction model is then optimized for a pressure range of 30–300 bar using high-pressure shock tube data from the literature. It is then validated with data obtained from the literature for methane combustion, H 2 S oxidation, and NO x effects on ignition delay. The effect of impurities on CH 4 combustion up to 16 atm is validated using NO x -doped methane studies obtained from the literature. In order to validate the model for high-pressure conditions, experiments are conducted at the UCF shock tube facility using natural gas identical mixtures with N 2 O as an impurity at ~100 bar. Current results show that there is a significant change in ignition delay with the presence of impurities. A comparison is made with experimental data using the developed model and predictions are found to be in good agreement. The model developed was used to study the effect of impurities on CO formation from sCO 2 combustors. It was found that NO x helps in reducing CO formation while the presence of H 2 S results in the formation of more CO. The reaction mechanism developed herein can also be used as a base mechanism to develop reduced mechanisms for use in CFD simulations.

Energy & Fuels↗

Modified Point-Kinetics Model for Neutron Precursors and Fission Product Behavior for Fluid-Fueled Molten Salt Reactors

Fluid-fueled nuclear reactors, such as molten salt reactors (MSRs), have recently gained significant interest. These advanced reactors represent a potential revolutionary shift in the implementation of nuclear power, and as a broad class of reactors, they have the potential to directly address many U.S. energy policy objectives. Fuel that is dissolved in the coolant requires methods to account for the birth, decay, and transport of fission products not only in the core but also throughout the loop and any auxiliary systems, such as off-gas, to which liquid fuel flows, gaseous products are carried, or solid particulates plate out. System models are particularly well suited to explore the wide range of phenomena that are associated with fluid-fueled systems, especially for safeguards analysis. However, before system dynamics can be explored, the compositions of fission products of the salt throughout the loop must be determined as they drive the dynamic behavior of a reactor.This paper describes the derivation of a modified point-kinetics model for obtaining a first-order approximation of the behavior of a salt-fueled system in which neutron precursors and fission products are born in the fuel-salt and transported outside the core. This paper also provides verification of the model using a steady-state analytic solution and provides additional cases exploring the response under transient cases. This model establishes a baseline model that can be used to explore the dynamic response of fluid-fueled reactors and to investigate important safeguards issues such as mass accountability of source terms. Here, the model is implemented in the Oak Ridge National Laboratory–developed, Modelica-based TRANSFORM library that was developed to investigate various aspects of advanced energy systems.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗