Summary
Here, we present an experimental method for decoupling the interdependent Coriolisforce and rotatingbuoyancy effects on fullfield heat transfer distributions of a rotating channel.
Abstract
An experimental method for exploring the heat transfer characteristics of an axially rotating channel is proposed. The governing flow parameters that characterize the transport phenomena in a rotating channel are identified via the parametric analysis of the momentum and energy equations referring to a rotating frame of reference. Based on these dimensionless flow equations, an experimental strategy that links the design of the test module, the experimental program and the data analysis is formulated with the attempt to reveal the isolated Coriolisforce and buoyancy effects on heat transfer performances. The effects of Coriolis force and rotating buoyancy are illustrated using the selective results measured from rotating channels with various geometries. While the Coriolisforce and rotatingbuoyancy impacts share several common features among the various rotating channels, the unique heat transfer signatures are found in association with the flow direction, the channel shape and the arrangement of heat transfer enhancement devices. Regardless of the flow configurations of the rotating channels, the presented experimental method enables the development of physically consistent heat transfer correlations that permit the evaluation of isolated and interdependent Coriolisforce and rotatingbuoyancy effects on the heat transfer properties of rotating channels.
Introduction
While thermodynamic laws dictate the improved specific power and thermal efficiency of a gas turbine engine by elevating the turbine entry temperature, several hot engine components, such as turbine blades, are prone to thermal damage. Internal cooling of a gas turbine rotor blade permits a turbine entry temperature in excess of the temperature limits of the creep resistance of the blade material. However, the configurations of the internal cooling channels must comply with the blade profile. In particular, the coolant rotates within the rotor blade. With such harsh thermal conditions for a running gas turbine rotor blade, an effective blade cooling scheme is crucial to ensure the structure's integrity. Thus, the local heat transfer properties for a rotating channel are important for the efficient usage of the limited coolant flow available. The acquisition of useful heat transfer data that are applicable to the design of the internal coolant passages at realistic engine conditions is of primary importance when an experimental method is developed for measuring the heat transfer properties of a simulated cooling passage inside a gas turbine rotor blade.
Rotation at a speed above 10,000 rpm considerably alters the cooling performance of a rotating channel inside a gas turbine rotor blade. The identification of engine conditions for such a rotating channel is permissible using the similarity law. With rotation, the dimensionless groups that control the transport phenomena inside a radially rotating channel can be revealed by deriving the flow equations relative to a rotating frame of reference. Morris^{1} has derived the momentum conservation equation of flow relative to a rotating frame of reference as:
(1)
In equation (1), the local fluid velocity, v̄, with the position vector, r̄, relative to a frame of reference rotating at the angular velocity, ω, is affected by the Coriolis acceleration in terms of 2(ω×v̄), the decoupled centripetal buoyancy force, β(TT_{ref})(ω×ω×r̄), the driven piezometric pressure gradient, , and the fluid dynamic viscosity, ν. The referenced fluid density, ρ_{ref}, is referred to a predefined fluid reference temperature T_{ref}, which is typical of the local fluid bulk temperature for experiments. If the irreversible conversion of mechanical energy into thermal energy is negligible, the energy conservation equation is reduced to:
(2)
The first term of equation (2) is obtained by treating the specific enthalpy to be directly related to the local fluid temperature, T, via the constant specific heat, C_{p}. As the perturbation of fluid density caused by the variation of fluid temperature in a heated rotating channel provides considerable influence on the motion of fluids when it links with the centripetal acceleration in equation (1), the fluid velocity and temperature fields in an axially rotating channel are coupled. Also, both Coriolis and centripetal accelerations vary simultaneously as the rotating speed is adjusted. Thus, the effects of Coriolis force and rotating buoyancy on the fields of fluid velocity and temperature are naturally coupled.
Equations (1) and (2) in the dimensionless forms disclose the flow parameters that govern the heat convection in a rotating channel. With a basically uniform heat flux imposed on a rotating channel, the local fluid bulk temperature, T_{b}, increases linearly in the streamwise direction, s, from the reference inlet level, T_{ref}. The local fluid bulk temperature is determined as T_{ref} + τs, where τ is the gradient of the fluid bulk temperature in the direction of flow. Substitutions of the following dimensionless parameters of:
(3)
(4)
(5)
(6)
(7)
into equations (1) and (2), where V_{mean}, N and d respectively stand for the mean flow through velocity, rotating velocity and channel hydraulic diameter, the dimensionless flow momentum and energy equations are derived as equations (8) and (9) respectively.
(8)
(9)
Evidently, η in equation (9) is a function of Re, Ro, and Bu = Ro^{2}βτdR, which are respectively referred to as Reynolds, rotation and buoyancy numbers. The Rossby number that quantifies the ratio between inertial and Coriolis forces is equivalent to the inverse rotation number in equation (8).
When T_{b} is calculated as T_{ref} + τs in a rotating channel subject to a uniform heat flux, the τ value can be alternatively evaluated as Q_{f}/(mC_{p}L) in which Q_{f}, m and L are the convective heating power, coolant mass flow rate and channel length, respectively. Thus, the dimensionless local fluid bulk temperature, η_{b}, is equal to s/d and the dimensionless temperature at channel wall, η_{w}, yields [(T_{w}T_{b})/Q_{f}][mC_{p}][L/d]+s/d. With the convective heat transfer rate defined as Q_{f}/(T_{w}T_{b}), the dimensionless walltofluid temperature difference, η_{w}η_{b}, is convertible into the local Nusselt number via equation (10) in which ζ is the dimensionless shape function of heating area and channel sectional area.
(10)
With a set of predefined geometries and the hydrodynamic and thermal boundary conditions, the dimensionless groups controlling the local Nusselt number of a rotating channel are identified as:
(11)
(12)
(13)
With experimental tests, the adjustment of rotating speed, N, for varying Ro to generate the heat transfer data at different strengths of Coriolis forces inevitably changes the centripetal acceleration, and thus, the relative strength of rotating buoyancy. Moreover, a set of heat transfer data collected from a rotating channel is always subject to a finite degree of rotating buoyancy effect. To disclose the individual effects of Coriolisforce and buoyancy on the heat transfer performance of a rotating channel requires the uncoupling of the Ro and Bu effects on Nu properties through the post data processing procedure that is inclusive in the present experimental method.
The engine and laboratory flow conditions for a rotating channel inside a gas turbine rotor blade can be specified by the ranges of Re, Ro and Bu. The typical engine conditions for the coolant flow through a gas turbine rotor blade, as well as the construction and commissioning of the rotating test facility that allowed experiments to be performed near the actual engine conditions was reported by Morris^{2}. Based on the realistic engine conditions summarized by Morris^{2}, Figure 1 constructs the realistic operating conditions in terms of Re, Ro and Bu ranges for a rotating coolant channel in a gas turbine rotor blade. In Figure 1, the indication of an engine's worst condition is referred to as the engine running condition at the highest rotor speed and the highest density ratio. In Figure 1, the lower limit and worst engine operating conditions respectively emerge at the lowest and highest engine speeds. It is extremely difficult to measure the fullfield Nu distribution of a rotating channel running at a real engine speed between 5000 and 20,000 rpm. However, based on the similarity law, laboratoryscale tests have been conducted at reduced rotating speeds but with several attempts to provide a full coverage of the realengine Re, Ro and Bu ranges. As an innovative experimental method, the NASA HOST program^{3}^{,}^{4}^{,}^{5}^{,}^{6} adopted the highpressure tests for increasing the fluid densities at the predefined Re in order to extend the Ro range by reducing the mean fluid velocity. In this regard, the specific relationships between Re, Ro and Bu for an ideal gas with a gas constant, R_{c}, and viscosity, μ, are related as:
(14)
(15)
To bring the laboratory conditions into the nominal correspondence with engine conditions seen in Figure 1, the rotating speed, N, coolant pressure, P, channel hydraulic diameter, d, rotating radius, R, and walltofluid temperature difference, T_{w}T_{b}, need to be controlled for matching the realistic Re, Ro and Bu ranges. Clearly, one of the most effective approaches to extend the Ro range is to increase channel hydraulic diameter, as Ro is proportional to d^{2}. As the laboratory heat transfer test at realistic N is extremely difficult, the coolant pressure, P, is technically easier to be raised for extending Ro range; even if Ro is only proportional to P. Based on this theoretical background, the design philosophy of the present experimental method is to increase Ro by pressurizing the rotating test channel using the maximum channel hydraulic diameter allowed to fit into the rotating rig. Having increased the Ro range, the range of Bu is accordingly extended as Bu is proportional to Ro^{2}. In Figure 1, the laboratory test conditions adopted to generate the heat transfer data of rotating channels are also included^{3}^{,}^{4}^{,}^{5}^{,}^{6}^{,}^{7}^{,}^{8}^{,}^{9}^{,}^{10}^{,}^{11}^{,}^{12}^{,}^{13}^{,}^{14}^{,}^{15}^{,}^{16}^{,}^{17}^{,}^{18}^{,}^{19}^{,}^{20}^{,}^{21}^{,}^{22}^{,}^{23}^{,}^{24}^{,}^{25}^{,}^{26}^{,}^{27}^{,}^{28}^{,}^{29}. As indicated in Figure 1, the coverage of realistic engine conditions by the available heat transfer data is still limited, especially for the required Bu range. The open and the colored solid symbols depicted in Figure 1 are the pointed and fullfield heat transfer experiments, respectively. As collected in Figure 1, most of the heat transfer data with cooling applications to gas turbine rotor blades^{1}^{,}^{2}^{,}^{3}^{,}^{4}^{,}^{5}^{,}^{6}^{,}^{7}^{,}^{8}^{,}^{9}^{,}^{10}^{,}^{11}^{,}^{12}^{,}^{13}^{,}^{14}^{,}^{15}^{,}^{16}^{,}^{17}^{,}^{18}^{,}^{20}^{,}^{21}^{,}^{22}^{,}^{23}^{,}^{24}^{,}^{25}^{,}^{26} are point measurements using the thermocouple method. The wall conduction effects on measuring the wall conductive heat flux and the temperatures at fluidwall interfaces undermine the quality of heat transfer data converted from the thermocouple measurements. Also, the heat transfer measurements^{1}^{,}^{2}^{,}^{3}^{,}^{4}^{,}^{5}^{,}^{6}^{,}^{7}^{,}^{8}^{,}^{9}^{,}^{10}^{,}^{11}^{,}^{12}^{,}^{13}^{,}^{14}^{,}^{15}^{,}^{16}^{,}^{17}^{,}^{18}^{,}^{20}^{,}^{21}^{,}^{22}^{,}^{23}^{,}^{24}^{,}^{25}^{,}^{26} using the thermocouple method cannot detect the twodimensional heat transfer variations over a rotating surface. With the present experimental method^{29}^{,}^{30}^{,}^{31}^{,}^{32}, the detection of fullfield Nusselt number distributions over the rotating channel wall is permissible. The minimization of wall conduction effect using 0.1 mm thick stainlesssteel foils with Biot numbers >>1 to generate the heating power by the present experimental method permits the onedimensional heat conduction from the heating foil to the coolant flow. In particular, the acquisition of fullfield heat transfer data involving both Ro and Bu effects is not permissible using the transient liquid crystal technique and the thermocouple method. With the current steadystate liquid crystal thermography method^{19}, the detectable temperature range of 3555 °C disables the generation of heat transfer data with realistic density ratios.
Using the flow parameters governing the heat convection in a rotating channel to demonstrate that the full coverage of realistic engine conditions seen in Figure 1 has not yet been achieved, so the need for acquiring the fullfield heat transfer data at realistic engine conditions has been continuously urged. The present experimental method enables the generation of fullfield heat transfer with both Coriolisforce and rotatingbuoyancy effects detected. The protocols are aimed at assisting the investigators to devise an experimental strategy relevant to the realistic fullfield heat transfer measurement of a rotating channel. Along with the method of parametric analysis that is unique to the present experimental method, the generation of heat transfer correlation for assessing the isolated and interdependent Ro and Bu effects on Nu is permitted.
The article illustrates an experimental method aimed at generating the twodimensional heat transfer data of a rotating channel with flow conditions similar to the realistic gas turbine engine conditions but operating at much lower rotating speeds in the laboratories. The method developed to select the rotating speed, the hydraulic diameter of test channel and the range of walltofluid temperature differences for acquiring the heat transfer data at realistic engine conditions is illustrated in the introduction. The calibration tests for the infrared thermography system, the heat loss calibration tests and the operation of the rotating heat transfer test rig are shown. The factors causing the significant uncertainties for heat transfer measurements and the procedures for decoupling the Coriolisforce and buoyancy effects on the heat transfer properties of a rotating channel are described in the article with the selective results to demonstrate the present experimental method.
Subscription Required. Please recommend JoVE to your librarian.
Protocol
NOTE: The details of rotating test facilities, data acquisition, data processing and the heat transfer test module emulating an internal cooling channel of a gas turbine rotor blade are in our previous works^{29}^{,}^{30}^{,}^{31}^{,}^{32}.
1. Preparation of Heat Transfer Tests
 Formulate the experimental conditions in terms of Re, Ro and Bu from the targeted operation conditions of a gas turbine rotor blade.
 Determine the N, P, d, R, and T_{w}  T_{b} needed for acquiring the tested Re, Ro and Bu using equations (14) and (15).
 Redefine the targeting Re, Ro and Bu if N, P, d, R, and T_{w}  T_{b} exceeds the limit of the experimental facilities.
 Design and construct the scaled heat transfer test module emulating a practical internal coolant channel in a gas turbine rotor blade^{2}.
