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,153 records · Page 64

Role of Macromolecular Structure in the Ultrafast Energy and Electron Transfer Dynamics of a Light-Harvesting Polymer

Ultrafast energy and electron transfer (EnT and ET, respectively) are characterized in a light-harvesting assembly based on a π-conjugated polymer (poly(fluorene)) functionalized with broadly absorbing pendant organic isoindigo (iI) chromophores using a combination of femtosecond transient absorption spectroscopy and large-scale computer simulation. Photoexcitation of the π-conjugated polymer leads to near-unity quenching of the excitation through a combination of EnT and ET to the iI pendants. The excited pendants formed by EnT rapidly relax within 30 ps, whereas recombination of the charge-separated state formed following ET occurs within 1200 ps. A computer model of the excited-state processes is developed by combining all-atom molecular dynamics simulations, which provides a molecular-level view of the assembly structure, with a kinetic model that accounts for the multiple excited-state quenching pathways. In conclusion, direct comparison of the simulations with experimental data reveals that the underlying structure has a dramatic effect on the partitioning between EnT and ET in the polymer assembly, where the distance and orientation of the pendants in relation to the backbone serve to direct the dominant quenching pathway.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

The tarantula toxin GxTx detains K + channel gating charges in their resting conformation

Allosteric ligands modulate protein activity by altering the energy landscape of conformational space in ligand–protein complexes. Here we investigate how ligand binding to a K + channel’s voltage sensor allosterically modulates opening of its K + -conductive pore. The tarantula venom peptide guangxitoxin-1E (GxTx) binds to the voltage sensors of the rat voltage-gated K + (Kv) channel Kv2.1 and acts as a partial inverse agonist. When bound to GxTx, Kv2.1 activates more slowly, deactivates more rapidly, and requires more positive voltage to reach the same K + -conductance as the unbound channel. Further, activation kinetics are more sigmoidal, indicating that multiple conformational changes coupled to opening are modulated. Single-channel current amplitudes reveal that each channel opens to full conductance when GxTx is bound. Inhibition of Kv2.1 channels by GxTx results from decreased open probability due to increased occurrence of long-lived closed states; the time constant of the final pore opening step itself is not impacted by GxTx. When intracellular potential is less than 0 mV, GxTx traps the gating charges on Kv2.1’s voltage sensors in their most intracellular position. Gating charges translocate at positive voltages, however, indicating that GxTx stabilizes the most intracellular conformation of the voltage sensors (their resting conformation). Kinetic modeling suggests a modulatory mechanism: GxTx reduces the probability of voltage sensors activating, giving the pore opening step less frequent opportunities to occur. This mechanism results in K + -conductance activation kinetics that are voltage-dependent, even if pore opening (the rate-limiting step) has no inherent voltage dependence. We conclude that GxTx stabilizes voltage sensors in a resting conformation, and inhibits K + currents by limiting opportunities for the channel pore to open, but has little, if any, direct effect on the microscopic kinetics of pore opening. The impact of GxTx on channel gating suggests that Kv2.1’s pore opening step does not involve movement of its voltage sensors.

Physiology↗

Phosphocreatine kinetics at the onset of contractions in skeletal muscle of MM creatine kinase knockout mice

Phosphocreatine (PCr) depletion during isometric twitch stimulation at 5 Hz was measured by (31)P-NMR spectroscopy in gastrocnemius muscles of pentobarbital-anesthetized MM creatine kinase knockout (MMKO) vs. wild-type C57B (WT) mice. PCr depletion after 2 s of stimulation, estimated from the difference between spectra gated to times 200 ms and 140 s after 2-s bursts of contractions, was 2.2 +/- 0.6% of initial PCr in MMKO muscle vs. 9.7 +/- 1.6% in WT muscles (mean +/- SE, n = 7, P < 0.001). Initial PCr/ATP ratio and intracellular pH were not significantly different between groups, and there was no detectable change in intracellular pH or ATP in either group after 2 s. The initial difference in net PCr depletion was maintained during the first minute of continuous 5-Hz stimulation. However, there was no significant difference in the quasi-steady-state PCr level approached after 80 s (MMKO 36.1 +/- 3.5 vs. WT 35.5 +/- 4.4% of initial PCr; n = 5-6). A kinetic model of ATPase, creatine kinase, and adenylate kinase fluxes during stimulation was consistent with the observed PCr depletion in MMKO muscle after 2 s only if ADP-stimulated oxidative phosphorylation was included in the model. Taken together, the results suggest that cytoplasmic ADP more rapidly increases and oxidative phosphorylation is more rapidly activated at the onset of contractions in MMKO compared with WT muscles.

