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 631 records · Page 35

Time-dependent numerical modeling of dust halo formation at comets

The evolution of gas and dust distributions following a spatially and temporally localized comet outburst was calculated using a hybrid kinetic-hydrodynamic method. In the inner coma the time-dependent continuity, momentum, and energy equations of the dusty gas flow were solved simultaneously using 12 dust sizes. Beyond 300 km a three-dimensional kinetic model was used to calculate the trajectory of each individual dust grain. It was found that following the onset of the comet outburst a gas-dust blast wave propagates outward in the inner coma. About 15 minutes after the increased gas and dust production was initiated at the nucleus, a new equilibrium was reached in the inner coma. The most important feature of this new steady state was the significant increase of the dust terminal velocities. These higher terminal velocity values resulted in larger apex distances for dust particles emitted during the outburst. The dust particles spend a relatively long time near their apex points; therefore, the outburst generates long-lasting distinct dust envelopes in front of the regular dust coma.

Gombosi, T. I.↗

On the structures and mapping of auroral electrostatic potentials

The mapping of magnetospheric and ionospheric electric fields in a kinetic model of magnetospheric-ionospheric electrodynamic coupling proposed for the aurora is examined. One feature is the generalization of the kinetic current-potential relationship to the return current region (identified as a region where the parallel drop from magnetosphere to ionosphere is positive); such a return current always exists unless the ionosphere is electrically charged to grossly unphysical values. A coherent phenomenological picture of both the low energy return current and the high energy precipitation of an inverted-V is given. The mapping between magnetospheric and ionospheric electric fields is phrased in terms of a Green's function which acts as a filter, emphasizing magnetospheric latitudinal spatial scales of order (when mapped to the ionosphere) 50 to 150 km. This same length, when multiplied by electric fields just above the ionosphere, sets the scale for potential drops between the ionosphere and equatorial magnetosphere.

Chiu, Y. T.↗

Investigation of low NOx staged combustor concept in high-speed civil transport engines

Levels of exhaust emissions due to high temperatures in the main combustor of high-speed civil transport (HSCT) engines during supersonic cruise are predicted. These predictions are based on a new combustor design approach: a rich burn/quick quench/lean burn combustor. A two-stage stirred reactor model is used to calculate the combustion efficiency and exhaust emissions of this novel combustor. A propane-air chemical kinetics model is used to simulate the fuel-rich combustion of jet fuel. Predicted engine exhaust emissions are compared with available experimental test data. The effect of HSCT engine operating conditions on the levels of exhaust emissions is also presented. The work described in this paper is a part of the NASA Lewis Research Center High-Speed Civil Transport Low NO(x) Combustor program.

Nguyen, Hung Lee↗

Investigation of low NOx staged combustor concept in high-speed civil transport engines

Levels of exhaust emissions due to high temperatures in the main combustor of high-speed civil transport (HSCT) engines during supersonic cruise are predicted. These predictions are based on a new combustor design approach: a rich burn/quick quench/lean burn combustor. A two-stage stirred reactor model is used to calculate the combustion efficiency and exhaust emissions of this novel combustor. A propane-air chemical kinetics model is used to simulate the fuel-rich combustion of jet fuel. Predicted engine exhaust emissions are compared with available experimental test data. The effect of HSCT engine operating conditions on the levels of exhaust emissions is also presented. The work described in this paper is a part of the NASA Lewis Research Center High-Speed Civil Transport Low NO(x) Combustor program.

Nguyen, Hung Lee↗

NUMERICAL INVESTIGATION OF FUEL PROPERTY EFFECTS ON MIXED-MODE COMBUSTION IN A SPARK-IGNITION ENGINE

