SEARCH · Engineering Papers

Results for “Kinetic modeling”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 577 records · Page 32

All sources NASA NTRS DOE OSTI

DOE OSTI2021

Experimental and computational investigations of ethane and ethylene kinetics with copper oxide particles for Chemical Looping Combustion

In this work, reaction pathways for the oxidation of methane, ethane, and ethylene with CuO was obtained by ReaxFF Molecular Dynamics (MD) simulations between temperatures of 1000 K and 2000 K. Experiments in a fixed-bed flow reactor were preformed with methane, ethane, and ethylene at temperatures ranging from 500 K to 1000 K with time-dependent species measurements from an Electron-Ionization Molecular Beam Mass Spectrometer (MBMS), and species validation with Gas Chromatography (GC) for detection of complete and intermediate combustion products. The MBMS and GC allow for the detection of oxygenated species and larger species produced from radical reformation. The simulation and experiment agree on the production of such species as CH 3 CHO, CH 2 O, CO, and H 2 O, which allow for the creation of simple C1 and C2 reaction pathways, which can be used in kinetic models of C2 species and larger fuels such as biofuels, which inherently depend on C1 and C2 kinetics and reaction pathway. The simulation and experiment disagree on the formation of C 2 H 2 , CH 3 OH, and CO 2 with large amounts of C 2 H 2 being measured in the ethylene oxidation simulations and CH 3 OH being formed in methane oxidation simulations, while neither species were experimentally found. In the case of CO 2 large amounts of CO 2 are rapidly produced in experiments with C2 fuels at 800 K, while little-to-no CO 2 was observed in simulations. This is believed to be resulting from the extremely short timescale of the simulations, preventing total oxidation of the fuel. Here, the differences in products produced between simulation and experiment allow for the potential to modify the ReaxFF potential functions to more accurately model the experimental products of Cu–H–O–C reaction kinetics.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2020

Mechanistic modelling of SARS-CoV-2 and other infectious diseases and the effects of therapeutics

Modern viral kinetic modeling and its application to therapeutics is a field that attracted the attention of the medical, pharmaceutical and modeling communities during the early days of the AIDS epidemic. Its successes lead to applications of modeling methods not only to HIV but a plethora of other viruses, such as hepatitis C virus (HCV), hepatitis B virus (HBV) and cytomegalovirus (CMV), which along with HIV cause chronic diseases, and viruses such as influenza, respiratory syncytial virus (RSV), West Nile virus (WNV), Zika virus, and SAR-CoV-2, which generally cause acute infections. Here we first review the historical development of mathematical models to understand HIV and HCV infections and the effects of treatment by fitting the models to clinical data. Furthermore, we then focus on recent efforts and contributions of applying these models towards understanding SARS-CoV-2 infection and highlight outstanding questions where modeling can provide crucial insights and help to optimize non-pharmaceutical and pharmaceutical interventions of the COVID-19 pandemic. The review is written from our personal perspective emphasizing the power of simple target cell limited models that provided important insights and then their evolution into more complex models that captured more of the virology and immunology. To quote Albert Einstein “Everything should be made as simple as possible, but not simpler”, and this idea underlies the modeling we describe below.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2022

Configuration Space Integration for Adsorbate Partition Functions: The Effect of Anharmonicity on the Thermophysical Properties of CO–Pt(111) and CH 3 OH–Cu(111)

