Elastoviscoplastic effects on liquid plug propagation and rupture occurring in airways are studied computationally using the Oldroyd-B and Saramito–Herschel–Bulkley models. The relevant parameters are selected from physiological values representative of the eighth-to-tenth generation branches of a typical adult lung. The respiration pushes the liquid plug, depositing a trailing film thicker than the leading film. As a result, the liquid plug gets drained and eventually ruptures. We model airway reopening considering a rigid axisymmetric tube whose inner surface is coated by a thin non-Newtonian liquid film. A critical elastic behaviour is revealed: for low Weissenberg number (subcritical), the viscoelastic stress is released in the liquid plug, while for high Weissenberg number (supercritical), the stretched polymeric chains release their stresses in the trailing film, giving rise to (i) hoop stress that increases the film thickness and (ii) axial stress that leads to a speed-up of the liquid plug. Under supercritical conditions, we also identify a resonance that amplifies the elastic stresses. A mechanical analogy is proposed to elucidate the resonance phenomenon. The occurrence of the resonance is robust upon a variation of Weissenberg number, Laplace number, reference solvent-to-total dynamic viscosity ratio, the surfactant elastoviscoplastic mucus. Our simulations confirm that a presence of surfactants do not significantly affect the results, except for the expected delay of airway reopening due to air–mucus surface contamination. Such a novel elastocapillary mechanism increases the risk of epithelial cell damage regardless of the occurrence of plug rupture.

Interface-resolved direct numerical simulations are performed to investigate bubble-induced transition from a laminar to elasto-inertial turbulent (EIT) state in a pressure-driven viscoelastic square channel flow. The Giesekus model is used to account for the viscoelasticity of the continuous phase, while the dispersed phase is Newtonian. Simulations are performed for both single- and two-phase flows for a wide range of Reynolds (${Re}$) and Weissenberg (${\textit{Wi}}$) numbers. In the absence of any discrete external perturbations, single-phase viscoelastic flow is transitioned to an EIT regime at a critical Weissenberg number ($Wi_{cr})$ that decreases with increasing ${Re}$. It is demonstrated that injection of bubbles into a laminar viscoelastic flow introduces streamline curvature that is sufficient to trigger an elastic instability leading to a transition to an EIT regime. The temporal turbulent kinetic energy spectrum shows a scaling of $-2$ for this multiphase EIT regime, and this scaling is found to be independent of size and number of bubbles injected into the flow. It is also observed that bubbles move towards the channel centreline and form a string-shaped alignment pattern in the core region at the lower values of ${Re}=10$ and ${\textit{Wi}}=1$. In this regime, there are disturbances in the core region in the vicinity of bubbles while flow remains essentially laminar. Unlike the solid particles, it is found that increasing shear-thinning effect breaks up the alignment of bubbles.

Cross-stream migration of a deformable fluid particle is investigated computationally in a pressure-driven channel flow of a viscoelastic fluid via interface-resolved simulations. Flow equations are solved fully coupled with the Giesekus model equations using an Eulerian–Lagrangian method and extensive simulations are performed for a wide range of flow parameters to reveal the effects of particle deformability, fluid elasticity, shear thinning and fluid inertia on the particle migration dynamics. Migration rate of a deformable particle is found to be much higher than that of a solid particle under similar flow conditions mainly due to the free-slip condition on its surface. It is observed that the direction of particle migration can be altered by varying shear thinning of the ambient fluid. With a strong shear thinning, the particle migrates towards the wall while it migrates towards the channel centre in a purely elastic fluid without shear thinning. An onset of elastic flow instability is observed beyond a critical Weissenberg number, which in turn causes a path instability even for a nearly spherical particle. An inertial path instability is also observed once particle deformation exceeds a critical value. Shear thinning is found to be suppressing the path instability in a viscoelastic fluid with a high polymer concentration whereas it reverses its role and promotes path instability in a dilute polymer solution. It is found that migration of a deformable particle towards the wall induces a secondary flow with a velocity that is approximately an order of magnitude higher than the one induced by a solid particle under similar flow conditions.

Direct numerical simulations are carried out to study the effect of finite Weissenberg number up to $Wi=16$ on laminar and turbulent channel flows of an elastoviscoplastic (EVP) fluid, at a fixed bulk Reynolds number of $2800$. The incompressible flow equations are coupled with the evolution equation for the EVP stress tensor by a modified Saramito model that extends both the Bingham viscoplastic and the finite extensible nonlinear elastic-Peterlin (FENE-P) viscoelastic models. In turbulent flow, we find that drag decreases with both the Bingham and Weissenberg numbers, until the flow laminarises at high enough elastic and yield stresses. Hence, a higher drag reduction is achieved than in the viscoelastic flow at the same Weissenberg number. The drag reduction persists at Bingham numbers up to 20, in contrast to viscoplastic flow, where the drag increases in the laminar regime compared with a Newtonian flow. Moreover, elasticity affects the laminarisation of an EVP flow in a non-monotonic fashion, delaying it at lower and promoting it at higher Weissenberg numbers. A hibernation phenomenon is observed in the EVP flow, leading to large changes in the unyielded regions. Finally, plasticity is observed to affect both low- and high-speed streaks equally, attenuating the turbulent dissipation and the fragmentation of turbulent structures.

We present numerical simulations of laminar and turbulent channel flow of an elastoviscoplastic fluid. The non-Newtonian flow is simulated by solving the full incompressible Navier–Stokes equations coupled with the evolution equation for the elastoviscoplastic stress tensor. The laminar simulations are carried out for a wide range of Reynolds numbers, Bingham numbers and ratios of the fluid and total viscosity, while the turbulent flow simulations are performed at a fixed bulk Reynolds number equal to 2800 and weak elasticity. We show that in the laminar flow regime the friction factor increases monotonically with the Bingham number (yield stress) and decreases with the viscosity ratio, while in the turbulent regime the friction factor is almost independent of the viscosity ratio and decreases with the Bingham number, until the flow eventually returns to a fully laminar condition for large enough yield stresses. Three main regimes are found in the turbulent case, depending on the Bingham number: for low values, the friction Reynolds number and the turbulent flow statistics only slightly differ from those of a Newtonian fluid; for intermediate values of the Bingham number, the fluctuations increase and the inertial equilibrium range is lost. Finally, for higher values the flow completely laminarizes. These different behaviours are associated with a progressive increases of the volume where the fluid is not yielded, growing from the centreline towards the walls as the Bingham number increases. The unyielded region interacts with the near-wall structures, forming preferentially above the high-speed streaks. In particular, the near-wall streaks and the associated quasi-streamwise vortices are strongly enhanced in an highly elastoviscoplastic fluid and the flow becomes more correlated in the streamwise direction.