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 721 records · Page 40
All sourcesNASA NTRSDOE OSTI
NASA NTRS1991

Propulsive lift noise

Propulsive life noise is the increase in noise that occurs when airframe surfaces are placed in the propulsive system's exhaust to increase their lift force. Increased local flow velocities and turbulence levels, due to the propulsive system exhaust gases passing along the airframe lifting surfaces, cause an increase in maximum lift coefficient. The airplane's flight speed required for takeoff and landing can then be significantly reduced, allowing operation from shorter runways than those of conventional commercial airports. Unfortunately, interaction of high velocity turbulent exhaust flow with the airframe's solid surfaces generates additional noise radiation. Aeroacoustic processes that cause propulsive lift noise also are present in airframe noise and propulsive system installation noise. Research studies of propulsive lift noise led to development of improved methods of predicting noise radiation from surfaces in turbulent flows. Noise reduction and prediction methods of aircraft noise are discussed.

Fink, Martin R.↗
NASA NTRS1989

Numerical computation of shock wave-turbulent boundary layer interaction in transonic flow over an axisymmetric curved hill

A control-volume based finite difference computation of a turbulent transonic flow over an axisymmetric curved hill is presented. The numerical method is based on the SIMPLE algorithm, and hence the conservation of mass equation is replaced by a pressure correction equation for compressible flows. The turbulence is described by a k-epsilon turbulence model supplemented by a near-wall turbulence model. In the method, the dissipation rate in the region very close to the wall is obtained from an algebraic equation and that for the rest of the flow domain is obtained by solving a partial differential equation for the dissipation rate. The other flow equations are integrated up to the wall. It is shown that the present turbulence model yields the correct location of the compression shock. The other computational results are also in good agreement with experimental data.

Kim, S.-W.↗
NASA NTRS1992

Laser-spectroscopic measurement techniques for hypersonic, turbulent wind tunnel flows

A review is given of the nature, present status, and capabilities of two laser spectroscopic methods for the simultaneous measurement of temperature, density, and their fluctuations owing to turbulence in high speed wind tunnel flows. One method is based on the two frequency excitation of nitric oxide seeded into a nitrogen flow, using tunable dye lasers. The second, more recent method relies on the excitation of oxygen in air flows using a tunable, ArF excimer laser. Signal are obtained from both the laser induced fluorescence and from Raman scattering of the same laser pulse. Measurements are demonstrated in the turbulent boundary layer of a Mach-2 channel flow.

Mckenzie, Robert L.↗
NASA NTRS1977

Supersonic separated turbulent boundary - layer over a wavy wall

A prediction method is developed for calculating distributions of surface heating rates, pressure and skin friction over a wavy wall in a two-dimensional supersonic flow. Of particular interest is the flow of thick turbulent boundary layers. The surface geometry and the flow conditions considered are such that there exists a strong interaction between the viscous and inviscid flow. First, using the interacting turbulent boundary layer equations, the problem is formulated in physical coordinates and then a reformulation of the governing equations in terms of Levy-Lees variables is given. Next, a numerical scheme for solving interacting boundary layer equations is adapted. A number of modifications which led to the improvement of the numerical algorithm are discussed. Finally, results are presented for flow over a train of up to six waves at various flow conditions.

Polak, A.↗
NASA NTRS1977

A critical study of higher-order numerical methods for solving the boundary-layer equations

A fourth-order box method is presented for calculating numerical solutions to parabolic, partial differential equations in two variables or ordinary differential equations. The method is the natural extension of the second-order Keller Box Scheme to fourth order and is demonstrated with application to the incompressible, laminar and turbulent boundary-layer equations. The efficiency of the present method is compared with other two-point and three-point higher-order methods; namely, the Keller Box Scheme with Richardson extrapolation, the method of deferred corrections, and the three-point spline methods. For equivalent accuracy, numerical results show the present method to be more efficient than the other higher-order methods for both laminar and turbulent flows.

Wornom, S. F.↗
NASA NTRS1979

Axial-flow compressor turning angle and loss by inviscid-viscous interaction blade-to-blade computation

A method for computation of the flow field around an arbitrary airfoil cascade on an axially symmetric blade-to-blade surface was developed which takes into account the development and separation of the blade surface boundary layers and mixing in the wake. The method predicts the overall fluid turning and total pressure loss in the context of an inviscid-viscous interaction scheme. The inviscid flow solution is obtained from a compressible flow matrix method. The viscous flow is obtained from a differential boundary layer method which calculates laminar, transitional and turbulent boundary layers. Provisions for the calculation of laminar and turbulent separation regions were added to the viscous scheme. The combined inviscid-viscous interaction scheme described yields results which are quantitatively consistent with experimental data. This suggests that the physical basis for the interactive system is correct and justifies continued exploration and use of the method.

Hansen, E. C.↗
NASA NTRS1982

Approximate method of predicting heating on the windward side of Space Shuttle Orbiter and comparisons with flight data