A method for computing anharmonic thermophysical properties for adsorbates on metal surfaces has been extended to include libration, or frustrated rotation. Classical phase space integration is used with Monte Carlo sampling of the configuration space to obtain the partition function of CO on Pt(111) and CH 3 OH on Cu(111). A minima-preserving neural network potential energy surrogate is used within the integration routines. Direct state counting using discrete variable representation is used to benchmark the results. We find that the phase space integration approach is in excellent agreement with the direct state counting results. Comparison with standard models such as the harmonic oscillator indicates that anharmonicity contributes significantly to the thermodynamic properties of CH 3 OH on Cu(111). We find that there is also a considerable difference between the harmonic oscillator and phase space integration for CO on Pt(111), although the discrepancy can largely be attributed to the presence of multiple binding sites within the unit cell. We demonstrate that a multisite harmonic oscillator model might be sufficient for CO-Pt(111). A more thorough description of the potential energy surface, which can be achieved with phase space integration, is necessary for weakly bound adsorbates such as CH 3 OH. In conclusion, the thermophysical properties were used to calculate free energies of adsorption on the respective metals, and subsequently the equilibrium constants and Langmuir isotherms in relevant temperature ranges. The results show that the choice of model to obtain partition functions greatly affects the resulting surface coverages in kinetic models.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2020

Characterizing magnetized plasmas with dynamic mode decomposition

Accurate and efficient plasma models are essential to understand and control experimental devices. Existing magnetohydrodynamic or kinetic models are nonlinear and computationally intensive and can be difficult to interpret, while often only approximating the true dynamics. In this work, data-driven techniques recently developed in the field of fluid dynamics are leveraged to develop interpretable reduced-order models of plasmas that strike a balance between accuracy and efficiency. In particular, dynamic mode decomposition (DMD) is used to extract spatio-temporal magnetic coherent structures from the experimental and simulation datasets of the helicity injected torus with steady inductive (HIT-SI) experiment. Three-dimensional magnetic surface probes from the HIT-SI experiment are analyzed, along with companion simulations with synthetic internal magnetic probes. A number of leading variants of the DMD algorithm are compared, including the sparsity-promoting and optimized DMD. Optimized DMD results in the highest overall prediction accuracy, while sparsity-promoting DMD yields physically interpretable models that avoid overfitting. These DMD algorithms uncover several coherent magnetic modes that provide new physical insights into the inner plasma structure. These modes were subsequently used to discover a previously unobserved three-dimensional structure in the simulation, rotating at the second injector harmonic. Finally, using data from probes at experimentally accessible locations, DMD identifies a resistive kink mode, a ubiquitous instability seen in magnetized plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2022

On the parameter combinations that matter and on those that do not: data-driven studies of parameter (non)identifiability

Abstract We present a data-driven approach to characterizing nonidentifiability of a model’s parameters and illustrate it through dynamic as well as steady kinetic models. By employing Diffusion Maps and their extensions, we discover the minimal combinations of parameters required to characterize the output behavior of a chemical system: a set of effective parameters for the model. Furthermore, we introduce and use a Conformal Autoencoder Neural Network technique, as well as a kernel-based Jointly Smooth Function technique, to disentangle the redundant parameter combinations that do not affect the output behavior from the ones that do. We discuss the interpretability of our data-driven effective parameters, and demonstrate the utility of the approach both for behavior prediction and parameter estimation. In the latter task, it becomes important to describe level sets in parameter space that are consistent with a particular output behavior. We validate our approach on a model of multisite phosphorylation, where a reduced set of effective parameters (nonlinear combinations of the physical ones) has previously been established analytically.

Evangelou, Nikolaos↗

DOE OSTI2024

Probing interspecies metabolic interactions within a synthetic binary microbiome using genome-scale modeling

