Engineering Papers⌕ Search

Engineering topics

Piomelli, Ugo

Publications and source records attributed to Piomelli, Ugo.

At least 19 records
All sourcesNASA NTRSDOE OSTI
NASA NTRS2004

Interaction of a Boundary Layer with a Turbulent Wake

The objective of this grant was to study the transition mechanisms on a flat-plate boundary layer interacting with the wake of a bluff body. This is a simplified configuration presented and designed to exemplify the phenomena that occur in multi-element airfoils, in which the wake of an upstream element impinges on a downstream one. Some experimental data is available for this configuration at various Reynolds numbers. The first task carried out was the implementation and validation of the immersed-boundary method. This was achieved by performing calculations of the flow over a cylinder at low and moderate Reynolds numbers. The low-Reynolds number results are discussed, which is enclosed as Appendix A. The high-Reynolds number results are presented in a paper in preparation for the Journal of Fluid Mechanics. We performed calculations of the wake-boundary-layer interaction at two Reynolds numbers, Re approximately equal to 385 and 1155. The first case is discussed and a comparison of the two calculations is reported. The simulations indicate that at the lower Reynolds number the boundary layer is buffeted by the unsteady Karman vortex street shed by the cylinder. This is shown: long streaky structures appear in the boundary layer in correspondence of the three-dimensionalities in the rollers. The fluctuations, however, cannot be self-sustained due to the low Reynolds-number, and the flow does not reach a turbulent state within the computational domain. In contrast, in the higher Reynolds-number case, boundary-layer fluctuations persist after the wake has decayed (due, in part, to the higher values of the local Reynolds number Re achieved in this case); some evidence could be observed that a self-sustaining turbulence generation cycle was beginning to be established. A third simulation was subsequently carried out at a higher Reynolds number, Re=3900. This calculation gave results similar to those of the Re=l155 case. Turbulence was established at fairly low Reynolds number, as a consequence of the high level of the free-stream perturbation. An instantaneous flow visualization for that case is shown. A detailed examination of flow statistics in the transitional and turbulent regions, including the evolution of the turbulent kinetic energy (TKE) budget and frequency spectra showed the formation and evolution of turbulent spots characteristic of the bypass transition mechanism. It was also observed that the turbulent eddies achieved an equilibrium, fully developed turbulent states first, as evidenced by the early agreement achieved by the terms in the TKE budget with those observed in turbulent flows. Once a turbulent Reynolds stress profile had been established, the velocity profile began to resemble a turbulent one, first in the inner region and later in the outer region of the wall layer. An extensive comparison of the three cases, including budgets, mean velocity and Reynolds stress profiles and flow visualization, is included. The results obtained are also presented.

Piomelli, Ugo↗
NASA NTRS2003

Numerical Simulations of Wake/Boundary Layer Interactions

Direct and large-eddy simulations of the interaction between the wake of a circular cylinder and a flat-plate boundary layer are conducted. Two Reynolds numbers are examined. The simulations indicate that at the lower Reynolds number the boundary layer is buffeted by the unsteady Karman vortex street shed by the cylinder. The fluctuations, however, cannot be self-sustained due to the low Reynolds-number, and the flow does not reach a turbulent state within the computational domain. In contrast, in the higher Reynolds-number case, boundary-layer fluctuations persist after the wake has decayed (due, in part, to the higher values of the local Reynolds number Re(sub theta) achieved in this case); some evidence could be observed that a self-sustaining turbulence generation cycle was beginning to be established.

Piomelli, Ugo↗
NASA NTRS1997

On the Computation of Sound by Large-Eddy Simulations

The effect of the small scales on the source term in Lighthill's acoustic analogy is investigated, with the objective of determining the accuracy of large-eddy simulations when applied to studies of flow-generated sound. The distribution of the turbulent quadrupole is predicted accurately, if models that take into account the trace of the SGS stresses are used. Its spatial distribution is also correct, indicating that the low-wave-number (or frequency) part of the sound spectrum can be predicted well by LES. Filtering, however, removes the small-scale fluctuations that contribute significantly to the higher derivatives in space and time of Lighthill's stress tensor T(sub ij). The rms fluctuations of the filtered derivatives are substantially lower than those of the unfiltered quantities. The small scales, however, are not strongly correlated, and are not expected to contribute significantly to the far-field sound; separate modeling of the subgrid-scale density fluctuations might, however, be required in some configurations.

Piomelli, Ugo↗
NASA NTRS1997

Large-Eddy Simulation of Transition to Turbulence in Boundary Layers

