The present thesis focuses on evaluating a conjugate heat transfer (CHT) simulation in a ribbed cooling passage with a fully developed flow assumption using LES with the immersed boundary method (IBM-LES-CHT). The IBM with the LES model (IBM-LES) and the IBM with CHT boundary condition (IBM-CHT) frameworks are validated prior to the main simulations by simulating purely convective heat transfer (iso-flux) in the ribbed duct, and a developing laminar boundary layer flow over a two-dimensional flat plate with heat conduction, respectively. For the main conjugate simulations, a ribbed duct geometry with a blockage ratio of 0.3 is simulated at a bulk Reynolds number of 10,000 with a conjugate boundary condition applied to the rib surface. The nominal Biot number is kept at 1, which is similar to the comparative experiment. As a means to overcome a large time scale disparity between the fluid and the solid regions, the use of a high artificial solid thermal diffusivity is compared to the physical diffusivity. It is shown that while the diffusivity impacts the instantaneous fluctuations in temperature, heat transfer and Nusselt numbers, it has an insignificantly small effect on the mean Nusselt number. The comparison between the IBM-LES-CHT and iso-flux simulations shows that the iso-flux case predicts higher local Nusselt numbers at the back face of the rib. Furthermore, the local Nusselt number augmentation ratio (EF) predicted by IBM-LES-CHT is compared to the body fitted grid (BFG) simulation, experiment and another LES conjugate simulation. Even though there is a mismatch between IBM-LES-CHT prediction and other studies at the front face of the rib, the area-averaged EF compares reasonably well in other regions between IBM-LES-CHT prediction and the comparative studies. / Master of Science / The present thesis focuses on the computational study of the conjugate heat transfer (CHT) investigation on the turbine internal ribbed cooling channel. Plenty of prior research on turbine internal cooling channel have been conducted by considering only the convective heat transfer at the wall, which assumes an iso-flux (constant heat flux) boundary condition at the surface. However, applying an iso-flux condition on the surface is far from the realistic heat transfer mechanism occurring in internal cooling systems. In this work, a conjugate heat transfer analysis of the cooling channel, which considers both the conduction within the solid wall and the convection at the ribbed inner wall surface, is conducted for more realistic heat transfer coefficient prediction at the inner ribbed wall. For the simulation, the computational mesh is generated by the immersed boundary method (IBM), which can ease the mesh generation by simply immersing the CAD geometry into the background volume grid. The IBM is combined with the conjugate boundary condition to simulate the internal ribbed cooling channel. The conjugate simulation is compared with the experimental data and another computational study for the validation. Even though there are some discrepancy between the IBM simulation and other comparative studies, overall results are in good agreement. From the thermal prediction comparison between the iso-flux case and the conjugate case v using the IBM, it is found that the heat transfer predicted by the conjugate case is different from the iso-flux case by more than 40 percent at the rib back face. The present study shows the potential of the IBM framework with the conjugate boundary condition for more complicated geometry, such as full turbine blade model with external and internal cooling system.

The desire to achieve higher heat transfer augmentation for turbine blades is fueled by the increased power output and efficiency that is achievable with high turbine inlet temperatures. The use of internal cooling channels fitted with pin fin arrays serves as one method of accomplishing this goal. Consequently, the addition of pin fin arrays comes at the expense of increased pressure drop. Therefore the pin fin geometry must be judiciously chosen to achieve the required heat transfer rate while minimizing the associated pressure drop. This project culminates in the measurement of both pin fin force and array pressure drop as they related to changes in the array geometry. Specifically, the effects of Reynolds number, spanwise pin spacing, streamwise pin spacing, pin aspect ratio, and flow incidence angle. Direct two-component force measurement is achieved with a cantilever beam force sensor that uses highly sensitive piezoresistive strain gauges, relating the strain at the base of the beam to the applied force. With proper characterization, forces as small as one-tenth the weight of a paper clip are successfully measured. Additionally, array pressure drop measurements are achieved using static pressure taps. Experiments were conducted over a range of Reynolds numbers between 7,500 and 35,000. Changes in the spanwise pin spacing were shown to substantially alter the pin fin drag and array pressure drop, while changes in the streamwise pin spacing were less influential. The experimental results also showed a dramatic reduction in the pin fin drag and array pressure drop for an inline flow incidence angle. Finally, changes in the pin aspect ratio were shown to have little effect on the array pressure drop. / Master of Science

