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 1,567 records · Page 87
All sourcesNASA NTRSDOE OSTI
DOE OSTI2019

Pilot-Scale CFP Commissioning: Creative Problem Solving and Lessons Learned

We have commissioned a pilot-scale, ex-situ catalytic fast pyrolysis (CFP) reactor system in the Thermal and Catalytic Process Development Unit (TCPDU) at the National Renewable Energy Laboratory. The FCC-style CFP system, including a riser and continuous regeneration, was added to the existing half-ton-per-day biomass fast pyrolysis unit and was designed to run a zeolite-type catalyst. The system is unique due to its scale and catalyst flow control scheme. As the design of this system was informed by early data from a bench-scale reactor, several technical scale-up challenges had to be overcome. Additionally, physical constraints and the complexities of biomass vapors required and inspired creative solutions. This talk will discuss the specific challenges we found and addressed, including an oxygen deficiency in the regenerator and plugging in the exit streams. We explored several options for increasing oxygen availability and selected a cyclone-catalyst return design on the exit line so we could increase air flow through the regenerator. We found that hot gas filtration is essential in any exit streams containing water vapor and trace amounts of catalyst; any combination of catalyst and condensate quickly makes mud that plugs equipment and severely limits run time. Duplicate filters are also required to extend run time. Most of the equipment and solutions were designed and installed by our team of engineers and technicians, working closely together in an iterative process to design, install, operate, then redesign. Collaborations with other national laboratories with computational fluid dynamics expertise helped inform designs and identify creative solutions to these scale-up challenges. Ultimately, data from the reactor system will be used to validate a comprehensive kinetic model of CFP that can be used to better design and inform new systems at various scales.

09 BIOMASS FUELS↗
DOE OSTI2021

Intra-Molecular Isotope Study of Light Hydrocarbons

Light hydrocarbons (C 1 -C 10 ) are the largest fraction of petroleum hydrocarbons and some of the most ubiquitous and mobile fluids within sedimentary basins. It has been increasingly recognized that hydrocarbons of varying chain-lengths can be stable even under mantle conditions and that their fluxes may play a key role in the carbon dynamics of the Earth. Major formation mechanisms of these light hydrocarbons include microbial, thermogenic, and abiogenic origins. Despite their proven utility, conventional approaches using the bulk or global stable isotopes of light hydrocarbons are often not sufficient to distinguish their sources (particularly, biogenic vs. abiogenic), maturation stages/temperatures, and migration-degradation processes. Most natural compounds, including light hydrocarbons, are composed of a set of diverse isotopic molecules that differ in the number of isotopic substitutions and/or positions of isotopic substitution within a given molecules. Many of light hydrocarbon molecules contain hydrogen and carbon in energetically non-equivalent positions: isotopomers or position-specific (intra-molecular) isotope fractionation. Position-specific isotopic compositions can differ due either to equilibrium or to path-dependent kinetic or biological processes. Thus, the position-specific isotopic fractionation between different carbon and hydrogen positions within a single molecular compound could be used regardless of their reactants and other products. We hypothesize that the formation, transport and degradation processes of light hydrocarbons in the subsurface induce and control significant, perhaps unique position-specific 13 C/ 12 C and 2 H/ 1 H isotope fractionation. In this proposal, we propose to develop a new dimension in stable isotope geochemistry of light hydrocarbons: intra-molecular isotope study of light hydrocarbons. The specific objectives of our proposal include: (a) develop robust and routine quantitative nuclear magnetic resonance (NMR) methods for high-precision/high-accuracy isotopomer analysis of propane and other light hydrocarbons (C 3 -C 10 ), (b) gain a first clear picture on position-specific δ 2 H and δ 13 C values of select light hydrocarbons of both natural and industrial sources, including samples of natural gases (C 3 and possibly C 4 -C 5 ) collected from select sedimentary basins with conventional and unconventional reservoirs, (c) acquire information on position-specific 13 C/ 12 C and 2 H/ 1 H isotope fractionation of light hydrocarbons in key natural processes (thermogenic formation and oxidative degradation) based on well-controlled laboratory experiments, and (d) establish a general framework for interpreting/modeling position-specific 13 C/ 12 C and 2 H/ 1 H isotope compositions of light hydrocarbons. Intra-molecular isotope study of light hydrocarbons is expected to open up great potentials that the bulk and position-specific isotope fractionation of light hydrocarbons, in combination with mathematical/theoretical kinetic models, could provide unprecedented wealth of information on the sources, transport, and sinks of light hydrocarbons in sedimentary basins and deep Earth.