Large-eddy simulation results for laminar-to-turbulent transition in a spatially developing boundary layer are presented. The disturbances are ingested into a laminar flow through an unsteady suction-and-blowing strip. The filtered, three-dimensional time- dependent Navier-Stokes equations are integrated numerically using spectral, high-order finite-difference, and three-stage low-storage Runge-Kutta methods. The buffer-domain technique is used for the outflow boundary condition. The localized dynamic model used to parameterize the subgrid-scale stresses begins to have a significant impact at the beginning of the nonlinear transition (or intermittency) region. The flow structures commonly found in experiments are also observed in the present simulation; the computed linear instability modes and secondary instability lambda-vortex structures are in agreement with the experiments, and the streak-like-structures and turbulent statistics compare with both the experiments and the theory. The physics captured in the present LES are consistent with the experiments and the full Navier-Stokes simulation (DNS), at a significant fraction of the DNS cost. A comparison of the results obtained with several SGS models shows that the localized model gives accurate results both in a statistical sense and in terms of predicting the dynamics of the energy-carrying eddies, without ad hoc adjustments.

Huai, Xiao-Li↗
NASA NTRS1997

Numerical Simulations Using the Immersed Boundary Technique

The immersed-boundary method can be used to simulate flows around complex geometries within a Cartesian grid. This method has been used quite extensively in low Reynolds-number flows, and is now being applied to turbulent flows more frequently. The technique will be discussed, and three applications of the method will be presented, with increasing complexity. to illustrate the potential and limitations of the method, and some of the directions for future work.

Piomelli, Ugo↗
NASA NTRS1993

High Reynolds number calculations using the dynamic subgrid-scale stress model

The dynamic subgrid-scale eddy viscosity model was used in the large-eddy simulation of the turbulent flow in a plane channel for Reynolds numbers based on friction velocity and channel half-width ranging between 200 and 2000, a range including values significantly higher than in previous simulations. The computed wall stress, mean velocity, and Reynolds stress profiles compare very well with experimental and direct simulation data. Comparison of higher moments is also satisfactory. Although the grid in the near-wall region is fairly coarse, the results are quite accurate: the turbulent kinetic energy peaks at y(+) of about 12, and the near-wall behavior of the resolved stresses is captured accurately. The model coefficient is o(0.001) in the buffer layer and beyond, where the cutoff wave numbers are in the decaying region of the spectra; in the near-wall region, the cutoff wave numbers are nearer the energy-containing range, and the resolved turbulent stresses become a constant fraction of the resolved stresses. This feature is responsible for the correct near-wall behavior of the model coefficient.

Piomelli, Ugo↗
NASA NTRS1993

Dynamic subgrid-scale modeling for high-speed transitional boundary layers

The subgrid scales are modeled dynamically in a large-eddy simulation of transitional boundary-layer flow along a hollow cylinder at a Mach number of 4.5. The behavior of the dynamic-model coefficients, which is determined from local information in the resolved field, is investigated through both an a priori test with direct numerical simulation data for the same case and a complete large-eddy simulation. Both contractions proposed by Germano et al. (1991) and Lilly (1992) are used for the unique determination of the coefficients of the dynamic model, and their results are compared and assessed. The behavior, as well as the energy cascade of the subgridscale field structure, is investigated at various stages of the transition process.

El-Hady, Nabil M.↗
NASA NTRS1993

Large-Eddy Simulation of Boundary Layer Transition on Swept Wings

The large-eddy simulation of the spatial evolution of a stationary crossflow vortex packet in a three-dimensional boundary layer was performed. Although a coarse grid was used (compared to that required by a direct numerical simulation) the essential features of the disturbance evolution, such as the spanwise disturbance spreading and the vortex rollover, were captured accurately. The eddy viscosity became significant only in the late nonlinear stages of the simulation.

Huai, Xiaoli↗
NASA NTRS1990

A dynamic subgrid-scale eddy viscosity model

One major drawback of the eddy viscosity subgrid-scale stress models used in large-eddy simulations is their inability to represent correctly with a single universal constant different turbulent field in rotating or sheared flows, near solid walls, or in transitional regimes. In the present work, a new eddy viscosity model is presented which alleviates many of these drawbacks. The model coefficient is computed dynamically as the calculation progresses rather than input a priori. The model is based on an algebraic identity (Germano 1990) between the subgrid-scale stresses at two different filtered levels and the resolved turbulent stresses. The subgrid-scale stresses obtained using the proposed model vanish in laminar flow and at a solid boundary, and have the correct asymptotic behavior in the near-wall region of a turbulent boundary layer. The results of large-eddy simulations of transitional and turbulent channel flow that use the proposed model are in good agreement with the direct simulation data.

Germano, Massimo↗
NASA NTRS1990

Subgrid-scale backscatter in transitional and turbulent flows

Most subgrid-scale (SGS) models for large-eddy simulations are absolutely dissipative (that is, they remove energy from the large scales at each point in the physical space). The actual SGS stresses, however, may transfer energy to the large scales (backscatter) at a given location. Direct numerical simulations of turbulent channel flow and compressible isotropic turbulence are used to study the backscatter phenomena. In all flows considered roughly 50 percent of the grid points were experiencing backscatter when a Fourier cutoff filter was used. The backscatter fraction was less with a Gaussian filter, and intermediate with a box filter in physical space. Moreover, the backscatter and forward scatter contributions to the SGS dissipation were comparable, and each was often much larger than the total SGS dissipation. The SGS dissipation (normalized by total dissipation) increased with filter width almost independently of filter type and Reynolds number. The amount of backscatter showed an increasing trend with Reynolds numbers. In the near-wall region of the channel, events characterized by strong Reynolds shear stress correlated fairly well with areas of high SGS dissipation (both forward and backward). In compressible isotropic turbulence similar results were obtained, independent of fluctuation Mach number.

