Open Journal of Fluid Dynamics
Vol.09 No.01(2019), Article ID:90537,27 pages
10.4236/ojfd.2019.91002

Flows of a Rarefied Gas between Coaxial Circular Cylinders with Nonuniform Surface Properties

Toshiyuki Doi

Department of Applied Mathematics and Physics, Graduate School of Engineering, Tottori University, Tottori, Japan

Copyright © 2019 by author(s) and Scientific Research Publishing Inc.

This work is licensed under the Creative Commons Attribution International License (CC BY 4.0).

http://creativecommons.org/licenses/by/4.0/

Received: December 12, 2018; Accepted: February 12, 2019; Published: February 15, 2019

ABSTRACT

Flows of a rarefied gas between coaxial circular cylinders with nonuniform surface properties are studied on the basis of kinetic theory. It is assumed that the outer cylinder is a diffuse reflection boundary and the inner cylinder is a Maxwell-type boundary whose accommodation coefficient varies in the circumferential direction. Three fundamental flows are studied: 1) a flow caused by the rotation of the outer cylinder (Couette flow), 2) a flow induced between the cylinders at rest kept at different temperatures (heat transfer problem), and 3) a flow induced by the circumferential temperature distribution along the cylindrical surfaces (thermal creep flow). The linearized ES-BGK model of the Boltzmann equation is numerically analyzed using a finite difference method. The time-independent behavior of the gas is studied over a wide range of the gas rarefaction degree, the radii ratio, and a parameter characterizing the distribution of the accommodation coefficient. Due to an effect of nonuniform surface properties, a local heat transfer occurs between the gas and the cylindrical surfaces in Couette flow; a local tangential stress arises in the heat transfer problem. However, the total heat transfer between the two cylinders in Couette flow and the total torque acting on the inner cylinder in the heat transfer problem vanish irrespective of the flow parameters. Two nondegenerate reciprocity relations arise due to the effect of nonuniform surface properties. The reciprocity relations among the above-mentioned three flows are numerically confirmed over a wide range of the flow parameters. The force on the inner cylinder, which also arises due to the effect of nonuniform surface properties in Couette flow and the heat transfer problems, is studied.

Keywords:

Rarefied Gas Flow, Micro Flow, Boltzmann Equation, Gas-Surface Interaction, Reciprocity Relations

1. Introduction

Internal flows of a rarefied gas have wide applications, e.g., in vacuum science and micro engineering, and have been studied extensively on the basis of kinetic theory (see Refs. [1] - [7] and the references therein). One of the characteristics of a rarefied gas flow is that the transports of mass and energy are strongly coupled due to the nonequilibrium effect. For example, in a flow between plane parallel walls caused by a pressure gradient (Poiseuille flow), or in a flow caused by a tangential motion of the wall (Couette flow), heat flows are induced in the gas simultaneously with the motion of the gas. These heat flows are, however, parallel to the wall surfaces, and do not contribute to a heat transfer between the gas and the wall surfaces.

On the other hand, another characteristic of a rarefied gas flow is that it is affected by the gas-surface interaction law (surface properties, in short). Therefore, a nonuniformness of the surface properties can cause special phenomena. Examples of the flows between plane parallel walls with unequal surface properties are studied, e.g., in Refs. [8] - [13]. Flows between parallel walls with nonuniform surface properties along the wall surface are studied, e.g., in Refs. [14] [15] [16] [17] [18]. In Refs. [15] [17] , it is assumed that one wall is a diffuse reflection boundary and the other wall is a Maxwell-type boundary whose accommodation coefficient varies periodically in the flow direction. The following four flows are studied: a flow induced by a temperature distribution along the wall surfaces (thermal creep flow), a flow induced between the walls at rest kept at different temperatures (heat transfer problem), and the above-mentioned Poiseuille and Couette flows. In this setting, the flow is no longer parallel, and thus, a heat transfer between the gas and the wall surfaces occurs due to the combined effects of gas rarefaction and nonuniform surface properties. In fact, a tangential motion of the wall causes a local heat transfer on the wall surface in Couette flow. In the heat transfer problem, on the other hand, a temperature difference between the walls induces a gas motion and a tangential stress on the wall surface. One of the interesting results in Ref. [17] is that the integrals of these heat flux and the tangential stress over one period of the variation of the accommodation coefficient vanish irrespective of the gas rarefaction degree and the distribution of the accommodation coefficient. To the best of the author’s knowledge, this vanishing property has not yet been derived from the Boltzmann equation and the theory of reciprocity [19] [20]. It is possible that this property may have some connection with “the nonexistence of one-way flow’’ in a pipe or channel by a periodic temperature distribution in rarefied gas dynamics [7]. Naturally, a question arises whether this vanishing property is a special characteristic limited to the plane channel flow or it has a universality over a wider class of flow problems.

In this paper, we study slow rarefied gas flows between coaxial circular cylinders with nonuniform surface properties on the basis of kinetic theory. In accordance with Refs. [15] [17] , we assume that the outer cylinder is a diffuse reflection boundary and the inner cylinder is a Maxwell-type boundary whose accommodation coefficient varies in the circumferential direction. The following three flows are studied: 1) a flow caused by the rotation of the outer cylinder (Couette flow), 2) a flow induced between the cylinders at rest kept at different temperatures (heat transfer problem), and 3) a flow induced by a circumferential temperature distribution along the cylindrical surfaces at rest (thermal creep flow). The linearized ES-BGK model [21] [22] of the Boltzmann equation is numerically analyzed using a finite difference method. The time-independent behavior of the gas is studied over a wide range of the gas rarefaction degree, the radii ratio, and a parameter characterizing the distribution of the accommodation coefficient. Special attention is focused on the heat flux and the tangential stress on the cylindrical surfaces that arise due to the effect of nonuniform surface properties, and on the resultant heat transfer between and the torque acting on the cylinders. The goals of this paper are as follows. First, it is clarified whether or not the net heat transfer in Couette flow and the torque in the heat transfer problem vanish irrespective of the flow parameters. Second, three reciprocity relations among the above-mentioned three flows are numerically confirmed over a wide range of the flow parameters. Note that in a conventional Couette flow and heat transfer problems in which the accommodation coefficient is uniform, it is obvious that the flow is axially symmetric and the cylinders are subject to no force. The third goal is therefore to study the force acting on the inner cylinder arising from the effect of nonuniform surface properties. To the best of the author’s knowledge, studies on relations among several flows between boundaries with tangentially nonuniform surface properties are not available except in Refs. [17] [18] and in the present paper. In this paper, we use the simplest model, i.e., the Maxwell-type boundary condition, to model the nonuniform surface property. For, this model is characterized by only one parameter (accommodation coefficient), and thus convenient to study over various distributions of the surface property and other parameters. Incidentally, the flow problems considered in this paper are simple models of micro lubrication and heat transfer between a journal and a bearing in which the journal is partially soiled or fabricated with various materials. Studies on these flows are therefore technically valuable in designing a micro lubrication and a heat transfer in micro-electro-mechanical systems (MEMS).

This paper is organized as follows. The problem, basic equations, and basic properties of the solutions are described in Section 2. The numerical method, computational condition, and the results of accuracy tests are summarized in Section 3. The results and the discussion are presented in Section 4. Finally, a conclusion is given in Section 5.

2. Problem and Basic Equation

2.1. Problem

Let us consider a rarefied gas between the coaxial circular cylinders placed at X 1 2 + X 2 2 = r 0 2 and X 1 2 + X 2 2 = r 1 2 , where ( X 1 , X 2 , X 3 ) is a Cartesian coordinate system and r 0 and r 1 are constants such that 0 < r 0 < r 1 (Figure 1(a)). The outer cylinder is a diffuse reflection boundary, whose temperature is denoted by T w 1 . The inner cylinder is a Maxwell-type boundary at rest whose accommodation coefficient α varies in the circumferential direction. To be specific, by introducing the cylindrical coordinates ( r , θ , z ) defined by

X 1 = r cos θ , X 2 = r sin θ , X 3 = z ( 0 θ < 2 π ) , (1)

α is the function of the azimuthal angle θ :

