Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “kinetic”

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 379 records · Page 21

Kinetic aspects of magnetotail dynamics - Observations

The observed kinetic (as opposed to fluid) aspects of particle behavior in the geomagnetic tail are discussed. Such behavior has been observed from thermal plasma energies to energies a few MeV or more. The importance of recognizing and analyzing the kinetic behavior of tail particles is that it provides a powerful diagnostic tool for inferring the state of the local and remote environment and the mechanisms responsible for observed tail dynamics, and a way of determining when fluid analyses (relying upon moments such as bulk velocity and temperature) may be misleading. Thus kinetic analyses provide a valuable complement to fluid-dynamics-based studies of the collective behavior of geomagnetic tail particles. Observed kinetic effects can be analyzed to examine spatial gradients; convection electric fields, distances to source regions, source-region characteristics, magnetic-field topology, heating/acceleration mechanism, boundary characteristics, and remote mechanisms affecting local particle distributions.

Mitchell, D. G.↗

Turbulent hydrocarbon combustions kinetics - Stochastic modeling and verification

Idealized reactors, that are designed to ensure perfect mixing and are used to generate the combustion kinetics for complex hydrocarbon fuels, may depart from the ideal and influence the kinetics model performance. A complex hydrocarbon kinetics model that was established by modeling a jet-stirred combustor (JSC) as a perfectly stirred reactor (PSR), is reevaluated with a simple stochastic process in order to introduce the unmixedness effect quantitatively into the reactor system. It is shown that the comparisons of the predictions and experimental data have improved dramatically with the inclusion of the unmixedness effect in the rich combustion region. The complex hydrocarbon kinetics is therefore verified to be mixing effect free and be applicable to general reacting flow calculations.

Wang, T. S.↗

Excitation of kinetic Alfven waves by resonant mode conversion and longitudinal heating of magnetized plasmas

The excitation of the kinetic Alfven wave by resonant mode conversion and longitudinal heating of the plasma by the kinetic Alfven wave were demonstrated on the basis of a macroscale particle simulation. The longitudinal electron current was shown to be cancelled by the ions. The kinetic Alfven wave produced an ordered motion of the plasma particles in the wave propagation direction. The electrons were pushed forward along the ambient magnetic field by absorbing the kinetic Alfven wave through the Landau resonance.

Tanaka, Motohiko↗

Ion beam generation at the plasma sheet boundary layer by kinetic Alfven waves

A two-dimensional quasi-linear numerical code was developed for studying ion beam generation at the plasma sheet boundary layer by kinetic Alfven waves. The model assumes that the central plasma sheet is the particle source, and that the last magnetic field lines on which kinetic Alfven waves exist and diffusion occurs can be either open or closed. As the possible source for the excitement of the kinetic Alfven waves responsible for ion diffusion, the resonant mode conversion of the surface waves to kinetic Alfven waves is considered. It is shown that, depending on the topology of the magnetic field at the lobe side of the simulation system, i.e., on whether field lines are open or closed, the ion distribution function may or may not reach a steady state.

Moghaddam-Taaheri, E.↗

Kinetic effects of Alfven wave nonlinearity. I - Ponderomotive density fluctuations

The Vlasov theory is used to study kinetic corrections to fluid descriptions of Alfven wave nonlinearity. The method is to obtain an expression for the second-order perturbed distribution function produced by a nonlinear Alfven wave. From this distribution function a kinetically correct expression is obtained for the plasma density perturbation associated with an envelope-modulated Alfven wave. This kinetic theory result differs substantially from the fluid expression when the plasma beta is greater than about 1, and the electron and ion temperatures are approximately equal. This result is of interest because density fluctuations are an observationally accessible indicator of wave nonlinearity in solar system Alfven waves. It also will assist in the determination of properties of Alfven waves in the interstellar medium. Finally, this analysis also yields a kinetically correct expression for u, the magnetic field-aligned component of the plasma fluid velocity.

Spangler, Steven R.↗

Critical evaluation of Jet-A spray combustion using propane chemical kinetics in gas turbine combustion simulated by KIVA-2

