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 235 records · Page 13

Integral boundary-layer models for turbulent separated flows

Two-dimensional, integral boundary-layer methods of the entrainment and dissipation type, with closure relationships which account for turbulent separated flows through the velocity and turbulence models, are developed. Kinematic shape factor relations used in the two methods are evaluated by comparison with recently obtained experimental data for separated flows. Prescribed-displacement-thickness computations are made with both methods and compared to experiment for a series of separated flows, including low speed and transonic shock-induced separation. Effects of upstream history on the turbulent shear stress development are shown to be important for the accurate prediction of the separated flows considered and the experimental trends are reasonably modeled by a rate equation derived from the turbulent kinetic energy equation.

Thomas, J. L.↗

Unsteady viscous-inviscid interaction procedures for transonic airfoils using Cartesian grids

A viscous-inviscid interaction procedure for transonic airfoils using an Euler/integral boundary layer formulation and Cartesian grids is presented. The approach is based on a time dependent formulation for both the integral boundary layer equations and the Euler equations. Effects of upstream history on the shear stress are modeled by a time dependent rate equation derived from the turbulent kinetic energy equation. Results are presented for two of the test cases reported by Cook et al. (1979) for the RAE 2822 supercritical airfoil and one of the cases reported by Harris (1981) for the NACA 0012 symmetric airfoil. In general, the results are in good agreement with experiment.

Fenno, Charles C., Jr.↗

Numerical prediction of turbulent flow over airfoil sections with a new nonequilibrium turbulence model

A new nonequilibrium turbulence closure model has been developed for computing wall bounded two-dimensional turbulent flows. This two-layer eddy viscosity model was motivated by the success of the Johnson-King model in separated flow regions. The influence of history effects are described by an ordinary differential equation developed from the turbulent kinetic energy equation. The performance of the present model has been evaluated by solving the flow around three airfoils using the Reynolds time-averaged Navier-Stokes equations. Excellent results were obtained for both attached and separated turbulent flows about the NACA 0012 airfoil, the RAE 2822 airfoil, and the Integrated Technology A 153W airfoil. Based on the comparison of the numerical solutions with the available experimental data, it is concluded that the new nonequilibrium turbulence model accurately captures the history effects of convection and diffusion on turbulence.

Ahmed, S.↗

Numerical computation of viscous flow about unconventional airfoil shapes

A new two-dimensional computer code was developed to analyze the viscous flow around unconventional airfoils at various Mach numbers and angles of attack. The Navier-Stokes equations are solved using an implicit, upwind, finite-volume scheme. Both laminar and turbulent flows can be computed. A new nonequilibrium turbulence closure model was developed for computing turbulent flows. This two-layer eddy viscosity model was motivated by the success of the Johnson-King model in separated flow regions. The influence of history effects are described by an ordinary differential equation developed from the turbulent kinetic energy equation. The performance of the present code was evaluated by solving the flow around three airfoils using the Reynolds time-averaged Navier-Stokes equations. Excellent results were obtained for both attached and separated flows about the NACA 0012 airfoil, the RAE 2822 airfoil, and the Integrated Technology A 153W airfoil. Based on the comparison of the numerical solutions with the available experimental data, it is concluded that the present code in conjunction with the new nonequilibrium turbulence model gives excellent results.

Ahmed, S.↗

Modeling Disturbance Dynamics in Transitional and Turbulent Boundary Layers

The dynamics of an ensemble of linear disturbances in boundary-layer flows at various Reynolds numbers is studied through an analysis of the transport equations for the mean disturbance kinetic energy and energy dissipation rate. Effects of adverse and favorable pressure-gradients on the disturbance dynamics are also included in the analysis. Unlike the fully turbulent regime where nonlinear phase scrambling of the fluctuations affects the flow field even in proximity to the wall, the early stage transition regime fluctuations studied here are influenced across the boundary layer by the solid boundary. In addition, the dominating dynamics in the disturbance kinetic energy equation is governed by the energy production, pressure-transport and viscous diffusion - also in contrast to the fully turbulent regime. For the disturbance dissipation rate, a dynamic balance exists between the destruction and diffusion of dissipation.

Grosch, C. E.↗

Comparison of numerical techniques for integration of stiff ordinary differential equations arising in combustion chemistry

