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 199 records · Page 11
All sourcesNASA NTRSDOE OSTI
DOE OSTI2016

WE-H-BRA-09: Application of a Modified Microdosimetric-Kinetic Model to Analyze Relative Biological Effectiveness of Ions Relevant to Light Ion Therapy Using the Particle Heavy Ion Transport System

Purpose: To evaluate the dose and biological effectiveness of various ions that could potentially be used for actively scanned particle therapy. Methods: The PHITS Monte Carlo code paired with a microscopic analytical function was used to determine probability distribution functions of the lineal energy in 0.3µm diameter spheres throughout a water phantom. Twenty million primary particles for 1H beams and ten million particles for 4He, 7Li, 10B, 12C, 14N, 16O, and 20Ne were simulated for 0.6cm diameter pencil beams. Beam energies corresponding to Bragg peak depths of 50, 100, 150, 200, 250, and 300mm were used and evaluated transversely every millimeter and radially in annuli with outer radius of 1.0, 2.0, 3.0, 3.2, 3.4, 3.6, 4.0, 5.0, 10.0, 15.0, 20.0 and 25.0mm. The acquired probability distributions were reduced to dose-mean lineal energies and applied to the modified microdosimetric kinetic model for five different cell types to calculate relative biological effectiveness (RBE) compared to 60Co beams at the 10% survival threshold. The product of the calculated RBEs and the simulated physical dose was taken to create biological dose and comparisons were then made between the various ions. Results: Transversely, the 10B beam was seen to minimize relative biological dose in both the constant and accelerated dose change regions, proximal to the Bragg Peak, for all beams traveling greater than 50mm. For the 50mm beam, 7Li was seen to provide the most optimal biological dose profile. Radially small fluctuations (<4.2%) were seen in RBE while physical dose was greater than 1% for all beams. Conclusion: Even with the growing usage of 12C, it may not be the most optimal ion in all clinical situations. Boron was calculated to have slightly enhanced RBE characteristics, leading to lower relative biological doses.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2021

Severity factor kinetic model as a strategic parameter of hydrothermal processing (steam explosion and liquid hot water) for biomass fractionation under biorefinery concept

Hydrothermal processes are an attractive clean technology and cost-effective engineering platform for bio-refineries based in the conversion of biomass to biofuels and high-value bioproducts under the basis of sustainability and circular bioeconomy. The deep and detailed knowledge of the structural changes by the severity of biomasses hydrothermal fractionation is scientifically and technological needed in order to improve processes effectiveness, reactors designs, and industrial application of the multi-scale target compounds obtained by stream explosion and liquid hot water systems. The concept of the severity factor [log 10 (R 0 )] established > 30 years ago, continue to be a useful index that can provide a simple description of the relationship between the operational conditions for biomass fractionation in second generation of biorefineries. Furthermore, this review develops a deep explanation of the hydrothermal severity factor based in lignocellulosic biomass fractionation with emphasis in research advances, pretreatment operations and the applications of severity factor kinetic model.

09 BIOMASS FUELS↗
DOE OSTI2017

A kinetic modeling study on octane rating and fuel sensitivity under HCCI conditions

With the development of advanced compression ignition engines such as the homogeneous charge compression ignition (HCCI) and the reactivity controlled compression ignition (RCCI), the long existing research and motor octane rating for fuels, i.e., RON and MON, might be insufficient to describe the fuel reactivity and performance with these new combustion strategies due to their distinctive engine operation conditions. It is also noted that with the stringent regulation for CO2 emission, a detailed knowledge of fuel reactivity is urgently needed for the selection of novel alternative bio-fuels for transportation, and optimization of their engine performance in the future. It is therefore very important to understand the reactivity and sensitivity for general classes of fuels as a function of the operating conditions in these advanced combustion engines. In this modeling study, using validated kinetic mechanism for ignition, we have simulated the performance of different classes of fuels including hydrogen, alkanes (straight and branched), alcohols and potential bio-fuels under typical HCCI engine conditions using detailed chemical kinetics. The crank angle corresponding to 50% total heat release (CA50) is utilized as an indicator for the fuel reactivity. Iso-contours of CA50 are then plotted in the engine operating parameter coordinate, such as the intake temperature, and pressure. It is shown that CA50 value can vary significantly with different engine operating conditions, and that fuels with the same RON/MON even exhibit qualitatively different CA50 iso-contours in specific operating regimes, implying differences in the intrinsic reactivity of different classes of fuels, as well as the insufficiency of RON/MON as comprehensive metrics of fuel reactivity in advanced compression ignition (ACI) strategies. In addition, the operating regime with high and low sensitivity, as well as potential non-monotonic change in fuel reactivity, has been identified by the CA50 iso-contour map for all studied fuels.

