Adding polymers to turbulent Newtonian fluid flows can have dramatic effects. A well-known example is a significant drag reduction by flexible polymers in turbulent wall-bounded flows. In numerical studies of polymer drag reduction, polymers are often modelled as dumbbells of two beads connected by a finitely extensible nonlinear elastic (FENE) spring. There are natural queries whether this highly simplified coarse-grained model is adequate for describing a polymer macromolecule in turbulent flows. By carrying out Eulerian–Lagrangian simulations of polymers, described by different models, in a turbulent pipe flow, Serafini et al. (J. Fluid Mech., vol. 987, 2024, R1) have demonstrated that the FENE dumbbell model can accurately capture polymer extension statistics as compared with the realistic Kuhn chain model. Their work further reveals the surprising result that increasing the number of beads in a FENE chain worsens its accuracy in characterizing polymer spatial conformations at large Weissenberg numbers.

Motivated by the recent numerical results of Khalid et al. (Phys. Rev. Lett., vol. 127, 2021, 134502), we consider the large-Weissenberg-number ($W$) asymptotics of the centre mode instability in inertialess viscoelastic channel flow. The instability is of the critical layer type in the distinguished ultra-dilute limit where $W(1-\beta )=O(1)$ as $W \rightarrow \infty$ ($\beta$ is the ratio of solvent-to-total viscosity). In contrast to centre modes in the Orr–Sommerfeld equation, $1-c=O(1)$ as $W \rightarrow \infty$, where $c$ is the phase speed normalised by the centreline speed as a central ‘outer’ region is always needed to adjust the non-zero cross-stream velocity at the critical layer down to zero at the centreline. The critical layer acts as a pair of intense ‘bellows’ which blows the flow streamlines apart locally and then sucks them back together again. This compression/rarefaction amplifies the streamwise-normal polymer stress which in turn drives the streamwise flow through local polymer stresses at the critical layer. The streamwise flow energises the cross-stream flow via continuity which in turn intensifies the critical layer to close the cycle. We also treat the large-Reynolds-number ($Re$) asymptotic structure of the upper (where $1-c=O(Re^{-2/3})$) and lower branches of the $Re$–$W$ neutral curve, confirming the inferred scalings from previous numerical computations. Finally, we remark that the viscoelastic centre-mode instability was actually first observed in viscoelastic Kolmogorov flow by Boffetta et al. (J. Fluid Mech., vol. 523, 2005, pp. 161–170).

Polymer chains in turbulent flows are generally modelled as dumbbells, i.e. two beads joined by a nonlinear spring. The dumbbell only maps a single spatial configuration, described by the polymer end-to-end vector, thus a multi-bead FENE (finitely extensible nonlinear elastic) chain seems a natural improvement for a more accurate characterisation of the polymer spatial conformation. At a large Weissenberg number, a comparison with the more accurate Kuhn chain reveals that the multi-bead FENE chain drastically overestimates the probability of folded configurations. Surprisingly, the dumbbell turns out to be the only meaningful bead-spring model to coarse-grain a polymer macromolecule in turbulent pipe flows.

We combine flow visualisation techniques and particle image velocimetry to experimentally investigate the higher-order transition to elastoinertial turbulence of Boger fluids ($El = 0.11\unicode{x2013}0.34$) in Taylor–Couette flows. The observed route to turbulence is associated with the appearance of chaotic inflow jets, termed flame patterns, for increasing inertia, and stable structures of solitons, known as diwhirls, for decreasing inertia. We also report for the first time spatially and temporally resolved flow fields in the meridional plane for the three characteristic viscoelastic flow regimes (diwhirls, flame patterns and elastoinertial turbulence). We observe in all cases coherent structures of dynamically independent solitary vortex pairs. The stability of these coherent structures is jet-dominated and can be mainly ascribed to the high extension of the polymer chains in the inflow boundaries in the $r$–$z$ plane. Solitary pairs are self-sustained when created through random events and do not split; instead, they merge when moving sufficiently close and annihilate when hoop stresses are not sufficient to sustain them. The highly localised and random events result in highly fluctuating, chaotic flow states. We estimate the decay exponent of spatial power spectral density, illustrating a universal scaling of $-2.5$ for elastoinertial turbulence. Based on our observations and in an effort to unify and combine precedent theories with our results, we suggest a mechanism for the origins of elastoinertial instabilities, accounting for both the effect of elasticity on the vortex formation and the effect of increasing/decreasing inertia on the flow dynamics.

