SEARCH · Engineering Papers

Results for “Turbulence Stabilization”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 181 records · Page 10

All sources NASA NTRS DOE OSTI

NASA NTRS1957

Base Pressure at Supersonic Speeds on Two-dimensional Airfoils and on Bodies of Revolution with and Without Fins Having Turbulent Boundary Layers

An analysis has been made of available experimental data to show the effects of most of the variables that are more predominant in determining base pressure at supersonic speeds. The analysis covers base pressures for two-dimensional airfoils and for bodies of revolution with and without stabilizing fins and is restricted to turbulent boundary layers. The present status of available experimental information is summarized as are the existing methods for predicting base pressure. A simple semiempirical method is presented for estimating base pressure. For two-dimensional bases, this method stems from an analogy established between the base-pressure phenomena and the peak pressure rise associated with the separation of the boundary layer. An analysis made for axially symmetric flow indicates that the base pressure for bodies of revolution is subject to the same analogy. Based upon the methods presented, estimations are made of such effects as Mach number, angle of attack, boattailing, fineness ratio, and fins. These estimations give fair predictions of experimental results. (author)

LOVE EUGENE S↗

NASA NTRS1987

Prediction of transition on a flat plate under the influence of free-stream turbulence using low-Reynolds-number two-equation turbulence models

The prediction characteristics and capabilities of two popular two-equation low-Reynolds-number turbulence models (Jones, Launder 1972; Lam, Bremhorst 1981) have been evaluated with respect to the prediction of transition of a flat plate under the influence of free-stream turbulence. The sensitivity of the predictions to free-stream turbulence intensity, initial starting location of the calculation, and the assumed initial starting profiles for k and epsilon has been determined and presented. Although both models predict the correct qualitative characteristics of transition, they also exhibit significant quantitative deficiencies with regard to both the predicted location and the length of transition. A modification to the production term in the turbulent kinetic energy equation is proposed which is based on a simple stability criterion and correlated to the free-stream turbulence level. The modification becomes inactive in the fully turbulent regime, but is shown to improve both the qualitative and quantitative characteristics of the transition predictions.

Schmidt, Rodney C.↗

NASA NTRS1980

A piloted simulator investigation of static stability and stability/control augmentation effects on helicopter handling qualities for instrument approach

A ground simulator experiment was conducted on the Flight Simulator for Advanced Aircraft at Ames Research Center to investigate the influence of several static stability and stability/control augmentation design parameters on helicopter flying qualities during terminal area operations in instrument conditions. Effects of light turbulence were included. Two levels of static stability in each rotational axis (pitch, roll, yaw) were examined for a hingeless rotor configuration. The variations in pitch and roll were: (1) stable and (2) neutral static stability; in yaw there were two stable levels. Four types of stability/control augmentation were also examined for the lower level of static stability in each axis. This latter investigation covered three helicopter rotor types: hingeless, articulated, and teetering. Four pilots performed a total of 105 evaluations of these parameters for a representative VOR instrument approach task. Pilot rating results indicate the acceptability of neutral static stability longitudinally and laterally and the need for pitch-roll attitude augmentation to achieve a satisfactory system.

Lebacqz, J. V.↗

NASA NTRS1986

Jeans criterion in a turbulent medium

According to the classical Jeans analysis, all the molecular clouds of mass larger than a few 100 M(solar), size larger than about 1pc and kinetic temperature Tk less than 30K are gravitationally unstable. We have shown that in clouds supported by internal supersonic motions, local gravitational instabilities may appear within molecular clouds which are globally stable. The argument is threefold: (1) when the turbulent kinetic energy is included into the internal energy term, the virial equilibrium condition shows that molecular clouds such as those observed, which are gravitationally unstable according to the Jeans criterion, are indeed globally stable if supported by a turbulent velocity field of power spectrum steeper than 3; (2) 2D compressible hydrodynamical simulations show that a supersonic turbulent velocity field generates a turbulent pressure within clouds, the gradients of which stabilize the unstable scales (i.e., the largest scales and the cloud itself) against gravitational collapse; (3) an analysis similar to the Jeans approach but including the turbulent pressure gradient term, gives basically the same results as those given in (1). Clouds of mean density lower than a critical value are found to be stable even though more massive than their Jeans mass. In clouds of mean density larger than that critical value, the gravitational instability appears only over a range of scales smaller than the cloud size, the largest scales being stable. In practice, the observed mean densities are lower than this critical value: the observation of a small number of cores and stars of a few solar masses embedded in clouds of several hundred solar masses can only be understood in terms of small scale density fluctuations of large amplitude generated by the supersonic turbulence which would occasionally overtake the limit of gravitational stability.