03 NATURAL GAS↗
DOE OSTI2019

Exploring Molecular Mechanisms of Cosolvent Enhanced Biomass Deconstruction: An Overview of Recent Progress

In this chapter, we explore the mechanisms driving efficient thermochemical deconstruction and fractionation of lignocellulosic biomass in cosolvent mixtures for the production of biofuels and bioproducts. New approaches that combine molecular dynamics simulations with experimental characterization highlight the important role of cosolvents and their influence on molecular processes relevant to overcoming the natural recalcitrance of lignocellulose. From localized phase behavior on cellulose fibers to changes in lignin conformation, insights into cosolvent effects on biomass deconstruction involve understanding complex interactions between the solvent–solvent, solvent–substrate, and solvent–catalyst. Recently coordinated research efforts at multiple scales using computational tools, nano-scale imaging, and kinetic models have significantly advanced identification of high-performance cosolvents, shedding light on the mechanisms behind solvent-biomass interactions. Specifically, multifunctional cosolvents have been employed in leading biomass “solvolysis” technologies that favorably coordinate with dilute acids to solubilize and depolymerize all three primary biomass fractions of biomass: cellulose, hemicellulose, and lignin. As a result, recent work in studying solvent effects of multifunctional cosolvents has led to a paradigm shift toward the use of first principles to design biomass deconstruction techniques with the aim of maximizing total carbon utilization in biomass. The synergistic pairing of cosolvents with biomass deconstruction will lead scientists and engineers to overcome key technical and economic barriers of biomass processing, enabling a successful transition of the current energy landscape toward a future bioeconomy based on biomass as an alternative to petroleum.

Mostofian, Barmak↗
DOE OSTI2021

VFA Biorefinery Design for Informed Production of Ground and Aviation Fuels

The rising demand for low-carbon intensity fuel options requires evaluation of a system which can support their production economically and sustainably. Volatile fatty acids (VFAs) ranging from C2 to C8 can be derived in high yield from arrested anaerobic digestion of biomass. This talk provides an overview of our work in which we upgrade wet waste-derived VFAs using ketonization to elongate their carbon backbone to reach a range relevant to jet or diesel fuels. Kinetic models were used to predict ketone profiles from varying VFA profiles, thereby informing potential carbon flow to either (a) mixed paraffins for use in diesel or aviation fuel, or (b) ethers for use in diesel fuel. To do this, the anticipated products were screened for critical fuel characteristics using predictive tools. The criteria for neat and blended bioblendstocks were chosen based on conventional petrofuel requirements, and they were applied to inform fuel targets, identify limiting characteristics, and guide conversion development. While paraffins can serve as either diesel or aviation fuels, the latter has stringent criteria which include a precise distillation curve range tied to carbon distribution. Fuels which fall outside this range typically fail to meet other property metrics, and as such flashpoint and viscosity were identified as potentially limiting properties of this VFA aviation fuel. However, paraffin fractions which fall outside aviation criteria may meet diesel fuel requirements, which are looser except for the flashpoint minimum. Flashpoint was also a concern for smaller VFA diesel ethers, which posed an additional oxygenate risk of high water solubility. This fuel-informed conversion design process was demonstrated through production of paraffins with reduced sooting as compared to petrodiesel (~34%), and increased renewable aviation fuel blend levels (>70%). VFA-derived ethers improved petrodiesel by increasing fuel autoignition quality by 56% and reducing sooting by 86%. These VFA aviation and diesel fuels could enable greenhouse gas (GHG) reductions of over 165% and 50%, respectively, relative to their purely fossil-derived counterparts. Collectively, this work provides (i) an overview of how we leveraged the flexibility of VFAs for fuel production; and (ii) insight into the potential of a VFA biorefinery to accommodate varied feed and fuel applications, while also considering the merit of tailored fuel production pathways and GHG reduction potential.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2024

Effect of Pressure on Structure and Extinction of Counterflow Premixed Ammonia/Air Flames

