COMSOL 6.4 - Slip-Flow Benchmark

As the absolute pressure in a gas flow is reduced, the mean free path of the gas molecules begins to approach the size of the vessel through which the flow occurs. Such rarefied flows are characterized by a parameter known as the Knudsen number, which is the ratio of the mean free path to the characteristic length of the geometry. In a rarefied flow the gas cannot be treated as a continuum and the kinetic nature of the flow must be considered. However, at moderate Knudsen numbers between approximately 0.01 and 0.1, the flow can be modeled using the Navier–Stokes equations, except in a thin region adjacent to the walls, termed the Knudsen layer. This regime is termed the slip-flow regime. In the slip-flow regime the Knudsen layer can be replaced by alternative boundary conditions for the Navier–Stokes equations. A slip velocity along the geometry walls develops, and the temperature at the walls of the structure becomes discontinuous at this level of approximation. Phenomena such as viscous slip and thermal creep or transpiration become important.

This model is a benchmark model for the Slip Flow interface. It is based on analytic and numeric calculations presented in Ref. 1. Air at atmospheric pressure flows through a conducting microchannel connecting two reservoirs maintained at different temperatures. A flow between the two reservoirs develops as a result of thermal creep along the channel wall, which in turn produces a pressure gradient. At steady state the net flow through the channel is zero, but a pressure gradient exists between the hot and cold sides of the reservoir and a circulating flow occurs.

Model Definition

The model consists of two chambers 7.5 μm by 15 μm, with a depth significantly larger than either of these dimensions (this means that the simulation can be performed using the 2D XY interface). The chambers are linked by a narrow channel of the same depth with width 1.5 μm and length 15 μm. The chamber and channel walls are fabricated from silicon and have a thickness of 1 μm. The walls of the two chambers are in thermal contact with heat sinks maintained at 300 K and 400 K, respectively. The channel walls are thermally insulated. Figure 1 shows the model geometry. The gas in the center of the channel is at a pressure of 1 atmosphere.

Figure 1: The model geometry.

For channels of micron scale dimensions the Knudsen number becomes greater than 0.01 even at atmospheric pressure. It is therefore necessary to use a slip condition on the surfaces of walls in the vicinity of the channel.

Slip Wall Boundary Condition

The slip velocity, uslip, along the walls of the channel is given by (Ref. 2):

(1)

where λ is the mean free path of the gas, n is the boundary normal, τ is the viscous stress tensor, Tw is the wall temperature, Tg is the temperature of the gas, μ is its viscosity, and ρ is its density. The slip coefficients: σs (the viscous slip coefficient), σT (the thermal slip coefficient), and ζT (the temperature jump coefficient) can be defined by material properties and the tangential momentum accommodation coefficient, av, within a generalized form of Maxwell’s original slip model (see Ref. 2 and Ref. 3). For a model in which the surface reflects some molecules diffusely and some specularly, av is the fraction of molecules which are reflected diffusely. Within Maxwell’s model the slip coefficients are:

where κ is the thermal conductivity of the gas.

The mean free path can be computed from the gas properties using the following equation (Ref. 2):

Results and Discussion

Figure 2 shows the mean free path of the gas throughout the simulation domain. Since the channel width is 1.5 μm, a mean free path of 75 nm corresponds to a Knudsen number of 0.05. This value is comparable to the Knudsen number of the gas in the channel and is within the slip-flow regime.

The velocity of the gas in the x direction is shown in Figure 3, together with the velocity streamlines. In the steady state there is no net flow through the channel, but a flow parallel to the walls, in the direction of the thermal gradient (cold to hot), develops due to thermal creep. To compensate for this flow a back flow develops in the center of the channel, which is driven by a pressure gradient in the gas.

The temperature distribution in the channels is shown in Figure 4. The thermal gradient is predominately parallel to the walls, however there are some normal thermal gradients in the vicinity of the opening to the channel, which result in a temperature jump between the gas and the wall, according to Equation 1. This temperature jump simulates the effect of the Knudsen layer, which is not captured by the continuum Navier–Stokes equations.

Figure 2: Gas mean free path. The channel width is 1.5 μm, so the Knudsen number varies between 0.064 and 0.045.

Figure 3: x-velocity (color) and streamlines (white) inside the channel. Flow is driven by thermal creep along the walls of the channel but a back-flow, driven by the pressure difference between the reservoirs, occurs in the center of the channel.

Figure 4: Temperature contours within the model. Note that a temperature jump occurs between the vessel walls and the gas when normal heat fluxes occur into the wall from the gas.

Figure 5: Pressure on the wall of the channel, as a function of position.

Figure 5 shows the relative pressure acting on the channel wall, as a function of position along the wall. The pressure is lowest in the coolest chamber, at 151 Pa below atmospheric, and highest in the warmer chamber, at 179 Pa above atmospheric pressure. The asymmetry of the pressure distribution about the center of the channel results from changes in the number density of the gas in the two reservoirs. The total pressure difference of 330 Pa is similar to the value predicted in Ref. 1 for a numerical simulation of an equivalent problem (336 Pa), and is also in good agreement with the value predicted by the analytic model of Ref. 1, which is (335 Pa). This model therefore serves as a benchmark for the Slip Flow interface.

References

1. G. Kariadakis, A. Beskok, and N. Aluru, Microflows and Nanoflows, Springer Science and Business Media, 2005.

2. E.H. Kennard, Kinetic Theory of Gases, McGraw Hill, New York, 1938.

3. J.C. Maxwell, “On stresses in rarefied gases arising from inequalities of temperature,” Phil. Trans. R. Soc. Lond., vol. 170, pp. 231–256, 1879.

Microfluidics_Module/Rarefied_Flow/slip_flow_benchmark

Modeling Instructions

From the menu, choose .

New

In the window, click