Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “turbulent flows”

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 1,117 records · Page 62
All sourcesNASA NTRSDOE OSTI
NASA NTRS1988

Holocinematographic Velocimeter

Holographic movies used to study turbulent flows. High-speed movie composed of single-exposure holograms of tracer particles in each of two simultaneous orthogonal views of flow. Time resolution set by movie frame rate, and total time limited by length of movie. Data from system used to verify results of theoretical models and numerical simulations of turbulence, visualize and measure coherent structures in "non-simple" turbulent flows, to examine mechanisms operative in various turbulence-control/drag-reduction concepts.

Weinstein, Leonard M.↗
NASA NTRS2013

Influence of Computational Drop Representation in LES of a Droplet-Laden Mixing Layer

Multiphase turbulent flows are encountered in many practical applications including turbine engines or natural phenomena involving particle dispersion. Numerical computations of multiphase turbulent flows are important because they provide a cheaper alternative to performing experiments during an engine design process or because they can provide predictions of pollutant dispersion, etc. Two-phase flows contain millions and sometimes billions of particles. For flows with volumetrically dilute particle loading, the most accurate method of numerically simulating the flow is based on direct numerical simulation (DNS) of the governing equations in which all scales of the flow including the small scales that are responsible for the overwhelming amount of dissipation are resolved. DNS, however, requires high computational cost and cannot be used in engineering design applications where iterations among several design conditions are necessary. Because of high computational cost, numerical simulations of such flows cannot track all these drops. The objective of this work is to quantify the influence of the number of computational drops and grid spacing on the accuracy of predicted flow statistics, and to possibly identify the minimum number, or, if not possible, the optimal number of computational drops that provide minimal error in flow prediction. For this purpose, several Large Eddy Simulation (LES) of a mixing layer with evaporating drops have been performed by using coarse, medium, and fine grid spacings and computational drops, rather than physical drops. To define computational drops, an integer NR is introduced that represents the ratio of the number of existing physical drops to the desired number of computational drops; for example, if NR=8, this means that a computational drop represents 8 physical drops in the flow field. The desired number of computational drops is determined by the available computational resources; the larger NR is, the less computationally intensive is the simulation. A set of first order and second order flow statistics, and of drop statistics are extracted from LES predictions and are compared to results obtained by filtering a DNS database. First order statistics such as Favre averaged stream-wise velocity, Favre averaged vapor mass fraction, and the drop stream-wise velocity, are predicted accurately independent of the number of computational drops and grid spacing. Second order flow statistics depend both on the number of computational drops and on grid spacing. The scalar variance and turbulent vapor flux are predicted accurately by the fine mesh LES only when NR is less than 32, and by the coarse mesh LES reasonably accurately for all NR values. This is attributed to the fact that when the grid spacing is coarsened, the number of drops in a computational cell must not be significantly lower than that in the DNS.

Bellan, Josette↗
NASA NTRS1988

A closure model of diffusion transport of the Reynolds-stress equations and its application to a turbulent step flow

A study is made for the development and computations of the separating and reattaching shear flows. The highlight of the study is an attempt to predict the third-moments of turbulent velocity which is responsible for the diffusion transport of the Reynolds stresses. The present computations show that the third-moments obtained by employing the low-Reynolds number model of transport equations improve the prediction of the third-moments.

Amano, R. S.↗
NASA NTRS1989

Numerical simulation of separated and vortical flows on bodies at large angles of attack

Navier-Stokes computations were carried out for high-angle-of-attack flow about bodies of revolution in supersonic and subsonic flow. Numerical and physical factors affecting flow predictions are discussed. Accurate prediction of the complex three-dimensional separated flow is found to depend on using grids sufficiently fine to resolve the details of both the viscous boundary layer and the off-surface separated flow structures, and, for turbulent flow cases, use of an eddy-viscosity turbulence model which accounts for the leeward vortical flow structures. An example of high-incidence flow computed about a more complex geometry, the F-18 fighter fuselage forebody and leading edge extension, is also presented.

Schiff, Lewis B.↗
NASA NTRS1993

Nongray radiative interactions in entry region flows with variable properties

The entrance region flow between parallel plates is investigated numerically and the variation of fluid properties with temperature is considered. The effect of variable properties for both laminar and turbulent flows is presented. Laminar flows are found to be affected to a larger extent than turbulent flows, as far as the effect of variable properties are concerned.

