Topic overview

physics.flu-dyn

2893 works7498 researchers0 institutions

Topic snapshot

What this area looks like now

2893works

7498authors

0experts visible

0communities

Next steps

Move from topic reading into action

The graph preview below keeps the nearby papers, people and communities visible in the same reading flow.

Inspect nearby papers, researchers, institutions and communities without opening a separate graph page.

Building this graph slice

BZPEER is loading the nearby papers, people, topics and institutions for this page.

preprint2016arXiv

Internal wave attractors examined using laboratory experiments and 3D numerical simulations

In the present paper, we combine numerical and experimental approaches to study the dynamics of stable and unstable internal wave attractors. The problem is considered in a classic trapezoidal setup filled with a uniformly stratified fluid. Energy is injected into the system at global scale by the small-amplitude motion of a vertical wall. Wave motion in the test tank is measured with the help of conventional synthetic schlieren and PIV techniques. The numerical setup closely reproduces the experimental one in terms of geometry and the operational range of the Reynolds and Schmidt numbers. The spectral element method is used as a numerical tool to simulate the nonlinear dynamics of a viscous salt-stratified fluid. We show that the results of three-dimensional calculations are in excellent qualitative and quantitative agreement with the experimental data, including the spatial and temporal parameters of the secondary waves produced by triadic resonance instability. Further, we explore experimentally and numerically the effect of lateral walls on secondary currents and spanwise distribution of velocity amplitudes in the wave beams. Finally, we test the assumption of a bidimensional flow and estimate the error made in synthetic schlieren measurements due to this assumption.

preprint2014arXiv

Electric-field-induced turbulent energy cascade in an oil-in-oil emulsion

We observe electro-hydrodynamically driven turbulent flows at low Reynolds numbers in a two-fluid emulsion consisting of micron-scale droplets. In the presence of electric fields, the droplets produce interacting hydrodynamic flows which result in a dynamical organization at a spatial scale much larger than the size of the individual droplets. We characterize the dynamics associated with these structures by both video imaging and a simultaneous, in situ, measurement of the time variation of the bulk Reynolds stress with a rheometer. The results display scale invariance in the energy spectra in both space and time.

preprint2006arXiv

The signature of laminar instabilities in the zone of transition to turbulence

We demonstrate that the space-time statistics of the birth of turbulent spots in boundary layers can be reconstructed qualitatively from the average behavior of macroscopic measures in the transition zone. The conclusion in \cite{vg04} that there exists a connection between the patterns in laminar instability and the birth of turbulent spots is strengthened. We examine why the relationship between instability and transition to turbulence is manifest in some cases and appears to be totally absent in others. Novel cellular automaton type simulations of the transition zone are conducted, and the pattern of spot birth is obtained from secondary instability analysis. The validity of the hypothesis of concentrated breakdown, according to which most turbulent spots originate at a particular streamwise location, is assessed. The predictions made lend themselves to straightforward experimental verification.

preprint2015arXiv

Triadic instability of a non-resonant precessing fluid cylinder

Flows forced by a precessional motion can exhibit instabilities of crucial importance, whether they concern the fuel of a flying object or the liquid core of a telluric planet. So far, stability analyses of these flows have focused on the special case of a resonant forcing. Here, we address the instability of the flow inside a precessing cylinder in the general case. We first show that the base flow forced by the cylinder precession is a superposition of a vertical or horizontal shear flow and an infinite sum of forced modes. We then perform a linear stability analysis of this base flow by considering its triadic resonance with two free Kelvin modes. Finally, we derive the amplitude equations of the free Kelvin modes and obtain an expression of the instability threshold and growth rate.

preprint2016arXiv

Numerical Study of the Sedimentation of Spheroidal Particles