Metabolic interactions within a microbial community play a key role in determining the structure, function, and composition of the community. However, due to the complexity and intractability of natural microbiomes, limited knowledge is available on interspecies interactions within a community. In this work, using a binary synthetic microbiome, a methanotroph-photoautotroph (M-P) coculture, as the model system, we examined different genome-scale metabolic modeling (GEM) approaches to gain a better understanding of the metabolic interactions within the coculture, how they contribute to the enhanced growth observed in the coculture, and how they evolve over time. Using batch growth data of the model M-P coculture, we compared three GEM approaches for microbial communities. Two of the methods are existing approaches: SteadyCom, a steady state GEM, and dynamic flux balance analysis (DFBA) Lab, a dynamic GEM. We also proposed an improved dynamic GEM approach, DynamiCom, for the M-P coculture. SteadyCom can predict the metabolic interactions within the coculture but not their dynamic evolutions; DFBA Lab can predict the dynamics of the coculture but cannot identify interspecies interactions. DynamiCom was able to identify the cross-fed metabolite within the coculture, as well as predict the evolution of the interspecies interactions over time. A new dynamic GEM approach, DynamiCom, was developed for a model M-P coculture. Constrained by the predictions from a validated kinetic model, DynamiCom consistently predicted the top metabolites being exchanged in the M-P coculture, as well as the establishment of the mutualistic N-exchange between the methanotroph and cyanobacteria. The interspecies interactions and their dynamic evolution predicted by DynamiCom are supported by ample evidence in the literature on methanotroph, cyanobacteria, and other cyanobacteria-heterotroph cocultures.

59 BASIC BIOLOGICAL SCIENCES↗

NASA NTRS1995

Dynamics and structure of turbulent premixed flames

In earlier work (Mantel & Bilger, 1994) the structure of the turbulent premixed flame was investigated using statistics based on conditional averaging with the reaction progress variable as the conditioning variable. The DNS data base of Trouve and Poinsot (1994) was used in this investigation. Attention was focused on the conditional dissipation and conditional axial velocity in the flame with a view to modeling these quantities for use in the conditional moment closure (CMC) approach to analysis of kinetics in premixed flames (Bilger, 1993). Two remarkable findings were made: there was almost no acceleration of the axial velocity in the flame front itself; and the conditional scalar dissipation remained as high, or higher, than that found in laminar premixed flames. The first finding was surprising since in laminar flames all the fluid acceleration occurs through the flame front, and this could be expected also for turbulent premixed flames at the flamelet limit. The finding gave hope of inventing a new approach to the dynamics of turbulent premixed flames through use of rapid distortion theory or an unsteady Bernoulli equation. This could lead to a new second order closure for turbulent premixed flames. The second finding was contrary to our measurements with laser diagnostics in lean hydrocarbon flames where it is found that conditional scalar dissipation drops dramatically below that for laminar flamelets when the turbulence intensity becomes high. Such behavior was not explainable with a one-step kinetic model, even at non-unity Lewis number. It could be due to depletion of H2 from the reaction zone by preferential diffusion. The capacity of the flame to generate radicals is critically dependent on the levels of H2 present (Bilger, et al., 1991). It seemed that a DNS computation with a multistep reduced mechanism would be worthwhile if a way could be found to make this feasible. Truly innovative approaches to complex problems often come only when there is the opportunity to work close at hand with the (in this case numerical) experimental data. Not only can one spot patterns and relationships in the data which could be important, but one can also get to know the limitations of the technique being used, so that when the next experiment is being designed it will address resolvable questions. A three-year grant from the Australian Research Council has enabled us to develop a small capability at the University of Sydney to work on DNS of turbulent reacting flow, and to analyze data bases generated at CTR. Collaboration between the University of Sydney and CTR is essential to this project and finding a workable modus operandum for this collaboration, given the constraints involved, has been a major objective of the past year's effort. The overall objectives of the project are: (1) to obtain a quantitative understanding of the dynamics of turbulent premixed flames at high turbulence levels with a view to developing improved second order closure models; and (2) to carry out new DNS experiments on turbulent premixed flames using a carefully chosen multistep reduced mechanism for the chemical kinetics, with a view to elucidating the laser diagnostic findings that are contrary to the findings for DNS using one-step kinetics. In this first year the objectives have been to make the existing CTR data base more accessible to coworkers at the University of Sydney, to make progress on understanding the dynamics of the flame in this existing CTR data base, and to carefully construct a suitable multistep reduced mechanism for use in a new set of DNS experiments on turbulent premixed flames.

Bilger, R. W.↗

DOE OSTI2020

Different Kinetic Reactivities of Electrons in Distinct TiO 2 Nanoparticle Trap States

