# Similarity Solution for Natural Convection on a Vertical Surface

### From Thermal-FluidsPedia

(Created page with '===6.4.1 Similarity Solution for Natural Convection on a Vertical Surface=== The boundary layer-type governing equations for the external convection problem shown in Fig. 6.1 ar…') |
|||

Line 1: | Line 1: | ||

- | |||

- | |||

The boundary layer-type governing equations for the external convection problem shown in Fig. 6.1 are eqs. (6.17), (6.20) and (6.21). Introducing the stream function <math>\psi </math>, and dimensionless temperature, <math>\theta </math>: | The boundary layer-type governing equations for the external convection problem shown in Fig. 6.1 are eqs. (6.17), (6.20) and (6.21). Introducing the stream function <math>\psi </math>, and dimensionless temperature, <math>\theta </math>: | ||

- | <center><math>u=\frac{\partial \psi }{\partial y},\text{ }v=-\frac{\partial \psi }{\partial x},\text{ }\theta =\frac{T-{{T}_{\infty }}}{{{T}_{w}}-{{T}_{\infty }}} \qquad \qquad() </math></center> | + | <center><math>u=\frac{\partial \psi }{\partial y},\text{ }v=-\frac{\partial \psi }{\partial x},\text{ }\theta =\frac{T-{{T}_{\infty }}}{{{T}_{w}}-{{T}_{\infty }}} \qquad \qquad(1) </math></center> |

(6.55) | (6.55) | ||

Line 11: | Line 9: | ||

- | <center><math>\frac{\partial \psi }{\partial y}\frac{{{\partial }^{2}}\psi }{\partial y\partial x}-\frac{\partial \psi }{\partial x}\frac{{{\partial }^{2}}\psi }{\partial {{y}^{2}}}=\nu \frac{{{\partial }^{3}}\psi }{\partial {{y}^{3}}}+g\beta ({{T}_{w}}-{{T}_{\infty }})\theta \qquad \qquad() </math></center> | + | <center><math>\frac{\partial \psi }{\partial y}\frac{{{\partial }^{2}}\psi }{\partial y\partial x}-\frac{\partial \psi }{\partial x}\frac{{{\partial }^{2}}\psi }{\partial {{y}^{2}}}=\nu \frac{{{\partial }^{3}}\psi }{\partial {{y}^{3}}}+g\beta ({{T}_{w}}-{{T}_{\infty }})\theta \qquad \qquad(2) </math></center> |

(6.56) | (6.56) | ||

- | <center><math>\frac{\partial \psi }{\partial y}\frac{\partial \theta }{\partial x}-\frac{\partial \psi }{\partial x}\frac{\partial \theta }{\partial y}=\alpha \frac{{{\partial }^{2}}\theta }{\partial {{y}^{2}}} \qquad \qquad() </math></center> | + | <center><math>\frac{\partial \psi }{\partial y}\frac{\partial \theta }{\partial x}-\frac{\partial \psi }{\partial x}\frac{\partial \theta }{\partial y}=\alpha \frac{{{\partial }^{2}}\theta }{\partial {{y}^{2}}} \qquad \qquad(3) </math></center> |

(6.57) | (6.57) | ||

Line 22: | Line 20: | ||

- | <center><math>\frac{\partial \psi }{\partial y}=\frac{\partial \psi }{\partial x}=0,\text{ }\theta =1\text{ at }y=0 \qquad \qquad() </math></center> | + | <center><math>\frac{\partial \psi }{\partial y}=\frac{\partial \psi }{\partial x}=0,\text{ }\theta =1\text{ at }y=0 \qquad \qquad(4) </math></center> |

(6.58) | (6.58) | ||

- | <center><math>\frac{\partial \psi }{\partial y}=\frac{\partial \psi }{\partial x}=0,\text{ }\theta =0\text{ at }y\to \infty \qquad \qquad() </math></center> | + | <center><math>\frac{\partial \psi }{\partial y}=\frac{\partial \psi }{\partial x}=0,\text{ }\theta =0\text{ at }y\to \infty \qquad \qquad(5) </math></center> |

(6.59) | (6.59) | ||

Line 33: | Line 31: | ||

