## 1. Introduction

The problem of the removal of aerosol particles from gas streams has become of increasing importance from the standpoint of public health and the recovery of valuable products. Technology of controlling the aerosol particles or improving the liquid phase of aerosol is very important in many industrial processes such as oil and petroleum, electronic, mining, and food, as well as waste products like noxious emission of aerosol in chemical plants. There are several ways for this purpose among which fibrous filters are more popular so that it is obvious to try to improve their efficiency. The efficiency of collection and the pressure drop are the most important practical considerations in the design of these fibrous filters [2]. Various experimental studies on liquid aerosol filtration have shown a filter clogging in several stages leading to a drainage stage with a constant pressure drop [3, 4]. Contal et al. [3] introduced four stages describing the clogging of fibrous filters by liquid droplets: In the first stage, the droplets impact the wet fibers. In the second stage, the amount of those droplets captured by fiber increases so that neighbor droplets coalesce. In the third stage, liquid shells form in the net that leads to a large increase of the pressure drop and to a dramatic decrease in the efficiency of the filter. Finally, there will be a pseudo-stationary state between the droplet collection and gravitational drainage of the liquid phase [1]. The first stage of clogging is studied in many researches experimentally. Patel et al. [5] investigated the drop breakup in a flow through fiber filters and the breakup probability for a drop. Hung and Yao [6] evaluated the impact of water droplets on a cylindrical object experimentally. Both these studies show that depending on the speed of the impacting droplets, the drops can be captured and broken into several pieces. Lorenceau et al. [1] investigated the threshold velocity and radius of drops captured by thin fibers. They have shown that the threshold impact velocity of drops decreases as the drop radius increases and vice versa. Following this study, Lorenceau et al. [7] investigated off-center impact of droplets on horizontal fibers.

Most of liquid droplets in aerosols produced in different industrial and natural processes show non-Newtonian behavior. Therefore, studying the phenomenon of impaction of a non-Newtonian droplet on thin fibers is of prime importance. All mentioned studies have been devoted to Newtonian fluids. Colliding of non-Newtonian droplets on thin fibers, to our knowledge, has not been investigated previously. However, colliding of non-Newtonian drops on flat plates is evaluated in some experimental and numerical studies. Haeri and Hashemabadi [8] studied spreading of power-law fluids on inclined plates both numerically and experimentally. Saïdi et al. [9] investigated experimentally the influence of yield stress on the fluid droplet impact control. Son and Kim [10] studied experimentally the spreading of inkjet droplet of non-Newtonian fluid on a solid surface with controlled contact angle at low Weber and Reynolds numbers. Kim and Baek [11] investigated the parameters that govern the impact dynamics of yield stress on the fluid droplets on solid surfaces.

In this chapter, we focus on the first stage of clogging and investigate the impaction of non-Newtonian power-law fluids on thin fibers. We aim to obtain the threshold radius of impacting droplets in different impact velocities. Effect of shear-thinning and shear-thickening behavior of droplets is evaluated and compared with corresponding Newtonian fluids. For this purpose, volume of fluid method is used and open source OpenFOAM software is applied for simulations.

The order of contents in this chapter is as follows: In Section 2 the governing equations and numerical methodology are given. Validation of the obtained results for Newtonian fluids in comparison with experimental observations presented by Lorenceau et al. [1] is provided in Section 3.1. A general description of the observations for non-Newtonian droplets captured by horizontal fiber is presented in Section 3.2, and the behavior of non-Newtonian drops affecting thin fibers is explained. Finally, a summary of the results and conclusions is provided.

## 2. Methodology

### 2.1. Volume of fluid method

Hirt and Nichols [12] demonstrated the volume of fluid (VOF) method and started a new trend in multiphase flow simulation. It relies on the definition of an indicator function

where *p* is pressure, and

where

In this study, a modified approach similar to the one proposed in [13] is used.

### 2.2. Power-law model

Generally, the form of the function relating

is adequate for many non-Newtonian fluids [16], where *K* is the fluid consistency coefficient and *n* is power-law index. This is a well-known power-law model, which is used in this study.

### 2.3. Computational method

All computations are performed using the code OpenFOAM [17], an open source computational fluid dynamics (CFD) toolbox, utilizing a cell-center-based finite volume method on a fixed unstructured numerical grid and employing the solution procedure based on the pressure implicit with splitting of operators (PISO) algorithm for coupling between pressure and velocity in transient flows [18].

A grid spacing equal to 1/20 of the droplet radius was used to discretize the computational domain. The mesh size was determined based on a mesh refinement study in which the grid spacing was progressively decreased until further reductions made no significant change in the predicted droplet shape during the impact. Bussmann et al. [19] have provided a detailed description of such a mesh refinement study earlier. The mesh size was uniform throughout the computational domain, which included the droplet and fiber.

