Effect of Ballistic Bouncing of Gas Particles across a Microchannel on Rarefied Gas Flows

This paper proposes a novel computationally efficient method of modeling rarefied gas flow in microchannels based on the newly discovered and mathematically proven Ballistic Principle of the Property Balance in Space (BPPBS). The mechanism of influence of the effect of rarefication on the gas flow is specifically investigated. Also, a differential form of the momentum balance equation governing gas flow in the channel between two parallel plates due to the pressure gradient along the channel and its exact implicit solution in the form of an integral equation have been derived. The theory does not use the generalized concept of viscosity based on the variable mean free path (MFP) in the Knudsen layer (KL). Comparing the normalized flow rate as a function of the inverse Knudsen number according to the current theory and the experimental data shows good agreement in the range of the inverse Knudsen number from 0.01 to about 40. The correlation factor is found to be about 0.995. The results show that our approach based on the BPPBS offers substantial and practical advantages in modeling and simulation of rarefied gases. The validity of the widely disseminated claim of the geometry-dependent MFP in the KL was analyzed.

and applying the second-order velocity slip boundary conditions (Maxwell-type assumption) in the following form [2] where η is the fluid dynamic viscosity, u is the gas slip velocity near the wall, u w is the tangential velocity of the wall, and f λ is the mean free path, C 1 and C 2 are constants. It is well-known from the gas kinetic theory the following relationship: where n is gas-particle density, m is gas-particle mass, v T is the magnitude of the thermal velocity, and η is the fluid dynamic viscosity in one-dimensional gas space. However, according to [3], the relationship above is "valid for flows in the region far away from the wall". Also, "within the KL, the flight paths of gas molecules are so significantly influenced by the wall that this relation becomes invalid [3]". Knudsen layer (KL), known from the literature as adjacent to a gas-solid interface region of several mean free paths (MFP) thick, and a variable MFP linked to defining a geometry-dependent "viscosity" in the adjacent to the gas-solid interface region are widely used abstractions that helped to explain gas microflow behavior. We continue the introduction by briefly describing some of wall, which is integrated into the existing models of a viscous gas flow and operated in a CFD arrangement based on the classical Navier-Stokes equations. As in the elementary kinetic theory, the viscosity is directly proportional to the molecular mean free path (MFP). The effective gas viscosity is defined similarly in the near-wall region [12] and [13].
In this approach, a wall is included in the system so that some molecules will hit the wall, and the wall will terminate their flight paths. Consequently, the MFP of all the molecules in the system may be smaller than the MFP in the unconfined gas space because of the boundary limiting effect [12] [13]. More specifically, in the most recent study, Abramov [12] introduced a new expression for the MFP of gases in the KL. First, he used a known from the kinetic theory for an ideal gas a probability distribution function ( ) Finally, Abramov [12] integrated Equation (5) over the azimuthal and polar angles of the direction of the flight of a suitable spherical coordinate system to obtain the dependence of the MFP as a function of the distance d from a wall, which has resulted in the following: where 1 is a standard notation for the exponential integral Hadjiconstantinou [14] noted that "for 0.1 Kn ≤ , kinetic effects are limited to the vicinity of the walls, one can obtain solutions of the linearized Boltzmann equation by superposing a Navier-Stokes solution with a kinetic boundary (Knudsen) layer correction, where the latter is significant only up to distances of the order of one mean free path from the wall". However, Guo et al. [3] extended the applicability of the Navier-Stokes equations to much higher Knudsen numbers (up to about 80, see Figure 4 in [3]). They claimed that because the data of their extended Navier-Stokes constitution (ENSC) model "agree well with the experimental data in the entire range of δ, and the Knudsen minimum in the flux is captured successfully", their model of the effective viscosity is valid. Still, it is difficult to imagine the applicability of the Navier-Stokes equations, which are designed to describe fluid flow in the infinite small volume, to the rarefied gas media having a gas-particle density low as one gas particle per 10 -80 distances between the walls. It is also disputed if such behavior of the effective MFP near the wall is physical or not [15]".
The authors of this study admitted "erroneous definition and calculation practices." They also highlighted that "variations in MFP as a function of the surface confinement, which disagrees with the kinetic theory and leads to wrong physical interpretations of nanoscale gas flows. This controversy occurs due to erroneous definition and calculation practices, such as consideration of gas wall collisions, using local bins smaller than a MFP, and utilizing time frames shorter than a mean collision time in the MFP calculations". They also stated that the collisions between the moving gas molecules and the solid wall should be ignored when evaluating the individual free paths. Finally, they highlighted that "[l]ocal MFP variations are physically admissible only if there are local density and temperature variations in the system".
Another approach consists in revisiting "the problem of micro-channel compressible gas flows and show[ing] that the axial diffusion of mass engendered by the density (pressure) gradient becomes increasingly significant with increased Knudsen number compared to the pressure-driven convection [16]." Dongari et al. [16] used "a recently proposed modification of the Navier-Stokes equations that include the diffusion of mass caused by the density and temperature gradients [17]". Dongari  Returning to the critical notes in [15], we also believe that the MFP is a function of the free flight distance that leads to a wrong physical interpretation of the MFP. In a "small-scale" view, a particle directed toward a wall is not aware of its possible collision with the wall during its free flight. Therefore, it may interact by collision with nearby particles with the same frequency as if it travels in the unbounded gas space unless a near-wall gas density is higher than the bulk gas density. Moreover, while recent literature presents a lot of theoretical and computational research dealing with the momentum transfer across the channel (the momentum associated with the gas flow along the channel), it is silent about the balance across the channel of both the mass and momentum associated with the thermal velocity component. Violation of these balances will lead to disagreement with the kinetic theory of gases and eventually produce an erroneous physical interpretation of microscale gas flows and, therefore, need to be verified.
Our recent publications [4] and [17] support some of the arguments in [15]. In these studies, we discuss a novel analytical molecular dynamics technique based on a newly discovered Ballistic Principle of the Property Balance in the Space (BPPBS) occupied by the gas. According to the BPPBS, in each nonmoving point of the gas space at a given time, the net rate of property (mass, momentum, or energy) influx per unit volume, formed by the converging ballistic particles (each traveling along a ballistic trajectory with certain survival probability) from the model gas system is equated to the temporal rate of property change per unit volume and the net rate of property efflux per unit volume, formed by the diverging ballistic particles. Moreover, each of the balances (mass, momentum, or energy) shall be maintained independently for any point of the gas space at any given time. In paper [4], we have validated the BPPBS by seven validation tests. The most important of them was demonstrating that, in the collision-dominated flow regime, the differential equations, which we converted from the derived integro-differential mass and momentum balance equations, are identical to the corresponding Navier-Stokes equations. This finding supports the assumption that the formulated integro-differential forms of the balance are exact implicit solutions for corresponding Navier-Stokes equations. Journal of Applied Mathematics and Physics One more of the validation tests included analysis and verification of the mass balance and momentum balance in a steady-state incompressible gas system confined by gas-solid interfaces with a purely diffuse scattering of particles at the uniform temperature. In such a system, all variables describing flow depend on the position in gas space and are not dependent on time so that the particle density n and associated with it the MFP, the mass m and the magnitude of the thermal velocity of the particle/molecule, v T , are constant. We have verified that, in each point of space between confining parallel plates, the mass and momentum balance are conserved. It implies that any deviation from the constant particle density or the MFP will violate the BPPBS and violate the physical principles.
Additionally, as we referenced above, the flow of a dilute gas near a solid surface is usually explained in various theories by non-continuum effects assumed to exist in the Knudsen layer. In contrast, our approach does not use the concept of the Knudsen layer at all. Particularly, using basic principles of Newton's mechanics and the kinetic theory of gases, we have derived an analytical representation of the gas flow velocity profile in the microchannel formed by two being at rest infinite parallel plates spaced at distance H (see Equation (116) (8) where n is the particle density, m is the mass, v T is the magnitude of the thermal velocity of the particle/molecule, σ is the momentum accommodation coefficient, which is the probability, for an incident particle, to accommodate momentum from the gas-solid interface and to scatter back in the model gas as a   Figure 7 in [24]). Again, the Knudsen minimum's appearance in the normalized volume flow rate results from neither usage of a geometry-dependent "viscosity" being generalized from variable local MFP near the gas-solid interface, nor the semi-empirical Maxwell-type assumptions but the result of the application of the BPPBS.
This paper further promotes the newly discovered BPPBS and its application to describe steady-state fluid flows in a microchannel. It specifically investigates the mechanism of influence of the Knudsen number on the rarefied gas flow. It offers a new way of computing a steady-state rarefied gas flow in a microchannel, which does not use the generalized concept of viscosity based on the variable MFP in the KL.
In Section 2 of this paper, we further analyze the validity of the currently widely disseminated claim of the MFP variation normal to the nano-channel surfaces as a function of the confinement level. We applied validation test 4.1, "Determining the total rate of collisions per unit area on a surface being in contact with the gas" of [4], for evaluation of how the MFP variation normal to the nano-channel surface claimed in [12] affects the rate of collisions per unit area.
It also provides a quantitative measure of how these results differentiate from corresponding the collisions' rate per unit area of an ideal gas according to the kinetic theory of gases.
Section 3 describes the properties and features of the model gas, the Ballistic Model's physical principles, and the general physical principles of constructing the property balance for one-dimensional steady-state gas flow.
In Section 4, we provide integro-differential forms of mass balance and momentum balance equations formulated for the gas space bounded by two infinite parallel plates.
In Section 5, we demonstrate the effect of ballistic bouncing of gas particles between confining plates on rarefied gas flows.
Section 6 shows a way to obtaining a differential form u x -momentum balance equation governing gas flow in the channel between two parallel plates due to the pressure gradient along the channel.
Finally, in Section 7, we present the discussion and conclusions and highlight the advantages of the proposed approach for modeling rarefied-gas flows or gas flows in MEMS.