An approximate method is developed for predicting laminar and turbulent heating rates on the windward side of the Space Shuttle Orbiter for both the wind-tunnel and flight environments. The method is based on a 'local infinite swept cylinder' analysis and includes both equilibrium-air chemistry and variable boundary-layer-edge entropy. The method is validated by comparing with data from wind-tunnel experiments and from the first and second Space Shuttle flights both along the windward-symmetry plane and in a lateral direction off the symmetry plane. Agreement with the flight data is good from approximately 67 km (peak heating) downward.

Hamilton, H. H., II↗
NASA NTRS1973

Unsteady turbulent boundary layers in two-dimensional, incompressible flow.

The paper presents the results of some calculations aimed at studying the characteristics of turbulent boundary layers in unsteady flow. The calculations were performed using the method of Singleton and Nash, which is based on the turbulent kinetic energy equation. Emphasis is placed on the approach to zero wall stress, and the results show the extent to which onset of this condition can be delayed by the effects of time dependence. The delay is not associated merely with alleviation of the pressure gradients, and quasi-steady methods cannot readily be used to calculate this type of flow.

Nash, J. F.↗
NASA NTRS1992

The Baldwin-Lomax model for separated and wake flows using the entropy envelope concept

Implementation of the Baldwin-Lomax algebraic turbulence model is difficult and ambiguous within flows characterized by strong viscous-inviscid interactions and flow separations. A new method of implementation is proposed which uses an entropy envelope concept and is demonstrated to ensure the proper evaluation of modeling parameters. The method is simple, computationally fast, and applicable to both wake and boundary layer flows. The method is general, making it applicable to any turbulence model which requires the automated determination of the proper maxima of a vorticity-based function. The new method is evalulated within two test cases involving strong viscous-inviscid interaction.

Brock, J. S.↗
NASA NTRS1973

Analysis of the injection of a heated turbulent jet into a cross flow

The development of a theoretical model is investigated of the incompressible jet injection process. The discharge of a turbulent jet into a cross flow was mathematically modeled by using an integral method which accounts for natural fluid mechanisms such as turbulence, entrainment, buoyancy, and heat transfer. The analytical results are supported by experimental data and demonstrate the usefulness of the theory for estimating the trajectory and flow properties of the jet for a variety of injection conditions. The capability of predicting jet flow properties, as well as two- and three-dimensional jet paths, was enhanced by obtaining the jet cross-sectional area during the solution of the conservation equations. Realistic estimates of temperature in the jet fluid were acquired by accounting for heat losses in the jet flow due to forced convection and to entrainment of free-stream fluid into the jet.

Campbell, J. F.↗
NASA NTRS1956

Turbulent shear spectra and local isotropy in the low-speed boundary layer

From measurements of turbulent shear spectra together with previously reported longitudinal turbulent energy spectra, the concept of local isotropy in a low-speed boundary layer was examined. Results of these measurements and measurements of the time derivatives of turbulent velocities in the x- and y- directions for various frequency bands showed no evidence of local isotropy in the boundary layer. Several methods (based on isotropy) of evaluating the turbulent dissipation term failed to give consistent answers, further emphasizing a lack of local isotropy.

BOUNDARY LAYER - COMPLETE WINGS↗
NASA NTRS1990

Computation of three-dimensional turbulent boundary layers with heat transfer in a plane of symmetry using embedded wall-layer functions

In the calculation of turbulent boundary layers, a large number of mesh points are required to adequately resolve the intense variation in the velocity and enthalpy in the near-wall region. A substantial reduction in computational effort may be realized by representing the velocity and enthalpy profiles in the wall layer by analytical embedded functions. The effectively inviscid flow in the outer part of the boundary layer may then be resolved by employing a relatively coarser mesh. To obtain complete profiles, the outer numerical solution is matched asymptotically to the inner wall-layer analytical solution. To date, this approach has been restricted to two-dimensional flows; in the present study, a method which may be utilized for turbulent boundary layers with heat transfer in a plane of symmetry is developed as a first step in the application of the embedded-function method to full three-dimensional flows. The present method uses only about half as many mesh points as that required in a conventional procedure, which calculates the flow all the way to the wall, but there is no degradation in accuracy of the computed results.

Degani, A. T.↗
NASA NTRS1985

Effects of wall curvature on turbulence statistics

A three-dimensional, time-dependent, direct numerical simulation of low-Reynolds number turbulent flow in a mildly curved channel was performed, and the results examined to determine the mechanism by which curvature affects wall-bounded turbulent shear flows. A spectral numerical method with about one-million modes was employed, and no explicit subgrid scale model was used. The effects of curvature on this flow were determined by comparing the concave and convex sides of the channel. The observed effects are consistent with experimental observations for mild curvature. The most significant difference in the turbulence statistics between the concave and convex sides is in the Reynolds shear stress. This is accompanied by significant differences in the terms of the Reynolds shear stress balance equations. In addition, it was found that stationary Taylor-Goertler vortices were present and that they had a significant effect on the flow by contributing to the mean Reynolds shear stress, and by enhancing the difference between the wall shear stresses.