With the growing interest in ammonia as a carbon-free fuel, numerous recent research efforts on fundamental combustion characteristics of ammonia have been reported. Despite this attention, there remain fundamental aspects that lack experimental investigation or validated chemical kinetic models. Specifically, the understanding of ammonia flame structure and the quantification of extinction limits at elevated pressures are incomplete. Furthermore, the chemistry of NOx in relation to ammonia flames has not been thoroughly examined, especially at pressure, and it is useful to explore the operating conditions that exhibit favorable flame stability with minimal NOx formation. In the current work, key characteristics of ammonia/air combustion, including strain-induced extinction limits and NO* chemiluminescence, are experimentally investigated over a wide range of equivalence ratios and pressures. This paper presents initial results in a study directed to improved comprehension of ammonia combustion chemistry, and acquisition of new datasets under elevated pressure and temperature conditions, to better refine rate expressions for important reactions in ammonia combustion.

08 HYDROGEN↗
DOE OSTI2018

Buffer Material Degradation in High Level Radioactive Waste Disposal - 18454

The concept of geological disposal for high-level radioactive waste (HLW) is based on a multi-barrier system comprising the natural geological barrier that repository host rock provides and its surroundings, as well as an engineered barrier system (EBS) that includes waste form, waste canisters, buffer materials, and backfill. The multi-barrier system is expected to perform its desired functions to isolate waste from the biosphere and achieve HLW disposal. Because of the interaction processes of radionuclide decay heat with hydraulic and water-mineral reactions in the radioactive waste repository, the bentonite buffer degradation resulting from the external temperature and pressure warrants investigation. However, seldom addressed in Taiwan is the influence of thermo-hydro-chemical (T-H-C) processes on buffer material degradation in the EBS of HLW disposal repositories as related to smectite clay dehydration. Therefore, we adopted the chemical kinetic model of smectite dehydration to calculate the amount of water expelled from smectite clay minerals caused by higher temperatures of waste decay heat. We discovered that smectite clay could cause volume shrinkage because of a loss of interlayer water, which can lead to bentonite buffer compression. The results of this study may be used in advanced research on the evolution of bentonite degradation for both performance assessments and safety analyses of final HLW disposal in Taiwan. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2024

Recovery of valuable metals from electronic waste using a novel ammonia-based hydrometallurgical process

The growing quantity of waste electrical and electronic equipment (WEEE), also known as electronic waste (E-waste), has been an area of growing public concern. The abundance of valuable metals contained in waste printed circuit boards (WPCBs) has made it a promising secondary resource, especially for Cu and Au. Although the recovery of metals from WPCBs via hydrometallurgical routes has been studied extensively over the past 20 years, most of the research has been limited in the laboratory. In current hydrometallurgical processes, strong acids and expensive oxidizers are often used to ensure a high recovery of metals without considering the sustainability aspects of the environment and economics. To improve upon current hydrometallurgical offerings, the current study seeks to develop an energy-saving, environmentally friendly, economic and sustainable process to efficiently recover the valuable metals from real-world end-of-life WPCBs. The new contributions presented in this study are 1) design and evaluation of a comprehensive hydrometallurgical flowsheet; 2) employment of real end-of-life PCBs as feed materials in an investigation on Cu-NH3 leaching kinetics; 3) further application of kinetic model on a counter flow process simulation; and 4) evaluation of the influences by co-existing metals in Au-S2O3 leaching and recommendation for favorable leaching conditions.

01 COAL, LIGNITE, AND PEAT↗
DOE OSTI2023

AN ACID BAKING APPROACH TO ENHANCE RARE EARTH ELEMENT RECOVERY FROM BITUMINOUS COAL SOURCES

