# Continuum flow limitations

### From Thermal-FluidsPedia

Line 1: | Line 1: | ||

- | The transport phenomena are usually modeled in continuum states for most applications – the materials are assumed to be continuous and the fact that matter is made of atoms is ignored. Recent development in fabrication and utilization of nanotechnology, micro devices, and microelectromechanical systems (MEMS) requires noncontinuum modeling of transport phenomena in nano- and microchannels. When the characteristic dimension, | + | The transport phenomena are usually modeled in continuum states for most applications – the materials are assumed to be continuous and the fact that matter is made of atoms is ignored. Recent development in fabrication and utilization of nanotechnology, micro devices, and microelectromechanical systems (MEMS) requires noncontinuum modeling of transport phenomena in nano- and microchannels. When the characteristic dimension, <math>L</math>, is small compared to the molecular mean free path, λ, which is defined as average distance between collisions for a molecule, the traditional Navier-Stokes equation and the energy equation based on the continuum assumption have failed to provide accurate results. The continuum assumption also fails when the gas is at very low pressure (rarefied). |

The continuum assumption may also not be valid in conventional sized systems – for example, the early stages of high-temperature heat pipe startup from a frozen state (Cao and Faghri, 1993) and microscale heat pipes (Cao and Faghri, 1994). During the early stage of startup of high-temperature heat pipes, the vapor density in the heat pipe core is very low and partly loses its continuum characteristics. The vapor flow in this condition is usually referred to as rarefied vapor flow. Because of the low density, the vapor in the rarefied state is somewhat different from the conventional continuum state. Also, the vapor density gradient is very large along the axial direction of the heat pipe. The vapor flow along the axial direction is caused mainly by the density gradient via vapor molecular diffusion. | The continuum assumption may also not be valid in conventional sized systems – for example, the early stages of high-temperature heat pipe startup from a frozen state (Cao and Faghri, 1993) and microscale heat pipes (Cao and Faghri, 1994). During the early stage of startup of high-temperature heat pipes, the vapor density in the heat pipe core is very low and partly loses its continuum characteristics. The vapor flow in this condition is usually referred to as rarefied vapor flow. Because of the low density, the vapor in the rarefied state is somewhat different from the conventional continuum state. Also, the vapor density gradient is very large along the axial direction of the heat pipe. The vapor flow along the axial direction is caused mainly by the density gradient via vapor molecular diffusion. | ||

Line 7: | Line 7: | ||

The continuum criterion is usually expressed in terms of the Knudsen number | The continuum criterion is usually expressed in terms of the Knudsen number | ||

- | <center><math>{\rm{Kn}} = \frac{\lambda }{L}</math> | + | <center><math>{\rm{Kn}} = \frac{\lambda }{L} \qquad\qquad(1)</math></center> |

Based on the degree of rarefaction of gas or the gases size, the flow regimes in various devices can be classified into four regimes: | Based on the degree of rarefaction of gas or the gases size, the flow regimes in various devices can be classified into four regimes: | ||

- | 1. ''Continuum regime'' (<math>{\rm{Kn}} < 0.001</math>). The Navier-Stokes and energy equations are valid (with no-slip/no jump boundary conditions). | + | 1.''Continuum regime'' (<math>{\rm{Kn}} < 0.001</math>). The Navier-Stokes and energy equations are valid (with no-slip/no jump boundary conditions). |

- | 2. ''Slip flow regime'' (<math>{\rm{0}}{\rm{.001 < Kn}} < 0.1</math>). The Navier-Stokes and energy equations can be used with the application of slip or jump boundary conditions, i.e., allowing non-zero axial fluid velocity near the wall of the object. | + | |

- | 3. ''Transition regime'' (<math>{\rm{0}}{\rm{.1 < Kn}} < 10</math>). The Navier-Stokes equation is not valid, and the flow must be solved using molecular based models such as the Boltzmann equation or Direct Simulation Monte Carlo (DSMC). | + | 2.''Slip flow regime'' (<math>{\rm{0}}{\rm{.001 < Kn}} < 0.1</math>). The Navier-Stokes and energy equations can be used with the application of slip or jump boundary conditions, i.e., allowing non-zero axial fluid velocity near the wall of the object. |

- | 4. ''Free molecular flow regime'' (<math>{\rm{Kn > 10}}</math>). The collision between molecules can be neglected and a collisionless Boltzmann equation can be used. | + | |

