Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Kinetic modeling”

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

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

At least 1,225 records · Page 68

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↗

Collisionless electrostatic interchange instabilities

The linear Vlasov dispersion equation for electrostatic plasma instabilities driven by gravity and weak density gradients perpendicular to a uniform magnetic field is derived and solved numerically. Two interchange instabilities emerge: the well-known fluid mode at long wavelengths and a kinetic model at wavelengths short compared with the ion gyroradius. The properties of both instabilities are studied, as well as the effects of gravity on the universal and lower-hybrid density drift instabilities. The results show that the kinetic interchange generally has a larger growth rate than the fluid interchange instability, indicating that, whenever the latter is present in a collisionless plasma, the former may also be found.

Gray, S. P.↗

CFD Analysis of Emissions for a Candidate N+3 Combustor

An effort was undertaken to analyze the performance of a model Lean-Direct Injection (LDI) combustor designed to meet emissions and performance goals for NASA's N+3 program. Computational predictions of Emissions Index (EINOx) and combustor exit temperature were obtained for operation at typical power conditions expected of a small-core, high pressure-ratio (greater than 50), high T3 inlet temperature (greater than 950K) N+3 combustor. Reacting-flow computations were performed with the National Combustion Code (NCC) for a model N+3 LDI combustor, which consisted of a nine-element LDI flame-tube derived from a previous generation (N+2) thirteen-element LDI design. A consistent approach to mesh-optimization, spraymodeling and kinetics-modeling was used, in order to leverage the lessons learned from previous N+2 flame-tube analysis with the NCC. The NCC predictions for the current, non-optimized N+3 combustor operating indicated a 74% increase in NOx emissions as compared to that of the emissions-optimized, parent N+2 LDI combustor.

nitrogen oxides↗

Analysis of Nitrogen Shockwave Radiation at Velocity Around 7km/s

In the last decade, the Electric Arc Shock Tube (EAST) facility has been widely used to provide experimental data aiming to validate radiative heating models used for atmospheric entry, mostly at velocities above 9km/s in air. In more recent years, another experimental campaign aimed to look at velocity below 9 km/s and its analysis resulted in a significant improvement of the air kinetic model using new rates found in the literature but discrepancies were still observed, especially for the N2 second positive system. Despite these previous studies, the overall uncertainties related to surface heating for lunar return remain around 50% To reduce these uncertainties, the last experimental campaign at EAST, named test 62, was conducted in pure nitrogen at velocity ranging from 6 to 11 km/s at a freestream pressure of 0.2 Torr to provide more fundamental data on nitrogen. A previous study analyzed the data obtained at a velocity of 10.3 km/s, where atomic radiation is dominant, and provided spatially resolved density measurements of different states of nitrogen. This study will focus on experimental data obtained at shock velocity of 6.88 and 7.15 km/s where molecular nitrogen is expected to be the dominant source of radiation.

Augustin Tibere-Inglesse↗

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↗

A Stochastic Reduced-Order Model for Statistical Microstructure Descriptors Evolution

Integrated computational materials engineering (ICME) models have been a crucial building block for modern materials development, relieving heavy reliance on experiments and significantly accelerating the materials design process. However, ICME models are also computationally expensive, particularly with respect to time integration for dynamics, which hinders the ability to study statistical ensembles and thermodynamic properties of large systems for long time scales. To alleviate the computational bottleneck, we propose to model the evolution of statistical microstructure descriptors as a continuous-time stochastic process using a non-linear Langevin equation, where the probability density function (PDF) of the statistical microstructure descriptors, which are also the quantities of interests (QoIs), is modeled by the Fokker–Planck equation. In this work, we discuss how to calibrate the drift and diffusion terms of the Fokker–Planck equation from the theoretical and computational perspectives. The calibrated Fokker–Planck equation can be used as a stochastic reduced-order model to simulate the microstructure evolution of statistical microstructure descriptors PDF. Considering statistical microstructure descriptors in the microstructure evolution as QoIs, we demonstrate our proposed methodology in three integrated computational materials engineering (ICME) models: kinetic Monte Carlo, phase field, and molecular dynamics simulations.

97 MATHEMATICS AND COMPUTING↗

Introducing the VRT gas turbine combustor

