Analogies and differences in different transport phenomena

From Thermal-FluidsPedia

(Difference between revisions)
Jump to: navigation, search
 
(6 intermediate revisions not shown)
Line 1: Line 1:
-
Most early work in theoretically predicting the heat and/or mass transfer in both laminar and turbulent flow cases was done using the analogy between momentum, heat, and mass and by predicting the approximate results for the heat and/or mass transfer coefficient from the momentum transfer or skin friction coefficient. Clearly there are severe limitations in using this simple approach; however, it is beneficial to understand the advantages and similarities for physical and mathematical modeling as well as the constraints involving this approach. We present this analogy for the classic problem of heat and mass transfer over a flat plate in this section.  Its application to more complicated geometries and boundary conditions as well as turbulent flow is not proven and caution should be taken in applying this approach to other cases.
+
Most early work in theoretically predicting the heat and/or mass transfer in both laminar and turbulent flow cases was done using the analogy between momentum, heat, and mass and by predicting the approximate results for the heat and/or mass transfer coefficient from the momentum transfer or skin friction coefficient. Clearly there are severe limitations in using this simple approach; however, it is beneficial to understand the advantages and similarities for physical and mathematical modeling as well as the constraints involving this approach.<ref>Faghri, A., Zhang, Y., and Howell, J. R., 2010, Advanced Heat and Mass Transfer, Global Digital Press, Columbia, MO.</ref> This analogy is presented here for the classic problem of heat and mass transfer over a flat plate in this section.  Its application to more complicated geometries and boundary conditions as well as turbulent flow is not proven and caution should be taken in applying this approach to other cases.
-
[[Image:Fig4.31.png|thumb|400 px|alt=Mass, momentum, and heat transfer in a laminar boundary layer |Figure 1: Mass, momentum, and heat transfer in a laminar boundary layer.]]
+
[[Image:Fig4.31.png|thumb|400 px|alt=Mass, momentum, and heat transfer in a laminar boundary layer | Mass, momentum, and heat transfer in a laminar boundary layer.]]
-
As presented before, a flat plate at constant wall temperature Tw is exposed to a free stream of constant velocity U∞, temperature T∞ and mass fraction ω1,∞.  Due to binary diffusion it can be presented with the following dimensionless conservation boundary layer equations and boundary conditions corresponding to Fig. 1.
+
As presented before, a flat plate at constant wall temperature Tw is exposed to a free stream of constant velocity ''U<sub>∞</sub>'', temperature ''T<sub>∞</sub>'' and mass fraction ''ω''<sub>1,∞</sub>.  Due to binary diffusion it can be presented with the following dimensionless conservation boundary layer equations and boundary conditions corresponding to the figure to the right.
Continuity
Continuity
 
 
Line 21: Line 21:
<math>\frac{\partial ^{2}u^{+}}{\partial y^{+2}}</math>
<math>\frac{\partial ^{2}u^{+}}{\partial y^{+2}}</math>
</center>
</center>
-
|{{EquationRef|(1)}}
+
|{{EquationRef|(2)}}
|}
|}
Energy
Energy
Line 30: Line 30:
<math>u^{+}\frac{\partial \theta }{\partial x^{+}}+v^{+}\frac{\partial \theta }{\partial y^{+}}=\frac{1}{\Pr \operatorname{Re}}\frac{\partial ^{2}\theta }{\partial y^{+2}}</math>
<math>u^{+}\frac{\partial \theta }{\partial x^{+}}+v^{+}\frac{\partial \theta }{\partial y^{+}}=\frac{1}{\Pr \operatorname{Re}}\frac{\partial ^{2}\theta }{\partial y^{+2}}</math>
</center>
</center>
-
|{{EquationRef|(1)}}
+
|{{EquationRef|(3)}}
|}
|}
Species
Species
Line 39: Line 39:
<math>u^{+}\frac{\partial \varphi }{\partial x^{+}}+v^{+}\frac{\partial \varphi }{\partial y^{+}}=\frac{1}{\text{Sc}\operatorname{Re}}\frac{\partial ^{2}\varphi }{\partial y^{+2}}</math>
<math>u^{+}\frac{\partial \varphi }{\partial x^{+}}+v^{+}\frac{\partial \varphi }{\partial y^{+}}=\frac{1}{\text{Sc}\operatorname{Re}}\frac{\partial ^{2}\varphi }{\partial y^{+2}}</math>
</center>
</center>
-
|{{EquationRef|(1)}}
+
|{{EquationRef|(4)}}
|}
|}
Line 47: Line 47:
<math>\text{at }y^{+}=0,\text{  u}^{\text{+}}=0,\text{  }\theta \text{=}\varphi \text{=0}</math>
<math>\text{at }y^{+}=0,\text{  u}^{\text{+}}=0,\text{  }\theta \text{=}\varphi \text{=0}</math>
</center>
</center>
-
|{{EquationRef|(1)}}
+
|{{EquationRef|(5)}}
|}
|}
Line 55: Line 55:
<math>y^{+}=0,\text{  }v^{+}=v_{w}^{+}</math>
<math>y^{+}=0,\text{  }v^{+}=v_{w}^{+}</math>
</center>
</center>
-
|{{EquationRef|(1)}}
+
|{{EquationRef|(6)}}
|}
|}
Line 63: Line 63:
<math>y^{+}=\infty ,\text{  }u^{+}=1,\text{  }\theta =\varphi =1</math>
<math>y^{+}=\infty ,\text{  }u^{+}=1,\text{  }\theta =\varphi =1</math>
</center>
</center>
-
|{{EquationRef|(1)}}
+
|{{EquationRef|(7)}}
|}
|}
where the dimensionless variables are defined as:
where the dimensionless variables are defined as:
Line 74: Line 74:
\end{align}</math>
\end{align}</math>
</center>
</center>
-
|{{EquationRef|(1)}}
+
|{{EquationRef|(8)}}
|}
|}
-
Let’s first consider the analogy between momentum and heat transfer. Equations (4.349) and (4.350) and appropriate boundary conditions are the same if Pr = 1.  The solutions for u+ and θ are, therefore, exactly the same if Pr = 1 and one expects to have a relationship between the skin friction coefficient, cf, and heat transfer coefficient, h.
+
Let’s first consider the analogy between momentum and heat transfer. Equations (2) and (3) and appropriate boundary conditions are the same if Pr = 1.  The solutions for ''u<sup>+</sup>'' and ''θ'' are, therefore, exactly the same if Pr = 1 and one expects to have a relationship between the skin friction coefficient, ''c<sub>f</sub>'', and heat transfer coefficient, <sub>h</sub>.
{| class="wikitable" border="0"
{| class="wikitable" border="0"
Line 83: Line 83:
<math>\frac{c_{f}}{2}=\frac{\tau _{w}}{\rho U_{\infty }^{2}}=\frac{\mu \left. \frac{\partial u}{\partial y} \right|_{y=0}}{\rho U_{\infty }^{2}}=\frac{\nu \left. \frac{\partial u^{+}}{\partial y^{+}} \right|_{y^{+}=0}}{U_{\infty }L}=\frac{\left. \frac{\partial u^{+}}{\partial y^{+}} \right|_{y^{+}=0}}{\operatorname{Re}}</math>
<math>\frac{c_{f}}{2}=\frac{\tau _{w}}{\rho U_{\infty }^{2}}=\frac{\mu \left. \frac{\partial u}{\partial y} \right|_{y=0}}{\rho U_{\infty }^{2}}=\frac{\nu \left. \frac{\partial u^{+}}{\partial y^{+}} \right|_{y^{+}=0}}{U_{\infty }L}=\frac{\left. \frac{\partial u^{+}}{\partial y^{+}} \right|_{y^{+}=0}}{\operatorname{Re}}</math>
</center>
</center>
-
|{{EquationRef|(1)}}
+
|{{EquationRef|(9)}}
|}
|}
Line 91: Line 91:
<math>Nu=\frac{hL}{k}=\frac{-k\left. \frac{\partial T}{\partial y} \right|_{y=0}L}{k\left( T_{\omega }-T_{\infty } \right)}=\left. \frac{\partial \theta }{\partial y^{+}} \right|_{y^{+}=0}</math>
<math>Nu=\frac{hL}{k}=\frac{-k\left. \frac{\partial T}{\partial y} \right|_{y=0}L}{k\left( T_{\omega }-T_{\infty } \right)}=\left. \frac{\partial \theta }{\partial y^{+}} \right|_{y^{+}=0}</math>
</center>
</center>
-
|{{EquationRef|(1)}}
+
|{{EquationRef|(10)}}
|}
|}
-
Since θ = u+ for Pr = 1, one also concludes that
+
Since ''θ'' = ''u<sub>+</sub>'' for Pr = 1, one also concludes that
{| class="wikitable" border="0"
{| class="wikitable" border="0"
Line 100: Line 100:
<math>\left. \frac{\partial \theta }{\partial y^{+}} \right|_{y^{+}=0}=\left. \frac{\partial u^{+}}{\partial y^{+}} \right|_{y^{+}=0}</math>
<math>\left. \frac{\partial \theta }{\partial y^{+}} \right|_{y^{+}=0}=\left. \frac{\partial u^{+}}{\partial y^{+}} \right|_{y^{+}=0}</math>
</center>
</center>
-
|{{EquationRef|(1)}}
+
|{{EquationRef|(11)}}
|}
|}
-
Therefore, combining eqs. (4.356) and (4.357) and using eq. (4.358) gives
+
Therefore, combining eqs. (9) and (10) and using eq. (11) gives
{| class="wikitable" border="0"
{| class="wikitable" border="0"
Line 109: Line 109:
<math>\frac{c_{f}}{2}=\frac{\text{Nu}}{\operatorname{Re}}</math>
<math>\frac{c_{f}}{2}=\frac{\text{Nu}}{\operatorname{Re}}</math>
</center>
</center>
-
|{{EquationRef|(1)}}
+
|{{EquationRef|(12)}}
|}
|}
-
This relationship between the friction coefficient and Nu is referred to as Reynolds analogy and is appropriate when Pr = 1. If Pr ≠ 1, we already concluded that <math>\frac{\delta _{T}}{\delta }=\Pr ^{-1/3}\text{  for  }0.5\le \Pr \le 10</math> from the similarity solution presented in Section 4.6. Using this information, one can generalize the result of the Reynolds analogy to Pr ≠ 1 by
+
This relationship between the friction coefficient and Nu is referred to as Reynolds analogy and is appropriate when Pr = 1. If Pr ≠ 1, we already concluded that <math>\frac{\delta _{T}}{\delta }=\Pr ^{-1/3}\text{  for  }0.5\le \Pr \le 10</math> from the [[similarity solutions|similarity solution]]. Using this information, one can generalize the result of the Reynolds analogy to Pr ≠ 1 by
{| class="wikitable" border="0"
{| class="wikitable" border="0"
Line 118: Line 118:
<math>\frac{c_{f}}{2}=\frac{\text{Nu}}{\operatorname{Re}}\Pr ^{-1/3}</math>
<math>\frac{c_{f}}{2}=\frac{\text{Nu}}{\operatorname{Re}}\Pr ^{-1/3}</math>
</center>
</center>
-
|{{EquationRef|(1)}}
+
|{{EquationRef|(13)}}
|}
|}
-
Now, we focus on similarities between the heat and mass transfer for comparison of eqs. (4.350) and (4.351) with their appropriate boundary conditions. It is clear that the solution for differential equations (4.350) and (4.351) for θ and <math>\varphi </math> are the same if Sc and Pr are interchanged appropriately.
+
Now, we focus on similarities between the heat and mass transfer for comparison of eqs. (3) and (4) with their appropriate boundary conditions. It is clear that the solution for differential equations (3) and (4) for ''θ'' and <math>\varphi </math> are the same if Sc and Pr are interchanged appropriately.
-
We already know the solution for eq. (4.350) for <math>v_{w}=0</math> from the similarity solution for 0.5 ≤ Pr ≤ 10
+
We already know the solution for eq. (3) for <math>v_{w}=0</math> from the similarity solution for 0.5 ≤ Pr ≤ 10
{| class="wikitable" border="0"
{| class="wikitable" border="0"
Line 128: Line 128:
<math>\text{Nu}_{x}=0.332\operatorname{Re}_{x}^{1/2}\Pr ^{1/3}</math>
<math>\text{Nu}_{x}=0.332\operatorname{Re}_{x}^{1/2}\Pr ^{1/3}</math>
</center>
</center>
-
|{{EquationRef|(1)}}
+
|{{EquationRef|(14)}}
|}
|}
-
Therefore, it can also be assumed that the solution of eq. (4.351) for <math>v_{w}=0</math> is
+
Therefore, it can also be assumed that the solution of eq. (4) for <math>v_{w}=0</math> is
{| class="wikitable" border="0"
{| class="wikitable" border="0"
Line 137: Line 137:
<math>\text{Sh}_{x}=0.332\operatorname{Re}_{x}^{1/2}\text{Sc}^{1/3}</math>
<math>\text{Sh}_{x}=0.332\operatorname{Re}_{x}^{1/2}\text{Sc}^{1/3}</math>
</center>
</center>
-
|{{EquationRef|(1)}}
+
|{{EquationRef|(15)}}
|}
|}
-
Combining eqs. (4.361) and (4.362) gives
+
Combining eqs. (14) and (15) gives
{| class="wikitable" border="0"
{| class="wikitable" border="0"
Line 146: Line 146:
<math>\frac{\text{Nu}}{\text{Sh}}\text{=}\frac{\text{Pr}^{\text{1/3}}}{\text{Sc}^{\text{1/3}}}</math>
<math>\frac{\text{Nu}}{\text{Sh}}\text{=}\frac{\text{Pr}^{\text{1/3}}}{\text{Sc}^{\text{1/3}}}</math>
</center>
</center>
-
|{{EquationRef|(1)}}
+
|{{EquationRef|(16)}}
|}
|}
-
It should be noted that the convective effect due to vertical velocity at the surface was neglected in predicting h and hm in the above analysis and therefore the analogy presented in (4.363) is for very low mass transfer at the wall. The effect is that the contribution of <math>v_{w}</math> in heat and mass flux was primarily due to diffusion.
+
It should be noted that the convective effect due to vertical velocity at the surface was neglected in predicting ''h'' and ''h''<sub>m</sub> in the above analysis and therefore the analogy presented in (16) is for very low mass transfer at the wall. The effect is that the contribution of <math>v_{w}</math> in heat and mass flux was primarily due to diffusion.
{| class="wikitable" border="0"
{| class="wikitable" border="0"
Line 155: Line 155:
<math>{q}''_{w}=\rho c_{p}v_{w}\left( T_{w}-T_{\infty } \right)-k\left. \frac{\partial T}{\partial y} \right|_{y=0}</math>
<math>{q}''_{w}=\rho c_{p}v_{w}\left( T_{w}-T_{\infty } \right)-k\left. \frac{\partial T}{\partial y} \right|_{y=0}</math>
</center>
</center>
-
|{{EquationRef|(1)}}
+
|{{EquationRef|(17)}}
|}
|}
Line 163: Line 163:
<math>{\dot{m}}''_{w}=\rho \left( \omega _{1,\omega }v_{w}-D_{12}\left. \frac{\partial \omega _{1}}{\partial y} \right|_{y=0} \right)</math>
<math>{\dot{m}}''_{w}=\rho \left( \omega _{1,\omega }v_{w}-D_{12}\left. \frac{\partial \omega _{1}}{\partial y} \right|_{y=0} \right)</math>
</center>
</center>
-
|{{EquationRef|(1)}}
+
|{{EquationRef|(18)}}
|}
|}
-
The analogy in momentum, heat, and mass transfer can also be applied to complex and/or coupled transport phenomena problems including phase change and chemical reactions. Obviously, in these circumstances, the simple relationship developed in eq. (4.363) is not applicable. To show the usefulness of this analogy to more complex transport phenomena problems we will apply it to sublimation with a chemical reaction for forced convective flow over a flat plate.
+
The analogy in momentum, heat, and mass transfer can also be applied to complex and/or coupled transport phenomena problems including phase change and chemical reactions. Obviously, in these circumstances, the simple relationship developed in eq. (16) is not applicable. To show the usefulness of this analogy to more complex transport phenomena problems we will apply it to sublimation with a chemical reaction for forced convective flow over a flat plate.
 +
 
 +
