Detached Eddy Simulations (DES) is a hybrid URANS-LES technique that was proposed to obtain computationally feasible solutions of high Reynolds number flows undergoing massive separation with reliable accuracy. Since its inception, DES has been applied to a wide variety of flow fields, but mostly limited to unbounded external aerodynamic flows. This is the first study to apply and validate DES to predict the internal flow and heat transfer in non-canonical flows of industrial relevance. The prediction capabilities of DES in capturing the effects of Coriolis forces, which are induced by rotation, and centrifugal buoyancy forces, which are induced by thermal gradients, are also authenticated. The accurate prediction of turbulent flows is sensitive to the level of turbulence predicted by the turbulence scheme. By treating the regions of interest in LES mode, DES allows the unsteadiness in these regions to develop and hence predicts the turbulence levels accurately. Additionally, this permits DES to capture the effects of system rotation and buoyancy. Computations on a rotating system (a sudden expansion duct) and a system subjected to thermal gradients (cavity with a heated wall) validate the prediction capability of DES. The application of DES is further extended to a non-canonical, internal flow which is of relevance in internal cooling of gas turbine blades. Computations of the fully developed flow and heat transfer shows that DES surpasses several shortcomings of the RANS model on which it is based. DES accurately predicts the primary and secondary flow features, the turbulence characteristics and the heat transfer in stationary ducts and in rotating ducts, where the effects of Coriolis forces and centrifugal buoyancy forces are dominant. DES computations are carried out at a computational cost that is almost an order of magnitude less than the LES with little compromise on the accuracy. However, the capabilities of DES in predicting the transition to turbulence are inadequate, as highlighted by the flow features and the heat transfer in the developing region of the duct. But once the flow becomes fully turbulent, DES predicts the flow physics and shows good quantitative agreement with the experiments and LES. / Ph. D.

Increasing the turbine inlet temperature of gas turbine engines significantly increases their power output and efficiency, but it also increases the likelihood of thermal failure. Internal passages with tiny ribs are typically cast into turbine blades to cool them, and the ability to accurately predict the flow and heat transfer within these channels leads to higher design reliability and prevention of blade failure resulting from local thermal loading. Prediction of the flow through these channels is challenging, however, because the flow is highly turbulent and anisotropic, and the presence of rotational body forces further complicates the flow. Large Eddy Simulations are used to study these flows because of their ability to predict the unsteady flow effects and anisotropic turbulence more reliably than traditional RANS closure models. Calculations in a stationary duct are validated with experiments in the developing flow, fully developed, and 180° bend regions to establish the accuracy and prediction capability of the LES calculations and to aid in understanding the major flow structures encountered in a ribbed duct. It is found that most flow and heat transfer calculations come to within 10-15% of the measurements, typically showing excellent agreement in all comparisons. In the developing flow region, Coriolis effects are found to destabilize turbulence and increase heat transfer along the trailing wall (pressure side), while decreasing leading wall heat transfer by stabilizing turbulence. Coriolis forces improve flow turning in the 180° bend by shifting the shape of the separated recirculation zone at the tip of the dividing wall and increasing the mainstream flow area. In addition, turbulence is attenuated near the leading wall throughout the bend, while Coriolis forces have little effect on trailing wall turbulence in the bend. Introducing and increasing centrifugal buoyancy in the developing flow region increases trailing wall heat transfer monotonically. Along the leading wall, buoyancy increases the size of the recirculation zones, shifting the peak heat transfer to a region upstream of the rib, which decreases heat transfer at low buoyancy parameters but increases it as the buoyancy parameter is increased beyond a value of 0.3. Centrifugal buoyancy in the 180° bend initially decreases the size of the recirculation zone at the tip of the dividing wall, increasing flow area and decreasing flow impingement. At high buoyancy, however, the recirculation zone shifts to the middle of the bend, increasing flow resistance and causing strong flow impingement on the back wall. The Boussinesq approximation is used in the buoyancy calculations, but the accuracy of the approximation comes into question in the presence of large temperature differences. A variable property algorithm is developed to calculate unsteady low speed flows with large density variations resulting from large temperature differences. The algorithm is validated against two test cases: Rayleigh-Bénard convection and Poiseuille-Bénard flow. Finally, design issues in rotating ribbed ducts are considered. The fully developed assumption is discussed with regard to the developing flow region, and controlling the recirculation zone in the 180° bend is considered as a way to determine the blade tip heat transfer and pressure drop across the bend. / Ph. D.

