Heat Transfer in Fully-Developed Internal Turbulent Flow

From Thermal-FluidsPedia

(Difference between revisions)
Jump to: navigation, search
Line 338: Line 338:
The temperature at the center of the tube, <math>{{\bar{T}}_{c}}</math>, can be obtained by letting  
The temperature at the center of the tube, <math>{{\bar{T}}_{c}}</math>, can be obtained by letting <math>{{y}^{+}}=y_{c}^{+}</math> in eq. (5.329), i.e.
in eq. (5.329), i.e.
{| class="wikitable" border="0"
{| class="wikitable" border="0"

Revision as of 03:08, 8 July 2010

Heat transfer in fully-developed turbulent flow in a circular tube subject to constant heat flux (q''w = const) will be considered in this subsection (Oosthuizen and Naylor, 1999). When the turbulent flow in the tube is fully developed, we have \bar{v}=0 and the energy eq. (5.268) becomes

\bar{u}\frac{\partial \bar{T}}{\partial x}=\frac{1}{r}\frac{\partial }{\partial r}\left[ r(\alpha +{{\varepsilon }_{H}})\frac{\partial \bar{T}}{\partial r} \right]


After the turbulent flow is hydrodynamically and thermally fully developed, the time-averaged temperature profile is no longer a function of axial distance from the inlet, i.e.,

\frac{\partial }{\partial x}\left( \frac{{{T}_{w}}-\bar{T}}{{{T}_{w}}-{{{\bar{T}}}_{c}}} \right)=0


where {{\bar{T}}_{c}} is the time-averaged temperature at the centerline of the tube, and Tw is the wall temperature. Thus, ({{T}_{w}}-\bar{T})/({{T}_{w}}-{{\bar{T}}_{c}}) is a function of r only, i.e.,

\frac{{{T}_{w}}-\bar{T}}{{{T}_{w}}-{{{\bar{T}}}_{c}}}=f\left( r \right)


where f is independent from x. Differentiating (5.295) yields

\frac{\partial \bar{T}}{\partial x}=\frac{d{{T}_{w}}}{dx}-\left( \frac{{{T}_{w}}-\bar{T}}{{{T}_{w}}-{{{\bar{T}}}_{c}}} \right)\left( \frac{d{{T}_{w}}}{dx}-\frac{d{{{\bar{T}}}_{c}}}{dx} \right)


At the wall, the contribution of eddy diffusivity on the heat transfer is negligible, and the heat flux at the wall becomes

{{{q}''}_{w}}=k{{\left. \frac{\partial \bar{T}}{\partial r} \right|}_{r={{r}_{0}}}}


Substituting eq. (5.296) into eq. (5.298), one obtains:



Since the heat flux is constant, q''w = const, it follows that ({{T}_{w}}-{{\bar{T}}_{c}})= Const, i.e.,



Therefore, eq. (5.297) becomes:

\frac{\partial \bar{T}}{\partial x}=\frac{d{{T}_{w}}}{dx}


For fully developed flow, the local heat transfer coefficient is:



where {{\bar{T}}_{m}} is the time-averaged mean temperature defined as:



Since q''w = const, it follows from eq. (5.302) that ({{T}_{w}}-{{\bar{T}}_{m}})= Const, i.e.,



Combining eqs. (5.300), (5.301) and (5.304), the following relationships are obtained:

\frac{\partial \bar{T}}{\partial x}=\frac{d{{T}_{w}}}{dx}=\frac{d{{{\bar{T}}}_{c}}}{dx}=\frac{d{{{\bar{T}}}_{m}}}{dx}


The time-averaged mean temperature, {{\bar{T}}_{m}} , changes with x as the result of heat transfer from the tube wall. By following the same procedure as that in Example 5.2, the rate of mean temperature change can be obtained as follows:

\frac{d{{{\bar{T}}}_{m}}}{dx}=\frac{4\pi {{{{q}''}}_{w}}}{\rho {{c}_{p}}D{{{\bar{u}}}_{m}}}


Substituting eq. (5.305) into eq. (5.294), the energy equation becomes:

\bar{u}\frac{d{{{\bar{T}}}_{m}}}{dx}=\frac{1}{{{r}_{0}}-y}\frac{\partial }{\partial y}\left[ ({{r}_{0}}-y)(\alpha +{{\varepsilon }_{H}})\frac{\partial \bar{T}}{\partial y} \right]


where y = r0r is the distance measured from the tube wall. Equation (5.307) is subject to the following two boundary conditions:

\frac{\partial \bar{T}}{\partial y}=0,\text{  }y={{r}_{0}} (axisymmetric condition)


\bar{T}={{T}_{w}}\text{ (unknown)},\text{ }y=0


Integrating eq. (5.307) in the interval of (r0, r) and considering eq. (5.308), we have:

({{r}_{0}}-y)(\alpha +{{\varepsilon }_{H}})\frac{\partial \bar{T}}{\partial y}=\frac{d{{{\bar{T}}}_{m}}}{dx}\int_{{{r}_{0}}}^{y}{({{r}_{0}}-y)\bar{u}dy}


which can be rearranged to

\frac{\partial \bar{T}}{\partial y}=\frac{I(y)}{({{r}_{0}}-y)(\alpha +{{\varepsilon }_{H}})}\frac{d{{{\bar{T}}}_{m}}}{dx}





Integrating eq. (5.311) in the interval of (0, y) and considering eq. (5.309), one obtains:

\bar{T}-{{T}_{w}}=\frac{d{{{\bar{T}}}_{m}}}{dx}\int_{0}^{y}{\frac{I(y)}{({{r}_{o}}-y)(\alpha +{{\varepsilon }_{H}})}dy}


If the profiles of axial velocity and the thermal eddy diffusivity are known, eq. (5.313) can be used to obtain the correlation for internal forced convection heat transfer. With the exception of the very thin viscous sublayer, the velocity profile in the most part of the tube is fairly flat. Therefore, it is assumed that the time-averaged velocity, \bar{u}, in eq. (5.312) can be replaced by {{\bar{u}}_{m}}, and I(y) becomes:

I(y)\simeq -\frac{\bar{u}_{m}^{2}}{2}{{({{r}_{o}}-y)}^{2}}


Substituting eqs. (5.314) and (5.306) into eq. (5.313) yields:

\bar{T}-{{T}_{w}}=-\frac{{{{{q}''}}_{w}}}{\rho {{c}_{p}}}\int_{0}^{y}{\frac{(1-y/{{r}_{o}})}{(\alpha +{{\varepsilon }_{H}})}dy}


which can be rewritten in terms of wall coordinate

\bar{T}-{{T}_{w}}=-\frac{{{{{q}''}}_{w}}}{\rho {{c}_{p}}}\sqrt{\frac{\rho }{{{\tau }_{w}}}}\int_{0}^{{{y}^{+}}}{\frac{(1-y/{{r}_{o}})}{[1/\Pr +({{\varepsilon }_{M}}/\nu )/{{\Pr }^{t}}]}d{{y}^{+}}}


where y+ is defined in eq. (5.277). To consider heat transfer in an internal turbulent flow, the entire turbulent boundary layer is divided into three regions: (1) inner region (y + < 5), (2) buffer region (5\le {{y}^{+}}\le 30), and (3) outer region (y + > 30). In the inner region {{\varepsilon }_{M}}={{\varepsilon }_{H}}=0 and eq. (5.316) becomes

\bar{T}-{{T}_{w}}=-\frac{{{{{q}''}}_{w}}}{\rho {{c}_{p}}}\sqrt{\frac{\rho }{{{\tau }_{w}}}}\Pr \int_{0}^{{{y}^{+}}}{(1-y/{{r}_{o}})d{{y}^{+}}}


Since the inner region is very thin, y/{{r}_{0}}\ll 1 and 1 − y / r0 is effectively equal to 1. Therefore, the temperature profile in the inner region becomes:

\bar{T}-{{T}_{w}}=-\left( \frac{{{{{q}''}}_{w}}}{\rho {{c}_{p}}} \right)\sqrt{\frac{\rho }{{{\tau }_{w}}}}\Pr {{y}^{+}}


The temperature at the boundary between the inner and buffer regions (y + = 5), {{\bar{T}}_{s}}, can be obtained from eq. (5.318) as

{{\bar{T}}_{s}}-{{T}_{w}}=-5\left( \frac{{{{{q}''}}_{w}}}{\rho {{c}_{p}}} \right)\sqrt{\frac{\rho }{{{\tau }_{w}}}}\Pr


In the buffer region where 5\le {{y}^{+}}\le 30, the eddy diffusivity in the buffer region is:

\frac{{{\varepsilon }_{M}}}{\nu }=\frac{{{y}^{+}}}{5}-1


Substituting eq. (5.320) into eq. (5.316) and assuming the turbulent Prandtl number {{\Pr }^{t}}=1, the following expression is obtained:

\bar{T}-{{\bar{T}}_{s}}=-\frac{{{{{q}''}}_{w}}}{\rho {{c}_{p}}}\sqrt{\frac{\rho }{{{\tau }_{w}}}}\int_{5}^{{{y}^{+}}}{\frac{(1-y/{{r}_{o}})}{[1/\Pr +{{y}^{+}}/5-1]}d{{y}^{+}}}


Since the buffer region is also very thin, 1 − y / r0 in eq. (5.321) is effectively equal to 1. Defining {{T}^{+}}={(\bar{T}-{{T}_{w}})}/{\left( -\frac{{{{{q}''}}_{w}}}{\rho {{c}_{p}}}\sqrt{\frac{\rho }{{{\tau }_{w}}}} \right)}\; and Integrating eq. (5.321) yields

\int_{5\Pr }^{{{T}^{+}}}{d{{T}^{+}}}=\int_{5}^{{{y}^{+}}}{\frac{d{{y}^{+}}}{1/\Pr +({{y}^{+}}-5)/(5{{\Pr }^{t}})}}



\bar{T}-{{\bar{T}}_{s}}=-5\left( \frac{{{{{q}''}}_{w}}}{\rho {{c}_{p}}} \right)\sqrt{\frac{\rho }{{{\tau }_{w}}}}\ln \left( \frac{{{y}^{+}}}{5}\Pr -\Pr +1 \right),\text{   }5<{{y}^{+}}<30


The temperature at the top of the buffer region where y + = 30, {{\bar{T}}_{b}} , becomes

{{\bar{T}}_{b}}-{{\bar{T}}_{s}}=-5\left( \frac{{{{{q}''}}_{w}}}{\rho {{c}_{p}}} \right)\sqrt{\frac{\rho }{{{\tau }_{w}}}}\ln (5\Pr +1)


For the outer region where {{\varepsilon }_{M}}\gg \nu \text{ and }{{\varepsilon }_{H}}\gg \alpha , eq. (5.316) becomes

\bar{T}-{{\bar{T}}_{b}}=-\frac{{{{{q}''}}_{w}}}{\rho {{c}_{p}}}\sqrt{\frac{\rho }{{{\tau }_{w}}}}\int_{30}^{{{y}^{+}}}{\frac{(1-y/{{r}_{o}})}{{{\varepsilon }_{M}}/\nu }d{{y}^{+}}}


where the turbulent Prandtl number is assumed to be equal to 1. It is assumed that the Nikuradse equation (5.276) is valid in the outer region and the velocity gradient in this region becomes:

\frac{\partial {{u}^{+}}}{\partial {{y}^{+}}}=\frac{2.5}{{{y}^{+}}}


The expression of apparent shear stress in this region, eq. (5.280) , can be non-dimensionalized using eqs. (5.277) and (5.278) as:

\frac{{{\tau }_{app}}}{{{\tau }_{w}}}=\frac{{{\varepsilon }_{M}}}{\nu }\frac{\partial {{u}^{+}}}{\partial {{y}^{+}}}


Substituting eqs. (5.275) and (5.326) into eq. (5.327), the eddy diffusivity in the outer region is obtained as:

\frac{{{\varepsilon }_{M}}}{\nu }=\left( 1-\frac{y}{{{r}_{o}}} \right)\frac{{{y}^{+}}}{2.5}


Substituting eq. (5.328) into eq. (5.325), the temperature distribution in this region becomes:

\bar{T}-{{\bar{T}}_{b}}=-2.5\left( \frac{{{{{q}''}}_{w}}}{\rho {{c}_{p}}} \right)\sqrt{\frac{\rho }{{{\tau }_{w}}}}\int_{30}^{{{y}^{+}}}{\frac{1}{{{y}^{+}}}d{{y}^{+}}}=-2.5\left( \frac{{{{{q}''}}_{w}}}{\rho {{c}_{p}}} \right)\sqrt{\frac{\rho }{{{\tau }_{w}}}}\ln \left( \frac{{{y}^{+}}}{30} \right)


which is valid from y + = 30 to the center of the tube where yc = r0 or

y_{c}^{+}=\frac{{{r}_{0}}}{\nu }\sqrt{\frac{{{\tau }_{w}}}{\rho }}


The temperature at the center of the tube, {{\bar{T}}_{c}}, can be obtained by letting {{y}^{+}}=y_{c}^{+} in eq. (5.329), i.e.

{{\bar{T}}_{c}}-{{\bar{T}}_{b}}=-2.5\left( \frac{{{{{q}''}}_{w}}}{\rho {{c}_{p}}} \right)\sqrt{\frac{\rho }{{{\tau }_{w}}}}\ln \left( \frac{{{r}_{0}}}{30\nu }\sqrt{\frac{{{\tau }_{w}}}{\rho }} \right)


The overall temperature change from the wall to the center of the tube can be obtained by adding eqs. (5.319), (5.324) and (5.331):

{{T}_{w}}-{{\bar{T}}_{c}}=\left( \frac{{{{{q}''}}_{w}}}{\rho {{c}_{p}}} \right)\sqrt{\frac{\rho }{{{\tau }_{w}}}}\left[ 2.5\ln \left( \frac{{{r}_{0}}}{30\nu }\sqrt{\frac{{{\tau }_{w}}}{\rho }} \right)+5\ln (5\Pr +1)+5\Pr  \right]


It follows from the definition of friction factor, eq. (5.282), that

{{\tau }_{w}}=\frac{1}{2}{{c}_{f}}\rho \bar{u}_{m}^{2}


Substituting eq. (5.333) into eq. (5.332) and considering the definition of Reynolds number, {{\operatorname{Re}}_{D}}={{\bar{u}}_{m}}D/\nu , eq. (5.332) becomes:

{{T}_{w}}-{{\bar{T}}_{c}}=\left( \frac{{{{{q}''}}_{w}}}{\rho {{c}_{p}}{{{\bar{u}}}_{m}}} \right)\sqrt{\frac{2}{{{c}_{f}}}}\left[ 2.5\ln \left( \frac{{{\operatorname{Re}}_{D}}}{60}\sqrt{\frac{{{c}_{f}}}{2}} \right)+5\ln (5\Pr +1)+5\Pr  \right] (5.334) In order to obtain the heat transfer coefficient, h={{{q}''}_{w}}/({{T}_{w}}-{{\bar{T}}_{m}}), the temperature difference {{T}_{w}}-{{\bar{T}}_{m}} must be obtained. If the velocity profile can be approximated by eq. (5.287), and the temperature and velocity can also be approximated by the one-seventh law, i.e.,

\frac{{{T}_{w}}-\bar{T}}{{{T}_{w}}-{{{\bar{T}}}_{c}}}={{\left( \frac{y}{{{r}_{o}}} \right)}^{1/7}}{{,}_{_{\text{        }}}}\frac{{\bar{u}}}{{{{\bar{u}}}_{c}}}={{\left( \frac{y}{{{r}_{o}}} \right)}^{1/7}}


it follows that

{{T}_{w}}-{{\bar{T}}_{m}}=\frac{\int_{0}^{{{r}_{o}}}{\bar{u}({{T}_{w}}-\bar{T})2\pi rdr}}{\int_{0}^{{{r}_{o}}}{\bar{u}2\pi rdr}}=\frac{5}{6}({{T}_{w}}-{{\bar{T}}_{c}})


Substituting eq. (5.334) into eq. (5.336) results in:

{{T}_{w}}-{{\bar{T}}_{m}}=\frac{5}{6}\left( \frac{{{{{q}''}}_{w}}}{\rho {{c}_{p}}{{{\bar{u}}}_{m}}} \right)\sqrt{\frac{2}{{{c}_{f}}}}\left[ 2.5\ln \left( \frac{{{\operatorname{Re}}_{D}}}{60}\sqrt{\frac{{{c}_{f}}}{2}} \right)+5\ln (5\Pr +1)+5\Pr  \right]


which can be rearranged to the following empirical correlation

\text{N}{{\text{u}}_{D}}=\frac{{{\operatorname{Re}}_{D}}\Pr \sqrt{\frac{{{c}_{f}}}{2}}}{\frac{5}{6}\left[ 2.5\ln \left( \frac{{{\operatorname{Re}}_{D}}}{60}\sqrt{\frac{{{c}_{f}}}{2}} \right)+5\ln (5\Pr +1)+5\Pr  \right]}


which can be used together with appropriate friction coefficient discussed in the previous subsection to obtain the Nusselt number.