Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Kinetic methods”

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 109 records · Page 6

CFD Modeling of Phase Change and Pressure Drop during Violent Sloshing of Cryogenic Fluid in a Small-Scale Tank

This paper presents the development of a two-phase CFD model that is used to study sloshing of a cryogenic fluid (LN2) undergoing phase change generated by lateral acceleration and its effects on the ensuing heat and mass transfer and pressure drop in a small-scale tank. In this model the interface is captured using the volume-of-fluid (VOF) method [1]. Kinetic-theory based Schrage relation [2] is used for calculating phase change mass transfer at the liquid-vapor interface. Computational results are compared to the data provided by a non-isothermal sloshing experiment with phase change conducted by the University of Tokyo and JAXA in 2018 using liquid Nitrogen in a transparent tank. The effect of different thicknesses of a liquid temperature stratification layer created at the interface on the pressure drop during sloshing was studied experimentally. CFD results for interface movement and tank pressure are presented and compared in this paper to experimental data for the cases with the smallest and largest liquid temperature stratification thickness.

phase change↗

Elementary solutions of coupled model equations in the kinetic theory of gases

The method of elementary solutions is employed to solve two coupled integrodifferential equations sufficient for determining temperature-density effects in a linearized BGK model in the kinetic theory of gases. Full-range completeness and orthogonality theorems are proved for the developed normal modes and the infinite-medium Green's function is constructed as an illustration of the full-range formalism. The appropriate homogeneous matrix Riemann problem is discussed, and half-range completeness and orthogonality theorems are proved for a certain subset of the normal modes. The required existence and uniqueness theorems relevant to the H matrix, basic to the half-range analysis, are proved, and an accurate and efficient computational method is discussed. The half-space temperature-slip problem is solved analytically, and a highly accurate value of the temperature-slip coefficient is reported.

Kriese, J. T.↗

A pressure correction method for the calculation of compressible chemical reacting flows

A recently developed noniterative method for the solution of the transient fluid flow equations at all speed is extended to handle chemical reacting flows. The species conservation equations are loosely coupled into the predictor/multicorrector sequence of the solution procedure. A split-operator method separates the chemical kinetics terms from the fluid-dynamical terms, as well as an implicit differencing method enhance the numerical stability. The method was applied for turbulent diffusion flame calculations and for the analyses of high pressure, axisymmetric turbulent hypersonic nozzle flows. The diffusion flame results were compared with a similar pressure method for fast chemistry integration scheme without operator-splitting. Simulations of the nozzle flow indicated that the nonideal intermolecular effects must be included in the analysis and design of high pressure hypersonic nozzle.

Chen, Z. J.↗

Development of an Innovative Algorithm for Aerodynamics-Structure Interaction Using Lattice Boltzmann Method

The lattice Boltzmann equation (LBE) is a kinetic formulation which offers an alternative computational method capable of solving fluid dynamics for various systems. Major advantages of the method are owing to the fact that the solution for the particle distribution functions is explicit, easy to implement, and the algorithm is natural to parallelize. In this final report, we summarize the works accomplished in the past three years. Since most works have been published, the technical details can be found in the literature. Brief summary will be provided in this report. In this project, a second-order accurate treatment of boundary condition in the LBE method is developed for a curved boundary and tested successfully in various 2-D and 3-D configurations. To evaluate the aerodynamic force on a body in the context of LBE method, several force evaluation schemes have been investigated. A simple momentum exchange method is shown to give reliable and accurate values for the force on a body in both 2-D and 3-D cases. Various 3-D LBE models have been assessed in terms of efficiency, accuracy, and robustness. In general, accurate 3-D results can be obtained using LBE methods. The 3-D 19-bit model is found to be the best one among the 15-bit, 19-bit, and 27-bit LBE models. To achieve desired grid resolution and to accommodate the far field boundary conditions in aerodynamics computations, a multi-block LBE method is developed by dividing the flow field into various blocks each having constant lattice spacing. Substantial contribution to the LBE method is also made through the development of a new, generalized lattice Boltzmann equation constructed in the moment space in order to improve the computational stability, detailed theoretical analysis on the stability, dispersion, and dissipation characteristics of the LBE method, and computational studies of high Reynolds number flows with singular gradients. Finally, a finite difference-based lattice Boltzmann method is developed for inviscid compressible flows.