- | <center><math>\eta =y\cdot H(x) \qquad \qquad() </math></center> | + | <center><math>\eta =y\cdot H(x) \qquad \qquad(6) </math></center> |

(6.60) | (6.60) | ||

- | where ''H(x)'' is an unspecified stretching function. The objective now is to reduce eqs. ( | + | where ''H(x)'' is an unspecified stretching function. The objective now is to reduce eqs. (2) and (3) to ordinary differential equations. The stream function, <math>\psi </math>, which is a function of ''x'' and ''y'', can be expressed as a function of ''x'' and <math>\eta </math>. If the similarity solution exists, one can express the stream function as |

- | <center><math>\psi (x,\eta )=\nu F(\eta )\cdot G(x) \qquad \qquad() </math></center> | + | <center><math>\psi (x,\eta )=\nu F(\eta )\cdot G(x) \qquad \qquad(7) </math></center> |

(6.61) | (6.61) | ||

Line 47: | Line 45: | ||

- | <center><math>\theta (x,\eta )=\theta (\eta ) \qquad \qquad() </math></center> | + | <center><math>\theta (x,\eta )=\theta (\eta ) \qquad \qquad(8) </math></center> |

(6.62) | (6.62) | ||

- | The derivatives in eqs. (6 | + | The derivatives in eqs. (6)–(8) can be obtained as shown below |

Line 75: | Line 73: | ||

- | where the primes for ''F'' and <math>\theta</math> denote the derivatives with respect to <math>\eta </math>, while the primes for ''G'' and ''H'' denote the derivatives with respect to ''x''. Substituting the above derivatives into eqs. ( | + | where the primes for ''F'' and <math>\theta</math> denote the derivatives with respect to <math>\eta </math>, while the primes for ''G'' and ''H'' denote the derivatives with respect to ''x''. Substituting the above derivatives into eqs. (2) and (3) and considering eq. (6), one obtains: |

- | <center><math>{F}'''+\frac{g\beta ({{T}_{w}}-{{T}_{\infty }})}{{{\nu }^{2}}{{H}^{3}}G}\theta +\frac{{{G}'}}{H}F{F}''-\left[ \frac{{H}'G}{{{H}^{2}}}+\frac{{{G}'}}{H} \right]{{({F}')}^{2}}=0 \qquad \qquad() </math></center> | + | <center><math>{F}'''+\frac{g\beta ({{T}_{w}}-{{T}_{\infty }})}{{{\nu }^{2}}{{H}^{3}}G}\theta +\frac{{{G}'}}{H}F{F}''-\left[ \frac{{H}'G}{{{H}^{2}}}+\frac{{{G}'}}{H} \right]{{({F}')}^{2}}=0 \qquad \qquad(9) </math></center> |

(6.63) | (6.63) | ||

- | <center><math>{\theta }''+\Pr \frac{{{G}'}}{H}F{\theta }'=0 \qquad \qquad() </math></center> | + | <center><math>{\theta }''+\Pr \frac{{{G}'}}{H}F{\theta }'=0 \qquad \qquad(10) </math></center> |

(6.64) | (6.64) | ||

- | In order to convert eqs. ( | + | In order to convert eqs. (9) and (10) to ordinary differential equations with <math>\eta </math> as the sole independent variable, all functions of ''x'' must be cancelled, which is possible only if the following combination of ''H'' and ''G'' and their derivatives are satisfied: |

- | <center><math>{{H}^{3}}G=A=\text{const} \qquad \qquad() </math></center> | + | <center><math>{{H}^{3}}G=A=\text{const} \qquad \qquad(11) </math></center> |

(6.65) | (6.65) | ||

- | <center><math>\frac{{{G}'}}{H}=B=\text{const} \qquad \qquad() </math></center> | + | <center><math>\frac{{{G}'}}{H}=B=\text{const} \qquad \qquad(12) </math></center> |

(6.66) | (6.66) | ||

- | <center><math>\frac{{H}'G}{{{H}^{2}}}=C=\text{const} \qquad \qquad() </math></center> | + | <center><math>\frac{{H}'G}{{{H}^{2}}}=C=\text{const} \qquad \qquad(13) </math></center> |

(6.67) | (6.67) | ||

- | Differentiation of eq. ( | + | Differentiation of eq. (11) and division of the resultant equation by <math>{H^4}</math> yields the following equation: |

- | <center><math>3\frac{{H}'G}{{{H}^{2}}}+\frac{{{G}'}}{H}=\text{0} \qquad \qquad() </math></center> | + | <center><math>3\frac{{H}'G}{{{H}^{2}}}+\frac{{{G}'}}{H}=\text{0} \qquad \qquad(14) </math></center> |

(6.68) | (6.68) | ||

- | which is satisfied if eqs. ( | + | which is satisfied if eqs. (12) and (13) are satisfied. Substituting eq. (14) into eq. (9), we have |

- | <center><math>{F}'''+\frac{g\beta ({{T}_{w}}-{{T}_{\infty }})}{{{\nu }^{2}}{{H}^{3}}G}\theta +\frac{{{G}'}}{H}F{F}''-\frac{2}{3}\frac{{{G}'}}{H}{{({F}')}^{2}}=0 \qquad \qquad() </math></center> | + | <center><math>{F}'''+\frac{g\beta ({{T}_{w}}-{{T}_{\infty }})}{{{\nu }^{2}}{{H}^{3}}G}\theta +\frac{{{G}'}}{H}F{F}''-\frac{2}{3}\frac{{{G}'}}{H}{{({F}')}^{2}}=0 \qquad \qquad(15) </math></center> |

(6.69) | (6.69) | ||

- | Therefore, satisfaction of eqs. ( | + | Therefore, satisfaction of eqs. (11) and (12) is sufficient to ensure that eqs. (9) and (10) become ordinary differential equations. Although any constants ''A'' and ''B'' in eqs. (11) and (12) will transform eqs. (15) and (10) into ordinary differential equations, the proper choice of values of these two constants will yield ordinary differential equations with simple forms. Let us choose |

Line 127: | Line 125: | ||

- | <center> <math>G(x)=4{{\left( \frac{1}{4}\text{G}{{\text{r}}_{x}} \right)}^{1/4}} \qquad \qquad() </math></center> | + | <center> <math>G(x)=4{{\left( \frac{1}{4}\text{G}{{\text{r}}_{x}} \right)}^{1/4}} \qquad \qquad(16) </math></center> |

(6.70) | (6.70) | ||

- | <center><math>H(x)=\frac{1}{x}{{\left( \frac{1}{4}\text{G}{{\text{r}}_{x}} \right)}^{1/4}} \qquad \qquad() </math></center> | + | <center><math>H(x)=\frac{1}{x}{{\left( \frac{1}{4}\text{G}{{\text{r}}_{x}} \right)}^{1/4}} \qquad \qquad(17) </math></center> |

(6.71) | (6.71) | ||

Line 138: | Line 136: | ||

- | <center><math>\text{G}{{\text{r}}_{x}}=\frac{g\beta ({{T}_{w}}-{{T}_{\infty }}){{x}^{3}}}{{{\nu }^{2}}} \qquad \qquad() </math></center> | + | <center><math>\text{G}{{\text{r}}_{x}}=\frac{g\beta ({{T}_{w}}-{{T}_{\infty }}){{x}^{3}}}{{{\nu }^{2}}} \qquad \qquad(18) </math></center> |

(6.72) | (6.72) | ||

Line 144: | Line 142: | ||

which is equivalent to the square of the Reynolds number based on the scale of the local velocity <math>{{u}_{0}}=\sqrt{g\beta ({{T}_{w}}-{{T}_{\infty }})x}</math>. | which is equivalent to the square of the Reynolds number based on the scale of the local velocity <math>{{u}_{0}}=\sqrt{g\beta ({{T}_{w}}-{{T}_{\infty }})x}</math>. | ||

- | Substituting eqs. ( | + | Substituting eqs. (16) and (17) into eqs. (15) and (10), the momentum and the energy equations become: |

- | <center><math>{F}'''+\theta +3F{F}''-2{{({F}')}^{2}}=0 \qquad \qquad() </math></center> | + | <center><math>{F}'''+\theta +3F{F}''-2{{({F}')}^{2}}=0 \qquad \qquad(19) </math></center> |

(6.73) | (6.73) | ||

- | <center><math>{\theta }''+3\Pr F{\theta }'=0 \qquad \qquad() </math></center> | + | <center><math>{\theta }''+3\Pr F{\theta }'=0 \qquad \qquad(20) </math></center> |

(6.74) | (6.74) | ||

Line 157: | Line 155: | ||

- | <center><math>u=\frac{\partial \psi }{\partial y}=\frac{\partial \psi }{\partial \eta }\frac{\partial \eta }{\partial y}=\frac{2\nu }{x}Gr_{x}^{1/2}{F}'(\eta ) \qquad \qquad() </math></center> | + | <center><math>u=\frac{\partial \psi }{\partial y}=\frac{\partial \psi }{\partial \eta }\frac{\partial \eta }{\partial y}=\frac{2\nu }{x}Gr_{x}^{1/2}{F}'(\eta ) \qquad \qquad(21) </math></center> |

(6.75) | (6.75) | ||

- | <center><math>v=-\frac{\partial \psi }{\partial x}=-\frac{\partial \psi }{\partial \eta }\frac{\partial \eta }{\partial x}=\frac{\nu }{x}{{\left( \frac{\text{G}{{\text{r}}_{x}}}{4} \right)}^{1/4}}[\eta {F}'(\eta )-3F(\eta ){{x}^{1/4}}] \qquad \qquad() </math></center> | + | <center><math>v=-\frac{\partial \psi }{\partial x}=-\frac{\partial \psi }{\partial \eta }\frac{\partial \eta }{\partial x}=\frac{\nu }{x}{{\left( \frac{\text{G}{{\text{r}}_{x}}}{4} \right)}^{1/4}}[\eta {F}'(\eta )-3F(\eta ){{x}^{1/4}}] \qquad \qquad(22) </math></center> |

(6.76) | (6.76) | ||

- | which are obtained from eqs. (6 | + | which are obtained from eqs. (6), (7), (16) and (17). Substituting the above expressions into eqs. (4) and (5), the boundary conditions for eqs. (19) and (20) can be obtained as follows |

- | <center><math>F(\eta )={F}'(\eta )=0\text{ and }\theta (\eta )=1\text{ at }\eta =0 \qquad \qquad() </math></center> | + | <center><math>F(\eta )={F}'(\eta )=0\text{ and }\theta (\eta )=1\text{ at }\eta =0 \qquad \qquad(23) </math></center> |

(6.77) | (6.77) | ||

- | <center><math>{F}'(\eta )=\theta (\eta )=0\text{ at }\eta \to \infty \qquad \qquad() </math></center> | + | <center><math>{F}'(\eta )=\theta (\eta )=0\text{ at }\eta \to \infty \qquad \qquad(24) </math></center> |

(6.78) | (6.78) | ||

- | Equations ( | + | Equations (19) and (20) are coupled nonlinear ordinary equations with boundary conditions specified at different <math>\eta</math> and they thus make external natural convection a boundary value problem. The original numerical solution was obtained by [[#References|Ostrach (1953)]] for a wide range of Prandtl numbers from 0.01 to 1000. Figure 1 shows the dimensionless velocity and temperature for various Prandtl numbers obtained by using a shooting method with the Range-Kutta method. As the Prandtl number increases, the maximum velocity in the boundary layer decreases, and the location at which the peak velocity occurs shifts to smaller <math>\eta</math>. As shown in Fig. 1(a), the thickness of the momentum boundary layer decreases with increasing Prandtl number. Similarly, the thickness of the thermal boundary layer also decreases with increasing Prandtl number, as indicated by Fig. 1(b). The absolute value of the dimensionless temperature gradient also increases with increasing Prandtl number. |

After the dimensionless temperature in the boundary layer is obtained, the local heat transfer coefficient at the surface of the vertical plate can be obtained from eq. (6.3) as follows: | After the dimensionless temperature in the boundary layer is obtained, the local heat transfer coefficient at the surface of the vertical plate can be obtained from eq. (6.3) as follows: | ||

- | <center><math>{{h}_{x}}=-\frac{k}{{{T}_{w}}-{{T}_{\infty }}}{{\left( \frac{\partial T}{\partial y} \right)}_{y=0}}=-k{\theta }'(0)\frac{1}{x}{{\left( \frac{\text{G}{{\text{r}}_{x}}}{4} \right)}^{1/4}} \qquad \qquad() </math></center> | + | <center><math>{{h}_{x}}=-\frac{k}{{{T}_{w}}-{{T}_{\infty }}}{{\left( \frac{\partial T}{\partial y} \right)}_{y=0}}=-k{\theta }'(0)\frac{1}{x}{{\left( \frac{\text{G}{{\text{r}}_{x}}}{4} \right)}^{1/4}} \qquad \qquad(15) </math></center> |

(6.79) | (6.79) | ||

Line 188: | Line 186: | ||

- | <center><math>\text{N}{{\text{u}}_{x}}=\frac{{{h}_{x}}x}{k}=-\frac{{\theta }'(0)}{\sqrt{2}}\text{Gr}_{x}^{1/4}=\phi (\Pr )\text{Gr}_{x}^{1/4} \qquad \qquad() </math></center> | + | <center><math>\text{N}{{\text{u}}_{x}}=\frac{{{h}_{x}}x}{k}=-\frac{{\theta }'(0)}{\sqrt{2}}\text{Gr}_{x}^{1/4}=\phi (\Pr )\text{Gr}_{x}^{1/4} \qquad \qquad(26) </math></center> |

(6.80) | (6.80) | ||

- | |||

- | |||

- | [[Image: chapter6g_(3).gif |thumb|400 px|alt= Velocity and temperature profile in the boundary layer for external natural convection over a vertical isothermal surface | Figure | + | [[Image: chapter6g_(3).gif |thumb|400 px|alt= Velocity and temperature profile in the boundary layer for external natural convection over a vertical isothermal surface | Figure 1: Velocity and temperature profile in the boundary layer for external natural convection over a vertical isothermal surface]] |

- | where <math>\phi (\Pr )=-{\theta }'(0)/\sqrt{2}</math> is a function of Prandtl number. The dependence of <math>\phi </math> on the Prandtl number is evidenced by eq. ( | + | where <math>\phi (\Pr )=-{\theta }'(0)/\sqrt{2}</math> is a function of Prandtl number. The dependence of <math>\phi </math> on the Prandtl number is evidenced by eq. (20) and by Fig. 1(b). The values of <math>\phi </math> for various Pr have been obtained numerically by [[#References|Ostrach (1953). Ede (1964)]] proposed the following function that correlated the numerical results: |

- | <center> <math>\phi (\Pr )=\frac{3}{4}{{\left[ \frac{2{{\Pr }^{2}}}{5(1+2{{\Pr }^{1/2}}+2\Pr )} \right]}^{1/4}}</math></center> | + | <center><math>\phi (\Pr )=\frac{3}{4}{{\left[ \frac{2{{\Pr }^{2}}}{5(1+2{{\Pr }^{1/2}}+2\Pr )} \right]}^{1/4}}</math></center> |

The local Nusselt number thus becomes: | The local Nusselt number thus becomes: | ||

- | <center><math>\text{N}{{\text{u}}_{x}}=\phi (\Pr )\text{Gr}_{x}^{1/4}=\frac{3}{4}{{\left[ \frac{2\Pr }{5(1+2{{\Pr }^{1/2}}+2\Pr )} \right]}^{1/4}}{{\text{(G}{{\text{r}}_{x}}\text{Pr)}}^{\text{1/4}}} \qquad \qquad() </math></center> | + | <center><math>\text{N}{{\text{u}}_{x}}=\phi (\Pr )\text{Gr}_{x}^{1/4}=\frac{3}{4}{{\left[ \frac{2\Pr }{5(1+2{{\Pr }^{1/2}}+2\Pr )} \right]}^{1/4}}{{\text{(G}{{\text{r}}_{x}}\text{Pr)}}^{\text{1/4}}} \qquad \qquad(27) </math></center> |

(6.81) | (6.81) | ||

Line 211: | Line 207: | ||

- | <center><math>\text{N}{{\text{u}}_{x}}=\frac{3}{4}{{\left[ \frac{2\Pr }{5(1+2{{\Pr }^{1/2}}+2\Pr )} \right]}^{1/4}}\text{Ra}_{x}^{\text{1/4}} \qquad \qquad() </math></center> | + | <center><math>\text{N}{{\text{u}}_{x}}=\frac{3}{4}{{\left[ \frac{2\Pr }{5(1+2{{\Pr }^{1/2}}+2\Pr )} \right]}^{1/4}}\text{Ra}_{x}^{\text{1/4}} \qquad \qquad(28) </math></center> |

(6.82) | (6.82) | ||

- | Equations ( | + | Equations (27) and (28) are valid for <math>0<\Pr <\infty </math>. |

+ | |||

As is demonstrated above, the similarity solution is obtained as a consequence of the geometrical similarity of the velocity and temperature profiles in the boundary layers, i.e., the velocity and temperature profiles vary with ''x'' according to the stretching functions, ''G(x) '' and ''H(x) ''. If the geometrical | As is demonstrated above, the similarity solution is obtained as a consequence of the geometrical similarity of the velocity and temperature profiles in the boundary layers, i.e., the velocity and temperature profiles vary with ''x'' according to the stretching functions, ''G(x) '' and ''H(x) ''. If the geometrical | ||

- | [[Image: chapter6g_(4).gif |thumb|400 px|alt= Velocity and temperature profiles in the boundary layer based on modified scale| Figure | + | [[Image: chapter6g_(4).gif |thumb|400 px|alt= Velocity and temperature profiles in the boundary layer based on modified scale| Figure 2: Velocity and temperature profiles in the boundary layer based on modified scale.]] |

- | similarity exists, the selection of the stretching functions is not unique, and the resulting solutions based on different choices, as long as they all satisfy eqs. ( | + | similarity exists, the selection of the stretching functions is not unique, and the resulting solutions based on different choices, as long as they all satisfy eqs. (11)–(13), are equivalent to the solution based on the stretching functions expressed in eqs. (16) and (17). However, a choice that better represents the physics of the problem will lead to the results presented in a physically more meaningful way. Figure 6.4 was obtained by modifying the numerical results shown in Fig. 1 by incorporating the modified stretching function for fluids of Pr > 1 [[#References|(Bejan, 2004)]] expressed in eqs. (6.43) and (6.44), which is equivalent to a selection of the stretching functions in the following form |

Line 229: | Line 226: | ||

- | It is clear in Fig. | + | It is clear in Fig. 2 that, in the limit Pr → ∞, the temperature profiles collapse onto a single curve, while the dimensionless velocity peak for fluids of Pr > 1 is consistently a number of order 1, showing that the velocity peak falls in the thermal boundary layer. Furthermore, as Pr increases, the velocity profile extends farther and farther into the isothermal fluid. All these features are anticipated by Fig. 6.1 and support the scale analysis, but cannot be seen from Fig. 1, in which the velocity and temperature profiles constantly shift as Pr changes, and the peak dimensionless velocity is not of order 1. |

- | <center><math>\Pr ({{{F}'''}_{1}}+{{\theta }_{1}})+{{\theta }_{1}}+{{F}_{1}}{{{F}''}_{1}}=\frac{2}{3}{F}'_{1}^{2} \qquad \qquad() </math></center> | + | <center><math>\Pr ({{{F}'''}_{1}}+{{\theta }_{1}})+{{\theta }_{1}}+{{F}_{1}}{{{F}''}_{1}}=\frac{2}{3}{F}'_{1}^{2} \qquad \qquad(29) </math></center> |

(6.83) | (6.83) | ||

- | <center><math>{{{\theta }''}_{1}}+{{F}_{1}}{{{\theta }'}_{1}}=0 \qquad \qquad() </math></center> | + | <center><math>{{{\theta }''}_{1}}+{{F}_{1}}{{{\theta }'}_{1}}=0 \qquad \qquad(30) </math></center> |

(6.84) | (6.84) | ||

- | Equations ( | + | Equations (29) and (30) allow simplification for the cases that <math>\Pr \to 0</math> or <math>\Pr \to \infty </math> (see Problem 6.4). Numerical solution of the simplified equation for these two cases yields the following results [[#References|(Le Fevre, 1956)]]: |

Line 247: | Line 244: | ||

0.600{{(\text{G}{{\text{r}}_{x}}{{\Pr }^{2}})}^{1/4}}=0.600{{(\text{R}{{\text{a}}_{x}}\Pr )}^{1/4}} & \text{as Pr}\to 0 \\ | 0.600{{(\text{G}{{\text{r}}_{x}}{{\Pr }^{2}})}^{1/4}}=0.600{{(\text{R}{{\text{a}}_{x}}\Pr )}^{1/4}} & \text{as Pr}\to 0 \\ | ||

0.503{{(\text{G}{{\text{r}}_{x}}\Pr )}^{1/4}}=0.503\text{Ra}_{x}^{1/4}\text{ } & \text{as Pr}\to \infty \\ | 0.503{{(\text{G}{{\text{r}}_{x}}\Pr )}^{1/4}}=0.503\text{Ra}_{x}^{1/4}\text{ } & \text{as Pr}\to \infty \\ | ||

- | \end{matrix} \right. \qquad \qquad() </math></center> | + | \end{matrix} \right. \qquad \qquad(31) </math></center> |

(6.85) | (6.85) | ||

- | Equations ( | + | Equations (28) and (31) indicate that <math>{{h}_{x}}\propto {{x}^{-1/4}}</math>. The average Nusselt number over the entire vertical wall, <math>{{\overline{\text{Nu}}}_{L}}</math>, is related to the local Nusselt number at <math>x=L</math> by the following (see Problem 6.6): |

- | <center> | + | <center><math>{{\overline{\text{Nu}}}_{L}}=\frac{4}{3}\text{N}{{\text{u}}_{L}} \qquad \qquad(32) </math></center> |

- | + | ||

(6.86) | (6.86) | ||

## Revision as of 15:36, 19 June 2010

The boundary layer-type governing equations for the external convection problem shown in Fig. 6.1 are eqs. (6.17), (6.20) and (6.21). Introducing the stream function ψ, and dimensionless temperature, θ:

(6.55)

the continuity equation (6.17) is satisfied and the momentum and energy equations (6.20) and (6.21) become:

(6.56)

(6.57)

with the following boundary conditions:

(6.58)

(6.59)

The idea behind the similarity solution is that the velocity and temperature profiles in different *x* in the boundary layers are geometrically similar, differing only by a stretching factor in the *x*-direction (Kays et al., 2005). Thus, the similarity variable should have the following form:

(6.60)

where *H(x)* is an unspecified stretching function. The objective now is to reduce eqs. (2) and (3) to ordinary differential equations. The stream function, ψ, which is a function of *x* and *y*, can be expressed as a function of *x* and η. If the similarity solution exists, one can express the stream function as

(6.61)

where *F*(η) and *G*(*x*) are the similarity function and the stretching function, respectively. The dimensionless temperature can be assumed as a function of η only, i.e.,

(6.62)

The derivatives in eqs. (6)–(8) can be obtained as shown below

where the primes for *F* and θ denote the derivatives with respect to η, while the primes for *G* and *H* denote the derivatives with respect to *x*. Substituting the above derivatives into eqs. (2) and (3) and considering eq. (6), one obtains:

(6.63)

(6.64)

In order to convert eqs. (9) and (10) to ordinary differential equations with η as the sole independent variable, all functions of *x* must be cancelled, which is possible only if the following combination of *H* and *G* and their derivatives are satisfied:

(6.65)

(6.66)

(6.67)

Differentiation of eq. (11) and division of the resultant equation by *H*^{4} yields the following equation:

(6.68)

which is satisfied if eqs. (12) and (13) are satisfied. Substituting eq. (14) into eq. (9), we have

(6.69)

Therefore, satisfaction of eqs. (11) and (12) is sufficient to ensure that eqs. (9) and (10) become ordinary differential equations. Although any constants *A* and *B* in eqs. (11) and (12) will transform eqs. (15) and (10) into ordinary differential equations, the proper choice of values of these two constants will yield ordinary differential equations with simple forms. Let us choose

The corresponding *G* and *H* functions then become:

(6.70)

(6.71)

where *G**r*_{x} is the *local* Grashof number defined as:

(6.72)

which is equivalent to the square of the Reynolds number based on the scale of the local velocity .

Substituting eqs. (16) and (17) into eqs. (15) and (10), the momentum and the energy equations become:

(6.73)

(6.74)

The velocity components in the *x*- and *y*-directions can be expressed in terms of similarity variables by the following equations:

(6.75)

(6.76)

which are obtained from eqs. (6), (7), (16) and (17). Substituting the above expressions into eqs. (4) and (5), the boundary conditions for eqs. (19) and (20) can be obtained as follows

(6.77)

(6.78)

Equations (19) and (20) are coupled nonlinear ordinary equations with boundary conditions specified at different η and they thus make external natural convection a boundary value problem. The original numerical solution was obtained by Ostrach (1953) for a wide range of Prandtl numbers from 0.01 to 1000. Figure 1 shows the dimensionless velocity and temperature for various Prandtl numbers obtained by using a shooting method with the Range-Kutta method. As the Prandtl number increases, the maximum velocity in the boundary layer decreases, and the location at which the peak velocity occurs shifts to smaller η. As shown in Fig. 1(a), the thickness of the momentum boundary layer decreases with increasing Prandtl number. Similarly, the thickness of the thermal boundary layer also decreases with increasing Prandtl number, as indicated by Fig. 1(b). The absolute value of the dimensionless temperature gradient also increases with increasing Prandtl number.

After the dimensionless temperature in the boundary layer is obtained, the local heat transfer coefficient at the surface of the vertical plate can be obtained from eq. (6.3) as follows:

(6.79)

The local Nusselt number is

(6.80)

where is a function of Prandtl number. The dependence of φ on the Prandtl number is evidenced by eq. (20) and by Fig. 1(b). The values of φ for various Pr have been obtained numerically by Ostrach (1953). Ede (1964) proposed the following function that correlated the numerical results:

The local Nusselt number thus becomes:

(6.81)

which can also be rewritten in terms of Rayleigh number

(6.82)

Equations (27) and (28) are valid for .

As is demonstrated above, the similarity solution is obtained as a consequence of the geometrical similarity of the velocity and temperature profiles in the boundary layers, i.e., the velocity and temperature profiles vary with *x* according to the stretching functions, *G(x) * and *H(x) *. If the geometrical

similarity exists, the selection of the stretching functions is not unique, and the resulting solutions based on different choices, as long as they all satisfy eqs. (11)–(13), are equivalent to the solution based on the stretching functions expressed in eqs. (16) and (17). However, a choice that better represents the physics of the problem will lead to the results presented in a physically more meaningful way. Figure 6.4 was obtained by modifying the numerical results shown in Fig. 1 by incorporating the modified stretching function for fluids of Pr > 1 (Bejan, 2004) expressed in eqs. (6.43) and (6.44), which is equivalent to a selection of the stretching functions in the following form

It is clear in Fig. 2 that, in the limit Pr → ∞, the temperature profiles collapse onto a single curve, while the dimensionless velocity peak for fluids of Pr > 1 is consistently a number of order 1, showing that the velocity peak falls in the thermal boundary layer. Furthermore, as Pr increases, the velocity profile extends farther and farther into the isothermal fluid. All these features are anticipated by Fig. 6.1 and support the scale analysis, but cannot be seen from Fig. 1, in which the velocity and temperature profiles constantly shift as Pr changes, and the peak dimensionless velocity is not of order 1.

**Failed to parse (PNG conversion failed;**

**check for correct installation of latex, dvips, gs, and convert): \Pr ({{{F}'''}_{1}}+{{\theta }_{1}})+{{\theta }_{1}}+{{F}_{1}}{{{F}''}_{1}}=\frac{2}{3}{F}'_{1}^{2} \qquad \qquad(29) **

(6.83)

(6.84)

Equations (29) and (30) allow simplification for the cases that or (see Problem 6.4). Numerical solution of the simplified equation for these two cases yields the following results (Le Fevre, 1956):

(6.85)

Equations (28) and (31) indicate that . The average Nusselt number over the entire vertical wall, , is related to the local Nusselt number at *x* = *L* by the following (see Problem 6.6):

(6.86)