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,027 records · Page 57
All sourcesNASA NTRSDOE OSTI
NASA NTRS2000

Study of Ring Current Dynamics During Geomagnetic Storms

This research program considered modeling the dynamical evolution of the ring current during several geomagnetic storms. The first year (6/01/1997-5/31/1998) of this successful collaborative research between the University of New Hampshire (UNH) and the University of California Los Angeles (UCLA) was supported by NASA grant NAG5-4680. The second and third years (6/01/1998-5/31/2000) were funded at UNH under NASA grant NAG5-7368. Research work at UNH concentrated on further development of a kinetic model to treat all of the important physical processes that affect the ring current ion population during storm conditions. This model was applied to simulate ring current development during several International Solar-Terrestrial Physics (ISTP) events, and the results were directly compared to satellite observations. A brief description of our major accomplishments and a list of the publications and presentations resulting from this effort are given.

Jordanova, Vania K.↗
DOE OSTI2021

Probing intermediate temperature heat release in autoignition of C3-C4 iso-alcohol/gasoline blends

In this work, we report an experimental and modeling study on the blending effects of C3–C4 iso-alcohols, namely iso-propanol and iso-butanol, on the characteristics of intermediate temperature heat release (ITHR) for a research-grade gasoline (FACE-F) in a rapid compression machine operating at diluted/stoichiometric fuel loading, compressed pressure of 43 bar, and compressed temperatures from 700 to 970 K. Changes in ITHR behavior are characterized through ITHR extent and evolution at 0 to 30 vol% iso-alcohol blending levels. Experimental observations reveal the strong promoting and decelerating effects on ITHR extent and evolution, respectively, within the low-temperature regime for both iso-alcohols, with stronger effects seen for iso-propanol; within the intermediate-temperature regime, the influence of both iso-alcohols diminishes greatly. Comprehensive chemical kinetic modeling is undertaken using a recently developed gasoline surrogate/alcohol model in conjunction with a five-component gasoline surrogate (FGF-LLNL), with good agreement obtained with the experiments for all the blends. Sensitivity, rate of production and flux analyses highlight the importance of chemical kinetic interactions via both fuel-specific and non-fuel-specific reactions. Disabling all the chemical kinetic interactions between iso-alcohol and FGF-LLNL in carbonated/non-fuel-specific species leads to somewhat extended ITHR duration and reduced ignition reactivity. Direct comparison between FGF-LLNL/iso-propanol and FGF-LLNL/iso-butanol blends is also made, where iso-propanol and iso-butanol are found to influence ITHR characteristics via significantly different chemical kinetic interactions.

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

Mathematical analysis of laboratory microbial experiments demonstrating deterministic chaotic dynamics

Presented is a system of four ordinary differential equations and a mathematical analysis of microbiological experiments in a four-component chemostat—nutrient n, rods r, cocci c, and predators p. The analysis is consistent with the conclusion that previous experiments produced features of deterministic chaotic and classical dynamics depending on dilution rate. The surrogate model incorporates as much experimental detail as possible, but necessarily contains unmeasured parameters. The objective is to understand better the differences between model simulations and experimental results in complex microbial populations. The key methodology for simulation of chaotic dynamics, consistent with the measured dilution rate and microbial volume averages, was to cause the preference of p for r vs. c to vary with the r and c concentrations, to make r more competitive for nutrient than c, and to recycle some dying p biomass, leading to a modified version of the Monod kinetics model. Our mathematical model demonstrated that the occurrence of chaotic dynamics requires a predator, p, preference for r versus c to increase significantly with increases in r and c populations. Also included is a discussion of several generalizations of the existing model and a possible involvement of the minimum energy dissipation principle. This principle appears fundamental to thermodynamic systems including living systems. Several new experiments are suggested.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

Investigation of the process of mixing reagents in the laser chamber of active medium generators of continuous-wave chemical lasers with promising nozzle arrays

It is shown that the use of a toothed nozzle array in the active medium generator (AMG) of continuous-wave chemical lasers (CWCLs) allows more efficient mixing of reagents to be achieved in the laser chamber than in an AMG with a slit nozzle array. The results of calculations of parameters of the active medium of a DF CWCL operating in the gain regime for the AMG with traditional and promising nozzle array designs are presented. A kinetic model of processes in the DF CWCL active medium is proposed, which allows obtaining a good agreement with the results of measurements of the small-signal gain. (paper)

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

Nucleation rate controlled grain boundary and lattice creep