Wu, Tong↗
DOE OSTI2021

Kinetic Monte Carlo Tool for Kinetic Modeling of Linear Step-Growth Polymerization: Insight into Recycling of Polyurethanes

A kinetic Monte Carlo model of polyurethane polymerization which explicitly tracks the polymer sequences is developed and shared. Here this model is benchmarked against theoretical and experimental polyurethane data and used to investigate the effect on oligomer distributions of unequal reactivity of the first and second isocyanate to react. The reverse reactions using thermodynamic consistency are then added to the framework, and analogous to the addition polymerization concept of ceiling temperature, equilibrium chain length distributions at various temperatures are calculated. For a mixture of three monomers AA, BB, and CC, where BB and CC do not react with one another, are present in stoichiometric proportions, and have different enthalpies of reaction with AA, an odd-even effect emerges. Odd length chains are more likely than even length chains for temperatures at which BB and CC have significantly different equilibrium conversions. The concept of ceiling temperature that is typically cited for addition polymers is extended here to provide a measure of conditions under which depolymerization for recycling is favored.

condensation polymerization↗
DOE OSTI2017

Reaction Kinetic Model of Dilute Acid-Catalyzed Hemicellulose Hydrolysis of Corn Stover under High-Solid Conditions

High solid conditions are desirable in pretreatment of lignocellulosic biomass. An advanced dilute-acid pretreatment reactor has been developed at National Renewable Energy Laboratory (NREL). It is a continuous auger-driven reactor that can be operated with high-solid charge at high temperature and with short residence time resulting high productivity and high sugar concentration. Here, we investigated the kinetics of the reactions associated with dilute-acid pretreatment of corn stover, covering the reaction conditions of the NREL reactor operation: 155-185 C, 1-2 wt% sulfuric acid concentration, and 1:2 solid to liquid ratio. The experimental data were fitted to a first-order biphasic model which assumes that xylan is comprised of two different fragments: fast and slow reacting fractions. Due to the high solid loading condition, significant amount of xylose oligomers was observed during the pretreatment. We also included the oligomers as an intermediate entity in the kinetic model. The effect of acid concentration was incorporated into the pre-exponential factor of Arrhenius equation. The kinetic model with bestfit kinetic parameters has shown good agreement with experimental data. The kinetic parameter values of the proposed model were noticeably different from those previously reported. The activation energies of xylan hydrolysis are lower and the acid exponents are higher than the average of literature values. The proposed model can serve as a useful tool for design and operation of pretreatment system pertaining to corn stover.

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

Kinetic model describing self-limiting CO 2 diffusion in supported amine adsorbents

A reaction–diffusion shrinking core model describing the decay in diffusivity of supported amine sorbents upon CO 2 sorption under both simulated direct air capture and point source capture conditions is described. The decay in CO 2 diffusivity is associated with crosslinking in the aminopolymer samples and general pore blockage in the amino-silane derived samples, which occurs as CO 2 is adsorbed. The model is used to extract four kinetic parameters that govern the CO 2 uptake kinetics and working capacity: an apparent reaction rate constant, an initial effective diffusivity, and two dimensionless decay parameters. Ideally, an initially reaction limited system would allow for direct determination of the intrinsic reaction rate constant; however, sorption experiments suggest mass transfer resistances related to gas mixing, external boundary layers and intraparticle diffusion are present. Reaction rate constants are determined and agree well with theoretical values predicted with the Eyring equation parameterized using density functional theory energies from literature sources. The kinetic performance is expressed as the average effective diffusivity as a function of average conversion, which can be correlated to the dispersion of sorption sites on the support and the morphology of the active sorbent phase. Four supports are impregnated or grafted with amines, SBA-15, single-walled zeolite nanotubes (ZNT), Syloid SiO 2 , and γ-Al 2 O 3 . Due to its pore structure, γ-Al 2 O 3 supported amines can be modeled at the μm scale or at the nm scale, where the shell balance is on the μm-sized macroporous particle aggregate or on the nm-sized amine film on the surface of the Al 2 O 3 nanoparticles, which comprise the spherical particle aggregates. Faster diffusion rates are maintained under 400 ppm rather than 10% CO 2 due to a slower reaction rate giving a slower decay in diffusivity. In conclusion, this work provides a first principles kinetic analysis of CO 2 sorption where previous models are semi-empirical and use arbitrary kinetic parameters.

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