It is well known that inertia-free shearing flows of a viscoelastic fluid with curved streamlines, such as the torsional flow between a rotating cone and plate or the flow in a Taylor–Couette geometry, can become unstable to a three-dimensional time-dependent instability at conditions exceeding a critical Weissenberg ($Wi$) number. However, the combined effects of fluid elasticity, shear thinning and finite inertia (as quantified by the Reynolds number $Re$) on the onset of elasto-inertial instabilities are not fully understood. Using a set of cone–plate geometries, we experimentally explore the entire $Wi$–$Re$ phase space for a series of nonlinear viscoelastic fluids (with the dependence on shear rate $\dot{\gamma}$ quantified using a shear-thinning parameter $\beta _P(\dot {\gamma })$). We tune $\beta _P(\dot {\gamma })$ by varying the dissolved polymer concentration in solution. This progressively reduces shear thinning but leads to finite inertial effects before the onset of elastic instability, and thus naturally results in elasto-inertial coupling. Time-resolved rheometric measurements and flow visualization experiments allow us to investigate the effects of flow geometry, and document the combined effects of varying $Wi, Re$ and $\beta _P(\dot {\gamma })$ on the emergence of secondary motions at the onset of instability. The resulting critical state diagram quantitatively depicts the competition between the stabilizing effects of shear thinning and the destabilizing effects of inertia. We extend the curved streamline instability criterion of Pakdel & McKinley (Phys. Rev. Lett., vol. 77, no. 12, 1996, p. 2459) for the onset of purely elastic instability in curvilinear geometries by using scaling arguments to incorporate shear thinning and finite inertial effects. The augmented condition facilitates predictions of the onset of instability over a broader range of flow conditions, and thus bridges the gap between purely elastic and elasto-inertial curved streamline instabilities.

A 4-component clay-polymer-salt-water system was studied by neutron scattering. The clay-salt-water system consisted of n-butylammonium vermiculite, n-butylammonium chloride and heavy water, and the volume fraction of clay in the system was held constant, at r = 0.01. Three polymers in the molecular weight range 10,000 to 30,000 were studied, poly(vinyl methyl ether) (PVME), poly (ethylene oxide) (PEO) and poly(acrylic acid) (PAA), at a polymer volume fraction of v = 0.01. The addition of PAA suppressed the clay swelling, irrespective of the salt concentration, c. The addition of the neutral polymers had no effect on the phase transition temperature, Tc, between the gel and tactoid phases of the system, its value remaining at 14 °C for c = 0.1 M and 30 °C for c = 0.01 M. At c = 0.01 M, the neutral polymers also had a negligible effect on the lattice constant d along the swelling axis of the clay colloid, but at c = 0.1 M, the d-value was significantly lower than in the system without added polymer. For a PVME sample of molecular weight 18,000, both d and Tc were measured as a function of ν, for volume fractions between 0 and 0.04. The addition of polymer, up to v = 0.04, had no effect on Tc. However, even for v values as low as 0.001, the vermiculite layers in the gel phase were more parallel and more regularly spaced than in the system without added polymer. In the gel phase, d decreased exponentially as a function of v, from 12 nm at v = 0 to 8 nm at v = 0.04. In the tactoid phase, at T < 14 °C, the d-value in the crystalline regions was equal to 1.94 nm at v = 0 and v = 0.04, showing that the spacing between the vermiculite layers is not affected by the added polymer when they are collapsed by an increase in temperature. The addition of a PVME sample of molecular weight 110,000, at v = 0.001, had no noticeable effect on either d or Tc.

