## Abstract

Understanding particle detachment from surfaces is necessary to better characterize dust generation and entrainment. Previous work has studied the detachment of particles from flat surfaces. This work generalizes this to investigate the aerodynamics of a particle attached to various locations on a model hill. This work serves as a model for dust aerosolization in a tube, as powder is injected into the Venturi dustiness tester (VDT). The particle is represented as a sphere in a parallel plate channel, or, in two dimensions, as a cylinder oriented perpendicular to the flow. The substrate is modified to include a conical hill (3D) or wedge (2D), and the test particle is located at various positions on this hill. The governing incompressible Navier–Stokes equations are solved using the finite volume fluent code. The coefficients of lift and drag are compared with the results on the flat substrate. Enhanced drag and significantly enhanced lift are observed as the test particle is situated near the summit of the hill.

## 1 Introduction

### 1.1 Introduction.

Aerosolization of particles from surfaces [1] is of particular interest to the study of the dustiness of powders [2,3]. Natural examples of this phenomenon are pollen grains or dust suspended in air, and sand particles carried away by the wind [47]. Airborne organic (e.g., soot, pollens, molds, bacteria, virus) and inorganic (e.g., silica, asbestos) dust may adversely affect health. Successful delivery of powdered pharmaceuticals presupposes efficient aerosolization of the powder [8]. Previous work [9,10] studied aerodynamics of particles located on flat surfaces.

The National Institute for Occupational Safety and Health (NIOSH) has a comprehensive program to study the dustiness of powders [2,3]. NIOSH and other laboratories [11] extensively utilize the Venturi dustiness tester (VDT) [12], which permits dustiness measurement at high Reynolds number of small quantities of expensive or potentially dangerous powders (e.g., pharmaceuticals, nanoparticles). In the VDT configuration [13], a small quantity (mg) of powder is placed in a slender tube; a high velocity air stream is forced through the tube and over the powder hill; individual particles are lifted from the powder hill and are aerosolized; the aerosol is swept into a containment chamber where it is then sampled. Dubey et al. [13] modeled the air flows in the VDT, assuming the aerosolization process to have already taken place. This work addresses the aerosolization process itself, i.e., the lift and drag of particles from the powder hill. Our numerical approach is restricted to flow in the laminar regime. Sharma et al. [14,15] have extended this work to higher Re, where the dynamics of vortex shedding from the hill becomes important.

### 1.2 Literature Review.

Many factors contribute to the detachment of particles from surfaces [16]. When a fully developed laminar boundary layer flows parallel to a substrate, the shear flow at the surface generates aerodynamic forces and moments. For the particle to detach, these aerodynamic forces must exceed the surface (e.g., capillary, adhesion) or body (e.g., gravity, van der Waals) forces which hold the particle onto the surface. The experimental measurement of these forces on isolated particles is challenging [1724]. Numerical modeling is appropriate to study such problems.

For fine and ultrafine particles (micron and nanometer diameter), the particle Reynolds number (based on the local flow velocity and particle radius) is typically less than unity. The Stokes drag (creeping flow) for a sphere in a uniform steady incompressible flow in the absence of any wall is Fdrag = 3 π μ u D, where u is the uniform freestream flow velocity, D is the sphere diameter, and μ is the fluid viscosity. Flift = 0. O'Neill [25] found that the presence of a wall increases the drag by a factor f = 1.7009. Leighton and Acrivos [26] included the first inertial correction to creeping flow and found Flift/Fdrag = 0.287 Re. This differs from the earlier Saffman [27] result, also a lowest inertial correction to creeping flow, valid for large shear flow and no wall; the corresponding 2D analog, namely, of a cylinder in linear shear flow, was studied by Bretherton [28]. Cherukat et al. [29] discuss the crossover between the Saffman and Leighton–Acrivos regimes. The effect of a finite gap between particle and wall for 2 < Re < 250 has been studied [3032]. Lee and Balachandar [33] studied the modification of lift and drag due to motion of the particle parallel to the boundary. The Finlay group [3436] examined the related problem of a sphere or cylinder attached to a wall in a Blasius boundary layer, where additional lift derives from the vertical velocity near the plate as the boundary layer develops.

Multiple laminar experiments have been conducted [3743]. The larger lift on an attached particle in a turbulent boundary layer has been studied [44,45]. The importance of aerodynamic forces on the migration and resuspension of small particles has been extensively studied [4657].

The primary focus of the earlier studies has been an understanding of the wall-induced forces on a particle. More relevant for the VDT [2,12,13] is a particle of a finely divided solid (a powder) situated in proximity to a collection of other similar particles, namely, on a powder hill. The aerodynamics of a particle when it is located on a hill has not previously been studied.

In the VDT, powder is preloaded at the base of a “tee”-shaped tube (Fig. 1), which is open to the atmosphere (at top and right) and which is inserted (towards the left) into a cubical dustiness sampling chamber. In a short (Δt = 1.5 s) dispersion phase, the pressure in the dustiness chamber is lowered, and an air stream is induced to flow over the powder hill. This air stream (Re ∼ 2 × 104) aerosolizes the powder particles from the hill, and the resulting aerosol enters the dustiness chamber (at left) as a jet. The earlier modeling study of the VDT [13] presupposed aerosolization—it considered an aerosol stream entering the dustiness chamber as a jet. In this work, we address the aerosolization process itself.

Fig. 1
Fig. 1
Close modal

As discussed later, in Sec. 2.2, the air flow in the VDT tube is turbulent (Retube ∼ 2 × 104). At these air flows, we expect vortices to be shed from the powder hill obstruction [5862]; this is seen in the work of Sharma et al. [14,15]. Vortex shedding is suppressed due to the confined geometry [6366]. Additional complications arise due to the porosity of the obstruction [67] and due to the dynamic disappearance of the obstacle (as particles are aerosolized from the hill). In the current, over-simplified, study, we restrict this investigation to low air flows (Resphere < 102) so that the flows may be treated as steady-state. We seek to understand whether the usual drag and lift of the powder particle are significantly altered from their planar values due to the obstruction geometry. In particular, we seek to determine whether it is any easier to aerosolize a particle from the hill than from a flat surface.

It is important to determine if the lift and drag forces, which contribute to the particle detachment process, and the torque, which determines the particle rotation upon detachment, are modified for a nonplanar substrate. Hence, this study considers flow over a particle attached to a hill and compares it to the idealized case of the isolated particle attached to the plane wall. This is not a true many-body study but represents a first attempt to take into account the effect of a nonplanar substrate on particle aerosolization. We note that effects of a nonplanar substrate will be important for any treatment of particle detachment from rough surfaces, which, at the microscopic scale, consist of protuberances and valleys. Length limitations restrict the discussion in this paper to lift and drag; our simulation results for torque are contained in Ref. [9].