2. Determination of Thermal Emissivity Coefficient for Infrared Thermography System
 Install the calibrated thermocouple on the back side of the scanned stainlesssteel heating foil.
 Spray a thin layer of black paint on the stainlesssteel heating foil scanned by the infrared camera.
 Create symmetrical flow fields on two sides of the stainlesssteel heating foil by placing a vertical thin stainlesssteel foil in a space with the free convective flows over the two sides of the vertical heating foil.
 Feed electrical heating power through the heating foil and measure temperatures simultaneously by thermocouple and infrared thermography system from the computer display at steady state.
 Repeat step 2.4 at least four times using elevated heater powers. Ensure that the wall temperatures corresponding to the heater powers used by steps 2.3 and 2.4 cover the T_{w} range determined by step 1.2.
 Calculate the T_{w} values scanned by the infrared thermography system using a number of selective thermal emissivity coefficients for the program that converts the infrared signals into temperature data.
 Compare the T_{w} data measured by the calibrated thermocouple and the infrared thermography system at the location corresponding to the thermocouple spot with the standard deviations evaluated.
 Select the thermal emissivity coefficient with the minimum standard deviation determined by step 2.7.
 Determine the maximum precision error for the infrared thermography system using the thermal emissivity coefficient determined by step 2.8.
3. Dynamic Balance of Rotating Rig
 Install the heat transfer test module, the infrared camera, the enveloping frame and all accessories on the rotating rig.
 Adjust the counterbalancing weight gradually until the running condition of the rotating rig satisfies the vibrational limitation for the infrared thermographic measurements to exhibit the stable thermal image on the computer display.