Kinetic Modelling and Experimental Studies for the Effects of Fe 2+ Ions on Xylan Hydrolysis with Dilute-Acid Pretreatment and Subsequent Enzymatic Hydrolysis

High-temperature (150-170 degrees C) pretreatment of lignocellulosic biomass with mineral acids is well established for xylan breakdown. Fe 2+ is known to be a cocatalyst of this process although kinetics of its action remains unknown. The present work addresses the effect of ferrous ion concentration on sugar yield and degradation product formation from corn stover for the entire two-step treatment, including the subsequent enzymatic cellulose hydrolysis. The feedstock was impregnated with 0.5% acid and 0.75 mM iron cocatalyst, which was found to be optimal in preliminary experiments. The detailed kinetic data of acid pretreatment, with and without iron, was satisfactorily modelled with a four-step linear sequence of first-order irreversible reactions accounting for the formation of xylooligomers, xylose and furfural as intermediates to provide the values of Arrhenius activation energy. Based on this kinetic modelling, Fe 2+ turned out to accelerate all four reactions, with a significant alteration of the last two steps, that is, xylose degradation. Consistent with this model, the greatest xylan conversion occurred at the highest severity tested under 170 ⁰C/30 min with 0.75 mM Fe 2+ , with a total of 8% xylan remaining in the pretreated solids, whereas the operational conditions leading to the highest xylose monomer yield, 63%, were milder, 150 degrees C with 0.75 mM Fe 2+ for 20 min. Furthermore, the subsequent enzymatic hydrolysis with the prior addition of 0.75 mM of iron(II) increased the glucose production to 56.3% from 46.3% in the control (iron-free acid). The detailed analysis indicated that conducting the process at lower temperatures yet long residence times benefits the yield of sugars. The above kinetic modelling results of Fe 2+ accelerating all four reactions are in line with our previous mechanistic research showing that the pretreatment likely targets multiple chemistries in plant cell wall polymer networks, including those represented by the C-O-C and C-H bonds in cellulose, resulting in enhanced sugar solubilization and digestibility.

09 BIOMASS FUELS↗
DOE OSTI2018

(Invited) Kinetic Models for the Degradation Mechanisms of PGM-Free ORR Catalysts

Platinum group metal-free (PGM-free) catalysts, based on earth-abundant elements such as carbon, nitrogen, and transition metals, have shown promising oxygen reduction reaction (ORR) activity in fuel cells. Intensive research activities have focused on the improvement of activity and understanding of active sites in PGM-free catalysts. However, the durability of PGM-free catalysts in fuel cell cathode under realistic operating conditions has been addressed to a much lesser degree than the initial ORR activity, with only qualitative degradation mechanisms proposed to account for the generally fast activity loss in the fuel cell. Herein, we introduce to PGM-free ORR electrocatalysis a quantitative logistic decay model, based on a two-step autocatalytic degradation mechanism to portray catalyst degradation behavior. This model describes well the current-time data measured in the kinetic range of the fuel cell (fuel cell cathode) operation. Here, the results further suggest that the formation of deactivation species during the ORR is responsible for the loss of catalytic active sites and ensuing cell performance loss.

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

Kinetic Model for Moisture-Controlled CO 2 Sorption

