Fully Developed Laminar Flow and Temperature Profile

From Thermal-FluidsPedia

Revision as of 02:51, 7 July 2010 by Founding Editor (Talk | contribs)
(diff) ← Older revision | Current revision (diff) | Newer revision → (diff)
Jump to: navigation, search

Fully Developed Velocity Distribution

The velocity distribution is derived from the continuity and momentum equations along with the slip condition. The velocity profile is expressed in terms of Knudsen number. The model for the analysis of velocity distribution can be considered as the flow of a fluid in a circular tube of radius ro, as shown in Fig. 5.22. The continuity equation in cylindrical coordinates is given as:

\frac{\partial \rho }{\partial t}+\frac{1}{r}\frac{\partial }{\partial r}\left( \rho rv \right)+\frac{\partial }{\partial x}\left( \rho u \right)=0

(1)

For steady and fully developed flow of an incompressible fluid, it becomes:

\frac{\partial u}{\partial x}=0

(1)

Therefore, rv = constant. Since the radial velocity at the wall is zero (impermeability condition), it can be concluded that v = 0 everywhere in the flow field. The steady parabolic momentum equation in the x-direction can be written in cylindrical coordinates as:

\rho u\frac{\partial u}{\partial x}+\rho v\frac{\partial u}{\partial r}=\frac{\mu }{r}\frac{\partial }{\partial r}\left( r\frac{\partial u}{\partial r} \right)-\frac{\partial p}{\partial x}

(1)

For fully developed, steady flow of incompressible fluid, it reduces to:

-\frac{dp}{dx}+\frac{\mu }{r}\frac{d}{dr}\left( r\frac{du}{dr} \right)=0

(1)

As the velocity is independent of x, a constant C1 can be defined as:

{{C}_{1}}=-\frac{1}{4\mu }\frac{dp}{dx}

(1)

Therefore, the momentum equation reduces to:

4{{C}_{1}}+\frac{1}{r}\frac{d}{dr}\left( r\frac{du}{dr} \right)=0

(1)

The above equation is integrated twice to yield:

u = C2lnr + C3C1r2

(1)

Since the velocity, u, is finite at the center of the tube (r = 0), we conclude C2=0. At the centerline of the tube (r = 0), the velocity is equal to uc, therefore, u(0) = uc = C3. At the inner surface of the tube wall (r = ro), the velocity is not zero due to the slip flow condition, but is equal to a finite velocity us:

u({{r}_{o}})={{u}_{s}}={{u}_{c}}-{{C}_{1}}r_{o}^{2}

(1)

Rearranging and solving for C1 we have:

{{C}_{1}}=\frac{{{u}_{c}}-{{u}_{s}}}{r_{0}^{2}}=-\frac{1}{4\mu }\frac{dp}{dx}

(1)

The velocity profile is obtained as follows:

u={{u}_{c}}-({{u}_{c}}-{{u}_{s}}){{(r/{{r}_{o}})}^{2}}={{u}_{c}}[1-{{(r/{{r}_{o}})}^{2}}]+{{u}_{s}}{{(r/{{r}_{o}})}^{2}}

(1)

The mean velocity in the tube, um, is now evaluated. The volumetric flow rate is:

\pi r_{0}^{2}{{u}_{m}}=\int_{0}^{{{r}_{0}}}{2\pi rudr}

(1)

therefore,

um = 2[uc / 2 − (ucus) / 4] = (uc + us) / 2

(1)

The centerline velocity can be written as:

uc = 2umus

(1)

Using the above relation, we get the following velocity profile in a microchannel:

u = 2(umus)(1 − r + 2) + us

(1)

where r+=r/ro. If the slip velocity is zero, the velocity distribution reduces to the Poiseuille flow distribution for conventional channels. The velocity slip is given by the following condition [see eq. (5.194)]:

{{u}_{s}}=-\frac{2-F}{F}\lambda {{\left. \frac{du}{dr} \right|}_{r={{r}_{o}}}}

(1)

For most applications, F has values near unity. Therefore,

{{u}_{s}}=-\lambda {{\left. \frac{du}{dr} \right|}_{r={{r}_{o}}}}

(1)


The slip velocity is derived using eq. (5.209) for velocity profile:

{{u}_{s}}=-\frac{\lambda }{{{r}_{o}}}{{\left. \frac{du}{d{{r}^{+}}} \right|}_{{{r}^{+}}=1}}=\frac{4\lambda ({{u}_{m}}-{{u}_{s}})}{{{r}_{o}}}=\frac{8\lambda ({{u}_{m}}-{{u}_{s}})}{D}

(1)

The Knudsen number, Kn = λ / D, is introduced in the above relation to obtain the following equation for the slip velocity:

\frac{{{u}_{s}}}{{{u}_{m}}}=\frac{8\text{Kn}}{1+8\text{Kn}}

(1)

Substituting the above into eq. (5.209), the non-dimensional expression for the velocity profile is obtained:

\frac{u}{{{u}_{m}}}=\frac{2(1-{{r}^{+2}})+8\text{Kn}}{1+8\text{Kn}}

(1)

The velocity distribution is shown in Fig. 5.23 for different Kn numbers. The slip velocity at the wall increases with an increasing Kn number.

Fully Developed Heat Transfer Coefficient in Microchannel The flow is considered to be steady, laminar, and fully developed both hydrodynamically and thermally. The conventional continuum approach is coupled with the two main characteristics of the microscale phenomena, the velocity slip and the temperature jump, as noted above. The energy equation including the viscous dissipation term for steady, fully developed flow in a pipe and neglecting axial conduction is:

u\frac{\partial T}{\partial x}=\frac{\alpha }{r}\frac{\partial }{\partial r}\left( r\frac{\partial T}{\partial r} \right)+\frac{\nu }{{{c}_{p}}}{{\left( \frac{\partial u}{\partial r} \right)}^{2}}

(1)

The variation of temperature in the r-direction at the center is zero due to axisymmetric conditions:

{{\left. \frac{\partial T}{\partial r} \right|}_{r=0}}=0

(1)

Also, the temperature jump condition for the fluid at the wall is written as follows:

{{T}_{s}}-{{T}_{w}}=-\frac{2-{{F}_{t}}}{{{F}_{t}}}\frac{4\gamma }{1+\gamma }\frac{\text{Kn}}{\Pr }{{\left. \frac{\partial T}{\partial r} \right|}_{r={{r}_{o}}}}

(1)

Both conventional boundary conditions of constant heat flux and constant temperature at the wall are solved analytically by Aydin and Avci (2006) for microchannel with circular cross-section, which are presented below.