Rare earth elements (REE) are a group of 17 elements typically classified as light and heavy rare earth elements, which play a crucial role in developing the latest technologies for energy, defense, and medical sectors. Even though REEs have been found in more than 200 minerals, only bastnaesite, monazite and xenotime are commercially exploited for REE extraction. However, the recent exponential increase in REE demand has spurred countries such as the United States into research for the extraction of REEs from secondary sources such as coal, acid mine drainage, and coal ash. Several coal sources (e.g., Fire Clay seam coal) across the United States have been identified to contain elevated concentrations of rare earth elements (>600 ppm), and various researchers have investigated the feasibility of both physical and hydrometallurgical extraction techniques for rare earth concentration and subsequent extraction. However, both physical and direct leaching were concluded to be inefficient for RE beneficiation and extraction due to low recoveries. Alternatively, thermal treatment provides promising means for RE recovery from bituminous coal sources. The positive impact of thermal treatment/calcination was established to be due to the decarbonization and dehydroxylation of the clays, which released entrapped rare earth elements within the dominant minerals and converted them into an acid-soluble form. Nonetheless, the improvement in recovery was limited to the light REEs (LREEs) with an insignificant increase in the heavy REEs (HREE). It was demonstrated that the light and heavy REEs in the material were associated with difficult-to-leach minerals such as monazite, xenotime, and zircon, which were not decomposed by simple calcination due to their high thermal stability. Hence, roasting in the presence of chemicals was necessary to ensure the decomposition of those REE-containing minerals. As such, this study was focused on the acid baking treatment of bituminous coals with an aim to enhance REE recovery, especially HREEs. Based on the presence of REE minerals like monazite and xenotime, three pre-leach treatment methods, i.e., 1) roasting, 2) direct acid baking, and 3) acid baking after roasting were investigated. Roasting tests at 600 ⁰C revealed that the recovery of light REEs (LREEs) was enhanced while the recovery of HREEs remained relatively unaffected. LREE and HREE recovery values of 38.3% and 21.3%, respectively, were achieved using a 50 g/L sulfuric acid solution at 5% solid concentration and a solution temperature of 75 ⁰C for 2 hours. Comparatively, direct acid baking at 250 °C provided substantial increased LREE and HREE recovery values to approximately 49.4% and 53.0%, respectively, using an equivalent acid dosage. Recoveries were maximized to 77.0% and 79.6% for LREE and HREE, respectively, by roasting followed by acid baking. Similar results were obtained from the treatment of a second bituminous coal source. Due to strong correlations between REE and Al recovery values, tests were performed on kaolinite and illite, which were prominent clay minerals within the source coals. These experiments revealed that the REE recovery improvements were likely a result of dehydroxylation of clays and subsequent release and decomposition of REE-bearing minerals such as monazite, xenotime and zircon. Subsequently, a parametric study was conducted to identify the impact of acid baking parameters on rare earth element recovery. The factors investigated using a three-level statistical experimental program were acid baking time, acid solution concentration, baking temperature, and acid solution-to-solids ratio which were found to significantly impact REE and contaminant element (Al, Fe, and Ca) recovery. An increase in baking temperature up to around 250 ⁰C improved the light and heavy REE recovery values by more than 50 absolute percentage points relative to performances achieved when direct leaching. As aforementioned, acid baking was needed to both decompose the clay minerals and liberate the REE minerals, which allowed access for the acid to solubilize the REEs. Acid concentration of the solution used for acid baking was studied as a means of minimizing the amount of acid needed to achieve a target REE recovery. However, thermo-gravimetric and differential scanning calorimetry analysis (TGA-DSC) of sulfuric acid under oxidizing atmosphere revealed that the addition of water decreased the evaporation temperature, which explains the lower REE recovery values obtained when using acid concentrations less than 100%. Using pure sulfuric acid at an acid-to-solid ratio of 0.8:1 resulted in recovery values of around 70% for both LREEs and HREEs. The decomposition reaction time was relatively quick with 65% of the TREEs recovered within the first 10 minutes. Following the identification of optimum operating conditions through the parametric study, a systematic leaching study was carried out to examine the impact of leaching parameters, such as solid-to-liquid (S/L) ratio, temperature, and time, on REE recovery using acid baking conditions of 1:1 acid to solids ratio at 250 °C for 30 minutes. The solid-to-liquid ratio was varied from 1-20% by weight at 25 °C, 50 °C, and 75 °C solution temperatures. The results indicated that reaction time and solution temperature considerably impacted the recovery of heavy and light REEs. Interestingly, LREE recovery reduced from 68% at 5% S/L to 58% at 20%S/L, whereas HREE recovery of 78% remained unaffected. The decrease in the LREE recovery was determined to be due to La and Ce precipitation, likely through isomorphic substitution with calcium in gypsum. The kinetic data indicated that 67% LREE and 77% HREE recovery could be obtained within the first 15 minutes of the reaction, suggesting fast reaction kinetics. Furthermore, raising the solution temperature from 25 °C to 75 °C increased the LREE and HREE recovery from 60% and 32% to 67% and 79%, respectively. The kinetic modeling results demonstrated that the rate-limiting step in the LREE dissolution was diffusion and chemical reaction, whereas the HREE extraction was controlled by only chemical reaction. The leaching study concluded that using 20% S/L at 75 °C for 15 minutes maximized LREE and HREE recovery. The lab-scale precipitation study showed that Fe and Al in solution could be removed at pH 4.5 followed by REE precipitation at pH 6.0 using 6 mol/L NaOH. Finally, the bench-scale data was used to develop a process flowsheet for REE recovery from low-grade bituminous coal sources using acid baking. Finally, based on the proposed flowsheet, a concentrated RE-cake obtained through selective precipitation at pH 6.5 was re-leached using HCl at pH 1.5. The resultant leachate was used to identify the impact of various operating parameters on REE recovery and purity with an aim to maximize REE precipitation efficiency while minimizing the oxalic acid dosage. The operating parameters for this investigation were oxalic acid dosage, iron (III) contamination, solution pH and temperature. The resultant model suggested that oxalic acid dosage and reaction pH are the most significant factors for the REE precipitation efficiency, followed by the interaction of oxalic dosage and Fe concentration. Test results indicated that increasing the oxalic acid concentration from 0g/L to 80g/L improved the REE precipitation efficiency from approximately 4.2% to 95.0%. Furthermore, raising the solution pH from 0.5 to 2.5 considerably enhanced the precipitation efficiency from 0.0% to 98.9%. A solution temperature elevation decreased REE recovery, which indicated an exothermic reaction between REEs and oxalate anions. Finally, a high level of Fe contamination adversely impacted REE precipitation efficiency. The speciation analysis revealed that the dominant iron species in the solution system were Fe-(C₂O₄)₃³⁻, Fe-(C₂O₄)²⁻, and Fe-(C₂O₄)⁺, which consumed the majority of the oxalate anions

