Search NASA⌕ Search

SEARCH · Search NASA

Results for “Kinetic Modelling”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

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

At least 1,081 records · Page 60
All sourcesNASA NTRSDOE OSTI
DOE OSTI2016

Heliospheric current sheet and effects of its interaction with solar cosmic rays

The effects of interaction of solar cosmic rays (SCRs) with the heliospheric current sheet (HCS) in the solar wind are analyzed. A self-consistent kinetic model of the HCS is developed in which ions with quasiadiabatic dynamics can present. The HCS is considered an equilibrium embedded current structure in which two main plasma species with different temperatures (the low-energy background plasma of the solar wind and the higher energy SCR component) contribute to the current. The obtained results are verified by comparing with the results of numerical simulations based on solving equations of motion by the particle tracing method in the given HCS magnetic field with allowance for SCR particles. It is shown that the HCS is a relatively thin multiscale current configuration embedded in a thicker plasma layer. In this case, as a rule, the shear (tangential to the sheet current) component of the magnetic field is present in the HCS. Taking into account high-energy SCR particles in the HCS can lead to a change of its configuration and the formation of a multiscale embedded structure. Parametric family of solutions is considered in which the current balance in the HCS is provided at different SCR temperatures and different densities of the high-energy plasma. The SCR densities are determined at which an appreciable (detectable by satellites) HCS thickening can occur. Possible applications of this modeling to explain experimental observations are discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Fuel Properties and Kinetics

The Co-Optimization of Fuels and Engines (CoOptima) program is focused on a two-pronged approach to develop new, lower net-carbon fuels and advanced engine technologies in tandem. This effort brings nine national labs and several university and industry partners together to simultaneously identify new fuel candidates and to further engine technology development in the light-, medium-, and heavy-duty vehicle sectors with a focus on fuel economy and emission reduction. The fuel properties and kinetics project is aimed at achieving a better understanding of fuel property kinetics and how they impact engine performance such that more efficient and lower emitting engines can be made a reality. Several promising biofuel candidates were identified through tiered screening approaches and were evaluated utilizing unique bench-scale measurement techniques. Kinetic models were developed to describe autoignition and soot pre-cursor formation mechanisms to identify the most promising candidates that advanced various combustion strategies.

ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2023

Kinetics and extinction of non-premixed cool and warm flames of dimethyl ether at elevated pressure

Here, the growing demand of clean and efficient propulsion and energy systems has sparked an interest in understanding low-temperature combustion at high pressure. Cool flame transition and extinction limits as well as oxygen concentration dependence at elevated pressures provide insights of the low-temperature and high-pressure fuel reactivity. A new experimental high-pressure counterflow burner platform was designed and developed to achieve the studies of high-pressure cool flames. Dimethyl ether (DME) was chosen to study its non-premixed cool flame in high-pressure counterflow burner at pressure up to 5 atm, perhaps for the first time. This paper investigates the effects of pressure on cool flame structure, extinction and transition limits, and oxygen concentration dependence as well as ozone assisted warm flames of DME in experiments and numerical simulations. The results show that the reignition transition from cool flame to hot flame occurs either with the decrease of the strain rate at a given fuel concentration and pressure or with the increase of fuel mole fraction or pressure at a given strain rate. Furthermore, it is shown that the higher pressure shifts the cool flame to higher strain rates and results in higher cool flame extinction strain rates. However, the existing kinetic model of DME fails in predicting the cool flame extinction limit at elevated pressures. Besides, the cool flame extinction limits are proportional to nth power of the oxygen concentration, [O 2 ] n , and the increase of pressure leads to stronger extinction limit dependence (larger n) on oxygen concentration. The present experiment and detailed kinetic analysis show clearly that increasing pressure promotes the low-temperature chemistry including the oxygen addition reactions. In addition, stable warm flame was first experimentally observed by using DME at elevated pressure with ozone sensitization.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2022