+ | 3.''Transition regime'' (<math>{\rm{0}}{\rm{.1 < Kn}} < 10</math>). The Navier-Stokes equation is not valid, and the flow must be solved using molecular based models such as the Boltzmann equation or Direct Simulation Monte Carlo (DSMC). | ||

+ | |||

+ | 4.''Free molecular flow regime'' (<math>{\rm{Kn > 10}}</math>). The collision between molecules can be neglected and a collisionless Boltzmann equation can be used. | ||

In the slip flow regime, the slip boundary condition refers to circumstance when the tangential velocity of the fluid at the wall is not the same as the wall velocity. Temperature jump is similarly defined as when the temperature of the fluid next to the wall is not the same as the wall temperature. | In the slip flow regime, the slip boundary condition refers to circumstance when the tangential velocity of the fluid at the wall is not the same as the wall velocity. Temperature jump is similarly defined as when the temperature of the fluid next to the wall is not the same as the wall temperature. | ||

- | The mean free path | + | The mean free path for dilute gases based on the kinetic theory can be rewritten in terms of temperature and pressure |

- | <center><math>\lambda = \frac{{1.051{k_b}T}}{{\sqrt 2 \pi {\sigma ^2}p}}</math> | + | <center><math>\lambda = \frac{{1.051{k_b}T}}{{\sqrt 2 \pi {\sigma ^2}p}} \qquad\qquad(2)</math></center> |

The transition density under which the continuum assumption is invalid can be obtained by combining eqs. (1) and (2) and <math>{\rm{Kn = }}0.001</math> i.e., | The transition density under which the continuum assumption is invalid can be obtained by combining eqs. (1) and (2) and <math>{\rm{Kn = }}0.001</math> i.e., | ||

- | <center><math>{\rho _{tr}} = \frac{{1.051{k_b}}}{{\sqrt 2 \pi {\sigma ^2}{R_g}D{\rm{Kn}}}}</math> | + | <center><math>{\rho _{tr}} = \frac{{1.051{k_b}}}{{\sqrt 2 \pi {\sigma ^2}{R_g}D{\rm{Kn}}}} \qquad\qquad(3)</math></center> |

- | where the ideal gas equation of state, <math>p = \rho {R_g}T</math> was used. | + | where the ideal gas equation of state, <math>p = \rho {R_g}T</math> was used. Assuming that the vapor is in the saturation state, the transition vapor temperature <math>{T_{tr}}</math> corresponding to the transition density can be obtained by using the Clausius-Clapeyron equation (see Chapter 2) combined with the equation of state: |

- | + | ||

- | Assuming that the vapor is in the saturation state, the transition vapor temperature | + | |

- | <center><math>{T_{tr}} = \frac{{{p_{sat}}}}{{\rho {R_g}}}\exp \left[ { - \frac{{{h_{\ell v}}}}{{{R_g}}}\left( {\frac{1}{{{T_{tr}}}} - \frac{1}{{{T_{sat}}}}} \right)} \right]</math> | + | <center><math>{T_{tr}} = \frac{{{p_{sat}}}}{{\rho {R_g}}}\exp \left[ { - \frac{{{h_{\ell v}}}}{{{R_g}}}\left( {\frac{1}{{{T_{tr}}}} - \frac{1}{{{T_{sat}}}}} \right)} \right] \qquad\qquad(4)</math></center> |

- | where | + | where <math>{P_{sat}}</math> and <math>{T_{sat}}</math> are the saturation pressure and temperature, <math>{h_{\ell v}}</math> is the latent heat of vaporization, and the vapor density <math>P</math> is given by eq. (3). Equation (4) can be rewritten as |

- | <center><math>\ln \left( {\frac{{{T_{tr}}\rho {R_g}}}{{{p_{sat}}}}} \right) + \frac{{{h_{\ell v}}}}{{{R_g}}}\left( {\frac{1}{{{T_{tr}}}} - \frac{1}{{{T_{sat}}}}} \right) = 0</math> | + | <center><math>\ln \left( {\frac{{{T_{tr}}\rho {R_g}}}{{{p_{sat}}}}} \right) + \frac{{{h_{\ell v}}}}{{{R_g}}}\left( {\frac{1}{{{T_{tr}}}} - \frac{1}{{{T_{sat}}}}} \right) = 0 \qquad\qquad(5)</math></center> |

