Search NASA⌕ Search

Engineering topics

Moin, P.

Publications and source records attributed to Moin, P..

At least 19 records
All sourcesNASA NTRSDOE OSTI
NASA NTRS2003

Large-Eddy Simulation of Conductive Flows at Low Magnetic Reynolds Number

In this paper we study the LES method with dynamic procedure in the context of conductive flows subject to an applied external magnetic field at low magnetic Reynolds number R(sub m). These kind of flows are encountered in many industrial applications. For example, in the steel industry, applied magnetic fields can be used to damp turbulence in the casting process. In nuclear fusion devices (Tokamaks), liquid-lithium flows are used as coolant blankets and interact with the surrounding magnetic field that drives and confines the fusion plasma. Also, in experimental facilities investigating the dynamo effect, the flow consists of liquid-sodium for which the Prandtl number and, as a consequence, the magnetic Reynolds number is low. Our attention is focused here on the case of homogeneous (initially isotropic) decaying turbulence. The numerical simulations performed mimic the thought experiment described in Moffatt in which an initially homogeneous isotropic conductive flow is suddenly subjected to an applied magnetic field and freely decays without any forcing. Note that this flow was first studied numerically by Schumann. It is well known that in that case, extra damping of turbulence occurs due to the Joule effect and that the flow tends to become progressively independent of the coordinate along the direction of the magnetic field. Our comparison of filtered direct numerical simulation (DNS) predictions and LES predictions show that the dynamic Smagorinsky model enables one to capture successfully the flow with LES, and that it automatically incorporates the effect of the magnetic field on the turbulence. Our paper is organized as follows. In the next section we summarize the LES approach in the case of MHD turbulence at low R(sub m) and recall the definition of the dynamic Smagorinsky model. In Sec. 3 we describe the parameters of the numerical experiments performed and the code used. Section 4 is devoted to the comparison of filtered DNS results and LES results. Conclusions are presented in Sec. 5.

Knaepen, B.↗
NASA NTRS1996

A new approach to turbulence modeling

A new approach to Reynolds averaged turbulence modeling is proposed which has a computational cost comparable to two equation models but a predictive capability approaching that of Reynolds stress transport models. This approach isolates the crucial information contained within the Reynolds stress tensor, and solves transport equations only for a set of 'reduced' variables. In this work, Direct Numerical Simulation (DNS) data is used to analyze the nature of these newly proposed turbulence quantities and the source terms which appear in their respective transport equations. The physical relevance of these quantities is discussed and some initial modeling results for turbulent channel flow are presented.

Perot, B.↗
NASA NTRS1994

Developments and applications of dynamic models for large eddy simulation of complex flows

The dynamic modeling procedure for large eddy simulation of turbulent flows is reviewed and recent developments in the theoretical aspects and applications are described. Methods for inclusion of backscatter of energy from small to large scale motions are presented. New formulations of the dynamic procedure are proposed which are optimized based on the subgrid scale flux vector or the energy dissipation rate instead of the subgrid scale stress tensor. Recent results from application of the model to forced isotropic turbulence with an inertial subrange, flow over a backward facing step at Reynolds number of 28000, and flow over a concave curved surface are presented.

Moin, P.↗
NASA NTRS1991

A dynamic subgrid-scale model for compressible turbulence and scalar transport

The dynamic subgrid-scale (SGS) model of Germano et al. (1991) is generalized for the large eddy simulation (LES) of compressible flows and transport of a scalar. The model was applied to the LES of decaying isotropic turbulence, and the results are in excellent agreement with experimental data and direct numerical simulations. The expression for the SGS turbulent Prandtl number was evaluated using direct numerical simulation (DNS) data in isotropic turbulence, homogeneous shear flow, and turbulent channel flow. The qualitative behavior of the model for turbulent Prandtl number and its dependence on molecular Prandtl number, direction of scalar gradient, and distance from the wall are in accordance with the total turbulent Prandtl number from the DNS data.

Moin, P.↗
NASA NTRS1991

Numerical study of axial turbulent flow over long cylinders

The effects of transverse curvature are investigated by means of direct numerical simulations of turbulent axial flow over cylinders. Two cases of Reynolds number of about 3400 and layer-thickness-to-cylinder-radius ratios of 5 and 11 were simulated. All essential turbulence scales were resolved in both calculations, and a large number of turbulence statistics were computed. The results are compared with the plane channel results of Kim et al. (1987) and with experiments. With transverse curvature the skin friction coefficient increases and the turbulence statistics, when scaled with wall units, are lower than in the plane channel. The momentum equation provides a scaling that collapses the cylinder statistics, and allows the results to be interpreted in light of the plane channel flow. The azimuthal and radial length scales of the structures in the flow are of the order of the cylinder diameter. Boomerang-shaped structures with large spanwise length scales were observed in the flow.

Neves, J. C.↗
NASA NTRS1990

Direct numerical simulation of a three-dimensional turbulent boundary layer