rare earth elements, acid baking, high-temperature↗
DOE OSTI2024

Influence of porous aluminosilicate grain size materials in experimental and modelling Cs + adsorption kinetics and wastewater column process

This paper focuses on the influence of the grain size of a geopolymer based adsorbent on its Cs + adsorption performances both in batch and fixed-bed process. The geopolymer phase was used as a binder to support NaY zeolite particle in a 20 wt% charged porous composite with 160 m 2 .g –1 of porous surface area. These samples were shaped with three grain sizes (50 /100/500 µm) to remove 80–90 mg/g of Cs + in batch and column operations. After their microstructural and porous characterizations, their efficiency and adsorption characteristics were investigated through adsorption isotherms and kinetic in the two processes. While the grain size has no influence on the maximal extraction capacity of the adsorbent, it strongly affects the sorption kinetic. By coupling experimental data and a modelling approach, the complex sorption mechanism was highlighted, suggesting a new insight of the contaminant sorption kinetic. Then, comparison of batch and column adsorption experiments illustrates the detailed explanation of various process parameters for column study. The results show challenges for fixed-bed column utilization by the choice of the appropriate grain size as a compromise between the material sorption kinetic and hydrodynamic considerations. Furthermore, this is of high importance to more accurately optimize the design of column adsorption to assess the transport of Cs+ in multi-porous tailored grain size materials.

36 MATERIALS SCIENCE↗
DOE OSTI2016

Development of Kinetics and Mathematical Models for High Pressure Gasification of Lignite-Switchgrass Blends