α ( θ ) = { α 2 + α 1 2 + α 2 α 1 2 cos θ 2 C ( 0 θ 2 π C ) , α 2 + α 1 2 + α 2 α 1 2 cos 2 π θ 2 ( 1 C ) ( 2 π C < θ < 2 π ) . (2)

Here, α 1 , α 2 , and C are constants such that 0 α 1 1 , 0 α 2 1 , 0 < C < 1 (see Figure 1(b)). The temperature of this wall is denoted by T w 0 . We study the following three flows: 1) a flow caused by the rotation of the outer cylinder (Couette flow), where the angular velocity is Ω and the cylinders are kept at the uniform temperature T w 0 = T w 1 = T 0 ; 2) a flow induced between the cylinders at rest kept at different temperatures T w 0 = T 0 and T w 1 = T 0 ( 1 + Δ HT ) (heat transfer problem); and 3) a flow induced by the temperature distribution T w 0 = T w 1 = T 0 ( 1 + Δ TC X 1 / r 0 ) along the cylindrical surfaces at rest, where T 0 , Δ HT , and Δ TC are constants. The time-independent behavior of the gas is studied on the basis of kinetic theory.

Figure 1. Schematic view of the system. (a) Schematic view from the X 3 direction and (b) distribution of the accommodation coefficient α ( θ ) along the surface of the inner cylinder (Equation (2)). In panel (a), the white part of the surface of the inner cylinder represents the interval in which α is close to the minimum value α 1 . In panel (b), the distribution α ( θ ) is plotted for C = 0.25 .

In this paper, we assume that the gas behavior is governed by the ES-BGK model of the Boltzmann equation with the Prandtl number 2/3. It is also assumed that the causes of the flows are so small that the equation and the boundary condition may be linearized. To be specific, we assume that r 1 | Ω | / V 0 1 , | Δ HT | 1 , and | r 1 r 0 1 Δ TC | 1 , where V 0 = ( 2 k B T 0 / m ) 1 / 2 , k B is the Boltzmann constant, and m is the mass of a molecule.

2.2. Basic Equation

Let us take r 0 , V 0 , and ρ 0 as the reference length, the reference speed, and the reference density of the system, where ρ 0 is the mean density of the gas over the annulus domain r 0 < r < r 1 , 0 < z < r 0 . We use the dimensionless variables r ^ and z ^ for the spatial coordinates r and z , and ζ ι ( ι = r , θ , z ) for the molecular velocity components ξ ι [7] :

r ^ = r r 0 , z ^ = z r 0 , ζ ι = ξ ι V 0 . (3)

The dimensionless molecular velocity is also denoted by ζ . The perturbation ϕ of the velocity distribution function f from the equilibrium state at rest with the density ρ 0 and the temperature T 0 is defined by

ϕ = f ρ 0 V 0 3 E 1 , (4)

where E ( ζ ) = π 3 / 2 exp ( ζ 2 ) with ζ 2 = ζ r 2 + ζ θ 2 + ζ z 2 . Correspondingly, the perturbations of the macroscopic variables ω , u ι , τ , P , P ι γ , Q ι ( ι = r , θ , z ; γ = r , θ , z ) for the density ρ , the flow velocity v ι , the temperature T, the pressure p, the stress tensor p ι γ , and the heat flow vector q ι are defined by

ω = ρ / ρ 0 1 , u ι = v ι / V 0 , τ = T / T 0 1 , P = p / p 0 1 , P ι γ = p ι γ / p 0 δ ι γ , Q ι = q ι / ( p 0 V 0 ) , (5)

where p 0 = k B ρ 0 T 0 / m .

The linearized ES-BGK model of the Boltzmann equation with the Prandtl number 2/3 in the time-independent and axially uniform ( ϕ / z ^ = 0 ) state is written in the cylindrical coordinate system as

ζ r ϕ r ^ + ζ θ r ^ ϕ θ + ζ θ 2 r ^ ϕ ζ r ζ θ ζ r r ^ ϕ ζ θ = 2 π K n [ ϕ + ω + 2 ζ r u r + 2 ζ θ u θ + ( ζ 2 3 2 ) τ 1 2 ( ζ r 2 P r r + ζ θ 2 P θ θ + ζ z 2 P z z + 2 ζ θ ζ r P θ r ζ 2 P ) ] , (6)

ω = ϕ E d ζ , (7)

u ι = ζ ι ϕ E d ζ , (8)

τ = 2 3 ( ζ 2 3 2 ) ϕ E d ζ , (9)

P = ω + τ , (10)

P ι γ = 2 ζ ι ζ γ ϕ E d ζ , (11)

Q ι = ζ ι ( ζ 2 5 2 ) ϕ E d ζ . (12)

In (6), K n = l / r 0 is the Knudsen number; l is the mean free path of the gas in the equilibrium state at rest with the density ρ 0 and the temperature T 0 . The mean free path l is related to the viscosity μ and the heat conductivity λ of the gas as [23]

μ = π 3 p 0 V 0 l , λ = 5 π 4 k B p 0 m V 0 l .

In (7)-(12) and in what follows, d ζ = d ζ r d ζ θ d ζ z , and the range of integration with respect to ζ ι is its whole space unless otherwise stated.

The boundary conditions on the cylindrical walls are

ϕ ( ζ r ) = ( 1 α ) ϕ ( ζ r ) + α [ ( ζ 2 2 ) τ w 0 2 π ζ r * < 0 ζ r * ϕ ( ζ * ) E ( ζ * ) d ζ * ] ( r ^ = 1 , ζ r > 0 ) , (13)

ϕ = 2 ζ θ u w + ( ζ 2 2 ) τ w 1 + 2 π ζ r * > 0 ζ r * ϕ ( ζ * ) E ( ζ * ) d ζ * ( r ^ = R , ζ r < 0 ) , (14)

where α = α ( θ ) , R = r 1 / r 0 is the radii ratio, and ϕ ( ± ζ r ) denotes the abbreviation for ϕ ( r ^ , θ , ± ζ r , ζ θ , ζ z ) . The τ w0 , τ w1 , and u w are defined by

τ w 0 = τ w 1 = 0, u w = R Δ CF ( Couetteflow ) , τ w 0 = 0, τ w 1 = Δ HT , u w = 0 ( heattransferproblem ) , τ w 0 = τ w 1 / R = Δ TC cos θ , u w = 0 ( thermalcreepflow ) , (15)

where Δ CF = r 0 Ω / V 0 . The periodic condition at the cross section θ = 2 π is

ϕ ( r ^ , 0 , ζ r , ζ θ , ζ z ) = ϕ ( r ^ , 2 π , ζ r , ζ θ , ζ z ) ( ζ θ > 0 ) , (16)

ϕ ( r ^ , 2 π , ζ r , ζ θ , ζ z ) = ϕ ( r ^ , 0 , ζ r , ζ θ , ζ z ) ( ζ θ < 0 ) . (17)

Finally, from the definition of the average density ρ 0 , the solution should satisfy

0 d θ 1 R d r ^ r ^ ϕ E d ζ = 0. (18)

From (6)-(18), it is easily seen that we can seek the solution ϕ that is symmetric with respect to ζ z .

From now on, we denote the solution for Couette flow, heat transfer, and the thermal creep problems, respectively, by the subscript CF, HT, and TC, e.g., ϕ CF , if necessary. The subscript is frequently represented by the dummy letter J, e.g., ϕ J with J = CF, HT, and TC.

The boundary value problem (6)-(18) is characterized by the following dimensionless parameters

K n = l r 0 , R = r 1 r 0 , α 1 , α 2 and C . (19)

We will solve the problem over a wide range of the parameters (19). Among the parameters (19), α 1 , α 2 , and C are those characterizing the distribution of the accommodation coefficient α ( θ ) in (2). The last parameter C represents a degree of asymmetry of the distribution α ( θ ) with respect to θ = π ; the α ( θ ) is symmetric when C = 1 / 2 and the degree of asymmetry increases as | 1 / 2 C | increases.

Multiplying (6) by E ,2 ζ θ E , and ζ 2 E , and integrating with respect to ζ , we obtain the conservation of mass, momentum, and energy:

1 r ^ r ^ u r r ^ + 1 r ^ u θ θ = 0 , (20)

1 r ^ 2 r ^ 2 P θ r r ^ + 1 r ^ P θ θ θ = 0 , (21)