[[Image:Fig4.32.png|thumb|400 px|alt=Sublimation with chemical reaction | Sublimation with chemical reaction. <ref name="K2001">Kaviany, M., 2001, Principles of Convective Heat Transfer, 2nd ed., Springer-Verlag, New York.</ref>]]
 +
 
During combustion involving a solid fuel, the solid fuel may either burn directly or be sublimated before combustion. In the latter case – which will be discussed in this subsection – gaseous fuel diffuses away from the solid-vapor interface. Meanwhile, the gaseous oxidant diffuses toward the solid-vapor interface. Under the right conditions, the mass flux of vapor fuel and the gaseous oxidant meet and the chemical reaction occurs at a certain zone known as the flame. The flame is usually a very thin region with a color dictated by the temperature of combustion.
During combustion involving a solid fuel, the solid fuel may either burn directly or be sublimated before combustion. In the latter case – which will be discussed in this subsection – gaseous fuel diffuses away from the solid-vapor interface. Meanwhile, the gaseous oxidant diffuses toward the solid-vapor interface. Under the right conditions, the mass flux of vapor fuel and the gaseous oxidant meet and the chemical reaction occurs at a certain zone known as the flame. The flame is usually a very thin region with a color dictated by the temperature of combustion.
-
Figure 4.32 shows the physical model of the problem under consideration.  The concentration of the fuel is highest at the solid fuel surface, and decreases as the location of the flame is approached.  The gaseous fuel diffuses away from the solid fuel surface and meets the oxidant as it flows parallel to the solid fuel surface. Combustion occurs in a thin reaction zone where the temperature is the highest, and the latent heat of sublimation is supplied by combustion. The combustion of solid fuel through sublimation can be modeled as a steady-state boundary layer type flow with sublimation and chemical reaction.
 