Jet-A spray combustion has been evaluated in gas turbine combustion with the use of propane chemical kinetics as the first approximation for the chemical reactions. Here, the numerical solutions are obtained by using the KIVA-2 computer code. The KIVA-2 code is the most developed of the available multidimensional combustion computer programs for application of the in-cylinder combustion dynamics of internal combustion engines. The released version of KIVA-2 assumes that 12 chemical species are present; the code uses an Arrhenius kinetic-controlled combustion model governed by a four-step global chemical reaction and six equilibrium reactions. Researchers efforts involve the addition of Jet-A thermophysical properties and the implementation of detailed reaction mechanisms for propane oxidation. Three different detailed reaction mechanism models are considered. The first model consists of 131 reactions and 45 species. This is considered as the full mechanism which is developed through the study of chemical kinetics of propane combustion in an enclosed chamber. The full mechanism is evaluated by comparing calculated ignition delay times with available shock tube data. However, these detailed reactions occupy too much computer memory and CPU time for the computation. Therefore, it only serves as a benchmark case by which to evaluate other simplified models. Two possible simplified models were tested in the existing computer code KIVA-2 for the same conditions as used with the full mechanism. One model is obtained through a sensitivity analysis using LSENS, the general kinetics and sensitivity analysis program code of D. A. Bittker and K. Radhakrishnan. This model consists of 45 chemical reactions and 27 species. The other model is based on the work published by C. K. Westbrook and F. L. Dryer.

Nguyen, H. L.↗

Interfacial reaction kinetics of coated SiC fibers with various titanium alloys

The kinetics of the reaction between the silicon carbide fibers and the titanium-based alloy matrix was investigated at temperatures from 800 to 1000 C for several titanium-based alloys (including Ti-1100 alloy and BETA 21S) and unalloyed Ti, reinforced with coated silicon carbide fiber SCS-6. The reaction zone growth kinetics was studied by exposing vacuum encapsulated samples to temperatures from 700 to 1000 C for times up to 150 hrs, followed by SAM observations of samples which were polished perpendicular to the fiber axis and etched. It was found that the reaction zone growth kinetics of the alpha (hcp) and beta (bcc) phases of unalloyed titanium reacting with SCS-6 fibers exhibited different values of the apparent activation energy and of the preexponential factor. Additions of other metals to Ti was found to slow down the reaction kinetics. Among the alloys studied, the Ti-1100 was the slowest reacting conventional alloy and the Ti-14Al-21Nb (in wt pct) was the slowest overall.

Gundel, D. B.↗

An upwind, kinetic flux-vector splitting method for flows in chemical and thermal non-equilibrium

We have developed new upwind kinetic difference schemes for flows with non-equilibrium thermodynamics and chemistry. These schemes are derived from the Boltzmann equation with the resulting Euler schemes developed as moments of the discretized Boltzmann scheme with a locally Maxwellian velocity distribution. Splitting the velocity distribution at the Boltzmann level is seen to result in a flux-split Euler scheme and is called Kinetic Flux Vector Splitting (KFVS). Extensions to flows with finite-rate chemistry and vibrational relaxation is accomplished utilizing nonequilibrium kinetic theory. Computational examples are presented comparing KFVS with the schemes of Van Leer and Roe for a quasi-one-dimensional flow through a supersonic diffuser, inviscid flow through two-dimensional inlet, and viscous flow over a cone at zero angle-of-attack. Calculations are also shown for the transonic flow over a bump in a channel and the transonic flow over an NACA 0012 airfoil. The results show that even though the KFVS scheme is a Riemann solver at the kinetic level, its behavior at the Euler level is more similar to the existing flux-vector splitting algorithms than to the flux-difference splitting scheme of Roe.

Eppard, W. M.↗

LSENS: A General Chemical Kinetics and Sensitivity Analysis Code for homogeneous gas-phase reactions. Part 3: Illustrative test problems

