The Knudsen Paradox in Micro-Channel Poiseuille Flows with a Symmetric Particle

Featured Application: The present work addresses a perceived knowledge gap in rareﬁed ﬂows, along with a relevant challenge while modeling gas–solid ﬂows in conﬁned geometries at the nano-scale, where simultaneous handling of local and non-local transport mechanisms over the particle surfaces must be realized. These phenomena are also of interest to several micro-ﬂuidic applications including (but not limited to) lab-on-a-chip devices, micro-total analytic systems ( µ TAS) and point-of-care diagnostics (POC). Abstract: The Knudsen paradox—the non-monotonous variation of mass-ﬂow rate with the Knudsen number—is a unique and well-established signature of micro-channel rareﬁed ﬂows. A particle which is not of insigniﬁcant size in relation to the duct geometry can signiﬁcantly alter the ﬂow behavior when introduced in such a system. In this work, we investigate the effects of a stationary particle on a micro-channel Poiseuille ﬂow, from continuum to free-molecular conditions, using the direct simulation Monte-Carlo (DSMC) method. We establish a hydrodynamic basis for such an investigation by evaluating the ﬂow around the particle and study the blockage effect on the Knudsen paradox. Our results show that with the presence of a particle this paradoxical behavior is altered. The effect is more signiﬁcant as the particle becomes large and results from a shift towards relatively more ballistic molecular motion at shorter geometrical distances. The need to account for combinations of local and non-local transport effects in modeling reactive gas–solid ﬂows in conﬁned geometries at the nano-scale and in nanofabrication of model pore systems is discussed in relation to these results.