4. Evaluation of Heat Loss Coefficients
 Fill the coolant channel of the heat transfer test module with thermal insulation material.
 Install the filled test module on the rotating test rig by fitting the test module on the rotating platform and connecting the heater power supply and all the instrumental cables.
 Activate the data acquisition system to scan the temporal T_{w} variation at a heating power until the steady state condition is satisfied. Ensure that the temporal T_{w} variations during several successive scans are less than +0.3 K at each steady state condition.
 Record the heater power, steadystate T_{w} data and the corresponding ambient temperature, T_{∞}.
 Repeat steps 4.3 and 4.4 at least five times using different heating powers at a fixed rotating speed.
 Repeat steps 4.2  4.4 with at least five rotating speeds. Ensure that the test range of the rotating speed covers all the N values determined by step 1.2.
 Repeat steps 4.3  4.6 with a reversed rotating direction.
 Construct the plots of heat loss flux against walltoambient temperature difference at each rotating speed.
 Correlate the heat loss coefficients as the functions of walltoambient temperature difference, rotating speed and direction of rotation.
 Incorporate the heat loss correlation into the post data process program for Nu accountancy.
5. Baseline Heat Transfer Tests
 Perform heat transfer tests at the targeting Reynolds numbers at zero rotating speed (Ro = N = 0) by feeding coolant flows and heater powers to the test module. Ensure the supplied coolant mass flow rate is constantly adjusted in order to control Reynolds number at the flow entry plane at the targeting value.
 Record all the relevant raw data, including the steady state wall temperatures, fluid temperatures, heater powers, flow pressures and ambient pressures and temperatures, for subsequent data processing.
 Evaluate the local and areaaveraged Nusselt numbers (Nu_{0}) over the scanned static channel walls.
