Turbulent flows present structures with a wide range of scales. The computation of the complete physics of a turbulent flow (termed DNS) is very expensive and is, for the time being, limited to low and medium Reynolds number flows. As a way to capture high Reynolds number flows, a part of the physics complexity has to be modeled. Large eddy simulation (LES) is a simulation strategy where the large turbulent eddies present on a given mesh are captured and the influence of the non-resolved scales onto the resolved ones is modeled. The present thesis reports on the development and validation of a methodology in order to apply LES for complex wall-bounded flows. Discretization methods and LES models, termed subgrid scale models (SGS), compatible with such a geometrical complexity are discussed. It is proved that discrete a kinetic energy conserving discretization of the convective term is an attractive solution to perform stable simulations without the use of an artificial dissipation, as upwinding. The dissipative effect of the SGS model is thus unaffected by any additional dissipation process. The methodology is first applied to a developed parallel fourth-order incompressible flow solver for cartesian non-uniform meshes. In order to solve the resulting Poisson equation, an efficient multigrid solver is also developed. The code is first validated using DNS (Taylor-Green vortex, channel flow, four-vortex system) and LES (channel flow), and finally applied to the investigation of an aircraft two-vortex system in ground effect. The methodology is then applied to improve a RANS-based industrial unstructured compressible flow solver, developed at CENAERO, to perform well for LES applications. The proposed modifications are tested successfully on the unsteady flow past a sphere at Reynolds of 300 and 10000, corresponding to the subcritical regime.

No description available.

The average and fluctuating components in a model ring-type gas turbine combustor are characterized using a Large Eddy Simulation at a Reynolds number of 11,000, based on the bulk velocity and the mean channel height. A spatial filter is applied to the incompressible Navier-Stokes equations, and a high pass filtered Smagorinsky model is used to model the sub-grid scales. Two cases are studied: one with only the swirler inlet active, and one with a single row of dilution jets activated, operating at a momentum flux ratio J of 100. The goal of both of these studies is to validate the capabilities of the solver NEK5000 to resolve important flow features inherent to gas turbine combustors by comparing qualitatively to the work of Jakirlic. Both cases show strong evidence of the Precessing Vortex Core, an essential flow feature in gas turbine combustors. Each case captures other important flow characteristics, such as corner eddies, and in general predicts bulk flow movements well. However, the simulations performed quite poorly in terms of predicting turbulence shear stress quantities. Difficulties in properly emulating the turbulent velocity entering the combustor for the swirl, as well as mesh quality concerns, may have skewed the results. Overall, though small length scale quantities were not accurately captured, the large scale quantities were, and this stress test on the HPF LES model will be built upon in future work that looks at more complex combustors.

In the context of Large Eddy Simulation (LES) of turbulent flows, there is a current need to compare and evaluate different proposed subfilter-scale models. In order to carefully compare subfilter-scale models and compare LES predictions to Direct Numerical Simulation (DNS) results (the latter would be helpful in the comparison and validation of models), there is a real need for a "grid-independent" LES capability and explicit filtering methods offer one means by which this may be achieved. Advantages of explicit filtering are that it provides a means for eliminating aliasing errors, allows for the direct control of commutation errors, and most importantly allows a decoupling between the mesh spacing and the filter width which is the primary reason why there are difficulties in comparing LES solutions obtained on different grids. This thesis considers the design and assessment of discrete explicit filters and their application to isotropic turbulence prediction.

In the context of Large Eddy Simulation (LES) of turbulent flows, there is a current need to compare and evaluate different proposed subfilter-scale models. In order to carefully compare subfilter-scale models and compare LES predictions to Direct Numerical Simulation (DNS) results (the latter would be helpful in the comparison and validation of models), there is a real need for a "grid-independent" LES capability and explicit filtering methods offer one means by which this may be achieved. Advantages of explicit filtering are that it provides a means for eliminating aliasing errors, allows for the direct control of commutation errors, and most importantly allows a decoupling between the mesh spacing and the filter width which is the primary reason why there are difficulties in comparing LES solutions obtained on different grids. This thesis considers the design and assessment of discrete explicit filters and their application to isotropic turbulence prediction.