Introduction
A sub-micron sized particle when transported in a micro-channel can significantly alter the rarefied behavior of the carrier gas. These effects emerge partly because the dispersed particle serves as an extra momentum source or sink in the system and partly due to the blockage induced by the particle, which further affects the effective hydrodynamic length of the bounding duct. These phenomena are of particular significance from a theoretical view-point due to the complex molecular interactions that prevail at the aforementioned spatial scales, necessitating novel approaches to understand the same. Moreover, these challenging systems are also of interest to several micro-fluidic applications including (but not limited to) lab-on-a-chip devices, micro-total analytic systems (µTAS) and point-of-care diagnostics (POC) [1,2]. The Knudsen paradox [3][4][5] is a unique signature of micro-channel rarefied flows. Usually, such flows are characterized by the Knudsen number, Kn = λ/a, where λ is the mean-free-path of the gas and a is a suitable system length scale (for instance, the channel height). Consequently, these flows can be classified into four regimes: continuum flows (Kn < 0.001), slip flows (0.001 ≤ Kn < 0.1), transition flows (0.1 ≤ Kn < 10) and free-molecular flows (Kn ≥ 10) [6,7]. The Knudsen paradox is characterized by the variation in the mass-flow rate (Q) of a molecular gas in micro-channels across these aforementioned regimes, where the ratio between the length and hydraulic diameter is large. Note that the flow in these micro-channels is driven by the inlet-outlet pressure gradient. Correspondingly, in such systems, Q decreases with increasing Kn, reaches a minimum at Kn ≈ O(1) (in the transition regime) and further increases with an increase in Kn (towards the free-molecular regime). For this reason, the Knudsen paradox is also referred to in literature as the Knudsen minimum. The paradox was first studied experimentally and theoretically by Knudsen [3] in experiments of Poiseuille flow driven by the identical pressure drop in channels with varying widths. The observation of this minimum was attributed, by Pollard and Present [4], to the imbalance between two counteracting molecular effects -the obstruction of long diffusion paths due to molecular collisions, leading to a net decrease in Q on the one hand and the development of drift transport leading to an increase in Q on the other hand. Thus, it follows that Q in a long tube must initially decrease with increasing pressure (when L > λ a, where L is the channel length) and pass through a minimum value due to the obstruction to molecular diffusion. As the pressure crosses this threshold (λ < a), drift transport increases (towards the Poiseuille form) leading to an increase in Q.
The Knudsen paradox has also been confirmed through numerical assessments, first by Cercignani and co-workers [5,8,9] and later by others (e.g., [10][11][12]). These authors solved the steady Boltzmann equation with a BGK (Bhatnagar-Gross-Krook) [13] collision term. For a detailed review of these approaches, the reader is referred to the book by Cercignani [14]. More recently, these numerical methods have been further extended and validated to handle a wide range of rarefied flows in a variety of geometries, as described in the reviews by Sharipov [15], Beskok [16], Karniadakis [17] and Colin [18]. In general, the reported methods are at the following levels of abstraction: (i) micro-scale resolved techniques that directly model molecular scale effects (such as inter-molecular and molecule-wall interactions); (ii) meso-scale resolved techniques (using the Boltzmann equation) that solve the molecular distribution function; and (iii) macro-scale (continuum) resolved techniques that incorporate velocity slip and temperature jump boundary conditions into the Navier-Stokes equations to model molecular phenomena. The continuum based methods are valid up to the transition regime [19], while the meso-and micro-scale methods can in principle be applied across the entire Kn range [6]. Thus, these two alternatives are preferred while evaluating the Knudsen paradox. Consequently, the lattice Boltzmann method (LBM) [20,21] and the direct simulation Monte-Carlo (DSMC) method [22,23] have gained prominence in rarefied micro-channel flow studies.
In general, research on single-phase rarefied Poiseuille flows are extensive (cf. [16][17][18]). However, related work on their multiphase counterparts (where a solid particle and a rarefied gas co-exist in the micro-channel) are relatively scarce. Currently, these reported multiphase studies are mostly limited to assessing stationary flow past a particle (e.g., [24][25][26][27]), indicative of a knowledge gap. In fact, investigations on the effect of an additional particulate phase on the Knudsen paradox noted in micro-channel Poiseuille flows (in the transition regime) have not been previously reported. An additional motive to study such rarefied particulate Poiseuille flows stems from the recent advancements in nanofabrication techniques. Novel nanofabricated reactor systems are increasingly used to study nanoscale reactive flows [28]. Such systems are often designed to mimic naturally occurring reactor geometries, such as individual pores inside a porous medium. Particles can either be dispersed (as in a gas-solid flow) or stationary (as in typical heterogeneous catalysis applications). A requirement for the existence of a Knudsen minimum is that the duct exhibits free paths for molecular motion that are much longer than the mean radius, a configuration which is typically not readily encountered in natural porous media [4]. However, due to restrictions imposed by nanofabrication as well as increasing demand for advanced experimental equipment, artificial model pores may well fulfill these requirements and thus exhibit Knudsen minima effects.
Furthermore, the spatial variability in the heterogeneous topology of natural porous media has been shown to reduce the relevance of rarefaction effects on the total flow rate in the slip flow regime [29], which could imply that effects observed in systems with regular pore patterns may differ from the systems they are meant to represent. As model pore manufacturing methods are swiftly advancing, possibilities to experimentally investigate reactive rarefied gas-solid flows in model pores will soon become available. It is then of great interest to understand the differences in flow and transport characteristics of artificially designed systems to their intended natural counterparts. Such studies are particularly relevant for understanding pressure-driven transport of nanoparticles through micro-channels as encountered during reactive particulate matter flows in catalytic automobile exhaust after-treatment devices [27]. The physical scales studied here are also relevant for nanoparticle and dust transport in ambient to low-pressure environments inside natural porous media (such as ceramic substrates) or man-made electronics (such as printed circuit boards) [30][31][32]. This paper aims to address these open questions in rarefied Poiseuille flows with the simultaneous presence of a particulate phase. Consequently, we investigate the Knudsen paradox in a regular micro-channel Poiseuille flow with a stationary obstacle placed at its center and discuss the implications of our observations for reactive particulate flows.
The obstacle is intended to represent a particle which is not small in relation to the duct geometry. For simplicity, the particle is assumed to be square shaped, occupying the entire depth of the micro-channel and immobile. We use the DSMC method developed by Bird [22,23] for the respective assessments. This is an efficient particle-based (micro-scale) method that simulates interactions between molecular parcels, which represents groups of molecules having the same properties such as velocity and temperature. The relevant macroscopic quantities of the flow are extracted by statistically averaging over these microscopic quantities. Consequently, we employ the thoroughly validated open source DSMC solver dsmcFoam+ (implemented within the OpenFOAM software framework [33]), developed by White and co-workers [34], for the simulations presented in this paper. We evaluate the variation of a normalized form of Q (referred to as Q norm ) across the entire Kn range in micro-channel Poiseuille flows with a particle. We investigate the effect of the particle size, in terms of a blockage ratio given as H/h, where H and h are the channel and particle heights, respectively, on the obtained Knudsen minimum. Additionally, the velocity along with the Mach number developed across the channel, are estimated to establish a hydrodynamic basis for the noted variation in Q norm .