NASA Discipline Cell Biology↗

Resting myocardial blood flow quantification using contrast-enhanced magnetic resonance imaging in the presence of stenosis: A computational fluid dynamics study

Purpose: The extent to which atherosclerotic plaques affect contrast agent (CA) transport in the coronary arteries and, hence, quantification of myocardial blood flow (MBF) using magnetic resonance imaging (MRI) is unclear. The purpose of this work was to evaluate the influence of plaque induced stenosis both on CA transport and on the accuracy of MBF quantification. Methods: Computational fluid dynamics simulations in a high-detailed realistic vascular model were employed to investigate CA bolus transport in the coronary arteries. The impact of atherosclerosis was analyzed by inserting various medium- to high-grade stenoses in the vascular model. The influence of stenosis morphology was examined by varying the stenosis shapes but keeping the area reduction constant. Errors due to CA bolus transport were analyzed using the tracer-kinetic model MMID4. Results: Dispersion of the CA bolus was found in all models and for all outlets, but with a varying magnitude. The impact of stenosis was complex: while high-grade stenoses amplified dispersion, mild stenoses reduced the effect. Morphology was found to have a marked influence on dispersion for a small number of outlets in the post-stenotic region. Despite this marked influence on the concentration–time curves, MBF errors were less affected by stenosis. In total, MBF was underestimated by −7.9% to −44.9%. Conclusions: The presented results reveal that local hemodynamics in the coronary vasculature appears to have a direct impact on CA bolus dispersion. Inclusion of atherosclerotic plaques resulted in a complex alteration of this effect, with both degree of area reduction and stenosis morphology affecting the amount of dispersion. This strong influence of vascular transport effects impairs the accuracy of MRI-based MBF quantification techniques and, potentially, other bolus-based perfusion measurement techniques like computed tomography perfusion imaging.

60 APPLIED LIFE SCIENCES↗

Comparison of semikinetic and generalized transport models of the polar wind

A rigorous comparison is made of a semikinetic model and the bi-Maxwellian-based 16-moment transport equations for polar wind conditions. The comparisons show remarkably close agreement in the corresponding predictions for the altitudinal variation of the various plasma parameters. These results indicate that the bi-Maxwellian-based transport equations are a powerful tool for modeling a wide range of thermal plasma flows throughout the solar-terrestrial environment. In the collisionless limit, the results indicate that the 16-moment transport equations are equivalent to the kinetic models commonly used in solar and polar wind studies.

Demars, H. G.↗

Progress in Validation of Wind-US for Ramjet/Scramjet Combustion

Validation of the Wind-US flow solver against two sets of experimental data involving high-speed combustion is attempted. First, the well-known Burrows- Kurkov supersonic hydrogen-air combustion test case is simulated, and the sensitively of ignition location and combustion performance to key parameters is explored. Second, a numerical model is developed for simulation of an X-43B candidate, full-scale, JP-7-fueled, internal flowpath operating in ramjet mode. Numerical results using an ethylene-air chemical kinetics model are directly compared against previously existing pressure-distribution data along the entire flowpath, obtained in direct-connect testing conducted at NASA Langley Research Center. Comparison to derived quantities such as burn efficiency and thermal throat location are also made. Reasonable to excellent agreement with experimental data is demonstrated for key parameters in both simulation efforts. Additional Wind-US feature needed to improve simulation efforts are described herein, including maintaining stagnation conditions at inflow boundaries for multi-species flow. An open issue regarding the sensitivity of isolator unstart to key model parameters is briefly discussed.

Engblom, William A.↗

New insights into fuel blending effects: Intermolecular chemical kinetic interactions affecting autoignition times and intermediate-temperature heat release

