Experimental studies of natural convection in yield stress fluids have revealed transient behaviours that contradict predictions from viscoplastic models. For example, at a sufficiently large yield stress, these models predict complete motionlessness; below a critical value, yielding and motion onset can be delayed in viscoplastic models. In both cases, however, experiments observe immediate motion onset. We present numerical simulations of the transient natural convection of elastoviscoplastic (EVP) fluids in a square cavity with differentially heated side walls, exploring the role of elasticity in reconciling theoretical predictions with experimental observations. We consider motion onset in EVP fluids under two initial temperature distributions: (i) a linear distribution characteristic of steady pure conduction, and (ii) a uniform distribution representative of experimental conditions. The Saramito EVP model exhibits an asymptotic behaviour similar to the Kelvin-Voigt model as, where material behaviour is primarily governed by elasticity and solvent viscosity. The distinction between motion onset and yielding, a hallmark of EVP models, is the key feature that bridges theoretical predictions with experimental observations. While motion onset is consistently immediate (as seen in experiments), yielding occurs with a delay (as predicted by viscoplastic models). Scaling analysis suggests that this delay varies logarithmically with the yield stress and is inversely proportional to the elastic modulus. The intensity of the initial pre-yield motion increases with higher yield stress and lower elastic modulus. The observed dynamics resemble those of under- and partially over-damped systems, with a power-law fit providing an excellent match for the variation of oscillation frequency with the elastic modulus.

The transport of particles in elastoviscoplastic (EVP) fluids is of significant interest across various industrial and scientific domains. However, the physical mechanisms underlying the various particle distribution patterns observed in experimental studies remain inadequately understood in the current literature. To bridge this gap, we perform interface-resolved direct numerical simulations to study the collective dynamics of spherical particles suspended in a pressure-driven EVP duct flow. In particular, we investigate the effects of solid volume fraction, yield stress, inertia, elasticity, shear-Thinning viscosity, and secondary flows on particle migration and formation of plug regions in the suspending fluid. Various cross-streamline migration patterns are observed depending on the rheological parameters of the carrier fluid. In EVP fluids with constant plastic viscosity, particles aggregate into a large cluster at the duct centre. Conversely, EVP fluids with shear-Thinning plastic viscosity induce particle migration towards the duct walls, leading to formation of particle trains at the corners. Notably, we observe significant secondary flows (compared to the mean velocity) in shear-Thinning EVP suspensions, arising from the interplay of elasticity, shear-Thinning viscosity and particle presence, which further enhances corner-ward particle migration. We elucidate the physical mechanism by which yield stress augments the first normal stress difference, thereby significantly amplifying elastic effects. Furthermore, through a comprehensive analysis of various EVP suspensions, we identify critical thresholds for elasticity and yield stress necessary to achieve particle focusing at the duct corners.

Optical coherence tomography (OCT) has become an indispensable tool for investigating mesoscopic features in soft matter and fluid mechanics. Its ability to provide high-resolution, non-invasive measurements in both spatial and temporal domains bridges critical gaps in experimental instrumentation, enabling the study of complex, confined, and dynamic systems. This review serves as both an introduction to OCT and a practical guide for researchers seeking to adopt this technology. A set of tutorials, complemented by Python scripts, is provided for both intensity- and Doppler-based techniques. The versatility of OCT is illustrated through case studies, including time-resolved velocimetry, particle-based velocity measurements, slip velocity characterization, detection of shear-induced structures, and analysis of fluid-fluid and fluid-structure interactions. Drawing on our experiences, we also present a set of practical guidelines for avoiding common pitfalls.

We introduce a novel unsteady shear protocol, which we name rotary shear (RS), where the flow and vorticity directions are continuously rotated around the velocity-gradient direction by imposing two out-of-phase oscillatory shears (OSs) in orthogonal directions. We perform numerical simulations of dense suspensions of rigid non-Brownian spherical particles at volume fractions between 0.40 and 0.55, subject to this new RS protocol, and compare with the classical OS protocol. We find that the suspension viscosity displays a similar non-monotonic response as the strain amplitude is increased: a minimum viscosity is found at an intermediate, volume-fraction-dependent strain amplitude. However, the suspension dynamics is different in the new protocol. Unlike the OS protocol, suspensions under RS do not show absorbing states at any and do not undergo the reversible-irreversible transition: the stroboscopic particle dynamics is always diffusive, which we attribute to the fact that the RS protocol is inherently irreversible due to its design. To validate this hypothesis, we introduce a reversible-RS (RRS) protocol, a combination of RS and OS, where we rotate the shear direction (as in RS) until it is instantaneously reversed (as in OS), and find the resulting rheology and dynamics to be closer to OS. Detailed microstructure analysis shows that both the OS and RRS protocols result in a contact-free, isotropic to an in-contact, anisotropic microstructure at the dynamically reversible-to-irreversible transition. The RS protocol does not render such a transition, and the dynamics remains diffusive with an in-contact, anisotropic microstructure for all strain amplitudes.