Evaluation of the Claim about Spatial Variation of MFP near a Wall
In our recent publications [4] and [18], we used several validation tests to evaluate the feasibility of the BM and the BPPBS. Particularly, we determined two gas properties: the rate of collisions per unit area and the pressure exerted on a solid wall contacting with the incompressible gas. We have found that both our Journal of Applied Mathematics and Physics results and the corresponding results from the kinetic theory of gases are correlated and analytically identical. These findings support the BM and the BPPBS.
Here we use the same approach in evaluating the claim about spatial variation of MFP near a wall claimed in [12]. Figure 1 is a perspective view of the geometry and coordinate system for determining the total rate of collisions per unit area on a planar flat surface being in contact with the gas. In a system with no gravitational force, the ballistic particles have straight-line trajectories. Here our consideration is limited to an incompressible model gas at the uniform temperature in a steady-state condition.

Determining the Total Rate of Collisions Per Unit Area on a Surface Being in Contact with the Gas in a Three-Dimensional Configuration
However, following [12], we assume the spatial variation of MFP near a wall expressed analytically by Equation (5). We plan to verify whether the notion of the variable MFP near a wall is physically justified.
In a semi-sphere filled with the incompressible gas over the being at rest surface A s having a directional vector  n , each particle having an instantaneous randomly directed vector of the thermal velocity of magnitude, v T may have the instant vector-velocity component directing a particle toward the surface A s . In Adopting Equation (34) of [4] to the conditions above and assigning 1 in = Ψ , the total rate of collisions per unit area on the surface A s , is given as: where Figure 1. The perspective view of the geometry and coordinate system for determining the total rate of collisions per unit area and the pressure exerted on the surface from the gas volume. and Z V is the rate of collisions per unit volume. In [4], we derived Z V by substitution of Equation (25) in Equation (77) of [4], so that Z V is expressed as where n is particle density, v T is the magnitude of the thermal velocity, and P c is the number of particles placed within a collision tube of a unit length in an unbounded gas space.
Recognizing the similarity of the geometry of computation in [4] and [12], we point out the identity of the following parameters In Equation (12), the left-hand terms are from [12], and the right-hand terms are from [4] and Figure 1. Since the number of particles placed within a collision tube of a unit length depends on the particle's trajectory [12], which, as we mentioned above, seems to be a wrong physical interpretation, we define ( ) , c P y θ as follows The number of particles placed within a collision tube of a unit length along trajectory 103 is a function of both the distance from the wall y and the angle of the incidence θ. In Equations (9) and (11), we will use ( ) Finally, using the geometry illustrated in Figure 1, considering that ( ) (15) and and substituting Equations (11), (13), (15) and (16) in Equation (9), in which , we obtain: where β is the model parameter. For the BM (compare with Equation (81) of Journal of Applied Mathematics and Physics [4]), 0 β = . One can easily perform integration of Equation (17), in which 0 β = , which is calculated as The result of integration in Equation (18) is identical to the rate of collisions per unit area of an ideal gas according to the kinetic theory of gases, thus supporting the Ballistic Model's validity. However, modification of the scheme of calculation by incorporation of the corresponding expression for the MFP derived in [12] has led to the dependence (see Equation (17) in which 1 which is not in agreement with the kinetic theory of gases.

Determining the Total Rate of Collisions Per Unit Area on a Surface Being in Contact with the Gas in a One-Dimensional Configuration
For clarity and simplicity of quantitative evaluation, we provide below the analysis of a one-dimensional gas system. Figure  Upon modification of Equation (9) to applying to the one-dimensional configuration, the total rate of collisions per unit area on plate 201 (schematically shown by arrow 203) is given as: Figure 2. The planar view of the geometry of a one-dimensional gas system for determining the total rate of collisions per unit area on the planar flat surface from the gas volume.
where Z V is the rate of collisions per unit volume and 0 i Q is the probability of free path traveling from the control volume CV ′ at y′ to the plate 201. In [4], we defined Z V for one-dimensional incompressible gas as follows where n is gas-particle density, v T is the magnitude of the thermal velocity, and P c is the number of particles placed within a collision tube of a unit length in an unbounded gas space. The survival probability that a particle will have traveled along the ballistic trajectory in incompressible model gas at the uniform temperature from ′ y to 0 y = is calculated as Substituting Equations (20) and (21) in Equation (19) and executing integration, we obtain The result of integration in Equation (22) is identical to the derivation of the rate of collisions per unit area of an ideal gas according to the kinetic theory of gases, thus supporting the Ballistic Model's validity.
We use the same approach to evaluate the rate of collisions per unit area of the gas, in which spatial variation of MFP near a wall [12] exists. Modifying Equation (13) so that it will be applicable for one-dimensional gas configuration, we obtain In Equations (20) and (21), we will use ( ) c P y′ defined in the equation above as a replacement for P c . Substituting Equations (20) and (21) Limiting integration within the KL of three MFP, the impact from each is the most significant and in which 1 c P y′  , the equation above is approximated as follows One can see that the integration above does not produce a meaningful value for the rate of collisions per unit area of the incompressible gas. We expect similar negative results of validation for the incompressible gas if we try other forms of MFP distribution functions, such as the power-law distribution function investigated by Dongari et al. [6] or Stops' exponential [13]. We support the claim Journal of Applied Mathematics and Physics in [15] that "[c]urrent claims of MFP variation normal to the nano-channel surfaces as a function of the confinement level are false".

Physical Principles of the Ballistic Model in One-Dimensional Gas Space
This section describes the BM's physical principles in one-dimensional incompressible steady-state gas flow at the uniform temperature, a lack of external field of the force, and low flow velocity (with the Mach number less than 0.3). These principles will be further applied to describe the steady-state gas flow in a micro-channel.

Properties and Features of the Model Gas
Here we assign the following properties of the model gas, which are adapted to the steady-state gas flow from the originally proposed in [25]: 1) The model gas enables a distant transport of one or more properties, including one or more of mass, momentum, and energy by particles being in a constant state of mostly random motion and interaction by collisions.
2) Each gas particle is assigned to travel by obeying a ballistic trajectory governed by a law of motion in free space. It overcomes a distance between any of the two points of the ballistic trajectory with a certain survival probability.
3) Each gas particle is adapted to transport a combination of one or more properties, comprising mass, momentum, and energy between a point of initial collision and a point of ending collision. 4) Each point within the space occupied by the model gas is treated as a point of collisions for converging particles, each following a ballistic trajectory with the same ending point simultaneously.
5) Each point of collisions is treated as either a point source for diverging ballistic particles or a point sink for converging ballistic particles. 6) Each of the particles moving from the point source to the point sink is treated as a property carrier. The property carrier is created in the point source by obtaining one or more properties of specific values being intrinsic to the model gas surrounding the point source. It is ended in the point sink by transferring one or more properties of specific values in the point sink.
7) The velocity of a point source equals the mass flow velocity of the model gas flow in a corresponding point of the initial collisions.
One can note from the above that the model gas properties differ from the properties typically assigned to the ideal gas.