6. Rotating Heat Transfer Tests
 Install the online monitoring program to monitor the test conditions at the targeting Re and Ro.
 Feed the measured coolant mass flow rate, airflow pressure, rotating speed and fluid temperature at channel entrance into the monitoring program to calculate the instant Re and Ro.
 Record all the relevant raw data, such as rotating speed, heater power, airflow and ambient pressures, as well as the wall and fluid temperatures for subsequent data processing after the predefined steadystate condition is satisfied.
 Repeat steps 6.2 and 6.3 with at least four ascending or descending heater powers at a set of fixed Re and Ro. Ensure that the test Re and Ro fall within ±1% differences from the targeting values by adjusting the rotating speed or the coolant mass flow rate or both.
 Ensure that the heat transfer tests at each set of fixed Re and Ro with different heater powers are continuously performed as the development of buoyancy induced flows is associated with the "history" of the flow development.
 Repeat steps 6.4 and 6.5 with four or five targeting Reynolds numbers (Re) at a fixed rotation number (Ro). Ensure the rotating speed is appropriately adjusted at each test Re to control both Re and Ro at the targeting values within ±1% differences.
 Repeat step 6.6 using four or five targeting rotation numbers (Ro).
 Repeat steps 6.2 to 6.7 with reversed rotating direction.
 Evaluate the local and areaaveraged Nusselt numbers (Nu) over the scanned rotating channel walls using a post data processing program.
7. Parametric Analysis
 Correlate the areaaveraged Nusselt numbers (Nu_{0}) collected from the static channel into the functions of Reynolds number.
 Evaluate the fullfield local Nu/Nu_{0} ratios at each fixed Re and Ro tested with the areaaveraged Nu/Nu_{0} ratios calculated.
 Verify the applicability of isolation Re effect by plotting the local and areaaveraged Nu/Nu_{0} ratios obtained with different Re but at identical Ro.
 Disclose the isolated impacts of rotating buoyancy on heat transfer properties of the rotating test channel by plotting the areaaveraged Nu/Nu_{0} ratios collected at the same Ro with different Re against Bu or density ratio (Δρ/ρ). Ensure the preferable selection of Bu or Δρ/ρ to construct this type of plot for obtaining the consistent data trend with a simple functional structure for heat transfer correlation.
 Extrapolate each Nu/Nu_{0} data trend collected at a fixed Ro but different Re into the limiting condition of Bu→0 or Δρ/ρ→0.
 Collect all the extrapolated Nu/Nu_{0} results with Bu→0 or Δρ/ρ→0 at all the tested Ro.
 Plot the extrapolated Nu/Nu_{0} results with vanished buoyancy interaction against Ro to disclose the uncoupled Coriolis force effects on the heat transfer properties.
 Correlate the test results collected by steps 7.4 and 7.7 into the functions of Ro and Bu.
