Search NASA⌕ Search

SEARCH · Search NASA

Results for “Kinetic Modelling”

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

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

At least 757 records · Page 42

A purely kinetic description of the evaporation of water droplets

The process of water evaporation, although deeply studied, does not enjoy a kinetic description that captures known physics and can be integrated with other detailed processes such as drying of catalytic membranes embedded in vapor-fed devices and chemical reactions in aerosol whose volumes are changing dynamically. In this work, we present a simple, three-step kinetic model for water evaporation that is based on theory and validated by using well-established thermodynamic models of droplet size as a function of time, temperature, and relative humidity as well as data from time-resolved measurements of evaporating droplet size. Additionally, the kinetic mechanism for evaporation is a combination of two limiting processes occurring in the highly dynamic liquid-vapor interfacial region: direct first order desorption of a single water molecule and desorption resulting from a local fluctuation, described using third order kinetics. The model reproduces data over a range of relative humidities and temperatures only if the interface that separates bulk water from gas phase water has a finite width, consistent with previous experimental and theoretical studies. The influence of droplet cooling during rapid evaporation on the kinetics is discussed; discrepancies between the various models point to the need for additional experimental data to identify their origin.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Hybrid Simulation of Low Temperature Plasmas: A Brief Tutorial

Modeling and simulation of low temperature plasmas is widely accepted as a useful tool from both physics and plasma reactor design points of view. Of the modeling approaches, hybrid models strive to combine the predictability of kinetic models with the lower computational burden of fluid models. This work presents a brief tutorial of hybrid plasma modeling and representative simulation results using hybrid models. Whenever possible, these results are compared to kinetic and/or fluid simulations as well as experimental data.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Laboratory Studies of the Role of Amines in Particle Formation, Growth and Climate

Organosulfur compounds have a variety of sources, particularly biological processes in the oceans. Their oxidation in air forms sulfur dioxide, which is further oxidized to sulfuric acid, as well as methanesulfonic acid (MSA). While sulfuric acid is a well known precursor to particles in air, MSA had not been regarded as a source of new particle formation. Laboratory studies were carried out under this project that showed MSA forms new particles quite efficiently in the presence of amines and water vapor. The data could be reproduced with a relatively simple kinetics model representing cluster formation and growth, which is promising for representing this chemistry in global climate models. The initial steps in the kinetics scheme are based on quantum chemical calculations of likely clusters. The organosulfur chemistry was introduced into an atmospheric model for southern California and used to predict the impact of going to a fossil-fuel free world in which anthropogenic emissions of SO 2 are removed, but the natural processes remain.

54 ENVIRONMENTAL SCIENCES↗

Dehydroxylation kinetics of kaolinite and montmorillonite examined using isoconversional methods

The use of calcined clays as supplementary cementitious materials (SCMs) in concrete is a promising strategy towards decarbonizing the cement and concrete industry. This is especially relevant considering the ever-increasing demand for concrete. Comprehensive understanding of the kinetics of calcination is essential towards maximizing the potential reactivity of clay minerals while ensuring energy efficiency. In this study, the kinetics of the dehydroxylation of kaolinite and montmorillonite are investigated under non-isothermal conditions at constant heating rate. Activation energies ( E a ) are determined via Friedman differential and advanced Vyazovkin incremental methods over the isoconversional range; these are devoid of computational approximations, thus allowing kinetic analysis without assuming a specific reaction model. Kinetic equations—in the differential form as well as a combination of differential and integral forms are compared against the experimentally determined reaction models to identify the most probable dehydroxylation mechanism for kaolinite and montmorillonite. A reaction order mechanism is established for dehydroxylation of kaolinite, while montmorillonite is noted to undergo dehydroxylation via a single-step reversible diffusion-controlled process. Kinetic triplet—comprising activation energy, reaction model and pre-exponential factor—is used to predict isothermal calcination conditions, which is further verified using analytical techniques. Heat release rates of clay-portlandite blends from isothermal calorimetry are used within a thermodynamic framework to quantify reactivity of the calcined clays. Here, the study demonstrates a general approach based on isoconversional methods to predict calcination conditions for different clays that can be used in efficient and optimized production of blended cements or SCMs.

36 MATERIALS SCIENCE↗