Problem Description
In this paper, we evaluate the pressure driven flow of rarefied nitrogen gas in a microchannel. This is a Cartesian 2D planar setup with the flow along the x-direction. A stationary square-shaped particle is placed at the center of the micro-channel. We consider three different configurations, given by particle blockage ratios (or H/h where H and h are the micro-channel and particle heights, respectively) 8, 4 and 2, which represent progressively increasing degree of channel obstruction (see Figure 1). The channel length (L) is five times its height (H), which is is 0.4 µm, to ensure that a fully developed flow with minimal entrance effects can be realized. The pressure drop (∆P) across the channel is varied so as to reasonably span over the entire Kn range. Such a fundamental setup is chosen since it perfectly mimics the flow configuration studied in [5,35,36] and is well suited to describe the Knudsen paradox in micro-channels.

DSMC Framework: dsmcFoam+
The direct simulation Monte-Carlo (DSMC), developed by Bird [23], is a stochastic particle-based technique for modeling real gases in which the molecular diameter is much smaller than the mean free path (λ). The primary objective of this method is to decouple molecular motion (modeled deterministically) and inter-molecular collision (modeled probabilistically) over small time intervals. The time interval over which the solution is sought is subdivided into smaller sub-intervals over which the particle motion and collisions are decoupled. Consequently, the trajectories are computed over a time interval which is lesser than the mean collision time. Further, DSMC relies on discrete parcels, which each contain a collection of molecules, to represent the total molecular number density (N) of the real gas. Thus, this method is computationally cheaper than a fully deterministic method like molecular dynamics (MD) and is ideally suited for studying rarefied gas dynamics [23].
We use the open-source solver dsmcFoam+ [34,36,37] for carrying out the assessments presented in this paper. This code-base is licensed under the GNU general public license with a publicly available software repository [38]. Further, it has been extensively validated against experimental and theoretical data across a wide variety of rarefied flows [34][35][36]. The algorithm employed in this solver can be summarized as follows. The flow domain is divided into cells (a Cartesian grid), which are in-turn populated with the DSMC parcels. Each of these parcels is assigned a position and a Maxwellian-distributed velocity based on the inlet conditions. Consequently, a Dirichlet type implicit boundary condition is used to assign the quantities associated with the molecules entering the computational domain (cf. [23,39]). The parcels are indexed into cells and arranged in an array to facilitate collisions. These parcels are then moved based on the type of interaction. Inter-molecular collisions are carried out probabilistically using the no-time counter (NTC) scheme proposed by Bird [23], while parcel-boundary interactions are simulated using a diffuse reflection model. The variable hard sphere (VHS) model [23] is used here to model the pair-wise collisions. The overall flow of molecules in the domain are driven by a downstream pressure outlet (Dirichlet type boundary condition for pressure, cf. [23,40]). These steps are repeated to increase the sample size until the statistical errors are small enough. Finally, the desired macroscopic flow properties are calculated by spatially and temporally averaging the sampled data. Note that, for the gas (di-atomic N 2 ) simulated in this study, the molecular collisions are considered elastic and they conserve both momentum and linear kinetic energy (of collision). For a more detailed account of the DSMC method and the dsmcFoam+ solver along with a detailed validation of the code, the reader is referred to [34,36,37].