An innovative annular combustor configuration is being developed for aircraft and other gas turbine engines. This design has the potential of permitting higher turbine inlet temperatures by reducing the pattern factor and providing a major reduction in NO(x) emission. The design concept is based on a Variable Residence Time (VRT) technique which allows large fuel particles adequate time to completely burn in the circumferentially mixed primary zone. High durability of the combustor is achieved by dual function use of the incoming air. The feasibility of the concept was demonstrated by water analogue tests and 3-D computer modeling. The computer model predicted a 50 percent reduction in pattern factor when compared to a state of the art conventional combustor. The VRT combustor uses only half the number of fuel nozzles of the conventional configuration. The results of the chemical kinetics model require further investigation, as the NO(x) predictions did not correlate with the available experimental and analytical data base.

Melconian, Jerry O.↗

Introducing the VRT gas turbine combustor

An innovative annular combustor configuration is being developed for aircraft and other gas turbine engines. This design has the potential of permitting higher turbine inlet temperatures by reducing the pattern factor and providing a major reduction in NO(x) emission. The design concept is based on a Variable Residence Time (VRT) technique which allows large fuel particles adequate time to completely burn in the circumferentially mixed primary zone. High durability of the combustor is achieved by dual function use of the incoming air. The feasibility of the concept was demonstrated by water analogue tests and 3-D computer modeling. The computer model predicted a 50 percent reduction in pattern factor when compared to a state of the art conventional combustor. The VRT combustor uses only half the number of fuel nozzles of the conventional configuration. The results of the chemical kinetics model require further investigation, as the NO(x) predictions did not correlate with the available experimental and analytical data base.

Melconian, Jerry O.↗

Experimental and modeling study of the pyrolysis and oxidation of an iso-paraffinic alcohol-to-jet fuel

Single pulse shock tube experiments were performed to determine the pyrolytic and oxidative decomposition products of an alcohol-to-jet fuel (ATJ) composed primarily of about 90% iso-dodecane. Experiments were conducted at 4 bar nominal pressure, 2 ms nominal reaction time, and temperatures ranging from 900 K to 1550 K. For oxidation, the equivalence ratio was 0.27. Gas chromatography (GC) was used to analyze the products of the chemical reactions and to determine their mole fractions. The experimental results were compared against mole fractions computed using a detailed iso-alkane kinetic model, developed for the present study. The model was able to accurately reproduce the experimental results of both ATJ pyrolysis and oxidation. The experimental data and the kinetic modelintroduced in the current study shall aid in progress in understanding the oxidation of the alternative jet fuels.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Understanding Automotive Exhaust Catalysts Using a Surface Science Approach: Model NOx Storage Materials

The structure-reactivity relationships of model BaO-based NOx storage/reduction catalysts were investigated under well controlled experimental conditions using surface science analysis techniques. The reactivity of BaO toward NO2, CO2, and H2O was studied as a function of BaO layer thickness [0\hBaO\30 monolayer (ML)], sample temperature, reactant partial pressure, and the nature of the substrate the NOx storage material was deposited onto. Most of the efforts focused on understanding the mechanism of NO2 storage either on pure BaO, or on BaO exposed to CO2 or H2O prior to NO2 exposure. The interaction of NO2 with a pure BaO film results in the initial formation of nitrite/nitrate ion pairs by a cooperative adsorption mechanism predicted by prior theoretical calculations. The nitrites are then further oxidized to nitrates to produce a fully nitrated surface. The mechanism of NO2 uptake on thin BaO films (\4 ML), BaO clusters (\1 ML) and mixed BaO/Al2O3 layers are fundamentally different: in these systems initially nitrites are formed only, and then converted to nitrates at longer NO2 exposure times. These results clarify the contradicting mechanisms presented in prior studies in the literature. After the formation of a nitrate layer the further conversion of the underlying BaO is slow, and strongly depends on both the sample temperature and the NO2 partial pressure. At 300 K sample temperature amorphous Ba(NO3)2 forms that then can be converted to crystalline nitrates at elevated temperatures. The reaction between BaO and H2O is facile, a series of Ba(OH)2 phases form under the temperature and H2O partial pressure regimes studied. Both amorphous and crystalline Ba(OH)2 phases react with NO2, and initially form nitrites only that can be converted to nitrates. The NO2 adsorption capacities of BaO and Ba(OH)2 are identical, i.e., both of these phases can completely be converted to Ba(NO3)2. In contrast, the interaction of CO2 with pure BaO results in the formation of a BaCO3 layer that prevents to complete carbonation of the entire BaO film under the experimental conditions applied in these studies. However, these ‘‘carbonated’’ BaO layers readily react with NO2, and at elevated sample temperature even the carbonate layer is converted to nitrates. The importance of the metal oxide/metal interface in the chemistry on NOx storage-reduction catalysts was studied on BaO(\1 ML)/Pt(111) reverse model catalysts. In comparison to the clean Pt(111), new oxygen adsorption phases were identified on the BaO/Pt(111) surface that can be associated with oxygen atoms strongly adsorbed on Pt atoms at the peripheries of BaO particles. A simple kinetic model developed helped explain the observed thermal desorption results. The role of the oxide/metal interface in the reduction of Ba(NO3)2 was also substantiated in experiments where Ba(NO3)2/O/Pt(111) samples were exposed to CO at elevated sample temperature. The catalytic decomposition of the nitrate phase occurred as soon as metal sites opened up by the removal of interfacial oxygen via CO oxidation from the O/Pt(111) surface. The temperature for catalytic nitrate reduction was found to be significantly lower than the onset temperature of thermal nitrate decomposition. We gratefully acknowledge the US Department of Energy (DOE), Office of Science, Division of Chemical Sciences, Geosciences, and Biosciences for the support of this work. The research described in this paper was performed at the Environmental Molecular Sciences Laboratory (EMSL), a national user facility sponsored by the DOE Office of Biological and Environmental Research and located at the Pacific Northwest National Laboratory (PNNL). PNNL is operated for the US DOE by Battelle under contract number DE-AC05-76RL01830.