Krishnamurthy, R.↗
NASA NTRS2018

Tetrahedral-Mesh Simulations of Shock-Turbulence Interaction

Despite decades of development of unstructured mesh methods, direct numerical simulations (DNS) of turbulent flows are still predominantly performed on structured or unstructured hexahedral meshes with high-order finite-difference methods, weighted essentially nonoscillatory (WENO) schemes, or hybrid schemes formed by their combinations. Tetrahedral meshes offer easy mesh generation and adaptation around complex geometries and the potential of an orientation-free grid that would benefit the isotropic nature of small-scale dissipation, as well as the solution accuracy of intermediate scales. To advance the state of the art of unstructured-mesh simulation capabilities for shock/turbulence interaction, DNS using pure tetrahedral meshes are carried out with the space-time conservation element, solution element (CESE) method in this research. By its design, the CESE method is constructed based on a non-dissipative scheme and is a genuinely multidimensional numerical framework that is free from the use of an approximate Riemann-solver. The numerical framework also provides the ability to add numerical dissipation (the nondissipative scheme acts as the reference state like that of the reversible state in thermodynamics) when needed (with justification from mathematics/physics). The above-mentioned features along with the CESE method's consistent shock-capturing approach and strong enforcement of flux conservation in spacetime offers a novel method to accurately simulate turbulent flows and their interaction with shocks using tetrahedral meshes. Two canonical problems, namely, isotropic turbulence interaction with a normal shock and a Mach 2.9 turbulent boundary layer flow over a 24deg compression corner are investigated in this study. Computational results show reasonably good agreement with experimental data and results from structured-mesh, high-order simulations available in the literature. Successful validation of these canonical problems demonstrated here paves the way for future high-fidelity supersonic flow simulations involving complex-geometries.

Venkatachari, Balaji Shankar↗
NASA NTRS1993

The evolution equation for the flame surface density in turbulent premixed combustion

The mean reaction rate in flamelet models for turbulent premixed combustion depends on two basic quantities: a mean chemical rate, called the flamelet speed, and the flame surface density. Our previous work had been primarily focused on the problem of the structure and topology of turbulent premixed flames, and it was then determined that the flamelet speed, when space-averaged, is only weakly sensitive to the turbulent flow field. Consequently, the flame surface density is the key quantity that conveys most of the effects of the turbulence on the rate of energy release. In flamelet models, this quantity is obtained via a modeled transport equation called the Sigma-equation. Past theoretical work has produced a rigorous approach that leads to an exact but unclosed formulation for the turbulent Sigma-equation. In the exact Sigma-equation, it appears that the dynamical properties of the flame surface density are determined by a single parameter, namely the turbulent flame stretch. Unfortunately, the turbulent flame stretch as well as the flame surface density is not available from experiments, and, in the absence of experimental data, little is known on the validity of the closure assumptions used in current flamelet models. Direct Numerical Simulation (DNS) is the alternative approach to get basic information on these fundamental quantities. In the present work, three-dimensional DNS of premixed flames in isotropic turbulent flow is used to estimate the different terms appearing in the Sigma-equation. A new methodology is proposed to provide the source and sink terms for the flame surface density, resolved both temporally and spatially throughout the turbulent flame brush. Using this methodology, our objective is to extract the turbulent flame stretch from the DNS data base and then perform extensive comparisons with flamelet models. Thanks to the detailed information produced by the DNS-based analysis, it is expected that this type of comparison will not only underscore the shortcomings of current models, but also suggest ways to improve them.

Trouve, Arnaud↗
NASA NTRS2022

Prediction of Turbulent Diffusing Flows Using FUN3D

The FUN3D code was used to perform Reynolds-averaged Navier-Stokes (RANS) simulations to compute subsonic flow in an S-duct diffuser and transonic flow in a two-dimensional diffuser using the Speziale-Sarkar-Gatski/Launder-Rodi-Reece (SSG/LRR) Reynolds stress model (RSM). For comparison purposes, additional simulations were run with the one-equation Spalart-Allmaras (SA) and the two-equation Menter Shear-Stress Transport (SST) turbulence models. Each model was run with and without the quadratic constitutive relation (QCR) for computing the turbulent stresses. It was shown that the simulations that utilized the RSM had better overall predictions of the diffusive flow fields compared to the simulations that utilized the one and two equation turbulence models.

Inlet Distortion↗
NASA NTRS1984

