This study employs a direct-forcing fictitious-domain method combined with the Oldroyd-B model to investigate the dynamics of spherical and oblate spheroidal particles in a viscoelastic shear flow. By varying the particle Reynolds number (0.1 ≤ Rep ≤ 10) and the Weissenberg number (1 ≤ Wi ≤ 5), we systematically examine the changes in particle equilibrium positions and rotational behaviours. For spherical particles, the equilibrium position shifts monotonically towards the wall with increasing Wi when Rep < 1 or Rep > 3, whereas it remains largely insensitive to elasticity for 1 ≤ Rep ≤ 3. For oblate particles, the rotational period increases steadily with Wi across low-to-moderate Rep, demonstrating that enhanced fluid elasticity significantly reduces the rotational speed. Owing to the interplay of fluid inertia and elasticity, oblate particles initialised with two representative initial orientations, denoted as np1 and np2 (specified by Euler angles as np1 = (0, 0, 0) and np2 = (π/2, π/4, π/2)) exhibit six distinct modes of lateral migration and rotational motion. Under weak fluid inertia, increasing elasticity induces a kayak-like rotation, while at moderate inertia, the motion is dominated by rolling within the shear plane. When Rep ≥ 7, the rotation may become locked, and this locked state can be released as Wi increases.

We investigate the presence of intermittency, a well-known phenomenon in Newtonian turbulence, in viscoelastic Taylor–Couette flow at fixed inertia, (Reynolds number, $\textit {Re}$ = 120). Using particle image velocimetry measurements in the r–z plane for a broad range of elasticities, $\textit{El}$ = 0–0.64, we reveal a dual role of viscoelasticity, even at moderate $\textit {Re}$: it can destabilise the flow and trigger turbulence at moderate elasticity but suppress turbulence intensity at higher values. Gradual turbulence suppression is evident through reduced velocity fluctuations, longer integral time scales proportional to the polymer relaxation time and weakened coherent structures as the flow approaches elasto-inertial turbulence. Higher-order statistics indicate that elasticity increases spatial and temporal intermittency while altering the dynamics: spatially, through a reduced viscous range and enhanced large-scale decorrelation; temporally, via a characteristic time scale consistent with polymer-driven state switching rather than a classical cascade. This spatio-temporal intermittency manifests as alternating active/hibernating turbulence cycles. During hibernation, the flow reorganises into a state distinct from active turbulence, characterised by low velocity fluctuations and Reynolds stresses, suppressing not only turbulence but also linearly unstable Newtonian states such as Taylor vortices. Rather than being featureless, this state exhibits weakened coherent structures in the bulk of the flow and sustained near-wall dynamics, where elastic stresses are likely generated. The present results demonstrate that hibernation is not exclusive to Newtonian turbulence, but also emerges in viscoelastic flows under conditions where inertia alone is insufficient to trigger or sustain turbulent motion.

Multiscale simulations offer a way to understand how polymer transformations proceed under realistic conditions. However, existing workflows rely on empirically tuned parameters. The MUltiscale SImulation frameworK (MUSIK) integrates conformer analysis, density functional theory, high-level coupled-cluster energetics, and graph-based lattice kinetic Monte Carlo (kMC) to generate chemically accurate, spatially and temporally resolved kinetic models for polymer transformations without parameter fitting. The present work introduces MUSIK as a reproducible workflow and provides a step-by-step tutorial that shows how to construct, run, and analyze a lattice kMC model for base-assisted dehydrochlorination of polyvinyl chloride (PVC), supported by post-processing scripts that automate kMC trajectory analysis and polymer-structure reconstruction. Tutorial results demonstrate how MUSIK links atomic energetics and local structure to macroscopic kinetics and structural evolution and reproduces key experimental observations for PVC dehydrochlorination. MUSIK applies broadly to arbitrary polymer chemistries and can provide practical recommendations for polymer recycling, upcycling, and materials design to address the plastic-waste challenge.