Subscription Required. Please recommend JoVE to your librarian.
Representative Results
Realistic operating conditions for the internal coolant flows inside a rotating gas turbine blade in terms of Re, Ro and Bu are compared with the emulated laboratory conditions in Figure 1. The data points fall in the realistic engine conditions using the present experimental method summarized in the protocols^{11}^{,}^{14}^{,}^{17}^{,}^{20}^{,}^{21}. Although the fullfield heat transfer data are more useful than the pointed heat transfer data measured from the rotating channels, most of the previous heat transfer experiments adopt the thermocouple method (Figure 1). The present infrared thermography method detects the fullfield heat transfer information from a rotating surface with the buoyancyinduced flows fully developed. With the free or forced convective external flows for a static or rotating test channel, the present protocols include the generation of heat loss correlations for post data processing (Figure 2). At the top of Figure 2, the construction of the heat transfer test module is also demonstrated. The correlative coefficients for all the fitted lines shown by Figure 2 fall between 0.950.98. In view of the h_{loss} correlation seen in the plot of h_{loss} against N in Figure 2, the error bars indicate the data range determined at each rotating speed.
Figure 3, Figure 4, and Figure 5 depict the selective heat transfer results measured from the static twopass Schannel with longitudinal wavy ribs, the rotating twopass Schannel^{31} and the rotating furrowed^{32} and pinfin channel^{33}. The estimated maximum uncertainties of the Nu measurements for the static Sribbed channel, the rotating Schannel^{31}, furrowed channel^{32} and pinfin channel^{33} are 7.9%, 8.8%, 9.2%, and 9.7%, respectively. To disclose the Re impact on the heat transfer properties of a coolant channel, the baseline fullfield heat transfer data detected from the static channel by the present infrared thermography method as typified by Figure 3 are essential. The diagram shown at the top of Figure 3 also depicts the channel configuration of the twopass Schannel with the longitudinal wavy ribs. The channel section is square with the semicircular sectioned longitudinal wavy ribs on two opposite heated walls of the inlet and outlet legs.
The applicability of isolated Re impact from Ro and Bu effects on local and regionally averaged heat transfer is permitted by presenting the heat transfer data in terms of Nu/Nu_{0} (Figure 4). Both patterns and levels of Nu/Nu_{0} at the same Ro with similar Bu seem to be weak functions of Re (Figure 4). The typical results from the protocol for disclosing the isolated Coriolis force effects on heat transfer properties are demonstrated in Figure 5. In Figure 5, the variations of Nu/Nu_{0} at each fixed Ro against Bu for two different rotating channels with wavy endwalls^{32} and diamond shaped pinfins^{33} tend to follow linearlike data trends. Thus, the linear extrapolation when Bu→0 is selected for the identified Nu/Nu_{0} levels at Bu = 0 and Ro>0. But, due to the different channel configurations, the Nu/Nu_{0} ratios measured from the rotating furrowed^{32} and pinfin^{33} channels as depicted in Figure 5 are respectively decreased and increased by raising Bu. In this regard, the depiction of Nu/Nu_{0} variations against density ratio (Δρ/ρ)^{3}^{,}^{4}^{,}^{5}^{,}^{6}^{,}^{34} has often led to the nonlinear Nu/Nu_{0} variations. Thus, the extrapolation of each Nu/Nu_{0} data trend at a fixed Ro toward the asymptotic limit of Δρ/ρ→0 with diminished buoyancy effect along a nonlinear data trend is often affected by the type of correlative function selected. Nevertheless, the data extrapolating procedure for the heat transfer results detected from the leading and trailing walls of the rotating channels^{32} demonstrates the applicability to unravel the isolated Coriolis force effects on heat transfer properties with vanished buoyancy interaction at Bu=0 (Figure 5).
The socalled zerobuoyancy Nu/Nu_{0} ratios are only controlled by Ro to reflect the isolated Coriolis force effects. The manner of heat transfer variations from the staticchannel references disclosed by steps 7.7 and 7.8 is typified by Figure 6. The separated Ro impact from the buoyancy effect on the heat transfer performances of a rotating channel is correlated as the Ro function to be a part of Nu/Nu_{0} correlation (Figure 6). The positive or negative ψ_{2} values in Figure 6 indicate the improving or impeding effects on heat transfer performances due to buoyancy interactions. The larger ψ_{2} magnitude, the higher degrees of rotating buoyancy impact are imposed on the heat transfer properties. The fitted lines indicated in Figure 6 are the plots of the correlative functions. The functional structures of the correlations for zerobuoyancy Nu/Nu_{0} ratios and ψ_{2} values are generally determined in accordance with the varying manners of the data trends emerged in Figure 6. As discussed previously, the different channel geometries between the furrowed^{32} and pinfin^{33} channels have respectively led to the negative and positive ψ_{2} values in Figure 6. But the common feature of the reduced magnitudes of ψ_{2} values caused by increasing Ro is observed for the two types of rotating channels^{32}^{,}^{33} in Figure 6. Having correlated the ψ_{2} values and the Nu/Nu_{0} ratios at zerobuoyancy conditions into the Ro functions, the heat transfer correlations, which permit the evaluation of the isolated and coupled Ro and Bu effects on Nu/Nu_{0}, is generated for the particular rotating channel.
Figure 1. Realistic operating Re, Ro and Bu ranges and the emulated laboratory conditions for a rotating coolant channel in a gas turbine rotor blade. The test conditions performed by NASA HOST program^{3}^{,}^{4}^{,}^{5}^{,}^{6} are indicated as the bar symbol. The open and solid symbols respectively signify the Bu, Ro, and Re test ranges for the pointed and fullfield heat transfer measurements. Numbers in brackets are references from which data are taken. Please click here to view a larger version of this figure.
Figure 2. Typical heat loss coefficients (h_{loss}) at various rotating speeds^{30} using the trapezoidal twinpass ribroughened rotating channel as an illustrative example. The diagram at the top depicts the constructional details of the rotating test module. The slope of each data trend constituted by the heat loss flux against the walltoambient temperature difference shown in the left lower portion reveals the heat loss coefficient at the specific rotating speed. By correlating the heat loss coefficients detected at all the rotating speed tested, the generated heat loss correlation typified by the right lower plot is incorporated into the data processing program for Nu accountancy. The error bars in the lower right plot indicate the ranges of h_{loss}^{30}. Please click here to view a larger version of this figure.
Figure 3. Local Nusselt number distribution of the static twinpass Schannel roughened by curly ribs at Re = 15,000 measured by present infrared thermography method. The top diagram depicts the endwall of the twopass wavy channel and the longitudinal Sribs. As indicated by the AA’ section view, the pair of longitudinal Sribs is arranged inline on two opposite channel endwalls. In the detailed distribution of Nusselt number over the twopass wavy endwall shown as the lower plot, the Nu data along the two longitudinal Sribs are discarded due to the wall conduction effects on the distributions of heatflux and walltemperature. Please click here to view a larger version of this figure.
Figure 4. Examples demonstrating the isolation of Re impact from Ro and Bu effect on local and regionallyaveraged heat transfer properties of rotating channel. The upper portion exhibits the detailed Nusselt number distributions at a fixed Ro of 0.15 with a different Re of 5000, 7500, and 12,500 to enlighten the impacts of Reynolds number on the heat transfer properties of the rotating endwall. The lower portion depicted the areaaveraged heat transfer properties over the rotational leading and trailing endwalls. The normalized Nu/Nu_{0} ratios highlight the heat transfer variations from the nonrotating scenarios by rotation. Adapted with permission from Chang et al. 2017^{31}. Please click here to view a larger version of this figure.
Figure 5. Examples demonstrating the uncoupled Ro effect from Bu impact on heat transfer properties of rotating channel^{32}^{,}^{33}. Each Budriven Nu/Nu_{0} variation is obtained at the fixed Ro and correlated as a linear function of Bu as indicated by the straight line in each plot. The correlation coefficients of these fitted lines fall between 0.96 and 0.98. The extrapolation of the Nu/Nu_{0} data trend toward Bu→0 along each fitted line reveals the Nu/Nu_{0} ratio at the tested Ro. The magnitude and slope of each Budriven Nu/Nu_{0} data trend disclose the manners of buoyancy effect on heat transfer performances. The magnitudes of the slopes represent the degrees of Bu impact on Nu/Nu_{0}. The positive and negative slopes respectively reflect the improving and impairing buoyancy impact on heat transfer levels. Numbers in brackets are references from which data are taken. Please click here to view a larger version of this figure.
Figure 6. Uncoupled Ro and Bu effects on regionally averaged heat transfer performances of the rotating wavy channel^{32}^{,}^{33}. The upper portion collects the heat transfer scenarios at various Ro but with vanished buoyancy effect at Bu = 0. Such Nu/Nu_{0} variations are solely caused by the various Coriolis forces at different Ro. The lower portion shows the variations of Bu impact on Nu/Nu_{0} at different Ro. The negative and positive ψ_{2} values indicate the respective impairing and improving Bu impacts on the heat transfer performances for the furrowed^{32} and pinfin^{33} channels. The dotted lines in this Figure are the correlation results for Nu/Nu_{0} at Bu = 0. Numbers in brackets are references from which data are taken. Please click here to view a larger version of this figure.
Subscription Required. Please recommend JoVE to your librarian.
Discussion
While the endwall temperatures of a rotating channel are detected by an infrared thermography system, the fluid temperatures are measured by thermocouples. As the alternative magnetic field of an AC motor that drives a rotating rig induces electrical potential to interfere the thermocouple measurements, the DC motor must be adopted to drive a rotating test rig.
The fluid temperature distribution over the exit plane of a heated channel is not uniform. At least five thermocouples on the existing plane of a rotating channel are recommended for measuring the local fluid exit temperatures. In particular, these thermocouples measuring the fluid temperatures installed in the flow passage are subject to centrifugal forces during the rotating tests. The thermocouple wires are easily bent toward the hot channel walls. Thus, a shielded thermocouple cable for measuring the fluid entry temperature is used. On the flow exit plane, a mesh with several thermocouple beads weaved on the mesh can be sandwiched between the exit flanges of a test channel to detect the fluid exit temperatures at the predefined locations under a rotating test condition.
With considerable rotation induced buoyancy effects on the flow and heat transfer characteristics of a rotating channel, the method selected to detect the fullfield heat transfer data needs to include both Coriolisforce and buoyancy effects. Using the transient liquid crystal method for measuring the fullfield heat transfer data, the thermal boundary layers are not yet fully developed as the temporal channelwall temperature variations are essential by this method for acquiring the convective heat transfer coefficients. As the centripetal acceleration could reach 10^{5 }x g in a coolant channel of a rotating gas turbine blade, the heat transfer data subject to the influences of the fully developed buoyancy flows, which are detectable by the present experimental method, are more practical for design activities.
The exposure of the scanned hot channel wall to an infrared camera inevitably incurs heat loss from the Joule heat generated by the heating foils. The protocols for conducting the heat loss calibration tests are critical for ensuring the quality of heat transfer data. Inheriting either from the free or forced convective external flows for a static or rotating test channel, the convective heat transfer coefficients can be correlated as the function of walltoambient temperature difference at a fixed rotating speed (Figure 2). It is preferable to envelop the entire rotating heat transfer test module with a shield for recovering the "freeconvective" like external flows during the rotating tests. The maximum experimental uncertainties of heat transfer data are generally reduced when the percentage of the heat loss flux from the supplied heat flux is reduced. Nevertheless, the heat loss coefficients are slightly increased as N increases even with the enveloped shield covering the entire heat transfer test module (Figure 3). The heat loss correlation is included in the post data processing program to evaluate the distribution of local heat loss flux for each set of heat transfer test results. As the thermal inertia of the heat transfer module filled by thermal insulation material is considerably increased, the time required for reaching the steadystate condition during each heat loss test is considerably extended from a heat transfer test with airflow.
It is essential to investigate the applicability of the isolating Re effect on heat transfer properties from those induced by rotation. As the Re effect on heat transfer performances depends on the channel configurations, it is not appropriate to customarily adopt the heat transfer correlations generated from other channel geometries as the staticchannel heat transfer references. The present experimental method isolates Re impact from Ro and Bu effects by presenting the heat transfer data in terms of Nu/Nu_{0}, in which the Nu_{0}_{ }data are measured for the static test channel. While the buoyancy effect in a rotation channel with centripetal acceleration about 10^{5 }x g is considerable, the gravitationdriven buoyancy effect on the heat transfer property of a static channel is generally negligible within the typical range of fluid density ratios examined for a static test channel.
During a heat transfer test after feeding heater power to generate the required temperature gradients for facilitating heat convection, a certain degree of buoyancy effect driven by the induced centripetal acceleration field in the rotating channel is inevitable. Such coupled Ro and Bu effects for a rotating channel at the realistic engine conditions are not negligible due to the extremely high centripetal accelerations. Thus, both Coriolis force and rotating buoyancy level are simultaneously altered when the rotating speed is adjusted. The simultaneous control of Ro and Re at the targeting values during the rotating experiment is essential for decoupling the Ro and Bu effects on heat transfer properties. Having fixed both Ro and Re, the heat transfer variations corresponding to the variation of heat flux, or buoyancy level, reflect the rotating buoyancy effect on heat transfer properties at the tested Ro. The Nu/Nu_{0} data converted from the data set generated in this manner permit the implementation of steps 7.4  7.8 for identifying the Coriolis force effect and rotating buoyance effect in isolation.
The Bu impact on the heat transfer property of a rotating channel is often Ro dependent as exemplified by Figure 6 in which the ψ_{2} values are varied as Ro changes. It is not appropriate to select the mathematic structure of the heat transfer correlation that treats the Ro and Bu as the independent parameters in the correlation.
In view of the Nu/Nu_{0} extrapolation toward the limiting condition of Bu→0, the linearlike Nu/Nu_{0} variations against the selected buoyancy parameter is preferable in order to reduce the uncertainty caused by the data extrapolation. In this regard, the fluid density ratio, Δρ/ρ or the buoyancy number, Bu, is recommended as the buoyancy parameter for disclosing the zerobuoyancy Nu/Nu_{0} level during such data extrapolating process.
With high pressure rotating tests, the deformations of heating foils and the constituent components of a rotating channel due to the thermal expansions at various patterns of temperature distribution often cause airflow leakage during the rotating test. Such small airflow leakage is difficult to be identified during the rotating test. Thus, immediate subsequent data processing is recommended for acquiring the heat transfer data of the rotating channel. By crossexamining the heat transfer results obtained from the previous rotating tests, the implication of any inconsistent data trend is the possible airflow leakage. The subsequent measures to detect and then prevent the airflow leakage are required.
We have demonstrated a method for generating the heat transfer data of a rotating channel at the realistic engine conditions with the Coriolisforce effect and rotating buoyancy effect uncoupled. The major limitation of the present experimental method for extending the test ranges of Ro and Bu is the sustainability of the infrared camera that rotates with the test channel. In general, 10 x g is the maximum sustainable centrifugal acceleration for an infrared camera. With respect to the existing method detecting the heat transfer rates of a rotating channel, the use of thin heating foil can minimize the effects of channelwall conduction on the distribution of local convective heat flux and the detection of temperatures at wallfluid interfaces. Also, the twodimensional fullfield heat transfer distribution over a rotating surface subject to the steadystate buoyancy effect are detectable using the present experimental technique. With the data analysis method developed, the influences of Coriolis force and rotating buoyancy on the fullfield heat transfer property of a rotating channel can be uncoupled. This method has already been applied to a wide range of rotating channel configurations. We expect that the present experimental strategy can lead to the design friendly heat transfer correlations and which will continue to extend for the full coverage of realistic engine conditions when the advancement of infrared camera technology permits its usages at the conditions with higher centrifugal accelerations.
Subscription Required. Please recommend JoVE to your librarian.
Disclosures
The authors have nothing to disclose.
Acknowledgments
The present research work was financially sponsored by the Ministry of Science and Technology of Taiwan under the grant NSC 942611E022001, NSC 952221E022018, NSC 962221E022015MY3 and NSC 972221E022013MY3.
Materials
Name  Company  Catalog Number  Comments 
Rotating test rig  Inhouse made  Design by this research group  
Heat transfer test module  Inhouse made  Design by this research group  
Mass flow meter  Eldride Product, Inc.  31003010101 3591007 