1 r ^ r ^ Q r r ^ + 1 r ^ Q θ θ = 0. (22)

Integrating (20)-(22) over 0 < θ < 2 π and applying the boundary condition (13) or (14), we obtain

0 u r d θ = 0, r ^ 2 0 P θ r d θ = c 1 , r ^ 0 Q r d θ = c 2 , (23)

where c 1 and c 2 are constants. The second relation of (23) is the conservation law of angular momentum.

2.3. Force, Torque, Heat Flow Rate, and the Reciprocity Relations

Once the macroscopic variables, e.g. P θ r ( r ^ , θ ) , are obtained, the force, torque, and the heat flow rate on the cylinders are obtained by the integral of the macroscopic variables as follows. The Cartesian components ( F 1 , F 2 ,0 ) of the force acting on the inner cylinder per unit length in the axial direction are given by

1 p 0 r 0 ( F 1 F 2 ) = 0 ( P r r cos θ P θ r sin θ P r r sin θ + P θ r cos θ ) d θ ( r ^ = 1 ) . (24)

For the convenience of the following discussion, we define the normalized torque T J ( r ^ ) acting on and the heat flow rate Q J ( r ^ ) onto the surface r ^ = r ^ , 0 θ < 2 π from the outer per unit length in the z ^ direction:

T J ( r ^ ) = 1 Δ J 0 r ^ 2 P θ r J ( r ^ , θ ) d θ , (25)

Q J ( r ^ ) = 1 Δ J 0 r ^ Q r J ( r ^ , θ ) d θ , (26)

where J = CF, HT, and TC. The T J and Q J are independent of r ^ due to the conservation laws (23), so that the position r ^ of evaluation may be omitted. The dimensional torque and the heat flow rate are, respectively, given by p 0 r 0 2 Δ J T J and p 0 r 0 V 0 Δ J Q J .

From the symmetric relation of the linearized Boltzmann equation [2] [19] [20] , the following three reciprocity relations hold among the three flow problems (J = CF, HT, and TC) for arbitrary values of the parameters (19):

Q CF ( R ) = T HT ( R ) , (27)

T TC ( R ) = H CF , (28)

Q TC ( R ) = H HT , (29)

where

H J = 1 Δ J 0 [ R 2 Q r J ( R , θ ) Q r J ( 1, θ ) ] cos θ d θ . (30)

Equations (27)-(29) will be numerically confirmed over a wide range of the parameters (19).

2.4. Cases of Uniform or Symmetric Distribution of the Accommodation Coefficient

When the accommodation coefficient α ( θ ) in the boundary condition (13) is uniform (conventional case, in short), the solutions of the boundary value problem (6)-(18) are simplified as follows. The solutions for the conventional cases are reported, e.g., in Refs. [24] [25] for Couette flow and Ref. [26] for the heat transfer problem.

Couette flow and the heat transfer problems reduce to axially symmetric problems, i.e., ϕ / θ = 0 , in the conventional case. In Couette flow problem, further, it is easily seen that the distribution function ϕ is an odd function with respect to ζ θ . Consequently from (7)-(12), we directly obtain

ω CF = u r CF = τ CF = P r r CF = P θ θ CF = P z z CF = Q r CF = 0, (31)

and thus

F 1 CF = F 2 CF = Q CF = H CF = 0. (32)

In the heat transfer problem, similarly, the distribution function ϕ is an even function with respect to ζ θ . By using the conservation laws (23) and the boundary condition (13) or (14), we obtain

u r HT = u θ HT = P θ r HT = Q θ HT = 0, (33)

and

F 1 HT = F 2 HT = T HT = H HT = 0. (34)

In the thermal creep flow problem, in contrast, the boundary value problem remains a spatially two dimensional one for r ^ and θ . Because of the special form of the temperature distribution τ w 0 and τ w 1 in (15), however, we find that the solution ϕ = ϕ e cos θ + ϕ o sin θ is compatible with the boundary value problem (6)-(18), where ϕ e and ϕ o are independent of θ and, respectively, even and odd functions with respect to ζ θ . From this form, we see that P θ r TC / sin θ , P r r TC / cos θ , and Q r TC / cos θ are independent of θ , and thus

F 2 TC = T TC = Q TC = 0. (35)

To summarize, when the accommodation coefficient is uniform, both sides of all the reciprocity relations (27)-(29) vanish; obviously, the cylinders are subject to no force in Couette and heat transfer problems. Probably because of this reason, the reciprocity relations are rarely discussed in the flow problems between coaxial circular cylinders. In the case of nonuniform accommodation coefficient to be studied in this paper, in contrast, (32), (34), and (35) do not generally hold, and all the three reciprocity relations (27)-(29) are of interest. The numerical confirmation of the reciprocity relations (27)-(29) is one of the important goals of this paper.

Special attention is focused on the quantities in (27). In Refs. [15] [17] , flow problems between plane parallel walls in which one of the walls has a periodic distribution of the accommodation coefficient were studied. In Ref. [17] , due to the effect of the nonuniform accommodation coefficient of the solid boundary, a local heat transfer across the wall surface occurs in plane Couette flow, and a local tangential stress on the surface occurs in the plane heat transfer problem. The numerical solution in Ref. [17] showed, however, that the integrals of the heat flux and the tangential stress over one period of the distribution of the accommodation coefficient vanish irrespective of the flow parameters (see Section 1). If this phenomenon has a universality beyond the flows between plane parallel walls, both sides of (27) may vanish. To clarify this point is an important motivation of this study.

Let us discuss another special case that the accommodation coefficient α ( θ ) is nonuniform but symmetric with respect to θ = π (symmetric case, in short). This is the case of C = 1 / 2 in (2). In this case, all the reciprocity relations (27)-(29) degenerate as follows. The derivation is quite similar to that in Ref. [15] and thus omitted. In Couette flow problem, in this case, the macroscopic variables u θ CF , P θ r CF , and Q θ CF are symmetric with respect to θ = π , whereas ω CF , u rCF , τ CF , P rrCF , and Q r CF are antisymmetric. Therefore,

F 1 CF = Q CF = H CF = 0. (36)

In the heat transfer and the thermal creep problems, ω J , u r J , τ J , P r r J , and Q r J (J = HT and TC) are symmetric with respect to θ = π , whereas u θ J , P θ r J , and Q θ J are antisymmetric. Therefore,

F 2 J = T J = 0 ( J = HT and TC ) . (37)

That is, all the reciprocity relations (27)-(29) degenerate also in the symmetric case. In this paper, therefore, the cases other than C = 1 / 2 are included.

3. Numerical Analysis

3.1. Numerical Method

A few preliminary processes are explained. First, the new variables ζ ρ and θ ζ are introduced in place of ζ r and ζ θ [7] :

ζ r = ζ ρ cos θ ζ , ζ θ = ζ ρ sin θ ζ ( ζ ρ > 0 , π < θ ζ π ) . (38)

Then, the number of derivative terms in (6) is reduced from four to three, which is convenient for the numerical analysis. Second, by introducing the marginal distribution functions Φ and Θ defined by [7]

Φ = 1 π ϕ exp ( ζ z 2 ) d ζ z , Θ = 2 π ζ z 2 ϕ exp ( ζ z 2 ) d ζ z , (39)

the molecular velocity component ζ z can be eliminated from the boundary value problem (6)-(18). Then, we are led to a boundary value problem for the four independent variables r ^ , θ , ζ ρ , and θ ζ .

In the numerical analysis, the infinite interval ζ ρ > 0 is replaced by a finite one 0< ζ ρ ζ D , where ζ D is a constant. The computational grid system is arranged in the four-dimensional space 1 r ^ R , 0 θ 2 π , 0 < ζ ρ ζ D , and π θ ζ π . Then, a finite difference scheme of the second order with respect to r ^ , θ , and θ ζ is constructed. The time-independent solution is obtained by an iteration method.

A few remarks on the numerical method are given below. First, the solution ϕ of the boundary-value problem is discontinuous across the surfaces

r ^ sin θ ζ = ± 1, | θ ζ | π / 2 , 0 θ 2 π , (40)