### 1.3 Outline of This Paper.

This paper is organized as follows: Section 1 introduces dustiness, the VDT, which motivates the problem of aerosolization of particles from a hill. Section 2 discusses the mathematical formulation underlying our simulations. We also discuss (Sec. 2.1,2) the peculiar geometry of a curved surface in contact with a flat surface, which necessitates displacing the particle slightly above the substrate. Grid refinement is discussed in Sec. 2.1.6. Conditions for the relevance of our simulations to the VDT are presented in Sec. 2.2. Section 3 presents the results of these simulations: 2D (cylinder on a wedge) in Sec. 3.1 and 3D (sphere on a cone) in Sec. 3.4. A major result of our study is that the aerodynamic lift on the particle is significantly enhanced at the summit of the hill. The origin of this effect is discussed in Sec. 3.3 and further in the Sec. S4 available in the Supplemental Materials on the ASME Digital Collection. Section 4 presents our conclusions and recommendations for further work on this problem.

## 2 Mathematical Formulation and Numerical Modeling

### 2.1 Numerical Modeling.

This work studies the aerodynamic lift and drag forces on a particle positioned on a conical hill and slightly elevated above the hill surface. The particle is approximated as a sphere (in 3D) or as a cylinder (in 2D), with the cylinder orientated perpendicular to the freestream direction. All the simulations are performed in the laminar regime 1 < Rechannel < 2000, corresponding to 0.07 < Rehill < 133 and to 0.01 < Recyl/sph < 13. We assume incompressible, steady flow, and we neglect the effects of additional surface and body forces. The parameters of standard temperature and pressure (STP) air are density ρ = 1.225 kg/m3 and dynamic viscosity μ = 1.78 × 10−5 kg/m-sec; lengths are measured in mm.

#### 2.1.1 Governing Equations.

Conservation of the mass is given by
$∂ρ∂t+∂ρ ui∂xi=0$
(3.1)

where ui is the fluid velocity.

Conservation of momentum is given by
$∂ρui∂t+∂(ujρui)∂xj−∂τij∂xj+∂p∂xi=0$
(3.2)
where p is the static pressure; the stress tensor, τ, is given by
$τij= μ(∂ui∂xj+∂uj∂xi−2∂uk3∂xkδij)$
(3.3)

#### 2.1.2 Geometry

##### 2.1.2.1 Elevated versus embedded particle.

Placing the particle in contact with the plane wall introduces a singular point at the contact point. This is numerically challenging for grid generation [35]. Because of this singular point, an improperly constructed grid may contain poor quality cells (high skewness, poor orthogonality, high aspect ratios). Potential solutions to this problem are to either slightly elevate the particle above the surface or slightly embed the particle into the surface. A brief discussion of this problem is provided in the Supplemental Material (Secs. S1 for 2D and S2 for 3D) on the ASME Digital Collection. In this study, we have taken the approach to elevate the test particle slightly above the wedge/conical substrate (perpendicular distance to the cylinder/sphere center x/R = 1.05).

##### 2.1.2.2 Two-dimensional study—flow over a cylinder on a wedge in a channel.

A cylinder, oriented perpendicular to the flow, is positioned (Fig. 2) at different locations, θ, on a wedge submerged in a laminar fully developed flow. The diameter of the cylinder is D; the wedge height is L. The channel height is 15 L (so there is minimal constriction of the flow by the wedge), and the inlet/outlet boundaries are placed 15 L/20 L upstream/downstream, respectively, of the wedge center.

Fig. 2
Fig. 2
Close modal

We have performed flow simulations and determined the lift, drag, and torque for the following cases: (a) three cylinder diameters, D = 2 mm, 1 mm, and 0.5 mm, at fixed particle location, θ, wedge half-angle, α, and height, L; (b) three wedge heights, L = 6.5 mm, 13 mm, and 19 mm, at fixed cylinder diameter, D, particle location, θ, and wedge half-angle, α; (c) three wedge half-angles, α = 30 deg, 45 deg, and 60 deg, at fixed cylinder diameter D, wedge height, L, and particle location, θ.

##### 2.1.2.3 Three-dimensional study—flow over a sphere on a cone in a channel.

A sphere is positioned at different locations on a cone, submerged in a laminar, fully developed flow. The boundaries are placed as in the 2D case (Sec. 2.1.2.2) so that the upstream effects of the cone and sphere have minimal effects on the solution. The lift and drag forces and torque acting on the sphere are computed for various locations of the sphere positioned axially (θ) on and azimuthally (φ) around the cone (the latter at fixed θ).

#### 2.1.3 Solver.

The commercial code used in this study is ansysfluent v.16, which employs a finite volume method, where the conservation laws are used in their integral form. The incompressible continuity and momentum equations are solved by the semi-implicit method for pressure-linked equations (SIMPLE) algorithm, which iteratively solves for the coupled velocity and pressure fields [68]. The residual for the convergence criterion was set at 10−6. The second-order accurate implicit upwind and second-order implicit central difference schemes are used to discretize the convection and diffusion terms, respectively.

#### 2.1.4 Boundary Conditions

• No slip boundary at the particle surface, on top and bottom plates, and on the surface of the wedge/cone.

• Slip boundary on the side walls.

• Fully developed parabolic velocity profile is initialized at the channel inlet.

• Pressure boundary condition (set to atmospheric pressure) at the channel outlet.

#### 2.1.5 Mesh Generation.

Ansys Mesh Modular is used to generate structured grids for the flow geometries considered. Care has been taken to resolve wall normal velocity gradients in the boundary layer on the no-slip surfaces.

#### 2.1.6 Grid Refinement Study

##### 2.1.6.1 2D Grid refinement–cylinder on a wedge.

A detailed grid refinement analysis [69] is carried out on successively refined grids to estimate the error in the numerical results due to the discretization. The grid refinement study is constructed for the case with cylinder diameter D = 2 mm; wedge height L = 13 mm; wedge half-angle α = 60 deg; the cylinder is located at θ = 60 deg on the wedge. The area of the domain is A = 8.1737 × 104 mm2. The refinement ratio is 1.5. For the medium grid, there are 36 cells around the cylinder and six cells in the gap. The absolute drag and lift forces (not normalized to their upstream values) acting on the cylinder are computed for successive grids, and the percent changes are presented in Table 1. As the changes in these forces with mesh refinement are less than 1%, the coarsest grid was used for the reported parametric study.

Table 1

