Ribbed steel pipes made of Type 2 aluminized steel are commonly used for culvert pipes for highway drainage. Typically aluminized steel pipes have shown good durability and are expected to have long service life, e.g. 75 years; also, they are used in a wide variety of soil and water conditions. However, early corrosion of aluminized steel pipes has been recently observed in some inland locations. Initial observations showed severe corrosion in forms of pits, both along the ribs and at the nearby flat portions of the pipes. It is critical to determine the cause of early deterioration and establish methods of durability prediction. The possibility of unusual environmental conditions is being investigated elsewhere, but this research focuses on possible mechanical factors aggravating corrosion, since it is prevalent near pipe rib deformations. While forming the rib bends in the pipe, the outer bend surface is exposed to extreme tensile stresses which would cause small coating cracks (microfissures) exposing base metal. Those may lead to early corrosion as galvanic protection from the surrounding aluminum may not be sufficient under certain environments. Electrochemical impedance spectroscopy was used to measure corrosion rate of both formed and flat aluminized steel samples in simulated natural waters. Initial findings show that specimens formed by spherical indentation were susceptible to early corrosion development in moderately aggressive simulated natural water, but not in a more benign, precipitating simulated natural water solution.

Les concepteurs de moteurs aéronautiques sont constamment sujets à la demande d’augmentation de puissance de la part des constructeurs d’aéronefs. Pour satisfaire à cette exigence, la température de sortie de la chambre de combustion peut être augmentée pour améliorer le rendement et la puissance de sortie du moteur. Cette élévation de température peut toutefois dépasser le point de fusion du matériau et, pour éviter les pannes de moteur, l’intégrité des aubes de la turbine repose notamment sur des systèmes de refroidissement internes,prélevant de l'air froid du compresseur. La conception de ces systèmes revient donc à maximiser l’amélioration du transfert de chaleur tout en minimisant le débit d’air via les pertes de charge afin d’éviter des pénalités de puissance du moteur. Or ces écoulements en canaux internes sont encore largement incontrôlés et mal compris. Dans le but de mieux comprendre ces écoulements en rotation se développant spatialement, ce travail porte sur l’étude via simulations numériques d’un canal de refroidissement droit, perturbé, en rotation. La configuration consiste en un canal carré équipé de 8 perturbateurs placés avec un angle de 90 degrés par rapport à l’écoulement principal. Pour les cas étudiés, des mesures PIV temporelles ont été effectuées à l'Institut VanKarman (VKI). Les conditions adiabatiques et isothermes ont été étudiées pour évaluer l’impact dela température de la paroi sur l’écoulement, en particulier dans les configurations en rotation. Les canaux statiques ainsi qu’en rotation positive et négative sont comparés avec, dans chaque cas,une prédiction d’écoulement adiabatique ou isotherme. Dans ce travail, les résultats de simulations aux grandes échelles (SGE) montrent que le modèle CFD haute fidélité est capable de reproduire les différences induites par la flottabilité sur la topologie de l'écoulement dans la région proche. Le modèle parvient également à prévoir l'augmentation (la diminution) de la turbulence autour des perturbateurs en rotation déstabilisante (stabilisante). Enfin et grâce à la SGE spatiale et temporelle complète, le développement spatial et l’instationnarité des écoulements secondaires sont analysés pour mieux comprendre leur origine et leurs différences potentielles entre les cas. Cette étude montre que la topologie du flux thermique en parois est déterminée par la structure des écoulements secondaires alors que l’intensité du flux thermique aux parois est déterminée par le niveau de fluctuations de l’écoulement dans l’espace interperturbateur / Aeronautical engine designers are constantly subject to increasing power demands from aircraft manufacturers. To satisfy this requirement, combustor outlet temperature can be increased to improve efficiency and output energy of the engine. This rise in temperature however can surpass the material melting point and to avoid engine failure, turbine blades rely on internal cooling systems. Turbine blade cooling often uses internal channels, taking cold air from the compressor flow. Design of these systems therefore resumes to maximizing heat transfer enhancement while minimizing airflow rate to avoid engine power penalties. However, such flows are still largely uncontrolled and miss-understood. In an attempt to better understand such spatially developing rotating flows, the present study deals with a computational investigation on a straight, rotating rib roughened cooling channel. The configuration consists in a squared channel equipped with 8 ribs turbulators placed with an angle of 90 degrees with respect to the flow direction. For the studied cases, time resolved two-dimensional Particle Image Velocimetry (PIV) measurements have been performed at the Van Karman Institute (VKI). Adiabatic as well as isothermal conditions have been investigated to evaluate the impact of the wall temperature on the flow, especially in the rotating configurations. Static as well as both positive and negative rotating channels are compared with, in each case, either an adiabatic or an isothermal flow prediction. In this work, Large Eddy Simulation (LES) results show that the high fidelity CFD model is able to reproduce the differences induced by buoyancy on the flow topology in the near rib region and resulting from an adiabatic or an isothermal flow in rotation. The model manages also to predict the turbulence increase (decrease) around the rib in destabilizing (stabilizing) rotation of the ribbed channels. Finally and thanks to the full spatial and temporal description produced by LES, the spatial development and the unsteadiness of secondary flows are analyzed to better understand their origin and potential differences in all a cases. This study shows that the wall heat flux topology is driven by the secondary flows structure and the wall heat flux intensity is driven by the level of flow fluctuations in the ribbed region

