Free boundary flow
From ThermalFluidsPedia
Yuwen Zhang (Talk  contribs) 

Line 1:  Line 1:  
  [[Image:Examples of free boundary flows.jpgthumb400 pxalt=Figure  +  [[Image:Examples of free boundary flows.jpgthumb400 pxalt=Figure 1 Examples of free boundary flows Figure 1 Examples of free boundary flows.]] 
  [[Image:Physical model of natural convection from a line heat source.jpgthumb400 pxalt=Figure  +  [[Image:Physical model of natural convection from a line heat source.jpgthumb400 pxalt=Figure 2 Physical model of natural convection from a line heat source Figure 2 Physical model of natural convection from a line heat source.]] 
  +  The common feature of external convection discussed in the preceding subsections is that the flows are always near a heated (or cooled) solid wall. For some other applications, however, the thermally induced flow occurs without presence of a solid wall and is referred to as free boundary flow. When a point or line heat source is immersed into a bulk fluid, the fluid near the heat source is heated, becomes lighter and rises to form a plume [see Fig. 1 (a)]. The plume generated by a point heat source is axisymmetric, but a line heat source will create a twodimensional plume. A thermal is a column of rising air near the ground due to the uneven solar heating of the ground surface. The lighter air near the ground rises and cools due to expansion [see Fig. 1(b)].  
  +  Let us consider now the twodimensional free boundary flow induced by a line heat source (see Fig. 2). Since the problem is symmetric about x = 0, only half of the domain (<math>x>0</math>) needs to be studied. In the coordinate system shown in Fig. 2, the continuity, momentum and energy equations are the same as those for natural convection near a [[Governing equations for natural convection on a vertical platevertical flat plate]]:  
  +  <center>  
+  <math>\frac{\partial u}{\partial x}+\frac{\partial v}{\partial y}=0</math>  
+  </center>  
+  
+  <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>  
+  <center>  
+  <math>u\frac{\partial T}{\partial x}+v\frac{\partial T}{\partial y}=\alpha \frac{{{\partial }^{2}}T}{\partial {{y}^{2}}}</math>  
+  </center>  
+  
+  The boundary conditions at <math>y\to \infty </math> can still be described by  
+  <center>  
+  <math>u=v=0,\text{ }y\to \infty </math>  
+  </center>  
+  <center>  
+  <math>T={{T}_{\infty }},\text{ }y\to \infty </math>  
+  </center>  
+  
+  However, the boundary conditions at the centerline, <math>x=0</math>, should be changed to:  
{ class="wikitable" border="0"  { class="wikitable" border="0"  
Line 32:  Line 51:  
    
 width="100%" <center>   width="100%" <center>  
  <math>{F}'''+\theta +(3+n)F{F}''(2n+2){{({F}')}^{2}}=0</math>  +  <big><math>{F}'''+\theta +(3+n)F{F}''(2n+2){{({F}')}^{2}}=0</math></big> 
</center>  </center>  
{{EquationRef(3)}}  {{EquationRef(3)}}  
Line 88:  Line 107:  
    
 width="100%" <center>   width="100%" <center>  
  <math>F(\eta )={F}''(\eta )=0,\text{ and }\theta (\eta )=1\text{ at }\eta =0</math>  +  <big><math>F(\eta )={F}''(\eta )=0,\text{ and }\theta (\eta )=1\text{ at }\eta =0</math></big> 
</center>  </center>  
{{EquationRef(9)}}  {{EquationRef(9)}}  
Line 101:  Line 120:  
}  }  
  Equations (7) – (10) can be solved numerically  +  Equations (7) – (10) can be solved numerically. With known intensity of the line heat source, the constant N in eq. (2) can be obtained from eq. (6): 
{ class="wikitable" border="0"  { class="wikitable" border="0"  
Line 112:  Line 131:  
where  where  
  
  
{ class="wikitable" border="0"  { class="wikitable" border="0"  
Line 124:  Line 141:  
is a function of Prandtl number. The values of I calculated at Pr = 0.7, 1.0, 6.7, and 10.0 by [[#ReferencesGebhart et al. (1970)]] are 1.245, 1.053, 0.407 and 0.328, respectively.  is a function of Prandtl number. The values of I calculated at Pr = 0.7, 1.0, 6.7, and 10.0 by [[#ReferencesGebhart et al. (1970)]] are 1.245, 1.053, 0.407 and 0.328, respectively.  
  
  
  
  
  
  
  
  
  
  
  
  
  
  
  
  
  
  
  
  
  
  
  
  
  
  
  
  
  
  
  
  
  
  
  
  
  
  
  
  
  
  
==References==  ==References==  
+  Faghri, A., Zhang, Y., and Howell, J. R., 2010, ''Advanced Heat and Mass Transfer'', Global Digital Press, Columbia, MO.  
Gebhart, B., Pera, L., and Schorr, A.W., 1970, “Steady Laminar Natural Convection Plumes above a Horizontal Line Heat Source,” International Journal of Heat and Mass Transfer, Vol. 13, pp. 161171  Gebhart, B., Pera, L., and Schorr, A.W., 1970, “Steady Laminar Natural Convection Plumes above a Horizontal Line Heat Source,” International Journal of Heat and Mass Transfer, Vol. 13, pp. 161171 
Current revision as of 09:28, 12 July 2010
The common feature of external convection discussed in the preceding subsections is that the flows are always near a heated (or cooled) solid wall. For some other applications, however, the thermally induced flow occurs without presence of a solid wall and is referred to as free boundary flow. When a point or line heat source is immersed into a bulk fluid, the fluid near the heat source is heated, becomes lighter and rises to form a plume [see Fig. 1 (a)]. The plume generated by a point heat source is axisymmetric, but a line heat source will create a twodimensional plume. A thermal is a column of rising air near the ground due to the uneven solar heating of the ground surface. The lighter air near the ground rises and cools due to expansion [see Fig. 1(b)].
Let us consider now the twodimensional free boundary flow induced by a line heat source (see Fig. 2). Since the problem is symmetric about x = 0, only half of the domain (x > 0) needs to be studied. In the coordinate system shown in Fig. 2, the continuity, momentum and energy equations are the same as those for natural convection near a vertical flat plate:
The boundary conditions at can still be described by
However, the boundary conditions at the centerline, x = 0, should be changed to:

which indicates that the velocity component in the xdirection and the temperature at the centerline are at their respective maximum. Gebhart et al. (1970) assumed that the difference between the centerline temperature and the bulk fluid temperature is of the following power law form:

where the index n will be determined using the overall energy balance. The continuity, momentum, and energy equations can be transformed to the following ordinary differential equations via the same similarity variables for natural convection over a vertical flat plate in Section 6.4.1 and considering eq. (2).
F''' + θ + (3 + n)FF'' − (2n + 2)(F')^{2} = 0 

Equations (3) and (4) are applicable to any natural convection problem that satisfies eq. (2). To get the ordinary differential equations that are specific for natural convection induced by a line heat source, proper values of N and n, as well as appropriate boundary conditions, must be specified. The energy convected across any horizontal plane in the plume is

which is identical to the intensity of the line heat source and it can be expressed in terms of similarity variables

Since the line heat source is the only source of heating, the above q' should be independent from x which can be true only if n is equal to − 3 / 5. Therefore, the ordinary differential equations for natural convection over a line heat source respectively become:


which are subject to the following boundary conditions:
F(η) = F''(η) = 0, and θ(η) = 1 at η = 0 

Equations (7) – (10) can be solved numerically. With known intensity of the line heat source, the constant N in eq. (2) can be obtained from eq. (6):

where

is a function of Prandtl number. The values of I calculated at Pr = 0.7, 1.0, 6.7, and 10.0 by Gebhart et al. (1970) are 1.245, 1.053, 0.407 and 0.328, respectively.
References
Faghri, A., Zhang, Y., and Howell, J. R., 2010, Advanced Heat and Mass Transfer, Global Digital Press, Columbia, MO.
Gebhart, B., Pera, L., and Schorr, A.W., 1970, “Steady Laminar Natural Convection Plumes above a Horizontal Line Heat Source,” International Journal of Heat and Mass Transfer, Vol. 13, pp. 161171