The motion of a freely rotating prolate spheroid in a simple shear flow of a dilute polymeric solution is examined in the limit of large particle aspect ratio, $\kappa$. A regular perturbation expansion in the polymer concentration, $c$, a generalized reciprocal theorem, and slender body theory to represent the velocity field of a Newtonian fluid around the spheroid are used to obtain the $O(c)$ correction to the particle's orientational dynamics. The resulting dynamical system predicts a range of orientational behaviours qualitatively dependent upon $c\, De$ ($De$ is the imposed shear rate times the polymer relaxation time) and $\kappa$ and quantitatively on $c$. At a small but finite $c\, De$, the particle spirals towards a limit cycle near the vorticity axis for all initial conditions. Upon increasing $\kappa$, the limit cycle becomes smaller. Thus, ultimately the particle undergoes a periodic motion around and at a small angle from the vorticity axis. At moderate $c\, De$, a particle starting near the flow–gradient plane departs it monotonically instead of spirally, as this plane (a limit cycle at smaller $c\, De$) obtains a saddle and an unstable node. The former is close to the flow direction. Upon further increasing $c\, De$, the saddle node changes to a stable node. Therefore, depending upon the initial condition, a particle may either approach a periodic orbit near the vorticity axis or obtain a stable orientation near the flow direction. Upon further increasing $c\, De$, the limit cycle near the vorticity axis vanishes, and the particle aligns with the flow direction for all starting orientations.

Polymers in a turbulent flow are stretched out by the fluctuating velocity gradient and exhibit a broad distribution of extensions $R$; the stationary probability density function (p.d.f.) of $R$ has a power-law tail with an exponent that increases with the Weissenberg number $\mathit {Wi}$, a non-dimensional measure of polymer elasticity. This study addresses the following questions. (i) What is the role of the non-Gaussian statistics of the turbulent velocity gradient on polymer stretching? (ii) How does the p.d.f. of $R$ evolve to its asymptotic stationary form? Our analysis is based on simulations of the dynamics of finitely extensible bead–spring dumbbells and chains, in the extremely dilute limit, that are transported in a homogeneous and isotropic turbulent flow, as well as in a Gaussian random flow. We show that while the turbulent flow is more effective at stretching small-$\mathit {Wi}$ stiff polymers, the Gaussian flow is more effective for high-$\mathit {Wi}$ polymers. This suggests that high-$\mathit {Wi}$ polymers (with large relaxation times) are stretched primarily by the cumulative effect of moderate strain rate events, rather than by short-lived extreme-valued strain rates. Next, we show that, beginning from a distribution of coiled polymers, the p.d.f. of $R$ exhibits two distinct regimes of evolution. At low to moderate $\mathit {Wi}$, the p.d.f. quickly develops a power-law tail with an exponent that evolves in time and approaches its stationary value exponentially. At high $\mathit {Wi}$, the rapid stretching of polymers first produces a peak in the p.d.f. near their maximum extension; a power law with a constant exponent then emerges and expands its range towards smaller $R$. The time scales of equilibration, measured as a function of $\mathit {Wi}$, point to a critical slowing down at the coil–stretch transition. Importantly, these results show no qualitative change when chains in a turbulent flow are replaced by dumbbells in a Gaussian flow, thereby supporting the use of the latter for reduced-order modelling.

Many key environmental, industrial and energy processes rely on controlling fluid transport within subsurface porous media. These media are typically structurally heterogeneous, often with vertically layered strata of distinct permeabilities – leading to uneven partitioning of flow across strata, which can be undesirable. Here, using direct in situ visualization, we demonstrate that polymer additives can homogenize this flow by inducing a purely elastic flow instability that generates random spatio-temporal fluctuations and excess flow resistance in individual strata. In particular, we find that this instability arises at smaller imposed flow rates in higher-permeability strata, diverting flow towards lower-permeability strata and helping to homogenize the flow. Guided by the experiments, we develop a parallel-resistor model that quantitatively predicts the flow rate at which this homogenization is optimized for a given stratified medium. Thus, our work provides a new approach to homogenizing fluid and passive scalar transport in heterogeneous porous media.