<p> The progress of technology in the electronic components industry has been rapidly growing. The evolution of various techniques has made it possible for this industry to grow and diversify with the market demand. Thus, the development of electronic component products over a short span of time requires having highly efficient tools for design and manufacturing. Advances in commercial Computational Fluid Dynamics (CFD) softwares and computational power have enabled modeling to a high level of architectural details. Nowadays, computer aided design becomes an essential design tool in the engineering environment. Computer analysis reduces both the time development cycle and the prototyping costs in the early to intermediate design phases. The accuracy of the computational prediction of heat transfer rates depends mostly on the correct choice of turbulent model. Although many turbulent models, rather than a universal turbulent model, have been developed during the last two decades, there is usually one model that performs better than others for certain flow conditions. </p> <p> In the present research, a turbulence model is selected from amongst a few candidates, namely standard k- 8, RNG k- 8, shear stress transport (SST), and Reynolds Stress Model (RSM), based on comparisons with experimental data and direct numerical simulation (DNS) results from previous work. The SST turbulence model shows excellent agreement with the DNS results and, hence, is considered an appropriate turbulence model for thermal analysis of electronic packages with elements that have almost the same heights. Moreover, the average Nusselt number of array of obstacles is obtained numerically using commercial code ANSYS-CFX 1 0.0. The effects upon the mean Nusselt number arising from parameteric changes in Reynolds number, element height, element width, and element-to-element distance are compared and discussed. Finally, the parametric study has offered a set of correlations for the mean Nusselt number of arrays of mounted obstacles in the channel flow. </p> / Thesis / Master of Applied Science (MASc)

[Truncated abstract] Shipbuilding is fast becoming a priority industry in Australia. With increasing demands to build fast vessels of lighter weight, shipbuilders are more concerned with noise and vibration problems in ships than ever. The objective of this thesis is to study the vibration response of coupled structures, in the hope that the study may shed some light in understanding the general features of ship vibration. An important feature characterizing the vibration in complex structures is the input mobility, as it describes the capacity of structures in accepting vibration energy from sources. The input mobilities of finite ribbed plate and plate/plate coupled structures are investigated analytically and experimentally in this study. It is shown that the input mobility of a finite ribbed plate is bounded by the input mobilities of the uncoupled plate and beam(s) that form the ribbed plate and is dependent upon the distance between the source location and the stiffened beam(s). Off-neutral axis loading on the beam (point force applied on the beam but away from the beam’s neutral axis) affects the input power, kinetic energy distribution in the component plates of the ribbed plate and energy flow into the plates from the beam under direct excitation ... solutions were then used to examine the validity of statistical energy analysis (SEA) in the prediction of vibration response of an L-shaped plate due to deterministic force excitations. It was found that SEA can be utilized to predict the frequency averaged vibration response and energy flow of L-shaped plates under deterministic force (moment) excitations providing that the source location is more than a quarter of wavelength away from the plate edges. Furthermore, a simple experimental method was developed in this study to evaluate the frequency dependent stiffness and damping of rubber mounts by impact test. Finally, analytical methods developed in this study were applied in the prediction of vibration response of a ship structure. It was found that input mobilities of ship hull structures due to machinery excitations are governed by the stiffness of the supporting structure to which the engine is mounted. Their frequency averaged values can be estimated from those of the mounting structure of finite or infinite extents. It was also shown that wave propagation in ship hull structures at low frequencies could be attenuated by irregularities imposed to the periodic locations of the ship frames. The vibration at higher frequencies could be controlled by modifications of the supporting structure.