Properties and Features of the Gas-Solid Interface
Here we assign the following additional unique properties to the model gas being in contact with a gas-solid interface revealing mixed diffuse and specular scattering of particles [4] [18]: 1) Each collision on a gas-solid interface of the model gas system, which has resulted in the scattering of the diffuse particles from the gas-solid interface, is Journal of Applied Mathematics and Physics treated as an act of interaction involving a property transport from the gas-solid interface to the scattered particle.
2) Each point of the diffuse particle scattering on the gas-solid interface is treated as a heterogeneous point source for each of the scattered particles.
3) The velocity of each of the heterogeneous point sources on the gas-solid interface equals the velocity of the gas-solid interface of corresponding points of diffuse particle scattering.
4) The point source strength of the heterogeneous point sources on the gas-solid interface is directly proportional to a property accommodation coefficient σ in a corresponding point of diffuse particle.
Note: Diffuse scattering is an act of interaction involving property transfer from a gas-solid interface to a scattered particle. Specular scattering does not involve property exchange between the gas-solid interface and a scattered particle. For diffuse scattering with no adsorption effects after colliding with the surface, the particle flux per unit time diffusively scattered back from the gas-solid The mentioned above trajectories were initially introduced in our publications [4] [18], and patent [25]. In this paper, we pay attention to a different group of trajectories such as 309 and 312, associated with the ballistic particles' bouncing between the being at rest confining plates, plate 1 and plate 2. These particles are not participated in property exchange by collisions within the gas space.

