Microfluid Mechanics/Models for Gases and Liquids

>back to Chapters of Microfluid Mechanics

>back to Chapters of Fluid Mechanics for Mechanical Engineers

Flowing medium

A fluid is composed of atoms and molecules. Depending on the phase of the fluid (gas,liquid or supercritical), the distance between the molecules shows orders of magnitude difference, being the largest in the gas phase and shortest in the liquid phase. As the distance between the molecules or the mean free path of the flowing medium approaches to the characteristic size of the flow device, the flow can not be treated as continuum.

In a solid, molecules form a regular lattice and oscillate around an equilibrium point. At this state, there is a strong attraction between the molecules and the kinetic energy of the molecules can not overcome this force in this phase of the matter. When enough energy is given to the molecules, e.g. by heating it, the matter melts and consequently becomes a liquid. The molecules gain kinetic energy as a result of added heat and start to move around in an irregular pattern. However, the density of liquids and solids, in other words the mean molecular distances at these two phases do not differ much from each other. When the liquid vaporizes and turns into the gas phase, the density drastically drops as the molecules starts to move freely between the intermolecular collisions.

Solid	Liquid	Gas
High density	High density	Low density
Low intermolecular distance (typical˜0.3 nm)	Low intermolecular distance (typical˜0.3 nm)	High intermolecular distance (typical˜3 nm)
Low kinetic kinetic energy of molecules	Higher kinetic energy of molecules	Highest kinetic kinetic energy of molecules
Molecules oscillates in a regular lattice arrangement.	Molecules build lattice forms over only short distances, but they move in an irregular pattern over longer distance.	Molecules moves freely between collisions.
Incompressible	Hardly compressible	Compressible

File:Phases of matter.svg

Cartoon showing the molecular difference between solids, liquids and gases.

Gas as a flowing medium

Macroscopic models for gas flow is better suited when the flow can be accepted as continuum (see continuum mechanics). Gas flows can not be treated as continuum, when their density drops or when they flow through micro or nano channels which have characteristic dimensions at the order of the mean free path of the gas molecules. Therefore, It might be necessary to model gas flows in microdevices in microscopic (molecular) level. Kinetic theory of gases Kinetic theory is a molecular model for gas flows. It recognizes the gas as a collection of many discrete molecules and, ideally, provides information on the position, velocity and the state of the molecules. This is done by following each molecule by using Newton's laws of motion and allowing binary collision of molecules with the other molecules and and the bounding walls. Macroscopic properties of gases, like density, pressure and temperature, can be best understood by utilizing this approach. Let us take a cubic container with a volume of $l^3$, which is full of $N$ simple gas (monatomic) molecules, and derive some of the macroscopic properties. We let the molecules to collide with each other and with the walls of the container. This concept is valid under the condition that collisions can be only of binary type and the molecules do not induce electromagnetic forces on each other. The mathematical expression for this conditions is as follows:

The temperature of an ideal monatomic gas is a measure of the average kinetic energy of its atoms. The size (Bohr radius) of helium atoms relative to their spacing is shown to scale under 1950 atmospheres of pressure. The atoms have a certain, average speed, slowed down here two trillion fold from room temperature.

The number density of the gas molecules in the container is: Template:Center top $n=\frac{N}{l^3}$ Template:Center bottom Hence, the average volume for one molecule in this container is $1/n$, so the mean molecular spacing $\delta$ reads Template:Center top $\delta =n^{-1/3}$ Template:Center bottom Gas molecules can be modeled by hard spheres, which has the effective molecular diameter $(d)$. Accordingly, the fraction of volume occupied by a molecule is in the order of $(d/\delta {)}^3$. For very low number densities, mean molecular spacing $(\delta )$ becomes relatively large compared to the effective diameter of the molecule $(d)$. Hence, a molecule can freely move outside region of influence of other molecules and almost all collisions of molecules can be accepted to be a binary one or collision with a wall. This situation can be written as: Template:Center top $\delta >>d$ Template:Center bottom and defines the circumstances in a dilute gas.

Density

The mass of each molecule can be calculated from the molar mass of the gas $M[kgmo{l}^{-1}]$ and the Avogadro constant $N_A=6.02214179(30)\times 10^{23}[mo{l}^{-1}]$ Template:Center top $m=\frac{M}{N_A}$ Template:Center bottom Hence, the density of the gas reads: Template:Center top $\rho =\frac{mN}{l^3}$ Template:Center bottom

Pressure

Pressure is explained by kinetic theory as arising from the force exerted by fluid molecules impacting on the walls of the container. When a gas molecule collides with the wall of the container perpendicular to the x coordinate axis and bounces off in the opposite direction with the same speed (an elastic collision), then the change of the momentum of the molecule ($\mathrm{\Delta }p$) upon wall collision is:

$\mathrm{\Delta }p=p_{final,x}-p_{initial,x}=-mu-mu=-2mu$

where $u$ is the x-component of the initial velocity of the molecule.

The molecule impacts one specific side wall once every

$\mathrm{\Delta }t=\frac{2l}{u}$

Hence the force exerted on the molecule by the wall is the rate of change of momentum per unit time:

$F^{*}=\frac{\mathrm{\Delta }p}{\mathrm{\Delta }t}=-\frac{m{u}^2}{l}.$

The force exerted by the molecule on the wall is in the negative direction:

$F=\frac{m{u}^2}{l}.$

Not all the molecules moves with the same speed. Say the speed of the molecules varies over the range $u_{max}$ and $u_{min}$, so the number of molecules moving over the range $u+{\mathrm{\Delta }}_u$ is denoted by $N_u$. Using this distribution, it is possible to calculate the total force exerted by the molecules on the wall in the x-direction is the sum over all molecules:

$\begin{array}{ccc}{F}_{total,x}& =& \mathrm{\Sigma }{N}_u\frac{m{u}^2}{l}\\ & & \\ & =& \frac{m}{l}N\frac{\mathrm{\Sigma }{N}_u{u}^2}{N}\\ & & \\ & =& \frac{m}{l}N\overline{u^2}\end{array}$

where $\overline{u^2}$ is the average of the square of x-component velocity of the molecules. As the pressure is the force per unit area, the x-component of the pressure reads

$P_x=\frac{m}{l}N\overline{u^2}\frac{1}{l^2}=\frac{Nm}{l^3}\overline{u^2}$

Other components of the pressures can be similarly derived:

$P_y=\frac{Nm}{l^3}\overline{v^2}$

$P_z=\frac{Nm}{l^3}\overline{w^2}$

When the cubic container is small enough, the pressure on the walls of the container can be approximated by the average of them

$\begin{array}{ccc}P& =& \frac{1}{3}(P_x+P_y+P_z)\\ & & \\ & =& \frac{Nm}{3l^3}(\overline{u^2}+\overline{v^2}+\overline{w^2})\\ & & \\ & =& \frac{Nm}{3l^3}\overline{c^2}\end{array}$

where $\overline{c^2}=\overline{u^2}+\overline{v^2}+\overline{w^2}$ is the mean square velocity of the molecules of the simple gas.

Rearranging the terms shows that Template:Center top

$P=\frac{2N}{3l^3}\frac{1}{2}m\overline{c^2}$ This is a first non-trivial result of the kinetic theory because it relates pressure, a macroscopic property, to the average (translational) kinetic energy per molecule $\frac{1}{2}m\overline{c^2}$ which is a microscopic property.

Since $Nm$ is the mass and $l^3$ is the volume, the pressure can also be written

Template:Center top $\begin{array}{ccc}P& =& \frac{1}{3}\rho \overline{c^2}\\ & & \\ & =& \frac{2}{3}E\end{array}$ Template:Center bottom where $E=\frac{1}{2}\rho \overline{c^2}$ is the kinetic energy of the molecules per unit volume.

Temperature and kinetic energy

After defining the pressure by the kinetic approach, it is possible to define temperature also by utilizing the ideal gas law, which is:

$P=\rho R_{\mathrm{s}\mathrm{p}\mathrm{e}\mathrm{c}\mathrm{i}\mathrm{f}\mathrm{i}\mathrm{c}}T\text{or}P=n{k}_{B}T$

where the specific gas constant is defined as

$R_{\mathrm{s}\mathrm{p}\mathrm{e}\mathrm{c}\mathrm{i}\mathrm{f}\mathrm{i}\mathrm{c}}=\frac{R}{M}\text{and}R=N_A{k}_B\to R_{\mathrm{s}\mathrm{p}\mathrm{e}\mathrm{c}\mathrm{i}\mathrm{f}\mathrm{i}\mathrm{c}}=\frac{k_B}{m}$

in which $R=8.314472(15)[Jmo{l}^{-1}{K}^{-1}]$, $M[kgmo{l}^{-1}]$, $m[kg]$ and $k_B=1.3806504(24)\times 10^{-23}[J{K}^{-1}]$ are the gas constant, the molar mass of the gas molecules, mass of one molecule and Boltzmann constant, respectively. It should be noted that ideal gas law is valid under ideal conditions, viz. at low pressure and at high temperatures, where $\delta >>d$, so that the molecular interaction is negligible.

Equating the pressure with the pressure derived by using kinetic approach

