LES and hybrid LES/RANS simulations for conjugate heat transfer over a matrix of cubes

Large Eddy Simulation (LES) and a novel k -l based hybrid LES/RANS approach have been applied to simulate a conjugate heat transfer problem involving flow over a matrix of surface mounted cubes. In order to assess the capability and reliability of the newly developed k -l based hybrid LES/RANS, numerical results are compared with new LES and existing RANS results. Comparisons include mean velocity profiles, Reynolds stresses and conjugate heat transfer. As well as for hybrid LES/RANS validation purposes, the LES results are used to gain insights into the complex flow physics and heat transfer mechanisms. Numerical simulations show that the hybrid LES/RANS approach is effective. Mean and instantaneous fluid temperatures adjacent to the cube surface are found to strongly correlate with flow structure. Although the LES captures more mean velocity field complexities, broadly time averaged wake temperature fields are found similar for the LES and hybrid LES/RANS. Copyright © 2005 by the American Institute of Aeronautics and Astronautics, Inc. All rights reserved.

Conjugate heat transfer and film cooling of a multi-stage cooling scheme

Cooling techniques play a key role in improving efficiency and power output of modern gas turbines. The conjugate technique of film and impingement cooling schemes is considered in this study. The Multi-Stage Cooling Scheme (MSCS) involves coolant passing from inside to outside turbine blade through two stages. The first stage; the coolant passes through first hole to internal gap where the impinging jet cools the external layer of the blade. Finally, the coolant passes through the internal gap to the second hole which has specific designed geometry for external film cooling. The effect of design parameters, such as, offset distance between two-stage holes, gap height, and inclination angle of the first hole, on upstream conjugate heat transfer rate and downstream film cooling effectiveness performance are investigated computationally. An Inconel 617 alloy with variable properties is selected for the solid material. The conjugate heat transfer and film cooling characteristics of MSCS are analyzed across blowing ratios of Br = 1 and 2 for density ratio, 2. This study presents upstream wall temperature distributions due to conjugate heat transfer for different gap design parameters. The maximum film cooling effectiveness with upstream conjugate heat transfer is less than adiabatic film cooling effectiveness by 24–34%. However, the full coverage of cooling effectiveness in spanwise direction can be obtained using internal cooling with conjugate heat transfer, whereas adiabatic film cooling effectiveness has narrow distribution.

Conjugate Heat Transfer Analysis of an Internally Cooled Gas Turbine Vane

A conjugate heat transfer (CHT) method was used to perform the aerothermal analysis of an internally cooled turbine vane, and was validated against experimental and empirical data.
Firstly, validation of the method with regard to internal cooling was done by reproducing heat transfer test data in a channel with pin fin heat augmenters, under steady constant wall temperature. The computed Nusselt numbers for the two tested configurations (full length circular pin fins attached to both walls and partial pin fins attached to one wall only) showed good agreement with the measurements. Sensitivity to mesh density was evaluated under this simplified case in order to establish mesh requirements for the analysis of the full component.
Secondly, the CHT method was applied onto a turbine vane test case from an actual engine. The predicted vane airfoil metal temperature was compared to the measured thermal paint data and the in-house empirical predictions. The CHT results agreed well with the thermal paint data and showed better prediction than the current empirical modeling approach.

Implementation of the conjugate heat transfer code in KIVA-4

KIVA is an open Computational Fluid Dynamics (CFD) source code that is capable to compute the transient two and three-dimensional chemically reactive fluid flows with spray. The latest version in the family of KIVA codes is the KIVA-4 which is capable of handling the unstructured mesh. This project focuses on the implementation of the Conjugate Heat Transfer code (CHT) in KIVA-4. The previous version of KIVA code with conjugate heat transfer code has been developed at Michigan Technological University by Egel Urip and is be used in this project. During the first phase of the project, the difference in the code structure between the previous version of KIVA and the KIVA-4 has been studied, which is the most challenging part of the project. The second phase involves the reverse engineering where the CHT code in previous version is extracted and implemented in KIVA-4 according to the new code structure. The validation of the implemented code is performed using a 4-valve Pentroof engine case. A solid cylinder wall has been developed using GRIDGEN which surrounds 3/4th of the engine cylinder and heat transfer to the solid wall during one engine cycle (0-720 Crank Angle Degree) is compared with that of the reference result. The reference results are nothing but the same engine case run in the previous version with the original code developed by Egel. The results of current code are very much comparable to that of the reference results which verifies that successful implementation of the CHT code in KIVA-4.

