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.
This paper will demonstrate the usefulness of standard non-linear statistical analysis on ICR and SIFT kinetic data. The specific systems used in the demonstration are the isotopic and change transfer reactions in the system of H2O+/D2O, H30+/D2O, and other permutations.
A novel, drag-diffusion Langevin method of hybrid, kinetic-fluid Coulomb scattering in plasmas is presented. Unlike previous methods, the frictional drag is always applied in the simulation frame of reference. The velocity-space diffusion is performed in the stationary-fluid frame of reference when anisotropic, and in the laboratory frame when isotropic. While the general method is mass-ratio independent, we focus on interactions of kinetic-ions and fluid-electrons to show first-order modifications to the electron velocity distribution function that are an important correction for the accurate calculation of electric resistivity. Inclusion of sub-cycling and a limit to the maximum collision frequency is shown to allow for arbitrarily large timesteps without numerical failure. Here the Langevin method is compared with a grid-based binary method and found to require a much less restrictive timestep in cases of ion–electron slowing and temperature equilibration; this finding differs from previous work and is dependent on the mass ratio.
A review of recent theoretical modeling of the effects of turbulence on acoustic wave propagation is presented. The Born approximation limit of the smoothing method of Keller and of the kinetic theory method of Howe are compared with the classical perturbation results of Chernov. Secular growth is observed in the higher order terms of the kinetic theory method of Howe when multiple scaling is ignored. The forward propagation - beam broadening model of Brown and Clifford is discussed as well as the single scatter model of DeLoach and of Rudd. Each model is examined with respect to its prediction of conservation of energy.
The reduction kinetics of hematite iron ore fines to metallic iron by hydrogen using a laboratory fluidized bed reactor were investigated in a temperature range between 873 K and 1073 K, by measuring the weight change of the sample portion during reduction. The fluidization conditions were checked regarding plausibility within the Grace diagram and the measured pressure drop across the material during experiments. The apparent activation energy of the reduction was determined against the degree of reduction and varied along an estimated two-peak curve between 11 and 55 kJ mol{sup −1}. Conventional kinetic analysis for the reduction of FeO to metallic iron, using typical models to describe gas–solid reactions, does not show results with high accuracy. Multistep kinetic analysis, using the Johnson–Mehl–Avrami model, shows that the initial stage of reduction from Fe{sub 2}O{sub 3} to Fe{sub 3}O{sub 4}, and partly to FeO, is controlled by diffusion and chemical reaction, depending on the temperature. Further reduction can be described by a combination of nucleation and chemical reaction, whereby the influence of nucleation increases with an increasing reduction temperature. The results of the kinetical analysis were linked to the shape of the curve from apparent activation energy against the degree of reduction.
An extension of the method of moments is developed for the numerical integration of the kinetic equations of droplet spectra evolution by condensation/evaporation and by coalescence/breakup processes. The number density function n sub k (x,t) in each separate droplet packet between droplet mass grid points (x sub k, x sub k+1) is represented by an expansion in orthogonal polynomials with a given weighting function. In this way droplet number concentrations, liquid water contents and other moments in each droplet packet are conserved and the problem of solving the kinetic equations is replaced by one of solving a set of coupled differential equations for the number density function moments. The method is tested against analytic solutions of the corresponding kinetic equations. Numerical results are obtained for different coalescence/breakup and condensation/evaporation kernels and for different initial droplet spectra. Also droplet mass grid intervals, weighting functions, and time steps are varied.
An efficient, decoupled direct method for calculating the first order sensitivity coefficients of homogeneous, batch combustion kinetic rate equations is presented. In this method the ordinary differential equations for the sensitivity coefficients are solved separately from , but sequentially with, those describing the combustion chemistry. The ordinary differential equations for the thermochemical variables are solved using an efficient, implicit method (LSODE) that automatically selects the steplength and order for each solution step. The solution procedure for the sensitivity coefficients maintains accuracy and stability by using exactly the same steplengths and numerical approximations. The method computes sensitivity coefficients with respect to any combination of the initial values of the thermochemical variables and the three rate constant parameters for the chemical reactions. The method is illustrated by application to several simple problems and, where possible, comparisons are made with exact solutions and those obtained by other techniques.
Transient pulse response experiments are used to construct rate/concentration (RC) kinetic dependencies—RC petals—and provide a new method to distinguish the timing and interplay of adsorption, surface reaction, and product formation on complex (industrial) materials. A petal shape arises as the dynamic “reaction-diffusion” experiment forces the concentration and reaction rate to return to zero. In contrast to the typical steady-state “Langmuir-type” RC dependence, RC petals have two branches, which arise as a result of decoupled gas and surface concentrations in the non-steady-state regime. To demonstrate this approach, the characteristics of petal shapes using ammonia decomposition as a probe reaction are presented. Ammonia, hydrogen, and nitrogen transformation rates are compared on three simple materials: Fe, Co, and a bimetallic CoFe preparation when ammonia is pulsed at 550°C in a low-pressure diffusion reactor. All materials demonstrate a two-branch kinetic RC dependence for ammonia adsorption, and rate constants are quantified in the low-coverage regime. We found that H 2 and N 2 product formation was dependent on the concentration of surface intermediates for all materials with one exception: for Co, an additional fast hydrogen generation process was observed; the rate of which coincided with ammonia adsorption. Nitrogen generation was only significant for CoFe and Co and on the CoFe catalyst, a self-inhibition property was observed. A method for estimating the number of active sites based on the RC petals is presented and was applied to the Fe and CoFe samples. The surface coverage and rate of formation/conversion of surface intermediates are interpreted from the examination of shape characteristics of the RC petals for each material.
A gas-kinetic solver is developed for the ideal magnetohydrodynamics (MHD) equations. The new scheme is based on the direct splitting of the flux function of the MHD equations with the inclusion of "particle" collisions in the transport process. Consequently, the artificial dissipation in the new scheme is much reduced in comparison with the MHD Flux Vector Splitting Scheme. At the same time, the new scheme is compared with the well-developed Roe-type MHD solver. It is concluded that the kinetic MHD scheme is more robust and efficient than the Roe- type method, and the accuracy is competitive. In this paper the general principle of splitting the macroscopic flux function based on the gas-kinetic theory is presented. The flux construction strategy may shed some light on the possible modification of AUSM- and CUSP-type schemes for the compressible Euler equations, as well as to the development of new schemes for a non-strictly hyperbolic system.
The efficiency of several algorithms used for numerical integration of stiff ordinary differential equations was compared. The methods examined included two general purpose codes EPISODE and LSODE and three codes (CHEMEQ, CREK1D and GCKP84) developed specifically to integrate chemical kinetic rate equations. The codes were applied to two test problems drawn from combustion kinetics. The comparisons show that LSODE is the fastest code available for the integration of combustion kinetic rate equations. It is shown that an iterative solution of the algebraic energy conservation equation to compute the temperature can be more efficient then evaluating the temperature by integrating its time-derivative.
This work describes the use of a cell phone camera and the L*a*b method (color space specified by the International Commission on Illumination) to characterize the color change in different vapochromic systems. In this study we have developed a semi-automatic color change analysis software that digitally analyzes images (e.g., video frames) collected while a vapochromic material is absorbing vapor. The advantages of using this method, as compared to reflectance spectroscopy or transmission spectroscopy through a thin film, include low cost, convenience, portability, ease of sample preparation, the absence of need for specialized equipment, and the ease of simultaneously collecting data on different samples under identical conditions. In addition, this method arguably provides direct insight into what a human would observe when monitoring these color changes by eye. Limitations of the method also are discussed in this paper.
A comparison of the efficiency of several algorithms recently developed for the efficient numerical integration of stiff ordinary differential equations is presented. The methods examined include two general-purpose codes EPISODE and LSODE and three codes (CHEMEQ, CREK1D, and GCKP84) developed specifically to integrate chemical kinetic rate equations. The codes are applied to two test problems drawn from combustion kinetics. The comparisons show that LSODE is the fastest code currently available for the integration of combustion kinetic rate equations. An important finding is that an iterative solution of the algebraic energy conservation equation to compute the temperature can be more efficient than evaluating the temperature by integrating its time-derivative.
A novel numerical method is employed to compute the integral form of the axi-symmetric Trubnikov-Rosenbluth potentials. Two methods for quadrature in pitch-angle are described and their convergence properties are studied. Careful attention is given to quadrature over a singular Green's function. Here it is shown that an infinite series representation of the Green's function can be used more efficiently than its closed form involving complete elliptic integrals. Then a collocation method in speed, with its associated quadrature scheme, is laid out and its convergence properties are studied. Using the proposed scheme, accurate low-order moments of the field collision operator are obtained using relatively few velocity space degrees of freedom. The scheme is showcased by solving for the equilibrium, axi-symmetric bootstrap current in tokamaks. A C 0 Gauss-Lobatto-Legendre finite element pitch-angle basis with vertex nodes at the trapped/passing boundary is shown, in the context of the integral methods used, to be much more efficient than the more common Legendre polynomial expansion.
The mass dependence of the transverse flow for Z = 1–5 fragments from the collisions of 40 Ar + 27 Al, 40 Ar + 48 Ti, and 40 Ar + 58 Ni at 47 MeV/nucleon is investigated experimentally in this article. The transverse flow values are determined using the in-plane components of the fragment transverse momenta, where three conventional methods, i.e., the kinetic flow tensor method, the transverse momentum analysis method, and the azimuthal correlation method, are applied to reconstruct the reaction plane in an event-by-event basis. It is demonstrated from the comparison of the present experimental mass dependent flow measurements and the model simulations using an improved antisymmetrized molecular dynamics model that the experimentally observed abnormal α transverse flow enhancement is closely related to the reaction plane reconstruction procedure in the flow extraction. We further investigate the physical existence of the abnormal α flow behavior using a two-particle azimuthal correlation method, which allows us to provide the relative flow magnitude information with an identification of fragment charge number without the knowledge of the reaction plane differing from the three conventional methods. It is found that the relative flow magnitudes deduced from the two-particle azimuthal correlation functions with an identification of Z, with the correction for the recoil effect imposed by the momentum conservation, show a monotonically increasing trend as a function of fragment charge number, with no exception of the α flow enhancement. These results, in addition to those from the improved antisymmetrized molecular dynamics model simulations, definitely provide experimental evidences for the inexistence of the abnormal α flow behavior in the heavy-ion collisions at the present incident energy region in nature.
A useful method for kinetic study of the hydrolysis of polyester elastomer is established which uses the number-average molecular weight. The reasonableness of this method is confirmed and the effect of magnetic particles on hydrolysis is considered.
Physical mixtures of catalysts can sometimes show synergistic activity which exceeds the sum of the activities of the same amount of catalysts employed separately. However, identifying such synergistic mixtures is non-trivial. Not all mixtures display synergy, and the number of combinations that are possible (even for binary mixtures of equal portions) scales very rapidly, and thus they would be costly to screen experimentally. In this work, we show that it is possible to predict synergistic mixtures using combinatorial kinetic simulations based on experimental data collected on individual catalysts. The data was first collected for conditions-of-interest under low conversions (also called near differential conditions, such that each condition approximates a small volume in a non-differential conditions reactor) to build a library of kinetic models (one model for each catalyst). This data was then used for combinatorial kinetic simulations of non-differential reactor scale conversions, and successfully predicted the qualitative behavior of two synergistic physical mixtures. The capability was utilized in the present work in the context of converting CO, C 3 H 6 , and NO species in simulated car exhaust, but is a general approach. We provide equations to estimate the costs of the screening method as well as the combinatorial kinetics method and show that the costs of identifying synergistic physical mixtures become much lower with the combinatorial kinetics method even before ten catalysts of interest. The cost of using prediction by the combinatorial kinetics method continue to become lower (relative to screening) when more catalysts are added -- and even more so if mixtures beyond two components, or consideration of multiple catalyst ratios, is of interest.
The American Society for Testing and Materials has published a new standard test method for characterizing time and temperature-dependence of material outgassing kinetics and the deposition kinetics of outgassed species on surfaces at various temperatures. This new ASTM standard, E 1559(1), uses the quartz crystal microbalance (QCM) collection measurement approach. The test method was originally developed under a program sponsored by the United States Air Force Materials Laboratory (AFML) to create a standard test method for obtaining outgassing and deposition kinetics data for spacecraft materials. Standardization by ASTM recognizes that the method has applications beyond aerospace. In particular, the method will provide data of use to the electronics, semiconductor, and high vacuum industries. In ASTM E 1559 the material sample is held in vacuum in a temperature-controlled effusion cell, while its outgassing flux impinges on several QCM's which view the orifice of the effusion cell. Sample isothermal total mass loss (TML) is measured as a function of time from the mass collected on one of the QCM's which is cooled by liquid nitrogen, and the view factor from this QCM to the cell. The amount of outgassed volatile condensable material (VCM) on surfaces at higher temperatures is measured as a function of time during the isothermal outgassing test by controlling the temperatures of the remaining QCM's to selected values. The VCM on surfaces at temperatures in between those of the collector QCM's is determined at the end of the isothermal test by heating the QCM's at a controlled rate and measuring the mass loss from the end of the QCM's as a function of time and temperature. This reevaporation of the deposit collected on the QCM's is referred to as QCM thermogravimetric analysis. Isothermal outgassing and deposition rates can be determined by differentiating the isothermal TML and VCM data, respectively, while the evaporation rates of the species can be obtained as a function of temperature by differentiating the QCM thermogravimetric analysis data.
Different functions are used to account for turbulence strength in the atmospheric boundary layer for different stability regimes. These functions are one of the sources for differences among different atmospheric models' predictions and associated biases. Also, turbulence strength is underrepresented in some of the resistance formulations. Here, to address these issues with dry deposition, firstly we take advantage of three‐dimensional (3‐D) turbulence information in estimating resistances by proposing and validating a 3‐D turbulence velocity scale that is relevant for different stability regimes of boundary layer. Secondly, we hypothesize and validate that friction velocity measured by 3‐D sonic anemometer can be effectively replaced by the new turbulence velocity scale multiplied by the von Karman constant. Finally, we (a) present a set of resistance formulations for ozone (O 3 ) based on the 3‐D turbulence velocity scale; (b) intercompare estimations of such resistances with those obtained using existing formulations; and, (c) evaluate simulated O 3 fluxes using a single‐point dry deposition model against long‐term observations of O 3 fluxes at the Harvard Forest (MA) site. Results indicate that the new resistance formulations work very well in simulating surface latent heat and O 3 fluxes when compared to respective existing formulations and measurements at a decadal time scale. Findings from this research may help to improve the capability of dry deposition schemes for better estimation of dry deposition fluxes and create opportunities for the development of a community dry deposition model for use in regional/global air quality models.
Presentation to be presented during the Site-Directed Research and Development (SDRD) program FY 2021 review meeting (Webex), September 22–23, 2021.