We show that vortex rings forming in suspensions and polymer solutions exhibit similar behaviour in terms of vortex-core growth, vorticity redistribution and circulation when compared at loadings (i.e. volume fraction in suspensions and concentration in polymer solutions) normalised by their respective critical loading. This suggests that the collective, loading-driven dynamics of dispersed particles or polymer chains can outweigh the unique characteristics of the dispersion and govern certain aspects of the flow response. Using a confined vortex ring as a canonical structure, we synthesise experimental data spanning dilute to concentrated regimes for both media, with our analysis indicating that similar collective dynamics exist in both systems. Specifically, the vortex-core growth rate relative to the unloaded case increases approximately linearly with normalised loading while circulation remains nearly unchanged over the measured parameter range. Shared collective dynamics are further suggested by radial vorticity profiles scaled according to loading, which reveal structural markers that align along consistent non-dimensional radii, demonstrating a topological self-similarity between flow fields despite fundamentally different disperse phases. Turbulent kinetic energy trends expose the limits of this correspondence: with increasing Reynolds number, turbulence rises in suspensions but remains largely unchanged in polymer solutions, reflecting differences in how semi-rigid particles and flexible chains interact with small-scale fluctuations. Altogether, these results uncover shared collective dynamics and support the cautious use of polymer analogues to study suspension flows within an equivalent loading regime. However, the limited test matrix motivates further experiments involving different flow media to assess the generality of these findings.

This study presents a combined analytical and experimental investigation into the rising dynamics of Newtonian droplets in viscoelastic fluids, focusing on discontinuous transitions in droplet velocity. Although velocity jumps have been observed in gas bubbles, this work is the first systematic investigation of the phenomenon in Newtonian droplets. Analytically, a triple-perturbation approach is developed using Reynolds (Re), Deborah (De) and capillary (Ca) numbers as expansion parameters. The surrounding viscoelastic medium is modelled using the nonlinear-Giesekus constitutive equation. The analysis reveals that increasing Deborah number intensifies and reorients normal stresses – particularly the polar ($\theta$-direction) and hoop components – along the droplet interface. These stresses migrate rearward with increasing elasticity, resulting in localised elongational flows and concentrated radial stresses near the rear stagnation point, which indicates the onset of elastic instabilities. Experimentally, Newtonian droplets are introduced into a variety of exterior fluids: (i) quasilinear Boger fluids, which are constant-viscosity viscoelastic solutions that isolate elastic effects; and (ii) a strongly nonlinear viscoelastic solution, which exhibits both high elasticity and significant shear-thinning behaviour. Equivalent Newtonian fluids of matched viscosity are also used as controls. Viscoelastic fluids are solutions of polyacrylamide in the water–glycerine mixture. Across all experiments, the dimensionless ranges were $\textit{Re}\approx 6\times 10^{-4}{-}4.59$, $\textit{De}\approx 0.057{-}0.802$ and $\textit{Ca}\approx 0.0026{-}0.638$, framing the scope of applicability. The experimental observations challenge the classical expectation of smooth droplet motion and demonstrate that Newtonian droplets in viscoelastic environments can exhibit sharp, weak or no velocity jumps, depending on the balance between elasticity, viscosity and deformability.

The addition of polymers to turbulent pipe flows induces significant drag reduction and fundamentally modifies turbulent flow structures. This study presents a fractal dimension analysis of polymeric turbulent pipe flows using velocity fields captured via two-dimensional particle image velocimetry in the streamwise-radial plane. Two experimental datasets were generated: one by varying the polymer concentration at a constant Reynolds number ($\textit{Re}$) and another by varying $\textit{Re}$ at a fixed polymer concentration. Friction factors were measured concurrently to quantify the extent of drag reduction. The two-dimensional fractal dimension was evaluated for isosurfaces of turbulent kinetic energy. While Newtonian turbulence exhibits a nearly constant fractal dimension at length scales exceeding a critical threshold, the introduction of polymers causes the fractal dimension to decrease monotonically with increasing concentration. Conversely, the fractal dimension remains insensitive to changes in the Reynolds number. The ratio of the critical length scale to the Kolmogorov scale varies according to both $\textit{Re}$ and polymer concentration; however, this scale ratio becomes independent of both parameters once the maximum drag reduction asymptote is reached. Spatial analysis of the one-dimensional fractal dimension across radial positions helps to further reveal the evolution of turbulence fractality. The results demonstrate that while flow inertia promotes the formation of space-filling structures, viscoelastic effects smooth these structures and transition them towards sheet-like or linear geometries. Finally, the correlation between the fractal dimension and turbulence intermittency is discussed.

