The laminar and turbulent regimes of a boundary layer on a flat plate are often represented with separate correlations under the assumption of a distinct “transition Reynolds number.” Average heat coefficients are then calculated by integrating across the “transition point.” Experimental data do not show an abrupt transition, but rather an extended transition region in which turbulence develops. The transition region may be as long as the laminar region. Although this transitional behavior has been known for many decades, few correlations have incorporated it. One attempt was made by Stuart Churchill in 1976. Churchill, however, based his curve fit on some doubtful assumptions about the data sets. In this paper, we develop different approximations through a detailed consideration of multiple data sets for , and varying levels of freestream turbulence for smooth, sharp-edged plates at zero pressure gradient. The result we obtain is in good agreement with the available measurements and applies smoothly over the full range of Reynolds number for either a uniform wall temperature or a uniform heat flux boundary condition. Fully turbulent air data are correlated to . Like Churchill's result, this correlation should be matched to the estimated transition condition of any particular flow. We also review the laminar analytical solutions for a uniform wall heat flux, and point out limitations of the classical Colburn analogy.
Simplified treatments of boundary layer heat transfer split the boundary layer into an upstream laminar section and downstream turbulent section. Correlations are then formed for these two sections, which are assumed to be separated by a distinct transition Reynolds number. An example is the correlation for average heat transfer coefficient proposed by Whitaker : the heat transfer coefficient is averaged over the length, using results from laminar theory and correlation for turbulent flow. A transition Reynolds number of 200,000 is assigned, and single equation depending on the overall Reynolds number is obtained.
The actual transition process has long been known to include a lengthy transition region, in which the heat transfer coefficient rises smoothly. Indeed, data show this region to have a length similar to the laminar region (see Sec. 4).
In Eq. (1), all four terms are included irrespective of the value of . Each term will become dominant in the appropriate range, as ensured by the exponents used and the terms' differing dependence on . Churchill's values of p = –5 and s = 2 produce fairly gradual transitions between the three indicated ranges. To account for the variability of transition, his result must be matched to any particular dataset by fitting the value of Reynolds number at the upper limit of the transition region, , into his formula for .
Churchill expressed the three variable terms as functions of a single parameter, , with different values of d for uniform wall temperature (UWT) or for uniform heat flux (UHF)
This function was based upon his and Ozoe's laminar flow curve-fits spanning all values of Pr [2,4]. As a consequence, he had to make an approximation for the turbulent Nusselt number. The full statement of Churchill's equation is given in Appendix A.
Churchill compared his model to aggregated data for the local and average heat transfer coefficient, with partial agreement. As seen in Fig. 1, the overall agreement with local data is not close. The transition away from laminar flow begins too soon, the slope in the transition region does not match, and the values in the fully turbulent range are slightly high. Therefore, a more accurate fit is desirable.
Here, we develop a correlation that accurately captures the lengthy transition region. The available data sets are reviewed together with key theoretical points so that differences among the measurements will be more clear. Churchill's basic concept is retained, but its primary elements are reëvaluated. The resultant equation is compared to a large number of independent measurements. The measurements clearly indicate that the transition Nusselt number rises more rapidly with Reynolds number when transition begins at a higher Reynolds number. We offer a straightforward formula for calculating the local (and average) heat transfer coefficient, which has good agreement with data in the literature.
In any situation for which both the laminar and turbulent regimes must be considered, so should be the transitional region. This fact is particularly important if local values of the heat transfer coefficient are required.
2 Theoretical Considerations and Assessment of Available Data Sets
The theory of boundary layer heat transfer is today very highly developed. To understand the differences in the data sets, a few fundamental points need to be in mind. These ideas are highlighted briefly, with no intention to give a comprehensive account.
2.1 Laminar Boundary Layers.
Most differences between the laminar measurements used herein and the laminar theory can be attributed to differences in the wall boundary condition and/or in the upstream initial conditions.
The heat transfer coefficient in a laminar boundary layer is strongly sensitive to the upstream history of the flow, including both the leading-edge configuration and any variation in wall temperature. Most of the experimental reports have taken care to describe leading-edge conditions. In some studies, upstream boundary layers were suctioned away through a slot ahead of the test piece. Most of the experiments had an unheated initial length or an initial section in which temperatures were not well controlled, and most of those experimentalists applied analytical correction factors to make their reported data appear as if the initial length were uniform (Blair  is notable for giving a clear account of the unheated region and reporting it without adjustment; see Sec. 4.).
The local heat transfer coefficient in a flat-plate boundary layer is higher for UHF than for UWT, by the factor of 0.453/0.332 = 1.36 . While that knowledge, as an experimental result with partial theoretical support, was available in the 1930s , many noteworthy heat transfer researchers were unaware of this fact during the 1940s and 1950s (see Appendix B). Instead, these researchers benchmarked laminar data on Pohlhausen's UWT solution .1 For example, Jakob and Dow  questioned the work of Fage and Falkner  on electrically heated platinum strips for being 27% above Pohlhausen's result; further, they reported that A.P. Colburn came to a similar opinion after they wrote to him. Likewise, Kestin et al.  were very critical of the transient measurements of Sugawara et al.  for this same reason. Yet, Sugawara's laminar measurements (and correlation) agree well with the analytical formula for a uniform wall heat flux; and the 1951 Japanese-language version of their work  explicitly considers a wall temperature varying as (as it would for a uniform heat flux), contrasting that result to Pohlhausen's formula.2
2.2 Freestream Turbulence.
Freestream turbulence is normally manipulated by placing a grid of bars upstream of the test plate; the formation and decay of grid turbulence have been very well understood in the wind-tunnel literature since the late 1930s [13,14].
Experimentalists concluded early on that freestream turbulence at levels up to 5% or so had little or no effect on the laminar heat transfer coefficient [7,15]. Experiments also clearly showed that increasing freestream turbulence strongly reduced the transition Reynolds number and that, with extremely low turbulence, the onset of transition could be delayed to . However, the effect of freestream turbulence on heat transfer in the turbulent boundary layer was extensively debated during the mid-twentieth century [7,10,11,15,17]. Mainly, these early investigations concluded that freestream turbulence had little effect at zero pressure gradient (the work of Sugawara et al.  is an exception). Eventually, more precise measurements encompassing higher turbulence levels showed that freestream turbulence can modestly increase the turbulent heat transfer coefficient. Blair , for example, measured up to an 18% increase for a nominal freestream turbulence level of 6%.
In what follows, we consider data with reported freestream turbulence levels below 5% (and mostly below 3%).
The various experimental reports differ significantly on how turbulence was characterized. Blair, for example, used multicomponent hot-wire anemometer measurements to fully describe the streamwise variation of the turbulence in each case. In other cases, experimentalists have reported only, say, “about 0.3%” without elaboration . The length and time scales of turbulence in the various experimental systems differ and cannot productively be compared using only the reported value of . However, variation of turbulence within a single experimental system is clearly comparable within that system. We note the level of freestream turbulence of the data sets whenever it has been reported.
2.3 Turbulent Boundary Layers.
Heat transfer coefficients in turbulent boundary layers are not sensitive to slow streamwise variations in the wall temperature, and the value of h for UHF is only about 4% greater than for UWT .
The near-wall velocity distribution in the fully turbulent boundary layer is essentially a function of the local wall shear stress, as embodied by the universal velocity profile or “law-of-the-wall” [20,21]. Similarly, the local heat transfer coefficient for nonmetallic fluids is mainly determined by the local shear stress and the Prandtl number, as shown by the Reynolds–Colburn analogy [22,23] and its more accurate generalizations based on boundary-layer structure [6,24,25]. For the internal flow case, shear stress is represented by the Darcy friction factor, proportional to , while the external flow case uses the skin friction coefficient, as .
Data for external, turbulent boundary-layer heat transfer have mainly been acquired in air flows. Water boundary layers have been measured, apparently, only by Žukauskas and Šlančiauskas  and by Hollingsworth . Data for other liquids are even more limited, amounting to measurements in “transformer oil,” again by Žukauskas and coworkers. This situation could seem discouraging until one recalls that a vast body of data exists for turbulent internal flow (in pipes), spanning an enormous range of Prandtl number. As a consequence of the universal velocity distribution, the generalized Reynolds–Colburn analogy for pipe flows is directly adaptable to the external boundary layers on a flat plate in nonmetallic, low-speed flows. Thus, the exhaustive correlations for internal flow in terms of can support the calculation of h from for those Pr at which external flow measurements have not been made (i.e., for most fluids other than air). We shall reinforce this point in this paper by comparison to data for air, water, and transformer oil.
Various authors have proposed values for the coefficients . For example, Prandtl  chose (pipe flow); Žukauskas and Šlančiauskas  proposed (flat plate); White  gave (flat plate); and Petukhov  gave (pipe flow). Gnielinski  examined thousands of data points in pipe flow, spanning , leading him to suggest . Gnielinski's values capture 90% of the liquid data to , with even better agreement for gases. In what follows, we adopt Gnielinski's coefficients.
2.4 Variable Property Effects.
For boundary layers with large temperature differences, property variation across the boundary may be important. For turbulent gas flows, property corrections are often made by evaluating a correlation at and multiplying the result by the absolute temperature ratio, , with n between 0.25 and 0.4 [26,32]. (For pipe flow, Gnielinski recommends n = 0.45 ). For the air data discussed here, temperature ratios are less than 1.11, and property variations are either negligible or small (). No corrections are applied in plotting the correlations for air data in what follows.
For liquids, corrections are generally based on a ratio of viscosities or Prandtl numbers, and most data are for turbulent pipe flows. For external boundary layers, the most comprehensive experimental study is due to Žukauskas and Šlančiauskas , who recommend evaluating Nusselt number correlations using properties at and multiplying the result by . (For turbulent pipe flow, Gnielinski recommends an exponent of 0.11). Typical values of the Prandtl number for the fluids discussed here are given in Appendix D.
For the water data used herein, the property corrections are modest, ranging from to 6%, with a single run (discussed below) reaching 10%. No corrections have applied to the water correlations in the plots shown.
For transformer oil, however, the corrections can exceed 20%; these will be addressed in context.
2.5 Data Sets Examined.
For this study, measurements of the local heat transfer coefficient through transition region are of primary interest. The data sets for local heat transfer coefficients are described in Tables 1 and 2. These data are used in examining the variation of the heat transfer coefficient in the transition region.
Data sets for average heat transfer coefficients are summarized in Table 3 because they were used by Churchill , Whitaker , and others in developing correlations. Apart from the laminar study of Fage and Falkner , those data sets are strongly affected by unheated starting lengths and other geometric differences, so we will not consider them in any detail.
3 Reconstituting Churchill's Formula
The term in Eq. (1) is intended to account for heat transfer at very low Reynolds numbers, as for creeping flow. Churchill selected on the basis of numerical simulations for by Dennis and Smith  and without experimental support. The contribution of is negligible for Reynolds numbers greater than 1000 or so in nonmetallic liquids. We shall restrict our attention to nonmetallic liquids, with and , and omit .
where a = 0.332 for UWT or 0.453 for UHF .
Churchill suggested a fixed exponent of c = 3/2, which seemed to fit some of Žukauskas and Šlančiauskas' data sets ; the latter authors had suggested 1.4. In contrast, our comparison to data in Sec. 4 shows c is usually larger than this and that the value c rises as rises. Our approach will be to treat c as a function of in fitting the data.
White recommends this equation for zero-pressure-gradient turbulent boundary layers at any .
The three component terms should be calculated from Eqs. (3) and (5), and either Eqs. (6) or (8), rather than Churchill's constructions. The value of should be fit to the dataset, and the value of c may either be fitted or estimated from Eq. (11), which is discussed in Sec. 4.4. When data are not available, Eq. (12) can provide a very rough estimate of .
Equation (9) can conveniently be used in a spreadsheet or other software… or by hand if one is old school.
4.1 Data for Air.
Figure 2 shows wind tunnel data including the transition region for UWT from Kestin et al.  and from Junkhan and Serovy . Values of c and have been fitted to each of these four data sets, and the corresponding curve from Eq. (9) (with Eqs. (3), (5), (6), and (7)) is plotted. The value of c trends higher as rises. Also shown, for reference, are two sets of fully turbulent measurements from Reynolds et al. ; the latter are in excellent agreement with both Eqs. (6) and (8). Reynolds et al. reported six additional data sets, which are substantially the same as those shown in Fig. 2.
Figure 3 shows wind tunnel measurements by Blair for several carefully characterized levels of freestream turbulence with UHF. Blair's test plate had an unheated starting length, cm, so that the laminar Nusselt number follows the theoretical prediction, (see Appendix B), which we use as the laminar term in Eq. (9) for this case. Agreement of the laminar theory with the data is excellent. Blair's fully turbulent data converge to Eq. (6). Values of c and were fitted as before, and again the agreement with data is excellent. The value of c becomes much larger for . As noted earlier, Blair concluded that freestream turbulence could increase the turbulent heat transfer coefficient, an effect that may be discernible in the data for .
Figure 4 shows fits through the transition region for four additional wind tunnel studies. Figure 4(a) shows data for seven different freestream speeds from Seban and Doughty . The trend of these data for Reynolds numbers above 106 is 7–8% below Eq. (6). In the laminar range, the data lie above Eq. (3); however, these data were analytically corrected for an unheated starting length, a process that the authors also note contributed to greater scatter in the transition range. Beyond these two differences, the fitted values of c and put Eq. (9) into good agreement with the transition data.
Figure 4(b) shows the one dataset in the study of Reynolds et al.  for which the boundary layer had not been intentionally tripped at the leading edge. Again using the fitted parameters, Eq. (9) has good agreement with the measurements. Equation (8), proposed by Reynolds et al. for fully turbulent flow, is also shown. The data can be observed to exceed the fully turbulent line just after transition. Reynolds et al. discussed and modeled this effect, which had also been described by Seban and Doughty. They attributed the effect to differences between the momentum and energy thicknesses of an already fully turbulent boundary layer and of the boundary layer formed at the end of transition (A similar, but less pronounced effect, might be perceived in Figs. 3 and 4(a).).
Figure 4(c) shows measurements from Žukauskas and Šlančiauskas  for a single plate at a fixed air speed with an increasing wall heat flux. The transition Reynolds number shifts slightly lower as the plate temperature increases. Equation (9) fits each case separately.
Figure 4(d) shows additional measurements from Žukauskas and Šlančiauskas from a different series of experiments. In the laminar range, the two data sets labeled Table 244 trend above laminar theory and have some significant scatter; the data labeled Table 25 have a large gap around the transition range are very scattered in the laminar range. The unusual laminar behavior might indicate a leading-edge influence of some unknown type (see Sec. 2.1). In the fully turbulent range, all three data sets lie in a fairly tight line about 15–25% above both Eq. (9) and measurements by other investigators; no reason for this difference is evident. However, the transition range of the Table 24 data sets is well represented by a slope c = 2.5, consistent with the other air data undergoing transition at similar Reynolds numbers.
together with the relationship proposed by Reynolds et al. . Equation (10) is substantially higher than other correlations and data for air, presumably because is it based on the data of Fig. 4(d), which are also higher. Similarly, the value proposed by Žukauskas and Šlančiauskas for the term a1 in Eq. (2) is lower than in other investigations.
Žukauskas and Šlančiauskas provide an additional dataset for high heat fluxes (not shown), with of 1.20–1.25, for which the variable properties correction would be 7–9%. Those data do not encompass the transition range. Their trend is similar to Fig. 4(d), but with somewhat more scatter.
To evaluate the performance of Eq. (6), all of the fully turbulent air data from Refs. , , , , and  are plotted in Fig. 5. This includes the data from Figs. 2–4(a) as well as six additional runs from Ref.  that are not in Fig. 2. The standard deviation of the data around Eq. (6) is , for a 95% confidence interval of . Figure 5 includes 328 data points; 326 of these (99.4%) are within of Eq. (6). The air data of Ref.  are not included because they are systematically high relative to all other studies (relative to Eq. (6), those data are 15–25% high, with a few points 30% high). Equation (8) is very close to Eq. (6) for air.
4.2 Data for Water.
Figure 6 shows Žukauskas and Šlančiauskas' data for room temperature water. In Fig. 6(a), Runs 7 and 8 have the same velocity and freestream temperature, but the temperature differences in Run 8 are significantly greater (see Table 2). For Run 7, the property-ratio correction averages 2%, and for Run 8 it averages 10%. Neither correction has been applied to the correlations plotted. Both data sets show good agreement with the laminar theory, Eq. (3), and follow a transition curve with c = 1.75 and . (The laminar correlation would be closer to the Run 8 data if the property correction were applied.)
Figure 6(b) shows data extending well into the turbulent range. The proposed correlation, Eq. (9), represents the transitional and turbulent data well. Also, shown here are the correlations of Hollingsworth  (see Table 2) and of Žukauskas and Šlančiauskas. The turbulent data for Run 18 and Hollingsworth's correlation are close one another, although about 10% below Eq. (6). The laminar data for these experiments are quite scattered. Property ratio corrections have not been applied. For Table 27, corrections average , and for the Table 25 data, corrections average .
Figure 7 shows Žukauskas and Šlančiauskas' data for , obtained using hot water in the fully turbulent regime. The expanded scale of this chart's ordinate makes the spread of data appear larger than on the other charts. Žukauskas and Šlančiauskas note that their measurement uncertainties were highest for the hot water experiments (reported as ). The data lie up to 20% below Eq. (6) at the highest Reynolds numbers (around ), and 10% below Hollingsworth's correlation. Variable properties effects on these data are negligible.
4.3 Data for Transformer Oil.
Figure 8 shows Žukauskas and Šlančiauskas' data for transformer oil at Pr of 109 and 257. For these data, property variations across the boundary layer are substantial. Therefore, the correction factor recommended by Žukauskas and Šlančiauskas, , has been applied to the proposed correlation, Eq. (9). A single average value (the average of Žukauskas and Šlančiauskas' pointwise values) has been used. With this adjustment, agreement between the data and Eq. (9) is remarkably good.
Figure 9(a) shows data for compared to the proposed correlation. For these two data sets, the variable properties correction recommended by Žukauskas and Šlančiauskas averages 12.7% and 17%. The proposed correlation is plotted both without adjustment and with a adjustment. The two curves bracket the data, in the turbulent regime and the later part of the transition regime; the data do not extend to laminar conditions.
Equation (6) is compared to Žukauskas and Šlančiauskas' data for fully turbulent flow of transformer oil at in Fig. 9(b); 13 of 16 points are within 10% of the equation (15 of 16 are within 15%). Those data do not extend to laminar and transitional conditions. The average variable properties correction for these data is and is not applied to correlation.
4.4 Recommended Values of and .
and 80% are within . The value c = 6 for Blair's low turbulence data does not lie on this curve, and more data for would be needed to extend the curve fit to higher .
Figures 4(b) and 6(a) display the intersection of the laminar correlation, Eq. (3), with the dotted line representing the transition correlation, Eq. (5). The intersection lies at . The data and the fit begin to rise above the laminar value upstream of ; and, at , the Nusselt number of Eq. (9) is 7% greater than the laminar Nusselt number.
for a laminar boundary starting at the leading edge under zero pressure gradient. However, the laminar momentum boundary layers measured mostly differ from that condition and differ between the various studies. This upstream variability together with unclear local values of leaves Mayle's equation with little predictive power for the heat transfer experiments summarized here. When experimental conditions are better defined, Mayle's equation has good experimental support .
In working scenarios with significant disturbances (), transition is likely to begin in the range . Only with extremely low levels of turbulence () can exceed 106; Schubauer and Skramstad  reported an asymptotic limit of for .
5 Average Nusselt Number
The transitional region contributes 29.2% of the total value. If transition were omitted (abruptly changing from laminar to turbulent flow, as often illustrated in textbooks), then: (a) if transition is taken to occur at (+11% error); but (b) if transition is taken to lie at (–17% error). Note also that the larger value exceeds the smaller value by 34%. The difference will be larger for smaller values of , and the difference will become smaller if increases.
New equations are proposed for calculating the Nusselt number of a flat plate boundary layer from the laminar regime, through the transition region, and into the fully turbulent region. These formulas are in good agreement with available data for air, water, and oil boundary layers spanning the ranges , with freestream turbulence levels up to 5%. The complete correlation is summarized in Table 4.
The transition region has a length similar to the laminar region, as shown in Fig. 11.
Data show that the Nusselt number in the transition region varies as where the exponent c increases with an increase in , the Reynolds number at the onset of transition. Through the examination of many independent experiments, values of c are fitted as a function of . The transition regime is described by Eqs. (5) and (11).
Churchill's multiregion correlation is modified to accurately predict the local Nusselt number in the transition and turbulent ranges. The new correlation, Eq. (9), shows good agreement with a large number of experimental data drawn from many independent studies.
Equation (6) correlates the fully turbulent air data to (two standard deviations); 99.4% of the data are within .
The fully turbulent air data of Žukauskas and Šlančiauskas are systematically higher than air data from other investigators. Similarly, their correlation, Eq. (10), is higher for air than other correlations.
The data of Fage and Falkner, for UHF laminar air flow, are in good agreement with corresponding analytical results, contrary to the suggestion of Jakob and Dow, down to an overall Reynolds number of 1400 (Appendix B).
Very few measurements are available for heat transfer in turbulent liquid boundary layers.
Older literature on convection contains a great deal of information that remains useful today.
In modern notation, Pohlhausen's equation is .
Sugawara et al. measured the transient heating and cooling of a plate in a wind tunnel at a dozen positions along the plate. The plate temperature response, locally, was treated as lumped. The plate would cool faster near the leading edge.
Reynolds et al. suggest multiplying by variable-properties correction factor of , which amounted to about 2% for their data. Their equation fits their data to a standard deviation of .
The table numbers for Žukauskas and Šlančiauskas's data refer to the data tables in the appendix of their book .
I am grateful for useful suggestions from Leon R. Glicksman (MIT), Je-Chin Han (Texas A&M University), D. Keith Hollingsworth (University of Alabama, Huntsville), and John H. Lienhard, IV (University of Houston).
- a =
coefficient in Eq. (3)
- b =
coefficient in Eq. (4)
- c =
exponent on Reynolds number, Eq. (4)
- Cf =
local skin friction coefficient
- d =
coefficient in Churchill's function
- f =
Darcy friction factor
- h =
local heat transfer coefficient, W m–2 K–1
average heat transfer coefficient, W m–2 K–1
integral defined by Eq. (B7)
- k =
thermal conductivity of fluid, W m–1 K–1
local Nusselt number, hx/k
average Nusselt number,
- p =
- Pr =
- qw =
wall heat flux, W m–2
- q0 =
a uniform heat flux, W m–2
local Reynolds number,
overall Reynolds number,
- s =
- St =
local Stanton number, Nux/Rex Pr
- Tw =
wall temperature, K
freestream temperature, K
- u0 =
mean freestream speed, m s–1
r.m.s. fluctuation of freestream speed, m s–1
- x =
distance from leading edge, m
- x0 =
unheated starting length, m
- ν =
kinematic viscosity of fluid, m2 s–1
- σ =
Churchill's function, Eq. (A2)
Churchill's Correlation in Full
and is a number typically between about 105 and 107. The actual value of must be fit to each specific dataset. If the Reynolds number at the end of the transition region is , an estimate is .
The equation is for uniform Tw. To adapt it to uniform qw, the constants 0.3387, 0.0468, and 2600 are replaced by 0.4637, 0.02052, and 7420, respectively.5 Churchill gave a similar equation for the average Nusselt number.
Laminar Boundary Layer With Uniform Heat Flux
The laminar similarity solution for uniform heat flux is not well known: even recently, some well-respected books report that there is no such solution . This appendix provides that background, as well as sketching the UHF solution with an unheated starting length, which is less well known than that for UWT.
B.1 Uniform Flux Over Entire Plate.
Imai reported 1% agreement with Levy's numerical values down to , which is within the accuracy of Levy's numerics.
Further, Fage and Falkner's experimental values for air fit well to , based on a value of that is 12% lower than modern property data. The temperature ratio, , is high enough that variable properties effects are not negligible. Without a variable properties correction, the experiments are within of Imai's formula. If a correction of is applied, the agreement is –6%.
B.2 Uniform Flux With Unheated Starting Length.
What is key here is that the approximation is bounded within 12% right up to x0.
This result is within 1% of Imai's formula.
Comments on the Classical Colburn Analogy
By combining this expression with an equation for Cf as power of the Reynolds number, the Nusselt number can be expressed as a product of powers of the Reynolds and Prandtl numbers. Such equations have commonly been used for turbulent boundary layers and turbulent pipe flows.
Unfortunately, these power laws cannot accurately span a wide range of Prandtl number: at high Reynolds numbers, the dependence of the Reynolds number exponent on the Prandtl number becomes significant so that a single, fixed exponent for Re and for Pr is inadequate to cover a wide range of conditions. Gnielinski  and Žukauskas and Šlančiauskas  discuss this point with reference to experimental data for either pipe or boundary layer flows. Transport models that take account of the structure of the turbulent boundary layer result in expressions having the general form of Eq. (2) [24,26,28,29]. Those expressions are capable of covering the full range of Prandtl number [29,30]. For a narrow range of Pr (say, just for gases or for some range of liquids), power-law formulas can be quite accurate, as shown by Gnielinski.
Colburn considered a wide spectrum of data in his work, but neither Colburn  nor Chilton and Colburn  had much data for external boundary layers. And the data they had were only for gases. Chilton and Colburn included three data points for evaporation of water into laminar air flow. In the turbulent regime, Colburn referred to two somewhat scattered sets of air data. In the laminar rage, Colburn used the UHF air data of Fage and Falkner , showing good agreement with Eq. (35); but he later confessed to Jakob and Dow that he had made a “slide-rule error” in plotting these data  and that the points should have been 21% higher. In fact, the Colburn analogy is incapable of discriminating the laminar thermal boundary condition.
Typical Values of Prandtl Number
Values of the Prandtl number for the fluids discussed in this paper are given in Table 5.