and solved iteratively for <math>{T_{tr}}</math> using the Newton-Raphson/secant method. The transition vapor temperature is the boundary between the continuum and noncontinuum regimes. | and solved iteratively for <math>{T_{tr}}</math> using the Newton-Raphson/secant method. The transition vapor temperature is the boundary between the continuum and noncontinuum regimes. |

## Revision as of 22:25, 5 October 2009

The transport phenomena are usually modeled in continuum states for most applications – the materials are assumed to be continuous and the fact that matter is made of atoms is ignored. Recent development in fabrication and utilization of nanotechnology, micro devices, and microelectromechanical systems (MEMS) requires noncontinuum modeling of transport phenomena in nano- and microchannels. When the characteristic dimension, *L*, is small compared to the molecular mean free path, λ, which is defined as average distance between collisions for a molecule, the traditional Navier-Stokes equation and the energy equation based on the continuum assumption have failed to provide accurate results. The continuum assumption also fails when the gas is at very low pressure (rarefied).

The continuum assumption may also not be valid in conventional sized systems – for example, the early stages of high-temperature heat pipe startup from a frozen state (Cao and Faghri, 1993) and microscale heat pipes (Cao and Faghri, 1994). During the early stage of startup of high-temperature heat pipes, the vapor density in the heat pipe core is very low and partly loses its continuum characteristics. The vapor flow in this condition is usually referred to as rarefied vapor flow. Because of the low density, the vapor in the rarefied state is somewhat different from the conventional continuum state. Also, the vapor density gradient is very large along the axial direction of the heat pipe. The vapor flow along the axial direction is caused mainly by the density gradient via vapor molecular diffusion.

The validity of the continuum assumption can also be violated in micro heat pipes. As the size of the heat pipe decreases, the vapor in the heat pipe may lose its continuum characteristics. The heat transport capability of a heat pipe operating under noncontinuum vapor flow conditions is very limited, and a large temperature gradient exists along the heat pipe length. This is especially true for miniature or micro heat pipes, whose dimensions may be extremely small.

The continuum criterion is usually expressed in terms of the Knudsen number

Based on the degree of rarefaction of gas or the gases size, the flow regimes in various devices can be classified into four regimes:

1.*Continuum regime* (Kn < 0.001). The Navier-Stokes and energy equations are valid (with no-slip/no jump boundary conditions).

2.*Slip flow regime* (0.001 < Kn < 0.1). The Navier-Stokes and energy equations can be used with the application of slip or jump boundary conditions, i.e., allowing non-zero axial fluid velocity near the wall of the object.

3.*Transition regime* (0.1 < Kn < 10). The Navier-Stokes equation is not valid, and the flow must be solved using molecular based models such as the Boltzmann equation or Direct Simulation Monte Carlo (DSMC).

4.*Free molecular flow regime* (Kn > 10). The collision between molecules can be neglected and a collisionless Boltzmann equation can be used.

In the slip flow regime, the slip boundary condition refers to circumstance when the tangential velocity of the fluid at the wall is not the same as the wall velocity. Temperature jump is similarly defined as when the temperature of the fluid next to the wall is not the same as the wall temperature.

The mean free path for dilute gases based on the kinetic theory can be rewritten in terms of temperature and pressure

The transition density under which the continuum assumption is invalid can be obtained by combining eqs. (1) and (2) and Kn = 0.001 i.e.,

where the ideal gas equation of state, *p* = ρ*R*_{g}*T* was used. Assuming that the vapor is in the saturation state, the transition vapor temperature *T*_{tr} corresponding to the transition density can be obtained by using the Clausius-Clapeyron equation (see Chapter 2) combined with the equation of state:

where *P*_{sat} and *T*_{sat} are the saturation pressure and temperature, is the latent heat of vaporization, and the vapor density *P* is given by eq. (3). Equation (4) can be rewritten as

and solved iteratively for *T*_{tr} using the Newton-Raphson/secant method. The transition vapor temperature is the boundary between the continuum and noncontinuum regimes.

## References

Cao, Y., and Faghri, A., 1993, “Simulation of the Early Startup Period of High Temperature Heat Pipes From the Frozen State by a Rarefied Vapor Self-Diffusion Model,” ASME *Journal of Heat Transfer*, Vol. 115, pp. 239-246.

Cao, Y., and Faghri, A., 1994, “Micro/Miniature Heat Pipes and Operating Limitations,” *Journal of Enhanced Heat Transfer*, Vol. 1, pp. 265-274.