Two-dimensional grid refinement results—flow over a cylinder on a wedge

No. of cells NA/N (mm2)(A/N)1/2 (mm)rDrag% differenceLift% difference
Coarse21,6723.77161.9422.02984 × 10−71.16822 × 10−7
Med48,8001.67491.2941.52.03285 × 10−70.1481.17489 × 10−70.571
Fine109,8000.744420.86281.52.03581 × 10−70.1461.18083 × 10−70.506
No. of cells NA/N (mm2)(A/N)1/2 (mm)rDrag% differenceLift% difference
Coarse21,6723.77161.9422.02984 × 10−71.16822 × 10−7
Med48,8001.67491.2941.52.03285 × 10−70.1481.17489 × 10−70.571
Fine109,8000.744420.86281.52.03581 × 10−70.1461.18083 × 10−70.506
##### 2.1.6.2 3D grid refinement–sphere on a cone.

A detailed grid refinement analysis [69] is carried out on successively refined grids to estimate the error in the numerical results due to the discretization. Sphere diameter D = 2 mm; cone height L = 13 mm; cone half-angle α = 60 deg; the sphere is located at θ = 80 deg on the cone. The volume of the domain is V = 1.5375 × 107 mm3. Two refinements are studied: r = 1.27 (for the coarse-to-medium grid refinement); r = 1.58 (for the medium-to-fine grid refinement). For the medium grid, there are 36 cells around the sphere and six cells in the gap. The absolute drag and lift forces (not normalized to their upstream values) acting on the sphere are computed for these successively refined grids, and the percent changes are presented in Table 2. As the changes in these forces with mesh refinement are less than 1%, the coarsest grid was used for the reported studies.

Table 2

Three-dimensional grid refinement results—flow over a sphere on a cone

No. of nodes NV/N (mm3)(V/N)1/3 (mm)rDrag force (μdyne)% differenceLift force (μdyne)% difference
Coarse941,07716.3382.537470.149222.3487
Medium1,935,9837.94171.99511.2770.21650.09622.40660.259
Fine7,685,4272.00051.26001.5869.6636−0.69222.3467−0.009
No. of nodes NV/N (mm3)(V/N)1/3 (mm)rDrag force (μdyne)% differenceLift force (μdyne)% difference
Coarse941,07716.3382.537470.149222.3487
Medium1,935,9837.94171.99511.2770.21650.09622.40660.259
Fine7,685,4272.00051.26001.5869.6636−0.69222.3467−0.009

### 2.2 Comment on the Applicability of These Simulations to Powder Aerosolization in the Venturi Dustiness Tester.

The diameter of the VDT inlet nozzle is Dtube = 0.44 cm. During injection, the volumetric flow rate, Q = 60.0 L/min, corresponds to an average flow velocity uav ∼ 65.8 m/s and Reynolds number Retube ∼ 19,900. While this is clearly turbulent flow throughout the bulk of the nozzle, a particle near the wall (i.e., within the viscous sublayer), experiences a much slower flow. Using a flow profile u ∼ umax (1 − r/R)1/n, with n ∼ 6 (appropriate for Re ∼ 19,900), umax/uav = (n + 1)(2n + 1)/2n2, whence umax ∼ 83.2 m/sec. For a d ∼ 1 μm particle resting on the tube surface, the velocity experienced at its center is uparticle ∼ umax (d/D)1/n ∼ 23.1 m/sec. The ratio of particle to tube Reynolds numbers is Reparticle/Retube = (uparticle/uav)*(d/D) ∼ 8 × 10−5, whence Reparticle ∼ 1.6. Modeling the drag, lift, and torque with a local laminar flow is thus a reasonable first approach to this problem.

## 3 Simulation Results and Discussion

### 3.1 Two-Dimensional Flow Over a Cylinder on a Wedge.

The cylinder (D = 2 mm) is positioned at various locations, θ, on the wedge (L = 13 mm, wedge half-angle α = 30 deg), as discussed in Sec. 2.1.2.2. Simulations are conducted for channel Reynolds numbers Rechannel < 1920, for which the flow remains laminar, consistent with pipe flow transition at Rec ∼ 2300 [70].

Figure 3(a) displays the normalized drag, Fdrag(θ)/Fdrag(0), and Fig. 4(a) displays the normalized lift Flift(θ)/Flift(0) normalized by their upstream values (i.e., at θ = 0).

Fig. 3
Fig. 3
Close modal
Fig. 4
Fig. 4
Close modal

The normalized drag increases as the cylinder is located higher up on the wedge, attaining a maximum drag just upstream of the wedge apex (θ ∼ 90 deg). For 120 deg < θ < 150 deg, the cylinder is “shielded” by the wedge and experiences no drag. The maximum normalized drag is in the range 5–7.5.

On the other hand, the normalized lift varies dramatically. As the test cylinder is positioned at the various locations, θ, on the wedge surface, the normalized lift increases as the test cylinder moves up the upstream side of the wedge, reaching a maximum at θ ∼ 80 deg, and then decreases on the downstream side of the wedge; there is an “undershoot” (the lift is minimum and even negative for θ ∼ 100 deg) and then a “rebound” as the lift increases further downstream of the wedge. The enhancement effect is not small (in the range 8–19) and increases with increasing Re. The undershoot effect is similarly not small (in the range −5 to −16) but decreases in magnitude with increasing Re. Recall that a negative lift means that the airflow actually presses the test particle against the surface.

The maximum normalized lift is 3–4 times larger than the maximum normalized drag. Finally, while the lift experiences an undershoot (negative lift at θ ∼ 100 deg), the normalized drag exhibits no such anomaly.

### 3.2 Wedge Angle Dependence.

We now discuss the effect of varying the wedge half-angle, α. The effects of varying the wedge height, L, and the cylinder diameter, D, are discussed in the Sec. S3 available in the Supplemental Materials on the ASME Digital Collection. For the wedge height study, the normalized drag (Fig. S5 available in the Supplemental Materials on the ASME Digital Collection) and normalized lift (Fig. S6 available in the Supplemental Materials on the ASME Digital Collection) are plotted as a function of channel Re.

Fig. 5
Fig. 5
Close modal
Fig. 6
Fig. 6
Close modal

Figures 3 and 4 present results for two additional wedge half-angles α = 45 deg and 60 deg, also at D = 2 mm and L = 13 mm. Both the normalized lift and drag appear to be relatively insensitive to the magnitude of the wedge half-angle, α. The increase in normalized drag, as the test cylinder is positioned at higher locations on the wedge, is preserved for these sharper wedges (Fig. 3). The general structure in the normalized lift (maximum lift at θ ∼ 80 deg, with minimum, negative lift at θ ∼ 100 deg) is preserved (Fig. 4).