Large-Eddy Simulations (LES or wall-resolved LES, WRLES) has been used extensively in capturing the physics of anisotropic turbulent flows. However, near wall turbulent scales in the inner layer in wall bounded flows makes it unfeasible for large Reynolds numbers due to grid requirements. This study evaluates the use of a wall model for LES (WMLES) on a channel with rotation at ã Reã _b = 34,000 from ã Roã _b = 0 to 0.38, non-staggered 90° ribbed duct with rotation at ã Reã _b = 20,000 from ã Roã _b = 0 to 0.70, stationary 45° staggered ribbed duct at ã Reã _b = 49,000, and two-pass smooth duct with a U-bend at ã Reã _b = 25,000 for ã Roã _b = 0 to 0.238 against WRLES and experimental data. In addition, for the two-pass smooth duct with a U-bend simulations, the synthetic eddy method (SEM) is used to artificially generate eddies at the inlet based on given flow characteristics. It is presented that WMLES captures the effects of Coriolis forces and predicts mean heat transfer augmentation ratios reasonably well for all simulations. The alleviated grid resolution for these simulations indicates significant reductions in resources, specifically, by a factor of 10-20 in non-staggered 90° ribbed duct simulations. The combined effects of density ratio, Coriolis forces, with SEM for the inlet turbulence, capture the general trends in heat transfer in and after the bend. / Master of Science

The desire for smaller, faster, and more efficient products places a strain on thermal management in components ranging from gas turbine blades to computers. Heat exchangers that utilize internal cooling flows have shown promise in both of these industries. Although pin fins are often placed in the cooling channels to augment heat transfer, their addition comes at the expense of increased pressure drop. Consequently, the pin fin geometry must be judiciously chosen to achieve the desired heat transfer rate while minimizing the pressure drop and accompanying pumping requirements. This project culminates in the construction of a new test facility and the development of a unique force measurement methodology. Direct force measurement is achieved with a cantilever beam force sensor that uses sensitive piezoresistive strain gauges to simultaneously measure aerodynamic lift and drag forces on a pin fin. After eliminating the detrimental environmental influences, forces as small as one-tenth the weight of a paper clip are successfully measured. Although the drag of an infinitely long cylinder in uniform cross flow is well documented, the literature does not discuss the aerodynamic forces on a cylinder with an aspect ratio of unity in channel flow. Measured results indicate that the drag coefficient of a cylindrical pin in a single row array is greater than the drag coefficient of an infinite cylinder in cross flow. This phenomenon is believed to be caused by an augmentation of viscous drag on the pin fin induced by the increased viscous effects inherent in channel flow. / Master of Science

The heat transfer characteristics through arrays of pin fins were studied for the further development of internal cooling methods for turbine airfoils. Low aspect ratio pin fin arrays were tested through a range of Reynolds numbers between 5000 and 30,000 to determine the effects of pin spacing as well as aspect ratio on pin and endwall heat transfer. Experiments were also conducted to determine the independent effects of pin spacing and aspect ratio on arrays with different flow incidence angles. The pin Nusselt numbers showed almost no dependence on pin spacing or flow incidence angle. Using an infrared thermogaphy technique, spatially-resolved Nusselt numbers were measured along the endwalls of each array. The endwall results showed that streamwise spacing had a larger effect than spanwise spacing on array-averaged Nusselt numbers. Endwall heat transfer patterns showed that arrays with flow incidence angles experienced less wake interaction between pins than arrays with perpendicular flow, which caused a slight decrease in heat transfer in arrays with flow incidence angles. The effect of flow incidence angle on array-average Nusselt number was greater at tighter pin spacings. Even though the pin Nusselt number was independent of pin spacing, the ratio of pin-to-endwall Nusselt number was dependent on flow conditions as well as pin spacing. The pin aspect ratio had little effect on the array-average Nusselt number for arrays with perpendicular flow; however, the effect of flow incidence angle on array-average Nusselt number increased as aspect ratio decreased. / Master of Science

The performance of many engineering devices from power electronics to gas turbines is limited by thermal management. Pin fins are commonly used to augment heat transfer by increasing surface area and increasing turbulence. The present research is focused on but not limited to internal cooling of turbine airfoils using pin fins. Although the pin fins are not limited to a single shape, circular cross-sections are most common. The present study examines heat transfer from a single row of circular pin fins with the row oriented perpendicular to the flow. The configurations studied have spanwise spacing to pin diameter ratios of two, four, and eight. Low aspect ratio pin fins were studied whereby the channel height to pin diameter was unity. The experiments are carried out for a Reynolds number range of 5000 to 30,000. Heat transfer measurements are taken on both the pin and on the endwall covering several pin diameters upstream and downstream of the pin row. The results show that the heat transfer augmentation relative to open channel flow is highest for the smallest spanwise spacing for the lowest Reynolds number flows. The results also indicate that the pin fin heat transfer is higher than on the endwall. / Master of Science

