Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “turbulence 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 397 records · Page 22

A comparative study of turbulence models for overset grids

The implementation of two different types of turbulence models for a flow solver using the Chimera overset grid method is examined. Various turbulence model characteristics, such as length scale determination and transition modeling, are found to have a significant impact on the computed pressure distribution for a multielement airfoil case. No inherent problem is found with using either algebraic or one-equation turbulence models with an overset grid scheme, but simulation of turbulence for multiple-body or complex geometry flows is very difficult regardless of the gridding method. For complex geometry flowfields, modification of the Baldwin-Lomax turbulence model is necessary to select the appropriate length scale in wall-bounded regions. The overset grid approach presents no obstacle to use of a one- or two-equation turbulence model. Both Baldwin-Lomax and Baldwin-Barth models have problems providing accurate eddy viscosity levels for complex multiple-body flowfields such as those involving the Space Shuttle.

Renze, Kevin J.↗

Local isotropy in high Reynolds number turbulent shear flows

This is a report on the continuation of experiments, which Dr. Srinivas Veeravalli and the present author started in 1991, to investigate the hypothesis of local isotropy in shear flows. This hypothesis, which states that at sufficiently high Reynolds numbers the small-scale structures of turbulent motions are independent of large-scale structures and mean deformations, has been used in theoretical studies of turbulence and computational methods like large-eddy simulation. The importance of Kolmogorov's ideas arises from the fact that they create a foundation for turbulence theory.

Saddoughi, Seyed G.↗

Status of Turbulence Modeling for Hypersonic Propulsion Flowpaths

This report provides an assessment of current turbulent flow calculation methods for hypersonic propulsion flowpaths, particularly the scramjet engine. Emphasis is placed on Reynolds-averaged Navier-Stokes (RANS) methods, but some discussion of newer meth- ods such as Large Eddy Simulation (LES) is also provided. The report is organized by considering technical issues throughout the scramjet-powered vehicle flowpath including laminar-to-turbulent boundary layer transition, shock wave / turbulent boundary layer interactions, scalar transport modeling (specifically the significance of turbulent Prandtl and Schmidt numbers) and compressible mixing. Unit problems are primarily used to conduct the assessment. In the combustor, results from calculations of a direct connect supersonic combustion experiment are also used to address the effects of turbulence model selection and in particular settings for the turbulent Prandtl and Schmidt numbers. It is concluded that RANS turbulence modeling shortfalls are still a major limitation to the accuracy of hypersonic propulsion simulations, whether considering individual components or an overall system. Newer methods such as LES-based techniques may be promising, but are not yet at a maturity to be used routinely by the hypersonic propulsion community. The need for fundamental experiments to provide data for turbulence model development and validation is discussed.

Georgiadis, Nicholas J.↗

Beam Propagation Through Atmospheric Turbulence Using an Altitude-Dependent Structure Profile with Non-Uniformly Distributed Phase Screens

For free-space optical communication systems, numerical wave optics simulations provide a useful technique for modeling turbulence-induced beam degradation when the analytical theory is insufficient for characterizing the atmospheric channel. Motivated by such applications we use a split-step method modeling the turbulence as a series of random phase screens using the Hufnagel-Valley turbulence profile. We employ a space-to-ground case study to examine the irradiance and phase statistics for uniformly and non-uniformly located screens and find better agreement with theory using a non-uniform discretization minimizing the contribution of each screen to the total scintillation. In this poster, we summarize the method and the results of the case study including a comparison to layered models used in astronomical imaging applications.

Chahine, Yousef K.↗

Flow tracing fidelity of scattering aerosol in laser Doppler velocimetry

An experimental method for the determination of the flow-tracing fidelity of a scattering aerosol used in laser Doppler velocimeters was developed with particular reference to the subsonic turbulence measurements. The method employs the measurement of the dynamic response of a flow-seeding aerosol excited by acoustic waves. The amplitude and frequency of excitation were controlled in order to simulate the corresponding values of fluid turbulence components. Experimental results are presented on the dynamic response of aerosols over the size range from 0.1 to 2.0 microns in diam and over the frequency range 100 Hz to 100 kHz. It was observed that unit-density spherical scatterers with diameters of 0.2 micron followed subsonic air turbulence frequency components up to 100 kHz with 98% fidelity.