-
To model the problem, the following assumptions are made:
 
-
1. The fuel is supplied by sublimation at a steady rate.
 
-
2. The Lewis number is unity, so the thermal and concentration boundary layers have the same thickness.
 
-
3. The buoyancy force is negligible.
 
-
[[Image:Fig4.32.png|thumb|400 px|alt=Sublimation with chemical reaction (Kaviany, 2001) |Figure 2: Sublimation with chemical reaction (Kaviany, 2001).]]
+
Figure to the right shows the physical model of the problem under consideration. The concentration of the fuel is highest at the solid fuel surface, and decreases as the location of the flame is approached. The gaseous fuel diffuses away from the solid fuel surface and meets the oxidant as it flows parallel to the solid fuel surface. Combustion occurs in a thin reaction zone where the temperature is the highest, and the latent heat of sublimation is supplied by combustion. The combustion of solid fuel through sublimation can be modeled as a steady-state boundary layer type flow with sublimation and chemical reaction.
 +
 
 +
To model the problem, the following assumptions are made:<br>
 +
1. The fuel is supplied by sublimation at a steady rate.<br>
 +
2. The Lewis number is unity, so the thermal and concentration boundary layers have the same thickness.<br>
 +
3. The buoyancy force is negligible.<br>
The conservation equations for mass, momentum, energy, and species in the boundary layer are  
The conservation equations for mass, momentum, energy, and species in the boundary layer are  
Line 182: Line 185:
<math>\frac{\partial (\rho u)}{\partial x}+\frac{\partial (\rho v)}{\partial y}=0</math>
<math>\frac{\partial (\rho u)}{\partial x}+\frac{\partial (\rho v)}{\partial y}=0</math>
</center>
</center>
-
|{{EquationRef|(1)}}
+
|{{EquationRef|(19)}}
|}
|}
Line 190: Line 193:
<math>u\frac{\partial u}{\partial x}+v\frac{\partial u}{\partial y}=\frac{\partial }{\partial y}\left( \nu \frac{\partial u}{\partial y} \right)</math>
<math>u\frac{\partial u}{\partial x}+v\frac{\partial u}{\partial y}=\frac{\partial }{\partial y}\left( \nu \frac{\partial u}{\partial y} \right)</math>
</center>
</center>
-
|{{EquationRef|(1)}}
+
|{{EquationRef|(20)}}
|}
|}
Line 198: Line 201:
<math>\frac{\partial }{\partial x}(\rho c_{p}uT)+\frac{\partial }{\partial y}(\rho c_{p}vT)=\frac{\partial }{\partial y}\left( k\frac{\partial T}{\partial y} \right)+{\dot{m}}'''_{o}\text{ }h_{c,o}</math>
<math>\frac{\partial }{\partial x}(\rho c_{p}uT)+\frac{\partial }{\partial y}(\rho c_{p}vT)=\frac{\partial }{\partial y}\left( k\frac{\partial T}{\partial y} \right)+{\dot{m}}'''_{o}\text{ }h_{c,o}</math>
</center>
</center>
-
|{{EquationRef|(1)}}
+
|{{EquationRef|(21)}}
|}
|}
-
 