### 3.3 Origin of the Enhanced Drag and Lift.

The origin of the enhanced drag and lift forces, as the cylinder is positioned at different locations on the wedge, is examined in the Sec. S4 available in the Supplemental Materials on the ASME Digital Collection, for the case of Re = 1. Both forces result from cancelations in the pressure and shear contributions from different regions of the cylinder surface; these cancelations become less complete at higher elevations on the wedge.

Figure 5 shows that the variation in the normalized drag closely follows the variation in the normalized velocity parallel to the substrate. Similarly, the variation in the normalized lift closely follows the variation in the normalized velocity normal to the substrate (Fig. 6). While the normalized drag is amplified over the normalized parallel velocity (the dotted line is above the solid line in Fig. 5), the normalized lift is suppressed below the normalized perpendicular velocity (the dotted line is below the solid line in Fig. 6).

When the fully developed flow encounters the wedge base, a new boundary layer begins to develop as the flow turns up the wedge; the incompressible fluid is accelerated (flow constriction), with a concomitant upward vertical acceleration. Similarly, when the newly established fully developed flow (going up the wedge) encounters the apex of the wedge (flow expansion), the flow is decelerated, with a concomitant downward vertical acceleration.

### 3.4 Three-Dimensional Flow Over a Sphere on a Cone.

The flat substrate is augmented by a cone, with half-angle α = 60 deg; the simulation flow remains laminar for Rechannel < 1920, consistent with pipe transition to turbulence at Rec ∼ 2300 [70]. Figures 7 and 8 display the normalized drag and lift, where these are normalized using the lift and drag forces experienced by the sphere at θ = 0 deg.

Fig. 7
Fig. 7
Close modal
Fig. 8
Fig. 8
Close modal

The normalized drag profile (Fig. 7) for the flow over a sphere on the cone, as function of location θ, is very similar to the 2D results (Fig. 3(a)). As in 2D, maximal drag is observed at θ = 80 deg for all Reynolds numbers. The drag increases as the sphere is moved from the bottom to the peak on the upstream side of the cone and then decreases gradually to minimal drag in the “sheltered” location downstream (θ = 180 deg). For Rechannel = 1920, the maximum normalized drag is 4.5, which is smaller than the maximum normalized drag (7.4) in the 2D simulations.

The normalized lift profile (Fig. 8) for the sphere placed on the cone is similar to the 2D results (Fig. 4(a)). A maximum normalized lift occurs (θ = 80 deg) for the sphere on the upstream side of the cone apex; a minimum negative lift occurs for the sphere located on the downstream side of the cone (θ = 100 deg). This suggests that the sphere located near the top of the cone can be more easily lifted by the fluid, compared to the identical sphere located near the bottom of the cone. The variation in normalized lift is considerably larger in 3D than in 2D, e.g., maximum normalized lift (for Re = 1920) is 104 in 3D but only 28 in 2D.

The behavior near the summit (θ = 80 deg, 100 deg) is examined as a function of channel Re. The normalized drag is shown in Fig. 9; the normalized drag is shown in Fig. 10. In both cases, the Re variation is less pronounced than the variation with location on the hill.

Fig. 9
Fig. 9
Close modal
Fig. 10
Fig. 10
Close modal

### 3.5 Azimuthal Variation.

We have also studied the azimuthal variation of the lift and drag forces on the sphere; we considered 6 azimuthal angles, φ, circumambulating the cone at θ = 60 deg.

Results for the normalized drag for flow over a sphere on a cone as function of azimuthal angle, φ, at θ = 60 deg, are plotted in Fig. 11. The reference drag used for normalization is the drag experienced by the sphere located at θ = 0 deg. There is only gradual variation in the normalized drag as the sphere moves azimuthally around the cone, with variation in the range 2.4–3.9

Fig. 11
Fig. 11
Close modal

We have also studied the behavior of the normalized lift for the flow over a sphere placed at various azimuthal angles around the cone for θ = 60 deg (Fig. 12). The maximum lift occurs for the sphere on the upstream side of the cone (φ = 0 deg), and the minimum occurs for the sphere located on the downstream side of the cone (φ = 180 deg). The normalized lift varies significantly when compared to the variations in the normalized drag. The net upwards force experienced by the sphere decreases as the sphere is moved azimuthally around the cone. This suggests that a particle located on the upstream side of the hill is more easily lifted by the fluid, compared to a similar particle located at other azimuthal positions on the hill.

Fig. 12
Fig. 12
Close modal

### 3.6 Enhanced Particle Detachment From the Hill Summit.

In order for the air flow to detach a particle from a flat surface, or from a hill, the aerodynamic lift force must exceed any adhesive forces (which considerably exceed the gravitational weight for fine and ultrafine particles). Measurement of particle adhesive forces is notoriously difficult, and the literature is extensive. For the following discussion, we consider two elegant experiments: (i) use of an atomic force microscope (AFM) tip to effect detachment [71]; (ii) detaching small particles via a Huntington bar shock pulse [72]. We look for data relevant to our simulations of silica particles.

Jones et al. [71] measured the adhesive force of three relevant systems: (i) small glass sphere (r = 20 μm) on the following flat plates: (a) hydrophilic glass: 4 μN < Fadh < 7 μN; (b) hydrophilic Si: 4 μN < Fadh < 9 μN; (ii) large glass sphere (r = 100 μm) on hydrophilic glass: 0.5 μN < Fadh < 2 μN. The ranges were the result of variations in the relative humidity. Wanka et al. [72] measured the adhesive force for silica spheres on a hydrophilic-treated polystyrene substrate: (i) d = 4 μm: Fadh = 70 nN; (ii) d = 10 μm: Fadh = 150 nN.

The aerodynamic lift force is given by Flift= Clift* (1/2 ρ u2) * πD2. The particle velocity is given by u = Resphereν/D, whence Flift= (π/2) * Clift* ρν2 Resphere2 = 0.42 nN Resphere2, using the density and kinematic viscosity of air.

Aerodynamic detachment from a flat surface should occur for the AFM system for the ranges: (i) small glass (20 μm): (a) hydrophilic glass: 100 < Resphere < 130; hydrophilic Si: 100 < Resphere < 150; (ii) large glass (100 μm) from hydrophilic glass: 35 < Resphere < 70. Similarly, aerodynamic detachment from a flat surface should occur for the Huntington shock systems: (i) 4 μm silica: Resphere ∼ 13; 10 μm silica: Resphere ∼ 19. Recall (Sec. 2.2) that a micron sized particle in the boundary layer of a flat surface (using the VDT flow conditions) experiences Resphere ∼ 2, which would seem to be insufficient to detach the particles. However, for a particle at the summit of a hill, the aerodynamic lift is enhanced by a factor of 102; the required detachment Resphere decreases by a factor of 10, which is then of the right order of magnitude for the VDT to effect detachment.