Bonazzola, S.↗

NASA NTRS1954

Wind-tunnel Tests of a Model of a Wingless Fin-controlled Missile to Obtain Static Stability and Control Characteristics Through a Range of Mach Numbers from 0.5 to 0.88

An investigation at medium to high subsonic speeds has been conducted in the Langley low-turbulence pressure tunnel to determine the static stability and control characteristics and to measure the fin normal forces and moments for a model of a wingless fin-controlled missile. The data were obtained at Reynolds number of 2.1 x 10(6) based on the missile maximum diameter or 17.7 x 10(6) based on missile length; this Reynolds number was found to be large enough to avoid any large scale effects between the test and the expected flight Reynolds number. With the horizontal-fin deflection limited to a maximum of 6 degrees, longitudinally stable and trimmed flight could not be maintained beyond an angle of attack of 17 degrees for a Mach number of 0.88 and beyond 20 degrees for a Mach number of 0.50 for any center-of-gravity location without the use of some auxiliary stability or control device such as jet vanes. Mach number had no appreciable effect on the center-of-pressure positions and only a slight effect on neutral-point position. There was a shift in neutral-point position of about 1 caliber as the angle of attack was varied through the range for which the neutral point could be determined. Yawing the model to angles of sideslip up to 7 degrees had little effect on the longitudinal stability at angles of attack up to 15 degrees; however, above 15 degrees, the effect of sideslip was destabilizing. With the vertical fins at a plus-or-minus 6 degree roll deflection, the rolling moment caused by yawing the model at high angles of attack could be trimmed out up to angles of sideslip of 6.5 degrees and an angle of attack of 26 degrees for a Mach number of 0.50; this range of sideslip angles was reduced to 3 degrees at a Mach number of 0.88. The data indicated that, at lower angles of attack, the trim range extended to higher angles of sideslip. The total normal-force and hinge-moment coefficients for both horizontal fins were slightly nonlinear with both angle-of-attack and fin deflection. The effect of Mach number was to reduce the slopes of the hinge-moment coefficient with angle of attack and deflection angle. In general, the effort of increasing the sideslip angle was to reduce the values of the fin normal-force and hinge-moment coefficients.

Burrows, Dale L↗

NASA NTRS1994

Modification of the Two-equation Turbulence Model in NPARC to a Chien Low Reynolds Number K-epsilon Formulation

This report documents the changes that were made to the two-equation k-epsilon turbulence model in the NPARC (National-PARC) code. The previous model based on the low Reynolds number model of Speziale, was replaced with the low Reynolds number k-epsilon model of Chien. The most significant difference was in the turbulent Prandtl numbers appearing in the diffusion terms of the k and epsilon transport equations. A new inflow boundary condition and stability enhancements were also implemented into the turbulence model within NPARC. The report provides the rationale for making the change to the Chien model, code modifications required, and comparisons of the performances of the new model with the previous k-epsilon model and algebraic models used most often in PARC/NPARC. The comparisons show that the Chien k-epsilon model installed here improves the capability of NPARC to calculate turbulent flows.

Georgiadis, Nicholas J.↗

NASA NTRS1997

Modeling of Turbulent Swirling Flows

Aircraft engine combustors generally involve turbulent swirling flows in order to enhance fuel-air mixing and flame stabilization. It has long been recognized that eddy viscosity turbulence models are unable to appropriately model swirling flows. Therefore, it has been suggested that, for the modeling of these flows, a second order closure scheme should be considered because of its ability in the modeling of rotational and curvature effects. However, this scheme will require solution of many complicated second moment transport equations (six Reynolds stresses plus other scalar fluxes and variances), which is a difficult task for any CFD implementations. Also, this scheme will require a large amount of computer resources for a general combustor swirling flow. This report is devoted to the development of a cubic Reynolds stress-strain model for turbulent swirling flows, and was inspired by the work of Launder's group at UMIST. Using this type of model, one only needs to solve two turbulence equations, one for the turbulent kinetic energy k and the other for the dissipation rate epsilon. The cubic model developed in this report is based on a general Reynolds stress-strain relationship. Two flows have been chosen for model evaluation. One is a fully developed rotating pipe flow, and the other is a more complex flow with swirl and recirculation.

Shih, Tsan-Hsing↗

NASA NTRS2000