In the present study, mixed-mode combustion of an E30 fuel in a direct-injection spark-ignition engine is numerically investigated at a fuel-lean operating condition using multidimensional computational fluid dynamics (CFD). A fuel surrogate matching Research Octane Number (RON) and Motor Octane Number (MON) of E30 is first developed using neural network based non-linear regression model. To enable efficient 3D engine simulations, a 164-species skeletal reaction mechanism incorporating NO, chemistry is reducedfrom a detailed chemical kinetic model. A hybrid approach that incorporates the G-equation model for tracking turbulentflame front, and the multi-zone well-stirred reactor model for predicting auto-ignition in the end gas, is employed to account for turbulent combustion interactions in the engine cylinder. Predicted in-cylinder pressure and heat release rate traces agree well with experimental measurements. The proposed modelling approach also captures moderated cyclic variability. Two different types of combustion cycles, corresponding to purely deflagrative and mixed-mode combustion, are observed. In contrast to the purely deflagrative cycles, mixed-mode combustion cycles feature early flame propagation followed by end-gas auto-ignition, leading to two distinctive peaks in heat release rate traces. The positive correlation between mixed-mode combustion cycles and early flame propagation is well captured by simulations. With the validated numerical setup, effects of NOx chemistry on mixed-mode combustion predictions are investigated. NOx chemistry is found to promote auto-ignition through residual gas recirculation, while the deflagrative flame propagation phase remains largely unaffected. Local sensitivity analysis is then performed to understand effects of physical and chemical properties of the fuel, i.e., heat of evaporation (HoV) and laminar flame speed (S-L). An increased HoV tends to suppress end-gas auto-ignition due to increased vaporization cooling, while the impact of HoV on flame propagation is insignificant. In contrast, an increased SL is found to significantly promote both flame propagation and auto-ignition. The promoting effect of SL on auto-ignition is not a direct chemical effect; it is rather caused by an advancement of the combustion phasing, which increases compression heating of the end gas.

detailed chemical kinetics↗

Pyrolysis of bio-derived dioxolane fuels: A ReaxFF molecular dynamics study

Alkyl-substituted 1,3-dioxolanes, including 4,5-dimethyl-2-pentan-3-yl-1,3-dioxolane (Fuel 1), 4,5-dimethyl-2-pentyl-1,3-dioxolane (Fuel 2), and 2-(heptan-3-yl)-4,5-dimethyl-1,3-dioxolane (Fuel 3), have been recently suggested as potential biodiesels. In this paper, we investigate the initial pyrolysis of the alkyl-substituted 1,3-dioxolanes at high temperatures using ReaxFF molecular dynamics (MD) simulations. We analyze the decomposition rate, reaction mechanism, and product distribution in the pyrolysis of the three alkyl 1,3-dioxlanes. The three fuels primarily decompose to 4,5-dimethyl-1,3-dioxolane radical and hydrocarbons derived from the alkyl side-chains. The further decomposition of 4,5-dimethyl-1,3-dioxolane radical primarily leads to 2-C 4 H 8 and CO 2 within a few decomposition steps. The hydrocarbon product distribution is significantly affected by the molecular structure of the alkyl side-chain, which would have a strong influence on the sooting tendency of these fuels. The ReaxFF simulations predict that the order of sooting tendency would be Fuel 3 > Fuel 1 > Fuel 2, which agrees with the measured sooting tendency trend. Based on the pyrolysis mechanism identified by ReaxFF, we propose a new alkyl dioxolane, 4-hexyl-5-methyl-1,3-dioxolane (Fuel 4), which might produce even less soot, by modifying the molecular structure of Fuel 2. Our ReaxFF simulation shows that Fuel 4 produce much less C 4 H 8 , an effective non-aromatic soot precursor, than Fuel 2. Moreover, more carbon atoms are bonded to each oxygen atom in Fuel 4 than Fuel 2, which would help reduce soot yield by removing more carbon atoms from the soot-producing pool of species. The major decomposition pathways identified in this work can be used to develop chemical kinetic models for 1,3-dioxolane based compounds, as biodiesel components, applicable to combustion engine simulations. We also demonstrate that the chemical kinetic insight offered by ReaxFF simulations can be used to design new fuel molecules with more desired properties.