Moser, R. D.↗
NASA NTRS2019

WRLES: Wave Resolving Large-Eddy Simulation Code, Theory and Usage

A computational fluid dynamics code has been developed for large-eddy simulations (LES) of turbulent flow. The code uses high-order of accuracy and high-resolution numerical methods to minimize solution error and maximize the resolution of the turbulent structures. Spatial discretization is performed using explicit central differencing. The central differencing schemes in the code include 2nd- to 12th-order standard central difference methods as well as 7-, 9-, 11- and 13-point dispersion relation preserving schemes. Solution filtering and high-order shock capturing are included for stability. Time discretization is performed using multistage Runge-Kutta methods that are up to 4th order accurate. Several options are available to model turbulence including: Baldwin-Lomax and Spalart-Allmaras Reynolds-averaged Navier-Stokes turbulence models, and Smagorinsky, Dynamic Smagorinsky and Vreman sub-grid scale models for LES. This report presents the theory behind the numerical and physical models used in the code and provides a user's manual to the operation of the code.

High-Order↗
NASA NTRS1986

Use of hot-wire anemometry for turbulence measurements in shock induced flows

A research program is currently being conducted with the aim to investigate the operating environment of future gas turbines. The present paper provides a description of the experimental methods which have been employed in performing turbulence intensity measurements in shock-induced flows. In a discussion of the instrumentation, attention is given to the employed low pressure shock tube, the hot-wire probe, the anemometer, the test facility, the experimental setup, the Kistler pressure transducer, and silicon piezoresistive gages. Aspects of instrumentation calibration are considered along with data corrections, experimental data, and data processing.

Hartung, L. C.↗
NASA NTRS1974

Theory for the effects of turbulence in a planetary atmosphere on radio occultation

Rytov's method is used to formulate the correlation functions for log-amplitude and phase fluctuations for both spherical and plane wave propagation in a turbulent medium whose correlation function for refractive index fluctuations is described by the product of a function of the average coordinate and a function of the difference coordinate. The results are applied to the study of radio occultation effects due to the atmosphere of Venus in the case of a flyby space probe. It is assumed that turbulence of isotropic and smoothly varying characteristics occurs in the Venusian atmosphere.

Woo, R.↗
NASA NTRS1986

Two-phase flow

An experimental program to characterize the spray from candidate nozzles for icing-cloud simulation is discussed. One canidate nozzle, which is currently used for icing research, has been characterized for flow and drop size. The median-volume diameter (MVD) from this air-assist nozzle is compared with correlations in the literature. The new experimental spray facility is discussed, and the drop-size instruments are discussed in detail. Since there is no absolute standard for drop-size measurements and there are other limitations, such as drop -size range and velocity range, several instruments are used and results are compared. A two-phase model was developed at Pennsylvania State University. The model uses the k-epsilon model of turbulence in the continous phase. Three methods for treating the discrete phase are used: (1) a locally homogeneous flow (LHF) model, (2) a deterministic separated flow (DSF) model, and (3) a stochastic separated flow (SSF) model. In the LHF model both phases have the same velocity and temperature at each point. The DSF model provides interphase transport but ignores the effects of turbulent fluctuations. In the SSF model the drops interact with turbulent eddies whose properties are determined by the k-epsilon turbulence model. The two-phase flow model has been extended to include the effects of evaporation and combustion.

Tacina, Robert R.↗
NASA NTRSDate not supplied

A Large-Eddy Simulation Of Turbulent Flow Over A Backward Facing Step

Large-eddy simulation (LES) was used to study the turbulent flow over a backward facing step. The simulations examine the details of the turbulent flow’s separated and reattaching regions created by the step geometry. An explicit high-order of accuracy/high resolution finite-difference computational fluid dynamics solver was used for both LES and Reynolds-averaged Navier-Stokes (RANS) simulations. The experimental data from the backward facing step experiment of Driver and Seegmiller was used for comparison. Data presented are derived from a RANS simulation using the Spalart-Allmaras turbulence model and LES simulations using the implicit LES assumption on three different grids of varying resolution. The RANS solution agrees with past work. It does an adequate job of predicting the flowfield and reattachment point, but misses details such as the strength of the recirculating region and magnitudes of the Reynolds stresses. The synthetic eddy method was used to provide a turbulent inflow to the quasi-2D LES domains. For the LES, the coarse grid was of insufficient resolution and under-predicted the skin friction coefficient and over-predicted the normal stresses. The finer grid LES results accurately predict the skin friction coefficient and provide more accurate Reynolds stresses. Triple product velocity correlations are also computed from the LES and the finer grid results did a good job predicting the magnitudes of the curves and the shapes were reasonably represented. It is recommended that future work include investigating the turbulent inflow domain length, width of the quasi-2D domain, and additional grid refinement

turbulence↗