Environmental Molecular Sciences Laboratory↗

Pulse current auxiliary forging of sintered TiAl alloys based on the investigation of dynamic recrystallization and phase transformation

Highlights: • DRX kinetics model of sintered near γ TiAl alloy was established. • DRX process was accompanied by twinning when compressed at 1125 °C/0.1 s{sup − 1}. • Hot deformation has a promotion to the phase transformation of γ → α. - Abstract: Pulse current auxiliary forging is a new and efficient way to fabricate high performance TiAl based alloys. Dynamic recrystallization (DRX) and phase transformation have a pronounced effect on the subsequent deformability of TiAl based alloys. Hence, the dynamic recrystallization kinetics were determined at the temperature range of 1125–1275 °C and strain rate range of 0.1–0.001 s{sup − 1}. In addition, the DRX mechanism and phase transformation of TiAl based alloys were investigated to model the microstructure evolution in order to provide guidance for the pulse current auxiliary forging process. When the specimens were compressed at 1125 °C/0.1 s{sup − 1}, twinning played an important role in the DRX process. With increasing deformation temperature, the microstructure transformed from near gamma into duplex microstructure. In addition, according to the kinetics of dynamic recrystallization and phase transformation, fully DRXed TiAl based alloys were fabricated through pulse current auxiliary forging. The as-forged TiAl-based alloys show good superplastic performance with an elongation over 410%.

36 MATERIALS SCIENCE↗

Chemical kinetics of octane sensitivity in a spark-ignition engine (Chemical Kinetics of Octane Sensitivity in a Spark Ignition Engine)

This article uses a chemical kinetic modeling approach to study the influences of fuel molecular structure on Octane Sensitivity (OS) in Spark Ignition (SI) engines. Octane Sensitivity has the potential to identify fuels that can be used in next-generation high compression, turbocharged SI engines to avoid unwanted knocking conditions and extend the range of operating conditions that can be used in such engines. While the concept of octane numbers of different fuels has been familiar for many years, the variations of their values and their role in determining Octane Sensitivity have not been addressed previously in terms of the basic structures of the fuel molecules. In particular, the importance of electron delocalization on low temperature hydrocarbon reactivity and its role in determining OS in engine fuel is described here for the first time. Finally, the role of electron delocalization on fuel reactivity and Octane Sensitivity is illustrated for a very wide range of engine fuel types, including n-alkane, 1-olefin, n-alcohol, and n-alkyl benzenes, and the unifying features of these fuels and their common trends, using existing detailed chemical kinetic reaction mechanisms that have been collected and unified to produce an overall model with unprecedented capabilities.

30 DIRECT ENERGY CONVERSION↗

Pyrolysis of Cyclohexane and 1-Hexene at High Temperatures and Pressures—A Photoionization Mass Spectrometry Study

