SEARCH · Engineering Papers
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.
Empirically based kinetic model describing synergistic inactivation of dry Bacillus subtilis by combined heat and gamma radiation environment
The BGK kinetic model for a binary mixture of gases was used to predict the backscattering return flux of molecules outgassing from a spherical spacecraft. For simplicity, the outgassing molecules were considered to be a single species, with the ambient atmospheric molecules being the second species in the binary mixture. Two basic collision mechanisms were considered: 1) collisions involving only outgas molecules (self-scattering); and 2) collisions between outgas molecules and molecules in the ambient atmosphere (ambient-scattering). Results were obtained in the form of simple equations relating outgas return flux to parameters based on outgas characteristics, spacecraft size, and the characteristics of the ambient flow.
Experiments are conducted in both rapid compression machine (RCM) and shock tube (ST) to better quantify autoignition behavior (e.g., ignition delay, heat release) and understand heteroatomic effects in heterocyclic compounds, which are important reference components for the combustion of biomass-derived liquid fuels. These tests focus on the nitrogen-containing, five-membered saturated ring, pyrrolidine, at diluted conditions covering pressures of 20 and 50 bar, temperatures of 720–1450 K and a range of stoichiometries (ϕ = 0.5–2). A chemical kinetic model is developed and coupled to an existing combustion kinetics framework describing key nitrogen containing intermediates (e.g. pyrrole, ammonia and NOx). H-abstraction reactions by OH, H, CH 3 and HO 2 , are determined using ab-initio transition state theory methods, while analogies to cyclopentane are adopted for many other reactions, such as ring-opening. Here, the autoignition measurements reveal the lack of negative temperature coefficient (NTC) behavior and low-temperature chemistry for pyrrolidine, as opposed to its saturated hydrocarbon analogue, cyclopentane. Interestingly, at the lowest temperatures (T < 750 K), the reactivity of cyclopentane is greater than pyrrolidine, while at higher temperatures, pyrrolidine becomes more reactive. Agreement between the experimental measurements and the model is good, and it is found that H-abstraction reactions by HO 2 and ensuing chemistry play key roles in controlling the reactivity of this cyclic amine. Most of the fluxes, i.e., >70 %, are predicted to move through 1- or 2-pyrroline (C 4 H 7 N) and then the cyclic C 4 H 6 N radical, at both lower and higher temperatures, to form either CH 2 CHCHCHNH via ring-opening or pyrrole via β-scission. It appears that the ring opens more easily at lower temperature whereas the C–H β-scission dominates at higher temperature and lower pressure, such that the reaction of the fuel radical intermediate carrying an unpaired electron on the nitrogen atom with HO 2 is the next most notable in promoting oxidation. When comparing pyrrolidine and cyclopentane, which exhibits distinct pathways in different temperature regimes, the pyrrolidine pathways and sensitivity analysis align more closely to the high temperature case of cyclopentane where the important role of HO 2 radicals is seen to provide chain branching through HO 2 reaction with the fuel, accompanied by H 2 O 2 formation and decomposition to OH. The formation of 5-membered diene rings and ring opening reactions are also found to be highly relevant. Of particular note, it is found that there is little influence of small molecule nitrogen-chemistry, e.g., NH 2 , HCN, NO/NO 2 on the reactivity of the pyrrolidine mixtures investigated here where no recirculated combustion gases are included.
The Mini-Magnetospheric Plasma Propulsion (M2P2), originally proposed by Winglee et al. [2000], is based on the two-fluid plasma model and requires a 15-km frontal standoff distance (or 20-km cross-sectional diameter) in order for the magnetic bubble to absorb sufficient momentum from the SW to accelerate a spacecraft to the unprecedented speeds of 50-80 km/s after an acceleration period of about three months. Winglee et al. [2000] derived the above size requirement based on an extrapolation of their simulated results in which a system much smaller than a M2P2 was used (p. 21,074 of their study). We submit, however, that a fluid model has no validity for such a small scale size-even in the region near the plasma source! It is assumed in the MHD fluid model, normally applied to the magnetosphere, that the characteristic scale-size is much greater than the Larmor radius and ion skin depth of the SW. In the case of the M2P2, however, the size of the magnetic bubble is actually less than or, at best, comparable to, the scale of these characteristic parameters and, therefore, a kinetic approach, which addresses the smallscale physical mechanisms involved, must be used. A fully three-dimensional version of the hybrid code is used in our M2P2 (Plasma Sails) studies was originally developed by Delamere et al. [1999]. The M2P2 plasma sail is an excellent application for this hybrid code. The primary advantage of this code is the seamless interface between fluid and kinetic descriptions of the ion populations. A kinetic description is not necessary for the dense inner regions of the magnetic bubble and tremendous computational savings can be realized by treating this dense, magnetized ion population with the fluid description. It is essential, however, that the outer bubble regions be treated kinetically as well as the SW protons. Comparison of full size M2P2 simulation based on 3D MHD and kinetic models show that kinetic treatment introduces much more asymmetry to the considering problem and the possibility of kinetic instabilities development.
In recent years, the space travel industry has grown exponentially, resulting in the need for a low-cost, efficient rocket engine fuel. As such, there has been renewed interest in utilizing liquid natural gas (LNG) as it is less likely to soot than kerosene-based fuels and has widespread availability and low cost of use as compared to the traditional RP-X and even liquid methane (LCH4) fuels due to the reduced need for refinement. However, current literature does not fully cover varying blends of natural gas and the effects that impurities have on natural gas/methane fuels, with most studies being confined to low-pressure applications. Therefore, new experimental ignition delay time measurements at rocket engine relevant pressures are essential to ensure current chemical kinetic models capture the ignition behavior at these elevated conditions and to determine the optimal purity levels for a reliable, cost-efficient, aerospace-grade LCH4/LNG fuel. The current work explores various blends of natural gas/methane fuel with various impurities, including higher hydrocarbons as well as nitrogen-carriers. A shock tube study was carried out to study the ignition delay times and carbon monoxide time histories of these natural gas blends utilizing chemiluminescence and laser absorption spectroscopy at 20, 50 and 100 atm over temperatures ranging from 1400-1700K and at an equivalence ratio of 1. The experimental data was then compared to two different chemical kinetic mechanisms: the industry-standard GRI 3.0 and the in-house UCF 2022. The data is used to improve chemical kinetic mechanisms and modeling of LNG fuels for rocket engine applications.
Heterogeneous olefin metathesis catalysts exhibit low active site densities and unpredictable kinetics due to dynamic active site formation and decay processes. Here, in this study, we establish a quantitative framework that captures active site generation, renewal, and decay in olefin metathesis over silica-supported molybdenum oxide catalysts, enabling a mechanistic explanation of catalytic behavior and strategies to achieve high, stable activity. Steady-state active site titrations reveal that 2,3-dimethyl-butene isomers (4MEs) cofeeding increases active site density by up to 4.3-fold, directly correlating with enhanced metathesis rates. Spectroscopic studies demonstrate that 4MEs facilitate Mo(VI) reduction to Mo(IV) and interact strongly with surface Si–OH groups, generating labile protons that promote active site formation via a 1,2-proton shift mechanism. Kinetic modeling indicates that ethylene acts as a decay promoter, shifting kinetic control away from the Chauvin cycle and suppressing metathesis activity. Comparative studies on catalysts with varying Mo loading reveal that promotion is most effective for dispersed molybdate species, with a decline at higher Mo loadings. These findings provide a unified mechanistic framework for heterogeneous olefin metathesis, offering new strategies to enhance active site accessibility, mitigate deactivation, and optimize catalyst design.
Dark particles with lepton-flavor-violating couplings to the tau lepton can induce rare neutrinoless $τ$ decays with large final state multiplicities. We study models where transitions of the type $τ^\pm\to \ell^\pm\,ϕ$, with $ϕ$ a light new particle, initiate a chain of decays in the dark sector that terminate with decays into electrons, muons, or pions. These decay cascades appear as rare five or even seven-body $τ$ decays with multiple reconstructable resonances. We survey several representative models: kinetically mixed dark photon, gauged $L_i-L_j$ models, and other more exotic charge assignments such as chiral $U(1)'$ extensions of the Standard Model. The main new ingredient is the possibility of flavor violation at very high scales. In these models, a number of channels that have not yet been searched for experimentally, such as $τ\to 5μ$, $τ\to 3μ\,2e$, $τ\to μ\,4e$, and hadronic channels like $τ\to μ\,4π$, typically dominate over the previously-considered signatures such as $τ\to 3μ$. While some of the models, such as the gauged $L_i-L_j$ ones, also contain more challenging channels with missing energy due to decays to neutrinos, they can still be searched for via fully visible channels.
Abstract. Ice growth from vapor deposition is an important process for the evolution of cirrus clouds, but the physics of depositional ice growth at the low temperatures (<235 K) characteristic of the upper troposphere and lower stratosphere is not well understood. Surface attachment kinetics, generally parameterized as a deposition coefficient αD, control ice crystal habit and also may limit growth rates in certain cases, but significant discrepancies between experimental measurements have not been satisfactorily explained. Experiments on single ice crystals have previously indicated the deposition coefficient is a function of temperature and supersaturation, consistent with growth mechanisms controlled by the crystal's surface characteristics. Here we use observations from cloud chamber experiments in the Aerosol Interactions and Dynamics in the Atmosphere (AIDA) aerosol and cloud chamber to evaluate surface kinetic models in realistic cirrus conditions. These experiments have rapidly changing temperature, pressure, and ice supersaturation such that depositional ice growth may evolve from diffusion limited to surface kinetics limited over the course of a single experiment. In Part 1, we describe the adaptation of a Lagrangian parcel model with the Diffusion Surface Kinetics Ice Crystal Evolution (DiSKICE) model (Zhang and Harrington, 2014) to the AIDA chamber experiments. We compare the observed ice water content and saturation ratios to that derived under varying assumptions for ice surface growth mechanisms for experiments simulating ice clouds between 180 and 235 K and pressures between 150 and 300 hPa. We found that both heterogeneous and homogeneous nucleation experiments at higher temperatures (>205 K) could generally be modeled consistently with either a constant deposition coefficient or the DiSKICE model assuming growth on isometric crystals via abundant surface dislocations. Lower-temperature experiments showed more significant deviations from any depositional growth model, with different ice growth rates for heterogeneous and homogeneous nucleation experiments.
The prevalent use of organic materials in manufacturing is a fire safety concern, and motivates the need for predictive thermal decomposition models. A critical component of predictive modeling is numerical inference of kinetic parameters from bench scale data. Currently, an active area of computational pyrolysis research focuses on identifying efficient, robust methods for optimization. This paper demonstrates that kinetic parameter calibration problems can successfully be solved using classical gradient-based optimization. We explore calibration examples that exhibit characteristics of concern: high nonlinearity, high dimensionality, complicated schemes, overlapping reactions, noisy data, and poor initial guesses. The examples demonstrate that a simple, non-invasive change to the problem formulation can simultaneously avoid local minima, avoid computation of derivative matrices, achieve a computational efficiency speedup of 10x, and make optimization robust to perturbations of parameter components. Techniques from the mathematical optimization and inverse problem communities are employed. By re-examining gradient-based algorithms, we highlight opportunities to develop kinetic parameter calibration methods that should outperform current methods.
Hydride vapor phase epitaxy (HVPE) is a promising technology that can aid in the cost reduction of III-V materials and devices manufacturing, particularly high-efficiency solar cells for space and terrestrial applications. However, recent demonstrations of ultrafast growth rates (~500 µm/h) via uncracked hydrides are not well described by present models for the growth. Therefore, it is necessary to understand the kinetics of the growth process and its coupling with transport phenomena, so as to enable fast and uniform epitaxial growth. In this work, we derive a kinetic model using experimental data and integrate it into a computational fluid dynamics simulation of an HVPE growth reactor. We also modify an existing hydride cracking model that we validate against numerical simulations and experimental data. Here, we show that the developed growth model and the improved cracking model are able to reproduce experimental growth measurements of GaAs in an existing HVPE system.
The necessity to reveal the formation mechanism of not only polycyclic aromatic hydrocarbons (PAHs) but also oxygenated PAHs (OPAHs) during combustion is increasing. Although many studies on PAHs have been conducted, fundamental studies investigating OPAH formation are still limited. Phenol, benzofuran, and dibenzofuran were selected as OPAHs in this study. We experimentally investigated fuel-rich oxidation in a flow reactor at atmospheric pressure, mean gas temperatures from 1050 to 1350 K, residence times from 0.2 to 1.5 s, and equivalence ratios from 3.0 to 12.0. Ethylene, toluene, and n-decane were used as fuels. Three kinds of OPAHs as well as 23 kinds of PAHs including monocyclic structures were quantitatively measured through a direct sampling using gas chromatography mass spectrometry. The results showed that concentrations of PAHs and OPAHs were strongly affected by the fuel type, with increasing concentrations in the following order: ethylene < n-decane < toluene. A chemical kinetic model for OPAH formation was developed based on a recent PAH growth model. The model showed that the predicted concentrations of OPAHs and PAHs were in reasonable agreement with the measured data in this study. Some modifications were made to the previous model based on recent literature studies and on the comparison of the simulated and measured results. The effect of the fuel type on the formation of PAHs and OPAHs was investigated through kinetic analysis using the model to discuss their reaction pathways.
Fuel-rich oxidation of three typical gasoline surrogate components, toluene, isooctane, and n-heptane, was investigated in an atmospheric-pressure flow reactor at mean gas temperatures from 1050 to 1350 K, equivalence ratio of 9.0, and residence times of 0.45 and 1.2 s. Not only polycyclic aromatic hydrocarbons (PAHs) up to 3 ring structure but also small intermediate products from C 1 to C 5 were quantified by a gas chromatograph mass spectrometry coupled with photon ionization and gas chromatograph with flame ionization detector, respectively. The kinetic model recently developed by Lawrence Livermore National Laboratory was revised to reflect the results of many recent investigations. Basically, the updated model could satisfactorily reproduce the experimental mole fractions of many species. The experimental and simulated results showed that PAH mole fractions produced were in the order of toluene, isooctane, and n-heptane. The kinetic analysis using the model was carried out to explore PAH formation pathways, especially focusing on naphthalene, acenaphthylene, and phenanthrene. Through rate of production analysis, it was found that the main formation pathways of many PAHs were affected by the fuels. Although resonantly stabilized radicals, such as benzyl and fulvenallenyl radicals, played a crucial role in the formation pathways of many PAHs in every fuel, they were produced through hydrogen elimination of toluene in toluene fuel, while they were formed from small intermediate products in isooctane and n-heptane fuels. Sensitivity analysis revealed the difference and similarity of the reactions with large positive coefficients according to the fuels studied here. In conclusion, the molecular growth reactions of aromatic species were influential in PAH production in every fuel, whereas the ring formation reactions from small species and the reactions involving toluene had large positive sensitivity coefficients in isooctane/n-heptane and toluene fuels, respectively.
The proposed project was aimed at assisting in the development of the Argon Power Cycle (APC), a zero emission, flexible, and efficient power generation technology. In the effort of bringing the APC to the market, NTS proposed in the scope of this project several partial but key activities to advance the development of the APC. These tasks encompassed 1) the analysis of the market, developing a techno-economic model with integrated market inputs, 2) further optimization of our chemical kinetics models for more accurate simulation of the process and 3) scale up of our original prototype to a larger, commercially viable platform. NTS expanded on this scope of work to secure additional R&D funding in amounts north of $7M, carry out customer discovery with the sponsorship of the NSF and grow the team’s engineering talent.
The chondrites are aggregates of components (e.g. chondrules, chondrule rims and matrix) that formed in the nebula but, at present, there is no consensus on how any of these components formed or whether their formation produced or post dated the chemical fractionations between the chondrites. Chondrites are, at present, the most primitive Solar System objects available for laboratory study and the conditions under which their principle components formed would provide the most direct constraints for models of nebula formation and evolution. The conditions under which chondrules formed is of particular importance because, if their relative abundance in chondrites approximates that in the nebula, they are the products of one of the most energetic and pervasive processes that operated in the early Solar System. The goal of this proposal was to combine theoretical modeling with a comprehensive study of the elemental and isotopic compositions of the major components in unequilibrated ordinary chondrites (UOCs), with the aim of determining the conditions in the nebula at the time of their formation. The isotopes of volatile and moderately volatile elements should be particularly revealing of conditions during chondrule formation, as evaporation under most conditions would lead to isotopic mass fractionation. Modeling of chondrule and matrix formation requires the development of a kinetic model of evaporation and condensation, and calibration of this model against experiments.
The chondrites are aggregates of components (e.g. chondrules, chondrule rims and matrix) that formed in the nebula but, at present, there is no consensus on how any of these components formed or whether their formation produced or post dated the chemical fractionations between the chondrites. Chondrites are, at present, the most primitive Solar System objects available for laboratory study and the conditions under which their principle components formed would provide the most direct constraints for models of nebula formation and evolution. The conditions under which chondrules formed is of particular importance because, if their relative abundance in chondrites approximates that in the nebula, they are the products of one of the most energetic and pervasive processes that operated in the early Solar System. The goal of this proposal was to combine theoretical modeling with a comprehensive study of the elemental and isotopic compositions of the major components in unequilibrated ordinary chondrites (UOCs), with the aim of determining the conditions in the nebula at the time of their formation. The isotopes of volatile and moderately volatile elements should be particularly revealing of conditions during chondrule formation, as evaporation under most conditions would lead to isotopic mass fractionation. Modeling of chondrule and matrix formation requires the development of a kinetic model of evaporation and condensation, and calibration of this model against experiments. Cosmic spherules present an opportunity to test our evaporation models under flash heating conditions that would be difficult to simulate experimentally. However, there is surprisingly little known about the isotopic compositions of silicate cosmic spherules, and a number of questions need to be addressed. Is the range of compositions they exhibit due to evaporation? If they are, are the relative volatilities consistent with the models/experiments and are the isotopic fractionations consistent with Rayleigh conditions? For instance, do the alkalis and S evaporate prior to significant melting so that conditions did not meet the Rayleigh criteria of rapid diffusion? If so, their isotopic fractionation might be considerably suppressed. Could this mechanism of K loss apply to chondrule formation? The Fe isotopic fractionation during evaporation of silicates has not been measured, so cosmic spherules might provide a clue to whether FeO diffusion is fast enough to maintain Rayleigh conditions during evaporation. And so on.
The risks of late effects from galactic cosmic rays (GCR) and solar particle events (SPE) are potentially a limitation to long-term space travel. The late effects of highest concern have significant lethality including cancer, effects to the central nervous system (CNS), and circulatory diseases (CD). For cancer and CD the use of age and gender specific models with uncertainty assessments based on human epidemiology data for low LET radiation combined with relative biological effectiveness factors (RBEs) and dose- and dose-rate reduction effectiveness factors (DDREF) to extrapolate these results to space radiation exposures is considered the current "state-of-the-art". The revised NASA Space Risk Model (NSRM-2014) is based on recent radio-epidemiology data for cancer and CD, however a key feature of the NSRM-2014 is the formulation of particle fluence and track structure based radiation quality factors for solid cancer and leukemia risk estimates, which are distinct from the ICRP quality factors, and shown to lead to smaller uncertainties in risk estimates. Many persons exposed to radiation on earth as well as astronauts are life-time never-smokers, which is estimated to significantly modify radiation cancer and CD risk estimates. A key feature of the NASA radiation protection model is the classification of radiation workers by smoking history in setting dose limits. Possible qualitative differences between GCR and low LET radiation increase uncertainties and are not included in previous risk estimates. Two important qualitative differences are emerging from research studies. The first is the increased lethality of tumors observed in animal models compared to low LET radiation or background tumors. The second are Non- Targeted Effects (NTE), which include bystander effects and genomic instability, which has been observed in cell and animal models of cancer risks. NTE's could lead to significant changes in RBE and DDREF estimates for GCR particles, and the potential effectiveness of radiation mitigator's. The NSRM- 2014 approaches to model radiation quality dependent lethality and NTE's will be described. CNS effects include both early changes that may occur during long space missions and late effects such as Alzheimer's disease (AD). AD effects 50% of the population above age 80-yr, is a degenerative disease that worsens with time after initial onset leading to death, and has no known cure. AD is difficult to detect at early stages and the small number of low LET epidemiology studies undertaken have not identified an association with low dose radiation. However experimental studies in mice suggest GCR may lead to early onset AD. We discuss modeling approaches to consider mechanisms whereby radiation would lead to earlier onset of occurrence of AD. Biomarkers of AD include amyloid beta (A(Beta)) plaques, and neurofibrillary tangles (NFT) made up of aggregates of the hyperphosphorylated form of the micro-tubule associated, tau protein. Related markers include synaptic degeneration, dentritic spine loss, and neuronal cell loss through apoptosis. Radiation may affect these processes by causing oxidative stress, aberrant signaling following DNA damage, and chronic neuroinflammation. Cell types to be considered in multi-scale models are neurons, astrocytes, and microglia. We developed biochemical and cell kinetics models of DNA damage signaling related to glycogen synthase kinase-3(Beta) (GSK3(Beta)) and neuroinflammation, and considered multi-scale modeling approaches to develop computer simulations of cell interactions and their relationships to A(Beta) plaques and NFTs. Comparison of model results to experimental data for the age specific development of A(Beta) plaques in transgenic mice will be discussed.
Data-driven approaches to construct reduced chemical kinetic models, that rely heavily on thermo-chemical datasets with full chemical kinetics, have been gaining popularity. Datasets from direct numerical simulations (DNS) under three-dimensional (3-D) realistic turbulent flow conditions are desirable but limited to carefully designed parametric conditions due to the computational cost. Constructing datasets from a large ensemble of zero-dimensional stirred reactors like perfectly stirred reactor (PSR) and partially stirred reactor (PaSR) is a computationally efficient solution to consider the turbulence-chemistry interactions and cover a broad range of parametric conditions. In this paper, we derive reduced chemistry models from solutions of a large number of PSR and PaSR reactors using autoencoder (AE) neural networks and principal component analysis (PCA), and conduct a priori examination of the reduced models in three temporally evolving 3-D DNS jet flames featuring local extinction and re-ignition. The results show that the reduced models derived from PaSR datasets, i.e., AE-PaSR and PCA-PaSR, generally show significant improvement over the ones derived from PSR datasets. Among all the reduced models, AE-PaSR shows the best agreement with DNS results on the reconstruction accuracy and the representation of temporally evolving local extinction and re-ignition events.
Abstract The Rayleigh–Taylor (RT) instability is ubiquitously observed, yet has traditionally been studied using ideal fluid models. Collisionality can vary strongly across the fluid interface, and previous work demonstrates the necessity of kinetic models to completely capture dynamics in certain collisional regimes. Where previous kinetic simulations used spatially and temporally constant collision frequency, this work presents five-dimensional (two spatial, three velocity dimensions) continuum-kinetic simulations of the RT instability using a more realistic spatially varying collision frequency. Three cases of collisional variation are explored for two Atwood numbers: low to intermediate, intermediate to high, and low to high. The low-to-intermediate case exhibits no RT instability growth, while the intermediate-to-high case is similar to a fluid-limit kinetic case with interface widening biased toward the lower-collisionality region. A novel contribution of this work is the low-to-high collisionality case that shows significantly altered instability growth through an upward movement of the interface and damped spike growth due to increased free-streaming particle diffusion in the lower region. Contributions to the energy flux from the non-Maxwellian portions of the distribution function are not accessible to fluid models and are greatest in magnitude in the spike and regions of low collisionality. Increasing the Atwood number results in greater RT instability growth and reduced upward interface movement. Deviation of the distribution function from Maxwellian is inversely proportional to collision frequency and concentrated around the fluid interface. The linear phase of RT instability growth is well described by theoretical linear growth rates accounting for viscosity and diffusion.