LSENS, the Lewis General Chemical Kinetics and Sensitivity Analysis Code, has been developed for solving complex, homogeneous, gas-phase chemical kinetics problems and contains sensitivity analysis for a variety of problems, including nonisothermal situations. This report is part 3 of a series of three reference publications that describe LSENS, provide a detailed guide to its usage, and present many example problems. Part 3 explains the kinetics and kinetics-plus-sensitivity analysis problems supplied with LSENS and presents sample results. These problems illustrate the various capabilities of, and reaction models that can be solved by, the code and may provide a convenient starting point for the user to construct the problem data file required to execute LSENS. LSENS is a flexible, convenient, accurate, and efficient solver for chemical reaction problems such as static system; steady, one-dimensional, inviscid flow; reaction behind incident shock wave, including boundary layer correction; and perfectly stirred (highly backmixed) reactor. In addition, the chemical equilibrium state can be computed for the following assigned states: temperature and pressure, enthalpy and pressure, temperature and volume, and internal energy and volume. For static problems the code computes the sensitivity coefficients of the dependent variables and their temporal derivatives with respect to the initial values of the dependent variables and/or the three rate coefficient parameters of the chemical reactions.

Bittker, David A.↗

LSENS, a general chemical kinetics and sensitivity analysis code for homogeneous gas-phase reactions. 2: Code description and usage

LSENS, the Lewis General Chemical Kinetics Analysis Code, has been developed for solving complex, homogeneous, gas-phase chemical kinetics problems and contains sensitivity analysis for a variety of problems, including nonisothermal situations. This report is part 2 of a series of three reference publications that describe LSENS, provide a detailed guide to its usage, and present many example problems. Part 2 describes the code, how to modify it, and its usage, including preparation of the problem data file required to execute LSENS. Code usage is illustrated by several example problems, which further explain preparation of the problem data file and show how to obtain desired accuracy in the computed results. LSENS is a flexible, convenient, accurate, and efficient solver for chemical reaction problems such as static system; steady, one-dimensional, inviscid flow; reaction behind incident shock wave, including boundary layer correction; and perfectly stirred (highly backmixed) reactor. In addition, the chemical equilibrium state can be computed for the following assigned states: temperature and pressure, enthalpy and pressure, temperature and volume, and internal energy and volume. For static problems the code computes the sensitivity coefficients of the dependent variables and their temporal derivatives with respect to the initial values of the dependent variables and/or the three rate coefficient parameters of the chemical reactions. Part 1 (NASA RP-1328) derives the governing equations describes the numerical solution procedures for the types of problems that can be solved by lSENS. Part 3 (NASA RP-1330) explains the kinetics and kinetics-plus-sensitivity-analysis problems supplied with LSENS and presents sample results.

Radhakrishnan, Krishnan↗

Direct numerical simulations of turbulent non-premixed methane-air flames modeled with reduced kinetics

Turbulent non-premixed stoichiometric methane-air flames modeled with reduced kinetics have been studied using the direct numerical simulation approach. The simulations include realistic chemical kinetics, and the molecular transport is modeled with constant Lewis numbers for individual species. The effect of turbulence on the internal flame structure and extinction characteristics of methane-air flames is evaluated. Consistent with earlier DNS with simple one-step chemistry, the flame is wrinkled and in some regions extinguished by the turbulence, while the turbulence is weakened in the vicinity of the flame due to a combination of dilatation and an increase in kinematic viscosity. Unlike previous results, reignition is observed in the present simulations. Lewis number effects are important in determining the local stoichiometry of the flame. The results presented in this work are preliminary but demonstrate the feasibility of incorporating reduced kinetics for the oxidation of methane with direct numerical simulations of homogeneous turbulence to evaluate the limitations of various levels of reduction in the kinetics and to address the formation of thermal and prompt NO(x).

Card, J. M.↗

Kinetic structure of rotational discontinuities: Implications for the magnetopause