09 BIOMASS FUELS↗

Experimental and theoretical studies of NO/x/ formation in a jet-stirred combustor.

As an approach to testing kinetic models which couple NO formation to combustion reactions, a jet-stirred combustor was used to study NO(x) formation under kinetically controlled combustion conditions. The experimental data were compared with kinetic calculations for a well-stirred reactor. The agreement between theory and experiment was quite good for the combustion of hydrogen/air and carbon monoxide/air, for which detailed mechanisms were used. A quasi-global mechanism for the combustion of propane/air resulted in substantial underprediction of NO(x) formation.

Engleman, V. S.↗

On The Development of Biophysical Models for Space Radiation Risk Assessment

Experimental techniques in molecular biology are being applied to study biological risks from space radiation. The use of molecular assays presents a challenge to biophysical models which in the past have relied on descriptions of energy deposition and phenomenological treatments of repair. We describe a biochemical kinetics model of cell cycle control and DNA damage response proteins in order to model cellular responses to radiation exposures. Using models of cyclin-cdk, pRB, E2F's, p53, and GI inhibitors we show that simulations of cell cycle populations and GI arrest can be described by our biochemical approach. We consider radiation damaged DNA as a substrate for signal transduction processes and consider a dose and dose-rate reduction effectiveness factor (DDREF) for protein expression.

Cucinotta, F. A.↗

Insights on keto-hydroperoxide formation from O 2 addition to the beta-tetrahydrofuran radical

Cyclic ethers provide an interesting case study of low-temperature oxidation chemistry, especially in relevance to biofuels. A recent experimental study revealed new questions regarding the low-temperature oxidation mechanism of tetrahydrofuran (THF) concerning the formation of keto-hydroperoxides. In particular, keto-hydroperoxides originating from the THF-β radical were not captured accurately by current literature models, motivating this work to calculate the energetics of the first and second O 2 -addition pathways for THF radicals. Electronic structure calculations at the CCSD(T)/cc-pV∞Z//M06-2X/cc-pVTZ level of theory were used to generate potential energy surfaces for the α-C and β-C THF radicals and subsequent pathways to the formation of the keto-hydroperoxide isomers. These are the first theoretical calculations of the second O 2 -addition radical pathways for the THF-β radical. Results from the theorical work provided further insight into the low-temperature oxidation of THF. This included identifying the pathways most likely to form the keto-hydroperoxide isomers observed in prior experimental work; and detecting that the shortcomings in prior models are likely due to uncertainties in R + THF abstraction reaction rates. Here, these conclusions will motivate future work for accurate THF kinetic model development.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Performance characterization of nanofiltration, reverse osmosis, and ion exchange technologies for acetic acid separation

A major obstacle to widespread implementation of bio-based fuels and chemicals is the lack of efficient and costeffective separation methods. To purify acetic acid produced by biochemical conversion of biomass via anaerobic digestion, this work employs two commonly used separation technologies: (1) ion-exchange (IX) resin and (2) pressure-driven membranes. This study tested five commercially available strong- and weak-base anion-exchange resins and five commercially available nanofiltration (NF) and reverse osmosis (RO) membranes. The pH of the feed solution significantly affected adsorption capacity. At pH 6.3, a strong-base IX resin (IRN-78) performed best (95.1% acetate removal). With strong-base IX resins, the Langmuir isotherm model fit well, whereas for weak-base IX resins, the Freundlich isotherm provided a better fit. A pseudo-second-order kinetic model fit well for both IRN-78 and IRA-67. Regarding membrane separation, RO (BW30XFR membrane) achieved the highest rejection (98.6% acetate rejection), whereas an NF membrane (NF*) achieved the best combination of permeate flux (105 L/(h∙m 2 )) and rejection (83.1% acetate rejection). For membrane performance, the experimental data were fit using the solution diffusion model. Increased pH in the feed solution lowered permeate flux but increased acetic acid rejection. Finally, then the acetic acid concentration in the feed solution increased, both permeate flux and acetic acid rejection decreased for membrane NF*.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Refined Outgassing Model for the LiH/LiOH System