Mei, Ren-Wei↗

Hydrogen chemistry - Perspective on experiment and theory

A review is presented of the advantages and limitations of various experimental methods for the investigation of the kinetics of hydrogen chemistry, including classic thermal and photochemical methods and the crossed molecular beam method. Special attention is given to the flash photolysis-resonance fluorescence apparatus developed by Braun et al, in which repetitive vacuum UV flashes result in the photolytic generation of the desired species, and to the discharge-flow technique. The use of various theoretical methods for the selection or elimination of kinetic data is considered in a brief discussion of the rate theory of two-body encounters and recombination-dissociation processes in neutral reactions. Recent kinetic studies of a series of OH reactions and of a major loss process for odd H in the stratosphere are summarized.

Kaufman, F.↗

Kinetic approach to the evaporation and condensation problem

In the paper, the Boltzmann equation governing the evaporation and condensation phenomena is solved by the Monte Carlo method. Based on the kinetic theory of gas the role of the non-equilibrium Knudsen layer and the growth of the hydrodynamic region outside the layer as time proceeds are simulated. Results show two possible types of transient developments in the vapor phase. The effects of the molecular absorption coefficient of the phase surface are examined. Except in the case of very strong evaporation the kinematic effects of binary collisions among vapor molecules on the mass flux rate are not serious. The limiting case of the quasi-steady evaporation and the maximal value of the evaporation rate are obtained.

Murakami, M.↗

Kinetic tetrazolium microtiter assay

A method for conducting an in vitro cell assay using a tetrazolium indicator is disclosed. The indicator includes a nonionic detergent which solubilizes a tetrazolium reduction product in vitro and has low toxicity for the cells. The incubation of test cells in the presence of zolium bromide and octoxynol (TRITON X-100) permits kinetics of the cell metabolism to be determined.

Pierson, Duane L.↗

Design, Construction and Testing of Annular Diffusers for High Speed Civil Transportation Combustor Applications

A theoretical and an experimental design study of subsonic flow through curved-wall annular diffusers has been initiated under this award in order to establish the most pertinent design parameters and hence performance characteristics for such devices, an the implications of their application in the design of engine components in the aerospace industries. The diffusers under this study are expected to contain appreciable regions of stall and the effects of swirl on their performance are being studied. The experimental work involves the application of Computer Aided Design software tool to the development of a suitable annular diffuse geometry and the subsequent downloading of such data to a CNC machine at Central State University (CSU). Two experimental run segments have been completed so far during FY-95 involving flow visualization and diffuser performance evaluation based on Kinetic Energy Dissipation. The method of calculation of the performance of diffusers based on pressure recovery coefficient has been shown to have some shortcomings and so the kinetic energy dissipation approach has been introduced in the run segment two with some success. The application of the discretized, full Navier Stokes and Continuity equations to the numerical study of the problem described above for the time-mean flow is expected to follow. Various models of turbulence are being evaluated for adoption throughout the study and comparisons would be made with experimental data where they exist. Assessment of diffuser performance based on the dissipated mechanical energy would also be made. The result of the investigations are expected to indicate that more cost effective component design of such devices as diffusers which normally contain complex flows can still be achieved.

Okhio, Cyril B.↗

Kinetic Monte Carlo Simulation of Oxygen and Cation Diffusion in Yttria-Stabilized Zirconia