The understanding of the sorption/desorption kinetics is essential for practical applications of moisture-controlled CO 2 sorption. We introduce an analytic model of the kinetics of moisture-controlled CO 2 sorption and its interpretation in two limiting cases. In one case, chemical reaction kinetics on pore surfaces dominates, in the other case, diffusive transport through the sorbent defines the kinetics. Here, we show that reaction kinetics, which is dominant in the first case, can be expressed as a linear combination of 1st and 2nd order kinetics in agreement with the static isotherm equation derived and validated in a previous paper. The interior transport kinetics can be described by non-linear diffusion equations. By combining all carbon species into a single equation, we can eliminate — in certain limits — the source terms associated with chemical reactions. In this case, the governing equation is ∂θ/∂t = –∇ · (–D eff ∇ θ ). For a sorbent in a form of a flat sheet or a membrane, one can maintain the same functional form of a diffusion equation by introducing a generalized effective diffusivity D M that combines contributions from both surface chemical reaction kinetics and interior diffusive transport kinetics. Experimental data of transient CO 2 flux in a preconditioned commercial anion exchange membrane fit well to the 1st order model as long as very dry states are avoided, validating the theory. The observed DM for a preconditioned commercial anion exchange membrane ranges from 6.6× 10 -14 to 7.1× 10 -14 m 2 s -1 at 35°C. These small values compared to typical ionic diffusivities imply a very slow kinetics, which will be the largest issue that needs to be addressed for practical application. The collected transient CO 2 flux data are used to predict the magnitude of a continuous CO 2 pumping flux in an active membrane that transports CO 2 against a CO 2 concentration gradient. The pumped CO 2 flux is supported by water flux due to a water concentration gradient.

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

Thermodynamic and Kinetic Modeling of Co-utilization of Glucose and Xylose for 2,3-BDO Production by Zymomonas mobilis

Prior engineering of the ethanologen Zymomonas mobilis has enabled it to metabolize xylose and to produce 2,3-butanediol (2,3-BDO) as a dominant fermentation product. When co-fermenting with xylose, glucose is preferentially utilized, even though xylose metabolism generates ATP more efficiently during 2,3-BDO production on a BDO-mol basis. To gain a deeper understanding of Z. mobilis metabolism, we first estimated the kinetic parameters of the glucose facilitator protein of Z. mobilis by fitting a kinetic uptake model, which shows that the maximum transport capacity of glucose is seven times higher than that of xylose, and glucose is six times more affinitive to the transporter than xylose. With these estimated kinetic parameters, we further compared the thermodynamic driving force and enzyme protein cost of glucose and xylose metabolism. It is found that, although 20% more ATP can be yielded stoichiometrically during xylose utilization, glucose metabolism is thermodynamically more favorable with 6% greater cumulative Gibbs free energy change, more economical with 37% less enzyme cost required at the initial stage and sustains the advantage of the thermodynamic driving force and protein cost through the fermentation process until glucose is exhausted. Glucose-6-phosphate dehydrogenase (g6pdh), glyceraldehyde-3-phosphate dehydrogenase (gapdh) and phosphoglycerate mutase (pgm) are identified as thermodynamic bottlenecks in glucose utilization pathway, as well as two more enzymes of xylose isomerase and ribulose-5-phosphate epimerase in xylose metabolism. Acetolactate synthase is found as potential engineering target for optimized protein cost supporting unit metabolic flux. Pathway analysis was then extended to the core stoichiometric matrix of Z. mobilis metabolism. Growth was simulated by dynamic flux balance analysis and the model was validated showing good agreement with experimental data. Dynamic FBA simulations suggest that a high agitation is preferable to increase 2,3-BDO productivity while a moderate agitation will benefit the 2,3-BDO titer. Taken together, this work provides thermodynamic and kinetic insights of Z. mobilis metabolism on dual substrates, and guidance of bioengineering efforts to increase hydrocarbon fuel production.

23-butanediol↗
DOE OSTI2016

Transitions From Stochastic to Point Kinetics Models In Fissile Solutions