## 3. Results and discussion

### 3.1. Validation of the Solution

In this section, we present the results obtained for water droplets impacting a fiber of radius

Defining

where *b* is the radius of the fiber, and also characteristic velocity of the drop as:

Variation of the dimensionless threshold radius versus dimensionless impact velocity is plotted in Figure 2, which shows a good agreement with the experimental and theoretical data presented by Lorenceau et al. [1].

Errors of the achieved results in comparison with experimental results of Lorenceau et al. (2004) are presented in Table 2.

### 3.2. Results for non-newtonian droplets

Four forces, capillary, gravity, inertia, and viscosity, are important during the impaction of droplet to the fiber. Capillary and gravity forces are in contrast. While the weight of the drop tends to detach it from the fiber, surface tension is responsible for the sustentation. In order to analyze the dynamics of the capture of a drop by a fiber, viscous and inertia effects must be considered. Here again, those effects are antagonistic. While inertia tends to make the drop cross the fiber, viscous effects will induce a dissipation and can thus possibly stop the drop [1]. Effects of variation of viscosity were not investigated in the previous works. In order to investigate these effects, both shear-thinning and shear-thickening fluids will be considered.

At first, we present general observations of threshold radiuses at different velocities for shear-thinning droplets,

Next, we present our observations for the two kinds of shear-thickening droplets, *n*=1.5, 2, impacting the fiber which are shown in Figures 4 and 5.

Results show that, on the one hand, in shear-thinning drops, the threshold radius decreased in a fixed velocity in comparison with the corresponding Newtonian drops. Moreover, instabilities were observed during the impaction of the drop. This is because viscosity in shear-thinning fluids decreases as it has an inverse relation with shear rate. That is, as a shear-thinning droplet impacts a thin fiber, due to high shear rates, viscosity decreases in comparison with the corresponding Newtonian fluid, so that the viscosity force reduces (i.e., a force that helps the drop to stick on the fiber is reduced). Thus, in order to balance the antagonistic forces, drop size has to be decreased, which means the threshold radius has to be decreased. On the other hand, in the case where the drop is a shear-thickening fluid, the threshold radius increased in a fixed velocity in comparison with the corresponding Newtonian fluid, since the viscosity of shear-thickening fluids has a direct relation to the shear rate. It means that, viscosity increases as the drop impacts the fiber compared to Newtonian drops. Thereby, in shear-thickening drops, viscous forces increase and the drop’s tendency to stick to the fiber increases. Therefore, the threshold drop size will also increase. This tendency is obvious in Figures 4 and 5. Threshold radius of Newtonian, shear-thinning and shear-thickening droplets versus impact velocity is plotted in Figure 6.

Figure 6 shows that the threshold radius will decrease generally with the increase of an impact velocity or vice versa. In shear-thinning droplets, viscosity has an inverse relation with shear rate, as mentioned before, so that it will decrease with the increase of impact velocity. It means inertia forces increase whereas viscosity forces decrease. That is why instabilities are observed in shear-thinning droplets impacting fiber at high speeds. The threshold radius of shear-thickening droplets, also, decreases with the increase of impact velocity. In shear-thickening droplets, although there is a direct relation between the viscosity and the shear rate, velocity increases with second power compared to capillary forces. That is, increase of inertia forces is higher than increase of viscosity forces so that the threshold radius will decrease to make the drop stick to the fiber. However, in higher velocities, this rate will decrease as it can be seen in Figure 6. In other words, rate of the decreasing of threshold radius with the increasing of impact velocity will decrease due to high growth of shear rate.

There are some general observations that are common to Newtonian and non-Newtonian fluids. For all kinds of fluids, the threshold radius of the droplet has a reverse relation with the impact velocity of the droplet. In all cases, at a fixed impact velocity, droplets with a radius greater than the threshold radius have passed the fiber without breakup, and drops with a radius lower than the threshold radius clung to the fiber entirely.

## 4. Conclusion

Impaction of non-Newtonian power-law droplet to the horizontal fiber of circular cross section is investigated in this study. Volume of fluid technique is employed, significantly reducing the computational cost. Outcomes are divided into three parts: First, it has been observed that the threshold radius of the droplets decreased with the increase of impact velocity for Newtonian, shear-thinning and shear-thickening fluids, due to the growth of inertia forces versus viscosity and capillary forces. Second, instabilities have been observed in shear-thinning droplets at high impact velocities due to severe reduction of viscosity. Finally, in shear-thickening droplets, at high impact speeds, rate of reduction of threshold radius has decreased due to increase of the viscosity forces.