Principles of Construction of the Property Balance in Micro-Channel
In the microscopic scale, the model gas flow is characterized by the group of particles of mass m, which move randomly and interact by collisions with effective collision cross-section σ c . In each of the points in space, the gas properties such as particle density n, the magnitude of thermal velocity v T , and the vector of mass flow velocity u quantify the model gas. According to the BM, each point in space occupied by the model gas may serve as both a sink accumulating property delivered by converging ballistic particles from the entire model gas system and a source dispersing into the surrounding property by diverging ballistic particles.
In the interests of simplicity, the particles are considered to have a unit mass unless otherwise stated.
Based on the BPPBS, we expect maintenance of an overall property balance in each collisions' points within the model gas system. In a steady-state gas flow, B (26) The table of the model parameters associated with defining the net rate of total property influx per unit volume and the net rate of total property efflux per unit volume is presented in Table 1.

Defining Components of a Net Rate of Total Property Influx Per Unit Volume in a Nonmoving Point of Gas Space in Micro-Channel
Referring to Figure 3 and Table 1, we formulate the following equality for a micro-channel with gas-solid interfaces revealing diffuse particle scattering:

Formulating the Net Rate of Property Influx from the Surrounding Gas in the General Nonmoving Point y in the Gas Space Confined between Two Parallel Plates
Modifying Equation (38) of [25] to applying to the one-dimensional steady-state gas system, we obtain: and where i n is a unit vector of arbitrary direction: i y y y y