Piomelli, Ugo↗
NASA NTRS1990

Large-eddy simulation of transitional channel flow

A large-eddy simulation (LES) of transition in plane channel flow was carried out. The LES results were compared with those of a fine direct numerical simulation (DNS), and with those of a coarse DNS that uses the same mesh as the LES, but does not use a residual stress model. While at the early stages of transition, LES and coarse DNS give the same results: the presence of the residual stress model was found to be necessary to predict accurately mean velocity and Reynolds stress profiles during the late stages of transition (after the second spike stage). The evolution of single Fourier modes is also predicted more accurately by the LES than by the DNS. As small scales are generated, the dissipative character of the residual stress starts to reproduce correctly the energy cascade. As transition progresses, the flow approaches its fully developed turbulent state, the subgrid scales tend towards equilibrium, and the model becomes more accurate.

Piomelli, Ugo↗
NASA NTRS1990

On the large-eddy simulation of transitional wall-bounded flows

The structure of the subgrid scale fields in plane channel flow has been studied at various stages of the transition process to turbulence. The residual stress and subgrid scale dissipation calculated using velocity fields generated by direct numerical simulations of the Navier-Stokes equations are significantly different from their counterparts in turbulent flows. The subgrid scale dissipation changes sign over extended areas of the channel, indicating energy flow from the small scales to the large scales. This reversed energy cascade becomes less pronounced at the later stages of transition. Standard residual stress models of the Smagorinsky type are excessively dissipative. Rescaling the model constant improves the prediction of the total (integrated) subgrid scale dissipation, but not that of the local one. Despite the somewhat excessive dissipation of the rescaled Smagorinsky model, the results of a large eddy simulation of transition on a flat-plate boundary layer compare quite well with those of a direct simulation, and require only a small fraction of the computational effort. The inclusion of non-dissipative models, which could lead to further improvements, is proposed.

Piomelli, Ugo↗
NASA NTRS1990

Application of renormalization group theory to the large-eddy simulation of transitional boundary layers

An eddy viscosity model based on the renormalization group theory of Yakhot and Orszag (1986) is applied to the large-eddy simulation of transition in a flat-plate boundary layer. The simulation predicts with satisfactory accuracy the mean velocity and Reynolds stress profiles, as well as the development of the important scales of motion. The evolution of the structures characteristic of the nonlinear stages of transition is also predicted reasonably well.

Piomelli, Ugo↗
NASA NTRS1990

On the application of large-eddy simulation to transitional wall-bounded flows

Large-eddy simulations of transition on a flat-plate boundary layer have been carried out. Standard residual stress models of the Smagorinsky type are excessively dissipative, but rescaling the model constant improves the prediction of the total (integrated) subgrid scale dissipation, leading to results that compare favorably with those of a direct simulation, at a small fraction of the computational effort.

Piomelli, Ugo↗
NASA NTRS1989

On the large-eddy simulation of transitional wall-bounded flows

The structure of the subgrid scale fields in plane channel flow has been studied at various stages of the transition process to turbulence. The residual stress and subgrid scale dissipation calculated using velocity fields generated by direct numerical simulations of the Navier-Stokes equations are significantly different from their counterparts in turbulent flows. The subgrid scale dissipation changes sign over extended areas of the channel, indicating energy flow from the small scales to the large scales. This reversed energy cascade becomes less pronounced at the later stages of transition. Standard residual stress models of the Smagorinsky type are excessively dissipative. Rescaling the model constant improves the prediction of the total (integrated) subgrid scale dissipation, but not that of the local one. Despite the somewhat excessive dissipation of the rescaled Smagorinsky model, the results of a large eddy simulation of transition on a flat-plate boundary layer compare quite well with those of a direct simulation, and require only a small fraction of the computational effort. The inclusion of non-dissipative models, which could lead to further improvements, is proposed.

Piomelli, Ugo↗
NASA NTRS1989

New approximate boundary conditions for large eddy simulations of wall-bounded flows

Two new approximate boundary conditions have been applied to the large eddy simulation of channel flow with and without transpiration. These new boundary conditions give more accurate results than those previously in use, and allow significant reduction of the required CPU time over simulations in which no-slip conditions are applied. Mean velocity profiles and turbulence intensities compare well both with experimental data and with the results of resolved simulations. The influence of the approximate boundary conditions remains confined near the point of application and does not affect the turbulence statistics in the core of the flow.

Piomelli, Ugo↗