Fissile solutions form an integral part of the fuel cycle, being used to store and process both fresh fuel and waste. Fissile solutions are also found within Aqueous Homogeneous Reactors (AHRs) including research reactors such as CRAC and SILENE and also potential future systems which could be used for the production of medical isotopes. In the case of fissile solutions in the fuel cycle, avoiding criticality is a priority but accidents have occurred in which a significant number of fissions occurred. In the case of AHRs, a sustained neutron population and fission rate is desired. In both cases the ability to simulate transients where changes in reactivity cause the neutron population to increase are important. A popular method for the simulation of such transients is the point kinetic method which allows the capture of the major features of a transient at a fraction of the computational runtime and code complexity of more comprehensive modelling technologies. A limitation of the point kinetics method and other mainstream methods of predicting neutron behaviour is that they assume the neutron population is large enough that the mean behaviour will be observed. When the neutron population is low, however, the deviation around the mean behaviour is expected to be relatively large as there is an aleatoric uncertainty arising from the uncertain fate of individual neutrons. Such limitations are of particular importance when the reactivity of a system is rapidly increasing while the neutron population is low, which is often true in accident scenarios involving fissile solutions. This provides the motivation for the coupling of point kinetics and a probabilistic models of fissile solutions. This summary aims to explore some of the behaviours resulting from a low neutron population and rising reactivity using a coupled probabilistic and point kinetics model of a fissile solution. As the effects of a low neutron population are the primary concern of this summary, the point kinetics and thermal hydraulic models used will be relatively simple but these models should be sufficient to capture the most important behaviours of the system in the first minute or so of a transient. This work is an extension of the work found in [1]. In that work the methodology was applied to fast pulses produced by systems such as Godiva and its results were found to compare favourably to experimental and analytical results, such as the burst wait time, in [2]. This work focuses on the application of this methodology to a fissile solution but the method could be applied to a number of criticality safety or reactor operation scenarios through suitable modification of the point kinetics model and would produce a more accurate prediction of the peak power obtained when a system first becomes critical compared to the use of a standard point kinetics model. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2022

Elucidating the temperature dependence of TRIP in Q&P steels using synchrotron X-Ray diffraction, constituent phase properties, and strain-based kinetics models

Understanding the deformation-induced martensitic transformation (DIMT) is critical for interpreting the structure-property relationships that govern the performance of transformation-induced plasticity (TRIP) assisted steels. However, modern TRIP-assisted steels often exhibit DIMT kinetics that are not easily captured by existing empirical models based on bulk tensile strain. We address this challenge by combined bulk uniaxial tensile tests and in-situ high energy synchrotron X-ray diffraction, which resolved the phase volume fractions, stress-strain response, and microstructure evolution of each constituent phase. A modification of the Olson-Cohen model is implemented, which describes the martensitic transformation kinetics as a function of the estimated partitioned strain in austenite, rather than the bulk tensile strain. This DIMT kinetic model is used as a framework to clarify the root cause of an insufficiently understood toughness trough reported for TRIP-assisted steels during deformation at elevated temperatures. Here, the importance of the temperature-dependent toughness is discussed, based on the opportunity to modify deformation processes to tailor the DIMT kinetics and mechanical properties during forming and in service.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Global kinetic modeling of NH 3 -SCR with two sites of NH 3 storage on Cu-SSZ-13

A comprehensive selective catalytic reduction (SCR) model is developed to detail the global kinetic reactions in the catalyst which is involved in two type of sites for high- and low-temperature NH 3 adsorption/desorption, as well as a low temperature H 2 O storage. The model accounts for the formation of ammonium nitrate in fast SCR and the decomposition formation of ammonium nitrate in the absence of NO, and addresses N 2 O relative to NO and NO 2 reactions with NH 3 . The model has been validated against commercial Cu-SSZ-13 measurements of NH 3 inventory, NH 3 oxidation, NO oxidation, standard SCR, fast SCR and NO 2 SCR at both comprehensive steady-state and transient conditions. In addition, the application of the SCR model for emissions control over a transient drive cycle at cold start was demonstrated based on diesel oxidation catalyst (DOC) and SCR aftertreatment system for a passenger car. In conclusion, the results show that the SCR model can provide the detailed emissions estimations for steady-state and transient conditions, can assist with understanding the impact of the SCR reaction mechanism, and can assist with the design and optimization of SCR catalysts to develop innovative technologies for co-optimization of engine fuel economy and emissions control technologies.

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

Kinetic Modeling of the Silicothermic Reduction of Manganese Oxide from Slag

