Search NASA⌕ Search

SEARCH · Search NASA

Results for “kinetic equation”

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 415 records · Page 23

Multifluids description of dynamics of upper atmosphere

A multifluids model to investigate ionospheric dynamics was established on kinetic theory. Its resultant equations are used to examine the following dynamic problems in the gamma region of 80-2000 Km of the ionosphere: (1) propagation of acoustic modes in the 500-2,000 Km of the ionosphere (two fluid model); (2) the relation between the cross field plasma drift instabilities and type I and type II ionospheric irregularities; and (3) time dependent neutral wind structure and horizontal pressure gradient.

Wu, S. T.↗

The NASA Ames Research Center one- and two-dimensional stratospheric models. Part 1: The one-dimensional model

A one-dimensional model of stratospheric trace constituents is described in detail. Specifically, the numerical solution of the species continuity equations, including a technique for treating the stiff differential equations representing the chemical kinetic terms, and an appropriate method for simulating the diurnal variations of the species concentrations, are discussed. A specialized treatment of atmospheric photodissociation rates is outlined in the text. The choice of a vertical eddy diffusivity profile and its success in predicting the vertical tracer distributions (carbon 14, methane, and nitrous oxide) are also discussed.

Turco, R. P.↗

Continuum models for static and dynamic analysis of repetitive lattices

A simple, rational approach is presented for developing continuum models for large repetitive lattice structures subjected to static and dynamic loadings. The procedure involves introducing kinematic assumptions to reduce the dimensionality of the lattice, and equating the strain and kinetic energies of the continuum model to those of the original lattice structure. The proposed procedure is applied to obtain effective elastic and dynamic characteristics of continuum plate models for double-layered grids. Numerical results are presented of stress and free vibration problems for two double-layered grids. These problems demonstrate the high accuracy of the continuum plate models developed. Also, an assessment is made of different approximations to these models.

Noor, A. K.↗

High-altitude gas releases - Transition from collisionless flow to diffusive flow in a nonuniform atmosphere

The paper is motivated by a need for describing high-altitude gas releases from rockets or the Space Shuttle. In the tenuous upper atmosphere, the injected gases expand from a collisionless to a collision-dominated state. The analysis presented extends the efforts of Baum (1973-1974) to include the effects of a nonuniform background atmosphere and a velocity-dependent collision frequency. The unsteady expansion of gas releases is analyzed using gas kinetic theory. The Boltzmann equation with the Krook collision integral is solved numerically. At late times (after many collisions with the background atmosphere) the solution is identical to the one given by diffusion in a nonuniform atmosphere (Bernhardt, 1979). The theoretical model predicts elongation and heating of the vapor trails due to collisions.

Bernhardt, P. A.↗

The molten carbonate carbon dioxide concentrator - Cathode performance at high CO2 utilization

An experimental study of the cathode performance of the molten carbonate carbon dioxide concentrator (MCCDC) at typical gas concentrations is reported, and the behavior is described in terms of standard electrochemical kinetic analysis. Theoretical aspects of the MCCDC are discussed, including the current-generating processes, the overpotential, and the reaction kinetics. A final working equation is derived, and the electrochemical cell is described. The parameters determined by the experiments include cathode overpotentials and removal and current efficiency as a function of current density, as well as carbon dioxide removal rate as a function of flow rate. The results are compared with data for other cells; the agreement with data for the molten carbonate fuel cell suggests that the major difference between the two cells is in the wetting characteristics. The MCCDC achieves high removal efficiencies at high current efficiencies. The lowest current efficiency recorded was 80 percent, and at most current densities, current efficiencies were above 100 percent.

Weaver, J. L.↗

A finite element computation of turbulent boundary layer flows with an algebraic stress turbulence model