Modeling study on the flow, heat transfer and energy conversion characteristics of low-power arc-heated hydrogen/nitrogen thrusters

A modeling study is conducted to investigate the effect of hydrogen content in propellants on the plasma flow, heat transfer and energy conversion characteristics of low-power (kW class) arc-heated hydrogen/nitrogen thrusters (arcjets). 1:0 (pure hydrogen), 3:1 (to simulate decomposed ammonia), 2:1 (to simulate decomposed hydrazine) and 0:1 (pure nitrogen) hydrogen/nitrogen mixtures are chosen as the propellants. Both the gas flow region inside the thruster nozzle and the anode-nozzle wall are included in the computational domain in order to better treat the conjugate heat transfer between the gas flow region and the solid wall region. The axial variations of the enthalpy flux, kinetic energy flux, directed kinetic-energy flux, and momentum flux, all normalized to the mass flow rate of the propellant, are used to investigate the energy conversion process inside the thruster nozzle. The modeling results show that the values of the arc voltage, the gas axial-velocity at the thruster exit, and the specific impulse of the arcjet thruster all increase with increasing hydrogen content in the propellant, but the gas temperature at the nitrogen thruster exit is significantly higher than that for other three propellants. The flow, heat transfer, and energy conversion processes taking place in the thruster nozzle have some common features for all the four propellants. The propellant is heated mainly in the near-cathode and constrictor region, accompanied with a rapid increase of the enthalpy flux, and after achieving its maximum value, the enthalpy flux decreases appreciably due to the conversion of gas internal energy into its kinetic energy in the divergent segment of the thruster nozzle. The kinetic energy flux, directed kinetic energy flux and momentum flux also increase at first due to the arc heating and the thermodynamic expansion, assume their maximum inside the nozzle and then decrease gradually as the propellant flows toward the thruster exit. It is found that a large energy loss (31-52%) occurs in the thruster nozzle due to the heat transfer to the nozzle wall and too long nozzle is not necessary. Modeling results for the NASA 1-kW class arcjet thruster with hydrogen or decomposed hydrazine as the propellant are found to compare favorably with available experimental data.

Single-phase fluid flow distribution and heat transfer in microstructured reactors

Single-phase microreactors and micro-heat-exchangers have been widely used in industrial and scientific applications over the last decade. In several cases, operation of microreactors has shown that their expected efficiency cannot be reached either due to non-uniform distribution of reactants between different channels or due to flow maldistribution between individual microreactors working in parallel. The latter problem can result in substantial temperature deviations between different microreactors resulting in thermal run away which could arise from an exothermicreaction. Thus advances in the understanding of heat transfer and fluid flow distribution continue to be crucial in achieving improved performance, efficiency and safety in microstructured reactors used for different applications. This paper presents a review of the experimental and numerical results on fluid flow distribution, heat transfer and combination thereof, available in the open literature. Heat transfer in microchannels can be suitably described by standard theory and correlations, but scaling effects (entrance effects, conjugate heat transfer, viscous heating, and temperature-dependent properties) have often to be accounted for in microsystems. Experiments with single channels are in good agreement with predictions from the published correlations. The accuracy of multichannel experiments is lower due to flow maldistribution. Special attention is devoted to theoretical and experimental studies on the effect of a flow maldistribution on the thermal and conversion response of catalytic microreactors. There view concludes with a set of design recommendations aimed at improving the reactor performance. (C) 2010 Elsevier Ltd. All rights reserved.

Natural convection and heat transfer in attics subject to periodic thermal forcing

The heat transfer through the attics of buildings under realistic thermal forcing has been considered in this study. A periodic temperature boundary condition is applied on the sloping walls of the attic to show the basic flow features in the attic space over diurnal cycles. The numerical results reveal that, during the daytime heating stage, the flow in the attic space is stratified; whereas at the night-time cooling stage, the flow becomes unstable. A symmetrical solution is seen for relatively low Rayleigh numbers. However, as the Ra gradually increases, a transition occurs at a critical value of Ra. Above this critical value, an asymmetrical solution exhibiting a pitchfork bifurcation arises at the night-time. It is also found that the calculated heat transfer rate at the night-time cooling stage is much higher than that during the daytime heating stage.