Nucleation versus diffusion rate-limited bicrystal and single crystal creep exhibit different scaling dependencies that enable the mechanisms to be isolated when measured as a function of sample size. It has recently been suggested that nucleation rate-limited kinetic models generally describe the non-Newtonian portion of the creep response well, but more direct evidence is required. This work analyzes the grain boundary creep response of UO 2 , a pyrochlore high entropy oxide, silver, and palladium, along with the lattice creep of silver using small-scale in situ loading in the transmission electron microscope. At small sizes, each system exhibits scale dependence associated with nucleation rate-limited kinetics. Fits of the data produce activation volumes on the order of a few Burgers vectors cubed with positive temperature coefficients as expected for nucleation kinetics. The activation enthalpies fall in the range of about 0.4 eV to 1.7 eV, being lower for the metals and higher for the oxides.

36 MATERIALS SCIENCE↗
DOE OSTI2018

Nonisothermal Crystallization Kinetics of Glassy Mold Fluxes

Crystallization of the solid glassy mold flux film occurring in the gap between the initial shell and mold wall is important, as it determines the in-mold heat transfer and mold lubrication during the process of continuous casting. In order to study the nonisothermal crystallization behavior of the glassy mold flux film in the continuous casting mold, the continuous heating transformation diagram, crystallization mechanism, and precipitate phases were investigated using the single hot thermocouple technique, kinetic models, a scanning electron microscope, and an energy-dispersive spectrometer (EDS). The results show that the initial crystallization temperature for CaO-SiO{sub 2} based flux A ranges from [1086 K to 1147 K (813 °C to 874 °C)], which is lower than the case of CaO-Al{sub 2}O{sub 3} based flux B ranging from [1205 K to 1245 K (932 °C to 972 °C)]. The crystallization kinetics for flux A are constant nucleation rate, two-dimensional growth, and control by diffusion. For flux B, they are constant nucleation rate, three-dimensional growth, and control by interface reaction. Besides, the EDS results indicate that the precipitate crystals in fluxes A and B are CaSiO{sub 3} and Ca{sub 2}AlSiO{sub 4}, respectively.

36 MATERIALS SCIENCE↗
NASA NTRS2006

Dynamic Response Testing in an Electrically Heated Reactor Test Facility

Non-nuclear testing can be a valuable tool in development of a space nuclear power or propulsion system. In a non-nuclear test bed, electric heaters are used to simulate the heat from nuclear fuel. Standard testing allows one to fully assess thermal, heat transfer, and stress related attributes of a given system, but fails to demonstrate the dynamic response that would be present in an integrated, fueled reactor system. The integration of thermal hydraulic hardware tests with simulated neutronic response provides a bridge between electrically heated testing and full nuclear testing. By implementing a neutronic response model to simulate the dynamic response that would be expected in a fueled reactor system, one can better understand system integration issues, characterize integrated system response times and response characteristics, and assess potential design improvements at a relatively small fiscal investment. Initial system dynamic response testing was demonstrated on the integrated SAFE-100a heat pipe cooled, electrically heated reactor and heat exchanger hardware, utilizing a one-group solution to the point kinetics equations to simulate the expected neutronic response of the system (Bragg-Sitton, 2005). The current paper applies the same testing methodology to a direct drive gas cooled reactor system, demonstrating the applicability of the testing methodology to any reactor type and demonstrating the variation in system response characteristics in different reactor concepts. In each testing application, core power transients were controlled by a point kinetics model with reactivity feedback based on core average temperature; the neutron generation time and the temperature feedback coefficient are provided as model inputs. Although both system designs utilize a fast spectrum reactor, the method of cooling the reactor differs significantly, leading to a variable system response that can be demonstrated and assessed in a non-nuclear test facility.

Bragg-Sitton, Shannon M.↗
DOE OSTI2017

Non-isothermal Crystallization Kinetics of Mold Fluxes for Casting High-Aluminum Steels