The Surface Pressure Response of a NACA 0015 Airfoil Immersed in Grid Turbulence: Characteristics of the Turbulence - Volume 1

Two grids have been developed for the Virginia Tech 6 ft x 6 ft Stability wind tunnel for the purpose of generating homogeneous isotropic turbulent flows for the study of unsteady airfoil response. The first, a square bi-planar grid with a 12" mesh size and an open area ratio of 69.4%, was mounted in the wind tunnel contraction. The second grid, a metal weave with a 1.2 in. mesh size and an open area ratio of 68.2% was mounted in the tunnel test section. Detailed statistical and spectral measurements of the turbulence generated by the two grids are presented for wind tunnel free stream speeds of 10, 20, 30 and 40 m/s. These measurements show the flows to be closely homogeneous and isotropic. Both grids produce flows with a turbulence intensity of about 4% at the location planned for the airfoil leading edge. Turbulence produced by the large grid has an integral scale of some 3.2 inches here. Turbulence produced by the small grid is an order of magnitude smaller. For wavenumbers below the upper limit of the inertial subrange, the spectra and correlations measured with both grids at all speeds can be represented using the von Karman interpolation formula with a single velocity and length scale. The spectra maybe accurately represented over the entire wavenumber range by a modification of the von Karman interpolation formula that includes the effects of dissipation. These models are most accurate at the higher speeds (30 and 40 m/s).

Bereketab, Semere↗

NASA NTRS1991

Multigrid solution of compressible turbulent flow on unstructured meshes using a two-equation model

The system of equations consisting of the full Navier-Stokes equations and two turbulence equations has been solved for in the steady-state using a multigrid strategy on unstructured meshes. The flow equations and turbulence equations are solved in a loosely coupled manner. The flow equations are advanced in time using a multistage Runge-Kutta time stepping scheme with a stability bound local time-step, while the turbulence equations are advanced in a point-implicit scheme with a time-step which guarantees stability and positivity. Low Reynolds number modifications to the original two-equation model are incorporated in a manner which results in well behaved equations for arbitrarily small wall distances. A variety of aerodynamic flows are solved for, initializing all quantities with uniform freestream values, and resulting in rapid and uniform convergence rates for the flow and turbulence equations.

Mavriplis, D. J.↗

NASA NTRS1996

Simulation and Modeling of the Elliptic Streamline Flow

Direct numerical simulation (DNS) was used to study homogeneous turbulence with a mean flow having elliptic streamlines. This flow combines the effects of rotation and strain on the turbulence. There are many important reasons for studying the elliptic streamline flow. This flow contains the effects of both rotation and strain and is therefore similar to the mean flow in a vortex strained in the plane perpendicular to its axis. Such flows provide insight into the fundamental vortical interactions within turbulence. A strained vortex also occurs in airplane wakes, in which each wingtip vortex induces a strain field on the other. The strain field can affect the stability of these vortices and thereby their turbulent structure downstream. The ability to understand and predict the turbulent structure of the vortices is important to the wake hazard problem which is of major concern for the safety of commercial aircraft.

Blaisdell, Gregory A.↗

NASA NTRS1981

On the OLP prediction of the unstable modes of the flat plate turbulent boundary layer

The Orr et al. (OLP) stability analysis is reformulated to yield profiles from hot-wire velocity covariance data in the boundary layer of a flow over a flat plate. The proper orthogonal decomposition theorem (PODT) is used to maximize the mean square of the inner product of the velocity vector with a deterministic candidate vector function in Hilbert space. The OLP formulation couples the PODT for the eddy viscosity with a global extremum principle for the disturbance kinetic energy. Derivation of the governing equations is described analytically, as are the extraction of the mean velocity profile and the eddy viscosity. Finally, techniques for characterizing the dominant mode of the velocity field are introduced and applied to three-dimensional flow over a flat plate. Good agreement is found between the predictions, using laboratory data, and the measured energy and three-dimensional structures of the dominant eddies.

Wadia, A. R.↗

NASA NTRS2001

A 7.5-Year Dataset of SSM/I-Derived Surface Turbulent Fluxes Over Global Oceans