Small interfering RNA (siRNA) is an emerging therapeutic modality for a varietyof diseases, including cancer, as siRNA can silence target genes in asequence-specific manner. The effective delivery of siRNA remains a majorchallenge due to rapid clearance by macrophages in the systemic environment.Nonspecific interactions with the serum proteins in the bloodstream contributeto macrophage uptake, limiting circulation time, thereby reducing the effectivedelivery of siRNA to the target site. Here, we report the efficient delivery ofsiRNA to cancer cells using hyaluronic acid (HA)-coated cationic polymericnanovector (PONI-Guan)/siRNA polyplexes. The guanidinium-functionalized polymersself-assemble with siRNA and enable cytosolic delivery. HA serves as anoninteracting protective shield on the polyplexes that prevent macrophageuptake in vitro. These nanovectors facilitate efficient siRNAdelivery to 4T1 triple-negative breast cancer cells in vitro,with a 4:1 selectivity relative to macrophages. Further, HA-coated polyplexesdemonstrated efficient STAT3 gene knockdown (~50%) in 4T1 cells. Intravenousadministration of HA-coated polyplexes in 4T1 tumor-bearing mice showedsignificantly (~50%) decreased accumulation in clearance organs, in comparisonto the PONI-Guan polyplexes. Collectively, HA-coated polyplexes provide aneffective strategy for selective siRNA delivery to tumor cells while avoidingmacrophage uptake.

To investigate how polymers influence energy transfer in three-dimensional turbulence, we conduct experiments in homogeneous bulk turbulence generated by a von Kármán swirling flow, using tomographic particle image velocimetry. A filtering approach is applied to the measured three-dimensional velocity fields to extract subgrid-scale (SGS) statistics, focusing on the filtered strain-rate tensor and SGS stress tensor. We find that polymer additives induce significant changes in the tensorial geometry: the strain-rate tensor shows a tendency towards an eigenvalue ratio of $1 : 0 : -1$, while the SGS stress tensor favours a $2 : -1 : -1$ configuration. The local energy flux – quantified by the inner product of the strain-rate and SGS stress tensors – is systematically suppressed by polymers and becomes increasingly intermittent. This suppression is linked to a reduced energy transfer efficiency, associated with the misalignment between the principal eigendirections of the two tensors. Anisotropic effects are also observed in the energy flux components, indicating that polymers affect vertical and horizontal energy transfer differently. Finally, the obtained SGS statistics allow for an a priori assessment of SGS models. Our results reveal that the nonlinear gradient model significantly outperforms the Smagorinsky model, particularly in polymer-laden turbulence. The diminished alignment between the strain-rate and SGS stress tensors may underlie the limitations of the Smagorinsky model, which assumes a scalar eddy-viscosity closure. These results provide new experimental insights into the SGS dynamics of polymeric turbulence and highlight the potential of nonlinear models for large-eddy simulations of viscoelastic flows.

We report an experimental study on the effects of polymer additives in the dissipative-scale flow field properties in turbulent Rayleigh–Bénard convection. The experiments were conducted in a cylindrical convection cell with a minute amount of polyacrylamide long-chain polymer. The local velocity gradient tensor was measured using an integrated home-made measurement system (J. Fluid Mech., 2024, vol. 984, p. A8). Although the single-roll large-scale circulation persists (owing to the slight tilt of the convection cell), polymers induce an anisotropic suppression of the dissipative-scale flow properties. The normal velocity gradient components are suppressed more than the shear components. The mean energy dissipation rate in both centre and side regions decreases, then levels off with increasing polymer concentration and the final reduction ratio exceeds 50 % in each region. In the side region, adding polymers has a stronger stabilising effect on the strain rate than the rotation. The anisotropic suppression of the velocity gradient tensor affects dissipation-rotation co-occurrence probability, velocity gradient triple decomposition and local streamline topology. Adding polymers also induces a deceleration effect and increases the contribution of local buoyancy in driving the flow. These results reveal that the addition of polymers can non-trivially manipulate dissipative-scale turbulence fields and energy cascades.