Cycloalkanes are important components of a wide range of fuels. However, there are few experimental data at simultaneously high temperatures and pressures similar to those found in practical systems. Such data are necessary for developing and testing chemical kinetic models. In this study, data relevant to cycloalkane pyrolysis were obtained from high repetition rate shock tube experiments coupled with synchrotron-based photoionization mass spectrometry diagnostics. The pyrolysis of cyclohexane was studied over 1270–1550 K and ~9 bar, while the more reactive primary decomposition product, 1-hexene, was studied at 1160–1470 K and ~5 bar. Insights into the decomposition of the parent molecules, the formation of primary products and the production of aromatic species were gained. Simulations were performed with models for cyclohexane and 1-hexene that were based on literature models. The results indicate that over several hundred microseconds reaction time at high pressures and temperatures the pyrolysis of cyclohexane is largely dominated by reactions initiated by cyclohexyl radicals. Furthermore, good agreement between the simulations and the experiments were observed for cyclohexane and 1-hexene with a modified version of the cyclohexane model. Conversely, the 1-hexene model did not reproduce the experimental observations.

1-alkenes↗

Validation of Wichner Predictive Model for Chronic Oxidation by Moisture of Nuclear Graphite

This document reports the results of multi-year extensive research on chronic oxidation of nuclear graphite by traces of moisture that might be present in the helium coolant of high-temperature gas-cooled reactors (HTGRs). The research focuses on four grades of nuclear graphite that are considered potential candidates for component fabrication in advanced nuclear reactors. The four grades are (1) PCEA (medium-fine grain, extruded, manufactured by GrafTech Inc., USA); (2) NBG-17 (medium-fine grain, vibrationally molded, manufactured by SGL Carbon, Germany/France); (3) IG-110 (microfine grain, isostatically molded, manufactured by Toyo Tanso, Japan); and (4) 2114 (microfine grain, isostatically pressed, manufactured by Mersen, USA). The experimental part of this research comprised three objectives that were consistently pursued for each graphite material: (1) accelerated oxidation tests in well-controlled conditions (temperature, gas composition) and analysis of chemical kinetic models and of their parameters for each graphite grade; (2) measurements of water vapor transport rates from H 2 O/He mixtures through thin graphite slabs for estimation of the water effective diffusivity specific to each material; (3) controlled oxidation of graphite specimens in H 2 O/He mixtures and characterization of the effects of oxidation temperatures on the density profiles of oxidized layers. These three groups of experimental results were interpreted using the model developed by Wichner et al. for predicting the penetration depth of oxidants (moisture or oxygen) as a function of oxidation conditions. Basically, starting from the kinetic and transport properties of each graphite grade, the ultimate goal was predicting the stable density profile in the oxidized layer and its dependence on the oxidation temperature. The profiles obtained this way were then compared against direct observations of oxidized layers in specimens obtained in well-controlled oxidation conditions. The general agreement between predictions and observations validated the general applicability of the Wichner model. The results obtained in this project were presented at annual Department of Energy–Advanced Reactor Technologies program reviews and at specialist meetings and conferences and were published in Oak Ridge National Laboratory (ORNL) reports and in professional journals. The results published in 2018 had an immediate impact in 2019, when a group from three US universities used ORNL kinetic data as input information in a complex simulation of chronic oxidation by moisture using the same graphite grades during long-time operation in a modular high-temperature gas-cooled reactor. Their conclusion, in line with the conclusions obtained at the end of this project, was that chronic graphite oxidation by moisture during normal operation conditions of an HTGR will be limited to a thin surface layer of graphite blocks at the reactor bottom, where the highest temperatures are. That will not affect the integrity of graphite components and reactor safety, provided the impurity concentration in the helium coolant is rigorously controlled and maintained at the design specifications.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Explicit least squares system parameter identification for exact differential input/output models

The equation error for a class of systems modeled by input/output differential operator equations has the potential to be integrated exactly, given the input/output data on a finite time interval, thereby opening up the possibility of using an explicit least squares estimation technique for system parameter identification. The paper delineates the class of models for which this is possible and shows how the explicit least squares cost function can be obtained in a way that obviates dealing with unknown initial and boundary conditions. The approach is illustrated by two examples: a second order chemical kinetics model and a third order system of Lorenz equations.

Pearson, A. E.↗