The efficiency and accuracy of several algorithms recently developed for the efficient numerical integration of stiff ordinary differential equations are compared. 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 interactive solution of the algebraic energy conservation equation to compute the temperature does not result in significant errors. In addition, this method is more efficient than evaluating the temperature by integrating its time derivative. Significant reductions in computational work are realized by updating the rate constants (k = at(supra N) N exp(-E/RT) only when the temperature change exceeds an amount delta T that is problem dependent. An approximate expression for the automatic evaluation of delta T is derived and is shown to result in increased efficiency.

Radhakrishnan, K.↗

Theory and modeling of atmospheric turbulence, part 1

The cascade transfer which is the only function to describe the mode coupling as the result of the nonlinear hydrodynamic state of turbulence is discussed. A kinetic theory combined with a scaling procedure was developed. The transfer function governs the non-linear mode coupling in strong turbulence. The master equation is consistent with the hydrodynamical system that describes the microdynamic state of turbulence and has the advantages to be homogeneous and have fewer nonlinear terms. The modes are scaled into groups to decipher the governing transport processes and statistical characteristics. An equation of vorticity transport describes the microdynamic state of two dimensional, isotropic and homogeneous, geostrophic turbulence. The equation of evolution of the macrovorticity is derived from group scaling in the form of the Fokker-Planck equation with memory. The microdynamic state of turbulence is transformed into the Liouville equation to derive the kinetic equation of the singlet distribution in turbulence. The collision integral contains a memory, which is analyzed with pair collision and the multiple collision. Two other kinetic equations are developed in parallel for the propagator and the transition probability for the interaction among the groups.

Source record↗

Filtering analysis of a direct numerical simulation of the turbulent Rayleigh-Benard problem

A filtering analysis of a turbulent flow was developed which provides details of the path of the kinetic energy of the flow from its creation via thermal production to its dissipation. A low-pass spatial filter is used to split the velocity and the temperature field into a filtered component (composed mainly of scales larger than a specific size, nominally the filter width) and a fluctuation component (scales smaller than a specific size). Variables derived from these fields can fall into one of the above two ranges or be composed of a mixture of scales dominated by scales near the specific size. The filter is used to split the kinetic energy equation into three equations corresponding to the three scale ranges described above. The data from a direct simulation of the Rayleigh-Benard problem for conditions where the flow is turbulent are used to calculate the individual terms in the three kinetic energy equations. This is done for a range of filter widths. These results are used to study the spatial location and the scale range of the thermal energy production, the cascading of kinetic energy, the diffusion of kinetic energy, and the energy dissipation. These results are used to evaluate two subgrid models typically used in large-eddy simulations of turbulence. Subgrid models attempt to model the energy below the filter width that is removed by a low-pass filter.

Eidson, T. M.↗

A new turbulence closure model for boundary layer flows with strong adverse pressure gradients and separation

A new turbulence closure model designed specifically to treat two-dimensional, turbulent boundary layers with strong adverse pressure gradients and attendant separation, is presented. The influence of history effects are modeled by using an ordinary differential equation (ODE) derived from the turbulence kinetic-energy equation, to describe the streamwise development of the maximum Reynolds shear stress in conjunction with an assumed eddy-viscosity distribution which has as its velocity scale the maximum Reynolds shear stress. In the outer part of the boundary layer, the eddy viscosity is treated as a free parameter which is adjusted in order to satisfy the ODE for the maximum shear stress. Because of this, the model s not simply an eddy-viscosity model, but contains features of a Reynolds-stress model. Comparisons with experiments are presented which clearly show the proposed model to be superior to the Cebeci-Smith model in treating strongly retarded and separated flows. In contrast to two-equation, eddy-viscosity models, it requires only slightly more computational effort than simple models like the Cebeci-Smith model.

Johnson, D. A.↗

Computational Fluid Mechanics

A one-equation turbulence model based on the turbulent kinetic energy equation is presented. The model is motivated by the success of the Johnson-King model and incorporates a number of features uncovered by Simpson's experiments on separated flows. Based on the results obtained, the model duplicates the success of algebraic models in attached flow regions and outperforms the two-equation models in detached flow regions.

Hassan, H. A.↗

Fast algorithms for combustion kinetics calculations: A comparison

To identify the fastest algorithm currently available for the numerical integration of chemical kinetic rate equations, several algorithms were examined. Findings to date are summarized. The algorithms examined include two general-purpose codes EPISODE and LSODE and three special-purpose (for chemical kinetic calculations) codes CHEMEQ, CRK1D, and GCKP84. In addition, an explicit Runge-Kutta-Merson differential equation solver (IMSL Routine DASCRU) is used to illustrate the problems associated with integrating chemical kinetic rate equations by a classical method. Algorithms were applied to two test problems drawn from combustion kinetics. These problems included all three combustion regimes: induction, heat release and equilibration. Variations of the temperature and species mole fraction are given with time for test problems 1 and 2, respectively. Both test problems were integrated over a time interval of 1 ms in order to obtain near-equilibration of all species and temperature. Of the codes examined in this study, only CREK1D and GCDP84 were written explicitly for integrating exothermic, non-isothermal combustion rate equations. These therefore have built-in procedures for calculating the temperature.

Radhakrishnan, K.↗

Verification of theory on weak turbulence relating to the sequence of diffuse plasma resonances in space.

An interpretation of the sequence of diffuse plasma resonances observed by space probes (Alouette 2 and ISIS-1 satellites) is developed in terms of wave-particle nonlinear interaction in a weakly turbulent plasma including the electrostatic electron cyclotron harmonic wave instability. The longest time duration of the center frequency of the diffuse plasma resonance is found to coincide with the most favorable condition for the electrostatic electron cyclotron harmonic wave instability which is obtained by solving the dispersion equation obtained for a linear approximation of the kinetic wave equation for the warm magnetoactive plasma. The electrostatic field due to the transmission of the intense rf pulse produces plasma turbulence involving nonlinear wave-wave interaction and temperature anisotropy which leads to instability. This instability supplies energy to the turbulence. The process can be thought of as a feedback system.

Oya, H.↗

Spectral study of wintertime kinetic energy of the Northern Hemisphere in the troposphere

Characteristics of the kinetic energy of wind fields at various pressure levels were analyzed, and significant wavenumbers in the wavenumber-frequency domain were identified. The nonlinear interaction terms of the kinetic energy equation were examined, and the distribution of the kinetic energy at the 850 mb, 500 mb, and 200 mb levels was calculated. A 5 deg latitude-longitude square grid was used, with NMC data for the 1975-1976 winter in the 20-60 deg N at 500 mb and 20-85 deg N for the 200 mb and 850 mb levels. The kinetic energy distribution was determined to be geography-dependent, with wavenumbers 6-9 westerly waves in the midfrequency range contributing significantly to kinetic energy maxima over the North Pacific and the east coast of North America. The contribution of the nonlinear interactions of these waves, which correspond to the longitudinal convergence of the kinetic energy flux, was found to be larger than the meridional convergence of the kinetic energy flux, and to occur mainly between 30-50 deg N. The nonlinear interactions were a negative contribution over the North Pacific at the 200 mb level.

Lee, H. N.↗

Analysis of Three-Dimensional, Nonlinear Development of Wave-Like Structure in a Compressible Round Jet

An analysis of the nonlinear development of the large-scale structures or instability waves in compressible round jets was conducted using the integral energy method. The equations of motion were decomposed into two sets of equations; one set governing the mean flow motion and the other set governing the large-scale structure motion. The equations in each set were then combined to derive kinetic energy equations that were integrated in the radial direction across the jet after the boundary-layer approximations were applied. Following the application of further assumptions regarding the radial shape of the mean flow and the large structures, equations were derived that govern the nonlinear, streamwise development of the large structures. Using numerically generated mean flows, calculations show the energy exchanges and the effects of the initial amplitude on the coherent structure development in the jet.

Dahl, Milo D.↗

Reduced Equations for Calculating the Combustion Rates of Jet-A and Methane Fuel

Simplified kinetic schemes for Jet-A and methane fuels were developed to be used in numerical combustion codes, such as the National Combustor Code (NCC) that is being developed at Glenn. These kinetic schemes presented here result in a correlation that gives the chemical kinetic time as a function of initial overall cell fuel/air ratio, pressure, and temperature. The correlations would then be used with the turbulent mixing times to determine the limiting properties and progress of the reaction. A similar correlation was also developed using data from NASA's Chemical Equilibrium Applications (CEA) code to determine the equilibrium concentration of carbon monoxide as a function of fuel air ratio, pressure, and temperature. The NASA Glenn GLSENS kinetics code calculates the reaction rates and rate constants for each species in a kinetic scheme for finite kinetic rates. These reaction rates and the values obtained from the equilibrium correlations were then used to calculate the necessary chemical kinetic times. Chemical kinetic time equations for fuel, carbon monoxide, and NOx were obtained for both Jet-A fuel and methane.

Molnar, Melissa↗

A method for calculating unsteady turbulent boundary layers in two- and three-dimensional flows.

The governing equations for an unsteady turbulent boundary layer on a swept infinite cylinder, composed of a continuity equation, a pair of momentum equations and a pair of rate equations for the shear stress, based on the turbulent kinetic-energy equation, are solved numerically. These rate equations for the shear-stress express a balance between the convection, production, dissipation and diffusion of the turbulent shear stress components. Due to the physical model of turbulence employed, this system of equations is hyperbolic. An explicit, second-order accurate, conditionally stable difference scheme for this system, together with appropriate boundary and initial conditions, is formulated and developed into a computer program. Calculations for oscillating free stream flows with no pressure gradient show significant unsteady effects on the turbulent boundary layer.

Singleton, R. E.↗

Assessment of a transitional boundary layer theory at low hypersonic Mach numbers

An investigation was carried out to assess the accuracy of a transitional boundary layer theory in the low hypersonic Mach number regime. The theory is based upon the simultaneous numerical solution of the boundary layer partial differential equations for the mean motion and an integral form of the turbulence kinetic energy equation which controls the magnitude and development of the Reynolds stress. Comparisions with experimental data show the theory is capable of accurately predicting heat transfer and velocity profiles through the transitional regime and correctly predicts the effects of Mach number and wall cooling on transition Reynolds number. The procedure shows promise of predicting the initiation of transition for given free stream disturbance levels. The effects on transition predictions of the pressure dilitation term and of direct absorption of acoustic energy by the boundary layer were evaluated.

Shamroth, S. J.↗