No description available.

The objective of this study is to develop a CFD tool for performing reliable large eddy simulation (LES) of the compressible evaporating two-phase turbulent flow in a gas turbine combustor. The KIVA-3V code originally developed by Los Alamos National Laboratory is used as a baseline code. The KIVA-3V code has been modified to facilitate LES calculations. Both the temporal and spatial accuracies of the original KIVA-3V code have been improved to second order. A one-equation subgrid scale (SGS) turbulence model is implemented to describe the unresolved turbulent subgrid effect. To ensure that there are sufficient particle numbers to capture the dynamic droplet dispersion process, the ETAB breakup model coupled with a new hybrid droplet-particle algorithm is also implemented into the code. Furthermore, the effect of the subgrid scale (SGS) velocity on the droplet dispersion is included. The SGS velocity is computed from the subgrid turbulent kinetic energy predicted by the one-equation SGS turbulence model. A new collision model based on the concept of "particle cloud" is proposed and implemented in the code. The new model greatly reduces the grid-dependence of the original O'Rourke model in a Cartesian mesh. The gas solver of the new LES version of KIVA-3V code, which will be referred as KIVA-LES hereby) is validated against large eddy simulations of natural and forced plane impinging jets. Predictions were carried out for different inflow conditions, which include a natural plane impinging jet with a random perturbation on the inflow plane and a forced plane impinging jet with a Strouhal number of 0.36, locked both in phase and laterally in space. The first simulation was performed to quantitatively study the mean flow and turbulence statistics. The computed field variables and turbulence intensity of streamwise velocity agreed well with the experimental results. The second simulation was performed to study the vortex structures of a forced plane impinging jet. The predictions captured the typical vortex structures of this kind of flow, such as spanwise rollers, successive ribs, cross ribs and wall ribs were reproduced by the simulation, which were also previously detected by the experiment of Sakakibara et al. (103) with digital particle image velocimetry (DPIV) system, but to our -best knowledge never wholly reproduced by numerical simulations to date. Moreover, the study has also led to some new findings related to the formation and evolution of successive ribs, cross ribs and wall ribs. The new collision model is tested against analytical solutions of simplified realistic collision problems in a box volume. The grid-dependence of the model is also checked against some spray test cases. The new collision scheme is computationally more efficient than the frequently used O'Rourke's (87) scheme since it abandons a sampling procedure to compute the collision number. The new model delivers sufficient accuracy in calculating the collision numbers in cases with uniformly distributed droplets although O'Rourke's model seems to perform better for these scenarios. However, for the prediction of a real spray in Cartesian gird, the new model has delivered much improved results. The predictions of the new model do not show any grid-dependent artefacts. KIVA-LES with the Lagrangian spray models is used to predict non-evaporating and evaporating diesel fuel sprays. The computed results are compared with the experimental data by Hiroyasu and Kadota (55) and Naber and Siebers (81), as well as the predictions of the original KIVA-3V. The predictions are in good agreement with the data. The large scale vortical structures are reproduced by the LES simulations, which cause "branch-like" spray shape and influence the spray penetration depth. The predictions have also captured the differences between the dense and dilute regions of the sprays. The LES analysis of diesel sprays has also demonstrated that SGS velocity has significant influence on the predicted spray angles. Most importantly, grid-convergent results, which were difficult to obtain with the original KIVA-3V, have been obtained in the present study. Finally, the validated code is used to study evaporating two-phase spray flow in a coaxial gas turbine model combustor. The predictions were compared with some published experimental data. This is a first step towards a more comprehensive numerical analysis of practical industrial combustors where multiple inlets and more complex combustor geometry are encountered. Good agreement with the data is achieved. The predictions have captured the "ring-like" vortex just downstream the annulus and "worm-like" streamwise vortical structure further downstream. The axial droplet mass flux and Sauter mean radius (SMR) are well predicted. Overall the present study has demonstrated the capability of KIVA-LES with the newly developed collision model to provide reasonably accurate predictions of evaporating two-phase flows in coaxial gas turbine combustors.

One microfilm reel (16 mm.) in pocket / Bibliography: leaves 194-201 / xiii, 212, 6 leaves, [9] leaves of plates : ill. (some col.) ; 30 cm. / Title page, contents and abstract only. The complete thesis in print form is available from the University Library. / Thesis (Ph.D.)--University of Adelaide, Dept. of Mechanical Engineering, 1987

The present research focuses on the development and computer implementation of a novel threedimensional, anisotropic turbulence model not only capable of handling complex geometries but also the turbulence driven secondary currents. The model equations comprise advanced algebraic Reynolds stress models in conjunction with Reynolds Averaged Navier-Stokes equations. In order to tackle the complex geometry of compound meandering channels, the body-fitted orthogonal coordinate system is used. The finite volume method with collocated arrangement of variables is used for discretization of the governing equations. Pressurevelocity coupling is achieved by the standard iterative SIMPLE algorithm. A central differencing scheme and upwind differencing scheme are implemented for approximation of diffusive and convective fluxes on the control volume faces respectively. A set of algebraic equations, derived after discretization, are solved with help of Stones implicit matrix solver. The model is validated against standard benchmarks on simple and compound straight channels. For the case of compound meandering channels with varying sinuosity and floodplain height, the model results are compared with the published experimental data. It is found that the present method is able to predict the mean velocity distribution, pressure and secondary flow circulations with reasonably good accuracy. In terms of engineering applications, the model is also tested to understand the importance of turbulence driven secondary currents in slightly curved channel. The development of this unique model has opened many avenues of future research such as flood risk management, the effects of trees near the bank on the flow mechanisms and prediction of pollutant transport.

This thesis considers the advantages of incorporating laminar kinetic energy modelling into turbulence modelling, in order to predict laminar-turbulent transition. The final aim is to implement an improved transition model into the industrial Finite-Volume code, Code Saturne. The literature review suggests that in order for a RANS-based model to predict transition accurately, modelling of complex, anisotropic phenomena is necessary. The Walters-Cokljat model is shown to compare very well to other transition modelling methods, including correlation-based modelling. The Walters-Cokljat model is a single-point RANS-based model that solves an additional transport equation for laminar kinetic energy. This transition model is especially desirable from an industrial stand-point, due to its single-point RANS basis, with only 3 transport equations. Although this method shows great promise as an industrial tool for transition prediction, results presented here show that there are aspects of the model that require modification. The definition of effective length-scale and the method of accounting for the effects of shear sheltering are the two main areas for consideration. The current definition of effective length-scale is found to be inappropriate for flows with large free-stream length-scales, which are common-place in turbomachinery applications. Another phenomenon commonly found in turbomachinery is separation-induced transition; however, the current function for shear sheltering effects inhibits transition when turbulence intensity is not the forcing factor. Additionally, when reviewed analytically, the definition and placement of the shear sheltering function does not match the observations of Jacobs and Durbin. Alternatives for the definitions of the effective length-scale and the shear sheltering function are proposed. The individual proposals are tested, and steps towards a full working implementation are documented.