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,009 records · Page 56

Kinetic approach to defect reduction in directed self-assembly

As a potential solution to next-generation nanolithography, directed self-assembly (DSA) of block copolymers (BCPs) is still restrained in high-volume manufacturing primarily due to its defectivity issue. Though defects possess greater free energies than aligned morphologies and are highly energetically unfavorable, they can be kinetically trapped by the energy barriers and persist for a long time during annealing. Therefore, understanding the kinetics of defect annihilation is crucial in revealing the mechanism of defect formation and in further reducing defectivity in DSA. We focus on two types of predominant defects in DSA- dislocation and bridge. Furthermore, a kinetic model of each defect type is developed through statistical analysis of experimental data, providing insight into possible approaches of further defect reduction. We also investigate the impact of annealing temperature and film thickness on annihilation kinetics and discuss the reasons behind the observed results. By simply optimizing annealing conditions and film thickness, we have successfully reduced the total defect density by 1 order of magnitude. Though these findings are based on polystyrene-b-poly(methyl methacrylate) (PS-b-PMMA), we anticipate they could be readily applied to other BCP platforms as well.

42 ENGINEERING↗

Methyl propyl ether combustion

Ethers are potential diesel blendstocks that have high cetane numbers and identified cost-effective production pathways. Among ethers, methyl propyl ether (MPE) is a small enough molecule to build a good detailed model for and a large enough molecule for studying key intramolecular hydrogen migrations involved in low-temperature feedback loops for larger ethers. In this study, a model was generated for the combustion of MPE from 600-900 K and 10-20 bar using the Reaction Mechanism Generator (RMG) package. Results from rapid compression machine (RCM) experiments were used to validate the kinetic model. At low temperatures, the generated model successfully captured the dominant low-temperature feedback loops. At higher temperatures, propanal is a major combustion intermediate for PME, suggesting that aldehydes likely play a significant role in larger ethers' combustion.

09 BIOMASS FUELS↗

Photoionization of seeded combustion products as a method of enhancing the efficiency of magnetohydrodynamic power generators

Here, in this study, we performed an experimental and computational investigation into the feasibility of utilizing photoionization to enhance the electrical conductivity of seeded oxy-fuel combustion products and improve the performance of magnetohydrodynamic (MHD) power generators. We applied a variety of optical and microwave diagnostics to study the ionization and recombination processes of potassium excited by an excimer laser in a high-velocity oxy-fuel free jet. Computational fluid dynamic (CFD) simulations were performed to model the thermophysical properties and species densities of the free jet. The CFD results were validated with position-dependent potassium concentration measurements. Electron recombination exponential lifetimes were measured through time-resolved microwave transmission. The experimental electron lifetimes were compared with lifetimes calculated from CFD-predicted species densities and literature recombination rates. It was determined that K + or O 2 are the most likely recombination partners for photoionized electrons. Time-resolved fluorescence measurements provided evidence of an ionization pathway involving a two-photon ionization of KOH . Finally, a zero-dimensional chemical kinetic model was developed to assess the fundamental viability of inducing a non-equilibrium electron population to provide a net energy return in combustion-driven MHD power generators. We determined that a high energy return is feasible for targeting electrode boundary layers with ultraviolet photoionization. We also found that photoionization could potentially lower the required temperature of the bulk gas flow.

20 FOSSIL-FUELED POWER PLANTS↗

Stochastic methods for aerosol chemistry: a compact molecular description of functionalization and fragmentation in the heterogeneous oxidation of squalane aerosol by OH radicals

The heterogeneous oxidation of organic aerosol by hydroxyl radicals (OH) can proceed through two general pathways: functionalization, in which oxygen functional groups are added to the carbon skeleton, and fragmentation, in which carbon–carbon bonds are broken, producing higher volatility, lower molecular weight products. An ongoing challenge is to develop a quantitative molecular description of these pathways that connects the oxidative evolution of the average aerosol properties (e.g. size and hygroscopicity) to the transformation of free radical intermediates. In order to investigate the underlying molecular mechanism of aerosol oxidation, a relatively compact kinetics model is developed for the heterogeneous oxidation of squalane particles by OH using free radical intermediates that convert reactive hydrogen sites into oxygen functional groups. Stochastic simulation techniques are used to compare calculated system properties over ten oxidation lifetimes with the same properties measured in experiment. The time-dependent average squalane aerosol mass, volume, density, carbon number distribution of scission products, and the average elemental composition are predicted using known rate coefficients. For functionalization, the calculations reveal that the distribution of alcohol and carbonyl groups is controlled primarily by the initial OH abstraction rate and to lesser extent by the branching ratio between secondary peroxy radical product channels. For fragmentation, the calculations reveal that the formation of activated alkoxy radicals with neighboring functional groups controls the molecular decomposition, particularly at high O/C ratios. Finally, this kinetic scheme provides a framework for understanding the oxidation chemistry of a model organic aerosol and informs parameterizations of more complex systems.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Numerical estimation of phase transformations in solid state during Yb:YAG laser heating of steel sheets