Formulating the Net Rate of Property Influx from the Confining Plates Having Gas-Solid Interfaces with a Purely Diffuse Scattering Of Particles in the General Nonmoving Point y in the Gas Space Confined between Two Parallel Plates Being at Rest
Modifying Equations (40), (41), (49), and (50) of [25] to applying to a one-dimensional steady-state gas system, we obtain:  Ψ are property values acquired by diffusively scattered particles from plate 1 and plate 2, respectively.

Formulating a Net Rate of Total Property Efflux Per Unit Volume from the General Nonmoving Point of Gas Space in Micro-Channel
We have recognized that the linear dimensions of the main control volume surrounding point y need to be sufficiently small to prevent two and more consecutive collisions of the same particle with other particles within the main control volume CV [4]. In Figure 3, the arrows 313 and 314 ballistic indicate the trajectory of gas particles escaping the CV positioned in point y.
Modifying Equation (73) of [25] to applying to the one-dimensional steady-state gas system, we obtain: out y y n Q y y y y y y for incompressible gas., and where + n is a unit vector of arbitrary direction:

Analytical Representation of a General Integral Form of the Property Balance Equation in the Micro-Channel
The integro-differential form of property balance equation is formulated by Equation (41) given below, which is obtained by substitution of Equations (28),

Analytical Representation of a General Integro-Differential Form of Mass Balance Equation in the Micro-Channel
To formulate a general integro-differential form of mass balance equation in each nonmoving point of space occupied by the gas, we modify Equation (42) by executing differentiation in the right-hand of the resulting equation followed by executing limit y ′ → y , and rearranging the terms. Then, we obtain the following general integro-differential form of the mass balance equation