{| class="wikitable" border="0"
{| class="wikitable" border="0"
Line 207: Line 209:
<math>\frac{\partial }{\partial x}(\rho u\omega _{o})+\frac{\partial }{\partial y}(\rho v\omega _{o})=\frac{\partial }{\partial y}\left( \rho D\frac{\partial \omega _{o}}{\partial y} \right)-{\dot{m}}'''_{o}</math>
<math>\frac{\partial }{\partial x}(\rho u\omega _{o})+\frac{\partial }{\partial y}(\rho v\omega _{o})=\frac{\partial }{\partial y}\left( \rho D\frac{\partial \omega _{o}}{\partial y} \right)-{\dot{m}}'''_{o}</math>
</center>
</center>
-
|{{EquationRef|(1)}}
+
|{{EquationRef|(22)}}
|}
|}
-
where <math>{\dot{m}}'''_{0}</math> is the rate of oxidant consumption (kg/m3*s), hc,0 is the heat released by combustion per unit mass consumption of the oxidant (J/kg), which is different from the combustion heat, and<math>\omega _{0}</math> is the mass fraction of the oxidant in the gaseous mixture.
+
where <math>{\dot{m}}'''_{0}</math> is the rate of oxidant consumption (kg/m<sup>3</sup>*s), ''h''<sub>c,0</sub> is the heat released by combustion per unit mass consumption of the oxidant (J/kg), which is different from the combustion heat, and<math>\omega _{0}</math> is the mass fraction of the oxidant in the gaseous mixture.
-
The corresponding boundary conditions of eqs. (4.366) to (4.369) are  
+
The corresponding boundary conditions of eqs. (19) to (22) are  
{| class="wikitable" border="0"
{| class="wikitable" border="0"
|-
|-
Line 220: Line 222:
\end{matrix}\_at\_y\to \infty </math>
\end{matrix}\_at\_y\to \infty </math>
</center>
</center>
-
| {{EquationRef|(1)}}
+
| {{EquationRef|(23)}}
|}
|}
{| class="wikitable" border="0"
{| class="wikitable" border="0"
Line 231: Line 233:
\end{matrix}\_at\_y=0</math>
\end{matrix}\_at\_y=0</math>
</center>
</center>
-
| {{EquationRef|(1)}}
+
| {{EquationRef|(24)}}
|}
|}
-
where <math>{\dot{m}}''_{f}</math> is the rate of solid fuel sublimation per unit area (kg/m2?s) and ρ is the density of the mixture.
+
where <math>{\dot{m}}''_{f}</math> is the rate of solid fuel sublimation per unit area (kg/m<sup>2</sup>*s) and ''ρ'' is the density of the mixture.
 +
 
