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 307 records · Page 17

Development and application of a non-Gaussian atmospheric turbulence model for use in flight simulators

A method is described for generating time histories which model the frequency content and certain non-Gaussian probability characteristics of atmospheric turbulence including the large gusts and patchy nature of turbulence. Methods for time histories using either analog or digital computation are described. A STOL airplane was programmed into a 6-degree-of-freedom flight simulator, and turbulence time histories from several atmospheric turbulence models were introduced. The pilots' reactions are described.

Reeves, P. M.↗

Electromagnetic precipitation and ducting of particles in turbulent boundary layers

A method for analyzing magnetic migration of particles in turbulent flows is applied to the prediction of particle trajectories and densities in turbulent aerodynamic boundary layers. Results for conditions typical of aircraft with 30-40 micron particles indicate a large upstream collection and a 5% loss of particles during one pass through the boundary layer. The capacity of the magnetic field to achieve a balance with turbulent diffusion in confining the particles to the boundary layer is discussed.

Davey, K. R.↗

A general method for calculating three-dimensional compressible laminar and turbulent boundary layers on arbitrary wings

The method described utilizes a nonorthogonal coordinate system for boundary-layer calculations. It includes a geometry program that represents the wing analytically, and a velocity program that computes the external velocity components from a given experimental pressure distribution when the external velocity distribution is not computed theoretically. The boundary layer method is general, however, and can also be used for an external velocity distribution computed theoretically. Several test cases were computed by this method and the results were checked with other numerical calculations and with experiments when available. A typical computation time (CPU) on an IBM 370/165 computer for one surface of a wing which roughly consist of 30 spanwise stations and 25 streamwise stations, with 30 points across the boundary layer is less than 30 seconds for an incompressible flow and a little more for a compressible flow.

Cebeci, T.↗

Atmospheric turbulence chamber for optical transmission experiment Characterization by thermal method

Consideration is given to an atmospheric turbulence chamber designed for optical wave propagation experiments. The chamber consists of ten small electric heater/blowers with an aluminum foil screen and three screens of 2-mm aluminum wire meshes. Calculations are made of the temperature structure constant squared on the basis of temperature structure function measurements derived from a differential microthermocouple system. Values are presented for the refractive-index structure constant squared. The average wind velocity and temperature are found to be, respectively, 0.41 m/sec and 53 C. The inner and outer scales of turbulence are 5.0 mm and 6.5 cm. It is shown that the measured temperature structure function and the power spectrum of temperature fluctuations satisfy, respectively, the 2/3 and -5/3 power similarity laws in the inertial subrange. Possible chamber improvements are discussed.

Gamo, H.↗

Prediction of Very High Reynolds Number Compressible Skin Friction

Flat plate skin friction calculations over a range of Mach numbers from 0.4 to 3.5 at Reynolds numbers from 16 million to 492 million using a Navier Stokes method with advanced turbulence modeling are compared with incompressible skin friction coefficient correlations. The semi-empirical correlation theories of van Driest; Cope; Winkler and Cha; and Sommer and Short T' are used to transform the predicted skin friction coefficients of solutions using two algebraic Reynolds stress turbulence models in the Navier-Stokes method PAB3D. In general, the predicted skin friction coefficients scaled well with each reference temperature theory though, overall the theory by Sommer and Short appeared to best collapse the predicted coefficients. At the lower Reynolds number 3 to 30 million, both the Girimaji and Shih, Zhu and Lumley turbulence models predicted skin-friction coefficients within 2% of the semi-empirical correlation skin friction coefficients. At the higher Reynolds numbers of 100 to 500 million, the turbulence models by Shih, Zhu and Lumley and Girimaji predicted coefficients that were 6% less and 10% greater, respectively, than the semi-empirical coefficients.

Carlson, John R.↗

Some problems of the calculation of three-dimensional boundary layer flows on general configurations

An accurate solution of the three-dimensional boundary layer equations over general configurations such as those encountered in aircraft and space shuttle design requires a very efficient, fast, and accurate numerical method with suitable turbulence models for the Reynolds stresses. The efficiency, speed, and accuracy of a three-dimensional numerical method together with the turbulence models for the Reynolds stresses are examined. The numerical method is the implicit two-point finite difference approach (Box Method) developed by Keller and applied to the boundary layer equations by Keller and Cebeci. In addition, a study of some of the problems that may arise in the solution of these equations for three-dimensional boundary layer flows over general configurations.

Cebeci, T.↗