Magnetic field rotations in the high ion beta magnetosheath that are part of the magnetopause structure are expected to have only a small normal component. We have studied the properties of rotational discontinuities (RDs) under these conditions, viewed as the limit of weak intermediate shocks (ISs), by performing hybrid simulations with a reflecting wall boundary condition (piston method). With this dynamic formation, the sense and size of rotation are not arbitrarily predetermined, but rather evolve from the given upstream (magnetosheath) and downstream (magnetospheric) boundary conditions, similar to what takes place at the magnetopause. This work focuses on several aspects: the observed minimum shear of RDs, their width, their internal signature, and their relation to ISs in isotropic plasmas. Our simulation results are in agreement with the minimum shear observations, that is, the RDs choose the sense of rotation that corresponds to the minimum angle between the upstream and downstream field vector. The RDs are stable, with a unique scale size. Typical gradient scale half widths are one to four ion inertial lengths with a total width up to ten times of that, in agreement with magnetopause observations. We develop a generalized fluid theory of RDs and discuss the characteristic internal signatures of the rotational layer, comparing the kinetic simulation results to predictions from the generalized fluid theory. The results show that ion inertia, anisotropic pressure, finite Larmor radius effects, nonzero ion heat flux, and reflected ions all contribute to the signatures of RDs on kinetic scales. The RDs may have upstream or downstream wave trains, which become weak for high ion beta and small normal components of the magnetic field. We explain the presence and direction of wave trains in terms of the kinetic properties of the Alfven/ion-cyclotron mode. Away from the RD limit there is a smooth transition to weak intermediate shocks, which have small jumps close to expected Rankine-Hugoniot values. Apart from that, there are few kinetic plasma signatures that distinguish RDs from their neighboring ISs. However, noncoplanar ISs evolve in time into thin RDs. Using the properties of RDs and ISs, we make specific suggestions how these discontinuities can be distinguished observationally in the case of an isotropic plasma.

Krauss-Varban, D.↗

Kinetic Study of Methyl Acetate Oxidation in a Pt/Al2O3 Fixed-Bed Reactor

To support technology development for future long-term missions, a metabolic simulator will be used in a closed chamber to test the functions of a Controlled Ecological Life Support System (CELSS). Methyl acetate (MA) was selected as the fuel because its metabolic respiratory quotient is near that of humans. A kinetic study of the catalytic oxidation of MA over Pt/Al203 was then conducted to support the design and operation of the simulator. Kinetic data were obtained as a conversion percentage of MA versus retention time. The reaction was studied at one atmosphere and temperatures from 220 to 340 deg. C. The inlet MA concentration was varied from 100 to 2000 ppm with retention times from 0.01 to 10 sec. A first-order rate law and a Langmuir-Hinshelwood rate equation were tested by nonlinear regression of the kinetic data to estimate rate constants in the rate law. Regression results of the L-H equation explain the kinetic data better than the results of the first-order rate law. A Taguchi experimental design was used to study the effects of temperature, retention time, and concentrations of MA, CO2, and O2 on the conversion of MA. Results indicate that temperature has greatest effect, followed by retention time, and finally MA concentration. It was further determined that the effects of CO2 and O2 concentrations, and the cross effects, are negligible.

Hoy, Michael↗

Metastable Solution Thermodynamic Properties and Crystal Growth Kinetics

The crystal growth rates of NH4H2PO4, KH2PO4, (NH4)2SO4, KAl(SO4)2 central dot 12H2O, NaCl, and glycine and the nucleation rates of KBr, KCl, NaBr central dot 2H2O, (NH4)2Cl, and (NH4)2SO4 were expressed in terms of the fundamental driving force of crystallization calculated from the activity of supersaturated solutions. The kinetic parameters were compared with those from the commonly used kinetic expression based on the concentration difference. From the viewpoint of thermodynamics, rate expressions based on the chemical potential difference provide accurate kinetic representation over a broad range of supersaturation. The rates estimated using the expression based on the concentration difference coincide with the true rates of crystallization only in the concentration range of low supersaturation and deviate from the true kinetics as the supersaturation increases.

Kim, Soojin↗

Kinetics and Mechanism of the Reaction of Hydoxyl Radicals with Acetonitrile under Atmospheric Conditions