Yttria-stabilized zirconia (YSZ) is of interest to the aerospace community, notably for its application as a thermal barrier coating for turbine engine components. In such an application, diffusion of both oxygen ions and cations is of concern. Oxygen diffusion can lead to deterioration of a coated part, and often necessitates an environmental barrier coating. Cation diffusion in YSZ is much slower than oxygen diffusion. However, such diffusion is a mechanism by which creep takes place, potentially affecting the mechanical integrity and phase stability of the coating. In other applications, the high oxygen diffusivity of YSZ is useful, and makes the material of interest for use as a solid-state electrolyte in fuel cells. The kinetic Monte Carlo (kMC) method offers a number of advantages compared with the more widely known molecular dynamics simulation method. In particular, kMC is much more efficient for the study of processes, such as diffusion, that involve infrequent events. We describe the results of kinetic Monte Carlo computer simulations of oxygen and cation diffusion in YSZ. Using diffusive energy barriers from ab initio calculations and from the literature, we present results on the temperature dependence of oxygen and cation diffusivity, and on the dependence of the diffusivities on yttria concentration and oxygen sublattice vacancy concentration. We also present results of the effect on diffusivity of oxygen vacancies in the vicinity of the barrier cations that determine the oxygen diffusion energy barriers.

Good, Brian↗

The PDF method for turbulent combustion

Probability Density Function (PDF) methods provide a means of calculating the properties of turbulent reacting flows. They have been successfully applied to many turbulent flames, including some with finite rate kinetic effects. Here the methods are reviewed with an emphasis on computational issues and their application to turbulent combustion.

Pope, S. B.↗

Kinetic Monte Carlo Simulation of Cation Diffusion in Low-K Ceramics

Low thermal conductivity (low-K) ceramic materials are of interest to the aerospace community for use as the thermal barrier component of coating systems for turbine engine components. In particular, zirconia-based materials exhibit both low thermal conductivity and structural stability at high temperature, making them suitable for such applications. Because creep is one of the potential failure modes, and because diffusion is a mechanism by which creep takes place, we have performed computer simulations of cation diffusion in a variety of zirconia-based low-K materials. The kinetic Monte Carlo simulation method is an alternative to the more widely known molecular dynamics (MD) method. It is designed to study "infrequent-event" processes, such as diffusion, for which MD simulation can be highly inefficient. We describe the results of kinetic Monte Carlo computer simulations of cation diffusion in several zirconia-based materials, specifically, zirconia doped with Y, Gd, Nb and Yb. Diffusion paths are identified, and migration energy barriers are obtained from density functional calculations and from the literature. We present results on the temperature dependence of the diffusivity, and on the effects of the presence of oxygen vacancies in cation diffusion barrier complexes as well.

Good, Brian↗

Growth kinetics of tetragonal lysozyme crystals

A method for immobilizing protein crystals in small volumes to determine growth rates on various faces is applied to study the growth kinetics of the (100) face of tetragonal hen-egg white lysozyme crystals at different degrees of bulk saturation. In normal gravity, transport is found to be dominated by convection for crystal sizes larger than a few microns, while in a microgravity environment, transport is diffusion-limited for sizes up to a few mm. It is found that convection can be significant even in microgravity for crystals approaching cm sizes, and that lysozyme growth is limited by surface kinetics in normal gravity.

Pusey, M.↗

Kinetics of Fast Atoms in the Terrestrial Atmosphere