The overall objective of the current project was to investigate the high pressure gasification characteristics of a feed containing both coal and biomass. The two feed types differ in their ash contents and ash composition, particularly the alkali content. Gasification of a combined feed of coal and biomass has the potential for considerable synergies that might lead to a dramatic improvement in process economics and flexibility. The proposed study aimed to develop a detailed understanding of the chemistry, kinetics, and transport effects during high pressure gasification of coal-biomass blend feed. Specifically, we studied to develop: (a) an understanding of the catalytic effect of alkali and other inorganic species present in the biomass and coal, (b) an understanding of processing conditions under which synergistic effects of the blending of coal and biomass might be observed. This included the role of particle size, residence time, and proximity of the two feed types, (c) kinetics of high pressure gasification of individual feeds as well as the blends, and (d) development of mathematical models that incorporate kinetics and transport models to enable prediction of gasification rate at a given set of operating conditions, and (e) protocols to extend the results to other feed resources. The goal was to provide a fundamental understanding of the gasification process and guide in optimizing the configurations and design of the next generation of gasifiers. The approach undertaken was centered on two basic premises: (1) the gasification for small particles without internal mass transfer limitations can be treated as the sum of two processes in series (pyrolysis and char gasification) , and (2) the reactivity of the char generated during pyrolysis not only depends on the pressure and temperature but is also affected by the heating rates. Thus low heating rates (10-50 °C/min) typical of PTGA fail to produce char that would typically be formed at high heating rates (~10 4 °C/sec), encountered in entrained flow gasifiers. The char morphology, also a function of the heating rate, would influence the transport rates during the char gasification phase. Thus, heating rate plays a critical role through which both, pyrolysis and char gasification, are interconnected. We utilized two complementary gasification experiments: PEFR (pressurized entrained flow gasifier) and PTGA (pressurized thermo-gravimetric analyzer). The PEFR allowed us to study gasification at pressures, temperatures, and heating rates relevant for coal-biomass gasifiers. The PTGA work was useful in understanding the basic chemistry of the evolution of various gaseous species during pyrolysis. These results helped improved our understanding of the chemistry and chemical changes during pyrolysis. The role alkali metals and other inorganics in char gasification using steam and/or CO 2 was investigated. Finally, the mathematical models for char gasification without the transport effects were developed at commercial operating conditions.

09 BIOMASS FUELS↗
DOE OSTI2025

Discovering CO Adsorption and Desorption Pathways from Chemical Reaction Neural Network Modeling of Transient Kinetics Spectroscopy

Here, we demonstrate a data-driven approach to interpreting surface reactions by combining time-resolved gas pulsing infrared spectroscopy with chemical reaction neural networks (CRNNs). Using CO adsorption and desorption on Pd(111) at 460–490 K as a model system, we show how transient kinetic data can reveal detailed reaction mechanisms. Starting with a simple one-species model, we systematically evaluate increasingly complex mechanisms involving hollow and bridge site adsorption. Despite the similar goodness of fit to the same experimental absorbance data, our models predict distinct coverage dynamics for different adsorption sites. Through analysis of spectral peak stability and predicted dynamics, we identify a mechanism in which CO primarily adsorbs on bridge sites followed by rapid conversion to hollow sites as being the most physically consistent with experimental observations. This work provides a framework for extracting mechanistic insights from limited experimental data, demonstrating how machine learning can bridge the gap between transient kinetic measurements and a molecular-level understanding of surface reactions.

36 MATERIALS SCIENCE↗
NASA NTRS1984

Micromolecular modeling

A reaction kinetics based model of the photodegradation process, which measures all important rate constants, and a computerized model capable of predicting the photodegradation rate and failure modes of a 30 year period, were developed. It is shown that the computerized photodegradation model for polyethylene correctly predicts failure of ELVAX 15 and cross linked ELVAX 150 on outdoor exposure. It is indicated that cross linking ethylene vinyl acetate (EVA) does not significantly change its degradation rate. It is shown that the effect of the stabilizer package is approximately equivalent on both polymers. The computerized model indicates that peroxide decomposers and UV absorbers are the most effective stabilizers. It is found that a combination of UV absorbers and a hindered amine light stabilizer (HALS) is the most effective stabilizer system.

Guillet, J. E.↗
DOE OSTI2017

PIC simulations of wave-particle interactions with an initial electron velocity distribution from a kinetic ring current model

Whistler wave-particle interactions play an important role in the Earth inner magnetospheric dynamics and have been the subject of numerous investigations. By running a global kinetic ring current model (RAM-SCB) in a storm event occurred on Oct 23–24 2002, we obtain the ring current electron distribution at a selected location at MLT of 9 and L of 6 where the electron distribution is composed of a warm population in the form of a partial ring in the velocity space (with energy around 15 keV) in addition to a cool population with a Maxwellian-like distribution. The warm population is likely from the injected plasma sheet electrons during substorm injections that supply fresh source to the inner magnetosphere. These electron distributions are then used as input in an implicit particle-in-cell code (iPIC3D) to study whistler-wave generation and the subsequent wave-particle interactions. Here, we find that whistler waves are excited and propagate in the quasi-parallel direction along the background magnetic field. Several different wave modes are instantaneously generated with different growth rates and frequencies. The wave mode at the maximum growth rate has a frequency around 0.62ω ce , which corresponds to a parallel resonant energy of 2.5 keV. Linear theory analysis of wave growth is in excellent agreement with the simulation results. These waves grow initially due to the injected warm electrons and are later damped due to cyclotron absorption by electrons whose energy is close to the resonant energy and can effectively attenuate waves. The warm electron population overall experiences net energy loss and anisotropy drop while moving along the diffusion surfaces towards regions of lower phase space density, while the cool electron population undergoes heating when the waves grow, suggesting the cross-population interactions.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2017