## 4 Conclusions

### 4.1 Conclusions.

In this work, we studied the variation in aerodynamic forces (lift and drag) on a test particle placed at various positions on a hill (2D–cylinder located on a wedge; 3D–sphere located on a cone). We have studied variations in lift and drag due to position on the hill (θ location), Reynolds number of the flow, and geometry of the hill (cone/wedge half-angle α); in 3D there is also azimuthal variation.

In the geometric models for the hill, the cylinder and sphere were minimally offset above the substrate, with the offset (measured normally from the particle center) required to obviate the otherwise encountered geometric singularity at the contact point. The incoming flow considered is steady laminar fully developed flow. The aerodynamic forces computed for each location (θ position) are normalized using the forces experienced by the particle upstream of the hill (i.e., at θ = 0 deg). In both 2D and 3D, the lift force experienced by the particle, near the peak of the hill upstream (wedge or cone), is greater than the force experienced by the particle on the planar substrate. This suggests that particles located upstream and near the peak of the obstacle tend to be lifted more easily.

A detailed parametric study was performed for different wedge angles, heights, and cylinder diameters. The results suggest that, as the cylinder diameter decreases, the normalized lift increases, whereas the normalized drag is relatively insensitive to cylinder size. The normalized lift and drag both increase with wedge height. In 3D, forces are also computed for the sphere at azimuthal locations around the cone.

### 4.2 Recommendations for Future Work.

An immediate extension of this work would be to modify the structure of the substrate to model a rough surface, which consists, at the microscopic scale, of hills and valleys. The dislodgement of particles from the rough surface (lift and drag) can then be modeled, taking into account the variation in lift and drag for the test particle situated upstream or downstream of the protuberance. Cheng et al. [73] have reported that a significant roughness reduces the normal pull-off force to a fraction of its smooth-surface value.

Any quantitative evaluation of particle detachment must take into account the physicochemical forces (surface, van der Waals, and electrostatic) between the particle and the surface. Extensive reviews of adhesion forces are available in Refs. [18] and [20]. Microparticles must overcome adhesion forces in order to detach from a surface. The adhesion force acting on the particle surface is distributed based on the surface irregularities present on the particle surface and its contacting surface. Particle shape and composition may be more important than the size, owing to their effect on the adhesion force [74].

This work has been limited to steady laminar flow conditions. Sharma et al. [14,15] have extended these results to the turbulent regime.

Finally, the powder hill has been modeled as a monolithic cone or wedge. A realistic powder hill has nonzero porosity and permeability and might be better modeled using Darcy flow through the hill.

## Acknowledgment

We thank J. Bennett, K.H. Dunn, and D.E. Evans (all of NIOSH) for helpful discussions and for a careful reading of the paper. This work was supported in part by the NIOSH Nanotechnology Research Center. This paper reports research contained in the thesis of Palakurthi [9], which was the last thesis defense that Professor Kirti Ghia adjudicated. We recall his acute yet generous constructive criticism which enhanced the quality of that thesis. We dedicate this paper to his memory.

The findings and conclusions in this paper are those of the authors and do not necessarily represent the views of the National Institute for Occupational Safety and Health. Mention of product or company name does not constitute endorsement by the Centers for Disease Control and Prevention. None of the authors has a financial relationship with a commercial entity that has an interest in the subject of this paper.

## References