An algebraic stress turbulence model and a computational procedure for turbulent boundary layer flows which is based on the semidiscrete Galerkin FEM are discussed. In the algebraic stress turbulence model, the eddy viscosity expression is obtained from the Reynolds stress turbulence model, and the turbulent kinetic energy dissipation rate equation is improved by including a production range time scale. Good agreement with experimental data is found for the examples of a fully developed channel flow, a fully developed pipe flow, a flat plate boundary layer flow, a plane jet exhausting into a moving stream, a circular jet exhausting into a moving stream, and a wall jet flow.

Kim, Sang-Wook↗

Kinetic theory for electrostatic waves due to transverse velocity shears

A kinetic theory in the form of an integral equation is provided to study the electrostatic oscillations in a collisionless plasma immersed in a uniform magnetic field and a nonuniform transverse electric field. In the low temperature limit the dispersion differential equation is recovered for the transverse Kelvin-Helmholtz modes for arbitrary values of K parallel, where K parallel is the component of the wave vector in the direction of the external magnetic field assumed in the z direction. For higher temperatures the ion-cyclotron-like modes described earlier in the literature by Ganguli, Lee and Plamadesso are recovered. In this article, the integral equation is reduced to a second-order differential equation and a study is made of the kinetic Kelvin-Helmholtz and ion-cyclotron-like modes that constitute the two branches of oscillation in a magnetized plasma including a transverse inhomogeneous dc electric field.

Ganguli, G.↗

The oblique behavior of low-frequency electromagnetic waves excited by newborn cometary ions

The free energy in oxygen or hydrogen ions freshly created in the solar wind stimulates low-frequency electromagnetic waves whose growth does not always maximize at parallel propagation. Exploration of the wave vector plane discloses the frequent occurrence of islets of oblique growth unconnected to the unstable parallel modes. Contour plots of the growth rate, real frequency, polarization, and magnetic compression characterize the oblique wave behavior for large values of the initial pitch angle of the cometary particles. Although wave-particle (Landau and cyclotron) resonances feed most of the surveyed oblique instabilities, some are seemingly fluidlike. The results, obtained from the numerical solution of the kinetic dispersion and wave equations, imply that newborn ions can easily excite significant oblique hydromagnetic wave activity. Cometary environments provide the adopted plasma model, but the study is helpful in the interpretation of other low-frequency wave observations in space.

Brinca, Armando L.↗

A review of near-wall Reynolds-stress

The advances made in second-order near-wall turbulence closures are summarized. All closures examined are based on some form of high Reynolds number models for the Reynolds stress and the turbulent kinetic energy dissipation rate equations. Consequently, most near-wall closures proposed to data attempt to modify the high Reynolds number models for the dissipation rate equation so that the resultant models are applicable all the way to the wall. The near-wall closures are examined for their asymptotic behavior so that they can be compared with the proper near-wall behavior of the exact equations. A comparison of the closure's performance in the calculation of a low Reynolds number plane channel flow is carried out. In addition, the closures are evaluated for their ability to predict the turbulence statistics and the limiting behavior of the structure parameters compared to direct simulation data.

So, R. M. C.↗

The saturation of the electron-cyclotron maser instability and the interpretation of solar millisecond spikes

A self-consistent numeric two-dimensional code of the kinetic wave-particle equations developed to investigate the maser dynamics in the solar context is applied to solar millisecond-spike observations in order to improve the diagnostic capabilities of the theory of the electron-cyclotron maser instablitity. Attention is given to the inhomogeneity of the magnetic field selecting magneto-ionic modes with relatively short saturation lengths and suppressing mechanisms such as collisional deflection, free-free absorption, and gyroresonance absorption. The time scales of maser saturation in respect to time scales of global particle changes in a magnetic loop are covered, relevant observations of solar millisecond spikes are described, and the interpretation in terms of physical parameters deduced from the quasi-linear maser simulations are presented. It is demonstrated that the quasi-linear simulations make it possible to constrain the physical parameters from the observed time scale and frequency.

Aschwanden, M. J.↗

Discontinuous Galerkin Methods for NonLinear Differential Systems