Fuel blending effects on chemically-dominated fuel properties, such as gasoline anti-knock quality, are influenced by fundamental chemical kinetic interactions between the blending agent and the base fuel. Historically, quantification of such interactions has focused on changes to the radical pool, including $\dot{O}$H and HO 2 , while intermolecular interactions pertaining to carbonated, non-fuel-specific intermediates are typically overlooked. In this regard, this work aims to derive new insight into intermolecular chemical kinetic interactions that are intrinsic to fuel blending effects via a case study on blends of 0–30% ethanol (by volume) into FGF-LLNL (a multi-component gasoline surrogate for FACE-F research gasoline) using a rapid compression machine at a diluted/stoichiometric fuel loading, compressed pressure of 40 bar and low- to intermediate-temperature regimes that are representative of boosted SI engine operation. Ethanol blending effects on the intermediate temperature heat release (ITHR) of FGF-LLNL are characterized using experimental measurements, where ethanol is found to promote the extent of ITHR and suppress the transition from ITHR to main ignition. Chemical kinetic modeling is undertaken using a recently updated gasoline surrogate model. Sensitivity analyses on ITHR characteristics further corroborate the ethanol blending effects, and highlight the significant dependence of ITHR on both fuel-specific and non-fuel-specific reactions. An approach allowing comprehensive characterization of the complex intermolecular chemical kinetic interactions between constitutes in a fuel blend is then proposed. Application of the approach to FGF-LLNL/E0–E30 reveals that ethanol perturbs the heat release and autoignition characteristics of FGF-LLNL not only by directly changing the $\dot{O}$H and HO 2 radical pools via fuel-specific reactions, but also through intermolecular interactions where participating intermediates can be produced and consumed by various sub-chemistries. Disabling the intermolecular interactions in carbonated species between ethanol and FGF-LLNL sub-chemistries leads to somewhat slower ITHR evolution and lower ignition reactivity. The role of the individual intermolecular interaction is also characterized using the proposed approach. Lastly, implications of intermolecular interactions for future studies that aim to improve model performance are highlighted, where it is found that further investigations on the core C0-C4 chemistry are needed for developing highly accurate chemistry models for complex fuel blends.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

High-Pressure Rate Rules for Ether Alkylperoxy Radical Isomerization

The first isomerization reaction of an alkylperoxy (RO 2 ) radical holds significant importance in low-temperature oxidation, as it governs the branching ratios of the hydroperoxyalkyl (QOOH) radicals, which influence the competition between the chain-propagation and chain-branching reactions. In this study, we systematically calculated high-pressure rate rules for the RO 2 isomerization reaction of monoethers, exploring 5-, 6-, 7-, and 8-membered ring transition states. Primary, secondary, and tertiary carbon sites, where both the abstracting peroxy group and the abstracted hydrogen are located, were considered, with particular emphasis on distinguishing between secondary carbons adjacent (alpha) and nonadjacent to the ether functional group. Using the G4//B3LYP/6-311++G(2df,2pd) level of theory and the transition state theory, we estimated the rate constants and the Arrhenius coefficient for over 120 possible isomerization reactions. We examined the effect of ring size and ring atoms, revealing that 6- and 7-membered ring isomerizations were generally the fastest. The impact of the ether functional group on transition states was investigated by comparing reactions with identical ring size, peroxy, and radical positions, but with the ether functional group positioned either outside (i.e., out) or inside (i.e., in) the transition state ring, leading to differences in the rate constants. When comparing to analogous alkane rate constants, differences of up to an order of magnitude were observed, underscoring the need for caution when assigning rate rules by analogy. We applied our rate constants in the di-iso-butyl ether kinetic model and evaluated their influence on low-temperature chemistry finding that they altered the branching ratios by up to a factor of 9, highlighting the significance of site-specific rate constants for more accurate low-temperature modeling.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Model for the Solid–Liquid Interfacial Free Energy at High Pressures

The free energy involved in the formation of an interface between two phases (e.g., a solid–liquid interface) is referred to as the interfacial free energy. For the case of solidification, the interfacial free energy dictates the height of the energy barrier required to nucleate stable clusters of the newly forming solid phase and is essential for producing an accurate solidification kinetics model using classical nucleation theory (CNT)-based methods. While various methods have been proposed for modeling the interfacial free energy for solid–liquid interfaces in prior literature, many of these formulations involve making restrictive assumptions or approximations, such as the system being at or near equilibrium (i.e., the system temperature is approximately equal to the melt temperature) or that the system is at pressures close to atmospheric. However, these approximations and assumptions may break down in highly non-equilibrium situations, such as in dynamic-compression experiments where metastable liquids that are undercooled by hundreds of kelvin or over-pressurized by several gigapascals or more are formed before eventually solidifying. We derive a solid–liquid interfacial free-energy model for such high-pressure conditions by considering the enthalpies of interactions between pairs of atoms or molecules. We also consider the contribution of interface roughness (disordering) by incorporating a multilayer interface model known as the Temkin n-layer model. Our formulation is applicable to a diverse variety of materials, and we demonstrate it by developing models specifically for two different materials: water and gallium. We apply our interfacial free-energy formulation to CNT-based kinetics simulations of several suites of dynamic-compression experiments that cause liquid water to solidify to the high-pressure solid polymorph ice VII and have found good agreement to the observed kinetics with only minor empirical fitting.