Performance of Polymer Electrolyte Fuel Cell Electrodes with Atomically Dispersed (AD) Fe-C-N ORR Catalyst

Activity of (AD)Fe-N-C catalyst with low Fe content is investigated in differential cells prepared by hot pressing anode gas diffusion electrodes (0.2 mgPt/cm2) to N211 membrane and brush painting the cathode catalyst ink. Polarization curves obtained in H2/O2 show double Tafel slopes which, in conjunction with the redox potential observed in cyclic voltammetry traces, forms the basis for a distributed ORR (oxygen reduction reaction) kinetic model with potential-dependent available sites. Application of this model to polarization data in H2/air provides the basis for formulating an oxygen transport model and leads to 7.7 ± 0.6 mA/cm2 catalyst activity at 0.9 V, 31.5-34.3 mA/cm2 cell performance at 0.8 V, and 0.8-2 s/cm O2 transport resistance. A coupled kinetic, O2 transport and proton transport illustrates the projected improvements needed in catalyst activity and electrode structure to approach the automotive target of 1000 mW/cm2 stack power density while meeting the 1.45 kW/oC heat rejection constraint at 1.5 cathode stoichiometry. The model indicates that we need twelve-fold higher mass activity for reducing the kinetic losses, doubling of active site density and an engineered electrode structure for 50% lower proton transport resistance, as well as a 50% reduction in electrode thickness to limit the O2 transport losses

ACTIVE-SITES↗

High-order continuum kinetic method for modeling plasma dynamics in phase space

Continuum methods offer a high-fidelity means of simulating plasma kinetics. While computationally intensive, these methods are advantageous because they can be cast in conservation-law form, are not susceptible to noise, and can be implemented using high-order numerical methods. Advances in continuum method capabilities for modeling kinetic phenomena in plasmas require the development of validation tools in higher dimensional phase space and an ability to handle non-cartesian geometries. To that end, a new benchmark for validating Vlasov-Poisson simulations in 3D (x,v x ,v y ) is presented. The benchmark is based on the Dory-Guest-Harris instability and is successfully used to validate a continuum finite volume algorithm. To address challenges associated with non-cartesian geometries, unique features of cylindrical phase space coordinates are described. Preliminary results of continuum kinetic simulations in 4D (r,z,v r ,v z ) phase space are presented.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Flame kinetics at scramjet-engine-relevant conditions: Role of prompt dissociation of weakly-bound radicals

Combustion in high-speed ram-based propulsion engines occurs under distinct thermodynamic conditions of high reactant temperatures (greater than 1000 K) and relatively low pressures (<5 atm). There is a lack of fundamental flame measurements at such conditions that result in adiabatic flame temperatures (T ad ) exceeding 2500 K. In this work, we have measured laminar flame speeds of oxygen-enriched CH 4 /oxidizer mixtures at sub-atmospheric conditions to probe kinetics at high T ad using the isobaric spherically expanding flame approach. Simulations with recent kinetic models revealed increasing differences between data and model predictions with increasing T ad , reaching up to 25 %. Kinetic analyses reveal that at the thermodynamic conditions in these O 2 -enriched flames, i.e., lower pressures and higher T ad , the effects of HCO prompt dissociation are accentuated. In addition to HCO, the prompt dissociations of CH 2 OH and C 2 H 5 are also considered. Here, the prompt dissociations of all three radicals were evaluated and their effects considered in flame speed simulations. Reaction path analysis for the present flames revealed that approximately half of the reaction flux for HCO formation undergoes prompt dissociation to H + CO. Furthermore, these analyses also revealed that the pathways and sensitive reactions are similar between oxygen-enriched fuel/oxidizer mixtures and preheated fuel/air mixtures, if both have similar T ad . Thus, flames of oxygen-enriched mixtures could be a surrogate to probe the flame chemistry of highly preheated mixtures at relatively low pressures that are often encountered in ram-based propulsion engine combustors.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Three-dimensional, off-lattice Monte-Carlo kinetics simulations of interstellar grain chemistry and ice structure