The role of collisional energy transfer in the thermal and prompt dissociation of 1-methyl allyl

Collisional energy transfer plays a pivotal role in a-priori calculations of the pressure-dependent rate constants obtained from a master equation. However, accurate determinations of collisional energy transfer parameters are rare and so most kinetic studies rely on best-fits of these parameters to experimental measurements (if available) or estimates guided from literature studies. In this work, we have quantified the effect of the uncertainty in energy transfer parameters on the thermal and prompt dissociation kinetics of a resonance stabilized radical, 1-methyl allyl (1MA), of relevance to the combustion of 1- and 2-butene isomers. Simulations using literature kinetics models were performed to assess the impact of these uncertainties on flame propagation and speciation data in laminar flames of 1- and 2-butene. Analyses of the uncertainty propagated by the energy transfer parameters for 1MA dissociation to the flame simulations, in particular the laminar flame speed, indicate an intricate coupling between the kinetics of 1MA dissociation (chain propagation) and its reaction with H-atoms (chain-termination). Ab-initio based theoretical calculations were also performed to obtain pressure-dependent kinetics for the reaction of 1MA with H-atoms. Lastly, theoretically calculated energy transfer parameters were used to best characterize the kinetics and branching between the chain propagating 1MA dissociation and the chain-terminating reaction 1MA + H.

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

Coarse-Grained Modeling Elucidates Differential Metabolism of Saccharomyces cerevisiae under Varied Nutrient Limitations

Microorganisms such as Saccharomyces cerevisiae have a native ability to adapt their metabolism to varying nutrient conditions. Understanding their responses to nutrient limitations is critical for decoding cellular physiology and designing strategies for metabolic engineering. While the influence of carbon availability on yeast metabolism has been extensively studied, the role of nitrogen availability remains relatively underexplored. In this study, we utilized a coarse-grained kinetic model to systematically analyze and compare the effects of carbon and nitrogen limitations on yeast metabolism. Our model successfully revealed the differential metabolic characteristics of S. cerevisiae under carbon- and nitrogen-limited chemostat conditions. It also highlighted the significance of protein activity regulation at varying carbon-to-nitrogen ratios, and elucidated distinct strategies employed to maintain ATP homeostasis. Furthermore, this study provides a computational tool for investigating yeast physiology under nutrient limitations and offers quantitative and mechanistic insights into yeast metabolism.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2017

System of kinetic equations describing large-scale processes in collisionless space plasma

A system of kinetic equations describing relatively slow large-scale processes in collisionless magnetoplasma structures with a spatial resolution of about the characteristic gyroradius is derived. Plasma is assumed to be quasineutral, while the magnetic and electric fields are determined by the instantaneous distributions of the particle and current densities and the stress tensor of all plasma components in the longrange instantaneous interaction approximation. A special version of equations is derived for the case of magnetized electrons described by the Vlasov equation in the drift approximation. The obtained system of equations can be used to develop a global numerical kinetic model of the Earth’s magnetosphere with a spatial resolution of about 100 km, as well as local models of certain regions of the Earth’s magnetosphere with a higher resolution.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2015

Partnership for Edge Physics Simulation (EPSI)