We consider the multidirectional swelling and drying of hydrogels formed from super-absorbent polymers and water, focusing on the elastic deformation caused by differential swelling. By modelling hydrogels as instantaneously incompressible, linear-elastic materials and considering situations in which there can be large isotropic strains (arising from swelling) while deviatoric strains remain small, it is possible to describe accurately a wider range of gel states than traditional linear elastic theories allow. An equation is derived relating the displacement field to the polymer fraction in such hydrogels, permitting the shape of the swelling gel to be determined as it evolves in time, using the formulation of Part 1 to find the local polymer fraction. We discuss the boundary conditions to be applied at the surfaces of a gel, both on the bulk elastic stress and on the pervadic (pore) pressure in the interstices. Similarities between the equation for the displacements and the equations of classical plate theory are investigated by considering a model problem of a slender cylinder with its base immersed in water drying by evaporation into the surrounding air. In this problem, there is differential drying along the axis of the cylinder, as the base remains swollen while the top dries. The results of our displacement formulation agree qualitatively with experiments that we have conducted, and provide a physical interpretation of the forced biharmonic equation describing the displacement field.

We introduce a new approach for modelling the swelling, drying and elastic behaviour of hydrogels, which leverages the tractability of classical linear-elastic theory whilst incorporating nonlinearities arising from large swelling strains. Relative to a reference state of a fully swollen gel, in which the polymer scaffold may only comprise less than $1\,\%$ of the total volume, a constitutive model for the Cauchy stress tensor is presented, which linearises around small deviatoric strains corresponding to shearing deformations of the material whilst allowing for a nonlinear relation between stress and isotropic strains. Such isotropic strains are considered only to be a consequence of losses and gains of water, while the hydrogel is taken to be instantaneously incompressible. The dynamics governing swelling and drying are described by coupling the interstitial flow of the water through the porous gel with the elastic response of the gel. This approach allows for a complete description of gel behaviour using only three macroscopic polymer-fraction-dependent parameters: an osmotic modulus, a shear modulus and a permeability. It is shown how these three material parameters can, in principle, be determined experimentally using a simple rheometry experiment in which a gel is compressed between two plates surrounded by water and the total force on the top plate is measured. To illustrate our approach, we solve for the swelling of a gel under horizontal confinement and for a partially dried hydrogel bead placed in water.

Recent studies of pressure-driven flows of dilute polymer solutions in straight channels demonstrated the existence of two-dimensional coherent structures that are disconnected from the laminar state and appear through a subcritical bifurcation from infinity. These travelling-wave solutions were suggested to organise the phase-space dynamics of purely elastic and elasto-inertial chaotic channel flows. Here, we consider a wide range of parameters, covering the purely elastic and elasto-inertial cases, and demonstrate that the two-dimensional travelling-wave solutions are unstable when embedded in sufficiently wide three-dimensional domains. Our work demonstrates that studies of purely elastic and elasto-inertial turbulence in straight channels require three-dimensional simulations, and no reliable conclusions can be drawn from studying strictly two-dimensional channel flows.

Thermo-responsive hydrogels are a promising material for creating controllable actuators for use in micro-scale devices, since they expand and contract significantly (absorbing or expelling fluid) in response to relatively small temperature changes. Understanding such systems can be difficult because of the spatially and temporally varying properties of the gel, and the complex relationships between the fluid dynamics, elastic deformation of the gel and chemical interaction between the polymer and fluid. We address this using a poro-elastic model, considering the dynamics of a thermo-responsive spherical hydrogel after a sudden change in the temperature that should result in substantial swelling or shrinking. We focus on two model examples, with equilibrium parameters extracted from data in the literature. We find a range of qualitatively different behaviours when swelling and shrinking, including cases where swelling and shrinking happen smoothly from the edge, and other situations that result in the formation of an inwards-travelling spherical front that separates a core and shell with markedly different degrees of swelling. We then characterise when each of these scenarios is expected to occur. An approximate analytical form for the front dynamics is developed, with two levels of constant porosity, that well-approximates the numerical solutions. This system can be evolved forward in time, and is simpler to solve than the full numerics, allowing for more efficient predictions to be made, such as when deciding dosing strategies for drug-laden hydrogels.

We examine the response of an inertialess viscoelastic channel flow to localized perturbations. We thus performed an experiment in which we perturbed the flow using a localized velocity pulse and probed the perturbed fluid packet downstream from the perturbation location. While for low Weissenberg numbers the perturbed fluid reaches the measurement location as a single velocity pulse, for sufficiently high Weissenberg numbers and perturbation strengths, a random number of pulses arrive at the measurement location. The average number of pulses observed increases with the Weissenberg number. This observation constitutes a transition to a novel elastic pulse-splitting regime. Our results suggest a possible new direction for studying the elastic instability of viscoelastic channel flows at high elasticity through the growth of localized perturbations.