This work concerns the numerical modeling of heat transfer and phase transformations in solid state occurring during the Yb:YAG laser beam heating process. The temperature field is obtained by the numerical solution into transient heat transfer equation with convective term. The laser beam heat source model is developed using the Kriging interpolation method with experimental measurements of Yb:YAG laser beam profile taken into account. Phase transformations are calculated on the basis of Johnson - Mehl - Avrami (JMA) and Koistinen - Marburger (KM) kinetics models as well as continuous heating transformation (CHT) and continuous cooling transformation (CCT) diagrams for S355 steel. On the basis of developed numerical algorithms 3D computer simulations are performed in order to predict temperature history and phase transformations in Yb:YAG laser heating process.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

A Numerical Model for the Dynamics of Pickup Ions Outside the Heliopause and IBEX “Ribbon” Observation

The leading mechanism for the origin of the Interstellar Boundary Explorer (IBEX) “ribbon” of enhanced energetic neutral atoms (ENAs) from the outer heliosphere is the so-called secondary ENA process. In this mechanism, neutralized ions escape the heliosphere, and after several times of charge exchange collision,s some propagate back toward Earth in directions nearly perpendicular to the local interstellar magnetic field. However, the physical processes governing the distribution of the pickup ions (PUIs) outside the heliopause are still unclear. In this study, we build a new global model of the IBEX ribbon, where a key component is to calculate the dynamics of the PUIs outside the heliopause by solving the gyrophase-averaged focused transport equations on top of the background magnetohydrodynamic-kinetic model for the heliosphere-interstellar medium interaction. We discuss how the properties of the simulated ribbon change with different scattering parameters and show simulation results of the ENA sky maps, ribbon centers, ribbon profiles at several azimuthal slices, the ribbon ENA source region, and the velocity distribution of PUIs in the ENA source region for different scattering frequencies. Our results show that the model can provide reasonable ribbon sky maps comparable to the observed ribbon when the PUI mean free path is on the order of several thousand astronomical units.

79 ASTRONOMY AND ASTROPHYSICS↗

On analysis and stochastic modeling of the particle kinetic energy equation in particle-laden isotropic turbulent flows

In this study, we analyze three-dimensional particle-laden, isotropic turbulence to develop an understanding of inertial particle dynamics from a kinetic energy perspective. Data trends implying inhomogeneous sampling of the flow by particles are identified and used to support a proposed particle behavior: particles appear to accumulate in regions of low flow kinetic energy over time because they lose kinetic energy and slow down in such regions, ultimately causing them to spend more time there. To elucidate this behavior, we derive a particle kinetic energy equation from the particle momentum equation, which incorporates inertial effects through the Schiller–Naumann drag correlation. Upon extracting fundamental physics from this equation, hypotheses regarding the role of the Stokes number in the temporal change of particle kinetic energy and the previously proposed particle behavior are evaluated using simulation data considering three Stokes numbers. Finally, a Fokker–Planck equation is used to derive the steady-state probability density function of the particle kinetic energy. The model fits the simulation data well and provides a tool for further investigation into understanding preferential concentration, as well as a reduced order model for predicting particle kinetic energy in turbulent flows.

42 ENGINEERING↗

Neutral interstellar helium parameters based on Ulysses/GAS and IBEX-LO observations: What are the reasons for the differences?