A model was developed to describe the mixed rate control kinetics in systems with MnO-SiO{sub 2}-CaO-Al{sub 2}O{sub 3} slags and Fe-Si metal droplets. During the reaction of Fe-Si droplets and slag, Mn{sup 2+} transport was found to be part of controlling the system but could not sufficiently describe all initial conditions (chiefly, changes to initial silicon content). The current model describes the kinetics of the stated system and offers answers to the question of rate control; the model has been fitted to nine datasets of varying initial conditions including initial [Si], initial (MnO), initial droplet size, and reaction temperature. The fitted mass transfer coefficients for metal and slag were 2.3 × 10{sup −4} and 6.7 × 10{sup −4} m/s, respectively; these values are constant across the nine datasets. Previous claims about the efficiency of silicon usage in reducing manganese have been modified; it appears that the formation of silica is favored throughout the reaction, but that the formation of a silicon monoxide gas layer on the metal surface dramatically impacts the rate of reaction. As a measure of overall fit of the model, the average of the root-mean-square errors for all datasets is 14 pct. Mass transport in slag is twice as influential on rate control as the metal phase. The simultaneous transport of both Mn{sup 2+} and silicate dimers controls mass transport in the slag. Both [Si] and [Mn] can control from the metal side but are dependent on the initial conditions.

36 MATERIALS SCIENCE↗
DOE OSTI2016

WE-H-BRA-01: BEST IN PHYSICS (THERAPY): Nano-Dosimetric Kinetic Model for Variable Relative Biological Effectiveness of Proton and Ion Beams

Purpose: Recent clonogenic cell survival and γH2AX studies suggest proton relative biological effectiveness (RBE) may be a non-linear function of linear energy transfer (LET) in the distal edge of the Bragg peak and beyond. We sought to develop a multiscale model to account for non-linear response phenomena to aid in the optimization of intensity-modulated proton therapy. Methods: The model is based on first-principle simulations of proton track structures, including secondary ions, and an analytical derivation of the dependence on particle LET of the linear-quadratic (LQ) model parameters α and β. The derived formulas are an extension of the microdosimetric kinetic (MK) model that captures dissipative track structures and non-Poissonian distribution of DNA damage at the distal edge of the Bragg peak and beyond. Monte Carlo simulations were performed to confirm the non-linear dose-response characteristics arising from the non-Poisson distribution of initial DNA damage. Results: In contrast to low LET segments of the proton depth dose, from the beam entrance to the Bragg peak, strong deviations from non-dissipative track structures and Poisson distribution in the ionization events in the Bragg peak distal edge govern the non-linear cell response and result in the transformation α=(1+c-1 L) α-x+2(c-0 L+c-2 L^2 )(1+c-1 L) β-x and β=(1+c-1 L)^2 β-x. Here L is the charged particle LET, and c-0,c-1, and c-2 are functions of microscopic parameters and can be served as fitting parameters to the cell-survival data. In the low LET limit c-1, and c-2 are negligible hence the linear model proposed and used by Wilkins-Oelfke for the proton treatment planning system can be retrieved. The present model fits well the recent clonogenic survival data measured recently in our group in MDACC. Conclusion: The present hybrid method provides higher accuracy in calculating the RBE-weighted dose in the target and normal tissues.

61 RADIATION PROTECTION AND DOSIMETRY↗
DOE OSTI2023

Advancement of hybrid fluid-kinetic modeling for HEDP and ICF science

We report on the development progress of a hybrid fluid-kinetic code for simulating fluids and plasmas in a wide range of environments, such as laser–matter interactions, inertial confinement fusion, magnetic confinement fusion, and pulsed power. The suite of numerical tools under development utilizes heterogeneous computer architectures and leverages the benefits of particle–based simulation techniques. By working to combine the kinetic particle-in-cell (PIC) model with a particle-based fluid simulation technique, such as smoothed particle hydrodynamics, we are developing a flexible framework capable of accurately modeling complex flows within and between kinetic and fluid regimes. The TriForce code is under development as a C++ framework for parallel, 3D, particle-based, hybrid fluid-kinetic plasma simulations. The fluid half of TriForce will be based upon the meshless smoothed-particle-hydrodynamics (SPH) approach, well-suited for shear, mixing, and turbulence, whereas the kinetic half resembles a traditional particle-in-cell (PIC) code; other particle-based approaches to fluid modeling that do use a mesh are also possible to use and are under investigation. Maxwell’s electromagnetic field equations are solved either via explicit or implicit algorithms or approximated via resistive magnetohydrodynamics (MHD) using an Ohm’s law and resulting induction equation (extended MHD is under development). A primary goal of enabling direct comparisons, from the same code, between results from the variants of MHD and implicit electromagnetic solutions is to improve our fundamental understanding of systems with magnetic fields. The code is under development to recover results from both radiation-MHD and fully kinetic codes in those limits, and is continuing to be developed from other follow-on grants to operate in between where both descriptions may co-exist and interact. For certain applications, it is desired for a simulation to contain fluid ions and electrons as well as kinetic ions and electrons. Typically, it is too computationally intensive to model a full-scale ICF or HEDP experiment fully kinetically since many cycles are expended with very small time steps on modeling the fluid part of a material that is well treated by the fluid approximation. In this case, many traditional PIC particles can be replaced with a single fluid particle representing the thermal part of the distribution function, and there are fewer needed kinetic particles, which describe the non-thermal part and can be sub-cycled relative to the fluid particle advance. Furthermore, a pure fluid code may, depending on the problem, simply lack many physically important details that are beyond the scope of its reduced approximations and assumptions. In this report, we summarize the objectives achieved in the development of the collisional and kinetic half of the code, and the physics problems to which the code has been applied in the areas of advanced and innovative fusion concepts, pulsed power, and magneto-inertial fusion.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Nucleation kinetics model for primary crystallization in Al–Y–Fe metallic glass

