Engineering topics

Meneveau, C.

Publications and source records attributed to Meneveau, C..

7 records

NASA NTRS1996

Dynamic Smagorinsky model on anisotropic grids

Large Eddy Simulation (LES) of complex-geometry flows often involves highly anisotropic meshes. To examine the performance of the dynamic Smagorinsky model in a controlled fashion on such grids, simulations of forced isotropic turbulence are performed using highly anisotropic discretizations. The resulting model coefficients are compared with a theoretical prediction (Scotti et al., 1993). Two extreme cases are considered: pancake-like grids, for which two directions are poorly resolved compared to the third, and pencil-like grids, where one direction is poorly resolved when compared to the other two. For pancake-like grids the dynamic model yields the results expected from the theory (increasing coefficient with increasing aspect ratio), whereas for pencil-like grids the dynamic model does not agree with the theoretical prediction (with detrimental effects only on smallest resolved scales). A possible explanation of the departure is attempted, and it is shown that the problem may be circumvented by using an isotropic test-filter at larger scales. Overall, all models considered give good large-scale results, confirming the general robustness of the dynamic and eddy-viscosity models. But in all cases, the predictions were poor for scales smaller than that of the worst resolved direction.

Scotti, A.↗

NASA NTRS1996

Dynamic model with scale-dependent coefficients in the viscous range

The standard dynamic procedure is based on the scale-invariance assumption that the model coefficient C is the same at the grid and test-filter levels. In many applications this condition is not met, e.g. when the filter-length, delta, approaches the Kolmogorov scale, and C(delta approaches eta) approaches O. Using a priori tests, we show that the standard dynamic model yields the coefficient corresponding to the test-filter scale (alpha delta) instead of the grid-scale (delta). Several approaches to account for scale dependence are examined and/or tested in large eddy simulation of isotropic turbulence: (a) take the limit alpha approaches 1; (b) solve for two unknown coefficients C(Delta) and C(alpha delta) in the least-square-error formulation; (c) the 'bi-dynamic model', in which two test-filters (e.g. at scales 2(delta) and 4(delta) are employed to gain additional information on possible scale-dependence of the coefficient, and an improved estimate for the grid-level coefficient is obtained by extrapolation, (d) use theoretical predictions for the ratio C(alpha delta)/C(delta) and dynamically solve for C(delta). None of these options is found to be entirely satisfactory, although the last approach appears applicable to the viscous range.

Meneveau, C.↗

NASA NTRS1994

A Lagrangian dynamic subgrid-scale model turbulence

A new formulation of the dynamic subgrid-scale model is tested in which the error associated with the Germano identity is minimized over flow pathlines rather than over directions of statistical homogeneity. This procedure allows the application of the dynamic model with averaging to flows in complex geometries that do not possess homogeneous directions. The characteristic Lagrangian time scale over which the averaging is performed is chosen such that the model is purely dissipative, guaranteeing numerical stability when coupled with the Smagorinsky model. The formulation is tested successfully in forced and decaying isotropic turbulence and in fully developed and transitional channel flow. In homogeneous flows, the results are similar to those of the volume-averaged dynamic model, while in channel flow, the predictions are superior to those of the plane-averaged dynamic model. The relationship between the averaged terms in the model and vortical structures (worms) that appear in the LES is investigated. Computational overhead is kept small (about 10 percent above the CPU requirements of the volume or plane-averaged dynamic model) by using an approximate scheme to advance the Lagrangian tracking through first-order Euler time integration and linear interpolation in space.

Meneveau, C.↗

NASA NTRS1992

On the local nature of the energy cascade

The local nature of the energy cascade in space and time is studied using direct numerical simulation of decaying and forced isotropic turbulence. To examine the concept that large scales evolve into smaller ones, we compute the Lagrangian correlation coefficient between local kinetic energy at different scales. This correlation is found to peak at a Lagrangian time-delay that increases with scale separation. The results show that, on average, the flow of energy to smaller scales is predominantly local in physical space and that the view of eddies decaying into smaller ones while transferring their kinetic energy appears to be, on average, quite realistic. To examine the spectral characteristics of the cascade under unsteady conditions, a pulse of large-scale energy is added to the large-eddy simulation of forced isotropic turbulence. As time progresses, the evolution of this pulse through bands of increasing wavenumbers is studied.

Meneveau, C.↗

NASA NTRS1992

Search for subgrid scale parameterization by projection pursuit regression

The dependence of subgrid-scale stresses on variables of the resolved field is studied using direct numerical simulations of isotropic turbulence, homogeneous shear flow, and channel flow. The projection pursuit algorithm, a promising new regression tool for high-dimensional data, is used to systematically search through a large collection of resolved variables, such as components of the strain rate, vorticity, velocity gradients at neighboring grid points, etc. For the case of isotropic turbulence, the search algorithm recovers the linear dependence on the rate of strain (which is necessary to transfer energy to subgrid scales) but is unable to determine any other more complex relationship. For shear flows, however, new systematic relations beyond eddy viscosity are found. For the homogeneous shear flow, the results suggest that products of the mean rotation rate tensor with both the fluctuating strain rate and fluctuating rotation rate tensors are important quantities in parameterizing the subgrid-scale stresses. A model incorporating these terms is proposed. When evaluated with direct numerical simulation data, this model significantly increases the correlation between the modeled and exact stresses, as compared with the Smagorinsky model. In the case of channel flow, the stresses are found to correlate with products of the fluctuating strain and rotation rate tensors. The mean rates of rotation or strain do not appear to be important in this case, and the model determined for homogeneous shear flow does not perform well when tested with channel flow data. Many questions remain about the physical mechanisms underlying these findings, about possible Reynolds number dependence, and, given the low level of correlations, about their impact on modeling. Nevertheless, demonstration of the existence of causal relations between sgs stresses and large-scale characteristics of turbulent shear flows, in addition to those necessary for energy transfer, provides important insight into the relation between scales in turbulent flows.

Meneveau, C.↗

NASA NTRS1990

Continuous wavelet analysis of coherent structures

We perform an analysis of planar cuts through three dimensional turbulent fields (planar channel flow and mixing layer) using the 2D continuous wavelet transform. We propose two new diagnostics: (1) a measure of intermittency I(r, vector x), which is the ratio of local energy and average energy at a given scale r; and (2) a local Reynolds number, defined on the local velocity contribution at a given scale, computed from the wavelet transform of the three velocity components, the scale of the transform, and molecular viscosity; this gives a representation of the local non-linearity of the flow viewed in both space and scale. We find, for the analyzed flows, strong small-scale intermittency located in the ejection regions for the channel flow and in the vortex core of the mixing layer.

Farge, M.↗

NASA NTRS1990

Turbulence dynamics in the wavelet representation

The phenomenon of small-scale intermittency is shown to motivate the decomposition of the velocity fields into modes that exhibit both localization in wavenumber and physical space. We review some basic properties of such a decomposition, called the wavelet transform. The wavelet-transformed Navier-Stokes equations are derived, and we define a new quantity Pi(r, vector-x, t), which is the flux of kinetic energy to scales smaller than r at position vector-x (at time t). The main goals of this research are also summarized.

Meneveau, C.↗