The first off-lattice Monte Carlo kinetics model of interstellar dust grain surface chemistry is presented. The positions of all surface particles are determined explicitly, according to the local potential minima resulting from the pair-wise interactions of contiguous atoms and molecules, rather than by a pre-defined lattice structure. The model is capable of simulating chemical kinetics on any arbitrary dust grain morphology, as determined by the user-defined positions of each individual dust grain atom. A simple method is devised for the determination of the most likely diffusion pathways and their associated energy barriers for surface species. The model is applied to a small, idealized dust grain, adopting various gas densities and using a small chemical network. Hydrogen and oxygen atoms accrete onto the grain to produce H{sub 2}O, H{sub 2}, O{sub 2}, and H{sub 2}O{sub 2}. The off-lattice method allows the ice structure to evolve freely; the ice mantle porosity is found to be dependent on the gas density, which controls the accretion rate. A gas density of 2 × 10{sup 4} cm{sup –3}, appropriate for dark interstellar clouds, is found to produce a fairly smooth and non-porous ice mantle. At all densities, H{sub 2} molecules formed on the grains collect within the crevices that divide nodules of ice and within micropores (whose extreme inward curvature produces strong local potential minima). The larger pores produced in the high-density models are not typically filled with H{sub 2}. Direct deposition of water molecules onto the grain indicates that amorphous ices formed in this way may be significantly more porous than interstellar ices that are formed by surface chemistry.

79 ASTRONOMY AND ASTROPHYSICS↗

Development of the FHR advanced natural circulation analysis code and application to FHR safety analysis

The University of California, Berkeley (UCB) is performing thermal hydraulics safety analysis to develop the technical basis for design and licensing of fluoride-salt-cooled, high-temperature reactors (FHRs). FHR designs investigated by UCB use natural circulation for emergency, passive decay heat removal when normal decay heat removal systems fail. The FHR advanced natural circulation analysis (FANCY) code has been developed for assessment of passive decay heat removal capability and safety analysis of these innovative system designs. The FANCY code uses a one-dimensional, semi-implicit scheme to solve for pressure-linked mass, momentum and energy conservation equations. Graph theory is used to automatically generate a staggered mesh for complicated pipe network systems. Heat structure models have been implemented for three types of boundary conditions (Dirichlet, Neumann and Robin boundary conditions). Heat structures can be composed of several layers of different materials, and are used for simulation of heat structure temperature distribution and heat transfer rate. Control models are used to simulate sequences of events or trips of safety systems. A proportional-integral controller is also used to automatically make thermal hydraulic systems reach desired steady state conditions. A point kinetics model is used to model reactor kinetics behavior with temperature reactivity feedback. The underlying large sparse linear systems in these models are efficiently solved by using direct and iterative solvers provided by the SuperLU code on high performance machines. Input interfaces are designed to increase the flexibility of simulation for complicated thermal hydraulic systems. In conclusion, this paper mainly focuses on the methodology used to develop the FANCY code, and safety analysis of the Mark 1 pebble-bed FHR under development at UCB is performed.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Experimental and modeling study of water time histories during H 2 S-N 2 O combustion in a shock tube

Here, with the interest in directly burning sour gas in gas turbines, and the fact that even small amounts of H 2 S or its combustion products can alter combustion characteristics, many research studies have been performed to better understand the combustion chemistry of H 2 S. In the present study, the water formation was followed by laser absorption with N 2 O as an oxidant, instead of O 2 . Nitrous Oxide being essentially consumed via N 2 O (+M) ⇌ N 2 + O (+M), the water formation via the H 2 S + O route can then be probed to further validate the models. Three H 2 S/N 2 O mixtures diluted in 98% Ar were studied to cover the following range of equivalence ratios: 0.5, 1.0, and 2.0, over a wide range of temperatures (1580–1940 K) around atmospheric pressure. A chemical kinetic model was then developed first by validating the base N 2 O kinetics mechanism, then by investigating the experimental results presented herein. The N 2 O kinetics mechanism was updated with recent work on the low- and high-pressure limits for N 2 O decomposition (N 2 O (+M) ⇌ N 2 + O (+M)) as well as a review of rate coefficients for N 2 O + H ⇌ N 2 + OH from the literature. Good agreement is shown for H 2 /N 2 O mixtures. Updates were then made to the H 2 S kinetics mechanism, specifically, an update from the literature on SO 2 + H (+M) ⇌ SO + OH (+M) and an adjustment to SO + SH ⇌ S 2 + OH. Additionally, reactions between SH and N 2 O were determined using W1BD data and transition state theory which required the addition of an NNS sub-mechanism. With these updates, the mechanism provides good agreement with the H 2 S/N 2 O experiments.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