Numerical Setup
The channel, considered in the DSMC simulations, has dimensions as shown in Figure  1 and is spatially discretized into 100 × 60 × 1 computational cells along the x, y and z directions, respectively. Note that this is a pseudo-3D setup where the z direction is not solved for (i.e., a 2D planar simulation). This is done so as to mimic the setups of Cercignani et al. [5,8] and other similar studies [34][35][36]. A stationary obstacle (as shown in Figure 1) is introduced in the center of the channel to obstruct the rarefied flow at varying degrees. This is given by particle blockage ratios (or H/h) 8, 4 and 2, which represent progressively increasing degree of channel obstruction. The surface and inlet temperatures are set to 300 K, while the inlet pressure is fixed based on the desired Kn (in our study, this varies from 5000 pa to 255,000 Pa). Further, the inlet to outlet pressure ratio is maintained at 3.0 for all the cases run (synonymous with the range prescribed in similar studies [34,35,37,41]). The time step size used in computing the DSMC parcel motion was determined smaller than the local mean collision time. The values of mass flow rate at the inlet and outlet are monitored until these are equivalent. The relevant macroscopic properties (e.g., temperature, density and velocity) are sampled after this steady state has been attained. The sampling is performed within each computational cell over a sufficiently large time step and over multiple realizations (to minimize statistical scattering) using relationships given by Hadjiconstantinou et al. [42]. These simulation conditions are summarized in Table 1. A typical DSMC simulation is usually constrained by the following discretization criteria: (i) the computational cell size must be smaller than the local mean free path (if possible collision partners are restricted to a parcel's current cell); (ii) the number of parcels per cell must be large enough to preserve collision statistics; and (iii) the simulation time step must be chosen so that the parcels only cross a fraction of the average cell length during each time step. Conditions (ii) and (iii) are satisfied by employing a time-step much smaller than the mean collision time, while (i) is met by utilizing a spatial discretization that has already been validated exhaustively for the current setup (cf. [34][35][36]). All the simulations were carried out serially on a single node with 20 cores using an Intel ® Xeon ® (E5-2650) processor (up to 3.0 GHz). The simulations contained around 200,000 DSMC parcels in the case at a pressure of 255,000 Pa (each parcel has a weighting factor of 10). Note that the parcel count reduces as the pressure (and consequently the mass flow of the gas) in the micro-channel reduces. Hence, to minimize the statistical scatter in the collected data, two realizations of each DSMC simulation were run for a total of 100,000 sampling intervals. The results presented in this paper are ensemble averaged over these two realizations (the scatter was noticeably reduced). Furthermore, the velocity extracted from the simulation (shown in Figure 2) is displayed both using the scattered raw data as well as a polynomial fitted using a weighted least squares estimate [43], so as to facilitate a simplified visualization of the data. Each DSMC simulation needed a physical time of 100 h.

Results and Discussion
We assess the Knudsen paradox in the numerical setup described above by varying the inlet pressure (P in ) and maintaining the inlet to outlet pressure ratio as 3. To compare our DSMC results to the results from Cercignani et al. [5,8], the mass flow rate Q (obtained from the simulations) is represented in a non-dimensional form. This dimensionless flow rate Q norm used is [34,41]: where L, H and w are the length, height and width of the micro-channel, respectively. Further, R and T are the specific gas constant and iso-thermal temperature of the gas. The variableṁ is the mass flux extracted by the dsmcFoam+ solver (using the mass flux measurement tool described in [34]) and P in and P o are the inlet and outlet pressures. Note that Q norm is independent of the local pressure gradient (depends only on its mean value), but depends mainly on the channel height represented in terms of a rarefaction parameter δ m (i.e., on an inverse Knudsen number-λ/H) [41] given by: Equations (1) and (2) can be used to describe the rarefied flow with the particle as well under the same constraints, i.e., the local pressure gradient (defined as H P dP dx ) is small (or 1). Correspondingly, we show the comparisons between dsmcFoam+ and the data of Cercignani et al. [5,8], along with the effect of particle blockage on the Knudsen paradox. Prior to this discussion, we first establish a hydrodynamic basis for the effect of the stationary particle in the channel. This is done by comparing the corresponding velocity profiles across the micro-channel setups with and without the particle in the transition regime (at Kn = 0.1). We also briefly evaluate the compressibility effects of the particle on the flow system by examining the Mach number developed in in this flow regime.

Hydrodynamics of a Stationary Particle in the Transition Regime
The presence of a stationary and inert particle in a 2D planar micro-channel Poiseuille flow of a rarefied gas in the transition regime alters the inherent response of the system. This is because of the additional momentum sink in the flow due to the particle. We investigate the hydrodynamic consequences of such a flow obstruction by comparing the velocity in the micro-channel with a benchmark case without the particle. Such an assessment has a two-fold objective: (i) deviations from the ideal micro-channel Poiseuille flow due to the presence of the particle can be identified; and (ii) the hydrodynamic basis for the observed Knudsen paradox can be deduced. We first discuss the impacts on the channel velocity profile (along the height) and follow this with a brief discussion on the compressibility of such particulate flows. Note that all these hydrodynamic assessments are done under the conditions defined in Table 1 and at an outlet pressure P o of 1 bar (Kn = 0.1). This Kn is chosen as it represents the unique hydrodynamics in the transition regime, where the Knudsen minimum has been observed in micro-channel Poiseuille flows [5,41].
The 2D planar micro-channel Poiseuille flow benchmark case without the particle is characterized by a parabolic velocity profile. Consequently, when the flow encounters a stationary particle at the center of the channel, its response is markedly altered. Figure 2 shows the variation of the co-axial velocity along the y-direction (U x (y)). As expected, the presence of the particle decelerates the overall flow, with the maximum velocity reduced due to the blockage. It follows that, at the center of the channel, the particle partitions the flow into two narrower channels each developing its own parabolic velocity profile for U x (y), with a maximum velocity lower than the benchmark case. Note that velocities at the wall boundaries are always nonzero due to a net slip-effect in the transition regime. The non-uniformities noted in the reported velocity profiles are attributed to the statistical scatter from the DSMC data. This scatter is more visible at the lowest blockage ratios as the largest particle drastically alters the local length scale of the system, thereby decreasing the probability of molecular collisions in this region, which in turn increases the need for sampling to converge the statistics. Nevertheless, the results obtained for the lowest blockage ratio are deemed reasonable enough as they follow the expected trend despite the scatter. These observations are further confirmed by the streamlines of the velocity along a x − y plane, shown in terms of the Mach number (Ma = u/c, where u is the local fluid velocity and c is the speed of sound) in Figure 3. The flow accelerates as it passes the particle, partly due to the pressure decrease but mostly due to continuity. Further, a symmetric low-velocity region is developed before and after the particle. The channeling effect created by the presence of the particle has a consequent impact on the Knudsen paradox, as will be explained later.
Furthermore, the mach number (Ma) is also used to assess compressibility effects in a flow system [44], where such effects are generally deemed negligible for Ma < 0.3. This criterion may however be violated in micro-channel Poiseuille flows. Figure 3 shows Ma in the micro-channel for all cases with and without the particle at the conditions reported in Table 1. The benchmark case shows that such Poiseuille flows are sub-sonic in nature with Ma < 0.8. The highest Ma is observed at the flow outlet as expected. The additional momentum sink, due to the particle, leads to a lower flow velocity of the rarefied gas at the outlet and as the particle size increases, there is a corresponding decrease in the developed Ma. Thus, compressibility effects can be disregarded while assessing micro-channel Poiseuille flows with a stationary particle with a non-negligible size in relation to the duct geometry (at the same driving pressure difference). This is an important observation for numerically representing such systems, as a constant density can be assumed while setting up the study.  Table 1) and the DSMC simulations at blockage ratios H/h of 8, 4 and 2, respectively.