Electrons added to TiO 2 and other semiconductors often occupy trap states, whose reactivity can determine the catalytic and stoichiometric chemistry of the material. We previously showed that reduced aqueous colloidal TiO 2 nanoparticles have two distinct classes of thermally equilibrated trapped electrons, termed Red/e – and Blue/e – . Presented here are parallel optical and electron paramagnetic resonance (EPR) kinetic studies of the reactivity of these electrons with solution-based oxidants. Optical stopped-flow measurements monitoring the reactions of TiO 2 /e – with substoichiometric oxidants showed a surprising pattern: an initial fast (seconds) decrease in TiO 2 /e – absorbance followed by a secondary, slow (minutes) increase in the broad TiO 2 /e – optical feature. The analysis revealed that the fast decrease is due to the preferential oxidation of the Red/e – trap states and the slow increase results from the re-equilibration of electrons from Blue/e - to Red/e - states. This kinetic model was confirmed by freeze-quench EPR measurements. Quantitative analysis of the kinetic data demonstrated that Red/e – react ~5 times faster than Blue/e – with the nitroxyl radical oxidant 4-methoxy-2,2,6,6-tetramethyl-1-piperidinyloxyl (4-MeO-TEMPO). Similar reactivity patterns were also observed in oxidations of TiO 2 /e – by O 2 , which like 4-MeO-TEMPO is a proton-coupled electron transfer (PCET) oxidant, and by the pure electron transfer (ET) oxidant potassium triiodide (KI 3 ). Furthermore, this suggests that the faster intrinsic reactivity of one trap state over another on the seconds–minutes time scale is likely a general feature of reduced TiO 2 reactivity. This differential trap-state reactivity is likely to influence the performance of TiO 2 in photochemical/electrochemical devices, and it suggests an opportunity for tuning catalysis.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

NASA NTRS1984

Theoretcial studies of solar-pumped lasers

A method of pumping a COhZ laser by a hot cavity was demonstrated. The cavity, heated by solar radiation, should increase the efficiency of solar pumped lasers used for energy conversion. Kinetic modeling is used to examine the behavior of such a COhZ laser. The kinetic equations are solved numerically vs. time and, in addition, steady state solutions are obtained analytically. The effect of gas heating filling the lower laser level is included. The output power and laser efficiency are obtained as functions of black body temperature and gas ratios (COhZ-He-Ar) and pressures. The values are compared with experimental results.

Harries, W. L.↗

NASA NTRS1989

Conservation equations and physical models for hypersonic air flows in thermal and chemical nonequilibrium

The conservation equations for simulating hypersonic flows in thermal and chemical nonequilibrium and details of the associated physical models are presented. These details include the curve fits used for defining thermodynamic properties of the 11 species air model, curve fits for collision cross sections, expressions for transport properties, the chemical kinetics models, and the vibrational and electronic energy relaxation models. The expressions are formulated in the context of either a two or three temperature model. Greater emphasis is placed on the two temperature model in which it is assumed that the translational and rotational energy models are in equilibrium at the translational temperature, T, and the vibrational, electronic, and electron translational energy modes are in equilibrium at the vibrational temperature, T sub v. The eigenvalues and eigenvectors associated with the Jacobian of the flux vector are also presented in order to accommodate the upwind based numerical solutions of the complete equation set.

Peter A. Gnoffo↗

NASA NTRS2000

Modeling of HiPco Process for Carbon Nanotube Production

High-pressure carbon monoxide (HiPco) reactor, developed at Rice University, is used to produce single-walled carbon nanotubes (SWNT) from gas-phase reactions of iron carbonyl and nickel carbonyl in carbon monoxide at high pressures (10 - 100 atm). Computational modeling is used to better understand the HiPco process. In the present model, decomposition of the precursor, metal cluster formation and growth, and carbon nanotube growth are addressed. Decomposition of precursor molecules is necessary to initiate metal cluster formation. The metal clusters serve as catalysts for carbon nanotube growth. Diameter of metal clusters and number of atoms in these clusters are some of the essential information for predicting carbon nanotube formation and growth, which is then modeled by Boudouard reaction (2CO ---> C(s) + CO2) with metal catalysts. The growth kinetic model is integrated with a two-dimensional axisymmetric reactor flow model to predict reactor performance.