Learning reduced kinetic Monte Carlo models of complex chemistry from molecular dynamics

Here, we propose a novel statistical learning framework for automatically and efficiently building reduced kinetic Monte Carlo (KMC) models of large-scale elementary reaction networks from data generated by a single or few molecular dynamics simulations (MD). Existing approaches for identifying species and reactions from molecular dynamics typically use bond length and duration criteria, where bond duration is a fixed parameter motivated by an understanding of bond vibrational frequencies. Conversely, we show that for highly reactive systems, bond duration should be a model parameter that is chosen to maximize the predictive power of the resulting statistical model. We demonstrate our method on a high temperature, high pressure system of reacting liquid methane, and show that the learned KMC model is able to extrapolate more than an order of magnitude in time for key molecules. Additionally, our KMC model of elementary reactions enables us to isolate the most important set of reactions governing the behavior of key molecules found in the MD simulation. We develop a new data-driven algorithm to reduce the chemical reaction network which can be solved either as an integer program or efficiently using L1 regularization, and compare our results with simple count-based reduction. For our liquid methane system, we discover that rare reactions do not play a significant role in the system, and find that less than 7% of the approximately 2000 reactions observed from molecular dynamics are necessary to reproduce the molecular concentration over time of methane. Furthermore, we describe a framework in this work that paves the way towards a genomic approach to studying complex chemical systems, where expensive MD simulation data can be reused to contribute to an increasingly large and accurate genome of elementary reactions and rates.

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

Kinetics and Reactor Modeling of VFA Ketonization for Sustainable Aviation Fuel Production

Ketonization of volatile fatty acids (VFAs) produced via arrested methanogenesis of wet waste represents the only unit operation of wet waste upgrading to sustainable (net-zero or negative life cycle CO2 emissions) aviation fuel not currently in industrial practice. Ketone product yields of close to 100% have been obtained during gas-phase reactions of VFAs over oxide catalysts at the laboratory scale, but design of an industrial ketonization reactor requires understanding of the impacts of reactant and product partial pressures, deactivation, and heat and mass transport phenomena on observed ketonization rates. Our work leverages rigorous kinetic analysis of ketonization processes to inform reactor scale-up efforts through packed-bed reactor modeling. We first present results of a kinetic study of ketonization of a model VFA, hexanoic acid, over an industrial ZrO2 catalyst performed in a packed-bed microreactor in conditions free of significant heat or mass transfer gradients. Major findings of the analysis include: (i) hexanoic acid saturates all catalyst active sites at relatively low partial pressure (~10 kPa) and (ii) ketonization products 6-undecanone, water, and CO2 inhibit reaction rates. Kinetic data are used to fit a rate expression quantifying the functional dependence of ketonization rate on partial pressures of VFA reactants and ketone, water, and CO2 products. A packed-bed reactor model describing vapor-phase hexanoic acid ketonization over ZrO2 extrudate catalysts is developed based on the kinetic model. The effects of (i) bed- and pellet-scale mass- and heat-transfer limitations and (ii) axial pressure drop guide our development of recommendations for optimal sizing, temperature, and influent flow composition of an industrial-scale ketonization reactor. The quantitative understanding of VFA ketonization developed in this study represents an advance toward derisking the VFA ketonization step of wet waste upgrading to sustainable aviation fuel.

biojet↗
DOE OSTI2019

Examining and Modeling Oxygen Uncoupling Kinetics of Cu-Based Oxygen Carriers for Chemical Looping with Oxygen Uncoupling (CLOU) in a Drop Tube Fluidized Bed Reactor