This report summarizes our investigations performed under NASA Grant NAG5-8058. The three-year research supported by the Geospace Sciences SR&T program (Ionospheric, Thermospheric, and Mesospheric Physics) has been designed to investigate fluxes of energetic oxygen and nitrogen atoms in the terrestrial thermosphere. Fast atoms are produced due to absorption of the solar radiation and due to coupling between the ionosphere and the neutral thermospheric gas. We have investigated the impact of hot oxygen and nitrogen atoms on the thermal balance, chemistry and radiation properties of the terrestrial thermosphere. Our calculations have been focused on the accurate quantitative description of the thermalization of O and N energetic atoms in collisions with atom and molecules of the ambient neutral gas. Upward fluxes of oxygen and nitrogen atoms, the rate of atmospheric heating by hot oxygen atoms, and the energy input into translational and rotational-vibrational degrees of atmospheric molecules have been evaluated. Altitude profiles of hot oxygen and nitrogen atoms have been analyzed and compared with available observational data. Energetic oxygen atoms in the terrestrial atmosphere have been investigated for decades, but insufficient information on the kinetics of fast atmospheric atoms has been a main obstacle for the interpretation of observational data and modeling of the hot geocorona. The recent development of accurate computational methods of the collisional kinetics is seen as an important step in the quantitative description of hot atoms in the thermosphere. Modeling of relaxation processes in the terrestrial atmosphere has incorporated data of recent observations, and theoretical predictions have been tested by new laboratory measurements.

Kharchenko, Vasili A.↗

Colloidal Oxide Properties and Kinetic Studies

Silica based oxide particles prepared by various methods are compared in terms of their gellation kinetics, densification kinetics, and microable structural evolutions begin at very low temperature. Additional work reported in literature further suggested that the colloidal behaviors and the dried gel structures of the same oxide from different preparation techniques differ fundamentally. These differences significantly affect the densification kinetics and the product microstructure. Gravitational segregation of particles of different sizes in a gravitational field is modelled as a diffusion/segregation problem. A structural free volume model was used to describe the activation process of jump statistics. The sizes of the particles, the interstices, and their packing topology are accounted for in the estimation of the free volume. An attempt frequency is ascribed to particles to account for the vibrational rearrangement during powder packing.

Chen, I. W.↗

Innovative methods for the measurement of I* quantum yields and kinetics by diode laser gain-versus-absorption

The quantum yields of a variety of candidate molecules for solar lasant materials to produce I* were tested. The absorption spectrum was measured for each compound and the I* yield determined by the diode laser or by infrared emission, using C3F7I as a standard. The results of these measurements are summarized. A GaAsInP diode laser system was developed to probe I and I* atoms to obtain yields and kinetics. A technique of gain-versus-absorption spectroscopy was investigated to measure quantum yields with high accuracy. The errors in the yield data were reduced to +/- 2% or less. In addition, experiments were set up to measure the rates of F-sublevel changing collisions in both the I ground state and the I* excited state. Finally, experiments and modelling were carried out to explore the possibility of measuring the recombination rates of I* with C3F7 radicals.

Leone, Stephen R.↗

Approximate method for predicting the permanent set in a beam in vacuo and in water subject to a shock wave

An approximate method to compute the maximum deformation and permanent set of a beam subjected to shock wave laoding in vacuo and in water was investigated. The method equates the maximum kinetic energy of the beam (and water) to the elastic plastic work done by a static uniform load applied to a beam. Results for the water case indicate that the plastic deformation is controlled by the kinetic energy of the water. The simplified approach can result in significant savings in computer time or it can expediently be used as a check of results from a more rigorous approach. The accuracy of the method is demonstrated by various examples of beams with simple support and clamped support boundary conditions.

Stiehl, A. L.↗

Supporting technology for the development of Controlled Ecological Life Support Systems (CELSS)

