# Similarity Solution for Natural Convection on a Vertical Surface

(Difference between revisions)
 Revision as of 03:13, 17 June 2010 (view source) (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…')← Older edit Current revision as of 20:16, 23 June 2010 (view source) (4 intermediate revisions not shown) Line 1: Line 1: - ===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. 1 from [[Governing Equations for Natural Convection]] are eqs. $\frac{\partial u}{\partial x}+\frac{\partial v}{\partial y}=0$, $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 }})$ and $u\frac{\partial T}{\partial x}+v\frac{\partial T}{\partial y}=\alpha \frac{{{\partial }^{2}}T}{\partial {{y}^{2}}}$ from [[External Natural Convection from Heated Vertical Plate]]. Introducing the stream function $\psi$, and dimensionless temperature, $\theta$: - + - 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 $\psi$, and dimensionless temperature, $\theta$: + -
$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()$
+
$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)$
- (6.55) + - the continuity equation (6.17) is satisfied and the momentum and energy equations (6.20) and (6.21) become: + the continuity equation $\frac{\partial u}{\partial x}+\frac{\partial v}{\partial y}=0$ is satisfied and the momentum and energy equations $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 }})$ and $u\frac{\partial T}{\partial x}+v\frac{\partial T}{\partial y}=\alpha \frac{{{\partial }^{2}}T}{\partial {{y}^{2}}}$ from [[External Natural Convection from Heated Vertical Plate]] become: -
$\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()$
+
$\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)$
- (6.56) + -
$\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()$
+
$\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)$
- (6.57) + Line 22: Line 19: -
$\frac{\partial \psi }{\partial y}=\frac{\partial \psi }{\partial x}=0,\text{ }\theta =1\text{ at }y=0 \qquad \qquad()$
+
$\frac{\partial \psi }{\partial y}=\frac{\partial \psi }{\partial x}=0,\text{ }\theta =1\text{ at }y=0 \qquad \qquad(4)$
- (6.58) + -
$\frac{\partial \psi }{\partial y}=\frac{\partial \psi }{\partial x}=0,\text{ }\theta =0\text{ at }y\to \infty \qquad \qquad()$
+
$\frac{\partial \psi }{\partial y}=\frac{\partial \psi }{\partial x}=0,\text{ }\theta =0\text{ at }y\to \infty \qquad \qquad(5)$
- (6.59) + Line 33: Line 30: -
$\eta =y\cdot H(x) \qquad \qquad()$
+
$\eta =y\cdot H(x) \qquad \qquad(6)$
- (6.60) + - where ''H(x)'' is an unspecified stretching function. The objective now is to reduce eqs. (6.56) and (6.57) to ordinary differential equations. The stream function, $\psi$, which is a function of ''x'' and ''y'', can be expressed as a function of ''x'' and $\eta$. If the similarity solution exists, one can express the stream function as + 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, $\psi$, which is a function of ''x'' and ''y'', can be expressed as a function of ''x'' and $\eta$. If the similarity solution exists, one can express the stream function as -
$\psi (x,\eta )=\nu F(\eta )\cdot G(x) \qquad \qquad()$
+
$\psi (x,\eta )=\nu F(\eta )\cdot G(x) \qquad \qquad(7)$
- (6.61) + Line 47: Line 44: -
$\theta (x,\eta )=\theta (\eta ) \qquad \qquad()$
+
$\theta (x,\eta )=\theta (\eta ) \qquad \qquad(8)$
- (6.62) + - The derivatives in eqs. (6.60)–(6.62) can be obtained as shown below + The derivatives in eqs. (6)–(8) can be obtained as shown below Line 75: Line 72: - where the primes for ''F'' and $\theta$ denote the derivatives with respect to $\eta$, while the primes for ''G'' and ''H'' denote the derivatives with respect to ''x''. Substituting the above derivatives into eqs. (6.56) and (6.57) and considering eq. (6.60), one obtains: + where the primes for ''F'' and $\theta$ denote the derivatives with respect to $\eta$, 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: -
${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()$
+
${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)$
- (6.63) + -
${\theta }''+\Pr \frac{{{G}'}}{H}F{\theta }'=0 \qquad \qquad()$
+
${\theta }''+\Pr \frac{{{G}'}}{H}F{\theta }'=0 \qquad \qquad(10)$
- (6.64) + - In order to convert eqs. (6.63) and (6.64) to ordinary differential equations with $\eta$ 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: + In order to convert eqs. (9) and (10) to ordinary differential equations with $\eta$ 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: -
${{H}^{3}}G=A=\text{const} \qquad \qquad()$
+
${{H}^{3}}G=A=\text{const} \qquad \qquad(11)$
- (6.65) + -
$\frac{{{G}'}}{H}=B=\text{const} \qquad \qquad()$
+
$\frac{{{G}'}}{H}=B=\text{const} \qquad \qquad(12)$
- (6.66) + -
$\frac{{H}'G}{{{H}^{2}}}=C=\text{const} \qquad \qquad()$
+
$\frac{{H}'G}{{{H}^{2}}}=C=\text{const} \qquad \qquad(13)$
- (6.67) + - Differentiation of eq. (6.65) and division of the resultant equation by ${H^4}$ yields the following equation: + Differentiation of eq. (11) and division of the resultant equation by ${H^4}$ yields the following equation: -
$3\frac{{H}'G}{{{H}^{2}}}+\frac{{{G}'}}{H}=\text{0} \qquad \qquad()$
+
$3\frac{{H}'G}{{{H}^{2}}}+\frac{{{G}'}}{H}=\text{0} \qquad \qquad(14)$
- (6.68) + - which is satisfied if eqs. (6.66) and (6.67) are satisfied. Substituting eq. (6.68) into eq. (6.63), we have + which is satisfied if eqs. (12) and (13) are satisfied. Substituting eq. (14) into eq. (9), we have -
${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()$
+
${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)$
- (6.69) + - Therefore, satisfaction of eqs. (6.65) and (6.66) is sufficient to ensure that eqs. (6.63) and (6.64) become ordinary differential equations. Although any constants ''A'' and ''B'' in eqs. (6.65) and (6.66) will transform eqs. (6.69) and (6.64) into ordinary differential equations, the proper choice of values of these two constants will yield ordinary differential equations with simple forms. Let us choose + 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 124: -
$G(x)=4{{\left( \frac{1}{4}\text{G}{{\text{r}}_{x}} \right)}^{1/4}} \qquad \qquad()$
+
$G(x)=4{{\left( \frac{1}{4}\text{G}{{\text{r}}_{x}} \right)}^{1/4}} \qquad \qquad(16)$
- (6.70) + -
$H(x)=\frac{1}{x}{{\left( \frac{1}{4}\text{G}{{\text{r}}_{x}} \right)}^{1/4}} \qquad \qquad()$
+
$H(x)=\frac{1}{x}{{\left( \frac{1}{4}\text{G}{{\text{r}}_{x}} \right)}^{1/4}} \qquad \qquad(17)$
- (6.71) + Line 138: Line 135: -
$\text{G}{{\text{r}}_{x}}=\frac{g\beta ({{T}_{w}}-{{T}_{\infty }}){{x}^{3}}}{{{\nu }^{2}}} \qquad \qquad()$
+
$\text{G}{{\text{r}}_{x}}=\frac{g\beta ({{T}_{w}}-{{T}_{\infty }}){{x}^{3}}}{{{\nu }^{2}}} \qquad \qquad(18)$
- (6.72) + which is equivalent to the square of the Reynolds number based on the scale of the local velocity ${{u}_{0}}=\sqrt{g\beta ({{T}_{w}}-{{T}_{\infty }})x}$. which is equivalent to the square of the Reynolds number based on the scale of the local velocity ${{u}_{0}}=\sqrt{g\beta ({{T}_{w}}-{{T}_{\infty }})x}$. - Substituting eqs. (6.70) and (6.71) into eqs. (6.69) and (6.64), the momentum and the energy equations become: + Substituting eqs. (16) and (17) into eqs. (15) and (10), the momentum and the energy equations become: -
${F}'''+\theta +3F{F}''-2{{({F}')}^{2}}=0 \qquad \qquad()$
+
${F}'''+\theta +3F{F}''-2{{({F}')}^{2}}=0 \qquad \qquad(19)$
- (6.73) + -
${\theta }''+3\Pr F{\theta }'=0 \qquad \qquad()$
+
${\theta }''+3\Pr F{\theta }'=0 \qquad \qquad(20)$
- (6.74) + Line 157: Line 154: -
$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()$
+
$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)$
- (6.75) + -
$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()$
+
$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)$
- (6.76) + - which are obtained from eqs. (6.60), (6.61), (6.70) and (6.71). Substituting the above expressions into eqs. (6.58) and (6.59), the boundary conditions for eqs. (6.73) and (6.74) can be obtained as follows + 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 -
$F(\eta )={F}'(\eta )=0\text{ and }\theta (\eta )=1\text{ at }\eta =0 \qquad \qquad()$
+
$F(\eta )={F}'(\eta )=0\text{ and }\theta (\eta )=1\text{ at }\eta =0 \qquad \qquad(23)$
- (6.77) + -
${F}'(\eta )=\theta (\eta )=0\text{ at }\eta \to \infty \qquad \qquad()$
+
${F}'(\eta )=\theta (\eta )=0\text{ at }\eta \to \infty \qquad \qquad(24)$
- (6.78) + - Equations (6.73) and (6.74) are coupled nonlinear ordinary equations with boundary conditions specified at different $\eta$ 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 6.3 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 $\eta$. As shown in Fig. 6.3(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. 6.3(b). The absolute value of the dimensionless temperature gradient also increases with increasing Prandtl number. + Equations (19) and (20) are coupled nonlinear ordinary equations with boundary conditions specified at different $\eta$ 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 $\eta$. 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. ${{h}_{x}}=-\frac{k}{{{T}_{w}}-{{T}_{\infty }}}{{\left( \frac{\partial T}{\partial y} \right)}_{y=0}}$ from [[Governing Equations for Natural Convection]] as follows: -
${{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()$
+
${{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)$
- (6.79) + Line 188: Line 185: -
$\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()$
+
$\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)$
- (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 6.3: Velocity and temperature profile in the boundary layer for external natural convection over a vertical isothermal surface]] + [[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 $\phi (\Pr )=-{\theta }'(0)/\sqrt{2}$ is a function of Prandtl number. The dependence of $\phi$ on the Prandtl number is evidenced by eq. (6.74) and by Fig. 6.3(b). The values of $\phi$ for various Pr have been obtained numerically by [[#References|Ostrach (1953). Ede (1964)]] proposed the following function that correlated the numerical results: + where $\phi (\Pr )=-{\theta }'(0)/\sqrt{2}$ is a function of Prandtl number. The dependence of $\phi$ on the Prandtl number is evidenced by eq. (20) and by Fig. 1(b). The values of $\phi$ for various Pr have been obtained numerically by [[#References|Ostrach (1953). Ede (1964)]] proposed the following function that correlated the numerical results: -
$\phi (\Pr )=\frac{3}{4}{{\left[ \frac{2{{\Pr }^{2}}}{5(1+2{{\Pr }^{1/2}}+2\Pr )} \right]}^{1/4}}$
+
$\phi (\Pr )=\frac{3}{4}{{\left[ \frac{2{{\Pr }^{2}}}{5(1+2{{\Pr }^{1/2}}+2\Pr )} \right]}^{1/4}}$
The local Nusselt number thus becomes: The local Nusselt number thus becomes: -
$\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()$
+
$\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)$
- (6.81) + Line 211: Line 206: -
$\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()$
+
$\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)$
- (6.82) + - Equations (6.81) and (6.82) are valid for $0<\Pr <\infty$. + Equations (27) and (28) are valid for $0<\Pr <\infty$. + 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 6.4: Velocity and temperature profiles in the boundary layer based on modified scale.]] + [[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. (6.65)–(6.67), are equivalent to the solution based on the stretching functions expressed in eqs. (6.70) and (6.71). 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. 6.3 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 + 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 2 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. ${{\delta }_{t}}\sim L\text{Ra}_{L}^{-1/4} \qquad \qquad(1)$ and $u\sim \frac{\alpha }{L}\text{Ra}_{L}^{1/2}$ from [[High Prandtl Number Fluids]], which is equivalent to a selection of the stretching functions in the following form Line 229: Line 225: - It is clear in Fig. 6.4 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. 6.3, in which the velocity and temperature profiles constantly shift as Pr changes, and the peak dimensionless velocity is not of order 1. + 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. 1 from [[Governing Equations for Natural Convection]] 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. -
$\Pr ({{{F}'''}_{1}}+{{\theta }_{1}})+{{\theta }_{1}}+{{F}_{1}}{{{F}''}_{1}}=\frac{2}{3}{F}'_{1}^{2} \qquad \qquad()$
+
$\Pr ({{{{F}'''}}_{1}}+{{\theta }_{1}})+{{\theta }_{1}}+{{F}_{1}}{{{{F}''}}_{1}}=\frac{2}{3}{{F}'}_{1}^{2} \qquad \qquad(29)$
- (6.83) + -
${{{\theta }''}_{1}}+{{F}_{1}}{{{\theta }'}_{1}}=0 \qquad \qquad()$
+
${{{\theta }''}_{1}}+{{F}_{1}}{{{\theta }'}_{1}}=0 \qquad \qquad(30)$
- (6.84) + - Equations (6.83) and (6.84) allow simplification for the cases that $\Pr \to 0$ or $\Pr \to \infty$ (see Problem 6.4). Numerical solution of the simplified equation for these two cases yields the following results [[#References|(Le Fevre, 1956)]]: + Equations (29) and (30) allow simplification for the cases that $\Pr \to 0$ or $\Pr \to \infty$ (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 243: 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() [/itex] + \end{matrix} \right. \qquad \qquad(31) [/itex] - (6.85) + - Equations (6.82) and (6.85) indicate that ${{h}_{x}}\propto {{x}^{-1/4}}$. The average Nusselt number over the entire vertical wall, + Equations (28) and (31) indicate that ${{h}_{x}}\propto {{x}^{-1/4}}$. The average Nusselt number over the entire vertical wall, ${{\overline{\text{Nu}}}_{L}}$, is related to the local Nusselt number at $x=L$ by the following (see Problem 6.6): -
${{\overline{\text{Nu}}}_{L}}$, is related to the local Nusselt number at $x=L$ by the following (see Problem 6.6): +
${{\overline{\text{Nu}}}_{L}}=\frac{4}{3}\text{N}{{\text{u}}_{L}} \qquad \qquad(32)$
- ${{\overline{\text{Nu}}}_{L}}=\frac{4}{3}\text{N}{{\text{u}}_{L}} \qquad \qquad()$
+ - (6.86) + ==References== ==References== + + Bejan, A., 2004, ''Convection Heat Transfer'', 3rd ed., John Wiley & Sons, New York. + + Ede, A.J., 1964, “Advances in Free Convection,” ''Advances in Heat Transfer'', Vol. 4, pp. 1-64, Academic Press, New York. + + Kays, W.M., Crawford, M.E., and Weigand, B., 2005, ''Convection Heat Transfer'', 4th ed., McGraw-Hill, New York, NY. + + Le Fevre, E.J., 1956, “Laminar Free Convection from a Vertical Plane Surface,” ''Proceedings of the 9th International Congress of Applied Mechanics'', Brussels, Vol. 4, pp. 168-174. + + Ostrach, S., 1953, An Analysis of Laminar Free Convection Flow and Heat Transfer about a Flat Plate Parallel to the Direction of the Generating Body Force, ''NASA Report'', 1111, pp. 63–79. ==Further Reading== ==Further Reading== ==External Links== ==External Links==

## Current revision as of 20:16, 23 June 2010

The boundary layer-type governing equations for the external convection problem shown in Fig. 1 from Governing Equations for Natural Convection are eqs. $\frac{\partial u}{\partial x}+\frac{\partial v}{\partial y}=0$, $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 }})$ and $u\frac{\partial T}{\partial x}+v\frac{\partial T}{\partial y}=\alpha \frac{{{\partial }^{2}}T}{\partial {{y}^{2}}}$ from External Natural Convection from Heated Vertical Plate. Introducing the stream function ψ, and dimensionless temperature, θ:

$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)$

the continuity equation $\frac{\partial u}{\partial x}+\frac{\partial v}{\partial y}=0$ is satisfied and the momentum and energy equations $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 }})$ and $u\frac{\partial T}{\partial x}+v\frac{\partial T}{\partial y}=\alpha \frac{{{\partial }^{2}}T}{\partial {{y}^{2}}}$ from External Natural Convection from Heated Vertical Plate become:

$\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)$

$\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)$

with the following boundary conditions:

$\frac{\partial \psi }{\partial y}=\frac{\partial \psi }{\partial x}=0,\text{ }\theta =1\text{ at }y=0 \qquad \qquad(4)$

$\frac{\partial \psi }{\partial y}=\frac{\partial \psi }{\partial x}=0,\text{ }\theta =0\text{ at }y\to \infty \qquad \qquad(5)$

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:

$\eta =y\cdot H(x) \qquad \qquad(6)$

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

$\psi (x,\eta )=\nu F(\eta )\cdot G(x) \qquad \qquad(7)$

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.,

$\theta (x,\eta )=\theta (\eta ) \qquad \qquad(8)$

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

$\frac{\partial \psi }{\partial y}=\nu {F}'(\eta )\frac{\partial \eta }{\partial y}G(x)=\nu {F}'(\eta )H(x)G(x)$

$\frac{\partial \psi }{\partial x}=\nu {F}'(\eta )y{H}'(x)G(x)+\nu F(\eta )\cdot {G}'(x)$

$\frac{{{\partial }^{2}}\psi }{\partial y\partial x}=\nu {F}'(\eta )[{H}'(x)G(x)+H(x){G}'(x)]+\nu {F}''(\eta )\eta {H}'(x)G(x)$

$\frac{{{\partial }^{2}}\psi }{\partial {{y}^{2}}}=\nu {{F}'}'(\eta )\cdot {{H}^{2}}(x)\cdot G(x)$

$\frac{{{\partial }^{3}}\psi }{\partial {{y}^{3}}}=\nu {F}'''(\eta )\cdot {{H}^{3}}(x)\cdot G(x)$

$\frac{\partial \theta }{\partial x}={\theta }'(\eta )\frac{\partial \eta }{\partial x}={\theta }'(\eta )y{H}'(x)$

$\frac{\partial \theta }{\partial y}={\theta }'(\eta )\frac{\partial \eta }{\partial y}={\theta }'(\eta )H(x)$

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:

${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)$

${\theta }''+\Pr \frac{{{G}'}}{H}F{\theta }'=0 \qquad \qquad(10)$

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:

${{H}^{3}}G=A=\text{const} \qquad \qquad(11)$

$\frac{{{G}'}}{H}=B=\text{const} \qquad \qquad(12)$

$\frac{{H}'G}{{{H}^{2}}}=C=\text{const} \qquad \qquad(13)$

Differentiation of eq. (11) and division of the resultant equation by H4 yields the following equation:

$3\frac{{H}'G}{{{H}^{2}}}+\frac{{{G}'}}{H}=\text{0} \qquad \qquad(14)$

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

${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)$

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

$\frac{g\beta ({{T}_{w}}-{{T}_{\infty }})}{{{\nu }^{2}}{{H}^{3}}G}=1$

$\frac{{{G}'}}{H}=3$

The corresponding G and H functions then become:

$G(x)=4{{\left( \frac{1}{4}\text{G}{{\text{r}}_{x}} \right)}^{1/4}} \qquad \qquad(16)$

$H(x)=\frac{1}{x}{{\left( \frac{1}{4}\text{G}{{\text{r}}_{x}} \right)}^{1/4}} \qquad \qquad(17)$

where Grx is the local Grashof number defined as:

$\text{G}{{\text{r}}_{x}}=\frac{g\beta ({{T}_{w}}-{{T}_{\infty }}){{x}^{3}}}{{{\nu }^{2}}} \qquad \qquad(18)$

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

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

${F}'''+\theta +3F{F}''-2{{({F}')}^{2}}=0 \qquad \qquad(19)$

${\theta }''+3\Pr F{\theta }'=0 \qquad \qquad(20)$

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

$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)$

$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)$

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

$F(\eta )={F}'(\eta )=0\text{ and }\theta (\eta )=1\text{ at }\eta =0 \qquad \qquad(23)$

${F}'(\eta )=\theta (\eta )=0\text{ at }\eta \to \infty \qquad \qquad(24)$

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. ${{h}_{x}}=-\frac{k}{{{T}_{w}}-{{T}_{\infty }}}{{\left( \frac{\partial T}{\partial y} \right)}_{y=0}}$ from Governing Equations for Natural Convection as follows:

${{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)$

The local Nusselt number is

$\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)$

Figure 1: Velocity and temperature profile in the boundary layer for external natural convection over a vertical isothermal surface

where $\phi (\Pr )=-{\theta }'(0)/\sqrt{2}$ 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:

$\phi (\Pr )=\frac{3}{4}{{\left[ \frac{2{{\Pr }^{2}}}{5(1+2{{\Pr }^{1/2}}+2\Pr )} \right]}^{1/4}}$

The local Nusselt number thus becomes:

$\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)$

which can also be rewritten in terms of Rayleigh number

$\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)$

Equations (27) and (28) are valid for $0<\Pr <\infty$.

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

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. (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 2 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. ${{\delta }_{t}}\sim L\text{Ra}_{L}^{-1/4} \qquad \qquad(1)$ and $u\sim \frac{\alpha }{L}\text{Ra}_{L}^{1/2}$ from High Prandtl Number Fluids, which is equivalent to a selection of the stretching functions in the following form

$G(x)=\frac{1}{\Pr }\text{Ra}_{x}^{1/4}$

$H(x)=\frac{1}{x}\text{Ra}_{x}^{1/4}$

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. 1 from Governing Equations for Natural Convection 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.

$\Pr ({{{{F}'''}}_{1}}+{{\theta }_{1}})+{{\theta }_{1}}+{{F}_{1}}{{{{F}''}}_{1}}=\frac{2}{3}{{F}'}_{1}^{2} \qquad \qquad(29)$

${{{\theta }''}_{1}}+{{F}_{1}}{{{\theta }'}_{1}}=0 \qquad \qquad(30)$

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

$\text{N}{{\text{u}}_{x}}=\left\{ \begin{matrix} 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 \\ \end{matrix} \right. \qquad \qquad(31)$

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

${{\overline{\text{Nu}}}_{L}}=\frac{4}{3}\text{N}{{\text{u}}_{L}} \qquad \qquad(32)$

## References

Bejan, A., 2004, Convection Heat Transfer, 3rd ed., John Wiley & Sons, New York.

Ede, A.J., 1964, “Advances in Free Convection,” Advances in Heat Transfer, Vol. 4, pp. 1-64, Academic Press, New York.

Kays, W.M., Crawford, M.E., and Weigand, B., 2005, Convection Heat Transfer, 4th ed., McGraw-Hill, New York, NY.

Le Fevre, E.J., 1956, “Laminar Free Convection from a Vertical Plane Surface,” Proceedings of the 9th International Congress of Applied Mechanics, Brussels, Vol. 4, pp. 168-174.

Ostrach, S., 1953, An Analysis of Laminar Free Convection Flow and Heat Transfer about a Flat Plate Parallel to the Direction of the Generating Body Force, NASA Report, 1111, pp. 63–79.