Recent analysis of the interstellar helium fluxes measured in 2009-2010 at Earth's orbit by the Interstellar Boundary Explorer (IBEX) has suggested that the interstellar velocity (both direction and magnitude) is inconsistent with that derived previously from Ulysses/GAS observations made in the period from 1990 to 2002 at 1.5-5.5 AU from the Sun. Both results are model dependent, and models that were used in the analyses are different. In this paper, we perform an analysis of the Ulysses/GAS and IBEX-Lo data using our state-of-the-art three-dimensional time-dependent kinetic model of interstellar atoms in the heliosphere. For the first time, we analyze Ulysses/GAS data from year 2007, the closest available Ulysses/GAS observations in time to the IBEX observations. We show that the interstellar velocity derived from the Ulysses 2007 data is consistent with previous Ulysses results and does not agree with the velocity derived from IBEX. This conclusion is very robust since, as is shown in the paper, it does not depend on the ionization rates adopted in theoretical models. We conclude that Ulysses data are not consistent with the new local interstellar medium (LISM) velocity vector from IBEX. In contrast, IBEX data, in principle, could be explained with the LISM velocity vector derived from the Ulysses data. This is possible for the models where the interstellar temperature increased from 6300 K to 9000 K. There is a need to perform further studies of possible reasons for the broadening of the helium signal core measured by IBEX, which could be an instrumental effect or could be due to unconsidered physical processes.

79 ASTRONOMY AND ASTROPHYSICS↗

Consistent Chemical Mechanism from Collaborative Data Processing

Numerical tool of Process Informatics Model (PrIMe) is mathematically rigorous and numerically efficient approach for analysis and optimization of chemical systems. It handles heterogeneous data and is scalable to a large number of parameters. The Boundto-Bound Data Collaboration module of the automated data-centric infrastructure of PrIMe was used for the systematic uncertainty and data consistency analyses of the H 2 /CO reaction model (73/17) and 94 experimental targets (ignition delay times). The empirical rule for evaluation of the shock tube experimental data is proposed. The initial results demonstrate clear benefits of the PrIMe methods for an evaluation of the kinetic data quality and data consistency and for developing predictive kinetic models.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Development of Kinetics for Soot Oxidation at High Pressures Under Fuel-Lean Conditions

The focus of the proposed research was to develop kinetic models for soot oxidation with the hope of developing a validated, predictive, multi-­scale, combustion model to optimize the design and operation of evolving fuels in advanced engines for transportation applications. The work focused on the relatively unstudied area of the fundamental mechanism for soot oxidation. The objectives include understanding of the kinetics of soot oxidation by O2 under high pressure which require: 1) development of intrinsic kinetics for the surface oxidation, which takes into account the dependence of reactivity upon nanostructure and 2) evolution of nanostructure and its impact upon oxidation rate and 3) inclusion of internal surface area development and possible fragmentation resulting from pore development and /or surface oxidation. These objectives were explored for a variety of pure fuel components and surrogate fuels. This project was a joint effort between the University of Utah (UU) and Pennsylvania State University (Penn State). The work at the UU focuses on experimental studies using a two-­stage burner and a high- pressure thermogravimetric analyzer (TGA). Penn State provided HRTEM images and guidance in the fringe analysis algorithms and parameter quantification for the images. This report focuses on completion done under supplemental funding.

02 PETROLEUM↗

Modeling a short dc discharge with thermionic cathode and auxiliary anode

A short dc discharge with a thermionic cathode can be used as a current and voltage stabilizer, but is subject to current oscillation. If instead of one anode two anodes are used, the current oscillations can be reduced. We have developed a kinetic model of such a discharge with two anodes, where the primary anode has a small opening for passing a fraction of the discharge current to an auxiliary anode. The model demonstrates that the current-voltage relationship of the discharge with two anodes is characterized everywhere by positive slope, i.e., positive differential resistance. Therefore, the discharge with two anodes is expected to be stable to the spontaneous oscillation in current that is induced by negative differential resistance. As a result, such a discharge can be used in an engineering application that requires stable plasma, such as a current and voltage stabilizer.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Analysis and prediction of reaction kinetics using the degree of rate control