Scaling of free convection heat transfer in a triangular cavity for Pr>1

Unsteady natural convection inside a triangular cavity subject to a non-instantaneous heating on the inclined walls in the form of an imposed temperature which increases linearly up to a prescribed steady value over a prescribed time is reported. The development of the flow from start-up to a steady-state has been described based on scaling analyses and direct numerical simulations. The ramp temperature has been chosen in such a way that the boundary layer is reached a quasi-steady mode before the growth of the temperature is completed. In this mode the thermal boundary layer at first grows in thickness, then contracts with increasing time. However, if the imposed wall temperature growth period is sufficiently short, the boundary layer develops differently. It is seen that the shape of many houses are isosceles triangular cross-section. The heat transfer process through the roof of the attic-shaped space should be well understood. Because, in the building energy, one of the most important objectives for design and construction of houses is to provide thermal comfort for occupants. Moreover, in the present energy-conscious society it is also a requirement for houses to be energy efficient, i.e. the energy consumption for heating or air-conditioning houses must be minimized.

A review of natural convection and heat transfer in attic-shaped space

This paper presents a review of studies on natural convection heat transfer in the triangular enclosure namely, in attic-shaped space. Much research activity has been devoted to this topic over the last three decades with a view to providing thermal comfort to the occupants in attic-shaped buildings and to minimising the energy costs associated with heating and air-conditioning. Two basic thermal boundary conditions of attic are considered to represent hot and cold climates or day and night time. This paper also reports on a significant number of studies which have been performed recently on other topics related to the attic space, for example, attics subject to localized heating and attics filled with porous media.

Free convection heat transfer through coupled thermal boundary layers in a partitioned triangular enclosure

A numerical investigation has been carried out for the coupled thermal boundary layers on both sides of a partition placed in an isosceles triangular enclosure along its middle symmetric line. The working fluid is considered as air which is initially quiescent. A sudden temperature difference between two zones of the enclosure has been imposed to trigger the natural convection. It is anticipated from the numerical simulations that the coupled thermal boundary layers development adjacent to the partition undergoes three distinct stages; namely an initial stage, a transitional stage and a steady state stage. Time dependent features of the coupled thermal boundary layers as well as the overall natural convection flow in the partitioned enclosure have been discussed and compared with the non-partitioned enclosure. Moreover, heat transfer as a form of local and overall average Nusselt number through the coupled thermal boundary layers and the inclined walls is also examined. The details results will be discussed in the full paper.

Analysis of heat transfer and fluid flow in a triangular enclosure in presence of an adiabatic fin

Natural convection thermal boundary layer adjacent to the heated inclined wall of a right angled triangle with an adiabatic fin attached to that surface is investigated by numerical simulations. The finite volume based unsteady numerical model is adopted for the simulation. It is revealed from the numerical results that the development of the boundary layer along the inclined surface is characterized by three distinct stages, i.e. a start-up stage, a transitional stage and a steady stage. These three stages can be clearly identified from the numerical simulations. Moreover, in presence of adiabatic fin, the thermal boundary layer adjacent to the inclined wall breaks initially. However, it is reattached with the downstream boundary layer next to the fin. More attention has been given to the boundary layer development near the fin area.

Unsteady buoyancy driven flows and heat transfer through coupled thermal boundary layers in a partitioned triangular enclosure

A numerical investigation has been carried out for the coupled thermal boundary layers on both sides of a partition placed in an isosceles triangular enclosure along its middle symmetric line. The working fluid is considered as air which is initially quiescent. A sudden temperature difference between two zones of the enclosure has been imposed to trigger the natural convection. It is anticipated from the numerical simulations that the coupled thermal boundary layers development adjacent to the partition undergoes three distinct stages; namely an initial stage, a transitional stage and a steady state stage. Time dependent features of the coupled thermal boundary layers as well as the overall natural convection flow in the partitioned enclosure have been discussed and compared with the non-partitioned enclosure. Moreover, heat transfer as a form of local and overall average Nusselt number through the coupled thermal boundary layers and the inclined walls is also examined.