Modern gas turbine engines operate at temperatures well above the melting point of the metal components. This has driven manufactures to develop sophisticated cooling methods which minimise the use of coolant to maximise engine efficiency by enabling further increases in operating temperature. This thesis investigates the cooling performance of engine representative mid-chord internal cooling passages for turbine blades. The work forms part of a larger E.C. FP7 project ERICKA (Engine Representative Internal Cooling Knowledge Applications).This thesis provides detailed maps of heat transfer coefficient (HTC) under a number of conditions, new experimental techniques, and has lead to a better understanding of the impact HTC distributions have on the thermal performance of a turbine blade at engine conditions. Transient liquid crystal experiments have been conducted on a large scale model of an engine representative internal cooling passage at three aspect ratios (width:height (chord length:spanwise length), 1:2, 1:3 and 1:4). Spatially resolved maps of Nusselt number have been produced for the full surface of the internal cooling passages. Little information exists in the literature for more engine representative geometries, and it is rare for spatial measurements to be presented over the full surface. The detailed maps provide validation data for CFD within the ERICKA programme. A novel method which produces spatially resolved maps in areas with highly non-one-dimensional heat transfer has been developed and validated. This method couples transient finite element analysis and data from transient liquid crystal experiments. Applied to the ribbed passage geometry, this produced spatially resolved maps of HTC over the rib surface. To the author’s best knowledge this is the first time spatial HTC maps have been presented for an engine representative rib. Industry best practice methods for internal cooling passage design typically apply averaged values of HTC, in part due to lack of spatially resolved data. To determine the significance of this approximation on blade design and life, experimental measurements have been applied to finite element (FEA) models at typical engine conditions. Application of a 3D HTC distribution to a FEA model of a section of ribbed wall demonstrated a significant under prediction (up to 58%) of localised thermal gradients when an average value is applied compared to a spatially resolved profile. This work demonstrated good agreement between distributions taken from experimental data and CFD predictions, indicating that CFD distributions may be more appropriate than bulk values in the design process. A 2D FEA study was undertaken to quantify the impact of HTC distribution approximations and aspect ratio on cooling of a generic turbine section. This study considered multiple adjacent internal cooling passages. It was confirmed that multi-pass arrangements offer greater heat removal for a given mass flow rate. Also a symmetric heat transfer profile with a higher HTC on the ribbed wall is the most desirable distribution. Use of average values significantly impacted the metal temperature, causing an underprediction up to 13◦C and 8◦C in the maximum and average values respectively. Based on the experimental HTC data, the 1:3 aspect ratio passage offered the lowest metal temperatures. Applying HTC distributions from CFD data (calculated with using the centreline temperature) showed, in general, good agreement, with the lowest metal temperatures (by up to 8◦C) in the 1:4 aspect ratio passage. Use of and HTC distribution provided by CFD prediction based on the mixed bulk temperature, produced average and peak metal temperatures 16◦C and 17◦C, respectively, lower in the 1:4 aspect ratio passage than the next best design. This highlights the need for appropriate and consistent method to be used in the analysis. As expected, reducing the passage aspect ratio led to increases in both thermal gradient and total pressure loss.