To support the development of Controlled Ecological Life Support Systems (CELSS) in the space program, a metabolic simulator has been selected for use in a closed chamber to test functions of the CELSS. This metabolic simulator is a catalytic reactor which oxidizes the methyl acetate to produce carbon dioxide and water vapor. In this project, kinetic studies of catalytic oxidation of methyl acetate were conducted using monolithic and pellet catalysts with 0.5% (by weight) platinum (Pt) on aluminum oxide (Al2O3). The reaction was studied at a pressure of one atmosphere and at temperatures varying from 160 C to 420 C. By-products were identified at the exit of the preheater and reactor. For the kinetic study with the monolithic catalyst, a linear regression method was used to correlate the kinetic data with zero-order, first-order and Langmuir-Hinshelwood models. Results indicate that the first-order model represents the data adequately at low concentrations of methyl acetate. For higher concentrations of methyl acetate, the Langmuir-Hinshelwood model best represents the kinetic data. Both rate constant and adsorption equilibrium constants were estimated from the regression. A Taguchi orthogonal array (L(sub 9)) was used to investigate the effects of temperature, flow rate, and concentration on the catalytic oxidation of methyl acetate. For the monolithic catalyst, temperature exerts the most significant effect, followed by concentration of methyl acetate. For the pellet catalyst, reaction temperature is the most significant factor, followed by gas flow rate and methyl acetate concentration. Concentrations of either carbon dioxide or oxygen were seen to have insignificant effect on the methyl acetate conversion process. Experimental results indicate that the preheater with glass beads can accomplish thermal cracking and catalytic reaction of methyl acetate to produce acetic acid, methanol, methyl formate, and 1-propanol. The concentration of all by-products was measured in ppmv (parts per million by volume). At higher temperatures, greater amounts of these products are produced, as expected. In all cases, methanol was the predominant concentration detected, followed by methyl formate. At temperatures lower than 320 C for the P-type monolithic catalyst, methanol, acetic acid, and acetone were detected, whereas, for the E-type monolithic catalyst, only methanol was detected at 160 C. Both P and E types of the monolithic catalyst were specified with the same substrates (ceramic), washcoat (Al2O3), and promoter (Pt). However, the manufacturing and treatment procedures were quite different. It was therefore concluded that the performance of the E-type monolithic catalyst is superior to that of the P-type for oxidation of methyl acetate. At higher reaction temperatures, e.g., above 420 C, all reactants and byproducts were completely oxidized using these two types of monolithic catalyst to produce carbon dioxide and water vapor. A complex heterogenous catalytic reaction mechanism was proposed to explain the formation of the byproducts (methanol, acetic acid, and methyl formate) as the methyl acetate traveled through the preheater packed with glass beads. The by-product, 1-propanol, may be formed only through a homogeneous reaction, since it is difficult to develop a reasonable sequence of heterogeneous reaction steps to explain its formation. The homogeneous thermal decomposition of methyl acetate to form free radicals was proposed to explain the formation of 1-propanol, and also methanol, in the preheater. A dual-site catalytic reaction mechanism was proposed for the oxidation of methyl acetate over Pt/Al2O3 monolithic catalyst. The dual-site mechanism describes the chemisorption of oxygen molecules as well as a physical adsorption of methyl acetate on the active sites. On the active sites, methyl acetate is oxidized rapidly to form carbon dioxide and water vapor. A rate equation derived from this mechanism gives the Langmuir-Hinshelwood rate formula which has been observed from the experimental data obtained in this project for high methyl acetate concentration (greater than 1000 ppmv) over a monolithic catalyst. If the oxygen concentration is very high and methyl acetate concentration is very low, the reaction rate equation is then reduced to a first-order with respect to methyl acetate concentration. The first-order model has also been observed from the experimental data obtained in this project for low methyl acetate concentration (less than 1000 ppmv).

Li, Ku-Yen↗

Assessment of three numerical methods for the computation of a low-density plume flow

Results from three numerical methods including one based on the Navier-Stokes equations, one based on kinetic theory using the DSMC method, and one based on the Boltzmann equation with a Krook-type collision term are compared to each other and to experimental data for a model problem of heated nitrogen flow in a conical nozzle expanding into a vacuum. The problem simulates flow in a resistojet, a low-thrust, electrothermal rocket. The continuum method is applied to both the internal flow and near-field plume. The DSMC and Boltzmann methods are applied primarily to the plume. Experimental measurements of Pitot pressure and flow angle, taken with an apparatus that duplicates the model nozzle flow, are used in the comparisons.

Penko, Paul F.↗