The gravity-driven motion of rigid particles in a viscous fluid is relevant in many natural and industrial processes, yet this has mainly been investigated for spherical particles. We therefore consider the sedimentation of non-spherical (spheroidal) isolated and particle pairs in a viscous fluid via numerical simulations using the Immersed Boundary Method. The simulations performed here show that the critical Galileo number for the onset of secondary motions decreases as the spheroid aspect ratio departs from 1. Above this critical threshold, oblate particles perform a zigzagging motion whereas prolate particles rotate around the vertical axis while having their broad side facing the falling direction. Instabilities of the vortices in the wake follow when farther increasing the Galileo number. We also study the drafting-kissing-tumbling associated with the settling of particle pairs. We find that the interaction time increases significantly for non-spherical particles and, more interestingly, spheroidal particles are attracted from larger lateral displacements. This has important implications for the estimation of collision kernels and can result in increasing clustering in suspensions of sedimenting spheroids.

preprint2006arXiv

The effect of wall heating on instability of channel flow

A comprehensive study of the effect of wall heating or cooling on the linear, transient and secondary growth of instability in channel flow is conducted. The effect of viscosity stratification, heat diffusivity and of buoyancy are estimated separately, with some unexpected results. From linear stability results, it has been accepted that heat diffusivity does not affect stability. However, we show that realistic Prandtl numbers cause a transient growth of disturbances that is an order of magnitude higher than at zero Prandtl number. Buoyancy, even at fairly low levels, gives rise to high levels of subcritical energy growth. Unusually for transient growth, both of these are spanwise-independent and not in the form of streamwise vortices. At moderate Grashof numbers, exponential growth dominates, with distinct Rayleigh-Benard and Poiseuille modes for Grashof numbers upto $\sim 25000$, which merge thereafter. Wall heating has a converse effect on the secondary instability compared to the primary, destabilising significantly when viscosity decreases towards the wall. It is hoped that the work will motivate experimental and numerical efforts to understand the role of wall heating in the control of channel and pipe flows.

preprint2016arXiv

Energy cascade in internal wave attractors

One of the pivotal questions in the dynamics of the oceans is related to the cascade of mechanical energy in the abyss and its contribution to mixing. Here, we propose internal wave attractors in the large amplitude regime as a unique self-consistent experimental and numerical setup that models a cascade of triadic interactions transferring energy from large-scale monochro-matic input to multi-scale internal wave motion. We also provide signatures of a discrete wave turbulence framework for internal waves. Finally, we show how beyond this regime, we have a clear transition to a regime of small-scale high-vorticity events which induce mixing. Introduction.

preprint2015arXiv

Collision model for fully-resolved simulations of flows laden with finite-size particles

We present a collision model for particle-particle and particle-wall interactions in interface-resolved simulations of particle-laden flows. Three types of inter-particle interactions are taken into account: (1) long- and (2) short-range hydrodynamic interactions, and (3) solid-solid contact. Long-range interactions are incorporated through an efficient and second-order accurate immersed boundary method (IBM). Short-range interactions are also partly reproduced by the IBM. However, since the IBM uses a fixed-grid, a lubrication model is needed for an inter-particle gap width smaller than the grid spacing. The lubrication model is based on asymptotic expansions of analytical solutions for canonical lubrication interactions between spheres in the Stokes regime. Roughness effects are incorporated by making the lubrication correction independent of the gap width for gap widths smaller than $\sim 1\%$ of the particle radius. This correction is applied until the particles reach solid-solid contact. To model solid-solid contact we use a variant of a linear soft-sphere collision model capable of stretching the collision time. This choice is computationally attractive because it allows to reduce the number of time steps required for integrating the collision force accurately and is physically realistic, provided that the prescribed collision time is much smaller than the characteristic timescale of particle motion. We verified the numerical implementation of our collision model and validated it against several benchmark cases for immersed head-on particle-wall and particle-particle collisions, and oblique particle-wall collisions. The results show good agreement with experimental data.

preprint2013arXiv

Attracting fixed points for heavy particles in the vicinity of a vortex pair

We study the behaviour of heavy inertial particles in the flow field of two like-signed vortices. In a frame co-rotating with the two vortices, we find that stable fixed points exist for these heavy inertial particles; these stable frame-fixed points exist only for particle Stokes number $St<St_{cr}$. We estimate $St_{cr}$ and compare this with direct numerical simulations, and find that the addition of viscosity increases the $St_{cr}$ slightly. We also find that the fixed points become more stable with increasing $St$ until they abruptly disappear at $St=St_{cr}$. These frame-fixed points are between fixed points and limit cycles in character.