Modeling Turbulence in Galactic Centers

Turbulence is a prevalent phenomenon in the interstellar medium, and in particular, the environment at the centers of galaxies. For example, detailed observations of the Milky Way’s Central Molecular Zone (CMZ) revealed that it has a complex and turbulent structure. Turbulence on galactic scales is often modeled using star formation and feedback. However, these effects do not appear to be sufficient for explaining the high-velocity dispersion observed in the CMZ, indicating that additional gas-stirring processes are likely to be operating. Here we introduce a proof-of-concept method to drive turbulence in gas that orbits under the influence of a galactic potential. Instead of relying on a particular physical mechanism, we have adopted a Fourier forcing module and have applied it using a smoothed particle hydrodynamics code. To test our method, we performed simulations of a simplistic model of the CMZ. Our turbulence injection method is capable of balancing the self-gravity of the gas, which allows us to run the simulations for long timescales and thereby follow the evolution of the CMZ. Our results show that turbulence induces a flocculent spiral pattern in our model, analogous to that found in galactic-scale simulations. Furthermore, we find that our turbulence injection method induces inward migration of gas, a result consistent with previous numerical simulations. We submit that this injection method is a promising new tool to simulate turbulence in galactic centers.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

On the use of wall functions as boundary conditions for two-dimensional separated compressible flows

A new and improved wall function method for compressible turbulent flows has been developed and tested. This method is applicable to attached and separated flows, to both high- and low-Reynolds number flows, and to flows with adiabatic and nonadiabatic surfaces. This wall function method has been applied to the Launder-Spalding k-epsilon two-equation model of turbulence. The tests consist of comparisons of calculated and experimental results for: (1) an axisymmetrical transonic shock-wave/boundary-wave interaction flow at low Reynolds number in an adiabatic tube, (2) an axisymmetrical high-Reynolds number transonic flow over a nonadiabatic bump, and (3) a two-dimensional supersonic high-Reynolds number flow on a nonadiabatic deflected flap. Each of these experiments had significant regions of flow separation. The calculations are performed with an implicit algorithm that solves the Reynolds-averaged Navier-Stokes equations. It is shown that the results obtained agree very well with the data for the complex compressible flows tested.

Viegas, J. R.↗

Discrimination of coherent features in turbulent boundary layers by the entropy method

Entropy in information theory is defined as the expected or mean value of the measure of the amount of self-information contained in the ith point of a distribution series x sub i, based on its probability of occurrence p(x sub i). If p(x sub i) is the probability of the ith state of the system in probability space, then the entropy, E(X) = - sigma p(x sub i) logp (x sub i), is a measure of the disorder in the system. Based on this concept, a method was devised which sought to minimize the entropy in a time series in order to construct the signature of the most coherent motions. The constrained minimization was performed using a Lagrange multiplier approach which resulted in the solution of a simultaneous set of non-linear coupled equations to obtain the coherent time series. The application of the method to space-time data taken by a rake of sensors in the near-wall region of a turbulent boundary layer was presented. The results yielded coherent velocity motions made up of locally decelerated or accelerated fluid having a streamwise scale of approximately 100 nu/u(tau), which is in qualitative agreement with the results from other less objective discrimination methods.

Corke, T. C.↗

Numerical experiments in homogeneous turbulence

The direct simulation methods developed by Orszag and Patternson (1972) for isotropic turbulence were extended to homogeneous turbulence in an incompressible fluid subjected to uniform deformation or rotation. The results of simulations for irrotational strain (plane and axisymmetric), shear, rotation, and relaxation toward isotropy following axisymmetric strain are compared with linear theory and experimental data. Emphasis is placed on the shear flow because of its importance and because of the availability of accurate and detailed experimental data. The computed results are used to assess the accuracy of two popular models used in the closure of the Reynolds-stress equations. Data from a variety of the computed fields and the details of the numerical methods used in the simulation are also presented.

Rogallo, R. S.↗

An explicit finite-volume time-marching procedure for turbulent flow calculations

A method was developed which calculates two-dimensional, transonic, viscous flow in ducts. The finite-volume, time-marching formulation is used to obtain steady flow solutions of the Reynolds-averaged form of the Navier-Stokes equations. The entire calculation is performed in the physical domain. Control volumes are chosen so that smoothing of flow properties, typically required for stability, is not required. Different time steps are used in the different governing equations. A new pressure interpolation scheme is introduced which improves the shock capturing ability of the method. A multi-volume method for pressure changes in the boundary layer allows calculations which use very long and thin control volumes (length/height - 1000). The method is compared with two test cases. Essentially incompressible turbulent boundary layer flow in an adverse pressure gradient is calculated and the computed distributions of mean velocity and shear are in good agreement with the measurements. Transonic viscous flow in a converging diverging nozzle is calculated; the Mach number upstream of the shock is approximately 1.25. The agreement between the calculated and measured shock strength and total pressure losses is good.