Nanoscale biological particles, such as lipoproteins (10–80 nm) or extracellular vesicles (30–200 nm), play pivotal roles in health and disease, including conditions like cardiovascular disorders and cancer. Their effective analysis is crucial for applications in diagnostics, quality control, and nanomedicine development. While elasto-inertial focusing offers a powerful method to manipulate particles without external fields, achieving consistent focusing of nanoparticles (<500 nm) has remained a challenge. In this study, elasto-inertial focusing of nanoparticles as small as 25 nm is experimentally demonstrated using straight high-aspect-ratio microchannels in a sheathless flow. Systematic investigations reveal the influence of channel width, particle size, viscoelastic concentration, and flow rate on focusing behavior. Additionally, through numerical simulations and experimental validation, insights are provided into particle migration dynamics and viscoelastic forces governing nanoparticle focusing. Finally, biological particles, including liposomes (90–140 nm), extracellular vesicles (100 nm), and lipoproteins (10–25 nm) is successfully focused, under optimized conditions, showcasing potential applications in medical diagnostics and targeted drug delivery. These findings mark a significant advancement toward size-based high-resolution particle separation, with implications for biomedicine and environmental sciences.

We study the role of the capillary number, Ca and of the surface wettability on the dynamics of the interface between an invading and a defending phase in a porous medium by means of numerical simulations. We employ a hybrid phase field-immersed boundary approach to successfully model the contact line dynamics over the solid objects. Using a phase-field method which naturally incorporates dynamic wetting we eliminate the need for empirical contact line models to address contact line singularity. We map the two dominant modes governing the motion of the interface, namely, capillary fingering, and stable penetration, in the (Ca - theta) plane, with theta the static contact angle prescribed at the solid pores. Capillary fingering dominates at lower values of Ca and pores hydrophobic to the invading phase, while a stable penetration is observed on hydrophillic surfaces. We present new measurements and analyses, including curvature probability density functions and average curvature. We also show that the pressure needed for the invading phase to advance at constant flow rate decreases with the capillary number, and increases with the contact angle at the capillary numbers considered. The latter is due to a significant increase in the length of the interface in the case of capillary fingering. Finally, we show that it is possible to identify the different interfacial modes by measuring the penetration length and velocity during the medium filling.

In this contribution, we present a robust and efficient computational framework capable of accurately capturing the dynamic motion and large deformation/deflection responses of highly-flexible rods interacting with an incompressible viscous flow. Within the partitioned approach, we adopt separate field solvers to compute the dynamics of the immersed structures and the evolution of the flow field over time, considering finite Reynolds numbers. We employ a geometrically exact, nonlinear Cosserat rod formulation in the context of the isogeometric analysis (IGA) technique to model the elastic responses of each rod in three dimensions (3D). The Navier–Stokes equations are resolved using a pressure projection method on a standard staggered Cartesian grid. The direct-forcing immersed boundary method is utilized for coupling the IGA-based structural solver with the finite-difference fluid solver. In order to fully exploit the accuracy of the IGA technique for FSI simulations, the proposed framework introduces a new procedure that decouples the resolution of the structural domain from the fluid grid. Uniformly distributed Lagrangian markers with density relative to the Eulerian grid are generated to communicate between Lagrangian and Eulerian grids consistently with IGA. We successfully validate the proposed computational framework against two- and three-dimensional FSI benchmarks involving flexible filaments undergoing large deflections/motions in an incompressible flow. We show that six times coarser structural mesh than the flow Eulerian grid delivers accurate results for classic benchmarks, leading to a major gain in computational efficiency. The simultaneous spatial and temporal convergence studies demonstrate the consistent performance of the proposed framework, showing that it conserves the order of the convergence, which is the same as that of the fluid solver.