preprint2005arXiv

Stability of flow through a slowly diverging pipe

Although the critical Reynolds number for linear instability of the laminar flow in a straight pipe is infinite, we show that it is finite for a divergent pipe, and approaches infinity as the inverse of the divergence angle. The velocity profile at the threshold of inviscid stability is obtained. A non-parallel analysis yields linear instability at surprisingly low Reynolds numbers, of about 150 for a divergence of 3 degrees, which would suggest a role for such instabilities in the transition to turbulence. A multigrid Poisson equation solver is employed for the basic flow, and an extended eigenvalue method for the partial differential equations describing the stability.

preprint2015arXiv

Oscillatory elastic instabilities in an extensional viscoelastic flow

Dilute polymer solutions are known to exhibit purely elastic instabilities even when the fluid inertia is negligible. Here we report the quantitative evidence of two consecutive oscillatory elastic instabilities in an elongation flow of a dilute polymer solution as realized in a T-junction geometry with a long recirculating cavity. The main result reported here is the observation and characterization of the first transition as a forward Hopf bifurcation resulted in a uniformly oscillating state due to breaking of time translational invariance. This unexpected finding is in contrast with previous experiments and numerical simulations performed in similar ranges of the $Wi$ and $Re$ numbers, where the forward fork-bifurcation into a steady asymmetric flow due to the broken spatial inversion symmetry was reported. We discuss the plausible discrepancy between our findings and previous studies that could be attributed to the long recirculating cavity, where the length of the recirculating cavity plays a crucial role in the breaking of time translational invariance instead of the spatial inversion. The second transition is manifested via time aperiodic transverse fluctuations of the interface between the dyed and undyed fluid streams at the channel junction and advected downstream by the mean flow. Both instabilities are characterized by fluid discharge-rate and simultaneous imaging of the interface between the dyed and undyed fluid streams in the outflow channel.

preprint2017arXiv

Instabilities of Internal Gravity Wave Beams

Internal gravity waves play a primary role in geophysical fluids: they contribute significantly to mixing in the ocean and they redistribute energy and momentum in the middle atmosphere. Until recently, most studies were focused on plane wave solutions. However, these solutions are not a satisfactory description of most geophysical manifestations of internal gravity waves, and it is now recognized that internal wave beams with a confined profile are ubiquitous in the geophysical context. We will discuss the reason for the ubiquity of wave beams in stratified fluids, related to the fact that they are solutions of the nonlinear governing equations. We will focus more specifically on situations with a constant buoyancy frequency. Moreover, in light of recent experimental and analytical studies of internal gravity beams, it is timely to discuss the two main mechanisms of instability for those beams. i) The Triadic Resonant Instability generating two secondary wave beams. ii) The streaming instability corresponding to the spontaneous generation of a mean flow.

preprint2017arXiv

Drag reduction in turbulent channel flow laden with finite-size oblate spheroids

We study suspensions of oblate rigid particles in a viscous fluid for different values of the particle volume fractions. Direct numerical simulations have been performed using a direct-forcing immersed boundary method to account for the dispersed phase, combined with a soft-sphere collision model and lubrication corrections for short-range particle-particle and particle-wall interactions. With respect to the single phase flow, we show that in flows laden with oblate spheroids the drag is reduced and the turbulent fluctuations attenuated. In particular, the turbulence activity decreases to lower values than those obtained by only accounting for the effective suspension viscosity. To explain the observed drag reduction we consider the particle dynamics and the interactions of the particles with the turbulent velocity field and show that the particle wall layer, previously observed and found to be responsible for the increased dissipation in suspensions of spheres, disappears in the case of oblate particles. These rotate significantly slower than spheres near the wall and tend to stay with their major axes parallel to the wall, which leads to a decrease of the Reynolds stresses and turbulence production and so to the overall drag reduction.

preprint2018arXiv

Mixing by internal waves quantified using combined PIV/PLIF technique