Small-amplitude viscous motion on arbitrary potential flows

This paper is concerned with small-amplitude, unsteady, vortical and entropic motion imposed on steady potential flows. It is restricted to the case where the spatial scale of the unsteady motion is small compared to that of the mean flow. Under such conditions, the unsteady motion may be influenced by viscosity even if the mean flow is not. An exact high-frequency (small-wavelength) solution is obtained for the small-amplitude viscous motion imposed on a steady potential flow. It generalizes the one obtained by Pearson (1959) for the homogeneous-strain case to the case of quasi-homogeneous strain. This result is used to study the effect of viscosity on rapidly distorted turbulent flows. Specific numerical results are given for a turbulent flow near a two-dimensional stagnation point.

Goldstein, M. E.↗
NASA NTRS1993

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.↗
NASA NTRSDate not supplied

Prediction of Turbulent Diffusing Flows Using FUN3D

The FUN3D code was used to perform Reynolds-averaged Navier-Stokes (RANS) simulations to compute subsonic flow in an S-duct diffuser and transonic flow in a two-dimensional diffuser using the Speziale-Sarkar-Gatski/Launder-Rodi-Reece (SSG/LRR) Reynolds stress model (RSM). For comparison purposes, additional simulations were run with the one-equation Spalart-Allmaras (SA) and the two-equation Menter Shear-Stress Transport (SST) turbulence models. Each model was run with and without the quadratic constitutive relation (QCR) for computing the turbulent stresses. It was shown that the simulations that utilized the RSM had better overall predictions of the diffusive flow fields compared to the simulations that utilized the one and two equation turbulence models.

Inlet Distortion↗
NASA NTRS1987

On the subgrid-scale modeling of compressible turbulence

A subgrid-scale model recently derived for use in the large-eddy simulation of compressible turbulent flows is examined from a fundamental theoretical and computational standpoint. It is demonstrated that this model, which is applicable only to compressible turbulent flows in the limit of small density fluctuations, correlates somewhat poorly with the results of direct numerical simulations of compressible isotropic turbulence at low Mach numbers. An alternative model, based on Favre-filtered fields, is suggested which appears to reduce these limitations.

Speziale, C. G.↗
NASA NTRS1988

The subgrid-scale modeling of compressible turbulence

A subgrid-scale model recently derived for use in the large-eddy simulation of compressible turbulent flows is examined from a fundamental theoretical and computational standpoint. It is demonstrated that this model, which is applicable only to compressible turbulent flows in the limit of small density fluctuations, correlates somewhat poorly with the results of direct numerical simulations of compressible isotropic turbulence at low Mach numbers. An alternative model, based on Favre-filtered fields, is suggested which appears to reduce these limitations.

Speziale, C. G.↗
NASA NTRS1995

Analysis of Aerothermal Heating Over Reusable Single-Stage to Orbit Vehicles

Using the General Aerodynamic Simulation Program (GASP) with finite-rate surface catalysis and radiative equilibrium surface temperature conditions, a computational study has been conducted to predict three-dimensional aero-convective heating over the proposed reusable single-stage to orbit vehicles. The solutions for three basic classes of usable single-stage-to-orbit vehicles, including wing body, lifting body and vertical lander, are presented and compared to provide information for vehicle design. The free stream conditions used in this study occur at an altitude of 57.5 km, ambient density of 4.129e-4 kg/cubic m, and velocity of 4,849 m/s. The convective heating over the vehicle surface in both laminar and turbulent flow regimes were computed. The Baldwin-Lomax algebraic model with a modified van Driest damping factor was applied in the turbulent flow calculations. A bifurcation model which guarantees the conservation of mass diffusive fluxes was adopted to account for the multi-species mass diffusion. The effects of turbulence and surface catalysis on the surface heating are discussed in detail The solutions indicate the degree of oxygen dissociation over the leeward side is sensitive to the shape of vehicle body. Predicted convective heating over the entire surface of lifting body is strongly affected by the presence of surface catalysis. For wing body and vertical lander, the surface catalysis effect does not appear to be important in the far downstream region, and turbulence is the dominant effect on the heating prediction in this region.

Chen, Y. K.↗
NASA NTRS1993

Investigation of the asymptotic state of rotating turbulence using large-eddy simulation