Infrared thermography system  NEC P384A8  310040104 3127A4 

Instrumentation slip ring  Michigan Scientific SR36M  310050662 3553372 
References
 Morris, W. D. Heat transfer and fluid flow in rotating coolant channels. John Wiley and Sons. ISBN 0471101214 (1981).
 Morris, W. D. A rotating facility to study heat transfer in the cooling passage of turbine rotor blades. Journal of Power and Energy. 210, (1), 5563 (1996).
 Wagner, J. H., Johnson, B. V., Graziani, R. A., Yeh, F. C. Heat transfer in rotating passages with smooth walls and radially outward flow. ASME Journal of Turbomachinery. 113, (1), 4251 (1991).
 Wagner, J. H., Johnson, B. V., Kopper, F. C. Heat transfer in rotating serpentine passages with smooth walls. ASME Journal of Turbomachinery. 113, (3), 321330 (1991).
 Wagner, J. H., Johnson, B. V., Steuber, G. D., Yeh, F. C. Heat transfer in rotating serpentine passages with trips normal to the flow. ASME Journal of Turbomachinery. 114, (4), 847857 (1992).
 Johnson, B. V., Wagner, J. H., Steuber, G. D., Yeh, F. C. Heat transfer in rotating serpentine passages with selected model orientations for smooth or skewed trip walls. ASME Journal of Turbomachinery. 116, (4), 738744 (1992).
 Hwang, G. J., Tzeng, S. C., Mao, C. P., Soong, C. Y. Heat transfer in a radially rotating fourpass serpentine channel with staggered halfv rib turbulators. ASME Journal of Heat Transfer. 123, (1), 3950 (2001).
 Azad, G. S., Uddin, M. J., Han, J. C., Moon, H. K., Glezer, B. Heat transfer in a twopass rectangular rotating channel with 45deg angled rib turbulators. ASME Journal of Turbomachinery. 124, (2), 251259 (2002).
 Griffith, T. S., AlHadhrami, L., Han, J. C. Heat transfer in rotating rectangular cooling channels (AR=4) with angled ribs. ASME Journal of Heat Transfer. 124, (4), 617625 (2002).
 AlHadhrami, L., Griffith, T. S., Han, J. C. Heat transfer in twopass rotating rectangular channels (AR=2) with five different orientations of 45 deg Vshaped rib turbulators. ASME Journal of Heat Transfer. 125, (2), 232242 (2003).
 Chang, S. W., Liou, T. M., Hung, J. H., Yeh, W. H. Heat transfer in a radially rotating squaresectioned duct with two opposite walls roughened by 45 deg staggered ribs at high rotation numbers. ASME Journal of Heat Transfer. 129, (2), 188199 (2007).
 Zhou, F., Lagrone, J., Acharya, S. Internal cooling in 4:1 AR passages at high rotation numbers. ASME Journal of Heat Transfer. 129, (12), 16661675 (2007).
 Liu, Y. H., Huh, M., Han, J. C., Chopra, S. Heat transfer in a twopass rectangular channel (AR=1:4) under high rotation numbers. ASME Journal of Heat Transfer. 130, (8), (2008).
 Chang, S. W., Liou, T. M., Chiou, S. F., Chang, S. F. Heat transfer in highspeed rotating trapezoidal duct with ribroughened surfaces and air bleeds from the wall on the apical side. ASME Journal of Heat Transfer. 130, (6), (2008).
 Wright, L. M., Liu, Y. H., Han, J. C., Chopra, S. Heat transfer in trailing edge, wedgeshaped cooling channels under high rotation numbers. ASME Journal of Heat Transfer. 130, (7), 111 (2008).
 Liou, T. M., Chen, M. Y., Tsai, M. H. Fluid flow and heat transfer in a rotating twopass square duct with inline 90deg ribs. ASME Journal of Turbomachinery. 124, (2), 260268 (2002).
 Chang, S. W., Liou, T. M., Yang, T. L., Hong, G. F. Heat transfer in radially rotating pinfin channel at high rotation numbers. ASME Journal of Turbomachinery. 132, (2), (2010).
 Rallabandi, A., Lei, J., Han, J. C., Azad, S., Lee, C. P. Heat transfer measurements in rotating bladeshape serpentine coolant passage with ribbed walls at high Reynolds numbers. ASME Journal of Turbomachinery. 136, (9), (2014).
 Mayo, I., Arts, T., Ahmed, E. H., Parres, B. Twodimensional heat transfer distribution of a rotating ribbed channel at different Reynolds numbers. ASME Journal of Turbomachinery. 137, (3), (2015).
 Chang, S. W., Yang, T. L., Liou, T. M., Fang, H. G. Heat transfer in rotating scaleroughened trapezoidal duct at high rotation numbers. Applied Thermal Engineering. 29, (8), 16821693 (2009).
 Liou, T. M., Chang, S. W., Chen, J. S., Yang, T. L., Lan, Y. A. Influence of channel aspect ratio on heat transfer in rotating rectangular ducts with skewed ribs at high rotation numbers. International Journal of Heat Mass Transfer. 52, (23), 53095322 (2009).
 Huh, M., Liu, Y. H., Han, J. C. Effect of rib height on heat transfer in a two pass rectangular channel (AR = 1:4) with a sharp entrance at high rotation numbers. International Journal of Heat Mass Transfer. 52, (19), 46354649 (2009).
 Xu, G., Li, Y., Deng, H. Effect of rib spacing on heat transfer and friction in a rotating twopass square channel with asymmetrical 90deg rib turbulators. Applied Thermal Engineering. 80, (5), 386395 (2015).
 Tao, Z., Yang, M., Deng, H., Li, H., Tian, S. Heat transfer study in a rotating ribbed twopass channel with enginesimilar cross section at high rotation number. Applied Thermal Engineering. 106, (5), 681696 (2016).
 Li, Y., Deng, H., Tao, Z., Xu, G., Chen, Y. Heat transfer characteristics in a rotating trailing edge internal cooling channel with two coolant inlets. International Journal of Heat Mass Transfer. 105, (2), 220229 (2017).
 Deng, H., Chen, Y., Tao, Z., Li, Y., Qiu, L. Heat transfer in a twoinlet rotating rectangular channel with sidewall fluid extraction. International Journal of Heat and Mass Transfer. 105, (2), 525534 (2017).
 You, R., Li, H., Tao, Z., Wei, K. Heat transfer investigation in a smooth rotating channel with thermography liquid crystal. ASME Turbo Expo. GT201656413, Turbomachinery Technical Conference and Exposition: Heat Transfer. 5, (B), V05BT16A006 1~10 (2016).
 Morris, W. D., Chang, S. W. An experimental study of heat transfer in a simulated turbine blade cooling passage. International Journal of Heat Mass Transfer. 40, (15), 37033716 (1997).
 Chang, S. W., Liou, T. M., Po, Y. Coriolis and rotating buoyancy effect on detailed heat transfer distributions in a twopass square channel roughened by 45° ribs at high rotation numbers. International Journal of Heat Mass Transfer. 53, (7), 13491363 (2010).
 Wang, W. J. Heat transfer in rotating twinpass trapezoidalsectioned passage with two opposite walls roughened by 45 degree ribs. Department of Marine Engineering, National Kaohsiung Marine University. Msc thesis (2006).
 Chang, S. W., Wu, P. S., Chen, C. S., Weng, C. C., Jiang, Y. R., Shih, S. H. Thermal performance of radially rotating twopass Sshaped zigzag channel. International Journal of Heat and Mass Transfer. 115, (B), 10111031 (2017).
 Chang, S. W., Lees, A. W., Liou, T. M., Hong, G. F. Heat transfer of a radially rotating furrowed channel with two opposite skewed sinusoidal wavy walls. International Journal of Thermal Sciences. 49, (5), 769785 (2010).
 Chang, S. W., Liou, T. M., Lee, T. H. Heat transfer of a rotating rectangular channel with a diamondshaped pinfin array at high rotation numbers. Journal of Turbomachinery Transactions of the ASME. 135, (4), 041007 1~10 (2013).
 Morris, W. D., Chang, S. W. Heat transfer in a radially rotating smoothwalled tube. The Aeronautical Journal. 102, (1015), 277285 (1998).