We present a novel characterization of mixing events associated with the propagation and overturning of internal waves, studied thanks to the simultaneous use of Particle Image Velocimetry (PIV) and Planar Laser Induced Fluorescence (PLIF) techniques. This combination of techniques had been developed earlier to provide an access to simultaneous velocity and density fields in two-layer stratified flows with interfacial gravity waves. Here, for the first time, we show how it is possible to implement it quantitatively in the case of a continuously stratified fluid where internal waves propagate in the bulk. We explain in details how the calibration of the PLIF data is performed by an iterative procedure, and we describe the precise spatial and temporal synchronizations of the PIV and PLIF measurements. We then validate the whole procedure by characterizing the Triadic Resonance Instability (TRI) of an internal wave mode. Very interestingly, the combined technique is then applied to a {precise} measurement of the turbulent diffusivity K_t associated with mixing events induced by an internal wave mode. Values up to K_t=15 mm^2\cdot s^{-1} are reached when TRI is present (well above the noise of our measurement, typically 1 mm^2\cdot s^{-1}), unambiguously confirming that TRI is a potential pathway to turbulent mixing in stratified flows. This work therefore provides a step on the path to new measurements for internal waves.

preprint2018arXiv

Highly accurate numerical computation of implicitly defined volumes using the Laplace-Beltrami operator

This paper introduces a novel method for the efficient and accurate computation of the volume of a domain whose boundary is given by an orientable hypersurface which is implicitly given as the iso-contour of a sufficiently smooth level-set function. After spatial discretization, local approximation of the hypersurface and application of the Gaussian divergence theorem, the volume integrals are transformed to surface integrals. Application of the surface divergence theorem allows for a further reduction to line integrals which are advantageous for numerical quadrature. We discuss the theoretical foundations and provide details of the numerical algorithm. Finally, we present numerical results for convex and non-convex hypersurfaces embedded in cuboidal domains, showing both high accuracy and thrid- to fourth-order convergence in space.

preprint2016arXiv

Linear stability analysis and direct numerical simulation of two layer channel flow

We study the stability of two-fluid flow through a plane channel at Reynolds numbers of a hundred to a thousand in the linear and nonlinear regimes. The two fluids have the same density but different viscosities. The fluids, when miscible, are separated from each other by a mixed layer of small but finite thickness, across which viscosity changes from that of one fluid to that of the other. When immiscible, the interface is sharp. Our study spans a range of Schmidt numbers, viscosity ratios and location and thickness of the mixed layer. A region of instability distinct from that of the Tollmien-Schlichting mode is obtained at moderate Reynolds numbers. We show that the overlap of the layer of viscosity-stratification with the critical layer of the dominant disturbance provides a mechanism for this instability. At very low values of diffusivity, the miscible flow behaves exactly like the immiscible in terms of stability characteristics. High levels of miscibility make the flow more stable. At intermediate levels of diffusivity however, in both linear and non-linear regimes, miscible flow can be more unstable than the corresponding immiscible flow without surface tension.

preprint2018arXiv

Internal Wave Attractors in 3D Geometries : trapping by oblique reflection

We study experimentally the propagation of internal waves in two different three-dimensional (3D) geometries, with a special emphasis on the refractive focusing due to the 3D reflection of obliquely incident internal waves on a slope. Both studies are initiated by ray tracing calculations to determine the appropriate experimental parameters. First, we consider a 3D geometry, the classical set-up to get simple, 2D parallelogram-shaped attractors in which waves are forced in a direction perpendicular to a sloping bottom. Here, however, the forcing is of reduced extent in along-slope, transverse direction. We show how the refractive focusing mechanism explains the formation of attractors over the whole width of the tank, even away from the forcing region. Direct numerical simulations confirm the dynamics,emphasize the role of boundary conditions and reveal the phase shifting in the transversal direction. Second, we consider a long and narrow tank having an inclined bottom, to simply reproduce a canal. In this case, the energy is injected in a direction parallel to the slope. Interestingly, the wave energy ends up forming 2D internal wave attractors in planes that are transverse to the initial propagation direction. This focusing mechanism prevents indefinite transmission of most of the internal wave energy along the canal.

preprint2018arXiv