This paper investigates the crystallization behavior of CaO-SiO{sub 2}- and CaO-Al{sub 2}O{sub 3}-based mold fluxes for casting high-aluminum steels using single hot thermocouple technology, developed kinetic models, and scanning electron microscope. The results showed that the crystallization ability of the typical CaO-SiO{sub 2}-based Flux A (CaO/SiO{sub 2} 0.62, Al{sub 2}O{sub 3} 2 mass pct) is weaker than that of CaO-Al{sub 2}O{sub 3}-based Flux B (CaO/SiO{sub 2} 4.11, Al{sub 2}O{sub 3} 31.9 mass pct) because of its higher initial crystallization temperature. The crystallization kinetics of Flux A was “surface nucleation and growth, interface reaction control” in the overall non-isothermal crystallization process, whereas that of Flux B was “constant nucleation rate, 1-dimensional growth, diffusion control, in the primary crystallization stage, and then transformed into constant nucleation rate, 3-dimensional growth, interface reaction control in the secondary crystallization stage.” The energy dispersive spectroscopy results for Flux B suggested that the variations in the crystallization kinetics for Flux B are due to different crystals precipitating in the primary (BaCa{sub 2}Al{sub 8}O{sub 15}) and secondary (CaAl{sub 2}O{sub 4}) crystallization periods during the non-isothermal crystallization process.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Isopentane Disproportionation in Lewis Acidic Chloroaluminate Ionic Liquid

Lewis acidic chloroaluminate ionic liquid diluted in dichloromethane catalyzes the disproportionation of alkanes, such as isopentane, via carbenium ions in a reaction readily initiated by carbenium ion precursors such as tert-butyl chloride (TBC). The carbenium ion-AlCl 4 – ion-pairs stabilized by the ionic liquid-dichloromethane solution are the key intermediates in two distinctive kinetic regimes, i.e., a transient regime (0–5 min) and a steady-state regime (after 5 min). The transient regime constitutes the majority of isopentane conversion and is governed by the initial carbenium ion concentration. In the steady-state regime, disproportionation occurs at a considerably lower rate, affected by the carbenium ion concentration, the concentration of the ionic liquid, and the reaction temperature. The formation of alkenes observed in the 1 H NMR spectra of the reacting substrates, along with the DFT calculations, suggests that deprotonation of carbenium ion-pairs reduces their concentration, decreasing, in turn, the reaction rate. Kinetic modeling indicates that the transient regime is significantly controlled by the hydride transfer (k HT ) and the deprotonation rate constants (k DP ), while the steady-state regime is additionally influenced by the alkene protonation rate constant (k –DP ). The overall activation energy of the reaction at the steady state, expressed as E a,steady-state regime = E a,HT – E a,DP + E a,–DP , was 54 kJ/mol. The reaction mechanism and the kinetics highlight the potential of Lewis acid-catalyzed conversions of hydrocarbons under remarkably mild conditions.

Deprotonation↗
DOE OSTI2023

Hydrophobic pockets built in polymer micelles enhance the reactivity of Cu 2+ ions