We propose to develop advanced simulation codes, based upon an extreme parallelism, first principles kinetic approach, to address the challenges associated with the edge region of magnetically confined plasmas. This work is relevant to both existing magnetic fusion facilities and essential for next-generation burning plasma experiments, such as ITER where success is critically dependent upon H-mode operation achieving an edge pedestal of sufficient height for good core plasma performance without producing deleterious large scale edge localized instabilities. The plasma edge presents a well-known set of multi-physics, multi-scale problems involving complex 3D magnetic geometry. Perhaps the greatest computational challenge is the lack of scale separation – temporal scales for drift waves, Alfven waves, ELM dynamics for example have strong overlap. Similar overlap occurs on the spatial scales for the ion poloidal gyro-radius, drift wave and pedestal width. The traditional approach of separating fusion problems into weakly interacting spatial or temporal domains clearly breaks down in the edge. A full kinetic model (full-f model) must be solved to understand and predict the edge physics including non-equilibrium thermodynamic issues arising from the magnetic topology (the open field lines producing a spatially sensitive velocity hole), plasma wall interactions, neutral and atomic physics. The plan here is to model these phenomena within a comprehensive first principles set of equations without the need for the insurmountable multiple-codes coupling issues by building on the XGC1 code developed under the SciDAC Proto-FSP Center for Plasma Edge Simulation (CPES). This proposal includes the critical participants in the XGC1 development. We propose enhancing the capability of XGC1 by including all the important turbulence physics contained in kinetic ion and electron electromagnetic dynamics, by extending the PIC technology to incorporate several positive features found in the Eulerian grid technology, combining the coarse grained XGC0 code into it with a multiscale projective time integration technique to be developed in the project for experimental time scale simulation of the edge multi-physics, and by implementing modern computational technologies to move forward with the extreme scale heterogeneous hardware/software platform development. Ultra-fast edge localized instabilities, involving compressional Alven modes, will be studied using a proper MHD or two-fluid codes employing kinetic closure information from XGC. The proposed research integrates the most important physical processes on overlapping temporal and spatial scales, to study (i) transitions and thresholds for enhanced confinement regimes, (ii) the predictive understanding of the edge pedestal formation, structure and dynamics, (iii) the effect of the edge plasma on the core confinement, (iv) the physics of Edge Localized Modes and their suppression or mitigation via external control techniques, and (v) the heat load on the material wall. Our approach will implement physical models that are valid in their relevant collisionality regimes for current experiments and future burning plasma devices. Comparisons with experimental data will be made at multiple levels including characterization of turbulence dynamics, flux-gradient responses in the energy, particle and momentum channels, and observations of global phenomena such as L-H transitions, momentum generation, ELM dynamics, ELM-free regimes and the tokamak density limit. The proposal team has established strong connections, through joint membership, with other organizations that will be critical for success. These include FASTMath, SUPER and QUEST Institutes respectively for meshing tools and solvers; code performance tuning and uncertainty quantification as well as a proposed Data Management and Visualization Institute.. We will also be teaming with ongoing fusion SciDAC’s: CSWPI for RF-edge plasmas interactions; CEMM for the M3D MHD simulations and with PEPSC for incorporation of energetic particle physics. Finally, the team includes representatives from each of the major U.S. experimental facilities who will coordinate the validation program and the development of synthetic diagnostics.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Scalable Implicit Solvers with Dynamic Mesh Adaptation for a Relativistic Drift-Kinetic Fokker–Planck–Boltzmann Model

In this work we consider a relativistic drift-kinetic model for runaway electrons along with a Fokker–Planck operator for small-angle Coulomb collisions, a radiation damping operator, and a secondary knock-on (Boltzmann) collision source. Here, we develop a new scalable fully implicit solver utilizing finite volume and conservative finite difference schemes and dynamic mesh adaptivity. A new data management framework in the PETSc library based on the p4est library is developed to enable simulations with dynamic adaptive mesh refinement (AMR), distributed memory parallelization, and dynamic load balancing of computational work. This framework and the runaway electron solver building on the framework are able to dynamically capture both bulk Maxwellian at the low-energy region and a runaway tail at the high-energy region. To effectively capture features via the AMR algorithm, a new AMR indicator prediction strategy is proposed that is performed alongside the implicit time evolution of the solution. This strategy is complemented by the introduction of computationally cheap feature-based AMR indicators that are analyzed theoretically. Numerical results quantify the advantages of the prediction strategy in better capturing features compared with nonpredictive strategies; and we demonstrate trade-offs regarding computational costs. The robustness with respect to model parameters, algorithmic scalability, and parallel scalability are demonstrated through several benchmark problems including manufactured solutions and solutions of different physics models. We focus on demonstrating the advantages of using implicit time stepping and AMR for runaway electron simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Multiscale Modeling of a Direct Nonoxidative Methane Dehydroaromatization Reactor with a Validated Model for Catalyst Deactivation

