Fluidized-bed gasification kinetics model development using genetic algorithm for biomass, coal, municipal plastic waste, and their blends
Not Available
Bhattarai, Ashish (ORCID:0009000924319431)↗
SEARCH · Engineering Papers
Results for “kinetic model”
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.
Microwave-Assisted Plastic Upcycling: Dynamic Data Reconciliation, Parameter Estimation, and Kinetic Modeling
Microwave (MW)-assisted catalytic pyrolysis offers a promising pathway for efficient plastic upcycling. This work develops an integrated modeling framework combining dynamic data reconciliation, a temperature-dependent rate model, and a yield model to represent the time-varying production rate of components in MW-assisted LDPE pyrolysis conducted in a batch reactor. An Arrhenius-type rate model with a temperature-dependent reaction order is developed. A biexponential correlation is proposed for the yield of gaseous products that enables to capture the evolving product formation behavior during conversion. In the yield correlation, one term is used to represent the initial increase in yield, reflecting the rapid formation of intermediate or primary products at the early stages of the reaction when a larger fraction of the reactant remains available. As conversion progresses, the influence of this term gradually diminishes. The other term accounts for the subsequent decrease in the predicted yield, representing secondary reactions such as further cracking or coke formation that reduce the concentration of certain products at higher conversion. The model is found to accurately represent reconciled experimental flow rate profiles from an in-house MW-assisted catalytic batch reactor for major products, including ethylene, ethane, 1-butene, and benzene, across 250−350 °C. Ethylene remains the dominant product but decreases from about 41.95% at 250 °C to 30.14% at 350 °C, while heavier products increase significantly, with 1-butene rising to nearly 8.37% and benzene reaching 2.17% at intermediate temperatures. The model shows that the ethylene production rate can be maximized at around 270 °C. The models developed in this work can be utilized for process optimization, reactor design and scale-up of microwave-assisted plastic conversion technologies, and economic analysis.
Damahe, Harish [West Virginia Univ., Morgantown, W↗
Kinetic modeling of hot tail runaway electron generation during plasma disruptions using the JOREK code
The generation of runaway electrons (REs) during disruptions poses a significant challenge for the operation of tokamaks. The production of these high-energy electrons can cause substantial damage, particularly when the plasma current is high, making it a critical concern for ITER. For the high-temperature plasmas anticipated in ITER, the primary generation of REs may be dominated by the hot tail mechanism, which consists of the acceleration of hot electrons from the pre-disruption population which have not yet thermalized with the bulk following the rapid cooling of the plasma. To account for the significant 3D effects on RE production, a hot tail modeling framework has been developed within the non-linear 3D extended MHD code JOREK. This paper presents the structure of this framework, which is based on test electrons evolving in MHD fields. The verification of the method shows good agreement with the reference DREAM code for 0D test cases, as well as for axisymmetric simulations of 15 MA ITER H-mode disruption scenarios. Furthermore, a proof-of-principle application to a DIII-D case demonstrates the framework’s capability to capture for the first time the hot tail generation in 3D MHD simulations in realistic geometry. Preliminary results suggest that the production of REs is significantly reduced by stochastic losses.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
A standardized workflow for kinetic metabolic model curation and dissemination
Kinetic metabolic models provide invaluable insights into cellular metabolism, supporting applications in synthetic biology, metabolic engineering, and systems biology. However, reproducibility and utility of these models hinge on clear and rigorous documentation, standardized annotation, and accessible visualization. This paper presents a workflow for building, annotating, visualizing, and sharing kinetic metabolic models. Our method integrates community standards and open-source tools to ensure reproducibility, interoperability, and user accessibility. This procedure enables researchers to produce reusable and well-documented kinetic models, advancing their role as powerful tools in metabolic research.
Cook, Margaret [Univ. of Washington, Seattle, WA (↗
Modeling MTS pyrolysis and SiC deposition kinetics using principal component analysis and neural networks
Accurate chemical kinetics modeling is crucial for improving the efficiency of chemical processing and synthesis of ceramic matrix composites. Detailed kinetic models are computationally expensive due to the large number of transported chemical species, while the simplified physics-based models, such as single-step global mechanisms, are efficient but often overlook key chemical intermediates and pathways. Recent deep learning approaches promise accurate and cost-effective models. Yet, they require additional closures for the transported nonlinear latent variables, complicating integration with existing solvers. In this work, we develop a hybrid linear—nonlinear reduced model for silicon carbide deposition from methyltrichlorosilane precursor by combining principal component analysis (PCA) and autoencoder (AE) neural network (NN) approaches. PCA is used to identify a smaller set of linear transport variables, enabling direct reuse of conventional transport solvers. NNs then reconstruct the full chemical state from these reduced variables. We demonstrate the method on a chemical vapor deposition reactor—comprising a gas-phase pyrolysis plug flow reactor and a heterogeneous surface reactor—over a wide range of temperatures, pressures, and residence times. Our PCA–AE model achieves high accuracy with only five transported scalars, achieving an eightfold cost reduction compared to detailed mechanisms, in both a priori (using data from the test set only) and a posteriori (coupled with a differential equation solver). In conclusion, notable errors arise primarily near training domain boundaries and for long residence times, indicating the need for domain shift indicators and better long-horizon predictions in future reduced chemistry model development.
autoencoder neural networks↗
Development of a kinetic-thermodynamic model for lime-stabilization of Na-bentonite
This study presents the first kinetic model to predict the solid and pore solution composition of Na-bentonite clay reacting with slaked lime over a period of 720 days. The model successfully accounts for most experimental data using a single kinetic rate constant. The following sequence of reactions was predicted by the model: initial rapid dissolution of portlandite within the first 7 days, leading to a decrease in pH and dissolved calcium, and concurrent formation of calcium silicate hydrates (C-S-H: jennite), calcium aluminate hydrate (C-A-H: C₄AH₁₃), calcium aluminosilicate hydrates (stratlingite) and hydrotalcite. After 7 days, jennite and stratlingite are predicted to transform into tobermorite-II, contributing to strength development up to 28 days. From 28 to 90 days, continued montmorillonite dissolution is predicted, along with minor formation of ettringite, partial tobermorite-II dissolution, and precipitation of secondary phases such as albite and talc. Experimentally, portlandite dissolution was confirmed by TGA and XRD and found to be complete within 7 days, in agreement with model predictions. However, other predicted solid-phase transformations (e.g., tobermorite-II formation and dissolution, ettringite, albite, and talc formation) could not be conclusively verified through experimental techniques. Aqueous phase measurements confirmed that the pH and Ca trends in solution, and that equilibrium was reached by 90 days.
Chemical kinetics↗
Towards UV-Models of Kinetic Mixing and Portal Matter VII: A Light Dark Photon in the 3c3L1A1B Model
The kinetic mixing (KM) portal, by which the Standard Model (SM) photon mixes with a light dark photon arising from a new U(1)DU(1)D gauge group, allows for the possibility of viable scenarios of sub-GeV thermal dark matter (DM) with appropriately suppressed couplings to the SM. This KM can only occur if particles having both SM and dark quantum numbers, here termed portal matter (PM), also exist. The presence of such types of states and the strong suggestion of a need to embed U(1)DU(1)D into a non-abelian gauge structure not too far above the TeV scale based on the RGE running of the U(1)DU(1)D gauge coupling is potentially indicative of an enlarged group linking together the visible and dark sectors. The gauge group G=SU(3)c×SU(3)L×U(1)A×U(1)B=3c3L1A1BG=SU(3)c×SU(3)L×U(1)A×U(1)B=3c3L1A1B is perhaps the simplest setup wherein the SM and dark interactions are partially unified in a non-abelian fashion that is not a simple product group of the form G=GSM×GDG=GSM×GD encountered frequently in earlier work. The present paper describes the implications and phenomenology of this type of setup.
Rizzo, Thomas G. [SLAC National Accelerator Labora↗
Toward UV models of kinetic mixing and portal matter. VII. A light dark photon in the 3 𝑐 3 𝐿 1 𝐴 1 𝐵 model
The kinetic mixing (KM) portal, by which the Standard Model (SM) photon mixes with a light dark photon arising from a new 𝑈(1) 𝐷 gauge group, allows for the possibility of viable scenarios of sub-GeV thermal dark matter with appropriately suppressed couplings to the SM. This KM can only occur if particles having both SM and dark quantum numbers, here termed portal matter, also exist. The presence of such types of states and the strong suggestion of a need to embed 𝑈(1) 𝐷 into a non-Abelian gauge structure not too far above the TeV scale based on the renormalization-group equation running of the 𝑈(1)𝐷 gauge coupling is potentially indicative of an enlarged group linking together the visible and dark sectors. The gauge group 𝐺 = 𝑆𝑈(3) 𝑐 × 𝑆𝑈(3) 𝐿 × 𝑈(1) 𝐴 × 𝑈(1) 𝐵 = 3 𝑐 3 𝐿 1 𝐴 1 𝐵 is perhaps the simplest setup wherein the SM and dark interactions are partially unified in a non-Abelian fashion that is not a simple product group of the form 𝐺 = 𝐺 SM × 𝐺 𝐷 encountered frequently in earlier work. The present paper describes the implications and phenomenology of this type of setup.
72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
A Kinetic Model-Driven Techno-Economic Analysis of Plastic Pyrolysis: Linking Process Dynamics to Economic Viability
This study employs a kinetic model integrated into Aspen Plus to predict pyrolysis product distribution under various conditions. A techno-economic assessment calculated the minimum selling price (MSP) of pyrolysis oil under different operating conditions for the baseline capacity of 100 kta, and across eight processing capacities ranging from 30 to 150 kta. The lowest MSP under the baseline capacity is estimated at $\$$420/ton, which is 33% lower than the 2023 average US crude oil price ($\$$74.6/bbl, equivalent to $\$$634/ton based on the density of pyrolysis oil). Under Monte Carlo simulation, accounting for variability in key economic and technical parameters, the mean MSP is estimated at $\$$1137/ton. The economic viability depends on feedstock price remaining below $\$$320/ton, defining the break-even feedstock price threshold. Sensitivity analysis further identifies capital investment and transportation cost as key economic drivers. Capacities beyond 90 kta show limited economies of scale benefits. Reducing product storage time cuts capital costs by 7% but raises operational risk. Uncertainty analysis suggests the economic feasibility of pyrolysis oil is unlikely to compete with crude oil without policy incentives.
petrochemicals↗
Real-fluid behavior in rapid compression machines: Does it matter?
Rapid compression machines (RCMs) have been extensively used to quantify fuel autoignition chemistry and validate chemical kinetic models at high-pressure conditions. Historically, the analyses of experimental and modeling RCM autoignition data have been conducted based on the adiabatic core hypothesis with ideal gas assumption, where real-fluid behavior has been completely overlooked, though this might be significant at common RCM test conditions. Here, this work presents a first-of-its-kind study that addresses two significant but overlooked questions for autoignition studies within RCMs in the fundamental combustion community: (i) experiment-wise, can unaccounted-for real-fluid behavior in RCMs affect the interpretation and analysis of RCM experimental data? and (ii) simulation-wise, can unaccounted-for real-fluid behavior in RCMs affect RCM autoignition modeling and the validation of chemical kinetic models? To this end, theories for real-fluid isentropic change are newly proposed and derived based on high-order Virial EoS, and are further incorporated into an effective-volume real-fluid autoignition modeling framework newly developed for RCMs. With detailed analyses, the strong real-fluid behavior in representative RCM tests is confirmed, which can greatly influence the interpretation of RCM autoignition experiments, particularly the determination of end-of-compression temperature and evolution of the adiabatic core in the reaction chamber. Furthermore, real-fluid RCM modeling results reveal that considerable error can be introduced into simulating RCM autoignition experiments when following the community-wide accepted effective-volume approach by assuming ideal-gas behavior, which can be as high as 64% in the simulated ignition delay time at compressed pressure of 125 bar and lead to contradictory validation results of chemical kinetic models. Therefore, we recommend the community to adopt frameworks with real-fluid behavior fully accounted for (e.g., the one developed in this study) to analyze and simulate past and future RCM experiments, so as to avoid misinterpretation of RCM autoignition experiments and eliminate the potential errors that can be introduced into the simulation results with the existing RCM modeling frameworks.
High-order Virial equation of state↗
Kinetic and modeling studies of the mechanism of the dehydrogenation of Mg(BH 4 ) 2 to Mg(B 3 H 8 ) 2
Since its discovery over 15 years ago, the reversible dehydrogenation of Mg(BH 4 ) 2 to Mg(B 3 H 8 ) 2 has remained one of the more intriguing hydrogen-cycling systems. While the mechanism of this reaction has been the subject of a good deal of speculation and computational studies, prior to this work it had not been probed through kinetic studies. Previous reports of the dehydrogenation of Mg(BH 4 ) 2 to Mg(B 3 H 8 ) 2 have not included kinetic studies. The present studies have shown that the dehydrogenation of Mg(BH 4 ) 2 to Mg(B 3 H 8 ) 2 is suppressed by hydrogen pressure indicating that the rate-limiting step in this process involves hydrogen elimination. Computational modeling of kinetic data obtained from monitoring the hydrogen elimination from Mg(BH 4 ) 2 to Mg(B 3 H 8 ) 2 under static vacuum over a range of temperatures supports that the dehydrogenation occurs through a reversible three-step process in which the elimination of hydrogen from the [B 3 H 10 ] − intermediate is rate limiting. A mechanism involving the low energy transfer of neighboring BH3 groups is proposed to account for the formation of [B 3 H 8 ] − at relatively low temperatures.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Kinetic analyses for solid-state phase transition of metastable amorphous-AlO x (2.5 < x ≤ 3.0) nanostructures into crystalline alumina polymorphs
Solid-solid phase change materials (SS-PCMs) hold promise for energy storage/dissipation in batteries and energetic materials. Yet, phase change kinetics for SS-PCMs undergoing metastable to semi-stable/stable phase transformations remain relatively ill-studied because trapping metastable phases remain challenging. Recently, we demonstrated the kinetic entrapment and stabilization of a highly disordered and amorphous Al-oxide phase m-AlO x @C (x~2.5-3.0) via laser ablation synthesis in solution (LASiS). We report here, to our knowledge, the first chemical kinetics analysis for S-S phase transition of the m-AlO 3 @C nanocomposites (< 5–8 nm sizes) into semi-stable equilibrium alumina phases (θ/γ-Al 2 O 3 ) via disproportionation reaction, while releasing excess trapped gases. Our results indicate the atomic density of the AlO 3 structures to be ~5–10 times less than that of the final Al 2 O 3 phases, which led to the hypothesis of a volume shrinkage process during their phase transition. Temperature-dependent X-ray diffraction studies reveal the high-temperature phase transition for m-AlO 3 → θ/γ-Al 2 O 3 to follow contracting volume kinetics model, thereby validating our earlier hypothesis. Using the geometric volume contraction model, reaction kinetics analyses from Arrhenius plots reveal the activation energy barrier for the phase transition to be ~270±11 kJ/mol. This makes the activation energy barrier nearly identical to the oxidation of micron-sized Al particles.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Numerical Simulation of Biogenic Fluid Catalytic Cracking (BFCC) Regenerators at Different Scales with MFIX-Exa
Catalytic Fast Pyrolysis (CFP) is a process that converts biomass into liquid intermediates suitable for transportation fuels by rapidly heating it in the presence of a catalyst, aiming to produce stable oils with reduced oxygen content. During CFP, the catalyst can become deactivated by the accumulation of coke, a carbon-rich deposit formed from the decomposition of biomass components. Unlike in petroleum refining, regenerating coked catalysts from biomass pyrolysis requires specific approaches due to the different chemical nature of the coke formed. An experimental technique, Temperature Programmed Oxidation (TPO), was used to study the de-coking process by gradually increasing temperature while monitoring the production of CO and CO2, which provides data for kinetic modeling. Utilizing data from TPO experiments, coke combustion kinetic model was developed to describe the rate of coke removal at different temperatures, allowing for simulation of regeneration processes. Then kinetic model is integrated into MFIX-Exa for the simulation of Biogenic Fluid Catalytic Cracker (BFCC) regenerator at different scales, enabling analysis of catalyst flow, temperature distribution, and regeneration efficiency under various operating conditions.
biogenic fluid catalytic cracking↗
A multiscale packed-bed reactor model for sustainable ethylene production via chemical looping oxidative coupling of methane
The rising global warming concerns and shale gas discovery have prompted research in the direction of greenhouse gas (GHG), such as methane, reduction and conversion. Oxidative coupling of methane (OCM) offers a pathway to low carbon-intense valorization of methane while producing ethylene, a chemical regarded as central to the petrochemical industry. Even after decades of OCM discovery, researchers keep understanding the process and underlying chemical reactions in a pursuit to achieve industrial viability for OCM. Here, in general, OCM suffers from low C 2 selectivity, yield and reactor temperature runaways due to highly exothermic nature of its reactions. Computational Fluid Dynamics (CFD) tools help analyze spatial gradients within the reactor to deeply understand the diffusion of species, mass and heat transfer phenomena. Furthermore, challenges associated with scaling up such as hot spot formation and parametric sensitivity can be addressed without having to expend on costly experiments. The current paper presents a multiscale packed-bed reactor CFD model coupled with a chemical kinetic model for the chemical looping OCM. The CFD model includes two scales i.e., macroscale for catalyst bed and microscale for individual pellets. Moreover, a chemical kinetic model based on 10 gas-phase reactions is integrated with the CFD model. An additional surface reaction for the formation of gas-phase oxygen from catalyst surface is added to account for the absence of feed oxygen. The model is calibrated against experimental results. The calibrated model captures trends in CH 4 conversion, C 2 selectivity and C 2 yield within a ± 4.35 % range across a temperature range of 700-900 °C. Moreover, model fidelity is evaluated by varying key computational parameters such as mesh resolution and time step size. The model is also verified by varying the inlet methane concentration and the gas hourly space velocity (GHSV) and comparing the results with literature. A sensitivity analysis and scale-up of the current model is undergoing.
Chemical looping↗
Closure models for the feedback of energetic particles on plasma turbulence
Energetic particles interact with the plasma surrounding them, resonating with certain types of plasma waves to stabilize them while destabilizing others, and changing the character of the background turbulence in ways that have not been fully quantified or understood. Interaction with the turbulent background plasma is key to the acceleration of many types of energetic particles including high-energy cosmic rays, solar energetic particles, and pick-up ions. The acceleration of particles is a process that would ideally be described by a kinetic model, a type of model that follows a probability distribution function (PDF) for all particles in 7-dimensional (x, y, z, v x , v y , v z , t) space. Because of the high dimensionality of a kinetic model, simulations that solve kinetic equations use the largest computational resources currently available, and are yet unable to simulate a realistic number of particles, reach the large scales necessary for astrophysical problems, and use high-precision numerical methods. Two available alternatives to kinetic plasma models have been explored for this problem, with limited success. One is a multi-fluid model produced by a cumulant discarding closure, which evolves coupled equations for the velocity, magnetic field, and internal energy for both the background plasma and the fluid of energetic particles. However, simulations that solve multi-fluid magnetohydrodynamic (MHD) equations are able to include the interaction with energetic particles only in crude ways, typically as an add-on pressure term. The second alternative is to use a hybrid method to couple a fluid description of the background plasma to a kinetic model or a Fokker–Planck model for the energetic particles. These methods are hampered by the physical modeling of the coupling. In this work, we develop a new model, which follows the PDF for all particles; this can be viewed as a step toward physical realism above a multi-fluid MHD model, while also being more computationally efficient than a kinetic model. The equations we develop model both the background plasma and the energetic particles self-consistently. Over the last decade, similar PDF methods have been developed to a high level of sophistication to model reactive flows and turbulent combustion for engineering applications. For treatment of the feedback of the energetic particles on a background plasma, a PDF closure approach should evaluate the mean characteristics, including the density, with better statistical quality than will particle-sampling procedures.
79 ASTRONOMY AND ASTROPHYSICS↗
Constraining the Low-Temperature Oxidation Mechanism of n -Hexanol through the Detection and Identification of C 6 Elusive Intermediates
Alcohol-based fuels are currently considered to be viable energy carriers for the transportation sector. Consequently, a comprehensive mechanistic understanding of the low-temperature oxidation of alcohols is essential for application in advanced low-temperature compression engines. Here, in this work, a multidimensional approach involving experimental investigations, kinetic modeling, and theoretical calculations was used to provide new insights into the low-temperature oxidation mechanism of a C 6 alcohol, n -hexanol (CH 3 (CH 2 ) 5 OH), through the detection and identification of elusive C 6 intermediates. The oxidation of n-hexanol was investigated in a jet-stirred reactor under stoichiometric conditions (ϕ = 1.0), an initial fuel concentration of 2%, a residence time of 2 s, a temperature range between 500 and 660 K, and a pressure of 700 Torr. The reactants, intermediates, and final products were detected and identified by means of molecular-beam mass spectrometry coupled with single-photon ionization employing tunable synchrotron-generated vacuum ultraviolet radiation. Chemical kinetic simulations were performed using a previously published kinetic model (Togbé et al., Energy Fuels 2010, 11, 5859−5875) to predict the reactivity of n-hexanol and elucidate the predominant formation pathways of the observed low-temperature species. Experimental photoionization efficiency curves in conjunction with ab initio calculations, enabled the identification of important low-temperature species, such as C 6 unsaturated alcohols, C 6 olefinic hydroperoxides, C 6 cyclic ethers, C 6 diones, and C 6 ketohydroperoxides. The results of this study provide valuable insight into the mechanism of the low-temperature oxidation chemistry of n -hexanol, contributing to the development of kinetic models for the low-temperature oxidation of n-hexanol and other long-chain linear alcohols.
alcohols↗
DME-Propane Blends Ignition Experiments and Modeling for Heavy-Duty Mixing Controlled Compression Engines
A blend of dimethyl ether (DME) and propane (C3H8) is being studied in a shock tube at heavy-duty engine conditions at 110 bar. Due to its intrinsic combustion properties, DME/propane blend can potentially replace diesel in mixing controlled compression ignition engines. A blend of DME/propane can reduce emissions in mixing controlled compression ignition in heavy-duty engines through modifications, which require simulations using a high-fidelity chemical kinetics model that can accurately predict the chemistry of the blend. An essential aspect of testing the chemical kinetics model is doing baseline fundamental chemistry studies on neat DME and propane, which include ignition delay time measurements. In this work, using a high-pressure shock tube, ignition delay times were gathered for DME/Propane blends at 110 bar diluted with AR to test chemical kinetic models published in the literature. These models include Aramco 3.0, NUIG V1.1, C3mech V3.3, and Dames et al. Comparisons with the experimental IDTs and models were conducted, and general agreement was observed. A sensitivity analysis was conducted, and important reactions were outlined.
Mohammed, Zuhayr Pasha [University of Central Flor↗