This talk considers simplified finite element discretization techniques for first-order systems of conservation laws equipped with a convex (entropy) extension. Using newly developed techniques in entropy symmetrization theory, simplified forms of the discontinuous Galerkin (DG) finite element method have been developed and analyzed. The use of symmetrization variables yields numerical schemes which inherit global entropy stability properties of the PDE (partial differential equation) system. Central to the development of the simplified DG methods is the Eigenvalue Scaling Theorem which characterizes right symmetrizers of an arbitrary first-order hyperbolic system in terms of scaled eigenvectors of the corresponding flux Jacobian matrices. A constructive proof is provided for the Eigenvalue Scaling Theorem with detailed consideration given to the Euler equations of gas dynamics and extended conservation law systems derivable as moments of the Boltzmann equation. Using results from kinetic Boltzmann moment closure theory, we then derive and prove energy stability for several approximate DG fluxes which have practical and theoretical merit.

Barth, Timothy↗

Budgets of divergent and rotational kinetic energy during two periods of intense convection

The derivations of the energy budget equations for divergent and rotational components of kinetic energy are provided. The intense convection periods studied are: (1) synoptic scale data of 3 or 6 hour intervals and (2) mesoalphascale data every 3 hours. Composite energies and averaged budgets for the periods are presented; the effects of random data errors on derived energy parameters is investigated. The divergent kinetic energy and rotational kinetic energy budgets are compared; good correlation of the data is observed. The kinetic energies and budget terms increase with convective development; however, the conversion of the divergent and rotational energies are opposite.

Buechler, D. E.↗

Development and Documentation of kL-Based Linear, Nonlinear, and Full Reynolds Stress Turbulence Models

The development and implementation of kL-based Reynolds Average Navier-Stokes (RANS) turbulence models are reported herein. The kL is based on Abdol-Hamid's closure and Menter's modi cation to Rotta's two-equation model. Rotta shows that a reliable transport equation can be formed from the turbulent length scale L, and the turbulent kinetic energy k. Rotta's kL equation is well suited for term-by-term modeling and displays useful features compared to other scale formulation. One of the important di erences is the inclusion of higher order velocity derivatives in the source terms of the scale equation. This can enhance the ability of RANS solvers to simulate unsteady ows in URANS mode. The present report documents the formulation of three model levels of turbulence models as implemented in the CFD code FUN3D. Methodology and calibration examples are shown in detail. The levels are the linear k-kL and the two-equation algebraic Reynolds stress model (ARSM) as well as the full Reynolds Stress Model (RSM). Attached, separated and corner ow cases are documented and compared with experimental, theoretical and other turbulence model data. The results show generally very good comparisons with canonical and experimental data. The results from this formulation are similar or better than results using the SST two- equation turbulence model. ARSM shows great promise with similar level of computational resources as general two equation turbulence models.

Abdol-Hamid, Khaled S.↗

kL-based Linear and Nonlinear Two-Equation Turbulence Models

The development and implementation of kL-based Reynolds Average Navier-Stokes (RANS) two-equation turbulence models are reported herein. The kL is based on Abdol-Hamid's closure and Menter's modification to Rotta's two-equation model. Rotta shows that a reliable transport equation can be formed from the turbulent length scale L, and the turbulent kinetic energy k. Rotta's kL equation is well suited for term-by-term modeling and displays useful features compared to other scale formulation. One of the important differences is the inclusion of higher order velocity derivatives in the source terms of the scale equation. This can enhance the ability of RANS solvers to simulate unsteady flows in URANS mode. The present report documents the formulation of two model levels of turbulence models as implemented in the computational fluid dynamics FUN3D code. The levels are the two-equation linear k-kL and the two-equation algebraic Reynolds stress model (ARSM). Free shear, separated and corner flow cases are documented and compared with experimental, and other turbulence model data. The results show generally very good comparisons with experimental data. The results from this formulation are similar or better than results using the SST two-equation turbulence model. ARSM shows great promise with similar level of computational resources as basic two-equation turbulence models.