Effect of the Particle on the Knudsen Paradox
We assess the Knudsen paradox in the 2D setup (described in Figure 1) by varying the inlet pressure (P in ) and maintaining the inlet to outlet pressure ratio as 3. As detailed above, the normalized mass flow rate in the channel Q norm (see Equation (1)) is extracted from the DSMC simulations across the entire Kn range, i.e., at different δ m (which is in turn related to Kn as shown in Equation (2)). Figure 4 shows the Knudsen minima estimated from the DSMC simulations along with inherent uncertainties in the sampled data (through error bars). In this figure, the solid lines represent cases without a particle, while the dashed lines represent the cases with one. The comparisons between the benchmark case (without a particle) and the results given by Cercignani et al. [5,8] are reasonable, given the statistical spread of the sampled data. The agreement is good at high δ m (or low Kn) values with a deviation noted at the lower ones. The result in this region is more sensitive to the scatter in the DSMC data (indicated by the wider error bars in this region) and is also known to be affected by the lateral boundary condition in pseudo-3D systems [45]. However, these obtained results clearly depict the Knudsen paradox (within the limits of certainty of the simulation setup) and are synonymous with those reported by White et al. [34] using the same solver. Moreover, the minima is noted at δ m = 0.97 (which is at around Kn ≈ 1, i.e in the transition regime), as reported by others (cf. [5,8,41]). The predicted Knudsen minima for all the cases simulated are summarized in Table 2.
When a particle is introduced in the channel, it serves as an additional momentum sink in the system and sections of high fluid velocity are created in the regions between the particle and wall (see Figure 3). The Knudsen minimum may have shifted towards higher Knudsen numbers (see Table 2), although our data does not provide enough statistical confidence to confirm this. Nevertheless, the noted trends could be due to a shift towards relatively more ballistic molecules (that do not collide with each other) at shorter geometrical distances-a trend similar to that observed when going from a channel flow to a duct flow [45] or with flows in a bent channel [46]. In other words, the presence of the particle reduces the governing local length scales of the system, with a decreased probability of molecular interactions along the diffusion paths of the gas molecules. Moreover, the momentum sink effect causes the total mass flow rate in the channel to be reduced. Such a reduction in Q norm is noticeable in Figure 4, with lower mass flows noted as the particle size increases, thereby effectively blocking the channel to a higher degree. It can also be noted that the relative changes in Q norm with increasing particle size in the transition regime diminishes, with the curve tending to flatten out. This occurs as beyond the continuum regime the frictional losses become relatively less important due to velocity slip. Moreover, as the particle becomes larger, there is an increased chance that the Knudsen paradox totally disappears, which means that Q norm monotonically decreases without any increase at higher Kn (or lower δ m ) similar to observations made in bent channels by Liu et al. [46]. Synonymous with the inferences drawn from the bent channel studies, we also attribute the observed behavior to the apparent increase in the local hydrodynamic resistances due to the particle. Thus, the global mass flow rate is now determined by the momentum loss imposed by the stationary particle and in this region the free paths for molecular motion are not longer than the mean channel height anymore. Note that on the basis of our data the blockage ratio (H/h) of 2 could represent a limiting case for the transition from a single to a multiple channel system (due to the partitioning effect of the flow from the particle). Such a channel system could possibly be characterized by a Knudsen number which is four times that of the upstream duct. Indeed, the transition is found (see Table 2) to occur at around Kn ≈ 5 which supports this hypothesis, given the scatter in the data.  ) with open squares represents the benchmark results (without a particle). The dashed lines represent the results obtained from the DSMC simulations of a micro-channel at blockage ratios H/h of 8, 4 and 2, respectively. The error bars indicate the statistical spread (±1 standard deviation) in the sampled data. Note that the error bars increase as H/h decreases (or the particle size increases) since the local Knudsen number in the restrictions created by the particle is higher. Table 2. Knudsen minimum in 2D planar micro-channel Poiseuille flows with a stationary particle in the middle. Note that the reported minima is subject to the variations in the sampled data shown in Figure 4.