Added mass: a complex face of tidal conversion

This paper revisits the problem of tidal conversion at a ridge in a uniformly stratified fluid of limited depth using measurements of complex-valued added mass. When the height of a sub-marine ridge is non negligible with respect to the depth of the water, the tidal conversion can be enhanced in the supercritical regime or reduced in the subcritical regime with respect to the large depth situation. Tidal conversion can even be null for some specific cases. Here, we study experimentally the influence of finite depth on the added mass coefficients for three diffierent ridge shapes. We first show that at low forcing frequency the tidal conversion is weakly enhanced by shallow depth for a semi-circular ridge. In addition, added mass coefficients measured for a vertical ridge show strong similarities with the ones obtained for the semi-circular ridge. Nevertheless, the enhancement of the tidal conversion at low forcing frequency for the vertical ridge has not been observed, in contrast with its supercritical shape. Finally, we provide the experimental evidence of a lack of tidal conversion due to the specific shape of a ridge for certain depth and frequency tuning.

preprint2017arXiv

Turbulent channel flow of a dense binary mixture of rigid particles

We study turbulent channel flow of a binary mixture of finite-size neutrally-buoyant rigid particles by means of interface-resolved direct numerical simulations. We fix the bulk Reynolds number and total solid volume fraction, $Re_b = 5600$ and $Φ=20\%$, and vary the relative fraction of small and large particles. The binary mixture consists of particles of two different sizes, $2h/d_l=20$ and $2h/d_s=30$ where $h$ is the half channel height and $d_l$ and $d_s$ the diameter of the large and small particles. While the particulate flow statistics exhibit a significant alteration of the mean velocity profile and turbulent fluctuations with respect to the unladen flow, the differences between the mono-disperse and bi-disperse cases are small. However, we observe a clear segregation of small particles at the wall in binary mixtures, which affects the dynamics of the near wall region and thus the overall drag. This results in a higher drag in suspensions with a larger amount of large particles. As regards bi-disperse effects on the particle dynamics, a non-monotonic variation of the particle dispersion in the spanwise (homogeneous) direction is observed when increasing the percentage of small/large particles. Finally, we note that particles of the same size tend to cluster more at contact whereas the dynamics of the large particles gives highest collision kernels due to a higher approaching speed.

preprint2017arXiv

Exact solitonic and periodic solutions of the extended KdV equation

The KdV equation can be derived in the shallow water limit of the Euler equations. Over the last few decades, this equation has been extended to include both higher order effects (KdV2) and an uneven river bottom. Although this equation is not integrable and has only one conservation law, exact periodic and solitonic solutions exist for the even bottom case. The method used to find them assumes the same function forms as for KdV solutions. KdV2 equation imposes more constraints on parameters of solutions. For soliton case KdV2 solution occurs for particular ratio of wave amplitude to water depth only. For periodic case physically relevant solutions are admissible only for two narrow intervals of elliptic parameter $m$. For a range of $m$ near one the cnoidal waves are upright as expected, but are inverted in $m$ region close to zero. Properties of exact solutions of KdV and KdV2 are compared.

preprint2018arXiv

Effects of the finite particle size in turbulent wall-bounded flows of dense suspensions

We use interface-resolved simulations to study finite-size effects in turbulent channel flow of neutrally-buoyant spheres. Two cases with particle sizes differing by a factor of 2, at the same solid volume fraction of 20% and bulk Reynolds number are considered. These are complemented with two reference single-phase flows: the unladen case, and the flow of a Newtonian fluid with the effective suspension viscosity of the same mixture in the laminar regime. As recently highlighted in Costa et al. (PRL 117, 134501), a particle-wall layer is responsible for deviations of the statistics from what is observed in the continuum limit where the suspension is modeled as a Newtonian fluid with an effective viscosity. Here we investigate the fluid and particle dynamics in this layer and in the bulk. In the particle-wall layer, the near wall inhomogeneity has an influence on the suspension micro-structure over a distance proportional to the particle size. In this layer, particles have a significant (apparent) slip velocity that is reflected in the distribution of wall shear stresses. This is characterized by extreme events (both much higher and much lower than the mean). Based on these observations we provide a scaling for the particle-to-fluid apparent slip velocity as a function of the flow parameters. We also extend the flow scaling laws in to second-order Eulerian statistics in the homogeneous suspension region away from the wall. Finite-size effects in the bulk of the channel become important for larger particles, while negligible for lower-order statistics and smaller particles. Finally, we study the particle dynamics along the wall-normal direction. Our results suggest that 1-point dispersion is dominated by particle-turbulence (and not particle-particle) interactions, while differences in 2-point dispersion and collisional dynamics are consistent with a picture of shear-driven interactions.