in the ( r ^ , θ , θ ζ ) space. Therefore, a conventional finite difference scheme cannot be applied directly. In this study, the hybrid scheme of finite difference and characteristic coordinate methods devised in Refs. [27] [28] is adopted. Incidentally, if the outer cylinder is the Maxwell-type boundary as well, the discontinuity staring from the inner cylinder reflects on the outer cylinder, and the discontinuous surfaces are present also in the region of | θ ζ | > π / 2 . In this case, the construction of the numerical scheme is more difficult. On the other hand, if the inner cylinder rotates, the problem is obviously a time-dependent one and elaborate to analyze. This is the reason why we chose the situation that the outer cylinder is a diffuse reflection boundary and the inner cylinder is at rest in this paper. Second, it is difficult to obtain an accurate solution at small Knudsen numbers because of the numerical error arising from the collision term [5]. To avoid this difficulty, we subtract the asymptotic solution [5] [23] for small Knudsen numbers from the numerical solution in calculating the collision term, or (7)-(12). In this paper, the asymptotic solution is constructed approximately using the numerical values of the macroscopic variables. This subtraction process is adopted when K n / ( R 1 ) 1 .

3.2. Computational Condition

The computational grid points are arranged as follows. For the coordinate r ^ , N r grid points are arranged in the interval 1 r ^ R , where grid size is uniform ( = d r 2 ( R 1 ) ) in | ( r ^ 1 ) / ( R 1 ) 1 / 2 | < 0.34 and nonuniform otherwise; the grid size is minimum ( = d r 1 ( R 1 ) ) at r ^ = 1 and R. For the coordinate θ , 101 grid points are arranged uniformly in the interval 0 θ 2 π . For the coordinate ζ ρ , 31 grid points are arranged uniformly in 0 ζ ρ ζ D , where we chose ζ D = 5 . For the coordinate θ ζ , N θ grid points are arranged uniformly in π θ ζ π . Here, N r , N θ , d r 1 , and d r 2 are constants chosen depending on the parameters (19) as follows: ( N r , N θ , d r 1 , d r 2 ) = ( 1281 , 401 , 2.7 × 10 5 , 1.3 × 10 3 ) when K n 0.2 and R = 5 , or ( N r , N θ , d r 1 , d r 2 ) = ( 641 , 641 , 5.4 × 10 5 , 2.5 × 10 3 ) otherwise.

3.3. Results of Accuracy Tests

The accuracy of the numerical solutions is tested as follows. The following tests are conducted over all the combinations of R = 1.2 , 2 , 5 , C = 0.1 , 0.25 , α 1 = 0.5 , α 2 = 1 , and K n = 0.1 , 1 , 10 of the parameters (19).

1) Conservation law. The solution of the boundary value problem (6)-(18) should satisfy the conservation laws (23). That is, the left-hand sides in (23) should be independent of r ^ . Therefore, the maximum variation of the left-hand side in the numerical solution in 1 r ^ R serves to estimate the numerical error. The maximum variations of the left-hand side of the second relation of (23) in Couette and thermal creep problems are, respectively, less than 0.63% and 0.45% relative to the value at r ^ = 1 . The maximum variations of the left-hand side of the third relation in heat transfer and thermal creep problems are, respectively, less than 0.072% and 2.4% relative to its value at r ^ = 1 . As we will see in Section 4, the left-hand side of the second relation in heat transfer problem and that of the third relation in Couette flow vanish ( c 1 = c 2 = 0 ), so that the relative error cannot be defined. Instead, the maximum variations of these two quantities relative to the maximum of the respective integrands r ^ 2 | P θ r HT | and r ^ | Q r CF | over the cylindrical surfaces are, respectively, 0.008% and 0.13%. The first relation of (23) is not used for the error estimate because the integrand u r vanishes identically on the cylindrical surfaces by (13) and (14).

2) Re-computation using a different computational grid system. For a test of accuracy, we also conducted re-computations using a coarser grid system. In the coarser system, the grid sizes in r ^ ,θ ,and θ ζ are two times coarser, that in ζ ρ is approximately 1.4 times coarser, and the upper limit ζ D of ζ ρ is reduced to ζ D = 4.5 , simultaneously. The error of the numerical solution using the grid system of Section 3.2 is estimated from the difference between the two numerical solutions. The estimated numerical error in the values of the normalized torque T CF and T TC (Equation (25)) at r ^ = 1 is 0.076% and 0.35%, respectively; the error in the normalized heat flow rate Q HT and Q TC (Equation (26)) at r ^ = 1 is 0.0072% and 0.36%, respectively. The other quantities Q CF and T HT will be shown to vanish in Section 4. The numerical error in the force ( F 1 , F 2 ) (Equation (24)) is estimated to be less than 0.14%, 0.32%, and 0.44% in Couette flow, heat transfer, and thermal creep problems, respectively.

The reciprocity relations (27)-(29) also serve to estimate the accuracy of the numerical solutions. The result will be presented in Section 4.4.

It may be noted that the accuracy in the numerical solution, e.g., the value of Q TC , in this paper is worse than that in our previous paper of the flow between plane parallel walls [17]. This is probably because that the basic Equation (6) is complicated (three derivative terms in place of two in Ref. [17] ), and that the substantial Knudsen number l / ( r 1 r 0 ) based on the gap size r 1 r 0 can be very small when r 1 / r 0 = R is large. At small Knudsen numbers, it is difficult to obtain an accurate numerical solution not only because the numerical error arising from the collision term is non-negligible [5] , but also because the quantities of interest are very small. For example, the worst case of the energy conservation law (2.4%) in the item 1) occurs in the thermal creep flow at ( K n , R , C ) = ( 0.1 , 5 , 0.1 ) . In this case, the normalized energy flow rate Q TC is 2.9 × 10 3 , which is extremely small. Note that the normalized distribution function ϕ TC / Δ TC is of the order of unity whereas the heat flow rate Q TC is extremely small due to a drastic cancellation in the successive integrals (12) and (26).

4. Results and Discussion

In this paper, the numerical computation is conducted for all the combinations of the parameters (19) as follows: R = 1.2 , 2 , 5 , C = 0.1 , 0.25 , 0.5 , α 1 = 0.5 , α 2 = 1 , and 15 Knudsen numbers K n ranging between 0.1 and 10. The results of the numerical analysis are presented in this section. In the following discussion, we also use the dimensionless Cartesian coordinates x 1 = X 1 / r 0 = r ^ cos θ and x 2 = X 2 / r 0 = r ^ sin θ , and the vector components

u 1 = u r cos θ u θ sin θ , u 2 = u r sin θ + u θ cos θ , (41)

and so on. Further, it is assumed that Δ CF , Δ HT , and Δ TC in (15) are positive without loss of generality.

4.1. Couette Flow

First, some examples of the flow field of Couette flow are presented in Figure 2. The flow velocity and the heat flow for the radii ratio R = 2 and the Knudsen numbers K n = 0.1 , 1 , and 10 are presented in Figures 2(a)-(c), respectively; those for R = 5 and the same K n ’s are presented in Figures 2(d)-(f), respectively ( α 1 = 0.5 , α 2 = 1 , C = 0.25 ). The arrow represents the normalized flow velocity ( u 1 CF , u 2 CF ) / Δ CF at its starting point, the magnitude of which is shown by the arrow at the upper right corner of each panel. The dotted line represents the “stream line’’ of the heat flow vector ( Q 1 CF , Q 2 CF ) , i.e., the isoline of the function Ψ defined by

Q 1 CF = Ψ x 2 , Q 2 CF = Ψ x 1 . (42)

The accommodation coefficient α ( θ ) is minimum at θ = π / 2 when C = 0.25 ; the white part along the surface of the inner cylinder represents the part in which α is close to the minimum ( 0.5 α 0.55 ). Now let us examine

Figure 2. Flow fields of Couette flow ( α 1 = 0.5 , α 2 = 1 , C = 0.25 ). (a, b, c) R = 2 , (d, e, f) R = 5 , (a, d) K n = 0.1 , (b, e) K n = 1 , and (c, f) K n = 10 . Arrow: the normalized flow velocity ( u 1 CF , u 2 CF ) / Δ CF at its starting point, the magnitude of which is represented by the arrow at the upper right corner of each figure. Dotted line: “stream line’’, or isoline Ψ = const (Equation (42)), of the heat flow vector; the arrow on the dotted line represents the direction of the heat flow vector. The white part along the surface of the inner cylinder represents the interval of 0.5 α 0.55 .