Due to the recent boom in shale gas production, aromatics production using direct nonoxidative methane dehydroaromatization (DHA) is being investigated extensively. However, due to rapid coke formation, catalysts in the nonoxidative methane DHA reactors get deactivated, which is one of the critical issues for the commercial success of the methane DHA process. In this paper, a model for catalyst deactivation is developed. Rate models for other DHA reactions are developed by considering the decrease in the catalyst activity with time. Due to the very fast coke formation rate on the fresh catalyst, there is coke formation immediately upon the introduction of the feed. Therefore, an algorithm is developed for estimation of the initial state of the reactor and the kinetic parameters by coupling an iterative direct substitution approach with an optimization approach. Transient experimental data from an in-house reactor are first reconciled and then used for developing the kinetic model including the coke formation model. Using the rate model, a dynamic, heterogeneous, multiscale reactor model with embedded heating is developed. Here, the model couples the catalyst pellet level model with a reactor level model. Impacts of temperature, L/D ratio, and scheduling of reactors on variability in conversion and yield with time are studied.

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

Origin of Multiple Peaks in the Potentiodynamic Oxidation of CO Adlayers on Pt and Ru-Modified Pt Electrodes

The research of the electrooxidation mechanism of CO ad on Pt based catalysts is very important for designing more effective CO-tolerant electrocatalysts for fuel cells. We have studied the origin of multiple peaks in the cyclic voltammograms of CO stripping from polycrystalline Pt and Ru modified polycrystalline Pt (Pt/Ru) surfaces in both acidic and alkaline media by differential electrochemical mass spectrometry (DEMS), DFT calculations, and kinetic Monte Carlo (KMC) simulations. A novel CO ad electrooxidation kinetic model on heterogeneous Pt and Pt/Ru catalysts is proposed to account for the multiple peaks experimentally observed. In this model, OH species prefer to adsorb at low-coordination sites or Ru sites and, thus, suppress CO repopulation from high-coordination sites onto these sites. Therefore, CO ad oxidation occurs on different facets or regions, leading to multiplicity of CO stripping peaks. This work provides a new insight into the CO electrooxidation mechanism and kinetics on heterogeneous catalysts.

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

NMR spin-rotation relaxation and diffusion of methane

The translational diffusion-coefficient and the spin-rotation contribution to the 1 H NMR relaxation rate for methane (CH 4 ) are investigated using MD (molecular dynamics) simulations, over a wide range of densities and temperatures, spanning the liquid, supercritical, and gas phases. Here, the simulated diffusion-coefficients agree well with measurements, without any adjustable parameters in the interpretation of the simulations. A minimization technique is developed to compute the angular velocity for non-rigid spherical molecules, which is used to simulate the autocorrelation function for spin-rotation interactions. With increasing diffusivity, the autocorrelation function shows increasing deviations from the single-exponential decay predicted by the Langevin theory for rigid spheres, and the deviations are quantified using inverse Laplace transforms. The 1 H spin-rotation relaxation rate derived from the autocorrelation function using the “kinetic model” agrees well with measurements in the supercritical/gas phase, while the relaxation rate derived using the “diffusion model” agrees well with measurements in the liquid phase. 1 H spin-rotation relaxation is shown to dominate over the MD-simulated 1 H- 1 H dipole-dipole relaxation at high diffusivity, while the opposite is found at low diffusivity. At high diffusivity, the simulated spin-rotation correlation time agrees with the kinetic collision time for gases, which is used to derive a new expression for 1 H spin-rotation relaxation, without any adjustable parameters.

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

Modeling Functional Organic Chemistry in Arctic Rivers: An Idealized Siberian System