preprint2018arXiv

Mixing and formation of layers by internal wave forcing

The energy pathways from propagating internal waves to the scales of irreversible mixing in the ocean are not fully described. In the ocean interior, the triadic resonant instability is an intrinsic destabilization process that may enhance the energy cascade away from topographies. The present study focuses on the integrated impact of mixing processes induced by a propagative normal mode-1 over long term experiments in an idealised setup. The internal wave dynamics and the evolution of the density profile are followed using the light attenuation technique. Diagnostics of the turbulent diffusivity $K_{T}$ and background potential energy $BPE$ are provided. Mixing effects result in a partially mixed layer colocated with the region of maximum shear induced by the forcing normal mode. The maximum measured turbulent diffusivity is 250 times larger than the molecular value, showing that diapycnal mixing is largely enhanced by small scale turbulent processes. Intermittency and reversible energy transfers are discussed to bridge the gap between the present diagnostic and the larger values measured in Dossmann et al, Experiments in Fluids, 57(8), 132 (2016). The mixing efficiency $η$ is assessed by relating the $BPE$ growth to the linearized $KE$ input. One finds a value of $Γ=12-19\%$ larger than the mixing efficiency in the case of breaking interfacial wave. After several hours of forcing, the development of staircases in the density profile is observed. This mechanism has been previously observed in experiments with weak homogeneous turbulence and explained by argument. The present experiments suggest that internal wave forcing could also induce the formation of density interfaces in the ocean.

preprint2018arXiv

Scale effects in internal wave attractors

As a necessary preliminary step toward geophysically significant extrapolations, we study the scale effects in internal wave attractors in the linear and nonlinear regimes. We use two geometrically similar experimental set-ups, scaled to factor 3, and numerical simulations (a spectral element method, based on the Nek5000 open solver) for a range of parameters that is typically accessible in laboratory. In the linear regime, we recover the classical viscous scaling for the beam width, which is not affected by variations of the amplitude of the input perturbation. In the nonlinear regime, we show that the scaling of the width-to-length ratio of the attractor branches is intimately related with the energy cascade from large-scale energy input to dissipation. We present results for the wavelength, amplitude and width of the beam as a function of time and as a function of the amplitude of the forcing.

preprint2017arXiv

A new analytical approach for modelling the added mass and hydrodynamic interaction of two cylinders subjected to large motions in a potential stagnant fluid

A potential theory is presented for the problem of two moving cylinders, with possibly different radii, large motions, immersed in an perfect stagnant fluid. We show that the fluid force is the superposition of an added mass term, related to the time variations of the potential, and a quadratic term related to its spatial variations. We provide new simple and exact analytical expressions for the fluid added mass coefficients, in which the effect of the confinement is made explicit. The self-added mass (resp. cross-added mass) is shown to decrease (resp. increase) with the separation distance and increase (resp. decreases) with the radius ratio. We then consider the case in which one cylinder translates along the line joining the centers with a constant speed. We show that the two cylinders are repelled from each other, with a force that diverges to infinity at impact. We extend our approach to the case in which one cylinder is imposed a sinusoidal vibration. We show that the force on the stationnary cylinder and the vibration displacement have opposite (resp. identical) axial (resp. transverse) directions. For large vibration amplitudes, this force is strongly altered by the nonlinear effects induced by the spatial variations of the potential. The force on the vibrating cylinder is in phase with the imposed displacement and is mainly driven by the added mass term. The results of this paper are of particular interest for engineers who need to grab the essential features associated to the vibration of a solid body in a still fluid.