SEARCH · Search NASA
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.
A survey designed to provide an introduction to the subject of turbulence modeling, and to explain the need for such models is given. The subject is developed along chronological lines since this provides a logical development plan and also because it then moves from relatively simple phenomenological models through more complicated procedures and ultimately to the subject of large-eddy simulation.
An investigation into the existence of hairpin vortices in turbulent channel flow is conducted using a database generated by the large eddy simulation technique. It is shown that away from the wall the distribution of the inclination angle of vorticity vector attains its maximum at about 45 deg to the wall. Two point correlations of velocity and vorticity fluctuations strongly support a flow model consisting of vortical structures inclined at 45 deg to the wall. The instantaneous vorticity vectors plotted in planes inclined at 45 deg show that the flow contains an appreciable number of hairpins. Vortex lines are used to display the three dimensional structure of hairpins, which are shown to be generated from deformation of transverse vortex filaments.
Several conditional sampling techniques are applied to a data base generated by large-eddy simulation of turbulent channel flow. It is shown that the bursting process is associated with well-organized horseshoe vortices inclined at about 45 deg. to the wall. These vortical structures are identified by examining the vortex lines of three-dimensional, ensemble averaged vorticity fields. Two distinct horseshoe-shaped vortices corresponding to the sweep and ejection events are detected. These vortices are associated with high Reynolds shear stress and hence make a significant contribution to turbulent energy production. The dependency of the ensemble averaged vortical structures on the detection criteria, and the question of whether this ensemble-averaged structure is an artifact of the ensemble averaging process are examined. The ensemble-averaged pattern of these vortical structures that emerge from the analysis could provide the basis for a hypothetical model of the organized structures of wall-bounded shear flows.
Models employed in the computation of turbulent flows are described and their application to internal flows is evaluated by examining the predictions of various turbulence models in selected flow configurations. The main conclusions are: (1) the k-epsilon model is used in a majority of all the two-dimensional flow calculations reported in the literature; (2) modified forms of the k-epsilon model improve the performance for flows with streamline curvature and heat transfer; (3) for flows with swirl, the k-epsilon model performs rather poorly; the algebraic stress model performs better in this case; and (4) for flows with regions of secondary flow (noncircular duct flows), the algebraic stress model performs fairly well for fully developed flow, for developing flow, the algebraic stress model performance is not good; a Reynolds stress model should be used. False diffusion and inlet boundary conditions are discussed. Countergradient transport and its implications in turbulence modeling is mentioned. Two examples of recirculating flow predictions obtained using PHOENICS code are discussed. The vortex method, large eddy simulation (modeling of subgrid scale Reynolds stresses), and direct simulation, are considered. Some recommendations for improving the model performance are made. The need for detailed experimental data in flows with strong curvature is emphasized.
The effects of prescribed wall motion on turbulent channel flow were examined with the objective of understanding the drag mechanisms and obtaining a possible means for drag reduction. A computer program for large-eddy simulation (LES) of turbulent flow in a channel was adapted to treat the case of time-varying wall motion using linearized boundary conditions. The code was applied to flow in a channel with prescribed wall motion on one wall. It was found that the results can be explained by a model based on a Stokes-like layer near the moving wall. This provides a theory which can be used to guide the selection of parameters. Considering both viscous and pressure drag components, it was found that net drag reduction is possible under certain conditions of wall motion. However, the work required to drive the wall motion can be greater than the reduced through-flow energy requirement. The major effects of the wall motion are restricted to the viscous sublayer of the turbulent flow.
Turbulence structures in a wall bounded shear layer during the bursting event detected by a conditional sampling technique are investigated using data obtained from large eddy simulation of turbulent channel flow. Streamlines are constructed from the ensemble-averaged velocity field to illustrate the flow patterns associated with the bursting event. The splatting motions during the sweep event and the existence of a pair of counterrotating streamwise vortices during the ejection process are exhibited.
A transport equation for the root-mean-square pressure fluctuations of turbulent flow is derived from the time-dependent momentum equation for incompressible flow. Approximate modeling of this transport equation is included to relate terms with higher order correlations to the mean quantities of turbulent flow. Three empirical constants are introduced in the model. Two of the empirical constants are estimated from homogeneous turbulence data and wall pressure fluctuations measurements. The third constant is determined by comparing the results of large eddy simulations for a plane channel flow and an annulus flow.
An investigation into the existence of hairpin vortices in turbulent channel flow is conducted using a database generated by the large eddy simulation technique. It is shown that away from the wall the distribution of the inclination angle of vorticity vector attains its maximum at about 45 deg to the wall. Two point correlations of velocity and vorticity fluctuations strongly support a flow model consisting of vortical structures inclined at 45 deg to the wall. The instantaneous vorticity vectors plotted in planes inclined at 45 deg show that the flow contains an appreciable number of hairpins. Vortex lines are used to display the three dimensional structure of hairpins, which are shown to be generated from deformation of transverse vortex filaments.
Several conditional sampling techniques are applied to a data base generated by large-eddy simulation of turbulent channel flow. It is shown that the bursting process is associated with well-organized horseshoe vortices inclined at about 45 deg to the wall. These vortical structures are identified by examining the vortex lines of three-dimensional, ensemble averaged vorticity fields. Two distinct horseshoe-shaped vortices corresponding to the sweep and ejection events are detected. These vortices are associated with high Reynolds shear stress and hence make a significant contribution to turbulent energy production. The dependency of the ensemble averaged vortical structures on the detection criteria, and the question of whether this ensemble-averaged structure is an artifact of the ensemble averaging process are examined. The ensemble-averaged pattern of these vortical structures that emerge from the analysis could provide the basis for a hypothetical model of the organized structures of wall-bounded shear flows.
The introduction of the supercomputer and recent advances in both Reynolds averaged, and large eddy simulation fluid flow approximation techniques to the Navier-Stokes equations, have created a robust environment for the exploration of problems of interest to the Navy in general, and the Naval Underwater Systems Center in particular. The nature of problems that are of interest, and the type of resources needed for their solution are addressed. The goal is to achieve a good engineering solution to the fluid-structure interaction problem. It is appropriate to indicate that a paper by D. Champman played a major role in developing the interest in the approach discussed.
New three-dimensional spectral multigrid algorithms are analyzed and implemented to solve the variable coefficient Helmholtz equation. Periodicity is assumed in all three directions which leads to a Fourier collocation representation. Convergence rates are theoretically predicted and confirmed through numerical tests. Residual averaging results in a spectral radius of 0.2 for the variable coefficient Poisson equation. In general, non-stationary Richardson must be used for the Helmholtz equation. The algorithms developed are applied to the large-eddy simulation of incompressible isotropic turbulence.
The requirements for computers for aerospace vehicle design problems, in particular fluid dynamics, are examined. The levels of approximation to Navier-Stokes equations are: (1) linearized inviscid form, (2) nonlinear inviscid equation, (3) Reynolds-averaged Navier-Stokes equation, (4) large eddy simulation with models for subgrid scale turbulence, and (5) full Navier-Stokes equations to calculate laminar/turbulent transition and turbulence dissipation. The Cray 2 with 256 megaword memory and 4.1 nsec clock and the Cray X-MP/48 with a 128 megaword memory and a 9.5 nsec clock of the Ames facility are described, and an example revealing their applicability is presented. The use of zonal methods to simulate aerodynamic flow fields numerically is discussed. Consideration is given to the development of a program for simulation of aerodynamic flow fields, the compilers, symbolic debuggers, the scheduling of multiple CPUs, and the analysis and storage of the data.
A survey designed to provide an introduction to the subject of turbulence modeling, and to explain the need for such models is given. The subject is developed along chronological lines since this provides a logical development plan and also because it then moves from relatively simple phenomenological models through more complicated procedures and ultimately to the subject of large-eddy simulation.
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.
From preliminary calculations it was concluded that the two-point Eddy Damping Quasi-Normal Markovian (EDQNM) closure accurately describes the behavior of second order moments. This closure can be applied as subgrid and supergrid models for Large Eddy Simulations at higher Reynolds numbers. In the case of homogeneous anisotropic turbulence, when the nonlinear terms are introduced the calculation becomes quite onerous but is still considerably less expensive than the calculation of a DS. The major merit of two-point closure models is that they can be easily applied to flows at Reynolds numbers that are unreachable by a DS. Work is in progress to derive expressions for the nonlinear terms that give good global conservation properties.
An alternative to the one-point closure model for turbulence, the large eddy simulation (LES), together with its more exact relative, direct numerical simulation (DNS) are discussed. These methods are beginning to serve as partial substitutes for turbulence experiments. The eddy damped quasi-normal Markovian (EDQNM) theory is reviewed. Angular distribution of the converted data was examined in relationship to EDQNM.
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.
New three-dimensional spectral multigrid algorithms are analyzed and implemented to solve the variable coefficient Helmholtz equation. Periodicity is assumed in all three directions which leads to a Fourier collocation representation. Convergence rates are theoretically predicted and confirmed through numerical tests. Residual averaging results in a spectral radius of 0.2 for the variable coefficient Poisson equation. In general, non-stationary Richardson must be used for the Helmholtz equation. The algorithms developed are applied to the large-eddy simulation of incompressible isotropic turbulence.