the results for R = 2 , focusing our attention on the heat flow field. The rotation of the outer cylinder causes the motion of the gas and the heat flow simultaneously. When the Knudsen number is very small, according to the asymptotic theory of the Boltzmann equation [5] [7] , the heat flow is present mainly in the Knudsen layers; the heat flow is a clockwise circulation flow along the outer cylinder and is a counterclockwise one along the inner cylinder. The heat flow is enhanced on passing the white part of small accommodation coefficient [29]. Then, the direction of the heat flow slightly shifts from the tangential direction due to the conservation on energy. As the Knudsen number increases (Figure 2(a)), the heat flow spreads into the bulk region, and forms a roll over the white part. This roll of the heat flow is not special to the present cylindrical Couette flow, but is observed also in the plane Couette flow in Ref. [17]. As the Knudsen number increases further (Figure 2(b)), the outer part of the roll is absorbed by the outer cylinder, and an outflow and an inflow of heat across the cylindrical surface occur. When the Knudsen number is sufficiently large (Figure 2(c)), the roll of the heat flow extends and merges at approximately θ = 7 π / 4 , and a one-way flow of heat from one part to the other part of the outer cylinder is established; the direction is approximately ( x 1 , x 2 ) = ( 1 , 1 ) . When the radii ratio is large ( R = 5 , Figures 2(d)-(f)), the behavior of the gas is qualitatively similar. For a small Knudsen number (Figure 2(d)), the heat flow is almost absent near the outer cylinder. This is because the breaking effect on the gas motion by the inner cylinder is small, and thus the gas flow near the outer cylinder is nearly a rigid-body rotation; therefore, the velocity distribution function is nearly Maxwellian, and thus the heat flow vanishes.

From the flow fields in Figure 2, one expects nonuniform distributions of the tangential stress and the heat flux along the cylindrical surfaces, which are examined in Figure 3 in detail. The dimensionless tangential stress P θ r CF along the cylindrical surfaces at r ^ = 1 and r ^ = R for the case of Figure 2(b) ( K n = 1 , R = 2 ) is shown in Figure 3(a) as a function of the azimuthal angle
θ . In this figure, r ^ 2 P θ r CF is plotted because P θ r decays approximately P θ r r ^ 2 due to the conservation law (23). Note that P θ r CF > 0 means that the tangential stress acts in the + θ ( θ ) direction on the inner (outer) cylinder. The magnitude of the tangential stress on the inner cylinder is minimum at approximately θ = π / 2 , i.e., the minimum position of α ( θ ) . The dotted line represents the average ( 2 π Δ CF ) 1 0 r ^ 2 P θ r CF d θ = ( ) 1 T CF , where T CF is the normalized torque defined by (25). The dotted lines at r ^ = 1 and r ^ = R are indistinguishable according to the conservation law (23). Incidentally from Figure 3(a), we can estimate that the inner cylinder is subject to a force (24) from the gas as follows. The tangential stress in the hemicircle of approximately π / 4 < θ < / 4 is less than the average. Therefore, the contribution to the force from this hemicircle and that from the other do not cancel out. Consequently, a force approximately in the direction of ( x 1 , x 2 ) = ( 1 , 1 ) , or θ = π / 4 , is expected. The force acting on the inner cylinder is discussed in detail in Section 4.6.

Similarly, the normal component Q r CF of the dimensionless heat flow vector on the cylinders for the same case (Figure 2(b)) is presented in Figure 3(b). Note that r ^ Q r CF is plotted (Equation (23)), and that Q r CF > 0 means a heat flow from the cylinder to the gas (from the gas to the cylinder) on the inner (outer) cylinder. Obviously, there is a distribution of inflow and outflow of heat along the cylindrical surfaces. The heat transfer is significant on the outer cylinder, which is easily seen from the flow field in Figure 2(b). The dotted line represents the average ( 2 π Δ CF ) 1 0 r ^ Q r CF d θ = ( ) 1 Q CF , where Q CF is the

Figure 3. Distribution of the dimensionless tangential stress P θ r CF and the normal component Q r CF of the heat flow along the cylindrical surfaces in Couette flow ( K n = 1 , R = 2 , α 1 = 0.5 , α 2 = 1 , C = 0.25 ; cf. Figure 2(b)). Solid line in (a): r ^ 2 P θ r CF / Δ CF ; solid line in (b): r ^ Q r CF / Δ CF . Dotted line in (a): the average ( 2 π Δ CF ) 1 0 r ^ 2 P θ r CF d θ at r ^ = 1 and r ^ = R ; dotted line in (b): the average ( 2 π Δ CF ) 1 0 r ^ Q r CF d θ at r ^ = 1 and r ^ = R ; the dotted lines at r ^ = 1 and r ^ = R are indistinguishable (Equation (23)).

normalized heat flow rate defined by (26). The dotted lines at r ^ = 1 and r ^ = R are indistinguishable according to the conservation law (23). Note that the average, or Q CF , is very small; to be specific, ( ) 1 Q CF at r ^ = 1 and r ^ = R is less than 0.004% of the maximum of the integrand | r ^ Q r CF / Δ CF | ( = 1.3 × 10 2 ) over the cylindrical surfaces. This property is common to every set of the values of the parameters (19); the average ( ) 1 Q CF is less than 0.15% of the maximum of | r ^ Q r CF / Δ CF | over all the cases of the parameters (19) conducted in this paper. In view of the accuracy tests in Section 3.3, these small values are of the same order as the inevitable numerical error. Therefore, we may conclude that

Q CF = 0 (43)

irrespective of the values of the parameters (19). That is, although a local heat transfer between the gas and the cylindrical surfaces occurs due to the effect of nonuniform accommodation coefficient, the total heat transfer between the two cylinders vanishes identically. Note that, nevertheless, a one-way heat flow from one part to the other part of the outer cylinder is present, as shown in Figure 2 and Figure 3(b).

4.2. Heat Transfer Problem

Next, some examples of the flow field of the heat transfer problem are presented in Figure 4. The parameters (19) in Figures 4(a)-(f) are the same as those in Figures 2(a)-(f). In the figures, the isothermal line τ HT = const is represented by the solid line. The meanings of other lines and symbols are similar to those in Figure 2. Due to an effect of the nonuniform accommodation coefficient, a flow is induced between the cylinders at rest kept at different temperatures. The mechanism of this flow is similar to that in the plane heat transfer problem [17].

Figure 4. Flow fields of the heat transfer problem ( α 1 = 0.5 , α 2 = 1 , C = 0.25 ). Arrow: the normalized flow velocity ( u 1 HT , u 2 HT ) / Δ HT ; dotted line: “stream line’’ of the heat flow vector ( Q 1 HT , Q 2 HT ) ; solid line: isothermal lines τ HT = 0.8 , 0.6 , from the outer. See the caption of Figure 2.

Because of the nonuniformness of the accommodation coefficient, the temperature jump of the gas at the cylindrical surfaces is nonuniform. Consequently, a temperature distribution of the gas is formed along the cylindrical surfaces; the temperature of the gas on the inner cylinder is the highest near the white part of small accommodation coefficient. Then, gas flows are induced along the cylindrical surfaces toward the hot part by the same mechanism as that of the thermal creep flow. Two major flows from both sides along the inner cylinder collide at approximately θ = 3 π / 4 , and finally two rolls are formed in the annulus; the one in approximately 3 π / 4 < θ < / 4 is counterclockwise and the other is clockwise. Incidentally, from Figure 4, it is expected that the inner cylinder is subject to a force approximately in the direction of ( x 1 , x 2 ) = ( 1 , 1 ) by the reaction of these thermal-creep type flows. The force on the inner cylinder will be presented in Section 4.6.

The distributions of the tangential stress P θ r HT and the normal component Q r HT of the heat flow vector along the cylindrical surfaces for the case of Figure 4(b) are presented in Figure 5(a) and Figure 5(b), respectively. The dotted lines represent the averages ( 2 π Δ HT ) 1 0 r ^ 2 P θ r HT d θ = ( 2 π ) 1 T HT and ( 2 π Δ HT ) 1 0 r ^ Q r HT d θ = ( 2 π ) 1 Q HT , respectively (see (25) and (26)). The average value of the tangential stress, or the normalized torque T HT , is very small, and this property is common to every set of the values of the parameters (19). To be specific, the value ( ) 1 T HT is less than 0.03% of the maximum of the integrand | r ^ 2 P θ r HT | over the cylindrical surfaces for every set of the values of the parameters (19). Therefore, following the same criterion as that given in the discussion on Figure 3(b), we concluded that