The global air-sea turbulent fluxes are needed for driving ocean models and validating coupled ocean-atmosphere global models. A method was developed to retrieve surface air humidity from the radiances measured by the Special Sensor Microwave/Imager (SSM/I) Using both SSM/I-retrieved surface wind and air humidity, they computed daily turbulent fluxes over global oceans with a stability-dependent bulk scheme. Based on this method, we have produced Version 1 of Goddard Satellite-Based Surface Turbulent Fluxes (GSSTF) dataset from the SSM/I data and other data. It provides daily- and monthly-mean surface turbulent fluxes and some relevant parameters over global oceans for individual F8, F10, and F11 satellites covering the period July 1987-December 1994. It also provides 1988-94 annual- and monthly-mean climatologies of the same variables, using only F8 and F1 1 satellite data. It has a spatial resolution of 2.0 degrees x 2.5 degrees lat-long and is archived at the NASA/GSFC DAAC. The purpose of this paper is to present an updated assessment of the GSSTF 1.0 dataset.

Chou, Shu-Hsien↗

NASA NTRS1976

Similarity theory of the buoyantly interactive planetary boundary layer with entrainment

A similarity model is developed for the vertical profiles of turbulent flow variables in an entraining turbulent boundary layer of arbitrary buoyant stability. In the general formulation the vertical profiles, internal rotation of the velocity vector, discontinuities or jumps at a capping inversion and bulk aerodynamic coefficients of the boundary layer are given by solutions to a system of ordinary differential equations in the similarity variable. To close the system, a formulation for buoyantly interactive eddy diffusivity in the boundary layer is introduced which recovers Monin-Obukhov similarity near the surface and incorporates a hypothesis accounting for the observed variation of mixing length throughout the boundary layer. The model is tested in simplified versions which depend only on roughness, surface buoyancy, and Coriolis effects by comparison with planetary-boundary-layer wind- and temperature-profile observations, measurements of flat-plate boundary layers in a thermally stratified wind tunnel and observations of profiles of terms in the turbulent kinetic-energy budget of convective planetary boundary layers. On balance, the simplified model reproduced the trend of these various observations and experiments reasonably well, suggesting that the full similarity formulation be pursued further.

Hoffert, M. I.↗

NASA NTRS1994

Roles of drizzle in a one-dimensional third-order turbulence closure model of the nocturnal stratus-topped marine boundary layer

This study focuses on the effects of drizzle in a one-dimensional third-order turbulence closure model of the nocturnal stratus-topped marine boundary layer. When the simulated drizzle rate is relatively small (maximum approximately equal to 0.6 mm/day), steady-state solutions are obtained. The boundary layer stabilizes essentially because drizzle causes evaporative cooling of the subcloud layer. This stabilization considerably reduces the buoyancy flux and turbulence kinetic energy below the stratus cloud. Thus, drizzle tends to decouple the cloud from the subcloud layer in the model, as suggested by many observational studies. In addition, the evaporation of drizzle in the subcloud layer creates small scattered clouds, which are likely to represent cumulus clouds, below the solid stratus cloud in the model. The sensitivity experiments show that these scattered clouds help maintain a coupled boundary layer. When the drizzle rate is relatively large (maximum approximately equal to 0.9 mm/day), the response of the model becomes transient with bursts in turbulent fluxes. This phenomenon is related to the formation of the scattered cloud layer below the solid stratus cloud. It appears that the model is inadequate to represent the heat and moisture transport by strong updrafts covering a small fractional area in cumulus convection.

Wang, Shouping↗

NASA NTRS1972

Instabilities and turbulence in highly ionized plasmas in a magnetic field

Physical mechanisms were considered which are responsible for plasma turbulence and the establishment of necessary conditions for energy exchange and transfer in the frequency spectrum. In addition, work was performed to better understand the drift instability in the highly inhomogeneous Rensselaer arc, and methods to suppress this instability using feedback stabilization techniques. Correlation techniques were refined to study plasma turbulence, the diffusion wave technique for monitoring cross-field diffusion was extended to include regimes of high turbulence levels, and a technique for coupling stabilizing RF power to the Rensselaer arc was developed.

Jennings, W. C.↗

NASA NTRS1999

Entropy Splitting and Numerical Dissipation