A gas bubble sitting at a liquid-gas interface can burst following the rupture of the thin liquid film separating it from the ambient, owing to the large surface energy of the resultant cavity. This bursting bubble forms capillary waves, a Worthington jet and subsequent droplets for a Newtonian liquid medium. However, rheological properties of the liquid medium like elastoviscoplasticity can greatly affect these dynamics. Using direct numerical simulations, this study exemplifies how the complex interplay between elasticity (in terms of elastic stress relaxation) and yield stress influences the transient interfacial phenomenon of bursting bubbles. We investigate how bursting dynamics depends on capillary, elastic and yield stresses by exploring the parameter space of the Deborah number ${{\textit {De}}}$ (dimensionless relaxation time of elastic stresses) and the plastocapillary number $\mathcal {J}$ (dimensionless yield-stress of the medium), delineating four distinct characteristic behaviours. Overall, we observe a non-monotonic effect of elastic stress relaxation on the jet development while plasticity of the elastoviscoplastic (EVP) medium is shown to affect primarily the jet evolution only at faster relaxation times (low ${{\textit {De}}}$). The role of elastic stresses on jet development is elucidated with the support of energy budgets identifying different modes of energy transfer within the EVP medium. The effects of elasticity on the initial progression of capillary waves and droplet formation are also studied. In passing, we study the effects of solvent-polymer viscosity ratio on bursting dynamics and show that polymer viscosity can increase the jet thickness apart from reducing the maximum height of the jet.

The combination of flow elasticity and inertia has emerged as a viable tool for focusing and manipulating particles using microfluidics. Although there is considerable interest in the field of elasto-inertial microfluidics owing to its potential applications, research on particle focusing has been mostly limited to low Reynolds numbers (Re<1), and particle migration toward equilibrium positions has not been extensively examined. In this work, we thoroughly studied particle focusing on the dynamic range of flow rates and particle migration using straight microchannels with a single inlet high aspect ratio. We initially explored several parameters that had an impact on particle focusing, such as the particle size, channel dimensions, concentration of viscoelastic fluid, and flow rate. Our experimental work covered a wide range of dimensionless numbers (0.05 < Reynolds number < 85, 1.5 < Weissenberg number < 3800, 5 < Elasticity number < 470) using 3, 5, 7, and 10 µm particles. Our results showed that the particle size played a dominant role, and by tuning the parameters, particle focusing could be achieved at Reynolds numbers ranging from 0.2 (1 µL/min) to 85 (250 µL/min). Furthermore, we numerically and experimentally studied particle migration and reported differential particle migration for high-resolution separations of 5 µm, 7 µm and 10 µm particles in a sheathless flow at a throughput of 150 µL/min. Our work elucidates the complex particle transport in elasto-inertial flows and has great potential for the development of high-throughput and high-resolution particle separation for biomedical and environmental applications. (Figure presented.)

Abstract: A numerical study of yield-stress fluids flowing in porous media is presented. The porous media is randomly constructed by non-overlapping mono-dispersed circular obstacles. Two class of rheological models are investigated: elastoviscoplastic fluids (i.e. Saramito model) and viscoplastic fluids (i.e. Bingham model). A wide range of practical Weissenberg and Bingham numbers is studied at three different levels of porosities of the media. The emphasis is on revealing some physical transport mechanisms of yield-stress fluids in porous media when the elastic behaviour of this kind of fluids is incorporated. Thus, computations of elastoviscoplastic fluids are performed and are compared with the viscoplastic fluid flow properties. At a constant Weissenberg number, the pressure drop increases both with the Bingham number and the solid volume fraction of obstacles. However, the effect of elasticity is less trivial. At low Bingham numbers, the pressure drop of an elastoviscoplastic fluid increases compared to a viscoplastic fluid, while at high Bingham numbers we observe drag reduction by elasticity. At the yield limit (i.e. infinitely large Bingham numbers), elasticity of the fluid systematically promotes yielding: elastic stresses help the fluid to overcome the yield stress resistance at smaller pressure gradients. We observe that elastic effects increase with both Weissenberg and Bingham numbers. In both cases, elastic effects finally make the elastoviscoplastic flow unsteady, which consequently can result in chaos and turbulence. Graphical abstract: (Figure presented.).