36 MATERIALS SCIENCE↗

Application of evolved gas analysis to cold-cap reactions of melter feeds for nuclear waste vitrification

The melter feed (a mixture of nuclear waste and glass forming and modifying additives) experiences multiple gas-evolving reactions in an electrical glass-melting furnace during the vitrification of nuclear wastes. Foams from the residual gases can significantly alter the melting rate through mass and heat transfers. Here, we employed the thermogravimetry-gas chromatography-mass spectrometry (TGA-GC-MS) combination to perform quantitative evolved gas analysis (EGA) and developed a simple calibration model which correlates the overall mass loss rate with the evolution rates for individual gases. The model parameters are obtained from the least squares analysis, assuming that the gas-evolving reactions are independent. Finally, the EGA adds the ‘chemical identity’ to the reactions indicated by the ‘phenomenological’ kinetic model.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Demonstration of a viable quantitative theory for interplanetary type II radio bursts

Between 29 November and 1 December 2013 the two widely separated spacecraft STEREO A and B observed a long lasting, intermittent, type II radio burst for the extended frequency range ≈ 4 MHz to 30 kHz, including an intensification when the shock wave of the associated coronal mass ejection (CME) reached STEREO A. We demonstrate for the first time our ability to quantitatively and accurately simulate the fundamental (F) and harmonic (H) emission of type II bursts from the higher corona (near 11 solar radii) to 1 AU. Our modeling requires the combination of data-driven three-dimensional magnetohydrodynamic simulations for the CME and plasma background, carried out with the BATS-R-US code, with an analytic quantitative kinetic model for both F and H radio emission, including the electron reflection at the shock, growth of Langmuir waves and radio waves, and the radiations propagation to an arbitrary observer. The intensities and frequencies of the observed radio emissions vary hugely by factors ≈ 10{sup 6} and ≈ 10{sup 3}, respectively; the theoretical predictions are impressively accurate, being typically in error by less than a factor of 10 and 20 %, for both STEREO A and B. We also obtain accurate predictions for the timing and characteristics of the shock and local radio onsets at STEREO A, the lack of such onsets at STEREO B, and the z-component of the magnetic field at STEREO A ahead of the shock, and in the sheath. Very strong support is provided by these multiple agreements for the theory, the efficacy of the BATS-R-US code, and the vision of using type IIs and associated data-theory iterations to predict whether a CME will impact Earth’s magnetosphere and drive space weather events.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Organic synthesis in the outer Solar System: Recent laboratory simulations for Titan, the Jovian planets, Triton and comets

We tabulate the most abundant gases and their radiation yields, for two experimental pressures: 0.24 mb, more relevant to upper atmosphere excitation, and 17 mb, more relevant to tropospheric, cosmic ray excitation. The yields computed in the 0.24 mb experiment combined with measured electronic fluxes and a simple, eddy diffusion model of Titan's atmosphere predict abundances of detected molecules in agreement with those found by Voyager and for heavier products, in somewhat better agreement with observation than photochemical absolute reaction rate kinetics models. All Voyager organics are accounted for and no detectable products are found that Voyager did not detect. A striking increase of products with multiple bonds is found with decreasing pressure. Hydrocarbon abundances decline slowly with increasing carbon number. Additionally, we list preliminary estimates for the yield of the heteropolymer, which seems to be produced in a quantity comparable (in moles of C+N consumed) to the total amount of gaseous product. The production rate required to sustain Titan's haze against sedimentation also indicates yields of this order. As can be seen from the table, over 10(exp 9) years substantial amounts of these products can accumulate on the surface -- ranging from cm thickness for the (C+N equals 4) species to a meter or more for HCN and C2H2; we also expect a meter or more of tholins. Similar analyses have been or are being done for the Jovian planets and Triton. Charged particle irradiation of hydrocarbon clathrates or mixed hydrocarbon/water ices produces a range of organic products, reddening and darkening of the ices and characteristic infrared spectra. From such spectra, the predicted emission by fine particles in cometary comae well-matches the observed 3.4 micron emission spectra of Comet Halley and other recent comets. Heliocentric evolution of organic emission features in comets is predicted. Organic products of such ice irradiation may account for colors and albedos on some of the satellites in the outer solar system, especially Triton and Pluto, where solid methane is known to exist.

Sagan, C.↗

Biophysics Representation of the Two-Hit Model of Alzheimer's Disease for the Exploration of Late CNS Risks from Space Radiation