Amorphous materials transition from solid-like to liquid-like behaviour (yield) under large stresses. Their constituent elements are caged in metastable configurations due to their neighbours. Microscale interactions between these elements lead to a large energy barrier to break the cages and trigger a plastic rearrangement. Thermal fluctuations can alter the yielding point as the elements hop to new configurations in anticipation. This work bridges the gap between molecular-scale physics and bulk rheology in thermal amorphous materials by connecting a classical density functional theory to a thermally activated elastoplastic model (EPM). We use a model system of solvent-free polymer-grafted nanoparticles which show rheological characteristics similar to those of soft glassy materials. We formulate the evolution of the free energy in a deforming array of polymer-grafted nanoparticles to obtain the energy landscape as an input to our EPM. We examine how the apparent yield stress depends on the shape of the energy landscape, thermal fluctuations and the rate of deformation. Our general scaling analyses reveal different regimes of structural relaxation governed by the applied shear rate and the inherent time scale for thermal hops. The complex interplay between mechanical loading and thermal fluctuations is further characterized by performing a variety of shear tests with different deformation history. The proposed framework provides an understanding of the yielding transition by integrating across a vast range of length and time scales.

The deposition of droplets onto a swollen polymer network induces the formation of a wetting ridge at the contact line. Current models typically consider either viscoelastic effects or poroelastic effects, while polymeric gels often exhibit both properties. In this study, we investigate the growth of the wetting ridge using a comprehensive large-deformation theory that integrates both dissipative mechanisms – viscoelasticity and poroelasticity. In the purely poroelastic case, following an initial instantaneous incompressible deformation, the growth dynamics exhibits scale-free behaviour, independent of the elastocapillary length or system size. A boundary layer of solvent imbibition between the solid surface (in contact with the reservoir) and the region of minimal chemical potential is created. At later times, the ridge equilibrates on the diffusion time scale given by the elastocapillary length. When viscoelastic properties are incorporated, our findings show that, during the early stages (prior to the viscoelastic relaxation time scale), viscoelastic effects dominate the growth dynamics of the ridge and solvent transport is significantly suppressed. Beyond the relaxation time, the late-time dynamics closely resembles that of the purely poroelastic case. These findings are discussed in light of recent experiments, showing how our approach offers a new interpretation framework for wetting of polymer networks of increasing complexity.

The effectiveness of polymer drag reduction by targeted injection is studied in comparison with that of a uniform concentration (or polymer ocean) in a turbulent channel flow. Direct numerical simulations are performed using a pseudo-spectral code to solve the coupled equations of a viscoelastic fluid using the finitely extensible nonlinear elastic dumbbell model with the Peterlin approximation. Light and heavy particles are used to carry the polymer in some cases, and polymer is selectively injected into specific flow regions in the other cases. Drag reduction is computed for a polymer ocean at a viscosity ratio of $\beta = 0.9$ for simulation validation, and then various methods of polymer addition at $\beta = 0.95$ are compared for their drag-reduction performance and general effect on the flow. It was found that injecting polymer directly into regions of high axial strain inside and around coherent vortical structures was the most effective at reducing drag, while injecting polymer very close to the walls was the least effective. The targeting methods achieved up to 2.5 % higher drag reduction than an equivalent polymer ocean, offering a moderate performance boost in the low drag-reduction regime.