$\frac{1}{3}\rho \overline{c^2}=\rho R_{\mathrm{s}\mathrm{p}\mathrm{e}\mathrm{c}\mathrm{i}\mathrm{f}\mathrm{i}\mathrm{c}}T,$

gives the temperature as:

$T=\frac{1}{3}\frac{\overline{c^2}}{R_{\mathrm{s}\mathrm{p}\mathrm{e}\mathrm{c}\mathrm{i}\mathrm{f}\mathrm{i}\mathrm{c}}}$ This equation reveals that at absolute zero temperature, the molecules will stop moving. Template:Center bottom

However, this result is not justified as gas ceases to behave as a perfect gas before reaching to the absolute zero of temperature.

Other forms can be obtained by replacing the specific gas constant with $R_{\mathrm{s}\mathrm{p}\mathrm{e}\mathrm{c}\mathrm{i}\mathrm{f}\mathrm{i}\mathrm{c}}=\frac{k_B}{m}$ as follows: Template:Center top

$T=\frac{1}{3}\frac{1}{k_B}m\overline{c^2}$

Template:Center bottom where $m$ is the mass of one gas molecule.

Template:Center top After grouping the terms properly, it can be seen that the kinetic energy of the molecules directly proportional to the temperature: $\frac{1}{2}m\overline{c^2}=\frac{3}{2}{k}_{B}T$ In other words, by heating a thermodynamic system, the kinetic energy of the molecules is increased. Template:Center bottom

Mean free path

The distance traveled by a molecule between two successive collision is called the mean free path. Under the assumptions that the molecules are hard spheres and the probability distribution of the molecular velocities follows that of Maxwell distribution law, the mean free path can be derived as: Template:Center top

$\lambda =\frac{1}{\sqrt{2}\pi d^{2}n}$

Template:Center bottom where $d$ and $n$ are the effective diameter of the gas molecule and the number density of the gas in the container.

Template:Center top While for an ideal gas $n=\frac{P}{k_{B}T}$, at a constant temperature $\lambda \propto \frac{1}{P}$ It means that mean free path increases with decreasing pressure. However, as the mean free path approaches to the size of the container, molecules starts to have collisions solely with the walls of the container. Template:Center bottom

Template:Center top

Template:Center top Hence, a key nondimensional parameter for microflows is the Knudsen number: Template:Center top $Kn=\frac{\lambda }{L}$ which characterizes the relative size of the mean free path to the characteristic size of the container $(L)$. Template:Center bottom Template:Center bottom

Molecular speeds and the speed of sound

Although the macroscopic velocity under adiabatic conditions in the container is zero, the molecules are still moving around with the mean square velocity of $\overline{c^2}$. The root mean square (rms) velocity of the molecules $c^{\prime }$ gives an idea of the speed of the molecules: Template:Center top $c^{\prime }=\sqrt{\overline{c^2}}=\sqrt{\frac{3P}{\rho }}=\sqrt{3R_{\mathrm{s}\mathrm{p}\mathrm{e}\mathrm{c}\mathrm{i}\mathrm{f}\mathrm{i}\mathrm{c}}T}$ Template:Center bottom

Two more statistically interesting speeds can be extracted from the Maxwell distribution of velocity:

$\begin{array}{cccc}\bar{c}& =& \sqrt{\frac{8}{\pi }{R}_{\mathrm{s}\mathrm{p}\mathrm{e}\mathrm{c}\mathrm{i}\mathrm{f}\mathrm{i}\mathrm{c}}T}& \text{:average molecular speed}\\ c_p& =& \sqrt{2R_{\mathrm{s}\mathrm{p}\mathrm{e}\mathrm{c}\mathrm{i}\mathrm{f}\mathrm{i}\mathrm{c}}T}& \text{:most probable molecular speed}\end{array}$

The speeds are ordered as follows:

$0.886\bar{c}=c_p<\bar{c}<c^{\prime }=1.085\bar{c}$

The speed of sound in an ideal gas is

$c_s=\sqrt{\gamma \frac{p}{\rho }}$

where $\gamma =\frac{C_p}{C_v}$ is the heat capacity ration and it takes the value 5/3 = 1.6667 for monoatomic and 7/5 = 1.4 for diatomic gas molecules. For ideal gas and air $\gamma$ is 1.4. Thus, in an ideal gas, the ratio between the speed of sound and rms velocity of the molecules is:

Template:Center top $\frac{c_s}{c^{\prime }}=\sqrt{\frac{\gamma }{3}}\approx 0.68$ Template:Center bottom which reveals that the speed of sound is slower than the all molecular speeds but it is at the same order of magnitude. This fact is not surprising while the waves are carried by molecular motion.