A rigorous stability estimate for arbitrary order of accuracy of spatial central difference schemes for initial-boundary value problems of nonlinear symmetrizable systems of hyperbolic conservation laws was established recently by Olsson and Oliger (1994) and Olsson (1995) and was applied to the two-dimensional compressible Euler equations for a perfect gas by Gerritsen and Olsson (1996) and Gerritsen (1996). The basic building block in developing the stability estimate is a generalized energy approach based on a special splitting of the flux derivative via a convex entropy function and certain homogeneous properties. Due to some of the unique properties of the compressible Euler equations for a perfect gas, the splitting resulted in the sum of a conservative portion and a non-conservative portion of the flux derivative. hereafter referred to as the "Entropy Splitting." There are several potential desirable attributes and side benefits of the entropy splitting for the compressible Euler equations that were not fully explored in Gerritsen and Olsson. The paper has several objectives. The first is to investigate the choice of the arbitrary parameter that determines the amount of splitting and its dependence on the type of physics of current interest to computational fluid dynamics. The second is to investigate in what manner the splitting affects the nonlinear stability of the central schemes for long time integrations of unsteady flows such as in nonlinear aeroacoustics and turbulence dynamics. If numerical dissipation indeed is needed to stabilize the central scheme, can the splitting help minimize the numerical dissipation compared to its un-split cousin? Extensive numerical study on the vortex preservation capability of the splitting in conjunction with central schemes for long time integrations will be presented. The third is to study the effect of the non-conservative proportion of splitting in obtaining the correct shock location for high speed complex shock-turbulence interactions. The fourth is to determine if this method can be extended to other physical equations of state and other evolutionary equation sets. If numerical dissipation is needed, the Yee, Sandham, and Djomehri (1999) numerical dissipation is employed. The Yee et al. schemes fit in the Olsson and Oliger framework.

Yee, H. C.↗

NASA NTRS2021

Boundary Layer Instabilities Over a Cone-Cylinder-Flare Model at Mach 6

Computations are performed to investigate the boundary-layer instabilities over a sharp cone-cylinder-flare model at zero degrees angle of attack. The model geometry and the flow conditions are selected to match the experiments conducted in the Boeing/AFOSR Mach 6 Quiet Tunnel (BAM6QT) at Purdue University. The geometry consists of a nominally sharp 5-degree half-angle cone, followed by a cylindrical segment and then a 10-degree flare. An axisymmetric separation bubble is generated as a result of the laminar shock/boundary-layer interaction in the cylinder-flare region. The comparison of the laminar flow solution and the schlieren images shows a remarkable agreement between the respective locations of both the boundarylayer edge and the reattachment shock. The predicted heat flux distribution is also in agreement with the measured values downstream of the reattachment location. The analysis of convective and global instabilities is performed for flare half angles equal to 8, 10, and 12 degrees and nosetip radii equal to 0.1, 1, and 5 mm. The linear amplification of first and second Mack mode instabilities that begin to amplify in the cone region are computed with a combination of the parabolized stability equations (PSE) and the harmonic linearized Navier-Stokes equations (HLNSE). The predicted frequency spectra of the surface pressure fluctuations associated with both planar and oblique instability waves are compared with the measured spectra at the various locations of the PCB and Kulite sensors. The comparison shows that the computational analysis captures the distinct lobes within the disturbance amplification spectra measured in the experiments, but some differences in amplification characteristics are noted at low frequencies. Overall, the oblique disturbances are found to be more amplified than the planar disturbances. To our knowledge, this represents the first successful comparison between convective instability analysis and measured surface pressure fluctuations for a hypersonic configuration with a separation bubble. Finally, the global instability analysis shows that the laminar flow becomes supercritical for flare half angles larger than 8 degrees. The unstable global mode for the experimental configuration of a 10 degrees flare and a sharp nosetip cone corresponds to a stationary three-dimensional disturbance that is concentrated in the recirculation region and achieves its maximum growth rate for an azimuthal wavenumber of 5.

Boundary layer transition↗

NASA NTRS2016

Nonlinear Transient Growth and Boundary Layer Transition

Parabolized stability equations (PSE) are used in a variational approach to study the optimal, non-modal disturbance growth in a Mach 3 at plate boundary layer and a Mach 6 circular cone boundary layer. As noted in previous works, the optimal initial disturbances correspond to steady counter-rotating streamwise vortices, which subsequently lead to the formation of streamwise-elongated structures, i.e., streaks, via a lift-up effect. The nonlinear evolution of the linearly optimal stationary perturbations is computed using the nonlinear plane-marching PSE for stationary perturbations. A fully implicit marching technique is used to facilitate the computation of nonlinear streaks with large amplitudes. To assess the effect of the finite-amplitude streaks on transition, the linear form of plane- marching PSE is used to investigate the instability of the boundary layer flow modified by spanwise periodic streaks. The onset of bypass transition is estimated by using an N- factor criterion based on the amplification of the streak instabilities. Results show that, for both flow configurations of interest, streaks of sufficiently large amplitude can lead to significantly earlier onset of transition than that in an unperturbed boundary layer without any streaks.

Laminar turbulent transition↗