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 883 records · Page 49

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↗

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↗

Numerical Analysis of Detonability Assessment in a Natural Gas-Air Fueled Rotating Detonation Engine

Rotating Detonation Engines (RDE) offer an alternative combustion strategy to replace conventional constant pressure combustion with a process that could produce a pressure gain without the use of a mechanical compressor. Recent numerical and experimental publications that consider air as the oxidizer have primarily focused on the ability of these annular combustors to sustain a stable continuous detonation wave when fueled by hydrogen. However, for this to be a viable consideration for the land-based power generation it is necessary to explore the ability to detonate natural gas and air within the confines of the annular geometry of an RDE. Previous studies on confined detonations have expressed the importance of permitting detonation cells to fully form within the combustor in order to achieve stability. This poses a challenge for natural gas–air fueled processes as their detonation cell size can be quite large even at moderate pressures. Despite the practical importance, only a few studies are available on natural gas detonations for air-breathing RDE applications. Moreover, the extreme thermodynamic condition (high temperature inside the combustor) allows limited accessibility inside the combustor for detailed experimental instrumentations, providing mostly single-point data. Recent experimental studies at the National Energy Technology Laboratory (NETL) have reported detonation failure at higher methane concentration in an air-breathing RDE fueled by natural gas-hydrogen fuel blends. This encourages to perform a detailed numerical investigation on the wave characteristics of detonation in a natural gas-air fueled RDE to understand the various aspects of instability associated with the natural gas-air detonation.This study is a numerical consideration of a methane-air fueled RDE with varying operating conditions to ascertain the ability to achieve a stable, continuous detonation wave. The simulations have been performed in a 2D unwrapped RDE geometry using the open-source CFD library “OpenFOAM” employing an unsteady pressure-based compressible reactive flow solver with a <em>k</em>–<em>ε</em> turbulence model in a structured rectangular grid system. Both reduced and detailed chemical kinetic models have been used to assess the effect of the chemistry on the detonation wave characteristics and the underlying flow features. A systematic grid sensitivity study has been conducted with various grid sizes to quantify the weakly stable overdriven detonation on a coarse mesh and oscillating features at fine mesh resolutions.The main focus of the current study is to investigate the effects of operating injection pressure on detonation wave characteristics of an air-breathing Rotating Detonation Engine (RDE) fueled with natural gas-hydrogen fuel blends. Wave speeds, peak pressures and temperatures, and dominant frequencies have been computed from the time histories. The flow structures were then visualized using 2D contours of temperature and species concentration.

Saha, Pankaj↗

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↗

The multi-assumption architecture and testbed (MAAT v1.0): R code for generating ensembles with dynamic model structure and analysis of epistemic uncertainty from multiple sources

Abstract. Here, computer models are ubiquitous tools used to represent systems across many scientific and engineering domains. For any given system, many computer models exist, each built on different assumptions and demonstrating variability in the ways in which these systems can be represented. This variability is known as epistemic uncertainty, i.e. uncertainty in our knowledge of how these systems operate. Two primary sources of epistemic uncertainty are (1) uncertain parameter values and (2) uncertain mathematical representations of the processes that comprise the system. Many formal methods exist to analyse parameter-based epistemic uncertainty, while process-representation-based epistemic uncertainty is often analysed post hoc, incompletely, informally, or is ignored. In this model description paper we present the multi-assumption architecture and testbed (MAAT v1.0) designed to formally and completely analyse process-representation-based epistemic uncertainty. MAAT is a modular modelling code that can simply and efficiently vary model structure (process representation), allowing for the generation and running of large model ensembles that vary in process representation, parameters, parameter values, and environmental conditions during a single execution of the code. MAAT v1.0 approaches epistemic uncertainty through sensitivity analysis, assigning variability in model output to processes (process representation and parameters) or to individual parameters. In this model description paper we describe MAAT and, by using a simple groundwater model example, verify that the sensitivity analysis algorithms have been correctly implemented. The main system model currently coded in MAAT is a unified, leaf-scale enzyme kinetic model of C 3 photosynthesis. In the Appendix we describe the photosynthesis model and the unification of multiple representations of photosynthetic processes. The numerical solution to leaf-scale photosynthesis is verified and examples of process variability in temperature response functions are provided. For rapid application to new systems, the MAAT algorithms for efficient variation of model structure and sensitivity analysis are agnostic of the specific system model employed. Therefore MAAT provides a tool for the development of novel or toy models in many domains, i.e. not only photosynthesis, facilitating rapid informal and formal comparison of alternative modelling approaches.

54 ENVIRONMENTAL SCIENCES↗

Temporal Analysis Of Product Computational Methods