Gas turbines for military applications, when operating in harsh environments like deserts often encounter unexpected operation faults. Such performance deterioration of the gas turbine decreases the mission readiness of the Air Force and simultaneously increases the maintenance costs. Some of the major factors responsible for the reduced performance are ingestion of debris during take off and landing, distorted intake flows during low altitude maneuvers, and hot gas ingestion during artillery firing. The focus of this thesis is to study ingestion of debris; specifically sand. The region of interest being the internal cooling ribbed duct of the turbine blade. The presence of serpentine passages and strong localized cross flow components makes this region prone to deposition, erosion, and corrosion (DEC) by sand particles. A Lagrangian particle tracking technique was implemented in a generalized coordinate multi-block Navier-Stokes solver in a distributed parallel framework. The developed algorithm was validated by comparing the computed particle statistics for 28 microns lycopodium, 50 microns glass, and 70 microns copper with available data [2] for a turbulent channel flow at Ret=180. Computations were performed for a particle-laden turbulent flow through a stationary ribbed square duct (rib pitch / rib height = 10, rib height / hydraulic diameter = 0.1) using an Eulerian-Lagrangian framework. Particle sizes of 10, 50, and 100 microns with response times (normalized by friction velocity and hydraulic diameter) of 0.06875, 1.71875, and 6.875 respectively are considered. The calculations are performed for a nominal bulk Reynolds number of 20,000 under fully developed conditions. The carrier phase was solved using Large Eddy Simulation (LES) with Dynamic Smagorinsky Model [1]. Due to low volume fraction of the particles, one-way fluid-particle coupling was assumed. It is found that at any given instant in time about 40% of the total number of 10 micron particles are concentrated in the vicinity (within 0.05 Dh) of the duct surfaces, compared to 26% of the 50 and 100 micron particles. The 10 micron particles are more sensitive to the flow features and are prone to preferential concentration more so than the larger particles. At the side walls of the duct, the 10 micron particles exhibit a high potential to erode the region in the vicinity of the rib due to secondary flow impingement. The larger particles are more prone to eroding the area between the ribs and towards the center of the duct. At the ribbed walls, while the 10 micron particles exhibit a fairly uniform propensity for erosion, the 100 micron particles show a much higher tendency to erode the surface in the vicinity of the reattachment region. The rib face facing the flow is by far the most susceptible to erosion and deposition for all particle sizes. While the top of the rib does not exhibit a large propensity to be eroded, the back of the rib is as susceptible as the other duct surfaces because of particles which are entrained into the recirculation zone behind the rib. / Master of Science

Gas turbine blades operate in very high temperatures to achieve a high thermal efficiency of the engine. For this reason, the blades have to be cooled to prevent degradation or even melting. The blades can be cooled using various techniques, both by cooling the inside of the blade with cooling channels, and by protecting the outside of the blade from the hot environment. One way to cool the blades from the inside is with rib turbulated channels. Straight square channels lined with 90◦, 45◦ and V-shaped ribs in a staggered configuration are investigated in this thesis.  Computational fluid dynamics (CFD), among other methods, can be used to predict important parameters such as heat transfer and pressure loss for different ribbed channel geometries. In this thesis a CFD model using RANS simulations with the turbulence models Lag Elliptic Blending k − ε, Realizable k − ε two-layer and SST k − ω is established and validated against experimental data by Taslim et al [1]. This is done by comparing the Nusselt number between a pair of ribs as well as the channel friction factor for 90◦, 45◦ and V-shape ribs. Different sensitivities are also investigated to get an understanding of the uncertainties found during the CFD implementation. These include the effect of mesh resolution, inlet turbulence intensity, rounded rib edges, wall roughness and temperature used for Reynolds number calculations. The Nusselt number and friction factor predictions of the turbulence models are also compared with existing empirical correlations.  The results of the investigation show that the CFD results for 90◦ ribs deviate the most from experimental results, while closer results are seen for the 45◦ and V-shape ribs.  In conclusion, the Lag Elliptic Blending k−ε model generally produces results closest to experimental data, especially for 90◦ ribs, but it shows some differences in Reynolds number trends. It proves to predict heat transfer and pressure loss closer to the experiment than the other models in flows where recirculation and reattachment has a significant impact. The Lag EB model is relatively stable and mesh independent. The SST k − ω model produces results rather similar to experimental data, but is unstable and sensitive to mesh resolution. The Realizable k − ε two-layer model produces results that are slightly less consistent with experimental data, but is very stable and insensitive to mesh resolution. The Nusselt number and friction factor from the investigated empirical correlations are closer to experimental results than the turbulence models for 90◦ inline ribs.