The effects of transverse strain on an initially two-dimensional turbulent boundary layer are studied in a direct numerical simulation of a planar channel flow with impulsively started transverse pressure gradient. Consistent with experiments in three-dimensional boundary layers, the simulation shows a decrease in the Reynolds shear stress with increasing transverse strain. Also, the directions of the Reynolds shear stress vector and the mean velocity gradient vector were found to differ. In addition, the simulation shows a drop in the turbulent kinetic energy. Terms in the Reynolds stress transport equations were computed. The balances indicate that the decrease in turbulent kinetic energy is a result of a decrease in turbulence production, along with an increase in turbulent dissipation. Intuitive reasoning and current turbulence models would predict an increase in kinetic energy along with increases in production and dissipation rates as a result of increased mean-flow strain rate. Later in the evolution of the flow, both turbulence production and dissipation increase.

Moin, P.↗
NASA NTRS1990

Active turbulence control in wallbounded flow using direct numerical simulations

An exploratory study of concepts for active control of turbulent boundary layers using the direct numerical simulation technique was performed. Significant drag reduction was achieved when the surface boundary condition was modified such that it could suppress the large-scale structures present in the wall region. This was achieved by prescribing the normal component of velocity at the wall to be 180 deg out of phase with the normal velocity slightly above the wall at each instant. The drag reduction was accompanied with significant reduction in the intensity of the wall-layer structures and reductions in the magnitude of Reynolds stresses throughout the flow. Suitability of wall-pressure and shear-stress fluctuations for detection of flow structures above the wall was examined. A preliminary result obtained by applying the present control strategy to a transitional flow is also briefly described, from which one can infer a possible linkage between the control strategy and flow stability.

Kim, J.↗
NASA NTRS1989

On the active control of wall-bounded turbulent flows

The direct numerical simulation technique is used to explore concepts for manipulation of turbulent boundary layers. Significant drag reduction was achieved when at each instant the normal component of velocity at the wall was prescribed to be 180 deg out of phase with the normal velocity slightly above the wall. The drag reduction is accompanied with significant reduction in the intensity of the wall-layer structures and reductions in the magnitude of Reynolds stresses throughout the flow. Suitability of wall pressure and shear stress fluctuations for detection of flow structures above the wall are examined.

Moin, P.↗
NASA NTRS1989

Numerical simulation of a three-dimensional turbulent boundary layer

The effects of transverse strain on an initially two-dimensional turbulent boundary layer are studied in a direct numerical simulation of a planar channel flow with impulsively started transverse pressure gradient. Consistent with experiments in three-dimensional boundary layers, the simulation shows a drop in the Reynolds shear stress with increasing transverse strain. Also, the directions of the Reynolds shear stress vector and the mean velocity gradient vector were found to differ. In addition, the simulation shows a drop in the turbulent kinetic energy. The terms in the budget of the Reynolds stress were computed. In general, the individual terms do not show a significant change. However, an increase in dissipation rate rather than reduction in total production appears to be responsible for the drop in turbulent kinetic energy.

Moin, P.↗
NASA NTRS1989

Cross correlation and length scales in turbulent flows near surfaces

Two kinds of length scales are used in turbulent flows; 'functional length scales' such as mixing length, dissipation length L(sub epsilon), etc., and 'flow-field length scales' derived from cross correlations of velocity, pressure, etc. in the flow. Some connection between these scales are derived here. We first consider the cross correlation R(sub vv)(y,y(sub 1)) of the normal components u at two heights y, y(sub 1) above a rigid surface, normalized by the velocity y(sub 1) (greater than y). For shear-free boundary layers it is found theoretically, and in field and numerical experiments that R(sub vv) approximately equals y/y(sub 1). For shear layers it is also found that R(sub vv) approximately equals f(y/y(sub 1)) less than or equal to y,y(sub 1). This function f differs slightly between low Reynolds number numerical simulations and field experiments. The lateral structure defined by R(sub vv)(y,r(sub 3); y(sub 1),0) is also self similar and shows that the eddies centered at about y(sub 1) appear to have constant lateral width a(sub 3) above and below y(sub 1), where a(sub 3, sup +) approximately equals 7+1/(1.4dU(sup +)/dy(sup +)), when normalized on u(sub *) and v, where U is the mean velocity. Results for L(sub epsilon, sup -1) from direct numerical simulation are found to compare well with the formula L(sub epsilon, sup -1) = A(sub B)/y + A(sub S)dU/dy/v, for unidirectional and reversing turbulent boundary layers and channel flow, except near where dU/dy approximately equals 0. The conclusion is that the large-scale eddy structure and length scales in these flows are determined by a combination of shear and blocking, and that the vertical component of turbulence has a self-similar structure in both kinds of boundary layer.

Hunt, J. C. R.↗
NASA NTRS1988

Dynamics of coherent structures in a plane mixing layer

An incompressible, time developing 3-D mixing layer with idealized initial conditions was simulated numerically. Consistent with the suggestions from experimental measurements, the braid region between the dominant spanwise vortices or rolls develops longitudinal vortices or ribs, which are aligned upstream and downstream of a roll and produce spanwise distortion of the rolls. The process by which this distortion occurs is explained by studying a variety of quantities of dynamic importance (e.g., production of enstrophy, vortex stretching). Other quantities of interest (dissipation, helicity density) are also computed and discussed. The currently available simulation only allows the study of the early evolution (before pairing) of the mixing layer. New simulations in progress will relieve this restriction.

