The interplay between viscoelasticity and inertia in dilute polymer solutions at high deformation rates can result in inertioelastic instabilities. The nonlinear evolution of these instabilities generates a state of turbulence with significantly different spatiotemporal features compared to Newtonian turbulence, termed elastoinertial turbulence (EIT). We ex- plore EIT by studying the dynamics of a submerged planar jet of a dilute aqueous polymer solution injected into a quiescent tank of water using a combination of schlieren imaging and laser Doppler velocimetry (LDV). We show how fluid elasticity has a nonmonotonic effect on the jet stability depending on its magnitude, creating two distinct regimes in which elastic effects can either destabilize or stabilize the jet. In agreement with linear stability analyses of viscoelastic jets, an inertioelastic shear-layer instability emerges near the edge of the jet for small levels of elasticity, independent of bulk undulations in the fluid column. The growth of this disturbance mode destabilizes the flow, resulting in a turbulence transition at lower Reynolds numbers and closer to the nozzle compared to the conditions required for the transition to turbulence in a Newtonian jet. Increasing the fluid elasticity merges the shear-layer instability into a bulk instability of the jet column. In this regime, elastic tensile stresses generated in the shear layer act as an “elastic membrane” that partially stabilizes the flow, retarding the transition to turbulence to higher levels of inertia and greater distances from the nozzle. In the fully turbulent state far from the nozzle, planar viscoelastic jets exhibit unique spatiotemporal features associated with EIT. The time-averaged angle of jet spreading, an Eulerian measure of the degree of entrainment, and the centerline velocity of the jets both evolve self-similarly with distance from the nozzle. The autocovariance of the schlieren images in the fully turbulent region of the jets shows coherent structures that are elongated in the streamwise direction, consistent with the suppression of streamwise vortices by elastic stresses. These coherent structures give a higher spectral energy to small frequency modes in EIT characterized by LDV measurements of the velocity fluctuations at the jet centerline. Finally, our LDV measurements reveal a frequency spectrum characterized by a −3 power-law exponent, different from the well-known −5/3 power-law exponent characteristic of Newtonian turbulence.
more »
« less
This content will become publicly available on December 1, 2024
Seagrass deformation affects fluid instability and tracer exchange in canopy flow
Abstract Monami is the synchronous waving of a submerged seagrass bed in response to unidirectional fluid flow. Here we develop a multiphase model for the dynamical instabilities and flow-driven collective motions of buoyant, deformable seagrass. We show that the impedance to flow due to the seagrass results in an unstable velocity shear layer at the canopy interface, leading to a periodic array of vortices that propagate downstream. Our simplified model, configured for unidirectional flow in a channel, provides a better understanding of the interaction between these vortices and the seagrass bed. Each passing vortex locally weakens the along-stream velocity at the canopy top, reducing the drag and allowing the deformed grass to straighten up just beneath it. This causes the grass to oscillate periodically even in the absence of water waves. Crucially, the maximal grass deflection is out of phase with the vortices. A phase diagram for the onset of instability shows its dependence on the fluid Reynolds number and an effective buoyancy parameter. Less buoyant grass is more easily deformed by the flow and forms a weaker shear layer, with smaller vortices and less material exchange across the canopy top. While higher Reynolds number leads to stronger vortices and larger waving amplitudes of the seagrass, waving amplitude is maximized at intermediate grass buoyancy. All together, our theory and computations develop an updated schematic of the instability mechanism consistent with experimental observations.
more »
« less
- Award ID(s):
- 1756279
- NSF-PAR ID:
- 10450541
- Date Published:
- Journal Name:
- Scientific Reports
- Volume:
- 13
- Issue:
- 1
- ISSN:
- 2045-2322
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Direct numerical simulations are performed to investigate a stratified shear layer at high Reynolds number ( $Re$ ) in a study where the Richardson number ( $Ri$ ) is varied among cases. Unlike previous work on a two-layer configuration in which the shear layer resides between two layers with constant density, an unbounded fluid with uniform stratification is considered here. The evolution of the shear layer includes a primary Kelvin–Helmholtz shear instability followed by a wide range of secondary shear and convective instabilities, similar to the two-layer configuration. During transition to turbulence, the shear layers at low $Ri$ exhibit a period of thickness contraction (not observed at lower $Re$ ) when the momentum and buoyancy fluxes are counter-gradient. The behaviour in the turbulent regime is significantly different from the case with a two-layer density profile. The transition layers, which are zones with elevated shear and stratification that form at the shear-layer edges, are stronger and also able to support a significant internal wave flux. After the shear layer becomes turbulent, mixing in the transition layers is shown to be more efficient than that which develops in the centre of the shear layer. Overall, the cumulative mixing efficiency ( $E^C$ ) is larger than the often assumed value of 1/6. Also, $E^C$ is found to be smaller than that in the two-layer configuration at moderate Ri . It is relatively less sensitive to background stratification, exhibiting little variation for $0.08 \leqslant Ri \leqslant 0.2$ . The dependence of mixing efficiency on buoyancy Reynolds number during the turbulence phase is qualitatively similar to homogeneous sheared turbulence.more » « less
-
more » « less
Abstract A new modeling methodology for ripple dynamics driven by oscillatory flows using a Eulerian two‐phase flow approach is presented in order to bridge the research gap between near‐bed sediment transport via ripple migration and suspended load transport dictated by ripple induced vortices. Reynolds‐averaged Eulerian two‐phase equations for fluid phase and sediment phase are solved in a two‐dimensional vertical domain with a
k ‐ε closure for flow turbulence and particle stresses closures for short‐lived collision and enduring contact. The model can resolve full profiles of sediment transport without making conventional near‐bed load and suspended load assumptions. The model is validated with an oscillating tunnel experiment of orbital ripple driven by a Stokes second‐order (onshore velocity skewed) oscillatory flow with a good agreement in the flow velocity and sediment concentration. Although the suspended sediment concentration far from the ripple in the dilute region was underpredicted by the present model, the model predicts an onshore ripple migration rate that is in very good agreement with the measured value. Another orbital ripple case driven by symmetric sinusoidal oscillatory flow is also conducted to contrast the effect of velocity skewness. The model is able to capture a net offshore‐directed suspended load transport flux due to the asymmetric primary vortex consistent with laboratory observation. More importantly, the model can resolve the asymmetry of onshore‐directed near‐bed sediment flux associated with more intense boundary layer flow speed‐up during onshore flow cycle and sediment avalanching near the lee ripple flank which force the onshore ripple migration. -
This study discusses turbulent suspension flows of non-Brownian, non-colloidal, neutrally buoyant and rigid spherical particles in a Newtonian fluid over porous media with particles too large to penetrate and move through the porous layer. We consider suspension flows with the solid volume fraction ${{\varPhi _b}}$ ranging from 0 to 0.2, and different wall permeabilities, while porosity is constant at 0.6. Direct numerical simulations with an immersed boundary method are employed to resolve the particles and flow phase, with the volume-averaged Navier–Stokes equations modelling the flow within the porous layer. The results show that in the presence of particles in the free-flow region, the mean velocity and the concentration profiles are altered with increasing porous layer permeability because of the variations in the slip velocity and wall-normal fluctuations at the suspension-porous interface. Furthermore, we show that variations in the stress condition at the interface significantly affect the particle near-wall dynamics and migration toward the channel core, thereby inducing large modulations of the overall flow drag. At the highest volume fraction investigated here, ${{\varPhi _b}}= 0.2$ , the velocity fluctuations and the Reynolds shear stress are found to decrease, and the overall drag increases due to the increase in the particle-induced stresses.more » « less
-
null (Ed.)We investigate the self-propulsion of an inertial swimmer in a linearly density stratified fluid using the archetypal squirmer model which self-propels by generating tangential surface waves. We quantify swimming speeds for pushers (propelled from the rear) and pullers (propelled from the front) by direct numerical solution of the Navier–Stokes equations using the finite volume method for solving the fluid flow and the distributed Lagrange multiplier method for modelling the swimmer. The simulations are performed for Reynolds numbers ( $Re$ ) between 5 and 100 and Froude numbers ( $Fr$ ) between 1 and 10. We find that increasing the fluid stratification strength reduces the swimming speeds of both pushers and pullers relative to their speeds in a homogeneous fluid. The increase in the buoyancy force experienced by these squirmers due to the trapping of lighter fluid in their respective recirculatory regions as they move in the heavier fluid is one of the reasons for this reduction. With increasing the stratification, the isopycnals tend to deform less, which offers resistance to the flow generated by the squirmers around them to propel themselves. This resistance increases with stratification, thus, reducing the squirmer swimming velocity. Stratification also stabilizes the flow around a puller keeping it axisymmetric even at high $Re$ , thus, leading to stability which is otherwise absent in a homogeneous fluid for $Re$ greater than $O(10)$ . On the contrary, a strong stratification leads to instability in the motion of pushers by making the flow around them unsteady and three-dimensional, which is otherwise steady and axisymmetric in a homogeneous fluid. A pusher is a more efficient swimmer than a puller owing to efficient convection of vorticity along its surface and downstream. Data for the mixing efficiency generated by individual squirmers explain the trends observed in the mixing produced by a swarm of squirmers.more » « less