The temporal analysis of products (TAP) reactor provides a wealth of data for transient studies of gas phase heterogeneous reactions, with the ability to provide a comprehensive picture of intrinsic reaction kinetics pertaining to complex catalytic mechanisms. A novel framework called the Y-Procedure was previously developed by Yablonsky and coworkers to extract information from this non-steady-state kinetic data. The Y-Procedure reconstructs the chemical transformation rate and gas composition in the catalytic zone with no a priori assumptions on the kinetic model and is used to relate changes in catalytic reaction rates and kinetic parameters to changes in a catalyst’s surface composition. This code provides the required steps to pre-process, analyze, and perform the Y-Procedure on experiments derived from the TAP reactor.

Kunz, Matthew R.↗

Verification of electromagnetic fluid-kinetic hybrid electron model in global gyrokinetic particle simulation

The fluid-kinetic hybrid electron model is verified in global gyrokinetic particle simulation of linear electromagnetic drift-Alfvenic instabilities in tokamak. In particular, we have recovered the {beta}-stabilization of the ion temperature gradient mode, transition to collisionless trapped electron mode, and the onset of kinetic ballooning mode as {beta}{sub e} (ratio of electron kinetic pressure to magnetic pressure) increases.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Effect of high energy electrons on H{sup −} production and destruction in a high current DC negative ion source for cyclotron

Recently, a filament driven multi-cusp negative ion source has been developed for proton cyclotrons in medical applications. In this study, numerical modeling of the filament arc-discharge source plasma has been done with kinetic modeling of electrons in the ion source plasmas by the multi-cusp arc-discharge code and zero dimensional rate equations for hydrogen molecules and negative ions. In this paper, main focus is placed on the effects of the arc-discharge power on the electron energy distribution function and the resultant H{sup −} production. The modelling results reasonably explains the dependence of the H{sup −} extraction current on the arc-discharge power in the experiments.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Kinetics of Recycling CF/EP Composites by Thermal Excitation of Cr{sub 2}O{sub 3}

Recycling high performance carbon fibers (CFs) by thermal excitation of Cr{sub 2}O{sub 3} semiconductor was proposed. The kinetic model of decomposition of epoxy resin was established and the method for solving kinetic parameters was proposed. The influences of temperature and treatment time on the decomposition rate of epoxy resin were investigated. The kinetic equation was established and the accuracy was verified experimentally. The morphology of the recycled CF was observed and monofilament tensile strength was determined. Results indicated that temperature and treatment time had positive correlation with decomposition rate of the epoxy resin. The reaction order n was 1.5, with excitation energy E of 85.03 kJ·mol{sup −1} and pre-exponential factor k{sub 0} of 1.26 × 10{sup 6} min{sup −1}. The established kinetic equation could be used to estimate the temperature and treatment time. Prolongation of treatment time could improve the monofilament tensile strength. The monofilament tensile strength of the recycled CF at 480 °C for 30 min was 102.72% of that of the original CF, indicating that surface heat treatment was performed on CF during the recycling process.

36 MATERIALS SCIENCE↗

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↗

Indirect determination of the electric field in plasma discharges using laser-induced fluorescence spectroscopy

The evaluation of electric fields is of prime interest for the description of plasma characteristics. In this work, different methods for determining the electric field profile in low-pressure discharges using one- and two-dimensional Laser-Induced Fluorescence (LIF) measurements are presented and discussed. The energy conservation, fluid, and kinetic approaches appear to be well-suited for the electric field evaluation in this region of the plasma flow. However, the numerical complexity of a two-dimensional kinetic model is penalizing due to the limited signal-to-noise ratio that can be achieved, making the computation of the electric field subject to large error bars. The ionization contribution which appears in the fluid model makes it unattractive on an experimental viewpoint. The energy conservation and 1D1V kinetic approaches should therefore be preferred for the determination of the electric field when LIF data are used.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

New risk metrics and mathematical tools for risk analysis: Current and future challenges

The current status of the food safety supply world wide, has led Food and Agriculture Organization (FAO) and World Health Organization (WHO) to establishing Risk Analysis as the single framework for building food safety control programs. A series of guidelines and reports that detail out the various steps in Risk Analysis, namely Risk Management, Risk Assessment and Risk Communication is available. The Risk Analysis approach enables integration between operational food management systems, such as Hazard Analysis Critical Control Points, public health and governmental decisions. To do that, a series of new Risk Metrics has been established as follows: i) the Appropriate Level of Protection (ALOP), which indicates the maximum numbers of illnesses in a population per annum, defined by quantitative risk assessments, and used to establish; ii) Food Safety Objective (FSO), which sets the maximum frequency and/or concentration of a hazard in a food at the time of consumption that provides or contributes to the ALOP. Given that ALOP is rather a metric of the public health tolerable burden (it addresses the total ‘failure’ that may be handled at a national level), it is difficult to be interpreted into control measures applied at the manufacturing level. Thus, a series of specific objectives and criteria for performance of individual processes and products have been established, all of them assisting in the achievement of FSO and hence, ALOP. In order to achieve FSO, tools quantifying the effect of processes and intrinsic properties of foods on survival and growth of pathogens are essential. In this context, predictive microbiology and risk assessment have offered an important assistance to Food Safety Management. Predictive modelling is the basis of exposure assessment and the development of stochastic and kinetic models, which are also available in the form of Web-based applications, e.g., COMBASE and Microbial Responses Viewer), or introduced into user-friendly softwares, (e.g., Seafood Spoilage Predictor) have evolved the use of information systems in the food safety management. Such tools are updateable with new food-pathogen specific models containing cardinal parameters and multiple dependent variables, including plate counts, concentration of metabolic products, or even expression levels of certain genes. Then, these tools may further serve as decision-support tools which may assist in product logistics, based on their scientifically-based and “momentary” expressed spoilage and safety level.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Scalar-Filtered Mass-Density-Function Simulation of Swirling Reacting Flows on Unstructured Grids

The scalar-filtered mass-density-function methodology is employed for large-eddy simulation of two swirl-stabilized nonpremixed flames. These are low-swirl (SM1) and high-swirl (SMA2) Sydney methane flames, both of which have been the subject of detailed laboratory measurements. Combustion chemistry is modeled via a flamelet model for the low-swirl flame, and a detailed finite-rate kinetics model for the high-swirl flame. The scalar-filtered mass-density-function is simulated by a Lagrangian Monte Carlo method on a domain discretized by unstructured grids. The simulated results are assessed via comparison with experimental data and show very good agreement. This demonstrates the capability of scalar-filtered mass-density-function for large-eddy simulation of complex flows and warrants future applications of the methodology for large-eddy simulation of practical combustor configurations.

42 ENGINEERING↗

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↗

Characterization and modelling of the boron-oxygen defect activation in compensated n-type silicon

A study of the activation of the light-induced degradation in compensated n-type Czochralski grown silicon is presented. A kinetic model is established that verifies the existence of both the fast and the slow components known from p-type and proves the quadratic dependence of the defect generation rates of both defects on the hole concentration. The model allows for the description of lifetime degradation kinetics in compensated n-type silicon under various intensities and is in accordance with the findings for p-type silicon. We found that the final concentrations of the slow defect component in compensated n-type silicon only depend on the interstitial oxygen concentration and on neither the boron concentration nor the equilibrium electron concentration n{sub 0}. The final concentrations of the fast defect component slightly increase with increasing boron concentration. The results on n-type silicon give new insight to the origin of the BO defect and question the existing models for the defect composition.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Toward the Development of a Fundamentally Based Chemical Model for Cyclopentanone: High-Pressure-Limit Rate Constants for H Atom Abstraction and Fuel Radical Decomposition

Theoretical aspects of the development of a chemical kinetic model for the pyrolysis and combustion of a cyclic ketone, cyclopentanone, are considered. We present calculated thermodynamic and kinetic data for the first time for the principal species including 2- and 3-oxo-cyclopentyl radicals, which are in reasonable agreement with the literature. Furthermore, these radicals can be formed via H atom abstraction reactions by H and Ö atoms and OH, HO 2 , and CH 3 radicals, the rate constants of which have been calculated. Abstraction from the β-hydrogen atom is the dominant process when OH is involved, but the reverse holds true for HO 2 radicals. We also determined the subsequent β-scission of the radicals formed, and it is shown that recent tunable VUV photoionization mass spectrometry experiments can be interpreted in this light. The bulk of the calculations used the composite model chemistry G4, which was benchmarked in the simplest case with a coupled cluster treatment, CCSD(T), in the complete basis set limit.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

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↗

The Degree of Rate Control: A Powerful Tool for Catalysis Research

The “degree of rate control” (DRC) is a mathematical approach for analysing multi-step reaction mechanisms that has proven very useful in catalysis research. It identifies the “ratecontrolling transition states and intermediates” (i.e., those whose DRCs are large in magnitude). Even in mechanisms with over 30 intermediates and transition states, these are generally just a few distinct chemical species whose energies, if they could be independently changed, would achieve a faster net reaction rate to the product of interest. For example, when there is a single “rate-determining step”, the DRC for its transition state (TS) is 1, which means (by definition) that if this TS’s energy could be decreased by kBT (where kB is Boltzmann’s constant and T is temperature), the net rate would increase by a factor of e. Since the (relative) energies of these key adsorbed intermediates and transition states can be adjusted by modifying the catalyst or solvent, or even a reactant’s molecular structure, the DRC values provide important ideas for catalyst improvement. The species with large DRCs are also the ones whose energetics must be most accurately measured or calculated to achieve an accurate kinetic model for any reaction mechanism. A tutorial on DRC analysis, the calculation of DRCs, and examples of the applications of DRCs in catalysis research is presented here. Applications of DRC analysis include: clarifying reaction kinetics, improving the accuracy of computational models, improving reaction conditions, improving choice of oxidant in selective oxidation, incorporation in algorithms which calculate net reaction rates of multistep mechanisms without solving the differential equations involved, and high-throughput computational screening of catalyst materials. Since DRC values can be determined experimentally, we conclude a full microkinetic model is not required to take advantage of DRC analysis.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