Gokcen, T.↗

NASA NTRS2002

Paralinear Oxidation of Silicon Nitride in a Water Vapor/Oxygen Environment

Three silicon nitride materials were exposed to dry oxygen flowing at 0.44 cm/s at temperatures between 1200 and 1400 C. Reaction kinetics were measured with a continuously recording microbalance. Parabolic kinetics were observed. When the same materials were exposed to a 50% H2O - 50% O2 gas mixture flowing at 4.4 cm/s, all three types exhibited paralinear kinetics. The material is oxidized by water vapor to form solid silica. The protective silica is in turn volatilized by water vapor to form primarily gaseous Si(OH)4. Nonlinear least squares analysis and a paralinear kinetic model were used to determine both parabolic and linear rate constants from the kinetic data. Volatilization of the protective silica scale can result in accelerated consumption of Si3N4. Recession rates under conditions more representative of actual combustors are compared to the furnace data.

Fox, Dennis S.↗

NASA NTRS1995

Initiation-promotion model of tumor prevalence in mice from space radiation exposures

Exposures in space consist of low-level background components from galactic cosmic rays (GCR), occasional intense-energetic solar-particle events, periodic passes through geomagnetic-trapped radiation, and exposure from possible onboard nuclear-propulsion engines. Risk models for astronaut exposure from such diverse components and modalities must be developed to assure adequate protection in future NASA missions. The low-level background exposures (GCR), including relativistic heavy ions (HZE), will be the ultimate limiting factor for astronaut career exposure. We consider herein a two-mutation, initiation-promotion, radiation-carcinogenesis model in mice in which the initiation stage is represented by a linear kinetics model of cellular repair/misrepair, including the track-structure model for heavy ion action cross-sections. The model is validated by comparison with the harderian gland tumor experiments of Alpen et al. for various ion beams. We apply the initiation-promotion model to exposures from galactic cosmic rays, using models of the cosmic-ray environment and heavy ion transport, and consider the effects of the age of the mice prior to and after the exposure and of the length of time in space on predictions of relative risk. Our results indicate that biophysical models of age-dependent radiation hazard will provide a better understanding of GCR risk than models that rely strictly on estimates of the initial slopes of these radiations.

NASA Discipline Radiation Health↗

DOE OSTI2023

From Hydrometeor Size Distribution Measurements to Projections of Wind Turbine Blade Leading-Edge Erosion

Wind turbine blade leading-edge erosion (LEE) is a cause of increased operation and maintenance costs and decreased annual energy production. Thus, detailed, site-specific quantification of likely erosion conditions are critically needed to inform wind plant owner/operator decisions regarding mitigation strategies. Estimating the damage potential at a wind plant site requires accurate measurement of precipitation intensity, phase, droplet size distributions, wind speeds and their joint statistics. The current work quantifies the effect of disdrometer type on the characterization of LEE potential at a site in the US Southern Great Plains. using observations from three co-located disdrometers (an optical, an impact and a video disdrometer), along with hub-height wind-speed observations from a Doppler lidar and two LEE models: a kinetic energy model and the Springer model. Estimates of total kinetic energy of hydrometeor impacts over the four-year study period vary by as much as 38%, and coating lifetime derived from accumulated distance-to-failure estimates from the Springer model differ by an even greater amount, depending on disdrometer type. Damage potential at this site is concentrated in time, with 50% of impact kinetic energy occurring in 6–12 h per year, depending on which set of disdrometer observations is used. Rotor-speed curtailment during the most erosive 0.1–0.2% of 10 min periods is found to increase blade lifetimes and lead to the lowest levelized cost of energy.