Large Eddy Simulation is a method of obtaining high accuracy computational results for modelling fluid flow. Unfortunately it is computationally expensive limiting it to users of large parallel machines. However, it may be that the use of LES leads to an over-resolution of the problem because the bulk of the computational domain could be adequately modelled using the Reynolds averaged approach. A study has been undertaken to assess the feasibility, both in accuracy and computational efficiency of using a parallel computer to solve both LES and RANS type turbulence models on the same domain for the problem flow over a circular cylinder at Reynolds number 3 900 To do this the domain has been created and then divided into two sub-domains, one for the LES model and one for the kappa - epsilon turbulence model. The hybrid model has been developed specifically for a parallel computing environment and the user is able to allocate modelling techniques to processors in a way which enables expansion of the model to any number of processors. Computational experimentation has shown that the combination of the Smagorinsky model can be used to capture the vortex shedding from the cylinder and the information successfully passed to the kappa - epsilon model for the dissipation of the vortices further downstream. The results have been compared to high accuracy LES results and with both kappa - epsilon and Smagorinsky LES computations on the same domain. The hybrid models developed compare well with the Smagorinsky model capturing the vortex shedding with the correct periodicity. Suggestions for future work have been made to develop this idea further, and to investigate the possibility of using the technology for the modelling of mixing and fast chemical reactions based on the more accurate prediction of the turbulence levels in the LES sub-domain.

A numerical study of the flow field generated by injection through diamondshaped orifices into a high-speed flow is presented in this document. Jet interaction flows have a wide range of applications in the field of engineering. These applications include the use of jets for fuel injection in scramjets, for reaction control of high-speed aerodynamic bodies and as cooling jets for skins of high-speed vehicles. A necessary requirement in the use of transverse jets for these and other applications is a thorough understanding of the physics of the interaction between the jet and freestream. This interaction generates numerous flow structures that include multiple shocks, vortices, recirculation regions and shear layers. This study involves diamond-shaped orifices that have the advantage of generating weaker or attached interaction shocks as compared to circular injectors. These injectors also negate the effects due to the recirculation region that is formed upstream of the injector. This study was undertaken in order to gain further understanding of the flow features generated by diamond-shaped injectors in a high-speed flow. Numerical simulations were performed using two different levels of turbulence models. Reynolds™ Averaged Navier-Stokes (RANS) simulations were performed using the GASP flow solver while Detached-Eddy Simulation (DES) runs were performed using the Cobalt flow solver. A total of fifteen diamond injector simulations were performed using the RANS model for a 15 half-angle diamond injector. The fifteen simulations spanned over five different injection angles and three jet total pressures. In addition to these, two circular injector simulations were also performed. In addition, low pressure normal injection through diamond and circular orifices simulations were performed using DES. Results obtained from CFD were compared to available experimental data. The resulting flow structure and the turbulent properties of the flow were examined in detail. The normal injection case through the diamond-shaped orifice at the lowest jet total pressure was defined as the baseline case and is presented in detail. In order to study the effect of different components of the vorticity transport equation, an in-house code was used post-process the results from the RANS runs.

No description available.

<p>The aim of the present study is to develop subgrid-scale models that are relevant for complex flows and combustion. A stochastic model based on a stochastic Smagorinsky constant with adjustable variance and time scale is proposed. The stochastic model is shown to provide for backscatter of both kinetic energy and scalar variance without causing numerical instabilities. A new subgrid-scale scalar flux model is developed using the same kind of methodology that leads to the explicit algebraic scalar flux model, EASFM, for RANS. The new model predicts the anisotropy of the subgrid-scales in a more realistic way than the eddy diffusion model. Both new models were tested in rotating homogeneous shear flow with a passive scalar. Rogallo’s method of moving the frame with the mean flow to enable periodic boundary conditions was used to simulate homogeneous shear flow.</p>

Optimal large eddy simulation (OLES) is an approach to LES sub-grid modeling that requires multi-point correlation data as input. Until now, this has been obtained by analyzing DNS statistics. In the finite-volume OLES formulation studied here, under the assumption of small-scale homogeneity and isotropy, these correlations can be theoretically determined from Kolmogorov inertial-range theory, small-scale isotropy, along with the quasi-normal approximation. These models are expressed as generalized quadratic and linear finite volume operators that represent the convective momentum flux. These finite volume operators have been analyzed to determine their characteristics as numerical approximation operators and as models of small-scale effects. In addition, the dependence of the model operators on the anisotropy of the grid and on the size of the stencils is analyzed to develop idealized general operators that can be used on general grids. The finite volume turbulence operators developed here will be applicable in a wide range of LES problems. / text