Previous models for H 2 and H 2 O outgassing from moisture exposed LiH specimens were based on temperature programmed desorption/decomposition (TPD) experiments and analyses that yielded kinetic descriptions which were specific to defined sets of initial conditions. Subsequently the inputs to the kinetic models were refined to correlate with infrared spectra of the adsorbed species. In this current enhancement, better agreement between prediction and experiment was obtained through: 1) refining the reaction mechanisms for some of the decomposition reactions; 2) developing the potential range of kinetic parameters by simultaneous fitting of all the TPD outgassing curves obtained at different heating rates; and, 3) limiting the scope of possible kinetic parameters to an optimized subset by matching the outgassing prediction curves at 3 different temperatures to the outcomes from the accurate but cumbersome, non-intuitive, and effort demanding isoconversional analysis.

36 MATERIALS SCIENCE↗

Interactions of Hydrazine and Blowby Gases

The interactions between hydrazine and blowby gases from pyrovalves was explored in this research project. Investigating the decomposition chemistry of hydrazine through detailed chemical kinetic modeling is a project started last summer while participating in the Summer Faculty Fellowship program. During the 1999-2000 academic year, the chemical kinetic mechanism for hydrazine decomposition developed while a SFF at NASA's White Sands Test Facility was further revised and validated against the limited experimental data in the literature. This mechanism was then used in assessing the effects of blowby gas species on hydrazine decomposition. The combustion products introduced into the fuel line by pyrovalve actuation consist primarily of hydrogen gas. Hydrogen is also a product of the decomposition of hydrazine. Additional gaseous chemical species are introduced into the fuel, as well as metals and metal salts that deposit onto the walls of the fuel line. The deposition process is undoubtedly very rapid, and exothermic. Therefore, the major focus of this summer's work was examining the effects of hydrogen presence on hydrazine decomposition, with some representative calculations including the remaining gaseous species found to exist in blowby gases. Since hydrogen is a product of hydrazine decomposition, all reactions necessary to evaluate its effect on hydrazine decomposition chemistry were in the original mechanism developed. However, the mechanism needed to be considerably expanded to include the reactions of the other gaseous blowby species with hydrazine, all the intermediate species formed in its decomposition, and each other. The expanded mechanism consists of 70 species interacting via a network of 452 reactions. Calculations with molecular hydrogen introduced into hydrazine gas in an inert bath gas indicate that H2 presence as an initial reactant in substantial amounts can dramatically impact the decomposition process for hydrazine. The other gaseous blowby species (CO, CO2, H2O, CH4, O2, and N2) were found to have little effect compared to the inclusion of hydrogen itself as an initial reagent. This result is undoubtedly due, in part, to the fact that the blowby gas used in these calculations consisted of 94.6% H2. A more rigorous examination of the behavior of the full detailed mechanism under a variety of conditions was not performed.

Meagher, Nancy E.↗

Elucidating the photodissociation fingerprint and quantifying the determination of organic hydroperoxides in gas-phase autoxidation

Hydroperoxides are formed in the atmospheric oxidation of volatile organic compounds, in the combustion autoxidation of fuel, in the cold environment of the interstellar medium, and also in some catalytic reactions. They play crucial roles in the formation and aging of secondary organic aerosols and in fuel autoignition. However, the concentration of organic hydroperoxides is seldom measured, and typical estimates have large uncertainties. In this work, we developed a mild and environmental-friendly method for the synthesis of alkyl hydroperoxides (ROOH) with various structures, and we systematically measured the absolute photoionization cross-sections (PICSs) of the ROOHs using synchrotron vacuum ultraviolet-photoionization mass spectrometry (SVUV-PIMS). A chemical titration method was combined with an SVUV-PIMS measurement to obtain the PICS of 4-hydroperoxy-2-pentanone, a typical molecule for combustion and atmospheric autoxidation ketohydroperoxides (KHPs). We found that organic hydroperoxide cations are largely dissociated by loss of OOH. This fingerprint was used for the identification and accurate quantification of the organic peroxides, and it can therefore be used to improve models for autoxidation chemistry. The synthesis method and photoionization dataset for organic hydroperoxides are useful for studying the chemistry of hydroperoxides and the reaction kinetics of the hydroperoxy radicals and for developing and evaluating kinetic models for the atmospheric autoxidation and combustion autoxidation of the organic compounds.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Probabilistic effects of porosity and chemical kinetics on the shock initiation of an octahydro-1,3,5,7-tetranitro-1,3,5,7-tetrazocine (HMX) based explosive