T HT = 0 (44)

irrespective of the values of the parameters (19). That is, although a local tangential stress is induced on the cylinders, the total torque acting on the inner cylinder vanishes identically irrespective of the parameters (19). Equation (44) is a counterpart of (43).

4.3. Thermal Creep Flow

Similarly, some examples of the flow field of the thermal creep flow are presented in Figure 6. Note that the distributions of the dimensionless wall temperature τ w 0 and τ w 1 are τ w 0 = τ w 1 = Δ TC x 1 in the Cartesian coordinates. Thermal creep flows are induced by this temperature distribution along the four parts of the cylindrical surfaces: the x 2 > 0 parts of the surfaces of the outer and inner cylinders, and the x 2 < 0 parts of those of the outer and inner cylinders. These flows interact, and a flow field of multi rolls is formed. Note

Figure 5. Distribution of the dimensionless tangential stress P θ r HT and the normal component Q r HT of the heat flow along the cylindrical surfaces in the heat transfer problem ( K n = 1 , R = 2 , α 1 = 0.5 , α 2 = 1 , C = 0.25 ; cf. Figure 4(b)). Solid line in (a): r ^ 2 P θ r HT / Δ HT ; solid line in (b): r ^ Q r HT / Δ HT . Dotted line: the average at r ^ = 1 and r ^ = R . See the caption of Figure 3.

Figure 6. Flow fields of the thermal creep flow ( α 1 = 0.5 , α 2 = 1 , C = 0.25 ). Arrow: the normalized flow velocity ( u 1 TC , u 2 TC ) / Δ TC ; dotted line: “stream line’’ of the heat flow vector ( Q 1 TC , Q 2 TC ) . See the caption of Figure 2.

that these thermal creep flows are induced even in the conventional case. The effect of the nonuniform accommodation coefficient works to make the flow asymmetric with respect to x 2 = 0 . The degree of asymmetry is, however, not so strong. The heat flow, denoted by the dotted line, is basically a one-way flow from the hotter side to the colder side of the outer cylinder.

The distributions of the tangential stress P θ r TC and the normal component Q r TC of the heat flow vector along the cylindrical surfaces for the case of Figure 6(b) are presented in Figure 7(a) and Figure 7(b), respectively. Note that P θ r TC / sin θ and Q r TC / cos θ are independent of θ in the conventional case (cf. Section 2.4). A weak deviation from this θ-dependence is observed in Figure 7(a) and Figure 7(b). The dotted lines represent the averages ( 2 π Δ TC ) 1 0 r ^ 2 P θ r TC d θ = ( ) 1 T TC and ( 2 π Δ TC ) 1 0 r ^ Q r TC d θ = ( 2 π ) 1 Q TC , respectively (see (25) and (26)). Note that the average values are small but do not vanish. To be specific, ( ) 1 T TC = 5.394 × 10 3 at r ^ = 1 and 5.392 × 10 3 at r ^ = R in Figure 7(a), and ( ) 1 Q TC = 2.112 × 10 2 at r ^ = 1 and 2.105 × 10 2 at r ^ = R in Figure 7(b). The magnitudes of these values are much larger than the estimated error (cf. Section 3.3). The values at r ^ = 1 and r ^ = R agree well; the relative differences are only 0.04% and 0.3%, respectively, which are of the same order as the numerical error. We also obtain nonzero values of T TC and Q TC for every set of the values of the parameters (19) except T TC at C = 0.5 (Equation (37)). That is, we find that both the torque T TC and the heat flow rate Q TC in the thermal creep flow do not generally vanish in contrast to (43) and (44). These two quantities appear in the reciprocity relations (28) and (29). In this way, we obtain two nondegenerate reciprocity relations by the effect of nonuniform accommodation coefficient. The reciprocity relations (28) and (29) are numerically confirmed in the next subsection. From Figure 7(a), incidentally, we can estimate the force (24) acting on the inner cylinder as follows. The contribution of the tangential stress from 0 < θ < π and that from π < θ < 2 π combine to produce a resultant force approximately in the x 1 direction. Note that, in contrast to Couette flow and heat transfer problems, a force acts even in the conventional thermal creep flow; this force is similar to the thermal force on a body with a nonuniform surface temperature, e.g., in Ref. [7]. From Figure 7(a), it is expected that the direction of the force is weakly affected by the nonuniformness of the accommodation coefficient.

4.4. Reciprocity Relations

The reciprocity relations (27), (28), and (29) are numerically confirmed now. From (43) and (44), which are numerically obtained in Sections 4.1 and 4.2, we directly have

Q CF ( R ) = T HT ( R ) = 0 (45)

irrespective of the parameters (19). We now examine the other two relations (28) and (29), which do not generally vanish, over the range of the parameters stated at the beginning of Section 4 except C = 0.5 . From the numerical results, the relative difference between the left-hand and right-hand sides of (28) is less than 1.0%; the relative difference between the left-hand and right-hand sides of (29) is less than 0.43%. In view of the accuracy tests in Section 3.3, these small differences are of the order of the numerical error. Therefore, we may conclude that the reciprocity relations (28) and (29) are numerically confirmed over the above-mentioned range of the parameters (19). It may be noted that the relative agreement in the nonvanishing equalities (28) and (29) is worse than the corresponding reciprocity equalities in the plane channel flows [17]. The reason may be attributed to the difficulty in obtaining accurate numerical solutions at small Knudsen numbers and a large radii ratio in the present paper (see the last paragraph in Section 3.3).

Figure 7. Distribution of the dimensionless tangential stress P θ r TC and the normal component Q r TC of the heat flow along the cylindrical surfaces in the thermal creep flow ( K n = 1 , R = 2 , α 1 = 0.5 , α 2 = 1 , C = 0.25 ; cf. Figure 6(b)). Solid line in (a): r ^ 2 P θ r TC / Δ TC ; solid line in (b): r ^ Q r TC / Δ TC . Dotted line: the average at r ^ = 1 and r ^ = R . See the caption of Figure 3.

When the accommodation coefficient α ( θ ) is uniform or symmetric with respect to π , both sides of all the reciprocity relations (27), (28), and (29) vanish identically (Section 2.4). In this paper, therefore, we especially considered the nonuniform and asymmetric distribution of the accommodation coefficient to avoid the degeneracy of the reciprocity relations. As a result, we obtained nondegenerate relations (28) and (29), although the magnitude is small. Nevertheless, both sides of (27) vanish irrespective of the values of the parameters (19). Note that the meaning of the second equality in (45) is completely different from that in (32), (34), (36), and (37) in the conventional or symmetric case; the integrands Q r CF and P θ r HT of Q CF and T HT do not vanish identically nor do not have any antisymmetry with respect to θ (cf. Figure 3(b) and Figure 5(a)). In the plane channel flow problem in Ref. [17] , a degeneracy corresponding to (45) was found between plane Couette flow and plane heat transfer problems. To clarify whether a similar property holds also in the flows between coaxial cylinders or not is one of the important goals of the present paper. The answer turned out to be affirmative. This result suggests that this degeneracy property may have a universality over a wider class of flows beyond the flows in a plane channel or an annulus.

In rarefied gas dynamics, a phenomenon referred to as “nonexistence of one-way flow’’ is known [7]. Consider a rarefied gas in a straight pipe with an arbitrary cross section, an arbitrary uniform accommodation coefficient, and an arbitrary periodic surface temperature distribution in the longitudinal direction. Then, a flow is induced due to a rarefaction effect. However, net mass flow is never induced through the pipe. The degeneracy properties observed in the present paper show a similarity to this phenomenon, and there may be some relation between them [17].

4.5. Torque on the Inner Cylinder and the Heat Flow Rate between the Cylinders

The torque acting on the inner cylinder and the heat flow rate between the two cylinders in the three flow problems are presented here over a wide range of the parameters. In the rest of Section 4, the Knudsen number