1.
Fuchs
,
N. A.
,
1964
,
The Mechanics of Aerosols
,
Pergamon
,
Oxford
, UK, Chap.
8
.
2.
Evans
,
D. E.
,
Turkevich
,
L. A.
,
Roettgers
,
C. T.
,
Deye
,
G. J.
, and
Baron
,
P. A.
,
2013
, “
Dustiness of Fine and Nanoscale Powders
,”
Ann. Occup. Hyg.
,
57
(
2
), pp.
261
277
.https://www.researchgate.net/publication/232248077_Dustiness_of_Fine_and_Nanoscale_Powders
3.
Evans
,
D. E.
,
Turkevich
,
L. A.
,
Roettgers
,
C. T.
, and
Deye
,
G. J.
,
2014
, “
Comment on Comparison of Powder Dustiness Methods
,”
Ann. Occup. Hyg.
,
58
(
4
), pp.
524
528
.10.1093/annhyg/met086
4.
Bagnold
,
R. A.
,
1954
,
The Physics of Blown Sand and Desert Dunes
,
Methuen
,
London
.
5.
Andreotti
,
B.
,
Claudin
,
P.
, and
Pouliquen
,
O.
,
2006
, “
Aeolian Sand Ripples: Experimental Study of Fully Developed States
,”
Phys. Rev. Lett.
,
96
(
2
), p.
02801
.10.1103/PhysRevLett.96.028001
6.
Duran
,
O.
,
Claudin
,
P.
, and
Andreotti
,
B.
,
2014
, “
Direct Simulation of Aeolian Sand Ripples
,”
PNAS
,
111
(
44
), pp.
15665
15668
.10.1073/pnas.1413058111
7.
Cheng
,
H.
,
Liu
,
C.
,
Li
,
J.
,
Liu
,
B.
,
Zheng
,
Z.
,
Zou
,
X.
,
Kang
,
L.
, and
Fang
,
Y.
,
2017
, “
Experimental Study of Aeolian Sand Ripples in a Wind Tunnel
,”
Earth Surf. Process. Landforms.
, 43(1), pp.
312
321
.10.1002/esp.4246
8.
Finlay
,
W. H.
,
2001
,
The Mechanics of Inhaled Pharmaceutical Aerosols: An Introduction
,
Academic
,
San Diego, CA
.
9.
Palakurthi
,
N. K.
,
2017
, “
Aerodynamics of Particle Detachment From Surfaces: A Numerical Study
,” M.Sc. thesis,
Univ. Cincinnati, Engineering & Applied Sciences, Mechanical Engineering, Cincinnati, OH
.
10.
Palakurthi
,
N. K.
,
Ghia
,
U.
, and
Turkevich
,
L. A.
,
2017
, “
Simulations of Particle Detachment From a Flat Surface
,”
ASME
Paper No. FEDSM 2017-69111.10.1115/FEDSM 2017-69111
11.
Bach
,
S.
,
Eickmann
,
U.
, and
Schmidt
,
E.
,
2013
, “
Comparison of Established Systems for Measuring the Dustiness of Powders With the UNC Dustiness Tester Developed Especially for Pharmaceutical Substances
,”
Ann. Occup. Hyg.
,
57
(
8
), pp.
1078
1086
.10.1093/annhyg/met022
12.
Boundy
,
M.
,
Leith
,
D.
, and
Polton
,
T.
,
2006
, “
Method to Evaluate the Dustiness of Pharmaceutical Powders
,”
Ann. Occup. Hyg.
,
50
(
5
), pp.
453
458
.10.1093/annhyg/mel004
13.
Dubey
,
P.
,
Ghia
,
U.
, and
Turkevich
,
L. A.
,
2017
, “
Computational Fluid Dynamics Analysis of the Venturi Dustiness Tester
,”
Powder Technol.
,
312
, pp.
310
320
.10.1016/j.powtec.2017.02.030
14.
Sharma
,
A.
,
Ghia
,
U.
, and
Turkevich
,
L. A.
,
2020
, “
Effects of Vortex Shedding on the Aerosolization of a Particle From a Hill Using Large Eddy Simulation
,”
AIAA
Paper No. 539120.10.1115/539120
15.
Sharma
,
A.
,
Ghia
,
U.
, and
Turkevich
,
L. A.
,
2020
, “
Large Eddy Simulation of Flow Over a Hemispherical Obstacle Within a Circular Tube
,”
ASME
Paper No. FEDSM2020-10596.10.1115/FEDSM2020-10596
16.
Sharma
,
M. M.
,
Chamoun
,
H.
,
Sarma
,
D. S. R.
, and
Schechter
,
R. S.
,
1992
, “
Factors Controlling the Hydrodynamic Detachment of Particles From Surfaces
,”
J. Coll. Int. Sci.
,
149
(
1
), pp.
121
134
.10.1016/0021-9797(92)90398-6
17.
Bowling
,
R. A.
,
1985
, “
An Analysis of Particle Adhesion on Semiconductor Surfaces
,”
J. Electrochem. Soc.
,
132
(
9
), pp.
2208
2214
.10.1149/1.2114320
18.
Bowling
,
R. A.
,
1988
, “
A Theoretical Review of Particle Adhesion in Particles on Surfaces 1: Detection, Adhesion
,”
Removal
,
K. L.
Mittal
, ed.,
Springer
,
Boston, MA
, pp.
129
142
.
19.
Willetts
,
B. B.
, and
Naddeh
,
K. F.
,
1986
, “
Measurements of Lift on Spheres Fixed in Low Reynolds Number Flows
,”
J. Hydraul. Res.
,
24
(
5
), pp.
425
435
.10.1080/00221688609499318
20.
Ranade
,
M. B.
,
1987
, “
Adhesion and Removal of Fine Particles on Surfaces
,”
Aerosol Sci. Tech.
,
7
(
2
), pp.
161
176
.10.1080/02786828708959155
21.
Yiantsios
,
S. G.
, and
Karabelas
,
A. J.
,
1995
, “
Detachment of Spherical Microparticles Adhering on Flat Surfaces by Hydrodynamic Forces
,”
J. Coll. Int. Sci.
,
176
(
1
), pp.
74
85
.10.1006/jcis.1995.0009
22.
Hatakentaka
,
K.
,
Tanaka
,
K.
,
Fuchiwaki
,
M.
, and
Shimizu
,
F.
,
2003
, “
Drag and Lift Acting on a Circular Cylinder in Laminar Boundary Layer: Part 1—Verification of Force Transducer for Measuring Small Force
,”
ASME
Paper No. FEDSM2003-45628.10.1115/FEDSM2003-45628
23.
Essawey
,
A. H. I.
,
2004
, “
Microparticle Detachment From Surfaces by Fluid Flow
,” Doctoral dissertation,
University of Notre Dame
, Notre Dame, IN.
24.
Ibrahim
,
A. H.
,
Dunn
,
P. F.
, and
Qazi
,
M. F.
,
2008
, “
Experiments and Validation of a Model for Microparticle Detachment From a Surface by Turbulent Air Flow
,”
J. Aerosol Sci.
,
39
(
8
), pp.
645
656
.10.1016/j.jaerosci.2008.03.006
25.
O'Neill
,
M. E.
,
1968
, “
A Sphere in Contact With a Plane Wall in a Slow Linear Shear Flow
,”
Chem. Eng. Sci.
,
23
(
11
), pp.
1293
1298
.10.1016/0009-2509(68)89039-6
26.
Leighton
,
D.
, and
Acrivos
,
A.
,
1985
, “
The Lift on a Small Sphere Touching a Plane in the Presence of a Simple Shear Flow
,”
Z. Angew. Math. Phys.
,
36
(
1
), pp.
174
178
.10.1007/BF00949042
27.
Saffman
,
P. G. T.
,
1965
, “
The Lift on a Small Sphere in a Slow Shear Flow
,”
J. Fluid Mech.
,
22
(
2
), pp.
385
400
(corrigendum Vol. 31, p. 624 (1968)).10.1017/S0022112065000824
28.
Bretherton
,
F. P.
,
1962
, “
Slow Viscous Motion Round a Cylinder in Simple Shear
,”
J. Fluid Mech.
,
12
(
4
), pp.
591
613
.10.1017/S0022112062000415
29.
Cherukat
,
P.
,
McLaughlin
,
J. B.
, and
Dandy
,
D. S.
,
1999
, “
A Computational Study of the Inertial Lift on a Sphere in a Linear Shear Flow Field
,”
Int. J. Multiphase Flow
,
25
(
1
), pp.
15
33
.10.1016/S0301-9322(98)00034-2
30.
Cherukat
,
P.
, and
McLaughlin
,
J. B.
,
1994
, “
The Inertial Lift on a Rigid Sphere in a Linear Shear Flow Field Near a Flat Wall
,”
J. Fluid Mech.
,
263
, pp.
1
18
.10.1017/S0022112094004015
31.
Cherukat
,
P.
,
McLaughlin
,
J. B.
, and
Graham
,
A. L.
,
1994
, “
The Inertial Lift on a Rigid Sphere Translating in a Linear Shear Flow Field
,”
Int. J. Multiphase Flow
,
20
(
2
), pp.
339
353
.10.1016/0301-9322(94)90086-8
32.
Zeng
,
L.
,
Najjar
,
F.
,
Balachandar
,
S.
, and
Fischer
,
P.
,
2009
, “
Forces on a Finite-Sized Particle Located Close to a Wall in a Linear Shear Flow
,”
Phys. Fluids
,
21
(
3
), p.
033302
.10.1063/1.3082232
33.
Lee
,
H.
, and
Balachandar
,
S.
,
2010
, “
Drag and Lift Forces on a Spherical Particle Moving on a Wall in a Shear Flow at Finite Re
,”
J. Fluid Mech.
,
657
, pp.
89
125
.10.1017/S0022112010001382
34.
Sweeney
,
L. G.
,
2006
, “
Blasius Boundary Layer Flow Over a Cylinder and Sphere Attached to a Wall: Characterization of Lift and Drag Forces
,” M.Sc. thesis,
University of Alberta
, Edmonton, AB.
35.
Sweeney
,
L. G.
, and
Finlay
,
W. H.
,
2007
, “
Lift and Drag Forces on a Sphere Attached to a Wall in a Blasius Boundary Layer
,”
J. Aerosol Sci.
,
38
(
1
), pp.
131
135
.10.1016/j.jaerosci.2006.09.006
36.
Martinez
,
R. C.
,
Sweeney
,
L. G.
, and
Finlay
,
W. H.
,
2009
, “
Aerodynamic Forces and Moment on a Sphere or Cylinder Attached to a Wall in a Blasius Boundary Layer
,”
Eng. App. Comp. Fluid Mech.
,
3
(
3
), pp.
289
295
.10.1080/19942060.2009.11015272
37.
Goldman
,
A. J.
,
Cox
,
R. G.
, and
Brenner
,
H.
,
1967
, “
Slow Viscous Motion of a Sphere Parallel to a Plane Wall: I—Motion Through a Quiescent Fluid
,”
Chem. Eng. Sci.
,
22
(
4
), pp.
637
651
.10.1016/0009-2509(67)80047-2
38.
Goldman
,
A. J.
,
Cox
,
R. G.
, and
Brenner
,
H.
,
1967
, “
Slow Viscous Motion of a Sphere Parallel to a Plane Wall: II—Couette Flow
,”
Chem. Eng. Sci.
,
22
(
4
), pp.
653
660
.10.1016/0009-2509(67)80048-4
39.
Bagnold
,
R. A.
,
1974
, “
Fluid Forces on a Body in Shear Flow; Experimental Use of Stationary Flow
,”
Proc. Roy. Soc. London A
,
340
(
1621
), pp.
147
171
.10.1098/rspa.1974.0145
40.
King
,
M. R.
, and
Leighton
, and
D. T.
Jr.
,
1997
, “
Measurement of the Inertial Lift on a Moving Sphere in Contact With a Plane Wall in a Shear Flow
,”
Phys. Fluids
,
9
(
5
), pp.
1248
1255
.10.1063/1.869264
41.
Dandy
,
D. S.
, and
Dwyer
,
H. A.
,
1990
, “
A Sphere in Shear Flow at Finite Reynolds Number: Effect of Shear on Particle Lift, Drag, and Heat Transfer
,”
J. Fluid Mech.
,
216
, pp.
381
410
.10.1017/S0022112090000477
42.
Burdick
,
G. M.
,
Berman
,
N. S.
, and
Beaudoin
,
S. P.
,
2001
, “
Describing Hydrodynamic Particle Removal From Surfaces Using the Particle Reynolds Number
,”
J. Nanopart. Res.
,
3
(
5/6
), pp.
453
465
.10.1023/A:1012593318108
43.
Muthanna
,
C.
,
Nieuwstadt
,
F. T. M.
, and
Hunt
,
J. C. R.
,
2005
, “
Measurement of the Aerodynamic Forces on a Small Particle Attached to a Wall
,”
Exp. Fluids
,
39
(
2
), pp.
455
463
.10.1007/s00348-005-0961-0
44.
Hall
,
D.
,
1988
, “
Measurement of the Mean Force on a Particle Near a Boundary in Turbulent Flow
,”
J. Fluid Mech.
,
187
, pp.
451
466
.10.1017/S0022112088000515
45.
Mollinger
,
A. M.
, and
Nieuwstadt
,
F. T. M.
,
1996
, “
Measurement of the Lift Force on a Particle Fixed to the Wall in the Viscous Sublayer of a Fully Developed Turbulent Boundary Layer
,”
J. Fluid Mech.
,
316
, pp.
285
306
.10.1017/S0022112096000547
46.
Segre
,
G.
, and
Silberberg
,
A.
,
1961
, “
Radial Particle Displacements in Poiseuille Flow of Suspensions
,”
Nature
,
189
(
4760
), pp.
209
210
.10.1038/189209a0
47.
Segre
,
G.
, and
Silberberg
,
A.
,
1962
, “
Behavior of Macroscopic Rigid Spheres in Poiseuille Flow I, II
,”
Fluid Mech.
,
14
, pp.
115
135
(136–157).10.1017/S002211206200110X
48.
Masironi
,
L. A.
, and
Fish
,
B. R.
,
1964
, “
Direct Observation of Particle Reentrainment From Surfaces
,”
Surface Contamination
,
B. R.
Fish
, ed.,
Pergamon
, Oxford, UK, pp.
55
59
.
49.
Cox
,
R. G.
, and
Brenner
,
H.
,
1968
, “
The Lateral Migration of Solid Particles in Poiseuille Flow: I. Theory
,”
Chem. Eng. Sci.
,
23
(
2
), pp.
147
173
.10.1016/0009-2509(68)87059-9
50.
Ho
,
B. P.
, and
Leal
,
R. G.
,
1974
, “
Inertial Migration of Rigid Spheres in Two-Dimensional Unidirectional Flows
,”
J. Fluid Mech.
,
65
(
2
), pp.
365
400
.10.1017/S0022112074001431
51.
Cox
,
R. G.
, and
Hsu
,
S. K.
,
1977
, “
The Lateral Motion of Solid Particles in a Laminar Flow Near a Plane
,”
Int. J. Multiphase Flow
,
3
(
3
), pp.
201
222
.10.1016/0301-9322(77)90001-5
52.
Reeks
,
M. W.
,
Reed
,
J.
, and
Hall
,
D.
,
1988
, “
On the Resuspension of Small Particles by a Turbulent Flow
,”
J. Phys. D
,
21
(
4
), pp.
574
589
.10.1088/0022-3727/21/4/006
53.
Schonberg
,
J. A.
, and
Hinch
,
E. J.
,
1989
, “
The Inertial Migration of a Sphere in Poiseuille Flow
,”
J. Fluid Mech.
,
203
, pp.
517
524
.10.1017/S0022112089001564
54.
McLaughlin
,
J. B.
,
1991
, “
Inertial Migration of a Small Sphere in Linear Shear Flows
,”
J. Fluid Mech.
,
224
, pp.
261
274
.10.1017/S0022112091001751
55.
Ziskind
,
G.
,
Fichman
,
M.
, and
Gutfinger
,
C.
,
1997
, “
Adhesion Moment Model for Estimating Particle Detachment From a Surface
,”
J. Aerosol Sci.
,
28
(
4
), pp.
623
634
.10.1016/S0021-8502(96)00460-0
56.
Ziskind
,
G.
,
Fichman
,
M.
, and
Gutfinger
,
C.
,
1998
, “
Effects of Shear on Particle Motion Near a surface—Application to Resuspension
,”
J. Aerosol Sci.
,
29
(
3
), pp.
323
338
.10.1016/S0021-8502(97)10009-X
57.
Zhang
,
F.
,
Reeks
,
M.
, and
Kissane
,
M.
,
2013
, “
Particle Resuspension in Turbulent Boundary Layers and the Influence of Non-Gaussian Removal Forces
,”
J. Aerosol Sci.
,
58
, pp.
103
128
.10.1016/j.jaerosci.2012.11.009
58.
Taneda
,
J.
,
1956
, “
Experimental Investigation of the Wake Behind a Sphere at Low Reynolds Number
,”
J. Phys. Soc. Jpn.
,
11
(
10
), pp.
1104
1108
.10.1143/JPSJ.11.1104
59.
Magarvey
,
R.
, and
Bishop
,
R. L.
,
1961
, “
Transition Ranges for Three-Dimensional Wakes
,”
Can. J. Phys.
,
39
(
10
), pp.
1418
1422
.10.1139/p61-169
60.
Tomboulides, A. G., Orszag
,
S.
, and
Karniadakis
,
G. E.
,
1993
, “
Direct and Large Eddy Simulations of Axisymmetric Wakes
,”
AIAA
Paper No. 93–9546.10.2514/6.1993-546
61.
Johnson
,
T. A.
, and
Patel
,
V. C.
,
1999
, “
Flow Past a Sphere Up to a Reynolds Number of 300
,”
J. Fluid Mech.
,
378
, pp.
19
70
.10.1017/S0022112098003206
62.
Przadka
,
W.
,
Miedzik
,
J.
,
Gumowski
,
K.
,
Goujon-Durand
,
S.
, and
Wesfreid
,
W.
,
2008
, “
The Wake Behind the Sphere: Analysis of Vortices During Transition From Steadiness to Unsteadiness
,”
Arch. Mech.
,
60
, pp.
467
474
.https://am.ippt.pan.pl/am/article/viewFile/v60p467/pdf
63.
Savory
,
E.
, and
Toy
,
N.
,
1986
, “
Hemispheres and Hemisphere-Cylinders in Turbulent Boundary Layers
,”
J. Wind Eng. Ind. Aerodyn.
,
23
, pp.
345
364
.10.1016/0167-6105(86)90054-1
64.
Savory
,
E.
, and
Toy
,
N.
,
1986
, “
The Flow Regime in the Turbulent Near Wake of a Hemisphere
,”
Exp. Fluids
,
4
(
4
), pp.
181
188
.10.1007/BF00717812
65.
Fedrizzi
,
M.
,
Giacobello
,
M.
,
Soria
,
J.
,
Atkinson
,
C.
, and
Jones
,
M.
,
2012
, “
Experimental Investigation of a Hemisphere in a Thin Flat Plate Boundary Layer
,”
Proceedings 18th Australasian Fluid Mechanics Conference
, Dec. 3–7, Launceston, Australia, pp.
18
21
.
66.
Wood
,
J. N.
,
De Nayer
,
G.
,
Schmidt
,
S.
, and
Breuer
,
M.
,
2016
, “
Experimental Investigation and Large-Eddy Simulation of the Turbulent Flow Past a Smooth and Rigid Hemisphere
,”
Flow, Turbul. Combust.
,
97
(
1
), pp.
79
119
.10.1007/s10494-015-9690-5
67.
Ledda
,
P. G.
,
Siconolfi
,
L.
,
Viola
,
F.
,
Gallaire
,
F.
, and
Camarri
,
S.
,
2018
, “
Suppression of Von Karman Vortex Streets Past Porous Rectangular Cylinders
,”
Phys. Rev. Fluids
,
3
(
10
), p.
103901
.10.1103/PhysRevFluids.3.103901
68.
Patankar
,
S. V.
, and
Spalding
,
D. B.
,
1972
, “
A Calculation Procedure for Heat, Mass and Momentum Transfer in Three-Dimensional Parabolic Flows
,”
Int. J. Heat Mass Transfer
,
15
(
10
), pp.
1787
1806
.10.1016/0017-9310(72)90054-3
69.
Celik
,
I. B.
,
Ghia
,
U.
, and
Roache
,
P. J.
,
2008
, “
Procedure for Estimation and Reporting of Uncertainty Due to Discretization in CFD Applications
,”
ASME J. Fluids Eng.
,
130
(
7
), p.
078001
.10.1115/1-2960953
70.
Schlichting
,
H.
,
1979
,
Boundary-Layer Theory
,
7
th ed.,
McGraw-Hill
,
New York
(tr. by J. Kestin from the German Grenzschlicht-Theorie (Karlsruhe: Braun)).
71.
Jones
,
R.
,
Pollock
,
H. M.
,
Cleaver
,
J. A. S.
, and
Hodges
,
C. S.
,
2002
, “
Adhesion Forces Between Glass and Silicon Surfaces in Air Studied by AFM: Effects of Relative Humidity Particle Size, Roughness, and Surface Treatment
,”
Langmuir
,
18
(
21
), pp.
8045
8055
.10.1021/la0259196
72.
Wanka
,
S.
,
Kappl
,
M.
,
Wolkenhauer
,
M.
, and
Butt
,
H.-J.
,
2013
, “
Measuring Adhesion Forces in Powder Collectives by Inertial Detachment
,”
Langmuir
,
29
(
52
), pp.
16075
16083
.10.1021/la404014d
73.
Cheng
,
W.
,
Dunn
,
P. F.
, and
Brach
,
R. M.
,
2002
, “
Surface Roughness Effects on Microparticle Adhesion
,”
J. Adhes.
,
78
(
11
), pp.
929
965
.10.1080/00218460214510
74.
Wu
,
Y.
,
Cliff
,
I. D.
, and
Russell
,
A. G.
,
1992
, “
Controlled Wind Tunnel Experiments for Particle Bounce-Off and Resuspension
,”
Aerosol Sci. Technol
,.
17
(
4
), pp.
245
262
.10.1080/02786829208959574