The high density of aluminum nanocrystals (>10 21 m -3 ) that develop during the primary crystallization in Al-based metallic glasses indicates a high nucleation rate (~10 18 m -3 s -1 ). Several studies have been advanced to account for the primary crystallization behavior, but none have been developed to completely describe the reaction kinetics. Recently, structural analysis by fluctuation electron microscopy has demonstrated the presence of the Al-like medium range order (MRO) regions as a spatial heterogeneity in as-spun Al 88 Y 7 Fe 5 metallic glass that is representative for the class of Al-based amorphous alloys that develop Al nanocrystals during primary crystallization. From the structural characterization, an MRO seeded nucleation configuration is established, whereby the Al nanocrystals are catalyzed by the MRO core to decrease the nucleation barrier. The MRO seeded nucleation model and the kinetic data from the delay time (τ) measurement provide a full accounting of the evolution of the Al nanocrystal density (N v ) during the primary crystallization under isothermal annealing treatments. Moreover, the calculated values of the steady state nucleation rates (J ss ) predicted by the nucleation model agree with the experimental results. Moreover, the model satisfies constraints on the structural, thermodynamic, and kinetic parameters, such as the critical nucleus size, the interface energy, and the volume-free energy driving force that are essential for a fully self-consistent nucleation kinetics analysis. In conclusion, the nucleation kinetics model can be applied more broadly to materials that are characterized by the presence of spatial heterogeneities.

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

Autoignition of trans -decalin, a diesel surrogate compound: Rapid compression machine experiments and chemical kinetic modeling

Decahydronaphthalene (decalin), with both cis and trans isomers, is a bicyclic alkane that is found in aviation fuels, diesel fuels, and alternative fuels from tar sands and oil shales. Between the two decalin isomers, trans-decalin has a lower cetane number, is energetically more stable, and has a lower boiling point. Moreover, trans-decalin has often been chosen as a surrogate component to represent two-ring naphthenes in transportation fuels. Recognizing the importance of understanding the chemical kinetics of trans-decalin in the development of surrogate models, an experimental and modeling study has been conducted. Experimentally, the autoignition characteristics of trans-decalin were investigated using a rapid compression machine (RCM) by using trans-decalin/O2/N2 mixtures at compressed pressures of PC=10–25 bar, low-to-intermediate compressed temperatures of TC=620–895 K, and varying equivalence ratios of Φ=0.5, 1.0, and 2.0. These new experimental data demonstrate the effects of pressure, fuel loading, and oxygen concentration on autoignition of trans-decalin. The current RCM data of trans-decalin at lower temperatures were also found to complement well with the literature shock tube data of decalin (mixture of cis+trans) at higher temperatures. Furthermore, a chemical kinetic model for the oxidation of trans-decalin has been developed with new reaction rates and pathways, including, for the first time, a fully-detailed representation of low-temperature chemical kinetics for transdecalin. This model shows good agreement with the overall ignition delay results of the current RCM experiments and the literature shock tube studies. Chemical kinetic analyses of the developed model were further conducted to help identify the fuel decomposition pathways and the reactions controlling the autoignition at varying conditions.

10 SYNTHETIC FUELS↗