Analytical Representation of a General Integro-Differential Form of Momentum Balance Equation in the Micro-Channel
To formulate a general integro-differential form of momentum (along with y-axis) balance equation in each nonmoving point of space occupied by the gas, we modify Equation (42) by assigning executing differentiation in the right-hand of the resulting equation followed by executing limit y ′ → y , and rearranging the terms. Then, we obtain the following general integro-differential form of the momentum balance equation

Solving a System of Mass-Balance and Momentum Balance Equations for Determining Zb
According to the BPPBS, the mass-balance and momentum-balance shall be maintained at every point of the gas space [4]. Recognizing that in both Equation (44) and Equation (46), there are two independent variables, u and Z b , one can solve a system of the above-mentioned equations and obtain a solution that The fact that the result of the derivation of Z b from the mass-and momentum-balance equations is identical to the result of the derivation of the rate of collisions per unit area of an ideal gas according to the kinetic theory of gases supports the Ballistic Model's validity. The functional relationship expressed by the equation above may also be obtained by analyzing fluxes formed by the ballistically traveling particles. Referring to Figure 3 and Equations (45) and (47) of [25] and considering symmetricity because of the identity of properties of gas-solid interfaces of both plate 1 and plate 2, we may write the following relationship: where 1 FS J and 2 FS J are the particle fluxes per unit time on the gas-solid interfaces of plate 1 and plate 2, respectively, 1 21 N is the particle flux per unit time on the gas-solid interface of plate 1, which is associated with the incident on the gas-solid interface of ballistic particles from the gas space following ballistic trajectories 307, and 1 2 N is the particle flux per unit time on the gas-solid interface of plate 1, which is associated with an incident on the gas-solid interface of Journal of Applied Mathematics and Physics ballistic particles from originated from diffuse scattering on the gas-surface interface of plate 4. Using Equations (43) and (42) of [25] as templates, we determine ( ) ( ) respectively. Referring again to Figure 3 and Equations (46) and (48) of [25], we formulate the following relationship: ( ) Substituting Equations (48)-(50) in Equation (51) and taking performing algebraic operations, we finally obtain which is identical to Equation (47). The above implies that Z b is constant at any Knudsen number.

Finding the Effect of Rarefication on the Effective Density of Particles in the Gas Space, Which Can Carry and Exchange by Collisions ux-Momentum
Here we shall admit that flux Resolving the equation above with respect to ef V Z , we obtain: Recognizing that the rate of collisions per unit length, P c , is constant for in-Journal of Applied Mathematics and Physics Kn − , we may define the normalized effective density, k ef , as follows: where Kn is the Knudsen number defined as the ratio of the mean free path 1 f c P λ = and the representative length scale H. Results for the normalized effective density profile between two parallel plates are presented in Figure 4. We intentionally provided a graphical interpretation of Equation (54) to highlight that the normalized effective density is constant within the gas space and depends only on the Knudsen number. As the Knudsen number's value increases from 0.2 to 20, the effective normalized density decreases from 1 to 0.048.