The shear stress at the solid fuel surface is
The shear stress at the solid fuel surface is
{| class="wikitable" border="0"
{| class="wikitable" border="0"
Line 243: Line 246:
\end{matrix}</math>
\end{matrix}</math>
</center>
</center>
-
| {{EquationRef|(1)}}
+
| {{EquationRef|(25)}}
-
|}
+
|}
The heat flux at the solid fuel surface is  
The heat flux at the solid fuel surface is  
Line 254: Line 257:
\end{matrix}</math>
\end{matrix}</math>
</center>
</center>
-
| {{EquationRef|(1)}}
+
| {{EquationRef|(26)}}
|}
|}
-
The exact solution of the heat and mass problem described by eqs. (4.366) to (4.369) can be obtained using numerical simulation.  It is useful here to introduce the results obtained by Kaviany (2001) using an analogy between momentum, heat, and mass transfer. Multiplying eq. (4.369) by hc,0 and adding the result to eq. (4.368), one obtains  
+
The exact solution of the heat and mass problem described by eqs. (19) to (22) can be obtained using numerical simulation.  It is useful here to introduce the results obtained by Kaviany <ref name="K2001"/> using an analogy between momentum, heat, and mass transfer. Multiplying eq. (22) by ''h''<sub>c,0</sub> and adding the result to eq. (21), one obtains  
{| class="wikitable" border="0"
{| class="wikitable" border="0"
Line 264: Line 267:
<math>\frac{\partial }{\partial x}\left[ \rho u(c_{p}T+\omega _{o}h_{c,o}) \right]+\frac{\partial }{\partial y}\left[ \rho v(c_{p}T+\omega _{o}h_{c,o}) \right]=\frac{\partial }{\partial y}\left[ k\frac{\partial T}{\partial y}+\rho Dh_{c,o}\frac{\partial \omega _{o}}{\partial y} \right]</math>
<math>\frac{\partial }{\partial x}\left[ \rho u(c_{p}T+\omega _{o}h_{c,o}) \right]+\frac{\partial }{\partial y}\left[ \rho v(c_{p}T+\omega _{o}h_{c,o}) \right]=\frac{\partial }{\partial y}\left[ k\frac{\partial T}{\partial y}+\rho Dh_{c,o}\frac{\partial \omega _{o}}{\partial y} \right]</math>
</center>
</center>
-
|{{EquationRef|(1)}}
+
|{{EquationRef|(27)}}
|}
|}
-
Considering the assumption that the Lewis number is unity, i.e. <math>Le=\nu /D=1</math>, eq. (4.374) can be rewritten as
+
Considering the assumption that the Lewis number is unity, i.e. <math>Le=\nu /D=1</math>, eq. (27) can be rewritten as
{| class="wikitable" border="0"
{| class="wikitable" border="0"
|-
|-
Line 275: Line 278:
\end{align}</math>
\end{align}</math>
</center>
</center>
-
| {{EquationRef|(1)}}
+
| {{EquationRef|(28)}}
|}
|}
which can be viewed as an energy equation with quantity <math>c_{p}T+\omega _{0}h_{c,0}</math>) as a dependent variable.  
which can be viewed as an energy equation with quantity <math>c_{p}T+\omega _{0}h_{c,0}</math>) as a dependent variable.  
-
Since <math>\partial \omega _{o}/\partial y=0</math> at y = 0, i.e., the solid fuel surface is not permeable for the oxidant, eq(4.373) can be rewritten as
+
Since <math>\partial \omega _{o}/\partial y=0</math> at ''y'' = 0, i.e., the solid fuel surface is not permeable for the oxidant, eq. (26) can be rewritten as
{| class="wikitable" border="0"
{| class="wikitable" border="0"
|-
|-
Line 287: Line 290:
\end{matrix}</math>
\end{matrix}</math>
</center>
</center>
-
| {{EquationRef|(1)}}
+
| {{EquationRef|(29)}}
|}
|}
Line 294: Line 297:
|-
|-
| width="100%" |<center>
| width="100%" |<center>
-
<math>\begin{align}
+
<math>{q}''_{w}=\frac{\tau _{w}}{U_{\infty }}\left[ (c_{p}T+\omega _{o}h_{c,o})_{w}-(c_{p}T+\omega _{o}h_{c,o})_{\infty } \right]=\frac{\tau _{w}}{U_{\infty }}[c_{p}(T_{w}-T_{\infty })+h_{c,o}(\omega _{o,w}-\omega _{o,\infty})]</math>
-
  & {q}''_{w}=\frac{\tau _{w}}{U_{\infty }}\left[ (c_{p}T+\omega _{o}h_{c,o})_{w}-(c_{p}T+\omega _{o}h_{c,o})_{\infty } \right] \\
+
-
& \begin{matrix}
+
-
  {} & =\frac{\tau _{w}}{U_{\infty }}\left[ c_{p}(T_{w}-T_{\infty })+h_{c,o}(\omega _{o,}_{w}-\omega _{o}_{,\infty }) \right] \\