Mazumder, M. K.↗

On motion of fluid in boundary layer near line of intersection of two planes

In the paper "The Mutual Interference of Boundary Layers," the authors investigated the problem of the interference of two planes intersecting at right angles on the boundary layers formed by the motion of fluid along the line of intersection of these planes. In the present paper, the results of the preceding one are generalized to the case of planes intersecting at any angle. The motion of a fluid in an angle less than 180 degrees is discussed and the enlargement of the boundary layers near the line of intersection of the planes, the limits of the interference effects of the boundary layers, and the corrections on the drag are determined. All computations are conducted by the Karman-Pohlhausen method for laminar and turbulent boundary layers. The results are reduced to tabulated form.

BOUNDARY LAYER - PLATES↗

Piloted Evaluation of a UH-60 Mixer Equivalent Turbulence Simulation Model

A simulation study of a recently developed hover/low speed Mixer Equivalent Turbulence Simulation (METS) model for the UH-60 Black Hawk helicopter was conducted in the NASA Ames Research Center Vertical Motion Simulator (VMS). The experiment was a continuation of previous work to develop a simple, but validated, turbulence model for hovering rotorcraft. To validate the METS model, two experienced test pilots replicated precision hover tasks that had been conducted in an instrumented UH-60 helicopter in turbulence. Objective simulation data were collected for comparison with flight test data, and subjective data were collected that included handling qualities ratings and pilot comments for increasing levels of turbulence. Analyses of the simulation results show good analytic agreement between the METS model and flight test data, with favorable pilot perception of the simulated turbulence. Precision hover tasks were also repeated using the more complex rotating-frame SORBET (Simulation Of Rotor Blade Element Turbulence) model to generate turbulence. Comparisons of the empirically derived METS model with the theoretical SORBET model show good agreement providing validation of the more complex blade element method of simulating turbulence.

Lusardi, Jeff A.↗

Analytical method for predicting the pressure distribution about a nacelle at transonic speeds

The formulation and development of a computer analysis for the calculation of streamlines and pressure distributions around two-dimensional (planar and axisymmetric) isolated nacelles at transonic speeds are described. The computerized flow field analysis is designed to predict the transonic flow around long and short high-bypass-ratio fan duct nacelles with inlet flows and with exhaust flows having appropriate aerothermodynamic properties. The flow field boundaries are located as far upstream and downstream as necessary to obtain minimum disturbances at the boundary. The far-field lateral flow field boundary is analytically defined to exactly represent free-flight conditions or solid wind tunnel wall effects. The inviscid solution technique is based on a Streamtube Curvature Analysis. The computer program utilizes an automatic grid refinement procedure and solves the flow field equations with a matrix relaxation technique. The boundary layer displacement effects and the onset of turbulent separation are included, based on the compressible turbulent boundary layer solution method of Stratford and Beavers and on the turbulent separation prediction method of Stratford.

Keith, J. S.↗

Numerical techniques for the solution of the compressible Navier-Stokes equations and implementation of turbulence models

The time-splitting explicit numerical method of MacCormack is applied to separated turbulent boundary layer flow problems. Modifications of this basic method are developed to counter difficulties associated with complicated geometry and severe numerical resolution requirements of turbulence model equations. The accuracy of solutions is investigated by comparison with exact solutions for several simple cases. Procedures are developed for modifying the basic method to improve the accuracy. Numerical solutions of high-Reynolds-number separated flows over an airfoil and shock-separated flows over a flat plate are obtained. A simple mixing length model of turbulence is used for the transonic flow past an airfoil. A nonorthogonal mesh of arbitrary configuration facilitates the description of the flow field. For the simpler geometry associated with the flat plate, a rectangular mesh is used, and solutions are obtained based on a two-equation differential model of turbulence.

Baldwin, B. S.↗

Calculation of three-dimensional turbulent subsonic flows in transition ducts