Chemical looping combustion with oxygen uncoupling (CLOU) is a promising carbon capture and storage (CCS) technology for conversion of gaseous and solid hydrocarbon fuels where the release of gaseous O 2 from an oxygen carrier is favored at high temperature and low O 2 partial pressure. One promising CLOU material is the copper oxide redox system (CuO–Cu 2 O). The primary objective of this study was to examine the use of a drop tube fluidized bed reactor (DT-FBR) for evaluating the kinetics of oxygen uncoupling with Cu-based oxygen carriers. Appropriate rate expressions from redox experiments are needed to model and scale up CLOU systems. Additionally, the determined oxygen uncoupling kinetic parameters were validated using a computational fluid dynamic model. Two oxygen carriers consisting of 20 wt % CuO/Al 2 O 3 and 9 wt % CuO/Al 2 O 3 were prepared by physical mixing. A third oxygen carrier sample was prepared at 9 wt % CuO/Al 2 O 3 by incipient wetness impregnation. The reaction rates of the copper oxide redox system are strongly dependent on thermodynamic effects of the oxygen partial pressure relative to the equilibrium oxygen partial pressure. The results presented in this paper offer an alternative and simplified kinetic analysis compared to that traditionally presented in the literature for the thermal reduction of copper oxide or other CLOU oxygen carriers. The reactor system used in this study allows for operating parameters to be adjusted, minimizing thermodynamic and mass transfer limitations, which eliminates the need for more complex kinetic/thermodynamic reaction models. The reaction kinetics measured in this study is compared based on the preparation method and the CuO weight percent loading. These results may aid in the development of CLOU technologies during reactor design and process modeling.

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

The Development of Kinetic and Radiation Hydrodynamics Modeling of Thermonuclear Burn Propagation in Isochoric p - 11 B Through the Support of the INFUSE Program

The report summarizes DOE INFUSE-supported work between HB11 Energy and the University of Rochester’s TriForce Institute to improve computational modeling of advanced fusion fuels, especially proton–boron-11 (p- 11 B). The project extended the TriForce particle-in-cell/Monte Carlo collision code to include physics needed for dense, high-temperature p- 11 B burn studies, including p- 11 B fusion reactions, three-alpha-particle reaction products, relativistic Coulomb collisions, large-angle nuclear scattering, bremsstrahlung radiation, inverse bremsstrahlung absorption, and photon transport. The upgraded models were verified against focused physics tests and against known deuterium–tritium burn behavior. The study then used one-dimensional spherical simulations to estimate the conditions required for thermonuclear burn propagation in isochoric p- 11 B fuel. The calculations found that burn propagation is possible in the model, but only under very extreme hot-spot conditions, such as about 7000 g/cm 3 at 500 keV or 9000 g/cm 3 at 300 keV for a 20-micron hot spot. These conditions are much more demanding than current demonstrated inertial confinement fusion hot spots. The report concludes that the INFUSE collaboration successfully advanced kinetic and radiation modeling capabilities for p- 11 B fusion and provided useful estimates of ignition requirements. However, the simulated fuel gains remain below what would be needed for practical inertial fusion energy, and further work is needed to reconcile differences among kinetic, radiation-hydrodynamic, and analytic models and to identify more achievable target designs.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Predicting synthetic mRNA stability using massively parallel kinetic measurements, biophysical modeling, and machine learning

Abstract mRNA degradation is a central process that affects all gene expression levels, though it remains challenging to predict the stability of a mRNA from its sequence, due to the many coupled interactions that control degradation rate. Here, we carried out massively parallel kinetic decay measurements on over 50,000 bacterial mRNAs, using a learn-by-design approach to develop and validate a predictive sequence-to-function model of mRNA stability. mRNAs were designed to systematically vary translation rates, secondary structures, sequence compositions, G-quadruplexes, i-motifs, and RppH activity, resulting in mRNA half-lives from about 20 seconds to 20 minutes. We combined biophysical models and machine learning to develop steady-state and kinetic decay models of mRNA stability with high accuracy and generalizability, utilizing transcription rate models to identify mRNA isoforms and translation rate models to calculate ribosome protection. Overall, the developed model quantifies the key interactions that collectively control mRNA stability in bacterial operons and predicts how changing mRNA sequence alters mRNA stability, which is important when studying and engineering bacterial genetic systems.

Cetnar, Daniel P.↗