Influence of the Effective Density of ux-Momentum Carriers on the Rarefied Gas Flow in the Microchannel
Considering the effect of the confining plates on the density of -momentum carriers in the channel's gas space, we analyze the velocity profile generated in the model gas due to the pressure gradient along the channel. We use a similar approach that has been initially introduced in our publications [4] and [18]. x y =∆ ∆ is at y within the model gas flow confined by two being at rest parallel plates at 0 y = and y H = . Plate 1 and Plate 2 confine across the y-axis a model gas flow along the x-axis and surfaces 301 and 302 positioned at distance ∆x bound a portion of the model gas system along the x-axis. Here we assume purely diffuse scattering of the particles from the walls 1 σ = , so that the probability, for an incident particle, to accommodate momentum from the gas-solid interface and to scatter back in the model gas as a diffuse particle equals to one. Also, because of the zero-velocity of the confining plates, the particles scattered from the walls do not carry u x -momentum, so that the fluxes associated with trajectories 308 and 311 are zeroed. Here we also consider that the magnitude of the thermal velocity of the particles is much higher than the magnitude of the mass flow velocity vector along the x-direction, To formulate a general integro-differential form of u x -momentum balance equation in each nonmoving point of space occupied by the gas, we modify Equation (42) by assigning executing differentiation in the right-hand of the resulting equation followed by executing limit y ′ → y , and rearranging the terms. Then, we obtain the following general integro-differential form of the momentum balance equation: One can note that the equation above, in which 1 eff k = and ef V V Z Z = , is identical to Equation (112) of our paper [4], in which 1 σ = . We will explain usage of k ef upon the following analysis. We have recognized that some of the particles initiated from collisions in the gas space and carrying u x momentum may travel by free path to a wall, then after accommodation with the gas-solid interface of a wall diffusively scatter back in the gas space. Here we shall recite that "the cause of the thermalization and "diffuse" escape of particles from the surface lies not in the trapping of molecules by the surface, as assumed by Maxwell himself and repeated in practically all the literature published on this question, but in the relaxation, in interaction with the phonon subsystem" [26]. Since the phonon-relaxation time is about 10 −11 -10 −10 sec, the thermalized particle scatters back from a gas-solid interface almost immediately. All the particles scattered from the walls do not carry u x -momentum. Some of these scattered particles may even travel to the opposite wall by free path. Such group of the particles will not participate in u x momentum transport anymore. This effect was not taken into consideration in our publications [4] and [18].
Here we also understand that the flow through the channel is motivated by the surface force exerted by the surroundings on the CV through the pressure force, and the pressure is directly proportional to gas particles density. We suggest N. Kislov Journal of Applied Mathematics and Physics that, in the gas of particle density n and the effective density n eff , the effective pressure force applied to the differential CV at y through the active gas particles may be expressed as Here we need to admit that while evaluating the effective pressure force applied to the differential CV, we used the effective density n ef of active particles, which is represented in the right-hand terms of Equation (59) (30) in Equation (59) is a non-dimensional velocity of the gas flow in the channel. The equation above's explicit solution is obtained numerically by sequential approximation described previously in [18].
Further integration of ( ) x y U along y-direction and normalization by H will result in finding a normalized volume flow rate xA U : (65) Then we can obtain the normalized velocity profile as: In Figure 5, we show the velocity profiles ( ) xN y U in the micro channel between two infinite parallel plates at different Knudsen numbers. We observe in Figure 5 that the velocity is nonzero at the channel walls, even at 0.001 Kn = .
It is also worth noting that the normalized volume flow rate xA U as a function of the inverse Knudsen number, according to the current theory, has the so-called "Knudsen minimum" at 0.9 δ ≅ . It is in good agreement with Knudsen's observation of the minimum at about ~1 Kn [29].

Deriving a Differential form ux Momentum Balance Equation Governing Gas Flow in the Channel between Two Parallel Plates Due to the Pressure Gradient along the Channel and Its Analysis
We have simulated gas flow in a plane channel using the BPPBS and obtained the gas flow velocity profile in the implicit form shown by Equation (64). It would also be interesting to obtain and analyze a differential form of the equation governing such gas flow in the channel. Below we show a way to obtaining a differential form u x momentum balance equation governing gas flow in the channel between two parallel plates due to the pressure gradient along the channel, like the presented in Step 2 of Section 4.4.3 of our paper [4].
Applying to Equation (64) The equation above is different from the classical Navier-Stokes equation applicable to incompressible Newtonian flow. Here we will not further proceed to obtain an explicit analytical solution of the differential equation above, as we did in Section 4.4.3 of our paper [4]. However, we point out that Equation (61) Remarkably, the normalized volume flow rate reaches a constant value for n → ∞ , which agrees with Knudsen's finding [29].