A method for computing three-dimensional turbulent subsonic flow in curved ducts is being developed. A set of tube-like surface oriented coordinates is employed for a general class of geometries applicable to subsonic diffusers with offset bends. The geometric formulation is complex and no previous treatment of this class of viscous flow problems is known to the authors. The duct centerline is a space curve specified by piecewise polynomials. A Frenet frame is located on the centerline at each axial location. The cross sections are described by superellipses imbedded in the Frenet frame. Duct surfaces are also coordinate surfaces, which greatly simplifies the boundary conditions. The resulting coordinates are nonorthogonal. An approximate set of governing equations is employed for viscous flows having strong flow in a primary flow direction. The derivation is coordinate invariant and the resulting equations are expressed in tensor form. These equations are solved by an efficient alternating direction implicit (ADI) method. This numerical method is generally stable and permits solution in difficult geometries using the general tensor formulation.

Levy, R.↗

A scheme for computing surface fluxes from mean flow observations

A computational method for estimating turbulent surface stress, sensible heat flux, and humidity flux from mean velocity, temperature, and humidity at a single height in the atmospheric surface layer is presented. The method is based on coupling a Monin-Obukhov similarity profile to a 'force-restore' formulation for the evolution of surface soil temperature to yield the local values of shear stress, heat flux, and surface temperature. In addition to reference-level mean flow properties, the parameters required to implement this method are the thermal heat capacity of the soil, the surface aerodynamic roughness, latitude, solar declination, surface albedo, surface emissivity, and atmospheric transmissivity. Sample calculations are presented for (1) constant atmospheric forcing at the reference level, and (2) variable atmospheric forcing corresponding to Kahle's (1977) measurements of windspeed, air temperature, and radiometer soil surface temperature under dry vegetatively sparse conditions in the Mohave Desert, Calif.

Hoffert, M. I.↗

Turbulent drag reduction research

Research in the area of turbulent drag reduction for attached flows is summarized. The most promising passive techniques utilize non-planar geometry. Of particular interest is the suitability of these devices for retrofit of existing vehicles. Five methods for reducing turbulent skin friction drag on bodies/fuselages are discussed. They are: (1) large-eddy breakup devices; (2) riblets; (3) slot injection optimization; (4) control of Emmons spot generation; and (5) relaminarization through massive suction. Except for the Emmons spot work these methods all indicate the possibility of sizable net reductions in skin friction for laboratory conditions.

Bushnell, D. M.↗

On Parametric Sensitivity of Reynolds-Averaged Navier-Stokes SST Turbulence Model: 2D Hypersonic Shock-Wave Boundary Layer Interactions

Examined is sensitivity of separation extent, wall pressure and heating to variation of primary input flow parameters, such as Mach and Reynolds numbers and shock strength, for 2D and Axisymmetric Hypersonic Shock Wave Turbulent Boundary Layer interactions obtained by Navier-Stokes methods using the SST turbulence model. Baseline parametric sensitivity response is provided in part by comparison with vetted experiments, and in part through updated correlations based on free interaction theory concepts. A recent database compilation of hypersonic 2D shock-wave/turbulent boundary layer experiments extensively used in a prior related uncertainty analysis provides the foundation for this updated correlation approach, as well as for more conventional validation. The primary CFD method for this work is DPLR, one of NASA's real-gas aerothermodynamic production RANS codes. Comparisons are also made with CFL3D, one of NASA's mature perfect-gas RANS codes. Deficiencies in predicted separation response of RANS/SST solutions to parametric variations of test conditions are summarized, along with recommendations as to future turbulence approach.

Hypersonic↗

A symmetrical laser Doppler velocity meter and its application to turbulence characterization

A symmetrical method of optical heterodyning of the Doppler shifted scattered laser radiation developed for velocity measurements with a minimal instrumental spectral broadening and a high signal-to-noise ratio. The method employs two laser beams incident on the moving scatterer and does not use any reference beam for heterodyning. The Doppler signal frequency is independent of the scattering angle and the signal possesses no receiving aperture broadening. Optical alignment is simple. Typical values of the instrumental spectral broadening were approximately 0.8 percent of the center frequency of the Doppler signal, and the signal-to-noise ratio was approximately 25 dB, obtained from an air flow system using submicron dioctylphthalate scattering aerosol. Experimental and theoretical studies were made on the characteristics of the Doppler signal and the effect of system parameters in turbulent flow measurement. The optimization process involved in the beam optics and in the use of a spatial filter is described. For localized flow measurement in any direction of the three-dimensional orthogonal coordinates, the system, using uncorrected optical components, had a sensing volume which can be described by a sensitive length of 600 microns and a diameter of 100 microns.