+
-
\end{matrix} \\
+
-
\end{align}</math>
+
</center>
</center>
-
| {{EquationRef|(1)}}
+
| {{EquationRef|(30)}}
|}
|}
Line 311: Line 309:
<math>-{q}''_{w}={\dot{m}}''_{f}h_{sv}+{q}''_{\ell }</math>
<math>-{q}''_{w}={\dot{m}}''_{f}h_{sv}+{q}''_{\ell }</math>
</center>
</center>
-
|{{EquationRef|(1)}}
+
|{{EquationRef|(31)}}
|}
|}
-
where the two terms on the right-hand side of eq. (4.378) represent the latent heat of sublimation, and the sensible heat required to raise the surface temperature of the solid fuel to sublimation temperature and heat loss to the solid fuel.  
+
where the two terms on the right-hand side of eq. (31) represent the latent heat of sublimation, and the sensible heat required to raise the surface temperature of the solid fuel to sublimation temperature and heat loss to the solid fuel.  
-
Combining eqs. (4.377) and (4.378) yields the rate of sublimation on the solid fuel surface
+
Combining eqs. (30) and (31) yields the rate of sublimation on the solid fuel surface
{| class="wikitable" border="0"
{| class="wikitable" border="0"
Line 321: Line 319:
<math>{\dot{m}}''_{f}=Z\frac{\tau _{w}}{U_{\infty }}</math>
<math>{\dot{m}}''_{f}=Z\frac{\tau _{w}}{U_{\infty }}</math>
</center>
</center>
-
|{{EquationRef|(1)}}
+
|{{EquationRef|(32)}}
|}
|}
-
where Z is the transfer driving force or transfer number defined as
+
where ''Z'' is the transfer driving force or transfer number defined as
{| class="wikitable" border="0"
{| class="wikitable" border="0"
Line 330: Line 328:
<math>Z=\frac{c_{p}(T_{\infty }-T_{w})+h_{c,o}(\omega _{o,\infty }-\omega _{o,w})}{h_{sv}+{q}''_{\ell }/\dot{{m}''}_{f}}</math>
<math>Z=\frac{c_{p}(T_{\infty }-T_{w})+h_{c,o}(\omega _{o,\infty }-\omega _{o,w})}{h_{sv}+{q}''_{\ell }/\dot{{m}''}_{f}}</math>
</center>
</center>
-
|{{EquationRef|(1)}}
+
|{{EquationRef|(33)}}
|}
|}
Using the friction coefficient gives
Using the friction coefficient gives
Line 339: Line 337:
<math>c_{f}=\frac{\tau _{w}}{\rho U_{\infty }^{2}/2}</math>
<math>c_{f}=\frac{\tau _{w}}{\rho U_{\infty }^{2}/2}</math>
</center>
</center>
-
|{{EquationRef|(1)}}
+
|{{EquationRef|(34)}}
|}
|}
-
eq. (4.379) therefore becomes
+
eq. (32) therefore becomes
{| class="wikitable" border="0"
{| class="wikitable" border="0"
Line 348: Line 346:
<math>{\dot{m}}''_{f}=\frac{c_{f}}{2}\rho U_{\infty }Z</math>
<math>{\dot{m}}''_{f}=\frac{c_{f}}{2}\rho U_{\infty }Z</math>
</center>
</center>
-
|{{EquationRef|(1)}}
+
|{{EquationRef|(35)}}
|}
|}
The surface blowing velocity of the gaseous fuel is then
The surface blowing velocity of the gaseous fuel is then
Line 357: Line 355:
<math>v_{w}=\frac{{\dot{m}}''_{f}}{\rho }=\frac{c_{f}}{2}U_{\infty }Z</math>
<math>v_{w}=\frac{{\dot{m}}''_{f}}{\rho }=\frac{c_{f}}{2}U_{\infty }Z</math>
</center>
</center>
-
|{{EquationRef|(1)}}
+
|{{EquationRef|(36)}}
|}
|}
-
where the friction coefficient, cf, can be obtained from the solution of the boundary layer flow over a flat plate with blowing on the surface.  The similarity solution of the boundary layer flow problem exists only if the blowing velocity satisfies  
+
where the friction coefficient, ''c<sub>f</sub>'', can be obtained from the solution of the boundary layer flow over a flat plate with blowing on the surface.  The similarity solution of the boundary layer flow problem exists only if the blowing velocity satisfies  
-
<math>v_{w}\propto x^{-1/2}</math>
+
<math>v_{w}\propto x^{-1/2}</math>. In this case, one can define a blowing parameter as
-
. In this case, one can define a blowing parameter as
+
{| class="wikitable" border="0"
{| class="wikitable" border="0"
Line 368: Line 365:
<math>B=\frac{(\rho v)_{w}}{(\rho u)_{\infty }}\operatorname{Re}_{x}^{1/2}</math>
<math>B=\frac{(\rho v)_{w}}{(\rho u)_{\infty }}\operatorname{Re}_{x}^{1/2}</math>
</center>
</center>
-
|{{EquationRef|(1)}}
+
|{{EquationRef|(37)}}
|}
|}
-
Combining eqs. (4.383) and (4.384) yields
+
Combining eqs. (36) and (37) yields
{| class="wikitable" border="0"
{| class="wikitable" border="0"
Line 377: Line 374:
<math>B=\frac{Z}{2}\operatorname{Re}_{x}^{1/2}c_{f}</math>
<math>B=\frac{Z}{2}\operatorname{Re}_{x}^{1/2}c_{f}</math>
</center>
</center>
-
|{{EquationRef|(1)}}
+
|{{EquationRef|(38)}}
|}
|}
-
Glassman (1987) recommended an empirical form of eq. (4.385) based on numerical and experimental results:
+
Glassman <ref name="G1987">Glassman, I.,  1987, Combustion, 2nd ed. Academic Press, Orlando, FL</ref> recommended an empirical form of eq. (38) based on numerical and experimental results:
{| class="wikitable" border="0"
{| class="wikitable" border="0"
Line 386: Line 383:
<math>B=\frac{\ln (1+Z)}{2.6Z^{0.15}}</math>
<math>B=\frac{\ln (1+Z)}{2.6Z^{0.15}}</math>
</center>
</center>
-
|{{EquationRef|(1)}}
+
|{{EquationRef|(39)}}
|}
|}
 +
 +
==References=
 +
{{Reflist}}

Current revision as of 19:49, 23 July 2010

Most early work in theoretically predicting the heat and/or mass transfer in both laminar and turbulent flow cases was done using the analogy between momentum, heat, and mass and by predicting the approximate results for the heat and/or mass transfer coefficient from the momentum transfer or skin friction coefficient. Clearly there are severe limitations in using this simple approach; however, it is beneficial to understand the advantages and similarities for physical and mathematical modeling as well as the constraints involving this approach.[1] This analogy is presented here for the classic problem of heat and mass transfer over a flat plate in this section. Its application to more complicated geometries and boundary conditions as well as turbulent flow is not proven and caution should be taken in applying this approach to other cases.

Mass, momentum, and heat transfer in a laminar boundary layer
Mass, momentum, and heat transfer in a laminar boundary layer.

As presented before, a flat plate at constant wall temperature Tw is exposed to a free stream of constant velocity U, temperature T and mass fraction ω1,∞. Due to binary diffusion it can be presented with the following dimensionless conservation boundary layer equations and boundary conditions corresponding to the figure to the right. Continuity

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

(1)

Momentum

u^{+}\frac{\partial u^{+}}{\partial x^{+}}+v^{+}\frac{\partial u^{+}}{\partial y^{+}}=\frac{1}{\operatorname{Re}} \frac{\partial ^{2}u^{+}}{\partial y^{+2}}

(2)

Energy

u^{+}\frac{\partial \theta }{\partial x^{+}}+v^{+}\frac{\partial \theta }{\partial y^{+}}=\frac{1}{\Pr \operatorname{Re}}\frac{\partial ^{2}\theta }{\partial y^{+2}}

(3)

Species

u^{+}\frac{\partial \varphi }{\partial x^{+}}+v^{+}\frac{\partial \varphi }{\partial y^{+}}=\frac{1}{\text{Sc}\operatorname{Re}}\frac{\partial ^{2}\varphi }{\partial y^{+2}}

(4)

\text{at }y^{+}=0,\text{  u}^{\text{+}}=0,\text{  }\theta \text{=}\varphi \text{=0}

(5)

y^{+}=0,\text{   }v^{+}=v_{w}^{+}

(6)

y^{+}=\infty ,\text{   }u^{+}=1,\text{   }\theta =\varphi =1

(7)

where the dimensionless variables are defined as:

