A fully Eulerian approach to particle inertial deposition in a physiologically realistic bifurcation

ABSTRACT The investigation of aerosol flows in bifurcations is significant in biomedical applications, because this geometry resembles the branching airways of the lung. Most of the existing computational studies regarding aerosol flow in branching airways use Lagrangian description for the particulate phase, when particle inertial effects are taken into consideration. \textcolor{blue}{In the present study a fully Eulerian model is used for studying the transport and deposition of inertial particles under steady state inspiratory flow in the single physiologically realistic bifurcation created by generations G3-G4 of the human lung}. In-house codes are used to create the geometry and the computational grid, obtain the fluid flow field and calculate particle concentration. Deposition fractions are calculated, concentration profiles are shown and deposition sites are indicated under different flow conditions (Reynolds number and flow asymmetry in the branches) and particles of various sizes (1-10$\mu$m). It is found that total particle deposition fraction is loosely dependent on fluid flow Reynolds number and our results are in agreement with experimental findings. Moreover, total deposition fraction does not change significantly between symmetric $Q_1/Q_2=1$ and asymmetric $Q_1/Q_2=2$ flow conditions, but is considerably lower for the totally obstructed case, $Q_2=0$. On the contrary, deposition sites and particle concentration profiles depend on the various flow conditions. Apart from deposition at the bifurcation, which is present at all cases, deposition sites downstream the bifurcation move towards the outer bifurcation wall and daughter tube exit as Reynolds number increases. For the $Q_2=0$, it is also shown that particles are trapped and deposit even in the obstructed daughter tube. In all studied cases, it is found that there is a particle-free region at the outer wall, at the beginning of the daughter tubes opposite the bifurcation. The straightforward formulation of deposition modelling and derivation of detailed concentration profiles along the bifurcation are major advantages of the fully Eulerian methodology used in the study.