“Degree of rate control” (DRC) analysis provides a quantitative approach for analysing the kinetics of multi-step reaction mechanisms that has been widely applied to both heterogeneous and homogeneous catalysis research, as well as electrocatalysis. The DRC of any given transition state or intermediate is defined as a partial derivative such that it approximately equals the fractional increase in net rate to the product of interest per differential decrease in its standard-state free energy for that species (÷RT), holding constant the standard-state free energies of all other transition states and intermediates. Even very complex mechanisms usually have only a few species with non-zero DRCs and are thus the species whose interactions with the catalyst most strongly affect the net rate. These key DRC values thus offer a simple and intuitive route to optimize catalyst materials, especially with the assistance of computational methods like density functional theory (DFT). These high-DRC species are also the species whose energetics must be most accurately measured or calculated to achieve an accurate kinetic model for any reaction mechanism. In simple cases with a single “rate-determining step”, the DRC for its transition state (TS) is + 1. Catalyst-bound intermediates, on the other hand, often have negative DRCs equal to a small integer times their fractional occupancy of catalyst sites. The apparent activation energy equals a weighted average of the standard-state enthalpies (relative to reactants) of all the species (intermediates, transition states and products) in the reaction mechanism, each weighted by its DRC (+RT). It has been shown that the apparent transfer coefficient in electrocatalysis, an inverted form of the Tafel slope, is a weighted average of the number of electrons transferred to generate each intermediate or product species in the mechanism, weighted again by the DRC. Quantitative analysis of kinetic isotope effects (KIEs, or the ratio of net rates for different reactant isotopes) in complex mechanisms has shown that the logarithm of the KIE equals the weighted average over all species in the mechanism of the difference between the two isotopes in their standard-state free energies (÷RT), again weighted by the DRC. The reaction orders with respect to fluid-phase concentrations of reactants, products and intermediates have also been proven to be directly related to DRCs. Thus, there are numerous experimental observables which equate to short linear combinations of DRCs, so that combinations of experimental measurements might provide access to DRC values. Since its invention for transition states in 1994 and its generalization to include intermediates in 2009, the DRC has thus far mainly been calculated for microkinetic models based either entirely on DFT or on DFT with the key energies (i.e., those for high-DRC species) being fine-tuned to match experiments. The relationships summarized in this work provide new opportunities for using experiments earlier in the development and optimization of microkinetic models that require input from computational catalysis.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Experimental Investigation of Uranium Volatility during Vapor Condensation

The predictive models that describe the fate and transport of radioactive materials in the atmosphere following a nuclear incident (explosion or reactor accident) assume that uranium-bearing particulates would attain chemical equilibrium during vapor condensation. Here, we show that kinetically driven processes in a system of rapidly decreasing temperature can result in substantial deviations from chemical equilibrium. This can cause uranium to condense out in oxidation states (e.g., UO 3 vs UO 2 ) that have different vapor pressures, significantly affecting uranium transport. To demonstrate this, we synthesized uranium oxide nanoparticles using a flow reactor under controlled conditions of temperature, pressure, and oxygen concentration. The atomized chemical reactants passing through an inductively coupled plasma cool from ~5000 to 1000 K within milliseconds and form nanoparticles inside a flow reactor. The ex situ analysis of particulates by transmission electron microscopy revealed 2–10 nm crystallites of fcc-UO 2 or α-UO 3 depending on the amount of oxygen in the system. α-UO 3 is the least thermodynamically preferred polymorph of UO 3 . The absence of stable uranium oxides with intermediate stoichiometries (e.g., U 3 O 8 ) and sensitivity of the uranium oxidation states to local redox conditions highlight the importance of in situ measurements at high temperatures. Therefore, we developed a laser-based diagnostic to detect uranium oxide particles as they are formed inside the flow reactor. Our in situ measurements allowed us to quantify the changes in the number densities of the uranium oxide nanoparticles (e.g., UO 3 ) as a function of oxygen gas concentration. Our results indicate that uranium can prefer to be in metastable crystal forms (i.e., α-UO 3 ) that have higher vapor pressures than the refractory form (i.e., UO 2 ) depending on the oxygen abundance in the surrounding environment. This demonstrates that the equilibrium processes may not dominate during rapid condensation processes, and thus kinetic models are required to fully describe uranium transport subsequent to nuclear incidents.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Effect of fulvic acid surface coatings on plutonium sorption and desorption kinetics on goethite

The rates and extent of plutonium (Pu) sorption and desorption onto mineral surfaces are important parameters for predicting Pu mobility in subsurface environments. The presence of natural organic matter, such as fulvic acid (FA), may influence these parameters. We investigated the effects of FA on Pu(IV) sorption/desorption onto goethite in two scenarios: when FA was (1) initially present in solution or (2) found as organic coatings on the mineral surface. A low pH was used to maximize FA coatings on goethite. Experiments were combined with kinetic modeling and speciation calculations to interpret variations in Pu sorption rates in the presence of FA. Our results indicate that FA can change the rates and extent of Pu sorption onto goethite at pH 4. Differences in the kinetics of Pu sorption were observed as a function of the concentration and initial form of FA. The fraction of desorbed Pu decreased in the presence of FA, indicating that organic matter can stabilize sorbed Pu on goethite. These results suggest that ternary Pu–FA–mineral complexes could enhance colloid-facilitated Pu transport. In conclusion, more representative natural conditions need to be investigated to quantify the relevance of these findings.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗

Experiments to populate and validate a processing model for polyurethane foam. BKC 44306 PMDI-10

We are developing computational models to elucidate the expansion and dynamic filling process of a polyurethane foam, PMDI. The polyurethane of interest is chemically blown, where carbon dioxide is produced via the reaction of water, the blowing agent, and isocyanate. The isocyanate also reacts with polyol in a competing reaction, which produces the polymer. Here we detail the experiments needed to populate a processing model and provide parameters for the model based on these experiments. The model entails solving the conservation equations, including the equations of motion, an energy balance, and two rate equations for the polymerization and foaming reactions, following a simplified mathematical formalism that decouples these two reactions. Parameters for the polymerization kinetics model are reported based on infrared spectrophotometry. Parameters describing the gas generating reaction are reported based on measurements of volume, temperature and pressure evolution with time. A foam rheology model is proposed and parameters determined through steady-shear and oscillatory tests. Heat of reaction and heat capacity are determined through differential scanning calorimetry. Thermal conductivity of the foam as a function of density is measured using a transient method based on the theory of the transient plane source technique. Finally, density variations of the resulting solid foam in several simple geometries are directly measured by sectioning and sampling mass, as well as through x-ray computed tomography. These density measurements will be useful for model validation once the complete model is implemented in an engineering code.

97 MATHEMATICS AND COMPUTING↗

Cloud based tool for analysis of chemical kinetic mechanisms

The size and complexity of chemical kinetic models used to simulate combustion of hydrocarbon fuels continue to increase with our improved understanding of the underlying physical mechanisms, along with the improved ability to utilize such models as the result of increased computational power and efficiency. As mechanisms grow beyond thousands of species and tens of thousands of reactions, it becomes increasingly difficult to manually check for errors and physical inconsistencies. In this work, we present several automated methods to check for such issues in the specification of chemical reaction models. First, we demonstrate how discontinuities in thermodynamic data can cause simulation difficulties manifested as long wall-clock times and convergence failures. To correct this type of problem, we describe an automated method for refitting thermodynamic parameters. We also outline several methods to check the timescales of reaction rate coefficients to ensure physical consistency. All the methods are made available through a web-based tool (https://combustiontools.llnl.gov) intended to aid mechanism developers and users to improve the accuracy and performance of fuel models used by the combustion community.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Self-consistent quantum kinetic theory of diatomic molecule formation

A quantum kinetic theory of molecule formation is presented which includes three-body recombination and radiative association for a thermodynamically closed system which may or may not exchange energy with its surrounding at a constant temperature. The theory uses a Sturmian representation of a two-body continuum to achieve a steady-state solution of a governing master equation which is self-consistent in the sense that detailed balance between all bound and unbound states is rigorously enforced. The role of quasibound states in catalyzing the molecule formation is analyzed in complete detail. The theory is used to make three predictions which differ from conventional kinetic models. These predictions suggest significant modifications may be needed to phenomenological rate constants which are currently in wide use. Implications for models of low and high density systems are discussed.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Effects of local order parameter dependent transport coefficient in diblock copolymers under applied electric fields

In this work, we present an approach for constructing thermodynamically consistent time-dependent models relevant to thin films of diblock copolymers in applied electric fields. The approach is based on the principles of linear irreversible thermodynamics, and, in this work, it is applied to study the effects of electric fields on thin films of incompressible diblock copolymers. Enforcement of local incompressibility constraint at all times leads to a local order parameter dependent transport coefficient in the model for the diblock copolymers. The dependence of the transport coefficient on the local order parameter is used to relate it with the diffusion constant of Rouse chains and leads to sensitivity of the model to initial conditions. In addition, transient behavior is found to be affected when compared with an ad hoc model assuming a constant transport coefficient. Numerical results such as electric field induced alignment of lamellae domains due to the field are found to be in qualitative agreement with experiments. This approach opens up a systematic way of developing kinetic models for simulating effects of electrolytes added to thin films containing diblock copolymers in the presence of applied electric fields.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