\begin{align}
  & u^{+}=\frac{u}{U_{\infty }},\text{   }\theta =\frac{T-T_{w}}{T_{\infty }-T_{w}}\text{,   }\varphi =\frac{\omega _{1}-\omega _{1,w}}{\omega _{1,\infty }-\omega _{1,w}}, \\ 
 & v^{+}=\frac{v}{U_{\infty }},\text{   }y^{+}=\frac{y}{L},\text{   }x^{+}=\frac{x}{L},\text{   }\operatorname{Re}=\frac{U_{\infty }L}{\nu } \\ 
\end{align}

(8)

Let’s first consider the analogy between momentum and heat transfer. Equations (2) and (3) and appropriate boundary conditions are the same if Pr = 1. The solutions for u+ and θ are, therefore, exactly the same if Pr = 1 and one expects to have a relationship between the skin friction coefficient, cf, and heat transfer coefficient, h.

\frac{c_{f}}{2}=\frac{\tau _{w}}{\rho U_{\infty }^{2}}=\frac{\mu \left. \frac{\partial u}{\partial y} \right|_{y=0}}{\rho U_{\infty }^{2}}=\frac{\nu \left. \frac{\partial u^{+}}{\partial y^{+}} \right|_{y^{+}=0}}{U_{\infty }L}=\frac{\left. \frac{\partial u^{+}}{\partial y^{+}} \right|_{y^{+}=0}}{\operatorname{Re}}

(9)

Nu=\frac{hL}{k}=\frac{-k\left. \frac{\partial T}{\partial y} \right|_{y=0}L}{k\left( T_{\omega }-T_{\infty } \right)}=\left. \frac{\partial \theta }{\partial y^{+}} \right|_{y^{+}=0}

(10)

Since θ = u+ for Pr = 1, one also concludes that

\left. \frac{\partial \theta }{\partial y^{+}} \right|_{y^{+}=0}=\left. \frac{\partial u^{+}}{\partial y^{+}} \right|_{y^{+}=0}

(11)

Therefore, combining eqs. (9) and (10) and using eq. (11) gives

\frac{c_{f}}{2}=\frac{\text{Nu}}{\operatorname{Re}}

(12)

This relationship between the friction coefficient and Nu is referred to as Reynolds analogy and is appropriate when Pr = 1. If Pr ≠ 1, we already concluded that \frac{\delta _{T}}{\delta }=\Pr ^{-1/3}\text{  for  }0.5\le \Pr \le 10 from the similarity solution. Using this information, one can generalize the result of the Reynolds analogy to Pr ≠ 1 by

\frac{c_{f}}{2}=\frac{\text{Nu}}{\operatorname{Re}}\Pr ^{-1/3}

(13)

Now, we focus on similarities between the heat and mass transfer for comparison of eqs. (3) and (4) with their appropriate boundary conditions. It is clear that the solution for differential equations (3) and (4) for θ and \varphi are the same if Sc and Pr are interchanged appropriately. We already know the solution for eq. (3) for vw = 0 from the similarity solution for 0.5 ≤ Pr ≤ 10

\text{Nu}_{x}=0.332\operatorname{Re}_{x}^{1/2}\Pr ^{1/3}

(14)

Therefore, it can also be assumed that the solution of eq. (4) for vw = 0 is

\text{Sh}_{x}=0.332\operatorname{Re}_{x}^{1/2}\text{Sc}^{1/3}

(15)

Combining eqs. (14) and (15) gives

\frac{\text{Nu}}{\text{Sh}}\text{=}\frac{\text{Pr}^{\text{1/3}}}{\text{Sc}^{\text{1/3}}}

(16)

It should be noted that the convective effect due to vertical velocity at the surface was neglected in predicting h and hm in the above analysis and therefore the analogy presented in (16) is for very low mass transfer at the wall. The effect is that the contribution of vw in heat and mass flux was primarily due to diffusion.

{q}''_{w}=\rho c_{p}v_{w}\left( T_{w}-T_{\infty } \right)-k\left. \frac{\partial T}{\partial y} \right|_{y=0}

(17)

{\dot{m}}''_{w}=\rho \left( \omega _{1,\omega }v_{w}-D_{12}\left. \frac{\partial \omega _{1}}{\partial y} \right|_{y=0} \right)

(18)

The analogy in momentum, heat, and mass transfer can also be applied to complex and/or coupled transport phenomena problems including phase change and chemical reactions. Obviously, in these circumstances, the simple relationship developed in eq. (16) is not applicable. To show the usefulness of this analogy to more complex transport phenomena problems we will apply it to sublimation with a chemical reaction for forced convective flow over a flat plate.

Sublimation with chemical reaction
Sublimation with chemical reaction. [2]

During combustion involving a solid fuel, the solid fuel may either burn directly or be sublimated before combustion. In the latter case – which will be discussed in this subsection – gaseous fuel diffuses away from the solid-vapor interface. Meanwhile, the gaseous oxidant diffuses toward the solid-vapor interface. Under the right conditions, the mass flux of vapor fuel and the gaseous oxidant meet and the chemical reaction occurs at a certain zone known as the flame. The flame is usually a very thin region with a color dictated by the temperature of combustion.

Figure to the right shows the physical model of the problem under consideration. The concentration of the fuel is highest at the solid fuel surface, and decreases as the location of the flame is approached. The gaseous fuel diffuses away from the solid fuel surface and meets the oxidant as it flows parallel to the solid fuel surface. Combustion occurs in a thin reaction zone where the temperature is the highest, and the latent heat of sublimation is supplied by combustion. The combustion of solid fuel through sublimation can be modeled as a steady-state boundary layer type flow with sublimation and chemical reaction.

To model the problem, the following assumptions are made:
1. The fuel is supplied by sublimation at a steady rate.
2. The Lewis number is unity, so the thermal and concentration boundary layers have the same thickness.
3. The buoyancy force is negligible.

The conservation equations for mass, momentum, energy, and species in the boundary layer are

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

(19)

u\frac{\partial u}{\partial x}+v\frac{\partial u}{\partial y}=\frac{\partial }{\partial y}\left( \nu \frac{\partial u}{\partial y} \right)

(20)

\frac{\partial }{\partial x}(\rho c_{p}uT)+\frac{\partial }{\partial y}(\rho c_{p}vT)=\frac{\partial }{\partial y}\left( k\frac{\partial T}{\partial y} \right)+{\dot{m}}'''_{o}\text{ }h_{c,o}

(21)

\frac{\partial }{\partial x}(\rho u\omega _{o})+\frac{\partial }{\partial y}(\rho v\omega _{o})=\frac{\partial }{\partial y}\left( \rho D\frac{\partial \omega _{o}}{\partial y} \right)-{\dot{m}}'''_{o}