La stabilité d’un fluide Newtonien diélectrique confiné dans un anneau cylindrique et soumis à un gradient radial de température et à un champ électrique est étudiée. Le gradient de température induit une stratification de la permittivité électrique du fluide et de sa masse volumique. Trois poussées thermiques rentrent alors en jeu : la gravité terrestre créée la poussée d’Archimède, la rotation des cylindres engendre la poussée centrifuge, et le champ électrique induit la poussée diélectrophorétique. L’effet de ces poussées est étudié dans différentes combinaisons, principalement à travers l’étude de la stabilité linéaire, mais également par la simulation numérique directe. / The stability of a Newtonian dielectric fluid confined in a cylindrical annulus and submitted to a radial temperature gradient and an electric field is studied. The temperature gradient induces a stratification of the electric permittivity and of the density. Thus three thermal buoyancies intervene: the Earth gravity creates the Archimedean buoyancy, the rotation of the cylinders generates the centrifugal buoyancy, and the electric field induces the dielectrophoretic buoyancy. The effect of these buoyancies is studied in different combination, principally through the linear stability analysis, but also by direct numerical simulation.

In this thesis, the ribbed ducts of the internal cooling passage in turbine blading are investigated to demonstrate the effects of high speed rotation. Rotation coupled with high temperature operating conditions alters the mean flow, turbulence, and heat transfer augmentation due to Coriolis and centrifugal buoyancy forces that arises from density stratification in the domain. Gas turbine engines operate in particle laden environments (sand, volcanic ash), and particulate matter ingested by the engine can make their way into the blade internal cooling passages over thousands of operating hours. These particulates can deposit on the walls of these cooling passages and degrade performance of the turbine blade. Large-Eddy Simulations (LES) with temperature dependent properties is used for turbulent flow and heat transfer in the ribbed cooling passages and Lagrangian tracking is used to calculate the particle trajectories together with a wall deposition model. The conditions used are Re=100,000, Rotation number, Ro = 0.0 and 0.2, and centrifugal Buoyancy parameters of Bo=0, 0.5, and 1.0. First, the independent effects of Coriolis and centrifugal buoyancy forces are investigated, with a focus on the additional augmentation obtained in heat transfer with the addition of centrifugal buoyancy. Coriolis forces are known to augment heat transfer at the trailing wall and attenuate the same at the leading wall. Phenomenological arguments stated that centrifugal buoyancy augments the effects of Coriolis forces in outward flow in the first pass while opposing the effect of Coriolis forces during inward flow in the second pass. In this study, it was found that in the first pass, centrifugal buoyancy had a greater effect in augmenting heat transfer at the trailing wall than in attenuating heat transfer at the leading wall. On the contrary, it aided heat transfer in the second half of the first pass at the leading wall by energizing the flow near the wall. Also, contrary to phenomenological arguments, inclusion of centrifugal buoyancy augmented heat transfer over Coriolis forces alone on both the leading and trailing walls of the second pass. Sand ingestion is then investigated, by injecting 200,000 particles in the size range of 0.5-175μm with 65% of the particles below 10 μm. Three duct wall temperatures are investigated, 950, 1000 and 1050 °C with an inlet temperature of flow and particles at 527 °C . The impingement, deposition levels, and impact characteristics are recorded as the particles move through the domain. It was found that the Coriolis force greatly increases deposition. This was made prevalent in the first pass, as 84% of the deposits in the domain occurred in the first pass for the rotating case, whereas only 27% of deposits occurred in the first pass for the stationary case with the majority of deposits occurring in the bend region. This was due to an increased interaction with the trailing wall in the rotating case whereas particles in the stationary case were allowed to remain in the mean flow and gain momentum, making rebounding from a wall during collision more likely than deposition. In contrast, the variation of wall temperatures caused little to no change in deposition levels. This was concluded to be a result of the high Reynolds number used in the flow. At high Reynolds numbers, the particles have a short residence times in the internal cooling circuit not allowing the flow and particles to heat up to the wall temperature. Overall, 87% of the injected particles deposited in the rotating duct whereas 58% deposited in the stationary duct. / Master of Science