A concern for long-term space travel outside the Earth s magnetic field is the late effects to the central nervous system (CNS) from galactic cosmic ray (GCR) or solar particle events (SPE). Human epidemiology data is severely limited for making CNS risk estimates and it is not clear such effects occur following low LET exposures. We are developing systems biology models based on biological information on specific diseases, and experimental data for proton and heavy ion radiation. A two-hit model of Alzheimer s disease (AD) has been proposed by Zhu et al.(1), which is the framework of our model. Of importance is that over 50% of the US population over the age of 75-y have mild to severe forms of AD. Therefore we recommend that risk assessment for a potential AD risk from space radiation should focus on the projection of an earlier age of onset of AD and the prevention of this possible acceleration through countermeasures. In the two-hit model, oxidative stress and aberrant cell cycle-related abnormalities leading to amyloid-beta plaques and neurofibrillary tangles are necessary and invariant steps in AD. We have formulated a stochastic cell kinetics model of the two-hit AD model. In our model a population of neuronal cells is allowed to undergo renewal through neurogenesis and is susceptible to oxidative stress or cell cycle abnormalities with age-specific accumulation of damage. Baseline rates are fitted to AD population data for specific ages, gender, and for persons with an apolipoprotein 4 allele. We then explore how low LET or heavy ions may increase either of the two-hits or neurogenesis either through persistent oxidative stress, direct mutation, or through changes to the micro-environment, and suggest possible ways to develop accurate quantitative estimates of these processes for predicting AD risks following long-term space travel.

Cucinotta, Francis A.↗

Application of soot carbonization kinetics to deduce meaningful soot formation rates in premixed flat flames

Soot formation rates measured in fuel-rich premixed flat flames are frequently used to calibrate or validate chemical kinetic models of soot formation. Unfortunately, these flames feature an extended region of soot precursor particle inception and carbonization that complicates interpretation of soot measurements and leads to a fundamental inconsistency in the nature of the soot material that is modeled versus what is being measured when using non-intrusive, optical techniques. In the work presented here, previously reported data on two canonical sooting ethylene-air premixed flames at 1 atm pressure are interpreted via a new analysis approach that combines soot optical dispersion coefficient measurements with soot carbonization kinetics. This analytical approach explicitly accounts for the production of poorly ordered soot precursor particle mass and its carbonization over time in the flames, providing a clear distinction between the formation rate of precursor particles and their transformation into ordered, solid soot particulate mass. In particular, the results of the analysis show that the precursor particles form much earlier in the flame than the majority of the carbonized soot and their formation rate is two to three times faster than that of ordered soot. The results also show that particle agglomeration begins when the particles are at an intermediate state of carbonization. In conclusion, these results offer a valuable new interpretation of these important datasets and should lead to substantial improvements in the development and calibration of quantitative soot models.

Soot formation↗

Evaluation of Water Injection Effect on NO(x) Formation for a Staged Gas Turbine Combustor

NO(x) emission control by water injection on a staged turbine combustor (STC) was modeled using the KIVA-2 code with modification. Water is injected into the rich-burn combustion zone of the combustor by a single nozzle. Parametric study for different water injection patterns was performed. Results show NO(x) emission will decrease after water being injected. Water nozzle location also has significant effect for NO formation and fuel ignition. The chemical kinetic model is also sensitive to the excess water. Through this study, a better understanding of the physics and chemical kinetics is obtained, this will enhance the STC design process.

Fan, L.↗

Evaluation of Water Injection Effect on NO(x) Formation for a Staged Gas Turbine Combustor

NO(x) emission control by water injection on a staged turbine combustor (STC) was modeled using the KIVA-2 code with modification. Water is injected into the rich-burn combustion zone of the combustor by a single nozzle. Parametric study for different water injection patterns was performed. Results show NO(x) emission will decrease after water being injected. Water nozzle location also has significant effect for NO formation and fuel ignition. The chemical kinetic model is also sensitive to the excess water. Through this study, a better understanding of the physics and chemical kinetics is obtained, this will enhance the STC design process.

Fan, L.↗

Theoretical study of production of unique glasses in space

Analytical functional relationships describing homogeneous nucleation and crystallization in various supercooled liquids were developed. The time and temperature dependent relationships of nucleation and crystallization (intrinsic properties) are being used to relate glass forming tendency to extrinsic parameters such as cooling rate through computer simulation. Single oxide systems are being studied initially to aid in developing workable kinetic models and to indicate the primary materials parameters affecting glass formation. The theory and analytical expressions developed for simple systems is then extended to complex oxide systems. A thorough understanding of nucleation and crystallization kinetics of glass forming systems provides a priori knowledge of the ability of a given system to form a glass.

Larsen, D. C.↗