In this work, we investigate the effects of porosity and chemical reactivity on the shock initiation of an HMX (octahydro-1,3,5,7-tetranitro-1,3,5,7-tetrazocine)-based polymer bonded explosive due to thin flyer impact using two-dimensional mesoscale simulations. Microstructural samples are computationally generated, and the calculated shock response is compared to experiments. The sample dimensions are 0.1×1.05 mm 2 , and each sample accounts for grain scale heterogeneities, including randomized grain morphology and pore spacing. Simulations are performed using the Lawrence Livermore National Laboratory (LLNL) arbitrary Lagrangian–Eulerian multi-physics code ALE3D. Chemistry and equation of states are handled with LLNL's thermochemical code Cheetah. We employ a modified Arrhenius-based chemical kinetics model that scales the temperature exponent fitting parameter with flyer shock conditions. This new model is calibrated to the mean response of experimental shock initiation thresholds for three separate Kapton flyer thicknesses. A simplified shock initiation criterion is proposed, which classifies the macroscale initiation response based on cross-sectional temperature profiles of the shocked sample. This criterion is used to further investigate the probabilistic response of the sample to variations in the mean porosity and chemical kinetics rate, and we find an inverse correlation between the shock pressure and the confidence interval of the experimental threshold velocity. Furthermore, the methodology presented in this work for quantifying the probabilistic nature of shock initiation is material agnostic. These mesoscale studies may be used to inform continuum scale models and expedite the development of future, tailored energetic formulations.

36 MATERIALS SCIENCE↗

Implementation of new mixture rules has a substantial impact on combustion predictions for H 2 and NH 3

Complex-forming reactions comprise a substantial fraction of all important combustion reactions and are central to combustion behavior. Despite being often called “pressure-dependent” reactions, their rate constants depend on not only the pressure but also the composition. While modern combustion codes allow arbitrarily high accuracy in treating pressure dependence, recent work has consistently demonstrated dramatic failures of essentially all available treatments of mixture dependence. In situations where mixture dependence is treated at all, it is inevitably treated through specification of pressure-dependent rate constants for a set of pure bath gases, which are then combined to estimate the rate constant in a mixture via a “mixture rule.” While there had been a generally unquestioning confidence in these mixture rules, they had, in reality, been scarcely tested until the last decade, when comparisons against master equation calculations revealed order-of-magnitude errors for important pressure-dependent reactions. New mixture rules, based on the reduced pressure, have recently been proposed and shown to reproduce master equation calculations for broad classes of complex-forming reactions very accurately. Here, in this work, we present an implementation of one such new mixture rule (“LMR-R”) in Cantera and then use it to enable simulations that use new high-accuracy ab initio data for individual bath gases (for the first time, since codes previously could not accommodate the complex bath gas dependence). Demonstrations focus on combustion of H 2 and NH 3 , where (1) high-accuracy ab initio data are available and (2) the impact is expected to be large due to the high fractions of efficient colliders (e.g., H 2 O and NH 3 ) in the burned and unburned gases. Indeed, we find the impact of this treatment to be substantial and may explain previous modeling difficulties for these important carbon-free fuels, particularly for NH 3 , whose extraordinarily high third-body efficiency (~20) is often omitted from kinetic models.

Ammonia↗