Mazumder, M. K.↗

Transonic viscous-inviscid interaction using Euler and inverse boundary-layer equations

The paper is concerned with the use of a zonal method for the computation of transonic viscous-inviscid interacting flow about airfoils. The inviscid portion of the flow is treated by using an Euler equation solution method, while an inverse integral compressible turbulent boundary-layer solution method is used for the viscous portion of the flow. The matching of the viscous and inviscid solutions is discussed, and some numerical results as well as comparisons with experimental data are presented.

Whitfield, D. L.↗

A three-dimensional incompressible flow simulation method and its application to the Space Shuttle main engine. II Turbulent flow

An implicit finite difference code cast in general curvilinear coordinates is further developed for three-dimensional incompressible turbulent flows. The code is based on the method of pseudocompressibility and utilizes the Beam and Warming implicit approximate factorization algorithm to achieve computational efficiency. A multiple-zone method is further extended to include composite-grids to overcome the excessive computer memory required for solving turbulent flows in complex three-dimensional geometries. A simple turbulence model is proposed for internal flows. The code is being used for the Space Shuttle Main Engine (SSME) internal flow analyses.

Chang, J. L. C.↗

Large-Eddy Simulation of Variable Speed Power Turbine Cascade with Inflow Turbulence

Numerical results are presented from theNASAGlenn Research Center’sin-house turbomachinery code, Glenn-HT applied to the Variable Speed Power Turbine (VSPT) experiment at the NASA Transonic Turbine Blade Cascade Facility.The main goal of this paper is to implement a digital filtering method to generate turbulence upstream and a sub-grid model (Localized dynamic k-equation model (LDKM)) in the framework of LES in order to investigate the effect of inflow turbulence on the transition seen in the VSPT experimental data at the cruise condition (incidence angle of 40oand Tu = 0.5%, 5%,10%, and 15%). Although the boundary layer on thesuction side and pressure side of the blades is initially laminar due to favorable pressure gradient, the laminar flow can transition to turbulent flow past a separation zoneon the suction sideor by natural or by-pass transition. This process determinesthe total-pressure lossesin the wake. Therefore, it is desirable to develop a reliable prediction tool to accurately capture the transition mechanism in blade rows operated under the conditions of low Reynolds number andat a variety offree stream turbulence conditions. Our numerical studies reveal that the location of separationis rather insensitive to the level of Tu, however the effect ofincreasingTu seems to be in reducing the size and ultimately suppressing the separation bubble. In addition, we performedspectral analysis to identify the peak frequencies in the region where the separation bubble is formed, which provides valuable insights into the transition/separation mechanism.

CFD↗

The nonlinear Galerkin method: A multi-scale method applied to the simulation of homogeneous turbulent flows

Using results of Direct Numerical Simulation (DNS) in the case of two-dimensional homogeneous isotropic flows, the behavior of the small and large scales of Kolmogorov like flows at moderate Reynolds numbers are first analyzed in detail. Several estimates on the time variations of the small eddies and the nonlinear interaction terms were derived; those terms play the role of the Reynolds stress tensor in the case of LES. Since the time step of a numerical scheme is determined as a function of the energy-containing eddies of the flow, the variations of the small scales and of the nonlinear interaction terms over one iteration can become negligible by comparison with the accuracy of the computation. Based on this remark, a multilevel scheme which treats differently the small and the large eddies was proposed. Using mathematical developments, estimates of all the parameters involved in the algorithm, which then becomes a completely self-adaptive procedure were derived. Finally, realistic simulations of (Kolmorov like) flows over several eddy-turnover times were performed. The results are analyzed in detail and a parametric study of the nonlinear Galerkin method is performed.

Debussche, A.↗