Submitted by Renan Luis Fragelli null (renan.fragelli@gmail.com) on 2017-03-28T15:54:09Z No. of bitstreams: 1 Dissertação de Mestrado - Renan Fragelli.pdf: 5490064 bytes, checksum: 80c3b38563331cc9ef7e88ff192c5c8b (MD5) / Approved for entry into archive by Luiz Galeffi (luizgaleffi@gmail.com) on 2017-03-29T20:46:29Z (GMT) No. of bitstreams: 1 fragelli_rl_me_bauru.pdf: 5490064 bytes, checksum: 80c3b38563331cc9ef7e88ff192c5c8b (MD5) / Made available in DSpace on 2017-03-29T20:46:29Z (GMT). No. of bitstreams: 1 fragelli_rl_me_bauru.pdf: 5490064 bytes, checksum: 80c3b38563331cc9ef7e88ff192c5c8b (MD5) Previous issue date: 2017-03-16 / Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES) / Este trabalho surgiu a partir da necessidade de produzir avanços em projeto que trata do desenvolvimento de um porta-ferramentas refrigerado internamente através de um fluido em mudança de fase e, na tentativa de minimizar a alta temperatura na ferramenta de corte através desse sistema de circulação. A utilização de nanofluidos surgiu como uma alternativa para a otimização da transferência térmica entre fluido e ferramenta de corte. A pesquisa consiste em avaliar a influência da adição de nanopartículas de prata numa solução de etilenoglicol e água deionizada, e também, da adição de nanopartículas de hematita (Fe2O3) no fluido refrigerante R141b. Em ambos os casos, as nanopartículas possuíam formato esférico, diâmetro médio de 30nm e foram avaliadas em concentrações. Além disso, as duas soluções foram submetidas a um campo elétrico na região de transferência térmica para analisar a influência do efeito eletrohidrodinâmico e, por fim, considerando as propriedades magnéticas da hematita, este nanofluido foi testado sob influência de um campo magnético. Os testes mostraram que as nanopartículas realmente influenciaram as propriedades dos fluidos e, por consequência, a quantidade de calor transferido. O nanofluido Ag/ETG+H2O(l) (0,023 vol%) resultou num incremento de 11% no valor do coeficiente de transferência térmica convectivo (h) quando sujeito ao campo elétrico. Para o caso do nanofluido Fe2O3/R141b, o valor de h aumentou em 30,3%, porém, quando sob efeito do campo magnético ou elétrico, o coeficiente foi prejudicado, resultando num valor menor que o do controle. Ao final, tem-se a proposta de um possível modelo desse porta-ferramentas. / This work arose from the need to produce advances in design development of an internally cooled toolholder through a phase change fluid. In order to minimize the high temperature in the cutting tool by this circulation system, using nanofluids emerged as an alternative to optimize heat transfer between the fluid and the cutting tool. The research consists in evaluate the influence of addition of silver nanoparticles in an ethylene glycol and deionized water solution, and also the addition of hematite nanoparticles (Fe2O3) in the refrigerant R141b. In both cases, nanoparticles had spherical shape, diameter of 30nm, and they were evaluated in different concentrations. Moreover, both nanofluids were subjected to an electric field in the heat transfer region to evaluate the influence of electrohydrodynamic effect and, finally, considering the magnetic properties of hematite, this nanofluid was tested under the influence of a magnetic field. The tests have shown that the nanoparticles really influence the properties of the fluids and, therefore, the amount of heat transferred. The nanofluid Ag/ETG+H2O(l) also presented a positive influence of the electric field, further enhancing the value of the convective heat transfer coefficient (h) in 11% (0,039 vol%). In the case of Fe2O3/R141b nanofluid, the h value increased 30.3%. However, when the nanofluid was under magnetic or electric effect, the value of h was deteriorated, resulting in a lesser value than the control. As conclusion, a new toolholder prototype is presented.

No description available.