Recently, the nature of viscoelastic drag-reducing turbulence (DRT), especially the maximum drag reduction (MDR) state, has become a focus of controversy. It has long been regarded as polymer-modulated inertial turbulence (IT), but is challenged by the newly proposed concept of elasto-inertial turbulence (EIT). This study is to repicture DRT in parallel plane channels by introducing dynamics of EIT through statistical, structural and budget analysis for a series of flow regimes from the onset of drag reduction to EIT. Some underlying mechanistic links between DRT and EIT are revealed. Energy conversion between velocity fluctuations and polymers as well as pressure redistribution effects are of particular concern, based on which a new energy self-sustaining process (SSP) of DRT is repictured. The numerical results indicate that at low Reynolds number ($Re$), weak IT flow is replaced by a laminar regime before the barrier of EIT dynamics is established with the increase of elasticity, whereas, at moderate $Re$, EIT-related SSP can get involved and survive from being relaminarized. This further explains the reason why relaminarization phenomenon is observed for low $Re$ while the flow directly enters MDR and EIT at moderate $Re$. Moreover, with the proposed energy picture, the newly discovered phenomenon that streamwise velocity fluctuations lag behind those in the wall-normal direction can be well explained. The repictured SSP certainly justifies the conjecture that IT nature is gradually replaced by that of EIT in DRT with the increase of elasticity.

This paper reports an exhaustive review of recent research progress in the development of millimeter-wave absorbers. With the advancement in technologies, microwave absorbers have shown their evolution in several applications such as defense, security identification, stealth technology, and several more. The importance of these absorbers is increasing due to electromagnetic interference (EMI) effects. Primarily, an abundant amount of absorbers has been developed in lower frequency bands and with the increasing demand for 5G technology, the usage of millimeter-wave bands has gained attention. This in turn requires the development of millimeter-wave absorbers to provide EMI shielding in several applications. Out of several materials, polymers have grabbed attention in the mitigation of EMI. An absorber that combines carbonaceous elements with polymers offers large design flexibility and tunability per the filler concentration. This paper focuses on the classification of these absorbers based on geometry as well as that by using polymers with their design challenges, merits, and demerits with their industrial applications. Comparative studies of geometrically based absorbers and polymeric-based absorbers are also shown.

It is recently found that viscoelastic pipe flows can be linearly unstable, leading to the possibility of a supercritical transition route, in contrast to Newtonian pipe flows. Such an instability is referred to as the centre mode, which was studied numerically by Chaudhary et al. (J. Fluid Mech., vol. 908, 2021, p. A11) based on an Oldroyd-B model. In this paper, we are interested in expanding the parameter space investigated and exploring the asymptotic scalings related to this centre instability in the Oldroyd-B viscoelastic pipe flow. It is found from the asymptotic analysis that the centre mode exhibits a three-layered asymptotic structure in the radial direction, a wall layer, a main layer and a central layer, which are driven by pure viscosity, axial and/or radial pressure gradient, and a combined effect of viscosity and elasticity, respectively. Depending on the relations of the control parameters, two regimes, the long-wavelength and short-wavelength centre instabilities, emerge, for which the central-layer thicknesses are of different orders of magnitude. Our large-Reynolds-number asymptotic predictions are compared to the numerical solutions of the original eigenvalue system, and favourable agreement is achieved, especially when the parameters approach their individual limits. In addition to revealing the dominant factors and their balances, the asymptotic analysis describes the instability system by reducing the number of control parameters, and furthermore explaining the collapse of the numerical results for different re-scalings.

We perform a linear stability analysis of viscoelastic plane Couette and plane Poiseuille flows with free-slip boundary conditions. The fluid is described by the Oldroyd-B constitutive model, and the flows are driven by a suitable body force. We find that both types of flow become linearly unstable, and we characterise the spatial structure of the unstable modes. By performing a boundary condition homotopy from the free-slip to no-slip boundaries, we demonstrate that the unstable modes are directly related to the least stable modes of the no-slip problem, destabilised under the free-slip situation. We discuss how our observations can be used to study recently discovered purely elastic turbulence in parallel shear flows.