Nicholson, Stephen↗

An explicit finite-volume time-marching procedure for turbulent flow calculations

A method was developed which calculates two-dimensional, transonic, viscous flow in ducts. The finite-volume, time-marching formulation is used to obtain steady flow solutions of the Reynolds-averaged form of the Navier-Stokes equations. The entire calculation is performed in the physical domain. Control volumes are chosen so that smoothing of flow properties, typically required for stability, is not required. Different time steps are used in the different governing equations. A new pressure interpolation scheme is introduced which improves the shock capturing ability of the method. A multi-volume method for pressure changes in the boundary layer allows calculations which use very long and thin control volumes (length/height - 1000). The method is compared with two test cases. Essentially incompressible turbulent boundary layer flow in an adverse pressure gradient is calculated and the computed distributions of mean velocity and shear are in good agreement with the measurements. Transonic viscous flow in a converging diverging nozzle is calculated; the Mach number upstream of the shock is approximately 1.25. The agreement between the calculated and measured shock strength and total pressure losses is good.

Nicholson, Stephen↗

A method for obtaining a statistically stationary turbulent free shear flow

The long-term goal of the current research is the study of Large-Eddy Simulation (LES) as a tool for aeroacoustics. New algorithms and developments in computer hardware are making possible a new generation of tools for aeroacoustic predictions, which rely on the physics of the flow rather than empirical knowledge. LES, in conjunction with an acoustic analogy, holds the promise of predicting the statistics of noise radiated to the far-field of a turbulent flow. LES's predictive ability will be tested through extensive comparison of acoustic predictions based on a Direct Numerical Simulation (DNS) and LES of the same flow, as well as a priori testing of DNS results. The method presented here is aimed at allowing simulation of a turbulent flow field that is both simple and amenable to acoustic predictions. A free shear flow is homogeneous in both the streamwise and spanwise directions and which is statistically stationary will be simulated using equations based on the Navier-Stokes equations with a small number of added terms. Studying a free shear flow eliminates the need to consider flow-surface interactions as an acoustic source. The homogeneous directions and the flow's statistically stationary nature greatly simplify the application of an acoustic analogy.

Timson, Stephen F.↗

Computation of turbulent flows over backward-facing step

A numerical method for computing incompressible turbulent flows is presented. The method is tested by calculating laminar recirculating flows and is applied in conjunction with a modified Kappa-epsilon model to compute the flow over a backward-facing step. In the laminar regime, the computational results are in good agreement with the experimental data. The turbulent flow study shows that the reattachment length is underpredicted by the standard Kappa-epsilon model. The addition of a term to the standard model that accounts for the effects of rotation on turbulent flow improves the results in the recirculation region and increases the computed reattachment length.

Mansour, N. N.↗

Methods for estimating pressure and thermal loads induced by elevon deflections on hypersonic-vehicle surfaces with turbulent boundary layers

Empirical anaytic methods are presented for calculating thermal and pressure distributions in three-dimensional, shock-wave turbulent-boundary-layer, interaction-flow regions on the surface of controllable hypersonic aircraft and missiles. The methods, based on several experimental investigations, are useful and reliable for estimating both the extent and magnitude of the increased thermal and pressure loads on the vehicle surfaces.

Kaufman, L. G., II↗

Estimating unsteady aerodynamic forces on a cascade in a three-dimensional turbulence field

An analytical method has been developed to estimate the unsteady aerodynamic forces caused by flow field turbulence on a wind tunnel turning vane cascade system (vane set). This method approximates dynamic lift and drag by linearly perturbing the appropriate steady state force equations, assuming that the dynamic loads are due only to free stream turbulence and that this turbulence is homogeneous, isotropic, and Gaussian. Correlation and unsteady aerodynamic effects are also incorporated into the analytical model. Using these assumptions, equations relating dynamic lift and drag to flow turbulence, mean velocity, and vane set geometry are derived. From these equations, estimates for the power spectra and rms (root mean squared value, delta) loading of both lift and drag can be determined.

Norman, T.↗