17 WIND ENERGY↗

NASA NTRS2005

Low NO(x) Combustor Development

The goal of these efforts was the development of an ultra-low emissions, lean-burn combustor for the High Speed Civil Transport. The HSCT Mach 2.4 FLADE C1 Cycle was selected as the baseline engine cycle. A preliminary compilation of performance requirements for the HSCT combustor system was developed. The emissions goals of the program, baseline engine cycle, and standard combustor performance requirements were considered in developing the compilation of performance requirements. Seven combustor system designs were developed. The development of these system designs was facilitated by the use of spreadsheet-type models which predicted performance of the combustor systems over the entire flight envelope of the HSCT. A chemical kinetic model was developed for an LPP combustor and employed to study NO(x) formation kinetics, and CO burnout. These predictions helped to define the combustor residence time. Five fuel-air mixer concepts were analyzed for use in the combustor system designs. One of the seven system designs, one using the Swirl-Jet and Cyclone Swirler fuel-air mixers, was selected for a preliminary mechanical design study.

Kastl, J. A.↗

NASA NTRS1988

The kinetic approach in magnetospheric plasma transport modeling

The need for a kinetic approach in magnetospheric plasma transport problems is reviewed, as are the trends in its recent applications. The need for kinetic modeling is particularly obvious when confronted with the astonishing variety of magnetospheric particle measurements that display compelling energy and pitch angle-related spatial and/or temporal dispersion, and various types of highly non-Maxwellian features in the distribution functions. Global problems in which the kinetic approach has recently been applied include solar wind plasma injection and dispersion over the cusp, substorm particle injection near synchronous orbit, synergistic energization of ionospheric ions into ring current populations by waves and induced electric field-driven convection, and ionospheric outflow from restricted source regions into the magnetosphere. Kinetic modeling can include efforts ranging from test-particle techniques to particle-in-cell studies, and this range is considered here. There are some areas where fluid and kinetic approaches have been combined or patched together, and these will be briefly discussed.

Horwitz, J. L.↗

DOE OSTI2023

A Comparative Study of Electrical Double Layer Effects for CO Reduction Reaction Kinetics

Solvation models describe how the interactions between the solutes and solvents affect the reactivity and selectivity in electrochemical processes. In this study, we developed a framework for evaluating the effects of applied potential and electrical double layer on CO reduction (COR), comparing fully explicit, implicit, and hybrid solvation models at the standard hydrogen electrode (SHE) scale. We analyzed all crucial intermediates leading to the production of C 1 and C 2 products and found good agreement across these models. Some notable differences were observed in the implicit description of *C and *CO adsorption at higher overpotentials and overall trends in the adsorption energies within the hybrid model. Using this unified SHE framework, we built comparable microkinetic models for COR kinetics and rates. Despite small differences in thermodynamics, solvation-model-based microkinetic simulations showed good agreement for onset potentials against the benchmark experiment. Only qualitative difference was observed for C 1+ versus hydrogen evolution at high overpotentials for the implicit model. Finally, we constructed a generalized C 2 selectivity map in descriptor space (U SHE , ΔG electrolyte CO ), which highlights the limitations of a copper-based COR catalyst and guides the search for optimal descriptor parameters to maximize C 2 selectivity. In conclusion, these findings demonstrate the importance of considering electrical double layer effects in reduction reactions and offer a useful framework for comparing solvation models and predicting optimal electrochemical conditions for specific applications.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

NASA NTRS1976

Combustion modeling in internal combustion engines

The fundamental assumptions of the Blizard and Keck combustion model for internal combustion engines are examined and a generalization of that model is derived. The most significant feature of the model is that it permits the occurrence of unburned hydrocarbons in the thermodynamic-kinetic modeling of exhaust gases. The general formulas are evaluated in two specific cases that are likely to be significant in the applications of the model.

Zeleznik, F. J.↗