Study of turbulent flows in rotating reference frames has long been an area of considerable scientific and engineering interest. Because of its importance, the subject of turbulence in rotating reference frames has motivated over the years a large number of theoretical, experimental, and computational studies. The bulk of these previous works has served to demonstrate that the effect of system rotation on turbulence is subtle and remains exceedingly difficult to predict. A rotating flow of particular interest in many studies, including the present work, is examination of the effect of solid-body rotation on an initially isotropic turbulent flow. One of the principal reasons for the interest in this flow is that it represents the most basic turbulent flow whose structure is altered by system rotation but without the complicating effects introduced by mean strains or flow inhomogeneities. The assumption of statistical homogeneity considerably simplifies analysis and computation. The principal objective of the present study has been to examine the asymptotic state of solid-body rotation applied to an initially isotropic, high Reynolds number turbulent flow. Of particular interest has been to determine the degree of two-dimensionalization and the existence of asymptotic self-similar states in homogeneous rotating turbulence.

Squires, Kyle D.↗
NASA NTRS1980

Three dimensional mean flow and turbulence characteristics of the near wake of a compressor rotor blade

The investigation was carried out using the rotating hot wire technique. Measurements were taken inside the end wall boundary layer to discern the effect of annulus and hub wall boundary layer, secondary flow, and tip leakage on the wake structure. Static pressure gradients across the wake were measured using a static stagnation pressure probe insensitive to flow direction changes. The axial and the tangential velocity defects, the radial component of velocity, and turbulence intensities were found to be very large as compared to the near and far wake regions. The radial velocities in the trailing edge region exhibited characteristics prevalent in a trailing vortex system. Flow near the blade tips found to be highly complex due to interaction of the end wall boundary layers, secondary flows, and tip leakage flow with the wake. The streamwise curvature was found to be appreciable near the blade trailing edge. Flow properties in the trailing edge region are quite different compared to that in the near and far wake regions with respect to their decay characteristics, similarity, etc. Fourier decomposition of the rotor wake revealed that for a normalized wake only the first three coefficients are dominant.

Ravindranath, A.↗
NASA NTRS2003

Numerical Simulation of High Drag Reduction in a Turbulent Channel Flow with Polymer Additives

The addition of small amounts of long chain polymer molecules to wall-bounded flows can lead to dramatic drag reduction. Although this phenomenon has been known for about fifty years, the action of the polymers and its effect on turbulent structures are still unclear. Detailed experiments have characterized two distinct regimes (Warholic et al. 1999), which are referred to as low drag reduction (LDR) and high drag reduction (HDR). The first regime exhibits similar statistical trends as Newtonian flow: the log-law region of the mean velocity profile remains parallel to that of the Newtonian ow but its lower bound moves away from the wall and the upward shift of the log-region is a function of drag reduction, DR. Although streamwise fluctuations are increased and transverse ones are reduced, the shape of the rms velocity profiles is not qualitatively modified. At higher drag reductions, of the order of 40-50%, the ow enters the HDR regime for which the slope of the log-law is dramatically augmented and the Reynolds shear stress is small (Warholic et al. 1999; Ptasinski et al. 2001). The drag reduction is eventually bounded by a maximum drag reduction (MDR) (Virk & Mickley 1970) which is a function of the Reynolds number. While several experiments report mean velocity profiles very close to the empirical profile of Virk & Mickley (1970) for MDR conditions, the observations regarding the structure of turbulence can differ significantly. For instance, Warholic et al. (1999) measured a near-zero Reynolds shear stress, whereas a recent experiment (Ptasinski et al. 2001) shows evidence of non-negligible Reynolds stress in their MDR flow. To the knowledge of the authors, only the LDR regime has been documented in numerical simulations (Sureshkumar et al. 1997; Dimitropoulos et al. 1998; Min et al. 2001; Dubief & Lele 2001; Sibilla & Baron 2002). This paper discusses the simulation of polymer drag reduced channel ow at HDR using the FENE-P (Finite Elastic non-linear extensibility-Peterlin) model which was used for the first LDR simulation by Sureshkumar et al. (1997). Flow and polymer parameters are close to realistic polymer drag reducing conditions. High drag reductions are achieved by using finite differences and a robust time stepping technique. A minimal channel flow is also used as a numerical experiment to investigate the effect of the outer region turbulent structures on the overall drag at HDR. The drag reducing action of the model is finally studied through the structure of energy transfers from the polymers to the velocity components. This investigation sheds some light on the details of polymer drag reduction.

Dubief, Yves↗