A theoretical explanation of the effect of oxygen poisoning on industrial Haber-Bosch catalysts

The Haber-Bosch process has been investigated extensively, however a low-temperature, low-pressure process remains elusive. As has been demonstrated many times, this stems in part from the difficulty of breaking the N-N triple bond. In this work, we highlight an additional reason for the lack of a low-temperature ammonia synthesis process: the effect of oxygen poisoning at low temperature. Using density functional theory (DFT), we have created a new model for the active site of industrial Haber-Bosch catalysts which explicitly includes the potassium promoter. Furthermore, we present a new micro-kinetic model for ammonia synthesis that includes the effect of oxygen poisoning due to trace water content in the input gas stream. Our model agrees well with previous experiments and shows that devising a strategy to avoid oxygen poisoning is crucial to creating a low-temperature Haber-Bosch process. Additionally, the model suggests that using a weaker-binding catalyst is one way to avoid oxygen poisoning if it is impractical to remove all water from the reactor.

08 HYDROGEN↗

Metabolic model guided strain design of cyanobacteria

Cyanobacteria are oxygenic photoautotrophs that serve as potential platforms for the production of biochemicals from cheap and renewable raw materials – sunlight, water, and carbon dioxide. Systems level analysis of the metabolic network of these organisms could enable the successful engineering of these organisms for the enhanced production of target chemicals. Metabolic modeling techniques including both stoichiometric and kinetic modeling with a genome-wide coverage enable a global assessment of metabolic capabilities. Recent studies guided by such modeling techniques have engineered strains for the enhanced production of valuable chemicals such as ethanol, n-butanol, 1,3-propanediol, glycerol, limonene, and isoprene from CO 2 .

59 BASIC BIOLOGICAL SCIENCES↗

Ultra-low concentration terbium (Tb) adsorption on garlic peels biosorbent and its application for $\mathrm{Nd}$-$\mathrm{Fe}$-$\mathrm{B}$ scraps recovery

Motivated by the strategic value of middle-heavy rare earth elements (MHREEs), we proposed a bio-adsorption method to recover ultra-low concentration terbium (Tb(Ⅲ)) from industrial wastewater that originated from the Nd-Fe-B scraps recovery process. It could be found that Tb(Ⅲ) ions as low as 5.34 ng·mL -1 (ppb) could be adsorbed onto Ca(Ⅱ)-modified garlic peels (Ca-GP) within 10 min with the adsorption efficiency of 99.2% at pH 3.5. The adsorption behavior of Tb(Ⅲ) onto Ca-GP conformed the Langmuir isotherm model and pseudo-second-order kinetic model. Attenuated total reflection Fourier Transform Infrared Spectroscopy (ATR-FTIR), X-ray photoelectron spectroscopy (XPS) and density-functional theory (DFT) calculations results showed that Tb(Ⅲ) ions could be ion-exchanged with cations in Ca-GP such as Ca(Ⅱ) and -COOH. In application, column experiment results showed that the maximum bed adsorption capacity for Tb, praseodymium (Pr), neodymium (Nd), dysprosium (Dy), and Total rare earth elements (REEs) ions calculated by the Thomas model was 0.06, 3.50, 9.65, 7.34, and 27.37 µg·g -1 , respectively, with the initial concentration of 0.37 ng·mL -1 Tb, 20 ng·mL -1 Pr, 44 ng·mL -1 Nd, 67 ng·mL -1 Dy and 170 ng·mL -1 total REEs ions of actual solutions (below the solubility of rare-earth hydroxide). Here this manuscript thus provided a novel cost-effective and efficient approach to the enrichment and recovery of MHREEs from ultra-low concentration REEs ions-containing solutions.

54 ENVIRONMENTAL SCIENCES↗

Controlled hydrogenation of a biomass-derived platform chemical formed by aldol-condensation of 5-hydroxymethyl furfural (HMF) and acetone over Ru, Pd, and Cu catalysts