We investigate the scale-by-scale transfers of energy, enstrophy and helicity in homogeneous and isotropic polymeric turbulence using direct numerical simulations. The study relies on the exact scale-by-scale budget equations, derived from the governing model equations, that fully capture the back-reaction of polymers on the fluid dynamics. Polymers act as dynamic sinks and sources and open alternative routes for interscale transfer whose significance is modulated by their elasticity, quantified through the Deborah number (${\textit{De}}$). Polymers primarily deplete the nonlinear energy cascade at small scales, by attenuating intense forward and inverse transfer events. At sufficiently high ${\textit{De}}$, a polymer-driven flux emerges and dominates at small scales, transferring on average energy from larger to smaller scales, while allowing for localised backscatter. For enstrophy, polymers inhibit the stretching of vorticity, with fluid–polymer interactions becoming the primary enstrophy source at high ${\textit{De}}$. Accordingly, an analysis of the small-scale flow topology reveals that polymers promote two-dimensional straining states and enhance the occurrence of shear and planar extensional flows, while suppressing extreme rotation events. Helicity, injected at large scales, exhibits a transfer mechanism analogous to energy, being dominated by nonlinear dynamics at large scales and by polymer-induced fluxes at small scales. Polymers enhance the breakdown of small-scale mirror symmetry, as indicated by a monotonic increase in relative helicity with ${\textit{De}}$ across all scales.

In 2019, in this journal, I discussed approaches for controlling the movement of molecules, in particular macromolecules, with an emphasis on how this enabled advances in the field of drug delivery – a field that has impacted billions of people worldwide. Since 2019, there have been advances in our work and this field including a striking demonstration in which drug delivery nanoparticles were crucial to the success of mRNA therapies and the Covid-19 vaccine. In this paper, I provide updates in such areas as i) developing new methods for oral drug delivery systems, ii) delivery of molecules to specific sites of the body, iii) new types of delivery systems, and iv) examples of machine learning/artificial intelligence in these areas. I also discuss advances in mRNA technology as it relates to drug delivery and the development of nanoparticles to protect and deliver vaccines, which saved and improved the lives of hundreds of millions of people throughout the world.

The Weissenberg effect, or rod-climbing phenomenon, occurs in non-Newtonian fluids where the fluid interface ascends along a rotating rod. Despite its prominence, theoretical insights into this phenomenon remain limited. In earlier work, Joseph & Fosdick (1973, Arch. Rat. Mech. Anal. vol. 49, pp. 321–380) employed domain perturbation methods for second-order fluids to determine the equilibrium interface height by expanding solutions based on the rotation speed. In this work, we investigate the time-dependent interface height through asymptotic analysis with dimensionless variables and equations using the Giesekus model. We begin by neglecting inertia to focus on the interaction between gravity, viscoelasticity and surface tension. In the small-deformation scenario, the governing equations indicate the presence of a boundary layer in time, where the interface rises rapidly over a short time scale before gradually approaching a steady state. By employing a stretched time variable, we derive the transient velocity field and corresponding interface shape on this short time scale, and recover the steady-state shape on a longer time scale. In contrast to the work of Joseph and Fosdick, which used the method of successive approximations to determine the steady shape of the interface, we explicitly derive the interface shape for both steady and transient cases. Subsequently, we reintroduce small but finite inertial effects to investigate their interaction with viscoelasticity, and propose a criterion for determining the conditions under which rod climbing occurs. Through numerical computations, we obtain the transient interface shapes, highlighting the interplay between time-dependent viscoelastic and inertial effects.

Shear-thinning fluids flowing through pipes are crucial in many practical applications, yet many unresolved problems remain regarding their turbulent transition. Using highly robust numerical tools for the Carreau–Yasuda model, we discovered that linear instability can arise when the power-law index falls below 0.35. This inelastic non-axisymmetric instability can universally arise in generalised Newtonian fluids that extend the power-law model. The viscosity ratio from infinite to zero shear rate can significantly impact instability, even if it is small. Two branches of finite-amplitude travelling-wave solutions bifurcate subcritically from the linear critical point. The solutions exhibit sublaminar drag reduction, a phenomenon not possible in the Newtonian case.