K n D = l r 1 r 0 = K n R 1 (46)

based on the gap size r 1 r 0 is used to present the results, because it is convenient to observe the similarity between the results for different radii ratios R. The normalized torque T CF of Couette flow is presented in Figure 8 as a function of K n D for various values of R and C. Note that the parameter C represents the degree of asymmetry of the accommodation coefficient α ( θ ) (see the explanation below (19)). The markers represent the numerical results. The torque is an increasing function of both the Knudsen number K n D and the radii ratio R. The torque is, however, almost independent of the parameter C; a similar property was observed also in plane Couette flow [17]. The solid line represents the estimate formula proposed in Ref. [14] :

T CF = η T CF 1 + ( 1 η ) T CF 2 , (47)

where T CF 1 and T CF 2 are normalized torque in the conventional Couette flow problem with α ( θ ) = const = α 1 and α ( θ ) = const = α 2 , respectively. The factor η is the constant defined by

η α 1 + ( 1 η ) α 2 = 1 0 α ( θ ) d θ . (48)

For (2), η = 1 / 2 irrespective of C. The estimate formula approximates the numerical results very well.

Next, the normalized heat flow rate Q HT between the cylinders (from the outer cylinder to the inner) in the heat transfer problem is presented in Figure 9. The heat flow rate is an increasing function of both K n D and R. A very weak dependence on the parameter C is observed also in this heat flow rate. The solid line represents the estimate formula, i.e., the counterpart of (47). The estimate formula approximates the numerical results very well.

Similarly, the normalized torque T TC and the heat flow rate Q TC in the thermal creep flow are presented in Figure 10(a) and Figure 10(b), respectively. The torque T TC and the heat flow rate Q TC vanish in the conventional case (Equation (35)). In other words, the nonzero values in Figure 10(a) and Figure 10(b) are essentially due to the effect of nonuniform accommodation coefficient. (Note that the dotted lines in Figure 7(a) and Figure 7(b) are divided by 2 π .) The torque and the heat flow rate increase with an increase in the Knudsen number. The dependence on the parameter C is perceptible. When C increases to approach 0.5, the torque T TC decreases to vanish (Equation (37)), whereas the heat flow rate Q TC increases and takes the maximum when C = 0.5 . The estimate formula, the counterpart of (47), is omitted because it vanishes identically.

4.6. Force Acting on the Inner Cylinder

Finally, the force acting on the inner cylinder in the three flow problems is presented in Figure 11. The figure shows the dimensionless magnitude | F ^ J | and the direction angle ϑ J (J = CF, HT, and TC) of the force defined by

F 1 J / ( p 0 r 0 ) = | F ^ J | cos ϑ J , F 2 J / ( p 0 r 0 ) = | F ^ J | sin ϑ J , 0.9 π < ϑ J 1.1 π (49)

Figure 8. Normalized torque T CF acting on the inner cylinder as a function of K n D , R , and C in Couette flow ( α 1 = 0.5 , α 2 = 1 ). Markers: numerical solutions; diamonds for R = 1.2 , circles for R = 2 , and squares for R = 5 ; closed symbols for C = 0.1 , symbols with plus for C = 0.25 , and open symbols for C = 0.5 (see the key in the figure). Solid line: estimate formula (47) for R = 1.2 , 2 , and 5 from the bottom.

Figure 9. Normalized heat flow rate Q HT between the cylinders as a function of K n D , R , and C in the heat transfer problem ( α 1 = 0.5 , α 2 = 1 ). Markers: numerical solutions (see the key in the figure). Solid line: estimate formula (cf. (47)) for R = 1.2 , 2 , and 5 from the bottom.

Figure 10. Normalized torque T TC acting on the inner cylinder and the heat flow rate Q TC between the cylinders as a function of K n D , R , and C in the thermal creep flow ( α 1 = 0.5 , α 2 = 1 ). (a) Normalized torque T TC and (b) the heat flow rate Q TC . Markers: numerical solution (see the key in the figure). The T TC vanishes when C = 0.5 (Equation (37)).

Figure 11. Dimensionless force acting on the inner cylinder as a function of K n D , R , and C ( α 1 = 0.5 , α 2 = 1 ). (a, b, c) The magnitude | F ^ J | and (d, e, f) the direction ϑ J (Equation (49)). (a, d) Couette flow, (b, e) heat transfer problem, and (c, f) thermal creep flow. Markers: numerical solution (see the key in Figure 11(a)). Solid line in Figure 11(c) and Figure 11(f): estimate formula (cf. (47)) for R = 1.2 , 2 , and 5 from the bottom.

as functions of the Knudsen number K n D . In the conventional Couette flow and the heat transfer problems, the flow is axially symmetric and, obviously, the inner cylinder is subject to no force. That is, the forces in these two problems, which are shown in Figures 11(a)-(e), arise essentially due to the effect of nonuniform accommodation coefficient. The magnitude | F ^ J | is an increasing function of the Knudsen number K n D ; the dependence on the parameter C is perceptible (Figure 11(a) and Figure 11(b)). In contrast, the direction ϑ J is only weakly dependent on K n D and R, and is determined almost solely by C. In Couette flow at C = 0.25 (Figure 11(d)), the direction is approximately ϑ CF = π / 4 , which agrees with the estimate given in the discussion on Figure 3(a). As C increases and approaches 0.5, the position θ of the minimum of the accommodation coefficient approaches π , and thus the direction ϑ CF of the force increases and approaches ϑ CF = π / 2 (Equation (36)). The direction ϑ HT in the heat transfer problem increases to approach zero from a negative value as C increases and approaches 0.5 (Equation (37)). This behavior may also be understood from Figure 4.

In the thermal creep flow in Figure 11(c) and Figure 11(f), in contrast, the cylinder is subject to a force even in the conventional case. The influence of the nonuniform accommodation coefficient on the direction ϑ TC of the force is weak (Figure 11(f)), in agreement with the estimate in the discussion on Figure 7(a). The direction ϑ TC is nearly equal to π , which is the direction in the conventional case (Equation (35)). The solid line represents the estimate formula corresponding to (47). The estimate formula approximates the numerical results very well.

Couette flow considered in this paper is a simple model of micro lubrication between a journal and a bearing in which the journal is partially soiled or fabricated with various materials. As shown in Figure 11, a force acting on the inner cylinder arises due to the effect of nonuniform surface properties of the journal. Note that the magnitude of the force is not weak. For example, the magnitude of the force is approximately 95% of the torque divided by the radius of the inner cylinder when K n D = 1 , R = 1.2 , and C = 0.25 ; the magnitude of the force increases as the clearance between the cylinders shrinks (or R decreases). This phenomenon poses a serious problem in micro engineering, because this force may cause an eccentricity of the journal and may affect the lubrication process. In this study, however, the model is not sufficiently good for some reasons, e.g., the inner cylinder is at rest, and the gap is not so narrow. A more realistic study focusing on the engineering applications, which is beyond the scope of this paper, will be conducted in a separate paper.

5. Conclusion

In this paper, we studied the flows of a rarefied gas between coaxial circular cylinders with nonuniform surface properties on the basis of kinetic theory. We assumed that the outer cylinder is a diffuse reflection boundary and the inner cylinder is a Maxwell-type boundary whose accommodation coefficient varies in the circumferential direction. Couette flow, the thermal creep flow, and the heat transfer problem were studied. The linearized ES-BGK model of the Boltzmann equation was numerically analyzed using the finite difference method. The time-independent behavior of the gas was studied over a wide range of the Knudsen number, radii ratio, and the parameter characterizing the distribution of the accommodation coefficient. Due to the effect of nonuniform surface properties, a local heat transfer between the gas and the cylindrical surface occurs in Couette flow; a local tangential stress arises in the heat transfer problem. However, the total heat transfer between the two cylinders in Couette flow and the total torque acting on the inner cylinder in the heat transfer problem vanish irrespective of the flow parameters. Two nondegenerate reciprocity relations arise due to the effect of nonuniform surface properties. The reciprocity relations among the above-mentioned three flows were numerically confirmed over a wide range of the flow parameters. The force acting on the inner cylinder, which also arises due to the effect of nonuniform surface properties in Couette flow and the heat transfer problems, was studied in detail.