Self-cleaning, anti-corrosion, anti-icing, dropwise-condensation, and drag-reduction are some applications in which superhydrophobic surfaces are implemented. To date, all the studies associated with superhydrophobic surfaces have been dedicated to understanding the liquid interaction with surfaces that are macroscopically smooth. The current study investigates the solid-liquid interaction of such surfaces which are fully decorated with macroscopic ridges (ribbed surfaces). In particular, the drop motion and impact on our newly designed non-wettable ribbed surface have been investigated in this work. Our experimental investigations have shown that liquid drops move faster on the ribbed surfaces due to lower friction induced by such a surface pattern. Moreover, an impacting droplet shows shorter contact time on ribbed surfaces. This concludes that ribbed surface pattern can be an efficient alternative design for the related applications. Besides the experimental studies, the theoretical analyses done in this work have led to, firstly a scaling model to predict descent velocity of a rolling viscous drops on an inclined non-wettable surface more accurately. Secondly, for curved superhydrophobic surfaces a scaling model which correlates the contact time of the impacting drop to its impact velocity has been developed. At the end, the knowledge obtained from this work has led to a special surface design which exhibits a contact time shorter than the inertial-capillary time scale, an unprecedented phenomenon.

<p>A first part of a ”Heat Transfer Handbook” about correlations for internal cooling of gas turbine vanes and blades has been created. The work is based on the cooling of vanes and blades 1 and 2 on different Siemens Gas Turbines. The cooling methods increase the heat transfer in the cooling channels by increasing the heat transfer coefficient and/or increasing the heat transfer surface area. The penalty paid for the increased heat transfer is higher pressure losses.</p><p>Three cooling methods, called rib turbulated cooling, matrix cooling and impingement cooling were investigated. Rib turbulated cooling and impingement cooling are typically used in the leading edge or mid region of the airfoil and matrix cooling is mostly applied in the trailing edge region.</p><p>Literature studies for each cooling method, covering both open literature and internal reports, were carried out in order to find correlations developed from tests. The correlations were compared and analyzed with focus on suitability for use in turbine conditions. The analysis resulted in recommendations about what correlations to use for each cooling method.</p><p>For rib turbulated cooling in square or rectangular ducts, four correlations developed by Han and his co-workers [3.5], [3.8], [3.9] and [3.6] are recommended, each valid for different channel and rib geometries. For U-shaped channels, correlations of Nagoga [3.4] are recommended.</p><p>Matrix cooling is relatively unknown in west, but has been used for many years in the former Soviet Union. Therefore available information in open literature is limited. Only one source of correlations was found. The correlations were developed by Nagoga [4.2] and are valid for closed matrixes. Siemens Gas Turbines are cooled with open matrixes, why further work with developing correlations is needed.</p><p>For impingement cooling on a flat target plate, a correlation of Florschuetz et al. [5.7] is recommended for inline impingement arrays. For staggered arrays, both the correlations of Florschuetz et al. [5.7] and Höglund [5.8] are suitable. The correlations for impingement on curved target plate gave very different results. The correlation of Nagoga is recommended, but it is also advised to consult the other correlations when calculating heat transfer for a specific case.</p><p>Another part of the work has been to investigate the codes of two heat transfer programs named Q3D and Multipass, used in the Siemens offices in Finspång and Lincoln, respectively. Certain changes in the code are recommended.</p>

A first part of a ”Heat Transfer Handbook” about correlations for internal cooling of gas turbine vanes and blades has been created. The work is based on the cooling of vanes and blades 1 and 2 on different Siemens Gas Turbines. The cooling methods increase the heat transfer in the cooling channels by increasing the heat transfer coefficient and/or increasing the heat transfer surface area. The penalty paid for the increased heat transfer is higher pressure losses. Three cooling methods, called rib turbulated cooling, matrix cooling and impingement cooling were investigated. Rib turbulated cooling and impingement cooling are typically used in the leading edge or mid region of the airfoil and matrix cooling is mostly applied in the trailing edge region. Literature studies for each cooling method, covering both open literature and internal reports, were carried out in order to find correlations developed from tests. The correlations were compared and analyzed with focus on suitability for use in turbine conditions. The analysis resulted in recommendations about what correlations to use for each cooling method. For rib turbulated cooling in square or rectangular ducts, four correlations developed by Han and his co-workers [3.5], [3.8], [3.9] and [3.6] are recommended, each valid for different channel and rib geometries. For U-shaped channels, correlations of Nagoga [3.4] are recommended. Matrix cooling is relatively unknown in west, but has been used for many years in the former Soviet Union. Therefore available information in open literature is limited. Only one source of correlations was found. The correlations were developed by Nagoga [4.2] and are valid for closed matrixes. Siemens Gas Turbines are cooled with open matrixes, why further work with developing correlations is needed. For impingement cooling on a flat target plate, a correlation of Florschuetz et al. [5.7] is recommended for inline impingement arrays. For staggered arrays, both the correlations of Florschuetz et al. [5.7] and Höglund [5.8] are suitable. The correlations for impingement on curved target plate gave very different results. The correlation of Nagoga is recommended, but it is also advised to consult the other correlations when calculating heat transfer for a specific case. Another part of the work has been to investigate the codes of two heat transfer programs named Q3D and Multipass, used in the Siemens offices in Finspång and Lincoln, respectively. Certain changes in the code are recommended.

Free liquid-jet impingement is well researched due to its high heat transfer ability and ease of implementation. This study considers both the steady state and transient heating of a patterned plate under slot-free-liquid jet impingement. The primary working fluid was water (H2O) and the plate material considered was silicon. Calculations were done for Reynolds number (Re) ranging from 500 to 1000 and indentation depths from 0.000125 to 0.0005 m for three different surface configurations. The effect of using different plate materials and R-134a as the working fluid were explored for the rectangular step case. The distributions of the local and average heat-transfer coefficient and the local and average Nusselt number were calculated for each case. A numerical model based in the FIDAP computer code was created to solve the conjugate heat transfer problem. The model used was developed for Cartesian coordinates for both steady state and transient conditions. Results show that the addition of surface geometry alters the fluid flow and heat transfer values. The highest heat-transfer coefficients occur at points where the fluid flow interacts with the surface geometry. The lowest heat-transfer coefficients are found in the indentations between the changes in geometry. The jet velocity has a large impact on the heat transfer values for all cases, with increasing jet velocity showing increased local heat-transfer coefficients and Nusselt number. It is observed that increasing the indentation depth for the rectangular and sinusoidal surfaces leads to a decrease in local heat transfer whereas for triangular patterns, a higher depth results in higher heat-transfer coefficient. The transient analysis showed that changing surface geometry had little effect on the time required to reach steady state. The selection of plate material has an impact on both the final maximum temperatures and the time required to reach steady state, with both traits being tied to the thermal diffusivity (α) of the material.

During the impinging heat transfer, a jet of working fluid, either gas or liquid, will besprayed onto the heat transfer surface. Due to the high turbulence of the fluid, the heat transfer coefficient between the wall and the fluid will be largely enhanced. Previously, an impinging type solar receiver with a cylindrical cavity absorber was designed for solar dish system. However, non-uniform temperature distribution in the circumferential direction was found on absorber surface from the numerical model, which will greatly limit receiver's working temperature and finally affect receiver's efficiency. One of the possible alternatives to solve the problem is through modifying the roughness of the target wall surface. This thesis work aims to evaluate the possibility and is focusing on the study of heat transfer characteristics. The simulation results will be used for future experimental impinging solar receiver optimization work. Computational Fluid Dynamics (CFD) is used to model the conjugate heat transfer phenomenon of atypical air impinging system. The simulation is divided into two parts. The first simulation was conducted with one rib arranged on the target surface where heat transfer coefficient is relatively low to demonstrate the effects of rib shape (triangular,rectangular, and semi-circular) and rib height (2.5mm, 1.5mm, and 0.5mm). The circular rib with 1.5mm height is proved to be most effective among all to acquirerelatively uniform temperature distribution. In the second part, the amount of ribs is taken into consideration in order to reach more uniform surface heat flux. The target wall thickness is also varied to assess its influence.