Pore-Scale Study on the Positive Feedback Between Stress and Porosity Caused by Pressure Solution in Porous Media

Pressure solution is an important process in the evolution of sedimentary rocks, which provide storage space for most of our petroleum resources. It directly influences the generation, migration, and storage of petroleum fluids in subsurface sedimentary rocks. Here, in this paper, we develop a pore-scale, mechanochemical model to demonstrate a possible positive feedback between the local porosity and pore surface stress, in which a higher local porosity causes a higher local pore surface stress, thus enhancing pressure solution and consequently further increasing the local porosity. Pore surface stress represents stress on a solid grain adjacent to a pore. Specifically, the pore-scale, mechanochemical model directly simulates the stress distribution over solid and pore surfaces using a finite element model. The dissolution of solids at the solid-pore interfaces under a far-from-equilibrium condition is simulated using a first-order kinetics model that accounts for the local stress distribution. The updated pore geometry, caused by pore surface dissolution, is then used in the stress simulation in the next numerical iteration. Two types of porous media, the Oriskany sandstone and an artificial porous medium with spherical pores, were tested in the mechanochemical simulation. The positive stress-porosity feedback during pressure solution was observed in both samples. In addition, the model quantitatively illustrated the distribution of local mineral dissolution rates on all pore surfaces, as well as its relation to the effective mineral dissolution rate of the entire sample. Based on the comparison between the two porous media, the local mineral dissolution was regulated by pore space distribution, geometry, and coalescence during pressure solution. This work is the first that uses direct, pore-scale numerical simulation to demonstrate the positive stress-porosity feedback during pressure solution, which has the potential to advance the understanding of the mechanical-chemical (MC) coupling in many geological processes that are relevant to subsurface energy systems, such as the recovery of petroleum hydrocarbons and geothermal energy.

CT scanning↗

Concepts Relevant for the Kinetic Analysis of Reversible Reaction Systems

The net rate of a reversible chemical reaction is the difference between unidirectional rates of traversal along forward and reverse reaction paths. In a multistep reaction sequence, the forward and reverse trajectories, in general, are not the microscopic reverse of one another; rather, each unidirectional route is comprised of distinct rate-controlling steps, intermediates, and transition states. Consequently, traditional descriptors of rate (e.g., reaction orders) do not reflect intrinsic kinetic information but instead conflate unidirectional contributions determined by (i) the microscopic occurrence of forward/reverse reactions (i.e., unidirectional kinetics) and (ii) the reversibility of reaction (i.e., nonequilibrium thermodynamics). This review aims to provide a comprehensive resource of analytical and conceptual tools which deconvolute the contributions of reaction kinetics and thermodynamics to disambiguate unidirectional reaction trajectories and precisely identify rate- and reversibility-controlling molecular species and steps in reversible reaction systems. The extrication of mechanistic and kinetic information from bidirectional reactions is accomplished through equation-based formalisms (e.g., De Donder relations) grounded in principles of thermodynamics and interpreted in the context of theories of chemical kinetics developed in the past 25 years. Here, the aggregate of mathematical formalisms detailed herein is general to thermochemical and electrochemical reactions and encapsulates a diverse body of scientific literature encompassing chemical physics, thermodynamics, chemical kinetics, catalysis, and kinetic modeling.

42 ENGINEERING↗

A comprehensive model to determine the effects of temperature and species fluctuations on reaction rates in turbulent reaction flows

The use of silane (SiH4) as an effective ignitor and flame stabilizing pilot fuel is well documented. A reliable chemical kinetic mechanism for prediction of its behavior at the conditions encountered in the combustor of a SCRAMJET engine was calculated. The effects of hydrogen addition on hydrocarbon ignition and flame stabilization as a means for reduction of lengthy ignition delays and reaction times were studied. The ranges of applicability of chemical kinetic models of hydrogen-air combustors were also investigated. The CHARNAL computer code was applied to the turbulent reaction rate modeling.

Magnotti, F.↗