An addition of polymers can significantly reduce drag in wall-bounded turbulent flows, such as pipes or channels. This phenomenon is accompanied by a noticeable modification of the mean-velocity profile. Starting from the premise that polymers reduce vortex stretching, we derive a theoretical prediction for the mean-velocity profile. After assessing this prediction by numerical experiments of turbulence with reduced vortex stretching, we show that the theory successfully describes experimental measurements of drag reduction in pipe flow.

The extensional rheology of dilute suspensions of spheres in viscoelastic/polymeric liquids is studied computationally. At low polymer concentration $c$ and Deborah number $\textit{De}$ (imposed extension rate times polymer relaxation time), a wake of highly stretched polymers forms downstream of the particles due to larger local velocity gradients than the imposed flow, indicated by $\Delta \textit{De}_{\textit{local}}\gt 0$. This increases the suspension’s extensional viscosity with time and $\textit{De}$ for $De \lt 0.5$. When $\textit{De}$ exceeds 0.5, the coil-stretch transition value, the fully stretched polymers from the far-field collapse in regions with $\Delta \textit{De}_{\textit{local}} \lt 0$ (lower velocity gradient) around the particle’s stagnation points, reducing suspension viscosity relative to the particle-free liquid. The interaction between local flow and polymers intensifies with increasing $c$. Highly stretched polymers impede local flow, reducing $\Delta \textit{De}_{\textit{local}}$, while $\Delta \textit{De}_{\textit{local}}$ increases in regions with collapsed polymers. Initially, increasing $c$ aligns $\Delta \textit{De}_{\textit{local}}$ and local polymer stretch with far-field values, diminishing particle–polymer interaction effects. However, beyond a certain $c$, a new mechanism emerges. At low $c$, fluid three particle radii upstream exhibits $\Delta \textit{De}_{\textit{local}} \gt 0$, stretching polymers beyond their undisturbed state. As $c$ increases, however, $\Delta \textit{De}_{\textit{local}}$ in this region becomes negative, collapsing polymers and resulting in increasingly negative stress from particle–polymer interactions at large $\textit{De}$ and time. At high $c$, this negative interaction stress scales as $c^2$, surpassing the linear increase of particle-free polymer stress, making dilute sphere concentrations more effective at reducing the viscosity of viscoelastic liquids at larger $\textit{De}$ and $c$.

Lubricant viscoelasticity arises due to a finite polymer relaxation time ($\lambda$) which can be exploited to enhance lubricant performance. In applications such as bearings, gears, biological joints, etc., where the height-to-length ratio ($H_0 / \ell _x$) is small and the shear due to the wall velocity ($U_0$) is high, a simplified two-dimensional computational analysis across the channel length and height reveals a finite increase in the load-carrying capacity of the film purely due to polymer elasticity. In channels with a finite length-to-width ratio, $a$, the spanwise effects can be significant, but the resulting mathematical model is computationally intensive. In this work, we propose simpler reduced-order models, namely via (i) a first-order perturbation in the Deborah number ($\lambda U_0 / \ell _x$) and (ii) the viscoelastic Reynolds approach extended from Ahmed & Biancofiore (J. Non-Newtonian Fluid Mech., vol. 292, 2021, 104524). We predict the variation in the net vertical force exerted on the channel walls (for a fixed film height) versus increasing viscoelasticity, modelled using the Oldroyd-B constitutive relation, and the channel aspect ratio. The models predict an increase in the net force, which is zero for the Newtonian case, versus both the Deborah number and the channel aspect ratio. Interestingly, for a fixed $\textit{De}$, this force varies strongly between the two limiting cases (i) $a \ll 1$, an infinitely wide channel, and (ii) $a \gg 1$, an infinitely short channel, implying a change in the polymer response. Furthermore, we observe a different trend (i) for a spanwise-varying channel, in which a peak is observed between the two limits, and (ii) for a spanwise-uniform channel, where the largest load value is for $a \ll 1$. When $a$ is O($1$), the viscoelastic response varies strongly and spanwise effects cannot be ignored.