Acknowledgements

This work was supported by JSPS KAKENHI Grant (18K03949).

Conflicts of Interest

The author declares no conflicts of interest regarding the publication of this paper.

Cite this paper

Doi, T. (2019) Flows of a Rarefied Gas between Coaxial Circular Cylinders with Nonuniform Surface Properties. Open Journal of Fluid Dynamics, 9, 22-48. https://doi.org/10.4236/ojfd.2019.91002

References

  1. 1. Kogan, M. (1969) Rarefied Gas Dynamics. Plenum Press, New York. https://doi.org/10.1007/978-1-4899-6381-9

  2. 2. Cercignani, C. (1988) The Boltzmann Equation and Its Applications. Springer, New York. https://doi.org/10.1007/978-1-4612-1039-9

  3. 3. Sharipov, F. and Seleznev, V. (1998) Data on Internal Rarefied Gas Flows. Journal of Physical and Chemical Reference Data, 27, 657-706. https://doi.org/10.1063/1.556019

  4. 4. Cercignani, C. (2000) Rarefied Gas Dynamics. Cambridge University Press, Cambridge.

  5. 5. Sone, Y. (2002) Kinetic Theory and Fluid Dynamics. Birkhäuser, New York. https://doi.org/10.1007/978-1-4612-0061-1

  6. 6. Karniadakis, G., Beskok, A. and Aluru, N. (2005) Microflows and Nanoflows: Fundamentals and Simulation. Springer, New York.

  7. 7. Sone, Y. (2007) Molecular Gas Dynamics. Birkhäuser, New York. https://doi.org/10.1007/978-0-8176-4573-1

  8. 8. Veijola, T., Kuisma, H. and Lahdenperä, J. (1998) The Influence of Gas-Surface Interaction on Gas-Film Damping in a Silicon Accelerometer. Sensors and Actuators A, 66, 83-92. https://doi.org/10.1016/S0924-4247(97)01732-9

  9. 9. Kang, S.C., Crone, R.M. and Jhon, M.S. (1999) A New Molecular Gas Lubrication Theory Suitable for Head-Disk Interface Modeling. Journal of Applied Physics, 85, 5594-5596. https://doi.org/10.1063/1.369810

  10. 10. Cercignani, C., Lampis, M. and Lorenzani, S. (2004) Plane Poiseuille Flow with Symmetric and Nonsymmetric Gas-Wall Interactions. Transport Theory and Statistical Physics, 33, 545-561. https://doi.org/10.1081/TT-200053939

  11. 11. Scherer, C.S., Prolo Filho, J.F. and Barichello, L.B. (2009) An Analytical Approach to the Unified Solution of Kinetic Equations in Rarefied Gas Dynamics. I. Flow problems. Journal of Applied Mathematics and Physics (ZAMP), 60, 70-115. https://doi.org/10.1007/s00033-008-7084-4

  12. 12. Doi, T. (2014) Plane Poiseuille Flow and Thermal Transpiration of a Highly Rarefied Gas between the Two Walls of Maxwell-Type Boundaries with Different Accommodation Coefficients: Effect of a Weak External Force. Journal of Applied Mathematics and Physics (ZAMP), 66, 1805-1820. https://doi.org/10.1007/s00033-014-0454-1

  13. 13. Doi, T. (2016) Transient Couette Flow of a Rarefied Gas between Plane Parallel Walls with Different Surface Properties. Physics of Fluids, 28, Article ID: 022006. https://doi.org/10.1063/1.4941738

  14. 14. Doi, T. (2015) Poiseuille Flow and Thermal Transpiration of a Rarefied Gas between Parallel Plates: Effect of Nonuniform Surface Properties of the Plates in the Transverse Direction. ASME Journal of Fluids Engineering, 137, Article ID: 101103. https://doi.org/10.1115/1.4030490

  15. 15. Doi, T. (2015) Poiseuille Flow and Thermal Transpiration of a Rarefied Gas between Parallel Plates II: Effect of Nonuniform Surface Properties in the Longitudinal Direction. Journal of Applied Mathematics and Physics (ZAMP), 66, 3405-3423. https://doi.org/10.1007/s00033-015-0580-4

  16. 16. Kumar, A., Datta, S. and Kalyanasundaram, D. (2016) Permeability and Effective Slip in Confined Flows Transverse to Wall Slippage Patterns. Physics of Fluids, 28, Article ID: 082002. https://doi.org/10.1063/1.4959184

  17. 17. Doi, T. (2017) Reciprocity Relations in Flows of a Rarefied Gas between Plane Parallel Walls with Nonuniform Surface Properties. Journal of Applied Mathematics and Physics, 68, 66. https://doi.org/10.1007/s00033-017-0811-y

  18. 18. Doi, T. (2018) Poiseuille, Thermal Transpiration and Couette Flows of a Rarefied Gas between Plane Parallel Walls with Nonuniform Surface Properties in the Transverse Direction and Their Reciprocity Relations. Fluid Dynamics Research, 50, Article ID: 025502. https://doi.org/10.1088/1873-7005/aa962d

  19. 19. Sharipov, F. (1994) Onsager-Casimir Reciprocity Relations for Open Gaseous Systems at Arbitrary Rarefaction I. General Theory for Single Gas. Physica A, 203, 437-456. https://doi.org/10.1016/0378-4371(94)90009-4

  20. 20. Takata, S. (2009) Symmetry of the Linearized Boltzmann Equation and Its Application. Journal of Statistical Physics, 136, 751-784. https://doi.org/10.1007/s10955-009-9793-2

  21. 21. Holway, J.H. (1966) New Statistical Models for Kinetic Theory: Methods of Construction. Physics of Fluids, 9, 1658-1673. https://doi.org/10.1063/1.1761920

  22. 22. Andries, P., Tallec, P.L., Perlat, J.-P. and Perthame, B. (2000) The Gaussian-BGK Model of Boltzmann Equation with Small Prandtl Number. European Journal of Mechanics B/Fluids, 19, 813-830. https://doi.org/10.1016/S0997-7546(00)01103-1

  23. 23. Takata, S., Hattori, M. and Hasebe, T. (2016) Slip/Jump Coefficients and Knudsen-Layer Corrections for the ES Model in the Generalized Slip-Flow Theory. AIP Conference Proceedings, 1786, Article ID: 040004. https://doi.org/10.1063/1.4967542

  24. 24. Cercignani, C. and Sernagiotto, F. (1967) Cylindrical Couette Flow of a Rarefied Gas. Physics of Fluids, 10, 1200-1024. https://doi.org/10.1063/1.1762263

  25. 25. Aoki, K., Yoshida, H., Nakanishi, T. and Garcia, A.L. (2003) Inverted Velocity Profile in the Cylindrical Couette Flow of a Rarefied Gas. Physical Review E, 68, Article ID: 016302. https://doi.org/10.1103/PhysRevE.68.016302

  26. 26. Bassanini, P., Cercignani, C. and Pagani, C.D. (1968) Influence of the Accommodation Coefficient on the Heat Transfer in a Rarefied Gas. International Journal of Heat and Mass Transfer, 11, 1359-1369. https://doi.org/10.1016/0017-9310(68)90181-6

  27. 27. Aoki, K., Sone, Y., Nishino, K. and Sugimoto, H. (1991) Numerical Analysis of Unsteady Motion of a Rarefied Gas Caused by Sudden Changes of Wall Temperature with Special Interest in the Propagation of a Discontinuity in the Velocity Distribution Function. In: Beylich, A.E., Ed., Rarefied Gas Dynamics, VCH, Weinheim, 222-231.

  28. 28. Sugimoto, H. and Sone, Y. (1992) Numerical Analysis of Steady Flows of a Gas Evaporating from Its Cylindrical Condensed Phase on the Basis of Kinetic Theory. Physics of Fluids A, 4, 419-440. https://doi.org/10.1063/1.858313

  29. 29. Wakabayashi, M., Ohwada, T. and Golse, F. (1996) Numerical Analysis of the Shear and Thermal Creep Flows of a Rarefied Gas over the Plane Wall of a Maxwell-Type Boundary on the Basis of the Linearized Boltzmann Equation for Hard-Sphere Molecules. European Journal of Mechanics B/Fluids, 15, 175-201.