We report the hydrophobicity-enhanced reactivity of Cu2+ ions as an ester hydrolase. Using a dipicolylamine (DPA) containing reversible addition–fragmentation chain transfer agent, the synthetic sequence, either hydrophobic or hydrophilic first in amphiphilic block copolymers of polystyrene-block-poly(N'N-dimethylacrylamide) (PS-b-PDMA), was varied to control the location of the binding motif, DPA, in the hydrophobic core or on the hydrated corona of polymer micelles. The hydrophobicity of Cu 2+ sites showed a significant impact (as large as 60 times more activity) on their catalytic efficiency towards ester hydrolase. With two different kinetic modes, including Michaelis–Menten and the reverse saturation kinetics models, the binding constant Kb of the substrates to Cu 2+ sites were quantitatively analyzed and we demonstrate that hydrophobicity favors the binding of the substrates to Cu2+ sites at polymer micelles with smaller sizes, however, K b decays exponentially with micellar diameters. Despite the diffusion barrier, hydrophobicity shows a profound impact on the catalytic rate constant kc that measures the single conversion rate of bound substrates to products. There is a 16–20 times kinetic enhancement in the hydrolase activity, completely endowed by the hydrophobic microenvironment of Cu 2+ sites compared to micelles with similar sizes. Our results indicate how the hydrophobicity of Cu 2+ -containing micelles can impact the catalytic efficiency and potentially illustrate a promising way toward the design of bioinspired catalysts.

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

Development of a Reduced-Order Model for Reacting Gas-Solids Flow using Proper Orthogonal Decomposition

This report summarizes the objectives, tasks and accomplishments made during the three year duration of this research project. The report presents the results obtained by applying advanced computational techniques to develop reduced-order models (ROMs) in the case of reacting multiphase flows based on high fidelity numerical simulation of gas-solids flow structures in risers and vertical columns obtained by the Multiphase Flow with Interphase eXchanges (MFIX) software. The research includes a numerical investigation of reacting and non-reacting gas-solids flow systems and computational analysis that will involve model development to accelerate the scale-up process for the design of fluidization systems by providing accurate solutions that match the full-scale models. The computational work contributes to the development of a methodology for obtaining ROMs that is applicable to the system of gas-solid flows. Finally, the validity of the developed ROMs is evaluated by comparing the results against those obtained using the MFIX code. Additionally, the robustness of existing POD-based ROMs for multiphase flows is improved by avoiding non-physical solutions of the gas void fraction and ensuring that the reduced kinetics models used for reactive flows in fluidized beds are thermodynamically consistent.

42 ENGINEERING↗
DOE OSTI2017

Experimental Study of Gas-Hydrate Dissociation Kinetics at Constant Temperatures Above 0°C

Recently, gas-hydrate storage technology has been drawing increasing attention as a new means of gas storage and transportation and may in the future become the optimal means of gas storage as compared with LNG technology because of its economic and safety advantages. The main purpose of the present work was to investigate experimentally the hydrate dissociation mechanism using a kinetic model at different temperatures and atmospheric pressure in the range 5-15°C. The obtained results have great significance for gas-hydrate storage technology and its application in industry.

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

The Reaction of HO2 and CH3O2: CH3OOH Formed from the Singlet Electronic State Surface

High-level coupled-cluster calculations in combination with two-dimensional master equation simulations were used to study the HO2 + CH3O2 reaction, which plays an important role in the oxidation of methane and hydrocarbons in the Earth’s atmosphere and low-temperature combustion. The main reaction pathways taking place on the lowest-lying triplet and singlet potential energy surfaces (PES) were characterized. Interestingly, methyl hydroperoxide (CH3OOH), the sole product, could be produced from both the triplet and singlet PESs, with a ratio of roughly 9:1. Formaldehyde is not made as a primary product, but can be formed via secondary chemistry. The formation of methyl tetraoxide (MTO) from the singlet PES is unimportant. The calculated reaction rate coefficients were found to be practically pressure-independent for p ≤ 760 Torr and can be given by k(T)=2.75×10−13×e+1.75 kcal mol−1/RT (in cm3/s), an expression useful for kinetics modeling over the range T = 200–800 K. The rate constant has a slight negative Arrhenius energy dependence of about −1.75 kcal mol–1, falling about a factor of 30 from 200 K to 800 K.

2DME↗
DOE OSTI2018

Shock interactions with heterogeneous energetic materials

The complex physical phenomenon of shock wave interaction with material heterogeneities has significant importance and nevertheless remains little understood. In many materials, the observed macroscale response to shock loading is governed by characteristics of the microstructure. Yet the majority of computational studies aimed at predicting phenomena affected by these processes, such as initiation and propagation of detonation waves in explosives, or shock propagation in geological materials, employ continuum material and reactive burn model treatment. In an effort to highlight the grain-scale processes that underlie the observable effects in an energetic system, a grain-scale model for hexanitrostilbene (HNS) has been developed. Measured microstructures were used to produce synthetic computational representations of the pore structure, and a density functional theory molecular dynamics (DFT-MD) derived equation of state (EOS) was used for the fully dense HNS matrix. The explicit inclusion of microstructure along with a fully-dense EOS resulted in close agreement with historical shock compression experiments. More recent experiments on dynamic reaction threshold were also reproduced by inclusion of a global kinetics model. The complete model was shown to reproduce accurately the expected response of this heterogeneous material to shock loading. Mesoscale simulations were shown to provide clear insight into the nature of threshold behavior, and are a way to understand complex physical phenomena.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Verification of an energy-conserving semi-implicit electrostatic particle-in-cell scheme for modeling high-density plasma at scale

A verification study of a semi-implicit energy-conserving electrostatic particle-in-cell algorithm is presented. The algorithm relaxes the time-step and mesh-size constraints that require resolution of the plasma period and Debye length associated with traditional explicit momentum-conserving particle-in-cell algorithms. Physical implications and applicability of using the semi-implicit scheme for modeling high-density plasmas are discussed. Where possible, numerical results are compared against analytical solutions. The simulation results indicate that the algorithm is stable at time steps larger than twice the inverse plasma frequency and cell sizes larger than the Debye length. It is found that the algorithm gives adequate results, provided that the distribution function and the spatiotemporal scales dictating the physics of the problem are resolved. As such, the algorithm may provide a robust method for kinetic modeling of high-density plasmas at scale.

Cyclotron resonance↗
DOE OSTI2023

Experiments on excitation of Alfvén eigenmodes by alpha-particles with bump-on-tail distribution in JET DTE2 plasmas

Dedicated experiments were performed in JET DTE2 plasmas for obtaining an α-particle bump-on-tail (BOT) distribution aiming at exciting Alfvén eigenmodes (AEs). Neutral beam injection-only heating with modulated power was used so that fusion-born α-particles were the only ions present in the MeV energy range in these DT plasmas. The beam power modulation on a time scale shorter than the α-particle slowing down time was chosen for modulating the α-particle source and thus sustaining a BOT in the α-particle distribution. High-frequency modes in the toroidicity-induced Alfven eigenmode (TAE) frequency range and multiple short-lived modes in a wider frequency range have been detected in these DT discharges with interferometry, soft x-ray cameras, and reflectometry. The modes observed were localised close to the magnetic axis and were not seen in the Mirnov coils. Analysis with the TRANSP and Fokker-Planck FIDIT codes confirms that α-particle distributions with BOT in energy were achieved during some time intervals in these discharges though no clear correlation was found between the times of the high-frequency mode excitation and the BOT time intervals. The combined magneto-hydrodynamic (MHD) and kinetic modelling studies show that the high-frequency mode in the TAE frequency range is best fitted with a TAE of toroidal mode number n= 9. This mode is driven mostly by the on-axis beam ions while the smaller drive due to the pressure gradient of α-particles allows overcoming the marginal stability and exciting the mode. The observed multiple short-lived modes in a wider frequency range are identified as the on-axis kinetic AEs predicted in Rosenbluth and Rutherford.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Reduction kinetics of hercynite redox materials for solar thermochemical water splitting

Solar thermochemical water splitting (STWS) using a hercynite redox cycle is a promising technology for producing renewable H 2 . Here, the reduction kinetics of hercynite (FeAl 2 O 4 ) is evaluated using thermogravimetric and XRD analyses. The results indicate that as-prepared hercynite materials undergo reduction via two different reaction mechanisms. The reaction first proceeds by a nucleation and growth mechanism (AE1), followed by a third-order kinetic model (F3). XRD analyses show the occurrence of superstoichiometric oxygen in the spinel structure of FeAl 2 O 4+δ in the second reaction mechanism, which indicates the formation of cationic vacancies. TGA and XRD analyses confirm that hercynite maintains its spinel structure when the materials are thermally reduced and oxidized with steam. Differential scanning calorimetry (DSC) performed on FeAl 2 O 4 particles shows two heats of reaction in agreement with the two observed mechanisms. This is the first kinetics study to report on the thermal reduction of hercynite.

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

An experimental investigation of flame and autoignition behavior of propane

We report autoignition delay time data are one important means to develop, quantify, and validate fundamental understanding of combustion chemistry at low temperatures (T<1200 K). However, low-temperature chemistry often has higher uncertainties and scatter in the experimental data compared with high-temperature ignition data (T>1200 K). In this study, autoignition properties of propane and oxygen mixtures were investigated using the University of Michigan rapid compression facility in order to understand the effects of ignition regimes on low-temperature ignition data. For the first time for propane, autoignition delay times were determined from pressure histories, and autoignition characteristics were simultaneously recorded using highspeed imaging of the test section through a transparent end-wall. Propane mixtures with fuel-to-O 2 equivalence ratios of φ = 0.25 and φ = 0.5 and O 2 -to-inert gas molar ratios of 1:3.76 were studied over the pressure range of 8.9 to 11.3 atm and the temperature range of 930 – 1070 K. The results showed homogeneous or strong autoignition occurred for all φ = 0.25 experiments, and inhomogeneous or mixed autoignition occurred for all φ = 0.5 experiments. While a limited temperature range is covered in the study, importantly the data span predicted transitions in autoignition behavior, allowing validation of autoignition regime hypotheses. Specifically, the results agree well with strong-autoignition limits proposed based on the Sankaran Criterion. The autoignition delay time data at the strong-ignition conditions are in excellent agreement with predictions using a well-validated detailed reaction mechanism from the literature and a zerodimensional modeling assumption. However, the experimental data at the mixed autoignition conditions were systematically faster than the model predictions, particularly at lower temperatures (T< ~970 K). The results are an important addition to the growing body of data in the literature that show mixed autoignition phenomena are important sources of the higher scatter observed in the low-temperature autoignition data for propane and other fuels. The results are discussed in terms of different methods to capture the effects of pre-autoignition heat release associated with mixed autoignition conditions and thereby address some of the discrepancies between kinetic modeling and experimental measurements.

42 ENGINEERING↗