Implications for Reactive Particulate Flows
The assessments presented in this paper are particularly relevant when studying local mass diffusivities in the transition regime. In such studies, the local mass diffusivity becomes influenced by the particle-to-wall distance and the determination of a local diffusivity (for example, in the region between the particle and the wall) requires a local Kn [47]. Both the continuum approximation and the assumption of thermodynamic quasiequilibrium break down at finite Kn [15,48], implying that completely different mechanisms may govern mass transfer on different sides of an asymmetrically positioned particle.
It has also been shown that the mass transfer coefficient in pressure-driven channel flows at varying degrees of rarefaction is essentially constant with Kn, whereas the heat transfer coefficient is strongly dependent on Kn, leading to a breakdown of the well-known heat and mass transfer analogies [49]. This difference in behavior originates from the non-local transfer phenomena: whereas a molecule coming from a distant region with a significantly different temperature is able to transfer a significantly different energy load (to the molecules in the near-particle region), the mass transfer capability of any given molecule is always the same, independent of the region from which the molecule originates [49]. Moreover, for nanofabricated systems with parallel channels the existence of Knudsen minima implies that the mass-flow rate through a given channel may vary non-monotonously with the degree of rarefaction in the transition regime. As the local rarefaction will depend on the local blockage and particle-wall distances for a gas-solid flow and as these may exhibit highly dynamic characteristics, oscillations in flow rate (and thus in retention time) may be expected to arise spontaneously under certain conditions. For reactive systems, where particle sizes change continuously due to chemical reactions, and which may be controlled by either mass or heat transfer or chemical kinetics, controlling the non-local transport characteristics may thus require accounting for highly complex Knudsen minima effects.
The local mass diffusivity in a reactive nano-scale system will thus depend on the local flow configuration. The present work demonstrates that, for an immobile particle, the Knudsen paradox may offer a possibility to indirectly assess or characterize these phenomena from the response of the global mass-flow rate to variations in the absolute pressure. Finally, for a suspended particle that is moving with the average gas flow velocity, the momentum sink effect may become insignificant, but the effects of the local length-scale reduction will remain in full. The latter effects are also the most relevant for nano-scale reactive flows.