The pulsed laser photolysis-pulsed laser induced fluorescence technique has been employed to determine absolute rate coefficients for the reaction OH + CH3CN (1) and its isotopic variants, OH + CD3CN (2), OD + CH3CN (3), and OD + CD3CN (4). Reactions 1 and 2 were studied as a function of pressure and temperature in N2, N2/O2, and He buffer gases. In the absence of O2 all four reactions displayed well-behaved kinetics with exponential OH decays and pseudo-first rate constants which were proportional to substrate concentration. Data obtained in N2 over the range 50-700 Torr at 298 K are consistent with k(sub 1), showing a small pressure dependence. The Arrhenius expression obtained by averaging data at all pressures in k(sub 1)(T) = (1.1(sup +0.5)/(sub -0.3)) x 10(exp -12) exp[(-1130 +/- 90)/T] cu cm /(molecule s). The kinetics of reaction 2 are found to be pressure dependent with k(sub 2) (298 K) increasing from (1.21 +/- 0.12) x 10(exp -14) to (2.16 +/- 0.11) x 10(exp -14) cm(exp 3)/ (molecule s) over the pressure range 50-700 Torr of N2 at 298 K. Data at pressures greater than 600 Torr give k(sub 2)(T) = (9.4((sup +13.4)(sub -5.0))) x 10(exp -13) exp[(-1180 +/- 250)/T] cu cm/(molecule s). The rates of reactions 3 and 4 are found to be independent of pressure over the range 50-700 Torr of N2 with 298 K rate coefficient given by k(sub 3) =(3.18 +/- 0.40) x 10(exp -14) cu cm/(molecule s) and k(sub 4) = (2.25 +/-0.28) x 10(exp -14) cu cm/(molecule s). In the presence of O2 each reaction shows complex (non-pseudo-first-order) kinetic behavior and/or an apparent decrease in the observed rate constant with increasing [O2], indicating the presence of significant OH or OD regeneration. Observation of regeneration of OH in (2) and OD in (3) is indicative of a reaction channel which proceeds via addition followed by reaction of the adduct, or one of its decomposition products, with O2. The observed OH and OD decay profiles have been modeled by using a simple mechanistic scheme to extract kinetic information about the adduct reations with O2 and branching ratios for OH regeneration. A plausible mechanism for OH regeneration in (2) involves OH addition to the nitrogen atom followed by O2 addition to the cyano carbon atom, isomeriazation and decomposition to D2CO + DOCN + OH. Our results suggest that the OH + CH3CN reaction occurs via a complex mechanism involving both bimolecular and termolecular pathways, analogous to the mechanisms for the the important atmospheric reactions of OH with CO and HNO3.

Hynes, A. J.↗

A Gas-Kinetic Scheme for Reactive Flows

In this paper, the gas-kinetic BGK scheme for the compressible flow equations is extended to chemical reactive flow. The mass fraction of the unburnt gas is implemented into the gas kinetic equation by assigning a new internal degree of freedom to the particle distribution function. The new variable can be also used to describe fluid trajectory for the nonreactive flows. Due to the gas-kinetic BGK model, the current scheme basically solves the Navier-Stokes chemical reactive flow equations. Numerical tests validate the accuracy and robustness of the current kinetic method.

Lian,Youg-Sheng↗

Time-Resolved O3 Chemical Chain Reaction Kinetics Via High-Resolution IR Laser Absorption Methods

Excimer laser photolysis in combination with time-resolved IR laser absorption detection of OH radicals has been used to study O3/OH(v = 0)/HO2 chain reaction kinetics at 298 K, (i.e.,(k(sub 1) is OH + 03 yields H02 + 02 and (k(sub 2) is H02 + 03 yields OH + 202). From time-resolved detection of OH radicals with high-resolution near IR laser absorption methods, the chain induction kinetics have been measured at up to an order of magnitude higher ozone concentrations ([03] less than or equal to 10(exp 17) molecules/cu cm) than accessible in previous studies. This greater dynamic range permits the full evolution of the chain induction, propagation, and termination process to be temporally isolated and measured in real time. An exact solution for time-dependent OH evolution under pseudo- first-order chain reaction conditions is presented, which correctly predicts new kinetic signatures not included in previous OH + 03 kinetic analyses. Specifically, the solutions predict an initial exponential loss (chain "induction") of the OH radical to a steady-state level ([OH](sub ss)), with this fast initial decay determined by the sum of both chain rate constants, k(sub ind) = k(sub 1) + k(sub 2). By monitoring the chain induction feature, this sum of the rate constants is determined to be k(sub ind) = 8.4(8) x 10(exp -14) cu cm/molecule/s for room temperature reagents. This is significantly higher than the values currently recommended for use in atmospheric models, but in excellent agreement with previous results from Ravishankara et al.

Kulcke, Axel↗