Her we report that we studied the hydrogenation at temperatures from 313–393 K of a biomass-derived platform molecule, 5-hydroxymethyl furfural (HMF)-acetone-HMF (HAH) over Pd, Ru, and Cu based catalysts. HAH was selectively hydrogenated to produce partially-hydrogenated monomers (PHAH) over Cu and Ru catalysts and to fully-hydrogenated HAH monomers (FHAH) over the Ru catalyst. Pd based catalysts yielded a mixture of partially and fully hydrogenated monomers. Lumped reaction kinetics models were employed to quantify the kinetic behavior for hydrogenation over Ru, Cu, and Pd catalysts. The 5-step pathway exhibited over Pd and Ru catalysts consists of both series and parallel reaction steps, where HAH is both converted to fully hydrogenated products sequentially via series reactions of partially hydrogenated intermediates, as well as converted directly in parallel reactions to form the fully hydrogenated products. In contrast, the 3-step pathway over the Cu catalyst consists only of the consecutive reaction steps, where the final product was formed via series reactions of intermediate products. Additionally, reaction over the Cu catalyst did not hydrogenate the furan rings of the HAH molecule and yielded a different final product than those hydrogenation over Pd and Ru catalysts. Batch conditions are determined for each hydrogenated product that give the highest yields in both batch and plug flow reactors.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Radiation from mixed multi-planar wire arrays

The study of radiation from different wire materials in wire array Z-pinch plasma is a very challenging topic because it is almost impossible to separate different plasmas at the stagnation. A new approach is suggested based on planar wire array (PWA) loads to assess this problem. Multi-planar wire arrays are implemented that consist of few planes, each with the same number of wires and masses but from different wire materials, arranged in parallel rows. In particular, the experimental results obtained with triple PWAs (TPWAs) on the UNR Zebra generator are analyzed with Wire Ablation Dynamics Model, non-local thermodynamic equilibrium kinetic model, and 2D radiation magneto-hydrodynamic to illustrate this new approach. In TPWAs, two wire planes were from mid-atomic-number wire material and another plane was from alloyed Al, placed either in the middle or at the edge of the TPWA. Spatial and temporal properties of K-shell Al and L-shell Cu radiations were analyzed and compared from these two configurations of TPWAs. Advantages of the new approach are demonstrated and future work is discussed.

74 ATOMIC AND MOLECULAR PHYSICS↗

Agent-based model forecasts aging of the population of people who inject drugs in metropolitan Chicago and changing prevalence of hepatitis C infections

People who inject drugs (PWID) are at high risk for blood-borne pathogens transmitted during the sharing of contaminated injection equipment, particularly hepatitis C virus (HCV). HCV prevalence is influenced by a complex interplay of drug-use behaviors, social networks, and geography, as well as the availability of interventions, such as needle exchange programs. To adequately address this complexity in HCV epidemic forecasting, we have developed a computational model, the Agent-based Pathogen Kinetics model (APK). APK simulates the PWID population in metropolitan Chicago, including the social interactions that result in HCV infection. We used multiple empirical data sources on Chicago PWID to build a spatial distribution of an in silico PWID population and modeled networks among the PWID by considering the geography of the city and its suburbs. APK was validated against 2012 empirical data (the latest available) and shown to agree with network and epidemiological surveys to within 1%. For the period 2010–2020, APK forecasts a decline in HCV prevalence of 0.8% per year from 44(±2)% to 36(±5)%, although some sub-populations would continue to have relatively high prevalence, including Non-Hispanic Blacks, 48(±5)%. The rate of decline will be lowest in Non-Hispanic Whites and we find, in a reversal of historical trends, that incidence among non-Hispanic Whites would exceed incidence among Non-Hispanic Blacks (0.66 per 100 per years vs 0.17 per 100 person years). APK also forecasts an increase in PWID mean age from 35(±1) to 40(±2) with a corresponding increase from 59(±2)% to 80(±6)% in the proportion of the population >30 years old. Our research highlight the importance of analyzing sub-populations in disease predictions, the utility of computer simulation for analyzing demographic and health trends among PWID and serve as a tool for guiding intervention and prevention strategies in Chicago, and other major cities.

59 BASIC BIOLOGICAL SCIENCES↗