Rivers of the Arctic will become ever more important for the global climate, since they carry a majority of continental dissolved organic carbon flux into the rapidly changing polar ocean. Aqueous organics comprise a wide array of functional groups, several of which are likely to impact coastal and open water biophysical properties. Light attenuation, interfacial films, aerosol formation, gas release and momentum exchange can all be cited. We performed Lagrangian kinetic modeling for the evolution of riverine organic chemistry as the molecules in question make their way from the highlands to Arctic outlets. Classes as diverse as the proteins, sugars, lipids, re-condensates, humics, bio-tracers and small volatiles are all included. Our reduced framework constitutes an idealized northward flow driving a major hydrological discharge rate and primarily representing the Russian Lena. Mountainous, high solute and tundra sources are all simulated, and they meet up at several points between soil and delta process reactors. Turnover rates are parameterized beginning with extrapolated coastal values imposed along a limited tributary network, with connections between different terrestrial sub-ecologies. Temporal variation of our total dissolved matter most closely resembles the observations when we focus on the restricted removal and low initial carbon loads, suggesting relatively slow transformation along the water course. Thus, channel combinations and mixing must play a dominant role. Nevertheless, microbial and photochemical losses help determine the final concentrations for most species. Chemical evolution is distinct for the various functionalities, with special contributions from pre- and post-reactivity in soil and delta waters. Several functions are combined linearly to represent the collective chromophoric dissolved matter, characterized here by its absorption. Tributaries carry the signature of lignin phenols to segregate tundra versus taiga sources, and special attention is paid to the early then marine behaviors of low molecular weight volatiles. Heteropolycondensates comprise the largest percentage of reactive carbon in our simulations due to recombination/accumulation, and they tend to be preeminent at the mouth. Outlet concentrations of individual structures such as amino acids and absorbers lie above threshold values for biophysical influence, on the monolayer and light attenuation. The extent of coastal spreading is examined through targeted regional box modeling, relying on salinity and color for calibration. In some cases, plumes reach the scale of peripheral arctic seas, and amplification is expected during upcoming decades. Conclusions are mapped from the Lena to other boreal discharges, and future research questions are outlined regarding the bonding type versus mass release as permafrost degrades. Dynamic aqueous organic coupling is recommended for polar system models, from headwaters to coastal diluent.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2016

Measurement of Kinetic Parameters by Oscillated Experiments

The availability of reliable values of the kinetic parameters (effective delayed neutron fraction {beta} and effective mean neutron generation time {Lambda}) is important for reactor kinetic analyses and for reactivity monitoring when adopting the point kinetic approach. The kinetic parameters may be obtained either by simulations (Monte Carlo or deterministic) of the system or by direct measurements. Several methods have been proposed for the direct experimental determination of the kinetic parameters, e.g. control rod drops and noise analysis approaches, e.g. In this work, the problem of the determination of kinetic parameters in subcritical systems is tackled. The system is subjected to an oscillated source and the value of the ratio {beta}/{Lambda} is related to a measurable quantity when the system response reaches an equilibrium condition. In the analysis presented, a pure frequency is considered: when analyzing more complicated source periodic behaviours, a Fourier analysis can allow to isolate the contributions of individual frequencies. The phase difference between the oscillating source and the system response is adopted as the experimental observable. The method for the interpretation of the measurements is based on the standard point kinetic model. The response of the system is obtained by localized neutron detectors and, therefore, spatial and spectral effects may also affect the results. This drawback needs some further discussion and will be considered in the future work. The phase shift between the source signal and the system response may be analytically retrieved by Laplace transforming the point kinetic equations in the presence of an oscillated source. An explicit formula relates the phase shift to the ratio {beta}/{Lambda} and to the frequency of the source. Based on this relationship, the ratio {beta}/{Lambda} can be determined. To validate the technique, the experimental data are simulated by a numerical code solving the time-dependent diffusion equations for a source-driven system. The work allows to assess the possibility of the technique and its limits of applicability. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2020