Discussion and Conclusions
The recently discovered and mathematically proven BPPBS [4] is used to determine the rarefied gas flow rate through the microchannel spaced between two infinite parallel plates in the current paper. Our theory does not use the generalized concept of viscosity based on the variable MFP in the KL. We formulated the Ballistic Model's physical principles for one-dimensional incompressible steady-state gas flow at the uniform temperature, a lack of external field of the force, and low flow velocity (with the Mach number less than 0.3). We have formulated the mass and momentum balance equations (momentum associated with the thermal component of velocity). Then we have solved them, revealing that the rate of collisions on the plates' gas-solid interface is independent of the rarefication (independent from Knudsen number). Since there is a chance for a particle diffusively scattered from one confining plate to reaching by the free path the opposite plate, we were able to derive an analytical function revealing the dependence of the effective density of particles currying u x -momentum from the Knudsen number. It has been found that the effective normalized density is Journal of Applied Mathematics and Physics constant within the gas space and depends only on the Knudsen number. Considering the effect of the confining plates on the density of u x -momentum carriers in the channel's gas space, we have formulated an integral form expressing implicitly the velocity profile generated in the model gas due to the pressure gradient along the channel. The implicit integral form of the velocity profile has been solved numerically by the sequential approximation. Comparing the normalized flow rate as a function of the inverse Knudsen number according to the current theory and the experimental data shows good agreement in the range of the inverse Knudsen number from 0.01 to about 40. The correlation factor is found to be about 0.995. Also, a differential form of the u x -momentum balance equation governing gas flow in the channel between two parallel plates due to the pressure gradient along the channel has been derived. It was noted that the differential form of the u x -momentum balance equation derived from the implicit integral form of the velocity profile is different from the classical Navier-Stokes equation applicable to incompressible Newtonian flow. However, the governing equation derived in the present theory is reduced, for 0 Kn → , to the Navier-Stokes equation for well-known general Couette flow. In turn, for Kn → ∞ , the governing equation is reduced to a constant value, thus representing the free-molecular regime. It implies that both the integral and the differential forms, according to our theory, are applicable to describe gas flow in the microchannel in the range of Knudsen numbers from zero to infinity.
We also analyzed the validity of the currently widely disseminated claim of the MFP variation in the KL as a function of the distance from the gas-solid interface. For this purpose, we compared the rate of collisions per unit area on a gas-solid for the incompressible gas at the uniform temperature in a semi-infinite space for our and other models based on the variable MFP in the KL region. First, we found that the collision rate per unit area obtained by applying the BM and the corresponding characteristic known from the kinetic theory of gases are correlated and analytically identical. These findings support our theory based on the BM and the BPPBS. Second, we incorporated in the equation expressing the total rate of collisions per unit area on the solid surface, which we have introduced earlier in [4], the function representing the variable MFP [12]. However, in the latter case, upon executing identical algebra operations, we could not produce a meaningful value for the collisions' rate per unit area of the incompressible gas.
Based on the results above, we may conclude the following: fictitious 1) "Knudsen layer" seems to be a useful abstraction to employ specific boundary conditions providing reasonable approximate solution beyond the Knudsen layer (by using Navier-Stokes equations or a lattice Boltzmann (LB) method to model a rarefied gas flow). However, gas particles' microscopic and physical behavior in the KL is not well understood, leading to a pure prediction of gas flow characteristics.
2) The assumption that the MFP varies as a function of the surface confinement (Knudsen number), widely disseminated in the literature, is erroneous. Journal of Applied Mathematics and Physics Such an assumption disagrees with the kinetic theory of gases. As we found, the interaction by collisions of the gas with variation in MFP does not produce a meaningful value for the rate of collisions per unit area. As a result, such an assumption leads to wrong physical interpretations of nanoscale gas flows.
3) We offer a new way of computing a steady-state rarefied gas flow in a microchannel based on recently discovered and mathematically proven by us the Ballistic Principle of the Property Balance in Space. We specifically investigated the mechanism of influence of the Knudsen number on the rarefied gas flow.

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