Scale analysis of natural convection

From Thermal-FluidsPedia

(Difference between revisions)
Jump to: navigation, search
Line 1: Line 1:
-
While scale analysis cannot provide the exact functions in eq. (6.36), the form of these functions can be provided [[#References|(Bejan, 2004)]]. Let us refer to Fig. 6.1 and consider the governing equations in the thermal boundary layer (<math>y\sim {{\delta }_{t}}</math>) for the entire flat plate (<math>x\sim L</math>). The thickness of the thermal boundary layer in which the effects of the heated wall are felt is much smaller than the length of the vertical plate, i.e. <math>{{\delta }_{t}}\ll L</math>. For the continuity equation (6.17) to be satisfied, the scales of the two terms must be the same:<math>\frac{u}{L}\sim \frac{v}{{{\delta }_{t}}}</math>
+
[[Image:Natural convection over a vertical flat plate.jpg|thumb|400 px|alt=Natural convection over a vertical flat plate (Pr >1) |Figure 6.1 Natural convection over a vertical flat plate (Pr >1).]]
 +
 
 +
While scale analysis cannot provide the exact functions in eq. (7), the form of these functions can be provided [[#References|(Bejan, 2004)]]. Let us refer to Fig. 6.1 and consider the governing equations in the thermal boundary layer (<math>y\sim {{\delta }_{t}}</math>) for the entire flat plate (<math>x\sim L</math>). The thickness of the thermal boundary layer in which the effects of the heated wall are felt is much smaller than the length of the vertical plate, i.e. <math>{{\delta }_{t}}\ll L</math>. For the continuity equation (8) to be satisfied, the scales of the two terms must be the same:<math>\frac{u}{L}\sim \frac{v}{{{\delta }_{t}}}</math>
The scale of the velocity component in the y-direction is therefore:
The scale of the velocity component in the y-direction is therefore:
Line 13: Line 15:
which indicates that <math>v\ll u</math> for flow in the thermal boundary layer. The scale of the velocity component in the x-direction, ''u'', is still unknown at this point.
which indicates that <math>v\ll u</math> for flow in the thermal boundary layer. The scale of the velocity component in the x-direction, ''u'', is still unknown at this point.
-
While the left hand side of the energy equation (6.21) shows the effect of advection, the right-hand side shows the effect of diffusion. The scales of the two advective terms on the left hand side of eq. (6.21) are:<math>u\frac{\partial T}{\partial x}\sim u\frac{\Delta T}{L},\text{ }v\frac{\partial T}{\partial y}\sim v\frac{\Delta T}{{{\delta }_{t}}}\sim u\frac{\Delta T}{L}</math>
+
While the left hand side of the energy equation (9) shows the effect of advection, the right-hand side shows the effect of diffusion. The scales of the two advective terms on the left hand side of eq. (9) are:<math>u\frac{\partial T}{\partial x}\sim u\frac{\Delta T}{L},\text{ }v\frac{\partial T}{\partial y}\sim v\frac{\Delta T}{{{\delta }_{t}}}\sim u\frac{\Delta T}{L}</math>
 
 
-
which indicates that the scale of the second term on the left hand side of eq. (6.21) is identical to the scale of the first term when the scale of the velocity component in the y-direction is given by eq. (6.37). The temperature difference <math>\Delta T={{T}_{w}}-{{T}_{\infty }}</math> in the above scale analysis represents the scale of the excess temperature, <math>T-{{T}_{\infty }}</math>. The scale of the right-hand side of eq. (6.21) is: <math>\alpha \frac{{{\partial }^{2}}T}{\partial {{y}^{2}}}\sim \alpha \frac{\Delta T}{\delta _{t}^{2}}</math>
+
which indicates that the scale of the second term on the left hand side of eq. (9) is identical to the scale of the first term when the scale of the velocity component in the y-direction is given by eq. (1). The temperature difference <math>\Delta T={{T}_{w}}-{{T}_{\infty }}</math> in the above scale analysis represents the scale of the excess temperature, <math>T-{{T}_{\infty }}</math>. The scale of the right-hand side of eq. (9) is: <math>\alpha \frac{{{\partial }^{2}}T}{\partial {{y}^{2}}}\sim \alpha \frac{\Delta T}{\delta _{t}^{2}}</math>
 
 
The scales of the two sides of the energy equation must be the same:<math>u\frac{\Delta T}{L}\sim \alpha \frac{\Delta T}{\delta _{t}^{2}}</math>
The scales of the two sides of the energy equation must be the same:<math>u\frac{\Delta T}{L}\sim \alpha \frac{\Delta T}{\delta _{t}^{2}}</math>
Line 41: Line 43:
where the scale of the thermal boundary layer thickness, δ<sub>t</sub>, is still unknown at this point.
where the scale of the thermal boundary layer thickness, δ<sub>t</sub>, is still unknown at this point.
-
The respective scales of the two inertial terms on the left-hand side of the momentum equation (6.20) are:
+
The respective scales of the two inertial terms on the left-hand side of the momentum equation (10) are:
 
 
<math>u\frac{\partial u}{\partial x}\sim \frac{{{u}^{2}}}{L}</math>  
<math>u\frac{\partial u}{\partial x}\sim \frac{{{u}^{2}}}{L}</math>  
Line 47: Line 49:
<math>v\frac{\partial u}{\partial y}\sim v\frac{u}{{{\delta }_{t}}}\sim \frac{{{u}^{2}}}{L}</math>
<math>v\frac{\partial u}{\partial y}\sim v\frac{u}{{{\delta }_{t}}}\sim \frac{{{u}^{2}}}{L}</math>
 
 
-
The respective scales of the viscosity and buoyancy terms on the right-hand side of eq. (6.20) are:
+
The respective scales of the viscosity and buoyancy terms on the right-hand side of eq. (10) are:
 
 
<math>\nu \frac{{{\partial }^{2}}u}{\partial {{y}^{2}}}\sim \nu \frac{u}{\delta _{t}^{2}}</math>
<math>\nu \frac{{{\partial }^{2}}u}{\partial {{y}^{2}}}\sim \nu \frac{u}{\delta _{t}^{2}}</math>
Line 96: Line 98:
|}
|}
-
is the Rayleigh number that is related to the Grashof number by  
+
is the Rayleigh number that is related to the Grashof number by <math>\text{R}{{\text{a}}_{L}}=\text{G}{{\text{r}}_{L}}\Pr </math>. Therefore, the relative importance of inertia and viscous forces depends on the Prandtl number, which is a property of the fluid. Thus, if the Prandtl number is high (<math>\Pr \gg 1</math>), the inertia term will be negligible and the viscosity term will balance the buoyancy term, whereas if the Prandtl number is low enough (<math>\Pr \ll 1</math>) – as for liquid metals – then the inertia term is considerable and balances the buoyancy term in steady state. The scale analysis for high- and low-Prandtl number fluids is presented in detail below.
-
<math>\text{R}{{\text{a}}_{L}}=\text{G}{{\text{r}}_{L}}\Pr </math>
+
 
-
. Therefore, the relative importance of inertia and viscous forces depends on the Prandtl number, which is a property of the fluid. Thus, if the Prandtl number is high (<math>\Pr \gg 1</math>), the inertia term will be negligible and the viscosity term will balance the buoyancy term, whereas if the Prandtl number is low enough (<math>\Pr \ll 1</math>) – as for liquid metals – then the inertia term is considerable and balances the buoyancy term in steady state. The scale analysis for high- and low-Prandtl number fluids is presented in detail below.
+
==Additional Equations==
 +
 
 +
{| class="wikitable" border="0"
 +
|-
 +
| width="100%" |<center>
 +
<math>{{\overline{\text{Nu}}}_{L}}=f(\text{G}{{\text{r}}_{L}},\Pr )</math>
 +
</center>
 +
|{{EquationRef|(7)}}
 +
|}
 +
 
 +
{| class="wikitable" border="0"
 +
|-
 +
| width="100%" |<center>
 +
<math>\frac{\partial u}{\partial x}+\frac{\partial v}{\partial y}=0</math>
 +
</center>
 +
|{{EquationRef|(8)}}
 +
|}
 +
 
 +
{| class="wikitable" border="0"
 +
|-
 +
| width="100%" |<center>
 +
<math>u\frac{\partial T}{\partial x}+v\frac{\partial T}{\partial y}=\alpha \frac{{{\partial }^{2}}T}{\partial {{y}^{2}}}</math>
 +
</center>
 +
|{{EquationRef|(9)}}
 +
|}
 +
 
 +
{| class="wikitable" border="0"
 +
|-
 +
| width="100%" |<center>
 +
<math>u\frac{\partial u}{\partial x}+v\frac{\partial u}{\partial y}=\nu \frac{{{\partial }^{2}}u}{\partial {{y}^{2}}}+g\beta (T-{{T}_{\infty }})</math>
 +
</center>
 +
|{{EquationRef|(10)}}
 +
|}

Revision as of 14:40, 7 July 2010

Natural convection over a vertical flat plate (Pr >1)
Figure 6.1 Natural convection over a vertical flat plate (Pr >1).

While scale analysis cannot provide the exact functions in eq. (7), the form of these functions can be provided (Bejan, 2004). Let us refer to Fig. 6.1 and consider the governing equations in the thermal boundary layer (y˜δt) for the entire flat plate (x˜L). The thickness of the thermal boundary layer in which the effects of the heated wall are felt is much smaller than the length of the vertical plate, i.e. {{\delta }_{t}}\ll L. For the continuity equation (8) to be satisfied, the scales of the two terms must be the same:\frac{u}{L}\sim \frac{v}{{{\delta }_{t}}}

The scale of the velocity component in the y-direction is therefore:

v\sim \frac{{{\delta }_{t}}}{L}u

(1)

which indicates that v\ll u for flow in the thermal boundary layer. The scale of the velocity component in the x-direction, u, is still unknown at this point.

While the left hand side of the energy equation (9) shows the effect of advection, the right-hand side shows the effect of diffusion. The scales of the two advective terms on the left hand side of eq. (9) are:u\frac{\partial T}{\partial x}\sim u\frac{\Delta T}{L},\text{ }v\frac{\partial T}{\partial y}\sim v\frac{\Delta T}{{{\delta }_{t}}}\sim u\frac{\Delta T}{L}

which indicates that the scale of the second term on the left hand side of eq. (9) is identical to the scale of the first term when the scale of the velocity component in the y-direction is given by eq. (1). The temperature difference \Delta T={{T}_{w}}-{{T}_{\infty }} in the above scale analysis represents the scale of the excess temperature, T-{{T}_{\infty }}. The scale of the right-hand side of eq. (9) is: \alpha \frac{{{\partial }^{2}}T}{\partial {{y}^{2}}}\sim \alpha \frac{\Delta T}{\delta _{t}^{2}}

The scales of the two sides of the energy equation must be the same:u\frac{\Delta T}{L}\sim \alpha \frac{\Delta T}{\delta _{t}^{2}}

The above equation can be used to estimate the scale of u as follows:

u\sim \alpha \frac{L}{\delta _{t}^{2}}

(2)

The scale of v can be obtained by substituting eq. (2) into eq. (1):

v\sim \frac{\alpha }{{{\delta }_{t}}}

(3)

where the scale of the thermal boundary layer thickness, δt, is still unknown at this point.

The respective scales of the two inertial terms on the left-hand side of the momentum equation (10) are:

u\frac{\partial u}{\partial x}\sim \frac{{{u}^{2}}}{L}

v\frac{\partial u}{\partial y}\sim v\frac{u}{{{\delta }_{t}}}\sim \frac{{{u}^{2}}}{L}

The respective scales of the viscosity and buoyancy terms on the right-hand side of eq. (10) are:

\nu \frac{{{\partial }^{2}}u}{\partial {{y}^{2}}}\sim \nu \frac{u}{\delta _{t}^{2}}


g\beta (T-{{T}_{\infty }})\sim g\beta \Delta T

We can very well see that three forces are at play in the boundary layer region, which are inertia, viscosity and buoyancy forces. Considering the scale of u obtained from eq. (2), the scales of these three forces are

\begin{matrix}
   \frac{{{\alpha }^{2}}L}{\delta _{t}^{4}}, & \frac{\nu \alpha L}{\delta _{t}^{4}}, & g\beta \Delta T  \\
   \text{Inertia} & \text{Viscous} & \text{Buoyancy}  \\
\end{matrix}

(4)

Among these three forces, the buoyancy force is never negligible because without it natural convection would not occur. Therefore, the scale of buoyancy force can be used to measure the importance of the inertial and viscous forces. Dividing eq. (4) by the scale of buoyancy force, gβΔT , one obtains the following:\begin{matrix}
   \frac{{{\alpha }^{2}}}{g\beta \Delta T{{L}^{3}}}{{\left( \frac{L}{{{\delta }_{t}}} \right)}^{4}}, & \frac{\nu \alpha }{g\beta \Delta T{{L}^{3}}}{{\left( \frac{L}{{{\delta }_{t}}} \right)}^{4}}, & 1  \\
   \text{Inertia} & \text{Viscous} & \text{Buoyancy}  \\
\end{matrix}

which can be expressed in terms of dimensionless parameters as (Bejan, 2004):

\begin{matrix}
   {{\left( \frac{L}{{{\delta }_{t}}} \right)}^{4}}\text{Ra}_{L}^{-1}{{\Pr }^{-1}}, & {{\left( \frac{L}{{{\delta }_{t}}} \right)}^{4}}\text{Ra}_{L}^{-1}, & 1  \\
   \text{Inertia} & \text{Viscous} & \text{Buoyancy}  \\
\end{matrix}

(5)

where

\text{R}{{\text{a}}_{L}}=\frac{g\beta ({{T}_{w}}-{{T}_{\infty }}){{L}^{3}}}{\nu \alpha }

(6)

is the Rayleigh number that is related to the Grashof number by \text{R}{{\text{a}}_{L}}=\text{G}{{\text{r}}_{L}}\Pr . Therefore, the relative importance of inertia and viscous forces depends on the Prandtl number, which is a property of the fluid. Thus, if the Prandtl number is high (\Pr \gg 1), the inertia term will be negligible and the viscosity term will balance the buoyancy term, whereas if the Prandtl number is low enough (\Pr \ll 1) – as for liquid metals – then the inertia term is considerable and balances the buoyancy term in steady state. The scale analysis for high- and low-Prandtl number fluids is presented in detail below.

Additional Equations

{{\overline{\text{Nu}}}_{L}}=f(\text{G}{{\text{r}}_{L}},\Pr )

(7)

\frac{\partial u}{\partial x}+\frac{\partial v}{\partial y}=0

(8)

u\frac{\partial T}{\partial x}+v\frac{\partial T}{\partial y}=\alpha \frac{{{\partial }^{2}}T}{\partial {{y}^{2}}}

(9)

u\frac{\partial u}{\partial x}+v\frac{\partial u}{\partial y}=\nu \frac{{{\partial }^{2}}u}{\partial {{y}^{2}}}+g\beta (T-{{T}_{\infty }})

(10)