Conclusions
In this paper, we present a DSMC investigation of the effect of a stationary particle on the Knudsen paradox observed in micro-channel Poiseuille flows. Such an investigation has previously not been reported in the literature. The objective with this study was to assess the effects from representing natural porous media through nanofabricated model pores in rarefied gas-solid systems and heterogeneous catalysis applications. The open source solver dsmcFoam+ was used to carry out the relevant studies. We show that the presence of a particle markedly alters the response of the rarefied flow system as it serves as an additional momentum sink in the flow. Further, as it becomes larger and subsequently blocks a larger portion of the duct (represented as the blockage ratios H/h), the mass-flow rate at the same nominal pressure difference decreases. Despite the noticeable scatter in the reported data, the general trends (in Q norm ) observed agree with conventional rarefied flow dynamics in micro-channels and parallelisms can be drawn from our study and that of Liu et al. [46] on micro-channels with bends. In both assessments, an increase in the local hydrodynamic resistances in the flow leads to noticeable deviations from the behavior derived for straight micro-channel Poiseuille flows by Cercignani et al. [5]. Consequently, with the presence of a particle, the Knudsen minimum is shifted towards higher Knudsen numbers with the curve flattening out towards the free-molecular flow regime. This effect is more pronounced as the particle becomes large and it results from a shift towards relatively more ballistic molecular motion at shorter geometrical distances. These observations have significant implications for the design of model pore systems for reactive rarefied gas-solid flows, as well as investigations of heterogeneous catalysis using stationary catalyst particles. Most notably, the local mass diffusivity, which determines the overall reaction rate at masstransfer limited conditions, will depend on the local flow configuration, which may thus be indirectly assessed from the response of the global mass-flow rate as a function of absolute pressure. Moreover, attention should be paid to the possibility for flow oscillations to manifest in nano-scale reactor systems with rarefied reactive gas-solid flow.
In conclusion, the present work addresses a perceived knowledge gap in rarefied flows along with a relevant challenge while modeling reactive gas-solid flows in confined geometries at the nano-scale, where simultaneous handling of local and non-local transport mechanisms over the particle surfaces must be realized. The observed trends can be extended towards asymmetric particles as well, since the presence of any particle in the flow increases the overall hydrodynamic resistances in it. Thus, we underline the need to account for these complex transition regime effects while assessing such rarefied gasparticle systems, particularly in the relevant micro-fluidic applications.

Data Availability Statement:
The original data is available from the authors upon request.

Acknowledgments:
The authors would like to thank the Centre for Scientific and Technical Computing at Chalmers University of Technology (C3SE), for providing the necessary computational resources.

Conflicts of Interest:
The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

Abbreviations
The following abbreviations are used in this manuscript: