This study performed detailed measurements on and around scaled up conducting and adiabatic airfoils with and without film cooling. The conducting vane was a matched Bi airfoil, which accurately scaled the convective heat transfer and conduction through the solid, in order to produce non-dimensional surface temperatures and thermal boundary layers that were representative of an actual engine. Measurements made on all vane models included surface temperature measurements and thermal profiles above the walls. Separate measurements on non-film cooled and film cooled conducting models allowed for the individual contributions of the internal convective cooling and external film cooling to the overall cooling scheme to be quantified. Surface temperature and thermal field measurements above the wall were also performed on a film cooled adiabatic model. For the conducting model with internal cooling only, strong streamwise temperature variations were seen. The surface temperature variations were highly dependent on the local external and internal heat transfer coefficients. Spanwise temperature variations also existed, but were modest in comparison to streamwise variations. Comparing the thermal fields above the film cooled adiabatic and conducting walls allowed for the assumption that the conducting wall would not significantly affect the thermal field in the film cooling jet to be tested. Near the edge of the film cooling jet the developing thermal boundary layer had a clear effect on the overlying gas temperature, suggesting that the common assumption that the adiabatic wall temperature is the appropriate driving temperature for heat transfer to a film cooled wall was invalid. On the jet centerline thermal boundary layer effects were less influential, due to the development of a new, thin boundary layer. This suggested that the adiabatic wall temperature as driving temperature for heat transfer was a reasonable assumption on the jet centerline for most cases tested. As film cooling momentum flux ratio increase, thermal boundary layer effects became more influential on the jet centerline. Additionally, the high resolution surface temperature measurements and thermal field measurements above the wall presented in the current study represent a significant improvement in the data available for validation of computational simulations of conducting turbine airfoils. / text

Dans un souci d’optimisation des fours industriels et de réduction des émissions de gaz à effet de serre,l’oxy-combustion est considérée comme l’une des solutions d’avenir. Les conditions existantes dans les chambres d’oxycombustion créent une interaction forte entre les différents phénomènes : Combustion,turbulence et transferts de chaleur. Pour mieux dimensionner les configurations futures il est nécessaire de pouvoir étudier la physique qui y règne, et ce pour un coût et un temps de retour raisonnables. De tels études nécessitent l’emploi d’outils de simulation de haute fidélité,et afin de modéliser les interactions inter-phénomènes à un coût acceptable le couplage de codes est utilisé. C’est avec cet objectif que les travaux présentés dans ce manuscrit se concentrent sur la mise au point d’une méthodologie de couplage entre codes d’écoulements réactifs et de transfert de chaleur dans les parois pour la réalisation de simulations de haute-fidélité massivement parallèles prédictives des chambres futures. / Oxycombustion is seen as one mean to attain the wished goals in terms of efficiency optimisation and Greenhouse Effect Gases emissions reduction for industrial furnaces. The extreme operating conditions, high pressure and temperature, lead to a strong interaction between the different phenomena which take place inside the combustion chambe r: Combustion, turbulence and heat transfer. To better design these futur oxyfuel processes, a mean to study the related physics with a reasonable computational cost and return time. Such studies require the use of high-fidelity numerical resolution tools, and in order to model the multi-physics interaction in a cost efficient way, code coupling. The operating conditions being extreme : High pressure and temperature, a strong interaction exists between the different phenomena occuring inside the chamber. To better understand the physics inside oxycombustion chambers,a multiphysics high-fidelity simulation methodology is developped.

Les simulations sont utilisées pour concevoir des chambres de combustion industrielles robustes et peu polluantes. Parmi les polluants, l’émission de particules de suies constitue une question sociétale et une priorité politico- industrielle, en raison de leurs impacts néfastes sur la santé et l'environnement. La taille des particules de suies joue un rôle important sur ces effets. Il est donc important de prévoir non seulement la masse totale ou le nombre de particules générées, mais également leur distribution en taille (PSD). De plus, les suies peuvent jouer un rôle important dans le rayonnement thermique. Dans des configurations confinées, la prédiction des transferts de chaleur est une question clé pour augmenter la robustesse des chambres de combustion. Afin de déterminer correctement ces transferts, les flux radiatifs et de conducto-convectifs aux parois doivent être pris en compte. Enfin, la température pariétale est aussi contrôlée par les transferts conjugués de chaleur entre les domaines fluides et solides. L’ensemble de ces transferts thermiques impactent la stabilisation de la flamme, la formation de polluants et la production de suies elle-même. Il existe donc un couplage complexe entre ces phénomènes et la simulation d'un tel problème multiphysique est aujourd'hui reconnu comme un important défi. Ainsi, l'objectif de cette thèse est de développer une modélisation multiphysique permettant la simulation de flammes suitées turbulentes avec le rayonnement thermique et les transferts conjugués de chaleur associés aux parois. Les méthodes retenues sont basées sur la Simulation aux Grandes Échelles (LES), une description en taille des suies, des transferts conjugués et un code Monte Carlo pour le rayonnement. La combinaison de telles approches est réalisable grâce aux ressources de calcul aujourd’hui disponibles afin d’obtenir des résultats de référence. Le manuscrit est organisé en trois parties. La première partie se concentre sur le développement d'un modèle détaillé pour la description de la production de suies dans les flammes laminaires. Pour cela, la méthode sectionnelle est retenue ici car elle permet la description de la PSD. La méthode est validée sur des flammes laminaires éthylène/air. Dans la deuxième partie, un formalisme LES spécifique à la méthode sectionnelle est développé et utilisé pour étudier deux flammes turbulentes : une flamme jet non-prémélangée et une flamme swirlée pressurisée confinée. Les champs de température et de fraction volumique de suies sont comparés aux données expérimentales. De bonnes prédictions sont obtenues et l’évolution des particules de suies dans de telles flammes est analysée à travers l'étude de l’évolution de leur PSD. Dans ces premières simulations, les pertes de chaleur aux parois reposent sur des mesures expérimentales de la température aux parois, et un modèle de rayonnement simple. Dans la troisième partie, une approche Monte Carlo permettant de résoudre l'équation de transfert radiatif avec des propriétés radiatives détaillées des phases gazeuse et solide est utilisée et couplée au solveur LES. Cette approche est appliquée à l'étude de la flamme jet turbulente. La prédiction des flux thermiques est comparée aux données expérimentales et la nature des transferts radiatifs est étudiée. Ensuite, une modélisation couplée de la combustion turbulente prenant en compte la production de suies, les transferts conjugués de chaleur et le rayonnement thermique est proposée en couplant les trois codes dédiés. Cette stratégie est appliquée pour la simulation du brûleur pressurisé confiné. L'approche proposée permet à la fois de prédire la température des parois et la bonne stabilisation de la flamme. Les processus de formation de suies se révèlent être affectés par la modélisation des transferts thermiques. Ceci souligne l’importance d’une description précise de ces transferts thermiques dans les développements futurs de modèles de production de suies et leur validation. / Numerical simulations are used by engineers to design robust and clean industrial combustors. Among pollutants, soot control is an urgent societal issue and a political-industrial priority, due to its harmful impact on health and environment. Soot particles size plays an important role in its negative effect. It is therefore important to predict not only the total mass or number of emitted particles, but also their population distribution as a function of their size. In addition, soot particles can play an important role in thermal radiation. In confined configurations, controlling heat transfer related to combustion is a key issue to increase the robustness and the life cycle of combustors by avoiding wall damages. In order to correctly determine these heat losses, radiative and wall convective heat fluxes must be accounted for. They depend on the wall temperature, which is controlled by the conjugate heat transfer between the fluid and solid domains. Heat transfer impacts the flame stabilization, pollutants formation and soot production itself. Therefore, a complex coupling exists between these phenomena and the simulation of such a multi-physics problem is today recognized as an extreme challenge in combustion, especially in a turbulent flow, which is the case of most industrial combustors. Thus, the objective of this thesis is to develop a multi-physics modeling enabling the simulation of turbulent sooting flames including thermal radiation and wall heat transfer. The retained methods based on Large-Eddy Simulation (LES), a soot sectional model, conjugate heat transfer, a Monte Carlo radiation solver are combined to achieve a stateof- the-art framework. The available computational resources make nowadays affordable such simulations that will yield present-day reference results. The manuscript is organized in three parts. The first part focuses on the definition of a detailed model for the description of soot production in laminar flames. For this, the sectional method is retained here since it allows the description of the particle size distribution (PSD). The method is validated on laminar premixed and diffusion ethylene/air flames before analyzing the dynamics of pulsed diffusion flames. In the second part, an LES formalism for the sectional method is developed and used to investigate two different turbulent flames: a non-premixed jet flame and a confined pressurized swirled flame. Predicted temperature and soot volume fraction levels and topologies are compared to experimental data. Good predictions are obtained and the different soot processes in such flames are analyzed through the study of the PSD evolution. In these first simulations, wall heat losses rely on experimental measurements of walls temperature, and a coarse optically-thin radiation model. In the third part, to increase the accuracy of thermal radiation description, a Monte Carlo approach enabling to solve the Radiative Transfer Equation with detailed radiative properties of gaseous and soot phases is used and coupled to the LES solver. This coupled approach is applied for the simulation of the turbulent jet flame. Quality of radiative fluxes prediction in this flame is quantified and the nature of radiative transfers is studied. Then, a whole coupled modeling of turbulent combustion accounting for soot, conjugate heat transfer and thermal radiation is proposed by coupling three dedicated codes. This strategy is applied for a high-fidelity simulation of the confined pressurized burner. By comparing numerical results with experimental data, the proposed approach enables to predict both the wall temperature and the flame stabilization. The different simulations show that soot formation processes are impacted by the heat transfer description: a decrease of the soot volume fraction is observed with increasing heat losses. This highlights the requirement of accurate description of heat transfer for future developments of soot models and their validation.

An average heat transfer coefficient, h_bar, is often used to solve heat transfer problems. It should be understood that this is an approximation and may provide inaccurate results, especially when the temperature field is of interest. The proper method to solve heat transfer problems is with a conjugate approach. However, there seems to be a lack of clear explanations of conjugate heat transfer in literature. The objective of this work is to provide a clear explanation of conjugate heat transfer and to determine the discrepancy in the temperature field when the interface boundary condition is approximated using h_bar compared to a local, or variable, heat transfer coefficient, h(x). Simple one-dimensional problems are presented and solved analytically using both h(x) and h_bar. Due to the one-dimensional assumption, h(x) appears in the governing equation for which the common methods to solve the differential equations with an average coefficient are no longer valid. Two methods, the integral equation and generalized Bessel methods are presented to handle the variable coefficient. The generalized Bessel method has previously only been used with homogeneous governing equations. This work extends the use of the generalized Bessel method to non-homogeneous problems by developing a relation for the Wronskian of the general solution to the generalized Bessel equation. The solution methods are applied to three problems: an external flow past a flat plate, a conjugate interface between two solids and a conjugate interface between a fluid and a solid. The main parameter that is varied is a combination of the Biot number and a geometric aspect ratio, A_1^2 = Bi*L^2/d_1^2. The Biot number is assumed small since the problems are one-dimensional and thus variation in A_1^2 is mostly due to a change in the aspect ratio. A large A_1^2 represents a long and thin solid whereas a small A_1^2 represents a short and thick solid. It is found that a larger A_1^2 leads to less problem conjugation. This means that use of h_bar has a lesser effect on the temperature field for a long and thin solid. Also, use of ¯ over h(x) tends to generally under predict the solid temperature. In addition is was found that A_2^2, the A^2 value for the second subdomain, tends to have more effect on the shape of the temperature profile of solid 1 and A_1^2 has a greater effect on the magnitude of the difference in temperature profiles between the use of h(x) and h_bar. In general increasing the A^2 values reduced conjugation.

In this thesis study, a least-squares finite element method (LSFEM) based conjugate heat transfer solver was developed. In the mentioned solver, fluid flow and heat transfer computations were performed separately. This means that the calculated velocity values in the flow calculation part were exported to the heat transfer part to be used in the convective part of the energy equation. Incompressible Navier-Stokes equations were used in the flow simulations. In conjugate heat transfer computations, it is required to calculate the heat transfer in both flow field and solid region. In this study, conjugate behavior was accomplished in a fully coupled manner, i.e., energy equation for fluid and solid regions was solved simultaneously and no boundary conditions were defined on the fluid-solid interface. To assure that the developed solver works properly, lid driven cavity flow, backward facing step flow and thermally driven cavity flow problems were simulated in three dimensions and the findings compared well with the available data from the literature. Couette flow and thermally driven cavity flow with conjugate heat transfer in two dimensions were modeled to further validate the solver. Finally, a microchannel conjugate heat transfer problem was simulated. In the flow solution part of the microchannel problem, conservation of mass was not achieved. This problem was expected since the LSFEM has problems related to mass conservation especially in high aspect ratio channels. In order to overcome the mentioned problem, weight of continuity equation was increased by multiplying it with a constant. Weighting worked for the microchannel problem and the mass conservation issue was resolved. Obtained results for microchannel heat transfer problem were in good agreement in general with the previous experimental and numerical works. In the first computations with the solver / quadrilateral and triangular elements for two dimensional problems, hexagonal and tetrahedron elements for three dimensional problems were tried. However, since only the quadrilateral and hexagonal elements gave satisfactory results, they were used in all the above mentioned simulations.

The present study considers the numerical simulation of the different flow characteristics involved in the conjugate heat transfer analysis of an internally cooled gas turbine blade. Conjugate simulations require full coupling of convective heat transfer in fluid regions to the heat diffusion in solid regions. Therefore, accurate prediction of heat transfer quantities on both external and internal surfaces has the uppermost importance and highly connected with the performance of the employed turbulence models. The complex flow on both surfaces of the internally cooled turbine blades is caused from the boundary layer laminar-to-turbulence transition, shock wave interaction with boundary layer, high streamline curvature and sequential flow separation. In order to discover the performances of different turbulence models on these flow types, analyses have been conducted on five different experimental studies each concerned with different flow and heat transfer characteristics. Each experimental study has been examined with four different turbulence models available in the commercial software (ANSYS FLUENT13.0) to decide most suitable RANS-based turbulence model. The Realizable k-&epsilon / model, Shear Stress Transport k-&omega / model, Reynolds Stress Model and V2-f model, which became increasingly popular during the last few years, have been used at the numerical simulations. According to conducted analyses, despite a few unreasonable predictions, in the majority of the numerical simulations, V2-f model outperforms other first-order turbulence models (Realizable k-&epsilon / and Shear Stress Transport k-&omega / ) in terms of accuracy and Reynolds Stress Model in terms of convergence.

This master thesis was provided by Scania AB. The objective of this thesis was to modify an application in the free Computational Fluid Dynamics software OpenFOAM to be able to handle spray and wall film modeling of a Urea Water Solution together with Conjugate Heat Transfer. The basic purpose is to widen the knowledge of the vaporization process of a Urea Water Solution in the exhaust gas after treatment system for a diesel engine by using OpenFOAM. First, urea has been modeled as a very viscous liquid at low temperature to mimic the solidication process of urea. Second, the development of the new application has been done. At last, test simulations of a simple test case are performed with the new application. The results are then compared with simplied hand calculations to verify a correct behavior of certain exposed source terms. The new application is working properly for the test case but to ensure the reliability, the results need to be compared with another Computational Fluid Dynamics software or more preferable, real experiments. For more advanced geometries, the continued development presented last in this thesis is highly recommended to follow.

This thesis presents a detailed, two-dimensional analysis of the performance of multi-element guard-heated hot-film wall shear stress microsensors for turbulent flows. Previous studies of conventional, single-element sensors show that a significant portion of heat generated in the hot-film travels through the substrate before reaching the fluid, causing spectral and phase errors in the wall shear stress signal and drastically reducing the spatial resolution of the sensor. Earlier attempts to reduce these errors have focused on reducing the effective thermal conductivity of the substrate. New guard-heated microsensor designs proposed to overcome the severe deficiencies of the conventional design are investigated in this thesis. Guard-heaters remove the errors associated with substrate heat conduction, by forcing zero temperature gradient at the edges and bottom face of the hot-film, and hence, block the indirect heat transfer to the flow. Air and water flow over the sensors are studied numerically to investigate design, performance and signal strength of the guard-heated sensors. Our results show, particularly for measurements in low-conductivity fluids such as air, that edge guard-heating needs to be supplemented by a sub-surface guard-heater, to make substrate conduction errors negligible. With this two-plane guard-heating, a strong non-linearity in the standard single-element designs can be corrected, and spectral and phase errors arising from substrate conduction can be eliminated. / Graduate / 0548 / etrati@uvic.ca

Design and analysis of conceptually different cooling systems for the human heart preservation are numerically investigated. A heart cooling container with required connections was designed for a normal size human heart. A three-dimensional, high resolution human heart geometric model obtained from CT-angio data was used for simulations. Nine different cooling designs are introduced in this research. The first cooling design (Case 1) used a cooling gelatin only outside of the heart. In the second cooling design (Case 2), the internal parts of the heart were cooled via pumping a cooling liquid inside both the heart’s pulmonary and systemic circulation systems. An unsteady conjugate heat transfer analysis is performed to simulate the temperature field variations within the heart during the cooling process. Case 3 simulated the currently used cooling method in which the coolant is stagnant. Case 4 was a combination of Case 1 and Case 2. A linear thermoelasticity analysis was performed to assess the stresses applied on the heart during the cooling process. In Cases 5 through 9, the coolant solution was used for both internal and external cooling. For external circulation in Case 5 and Case 6, two inlets and two outlets were designed on the walls of the cooling container. Case 5 used laminar flows for coolant circulations inside and outside of the heart. Effects of turbulent flow on cooling of the heart were studied in Case 6. In Case 7, an additional inlet was designed on the cooling container wall to create a jet impinging the hot region of the heart’s wall. Unsteady periodic inlet velocities were applied in Case 8 and Case 9. The average temperature of the heart in Case 5 was +5.0oC after 1500 s of cooling. Multi-objective constrained optimization was performed for Case 5. Inlet velocities for two internal and one external coolant circulations were the three design variables for optimization. Minimizing the average temperature of the heart, wall shear stress and total volumetric flow rates were the three objectives. The only constraint was to keep von Mises stress below the ultimate tensile stress of the heart’s tissue.

La problématique associée à ce travail de thèse concerne la détermination des causes de fatigue thermique dans des conduites hydrauliques. Ce travail porte plus particulièrement sur l’étude d’une zone de mélange orthogonale dans laquelle débouchent deux écoulements à différentes températures et présentant des régimes dynamiques turbulents. Ce type de configuration est utilisé afin de favoriser la baisse de température de l’écoulement (ou l’augmentation de température selon les cas). Pour des valeurs élevées du nombre de Reynolds et des écarts de températures des écoulements en confluence, la charge thermique engendrée par les fluctuations pariétales de température peut entraîner la formation de fissures dans les canalisations.Pour faire face à ce problème, les méthodes de dimensionnement des conduites hydrauliques mettent en évidence la nécessité de connaître la charge thermique imposée sur la paroi sensible. Le besoin attendu en termes de caractérisation dynamique et thermique de ces écoulements est à l’origine de ce projet de recherche qui s’attache ainsi à déterminer pour une configuration industrielle simplifiée, les phénomènes dynamiques et thermiques instationnaires et tridimensionnels à proximité de la paroi. Pour cela, l’étude s’organise en deux phases. Une approche expérimentale, faisant appel à la vélocimétrie par images de particules et la thermographie infrarouge, est mise en oeuvre dans un écoulement turbulent à la jonction entre une conduite de section rectangulaire et un piquage cylindrique. L’étude expérimentale permet de mesurer l’effet de différents paramètres de l’écoulement (régime dynamique, rapport des vitesses de l’écoulement, écart de température). Une des configurations expérimentales est de plus modélisée par simulations des grandes échelles dans laquelle le transfert thermique conjugué entre le fluide et la paroi d’intérêt est pris en compte. Les champs statistiques de vitesse et température issus du calcul numérique sont comparés à la base de données expérimentale. L’analyse des mécanismes thermodynamiques près de la paroi et dans la paroi, est de plus réalisée au moyen des grandeurs statistiques de l’écoulement, des corrélations de vitesse-température et de l’enregistrement temporel de sondes numériques placées dans l’écoulement et dans la paroi. / This research work originates from thermal stress issues in hydraulic pipes. It specifically focuses on determination of thermal load created by mixing flows in orthogonal pipe junctions. Such mixing zones are commonly used to cool or heat the main flow. For high Reynolds number values and high difference of temperature of the two flows, thermal load created by strong temperature fluctuations at the wall can initiate and propagate cracks through the pipe. Standards in pipe design pointed out the need in knowing these thermal loads applied at the downstream junction wall. Thus, this study aims to characterize three-dimensional and unsteady dynamical and thermal phenomena in a simplified industrial configuration. First, an experimental investigation based on particle image velocimetry and infrared thermography has been carried out on a cylindrical jet pipe flowing through a main rectangular channel flow. Different parameters have been studied such as Reynolds number influence, velocity ratio between secondary and main flow and temperature difference. One of the experimental configurations has also been investigated using numerical simulations. LargeEddy Simulations with conjugate heat transfer model has been developed and compared to the experimental database. Dynamics and Thermal mechanisms have been analysed near the wall and within the wall with obtained statistical flow fields, velocity-temperature moments and numerical probes extracted from fluid and solid parts.