Abdol-Hamid, K. S.↗

Juncture Flow Computations using kL-Based Turbulence Models

The development and implementation of kL-based Reynolds-Average Navier-Stokes (RANS) two-equation turbulence models are reported in this paper. The kL model is based on Abdol-Hamid's closure and Menter's modification to Rotta's two-equation model. Rotta showed that a reliable transport equation can be formed from the turbulence length scale L, and the turbulence kinetic energy k. Rotta's kL equation is well suited for a term-by-term modeling and displays useful features compared to other scale formulations. One of the important differences is the inclusion of higher-order velocity derivatives in the source terms of the scale equation. This can enhance the ability of RANS solvers to simulate unsteady flows in URANS mode. The k-kL scheme has been modified to include the Algebraic Reynolds Stress Model (ARSM) and the Quadratic Constitutive Relation (QCR) as nonlinear models. These models have gone through extensive validations using two-dimensional benchmark problems. The present study documents the application of the k-kL schemes to simulate flow around the Juncture Flow Model (JFM). The k-kL prediction results show generally good comparisons with measurements. The results from this formulation are similar to, or better than results using the SA-RC-QCR2000 two-equation turbulence model. The ARSM and the QCR formulations of k-kL show promise with a similar level of computational resources as basic two-equation turbulence models.

K S Abdol-Hamid↗

Investigation of the hydrochlorination of SiCl4

The hyrochlorination of SiC14 and m.g. silicon metal to produce SiHC13, was investigated. Reaction kinetic measurements were carried out to collect additional rate data at 525 C and 550 C. A theoretical study was carried out to provide a kinetic model and a rate equation for the hydrochlorination reaction. Results of this preliminary study show that the rate of formation of SiHC13 follows a pseudo first order kinetics. The rate constants were measured at three temperatures, 550 C, 500 C and 450 C, respectively. The activation energy was determined from the Arrhenius plot to give a value of 13.2 Kcal/mole. The design of a quartz reactor to measure reaction rates and equilibrium conversion of SiHC13 at reaction temperature up to 650 C was completed.

Mui, J. Y. P.↗

Implementation of a Transition Model in a NASA Code and Validation Using Heat Transfer Data on a Turbine Blade

The purpose of this report is to summarize and document the work done to enable a NASA CFD code to model laminar-turbulent transition process on an isolated turbine blade. The ultimate purpose of the present work is to down-select a transition model that would allow the flow simulation of a variable speed power turbine to be accurately performed. The flow modeling in its final form will account for the blade row interactions and their effects on transition which would lead to accurate accounting for losses. The present work only concerns itself with steady flows of variable inlet turbulence. The low Reynolds number k- model of Wilcox and a modified version of the same model will be used for modeling of transition on experimentally measured blade pressure and heat transfer. It will be shown that the k- model and its modified variant fail to simulate the transition with any degree of accuracy. A case is thus made for the adoption of more accurate transition models. Three-equation models based on the work of Mayle on Laminar Kinetic Energy were explored. The three-equation model of Walters and Leylek was thought to be in a relatively mature state of development and was implemented in the Glenn-HT code. Two-dimensional heat transfer predictions of flat plate flow and two-dimensional and three-dimensional heat transfer predictions on a turbine blade were performed and reported herein. Surface heat transfer rate serves as sensitive indicator of transition. With the newly implemented model, it was shown that the simulation of transition process is much improved over the baseline k- model for the single Reynolds number and pressure ratio attempted; while agreement with heat transfer data became more satisfactory. Armed with the new transition model, total-pressure losses of computed three-dimensional flow of E3 tip section cascade were compared to the experimental data for a range of incidence angles. The results obtained, form a partial loss bucket for the chosen blade. In time the loss bucket will be populated with losses at additional incidences. Results obtained thus far will be discussed herein.

Ameri, Ali A.↗