Viscosity

The viscosity of the fluid is responsible for the momentum transport in molecular level. The dynamic viscosity of an ideal gas, which is derived in the frame of kinetic theory, is: Template:Center top $\mu =\frac{1}{3}\rho \lambda \bar{c}$ Template:Center bottom substituting the $\lambda$ and $\rho$ in the above equation results in Template:Center top $\mu =\frac{1}{3}\frac{m}{\sqrt{2}\pi d^2}\bar{c}=\frac{1}{3}\frac{m}{\sqrt{2}\pi d^2}\sqrt{\frac{8}{\pi }{R}_{\mathrm{s}\mathrm{p}\mathrm{e}\mathrm{c}\mathrm{i}\mathrm{f}\mathrm{i}\mathrm{c}}T}$ Template:Center bottom which states that the dynamic viscosity is only a function of the mass, the diameter and the average speed of the molecules, which is in turn a function of measured temperature. In other words, it is supposedly independent of pressure and density, which is in contradiction to experimental findings showing that viscosity is pressure dependent. In addition to that , this relation does not hold at very low pressures at which the mean free path reaches to the size of the container or flow apparatus and at very high pressures at which intermolecular forces can not be neglected. Furthermore, viscosity was observed to increase much faster than the $\sqrt{T}$ proportionality suggests. Nevertheless, this relation holds from a few mbars to several bars of pressure. A more accurate equation had been derived by using the intermolecular forces .

Liquid as a flowing medium

The molecules are too close to each other in a liquid, so that strong intermolecular interactions dominates the motion of the molecules. Moreover, the mean free path of the liquid molecules $(\lambda )$ are in the order of the mean molecular spacing $(\delta )$.

As the distance between the molecules becomes shorter, the forces between the molecules should be accounted properly. The intermolecular interaction between two molecules separated by a distance $r$ can be modeled by the Lennard-Jones pair pottential. Template:Center top $w(r)=\underset{negative}{\underbrace{-\frac{A}{r^6}}}+\underset{positive}{\underbrace{\frac{B}{r^{12}}}}[J]$ Template:Center bottom where $A$ and $B$ are molecule dependent constants and they are equal to Template:Center top $A=4\epsilon {\sigma }^6\text{and}B=4\epsilon {\sigma }^{12}$ Template:Center bottom in which $\epsilon$ is the depth of the potential well (maximum attraction energy) and $\sigma$ is the finite distance at which the inter-molecule potential is zero. $\sigma$ is also called as the collision diameter. Inserting $A$ and $B$ in the above interaction potential reveals the other form of the Lennard-Jones pair pottential: Template:Center top $w(r)=4\epsilon [-{\left(\frac{\sigma }{r}\right)}^6+{\left(\frac{\sigma }{r}\right)}^{12}]$ The van der Waals potential is the negative contribution, which varies with the inverse 6th power of the distance $r$. The intermolecular force is the derivative of the pair potential w.r.t. the distance. Template:Center top $F(r)=-\frac{dw(r)}{dr}=\underset{attractive}{\underbrace{-6\frac{A}{r^7}}}+\underset{repulsive}{\underbrace{12\frac{B}{r^{13}}}}[N]$ Template:Center bottom Hence, it can be seen that the first term is the attractive force, i.e. the van der Waals force. The second term is the repulsive force and it hinders that two atoms and / or molecules do not exist at the same location at an instant. Note that, force and potential energy are directly related. The potential energy (Lennard-Jones potential) is the negative of the work done by the net electromagnetic force moving the molecule towards another molecule to its given position in space from infinity.

File:Lennard jones potential force.png Lennard-Jones potential and the intermolecular force derived from it. Selected constants are A=1e-77 [J m6] & B=1e-134 [J m12].

Definition of pressure

The pressure on a wall in liquids is a result of molecular collisions on a wall and the gravitation. In microfluidic systems, the effect of gravity is insignificant. However, the molecular collisions can not be easily modeled as was done for dilute gases by kinetic theory, because liquid molecules are densely packed and their motion highly influenced by the intermolecular interactions.

Viscosity

The derivation of viscosity using kinetic approach is not straightforward as the electromagnetic molecular interactions play dominant role in the liquids. One rough approximation was derived by Eyring and coworkers , which is: Template:Center top

$\mu \approx \frac{N_{A}h\rho }{M}\exp (3.8T_B/T)$

Template:Center bottom where $h$ and $T_B$ are the planck constant and the boiling temperature of the liquid under one atmosphere, respectively. This equation, though being only an approximation, gives the observed decrease of viscosity with increasing temperature.

External Links

• This applet shows basic features of gas flow as the kinetic theory suggests.