Protein folding from heterogeneous unfolded state revealed by time-resolved X-ray solution scattering

One of the most challenging tasks in biological science is to understand how a protein folds. In theoretical studies, the hypothesis adopting a funnel-like free-energy landscape has been recognized as a prominent scheme for explaining protein folding in views of both internal energy and conformational heterogeneity of a protein. Despite numerous experimental efforts, however, comprehensively studying protein folding with respect to its global conformational changes in conjunction with the heterogeneity has been elusive. Here we investigate the redox-coupled folding dynamics of equine heart cytochrome c (cyt-c) induced by external electron injection by using time-resolved X-ray solution scattering. A systematic kinetic analysis unveils a kinetic model for its folding with a stretched exponential behavior during the transition toward the folded state. With the aid of the ensemble optimization method combined with molecular dynamics simulations, we found that during the folding the heterogeneously populated ensemble of the unfolded state is converted to a narrowly populated ensemble of folded conformations. These observations obtained from the kinetic and the structural analyses of X-ray scattering data reveal that the folding dynamics of cyt-c accompanies many parallel pathways associated with the heterogeneously populated ensemble of unfolded conformations, resulting in the stretched exponential kinetics at room temperature. This finding provides direct evidence with a view to microscopic protein conformations that the cyt-c folding initiates from a highly heterogeneous unfolded state, passes through still diverse intermediate structures, and reaches structural homogeneity by arriving at the folded state.

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

Understanding fundamental effects of biofuel structure on ignition and physical fuel properties

The development of biochemicals and biofuels with advantaged properties of higher combustion efficiency at lower emission levels can be aided by a priori prediction of critical global characteristics based on molecular structure. Existing predictive techniques are incomplete for comprehensive ignition and physical fuel property prediction. Here this investigation focuses on filling this knowledge gap by developing a priori prediction techniques for predicting important ignition as well as physical fuel properties of promising advantaged biofuels. We demonstrate the utility of this a priori prediction technique on saturated (≤C 5 ) alcohols, both linear and branched in nature by correlating their molecular structure with both important ignition and physical properties of the fuels. This work utilizes a semi-automated method of predicting ignition characteristics via accurate fundamental-based calculations of the fuel radical cascades that govern low-temperature combustion. This numerical methodology is supported by the relevant fuel ignition delay time (IDT) and research octane number (RON) data experimentally measured by the Advanced Fuel Ignition Delay Analyzer (AFIDA). Phyiscal fuel property trends are also analyzed by correlating the molecular structure of the fuel to its physical characteristics. Based on their molecular structure, recommendations are made on an a priori method for selecting biofuels that have advantaged properties for the selected fuel application. Lastly, the utility of this investigation is extended beyond the subset of alcohols considered by analyzing the IDT trends of i-Pentanol and n-Hexanol, revealing the fidelity of the novel kinetic modelling approach and the accuracy of the predicted trends in this study.

09 BIOMASS FUELS↗
DOE OSTI2023

Establishing a Design Strategy for Corrosion Resistant Structural Materials in Molten Salt Technologies

To address the gap in material design strategies for structural components in molten salt-powered technologies, the current work aimed to identify key parameters governing corrosion induced degradation of materials during exposures in molten halide salts. The applicability of the data on dissolution rates of pure metals (e.g., Cr) in a given molten salt to enable predictions of corrosion behavior of multicomponent binary and ternary Fe-based model and Ni-based commercial alloys isothermally exposed in the purified KCl-MgCl 2 in quartz capsules at 700 and 800°C was evaluated. The influence of alloy chemistry on initial dissolution rates and the consequent impact on the time-dependent corrosion induced microstructural evolution in the alloy subsurface was predicted with a coupled thermodynamic-kinetic model. Finally, the model was well able to predict the Cr depletion in all investigated alloys demonstrating the ability of the proposed approach to enable a unique design strategy for structural materials in molten salts.

36 MATERIALS SCIENCE↗