(22)

where {\dot{m}}'''_{0} is the rate of oxidant consumption (kg/m3*s), hc,0 is the heat released by combustion per unit mass consumption of the oxidant (J/kg), which is different from the combustion heat, andω0 is the mass fraction of the oxidant in the gaseous mixture. The corresponding boundary conditions of eqs. (19) to (22) are

u\to U_{\infty }\begin{matrix}
   , & T\to T_{\infty }\begin{matrix}
   , & \omega _{o}\to \omega _{o,\infty }  \\
\end{matrix}  \\
\end{matrix}\_at\_y\to \infty

(23)

u=0\begin{matrix}
   , & v=\frac{{\dot{m}}''_{f}}{\rho }\begin{matrix}
   , & \frac{\partial \omega _{o}}{\partial y}=0  \\
\end{matrix}  \\
\end{matrix}\_at\_y=0

(24)

where {\dot{m}}''_{f} is the rate of solid fuel sublimation per unit area (kg/m2*s) and ρ is the density of the mixture.

The shear stress at the solid fuel surface is

\tau _{w}=\mu \frac{\partial u}{\partial y}\begin{matrix}
   , & y=0  \\
\end{matrix}

(25)

The heat flux at the solid fuel surface is

{q}''_{w}=-k\frac{\partial T}{\partial y}\begin{matrix}
   , & y=0  \\
\end{matrix}

(26)

The exact solution of the heat and mass problem described by eqs. (19) to (22) can be obtained using numerical simulation. It is useful here to introduce the results obtained by Kaviany [2] using an analogy between momentum, heat, and mass transfer. Multiplying eq. (22) by hc,0 and adding the result to eq. (21), one obtains

\frac{\partial }{\partial x}\left[ \rho u(c_{p}T+\omega _{o}h_{c,o}) \right]+\frac{\partial }{\partial y}\left[ \rho v(c_{p}T+\omega _{o}h_{c,o}) \right]=\frac{\partial }{\partial y}\left[ k\frac{\partial T}{\partial y}+\rho Dh_{c,o}\frac{\partial \omega _{o}}{\partial y} \right]

(27)

Considering the assumption that the Lewis number is unity, i.e. Le = ν / D = 1, eq. (27) can be rewritten as

\begin{align}
  & \frac{\partial }{\partial x}\left[ \rho u(c_{p}T+\omega _{o}h_{c,o}) \right]+\frac{\partial }{\partial y}\left[ \rho v(c_{p}T+\omega _{o}h_{c,o}) \right] \\ 
 & =\frac{\partial }{\partial y}\left[ \rho \alpha \frac{\partial }{\partial y}(c_{p}T+\omega _{o}h_{c,o}) \right] \\ 
\end{align}

(28)

which can be viewed as an energy equation with quantity cpT + ω0hc,0) as a dependent variable. Since \partial \omega _{o}/\partial y=0 at y = 0, i.e., the solid fuel surface is not permeable for the oxidant, eq. (26) can be rewritten as

{q}''_{w}=-\rho \alpha \frac{\partial }{\partial y}(c_{p}T+\omega _{o}h_{c,o})\begin{matrix}
   , & y=0  \\
\end{matrix}

(29)

Analogy between surface shear stress and the surface energy flux yields

{q}''_{w}=\frac{\tau _{w}}{U_{\infty }}\left[ (c_{p}T+\omega _{o}h_{c,o})_{w}-(c_{p}T+\omega _{o}h_{c,o})_{\infty } \right]=\frac{\tau _{w}}{U_{\infty }}[c_{p}(T_{w}-T_{\infty })+h_{c,o}(\omega _{o,w}-\omega _{o,\infty})]

(30)

The energy balance at the surface of the solid fuel is

-{q}''_{w}={\dot{m}}''_{f}h_{sv}+{q}''_{\ell }

(31)

where the two terms on the right-hand side of eq. (31) represent the latent heat of sublimation, and the sensible heat required to raise the surface temperature of the solid fuel to sublimation temperature and heat loss to the solid fuel. Combining eqs. (30) and (31) yields the rate of sublimation on the solid fuel surface

{\dot{m}}''_{f}=Z\frac{\tau _{w}}{U_{\infty }}

(32)

where Z is the transfer driving force or transfer number defined as

Z=\frac{c_{p}(T_{\infty }-T_{w})+h_{c,o}(\omega _{o,\infty }-\omega _{o,w})}{h_{sv}+{q}''_{\ell }/\dot{{m}''}_{f}}

(33)

Using the friction coefficient gives

c_{f}=\frac{\tau _{w}}{\rho U_{\infty }^{2}/2}

(34)

eq. (32) therefore becomes

{\dot{m}}''_{f}=\frac{c_{f}}{2}\rho U_{\infty }Z

(35)

The surface blowing velocity of the gaseous fuel is then

v_{w}=\frac{{\dot{m}}''_{f}}{\rho }=\frac{c_{f}}{2}U_{\infty }Z

(36)

where the friction coefficient, cf, can be obtained from the solution of the boundary layer flow over a flat plate with blowing on the surface. The similarity solution of the boundary layer flow problem exists only if the blowing velocity satisfies v_{w}\propto x^{-1/2}. In this case, one can define a blowing parameter as

B=\frac{(\rho v)_{w}}{(\rho u)_{\infty }}\operatorname{Re}_{x}^{1/2}

(37)

Combining eqs. (36) and (37) yields

B=\frac{Z}{2}\operatorname{Re}_{x}^{1/2}c_{f}

(38)

Glassman [3] recommended an empirical form of eq. (38) based on numerical and experimental results:

B=\frac{\ln (1+Z)}{2.6Z^{0.15}}

(39)

=References

  1. Faghri, A., Zhang, Y., and Howell, J. R., 2010, Advanced Heat and Mass Transfer, Global Digital Press, Columbia, MO.
  2. 2.0 2.1 Kaviany, M., 2001, Principles of Convective Heat Transfer, 2nd ed., Springer-Verlag, New York.
  3. Glassman, I., 1987, Combustion, 2nd ed. Academic Press, Orlando, FL