Hussain, Fazle↗
NASA NTRS1988

Sampling inhomogeneous turbulent fields

The reconstruction of an inhomogeneous random process from a finite number of discrete samples can be performed in terms of the Karhunen-Loeve (KL) expansion for that process. The n(th) eigenfunction has n - 1 zero crossings which are the sampling points for the inhomogeneous process. The rapid variation of the KL eigenfunctions makes it unnecessary to have a high density of sampling (or grid points) near the wall. However, this result should not be construed to indicate that with spectral simulations significantly fewer grid points are required with the KL expansion as compared to other orthogonal expansions. Moin and Moser (1989) have shown that the advantage of the KL expansion over Chebychev expansion rapidly diminishes when high percentage (say 90 percent) energy recovery is demanded.

Adrian, R. J.↗
NASA NTRS1988

Dynamical interpretation of conditional patterns

While great progress is being made in characterizing the 3-D structure of organized turbulent motions using conditional averaging analysis, there is a lack of theoretical guidance regarding the interpretation and utilization of such information. Questions concerning the significance of the structures, their contributions to various transport properties, and their dynamics cannot be answered without recourse to appropriate dynamical governing equations. One approach which addresses some of these questions uses the conditional fields as initial conditions and calculates their evolution from the Navier-Stokes equations, yielding valuable information about stability, growth, and longevity of the mean structure. To interpret statistical aspects of the structures, a different type of theory which deals with the structures in the context of their contributions to the statistics of the flow is needed. As a first step toward this end, an effort was made to integrate the structural information from the study of organized structures with a suitable statistical theory. This is done by stochastically estimating the two-point conditional averages that appear in the equation for the one-point probability density function, and relating the structures to the conditional stresses. Salient features of the estimates are identified, and the structure of the one-point estimates in channel flow is defined.

Adrian, R. J.↗
NASA NTRS1988

Eddies, streams, and convergence zones in turbulent flows

Recent studies of turbulent shear flows have shown that many of their important kinematical and dynamical properties can be more clearly understood by describing the flows in terms of individual events or streamline patterns. These events or flow regions are studied because they are associated with relatively large contributions to certain average properties of the flow, for example kinetic energy, Reynolds stress, or to particular processes in the flow, such as mixing and chemical reactions, which may be concentrated at locations where streamlines converge for fast chemical reactions (referred to as convergence or C regions), or in recirculating eddying regions for slow chemical reactions. The aim of this project was to use the numerical simulations to develop suitable criteria for defining these eddying or vortical zones. The C and streaming (S) zones were defined in order to define the whole flow field. It is concluded that homogeneous and sheared turbulent flow fields are made up of characteristic flow zones: eddy, C, and S zones. A set of objective criteria were found which describe regions in which the streamlines circulate, converge or diverge, and form high streams of high velocity flow.

Hunt, J. C. R.↗
NASA NTRS1988

Karhunen-Loeve expansion of Burgers' model of turbulence

The properties of the Karhunen-Loeve expansion of a strongly inhomogeneous random process are examined with emphasis on applications to turbulent flow fields. The ability of the KL expansion to represent functions that have both slow and rapid variations in a relatively small number of expansion terms is tested on a one-dimensional model based on the forced Burgers' equation. The rate of the convergence of the expansion is evaluated, and its dependence on the Reynolds number is determined. It is shown that the KL eigenfunctions possess wall boundary layers attached to outer structures that are independent of the Reynolds number (at high Reynolds numbers). It is also shown that the spectrum of eigenvalues is broad at large Reynolds numbers, requiring many terms to represent higher-order derivatives of the function.

Chambers, D. H.↗
NASA NTRS1988

Reynolds-stress and dissipation-rate budgets in a turbulent channel flow

The budgets for the Reynolds stresses and for the dissipation rate of the turbulence kinetic energy are computed using direct simulation data of a turbulent channel flow. The budget data reveal that all the terms in the budget become important close to the wall. For inhomogeneous pressure boundary conditions, the pressure-strain term is split into a return term, a rapid term, and a Stokes term. The Stokes term is important close to the wall. The rapid and return terms play different roles depending on the component of the term. A split of the velocity pressure-gradient term into a redistributive term and a diffusion term is proposed, which should be simpler to model. The budget data is used to test existing closure models for the pressure-strain term, the dissipation rate, and the transport rate. In general, further work is needed to improve the models.

Mansour, N. N.↗
NASA NTRS1988

Ejection mechanisms in the sublayer of a turbulent channel

The structure of the vorticity field in the viscous wall layer of a turbulent channel is studied by examining the results of a fully resolved direct numerical simulation. It is shown that this region is dominated by intense three-dimensional shear layers in which the dominant vorticity component is spanwise. The advection and reproduction processes of these structures are examined and shown to be consistent with the classical generation mechanism for two-dimensional Tollmien-Schlichting waves. This process is fundamentally different from the usually accepted mechanism involving hairpin vortices.

Jimenez, Javier↗