JFM Papers
The optimal kinematic dynamo driven by steady flows in a sphere
- L. Chen, W. Herreman, K. Li, P. W. Livermore, J. W. Luo, A. Jackson
-
- Published online by Cambridge University Press:
- 25 January 2018, pp. 1-32
-
- Article
- Export citation
-
We present a variational optimization method that can identify the most efficient kinematic dynamo in a sphere, where efficiency is based on the value of a magnetic Reynolds number that uses enstrophy to characterize the inductive effects of the fluid flow. In this large-scale optimization, we restrict the flow to be steady and incompressible, and the boundary of the sphere to be no-slip and electrically insulating. We impose these boundary conditions using a Galerkin method in terms of specifically designed vector field bases. We solve iteratively for the flow field and the accompanying magnetic eigenfunction in order to find the minimal critical magnetic Reynolds number $Rm_{c,min}$ for the onset of a dynamo. Although nonlinear, this iteration procedure converges to a single solution and there is no evidence that this is not a global optimum. We find that $Rm_{c,min}=64.45$ is at least three times lower than that of any published example of a spherical kinematic dynamo generated by steady flows, and our optimal dynamo clearly operates above the theoretical lower bounds for dynamo action. The corresponding optimal flow has a spatially localized helical structure in the centre of the sphere, and the dominant components are invariant under rotation by $\unicode[STIX]{x03C0}$.
Mechanisms in the hypersonic laminar near wake of a blunt body
- W. Schuyler Hinman, Craig T. Johansen
-
- Published online by Cambridge University Press:
- 25 January 2018, pp. 33-75
-
- Article
- Export citation
-
A new theoretical framework, based on the analysis of Navier–Stokes solutions for the hypersonic laminar near wake of two-dimensional and axisymmetric blunt bodies, is presented. A semi-empirical relationship is derived between the free-stream Mach and Reynolds numbers and a characteristic wake Reynolds number. A control volume analysis was performed to assess the validity of some common assumptions used in the literature. Analysis of the momentum and vorticity equations is used to assess the dominant mechanisms of momentum transfer along and across the dividing streamline and centreline which enclose the near wake. An observed stagnation pressure gain along the dividing streamline is explained using the entropy transport equation, demonstrating an unbalance between entropy generation due to viscous dissipation and entropy diffusion. The rear-stagnation point flow is analysed using an analogy to a reversed flow jet which allows for the centreline Mach number to be solved. A new viscous–inviscid interaction theory is presented for the reattachment shock formation process for both planar and axisymmetric wakes. Finally, all of the sub-mechanisms are combined into an overall wake mechanism. The resulting equations constitute the first overall theoretical framework of the laminar near-wake mechanism including separation, reattachment, rear-stagnation point flow and dividing streamline stagnation pressure gain for both planar and axisymmetric near wakes. Scaling arguments are presented throughout the work for each of the key sub-mechanisms. Recommendations are made for how experimental and numerical results for the near wake should be presented. The equations and recommendations presented here are then used to perform a detailed disambiguation of laminar capsule studies in the literature.
The critical point of the transition to turbulence in pipe flow
- Vasudevan Mukund, Björn Hof
-
- Published online by Cambridge University Press:
- 25 January 2018, pp. 76-94
-
- Article
- Export citation
-
In pipes, turbulence sets in despite the linear stability of the laminar Hagen–Poiseuille flow. The Reynolds number ($Re$) for which turbulence first appears in a given experiment – the ‘natural transition point’ – depends on imperfections of the set-up, or, more precisely, on the magnitude of finite amplitude perturbations. At onset, turbulence typically only occupies a certain fraction of the flow, and this fraction equally is found to differ from experiment to experiment. Despite these findings, Reynolds proposed that after sufficiently long times, flows may settle to steady conditions: below a critical velocity, flows should (regardless of initial conditions) always return to laminar, while above this velocity, eddying motion should persist. As will be shown, even in pipes several thousand diameters long, the spatio-temporal intermittent flow patterns observed at the end of the pipe strongly depend on the initial conditions, and there is no indication that different flow patterns would eventually settle to a (statistical) steady state. Exploiting the fact that turbulent puffs do not age (i.e. they are memoryless), we continuously recreate the puff sequence exiting the pipe at the pipe entrance, and in doing so introduce periodic boundary conditions for the puff pattern. This procedure allows us to study the evolution of the flow patterns for arbitrary long times, and we find that after times in excess of $10^{7}$ advective time units, indeed a statistical steady state is reached. Although the resulting flows remain spatio-temporally intermittent, puff splitting and decay rates eventually reach a balance, so that the turbulent fraction fluctuates around a well-defined level which only depends on $Re$. In accordance with Reynolds’ proposition, we find that at lower $Re$ (here 2020), flows eventually always resume to laminar, while for higher $Re$ (${\geqslant}2060$), turbulence persists. The critical point for pipe flow hence falls in the interval of $2020<Re<2060$, which is in very good agreement with the recently proposed value of $Re_{c}=2040$. The latter estimate was based on single-puff statistics and entirely neglected puff interactions. Unlike in typical contact processes where such interactions strongly affect the percolation threshold, in pipe flow, the critical point is only marginally influenced. Interactions, on the other hand, are responsible for the approach to the statistical steady state. As shown, they strongly affect the resulting flow patterns, where they cause ‘puff clustering’, and these regions of large puff densities are observed to travel across the puff pattern in a wave-like fashion.
Geometric generalised Lagrangian-mean theories
- Andrew D. Gilbert, Jacques Vanneste
-
- Published online by Cambridge University Press:
- 25 January 2018, pp. 95-134
-
- Article
- Export citation
-
Many fluctuation-driven phenomena in fluids can be analysed effectively using the generalised Lagrangian-mean (GLM) theory of Andrews & McIntyre (J. Fluid Mech., vol. 89, 1978, pp. 609–646) This finite-amplitude theory relies on particle-following averaging to incorporate the constraints imposed by the material conservation of certain quantities in inviscid regimes. Its original formulation, in terms of Cartesian coordinates, relies implicitly on an assumed Euclidean structure; as a result, it does not have a geometrically intrinsic, coordinate-free interpretation on curved manifolds, and suffers from undesirable features. Motivated by this, we develop a geometric generalisation of GLM that we formulate intrinsically using coordinate-free notation. One benefit is that the theory applies to arbitrary Riemannian manifolds; another is that it establishes a clear distinction between results that stem directly from geometric consistency and those that depend on particular choices. Starting from a decomposition of an ensemble of flow maps into mean and perturbation, we define the Lagrangian-mean momentum as the average of the pull-back of the momentum one-form by the perturbation flow maps. We show that it obeys a simple equation which guarantees the conservation of Kelvin’s circulation, irrespective of the specific definition of the mean flow map. The Lagrangian-mean momentum is the integrand in Kelvin’s circulation and distinct from the mean velocity (the time derivative of the mean flow map) which advects the contour of integration. A pseudomomentum consistent with that in GLM can then be defined by subtracting the Lagrangian-mean momentum from the one-form obtained from the mean velocity using the manifold’s metric. The definition of the mean flow map is based on choices made for reasons of convenience or aesthetics. We discuss four possible definitions: a direct extension of standard GLM, a definition based on optimal transportation, a definition based on a geodesic distance in the group of volume-preserving diffeomorphisms, and the ‘glm’ definition proposed by Soward & Roberts (J. Fluid Mech., vol. 661, 2010, pp. 45–72). Assuming small-amplitude perturbations, we carry out order-by-order calculations to obtain explicit expressions for the mean velocity and Lagrangian-mean momentum at leading order. We also show how the wave-action conservation of GLM extends to the geometric setting. To make the paper self-contained, we introduce in some detail the tools of differential geometry and main ideas of geometric fluid dynamics on which we rely. These include variational formulations which we use for alternative derivations of some key results. We mostly focus on the Euler equations for incompressible inviscid fluids but sketch out extensions to the rotating–stratified Boussinesq, compressible Euler, and magnetohydrodynamic equations. We illustrate our results with an application to the interaction of inertia-gravity waves with balanced mean flows in rotating–stratified fluids.
Revisiting slope influence in turbulent bedload transport: consequences for vertical flow structure and transport rate scaling
- Raphael Maurin, Julien Chauchat, Philippe Frey
-
- Published online by Cambridge University Press:
- 25 January 2018, pp. 135-156
-
- Article
- Export citation
-
Gravity-driven turbulent bedload transport has been extensively studied over the past century in regard to its importance for Earth surface processes such as natural riverbed morphological evolution. In the present contribution, the influence of the longitudinal channel inclination angle on gravity-driven turbulent bedload transport is studied in an idealised framework considering steady and uniform flow conditions. From an analytical analysis based on the two-phase continuous equations, it is shown that: (i) the classical slope correction of the critical Shields number is based on an erroneous formulation of the buoyancy force, (ii) the influence of the slope is not restricted to the critical Shields number but affects the whole transport formula and (iii) pressure-driven and gravity-driven turbulent bedload transport are not equivalent from the slope influence standpoint. Analysing further the granular flow driving mechanisms, the longitudinal slope is shown to not only influence the fluid bed shear stress and the resistance of the granular bed, but also to affect the fluid flow inside the granular bed – responsible for the transition from bedload transport to debris flow. The relative influence of these coupled mechanisms allows us to understand the evolution of the vertical structure of the granular flow and to predict the transport rate scaling law as a function of a rescaled Shields number. The theoretical analysis is validated with coupled fluid–discrete element simulations of idealised gravity-driven turbulent bedload transport, performed over a wide range of Shields number values, density ratios and channel inclination angles. In particular, all the data are shown to collapse onto a master curve when considering the sediment transport rate as a function of the proposed rescaled Shields number.
Propagation of a strong shock over a random bed of spherical particles
- Y. Mehta, C. Neal, K. Salari, T. L. Jackson, S. Balachandar, S. Thakur
-
- Published online by Cambridge University Press:
- 25 January 2018, pp. 157-197
-
- Article
- Export citation
-
Propagation of a strong incident shock through a bed of particles results in complex wave dynamics such as a reflected shock, a transmitted shock, and highly unsteady flow inside the particle bed. In this paper we present three-dimensional numerical simulations of shock propagation in air over a random bed of particles. We assume the flow is inviscid and governed by the Euler equations of gas dynamics. Simulations are carried out by varying the volume fraction of the particle bed at a fixed shock Mach number. We compute the unsteady inviscid streamwise and transverse drag coefficients as a function of time for each particle in the random bed for different volume fractions. We show that (i) there are significant variations in the peak drag for the particles in the bed, (ii) the mean peak drag as a function of streamwise distance through the bed decreases with a slope that increases as the volume fraction increases, and (iii) the deviation from the mean peak drag does not correlate with local volume fraction. We also present the local Mach number and pressure contours for the different volume fractions to explain the various observed complex physical mechanisms occurring during the shock–particle interactions. Since the shock interaction with the random bed of particles leads to transmitted and reflected waves, we compute the average flow properties to characterize the strength of the transmitted and reflected shock waves and quantify the energy dissipation inside the particle bed. Finally, to better understand the complex wave dynamics in a random bed, we consider a simpler approximation of a planar shock propagating in a duct with a sudden area change. We obtain Riemann solutions to this problem, which are used to compare with fully resolved numerical simulations.
Two-level, two-phase model for intense, turbulent sediment transport
- Jose M. Gonzalez-Ondina, Luigi Fraccarollo, Philip L.-F. Liu
-
- Published online by Cambridge University Press:
- 26 January 2018, pp. 198-238
-
- Article
- Export citation
-
The study of sediment transport requires in-depth investigation of the complex effects of sediment particles in fluid turbulence. In this paper we focus on intense sediment transport flows. None of the existing two-phase models in the literature properly replicates the liquid and solid stresses in the near bed region of high concentration of sediment. The reason for this shortcoming is that the physical processes occurring at the length scale of the particle collisions are different from those occurring at larger length scales and therefore, they must be modelled independently. We present here a two-level theoretical derivation of two-phase, Favre averaged Navier–Stokes equations (FANS). This approach treats two levels of energy fluctuations independently, those associated with a granular spatial scale (granular temperature and small-scale fluid turbulence) and those associated with the ensemble average (turbulent kinetic energy for the two phases). Although similar attempts have been made by other researchers, the two level approach ensures that the two relevant length scales are included independently in a more consistent manner. The model is endowed with a semi-empirical formulation for the granular scale fluid turbulence, which is important even in the dense collisional shear layer, as has been recently recognized. As a result of the large and small scale modelling of the liquid and solid fluctuations, predictions are promising to be reliable in a wide range of flow conditions, from collisional to turbulent suspensions. This model has been validated for steady state flows with intense, collisional or mixed collisional–turbulent sediment transport, using various sources of detailed experimental data. It compares well with the experimental results in the whole experimental range of Shields parameters, better than previous models, although at the cost of increased complexity in the equations. Further experiments on turbulent suspensions would be necessary to definitely assess the model capabilities.
Rapidly rotating precessing cylinder flows: forced triadic resonances
- Juan M. Lopez, Francisco Marques
-
- Published online by Cambridge University Press:
- 25 January 2018, pp. 239-270
-
- Article
- Export citation
-
Rapidly rotating cylinder flows subjected to low-amplitude precessional forcing are studied numerically over a range of cylinder and precessional rotation rates. For sufficiently small rotation rates, viscous effects lead to a forced overturning flow that is steady in the precession (table) frame of reference. Increasing the rotation rates, this forced flow loses stability in a Hopf bifurcation, which can be either supercritical or subcritical, and may preserve or break the symmetry of the system, depending on the parameter regime studied. Regardless of these details of the Hopf bifurcation, it is found that the Hopf instability is associated with a slightly detuned triadic resonance between the forced overturning flow and two free Kelvin modes (inviscid eigenmodes of the rotating cylinder). Further increases in rotation rates lead to a sequence of secondary instabilities which also follow a generic pattern irrespective of the parameter regime investigated. The relationship between this sequence of instabilities and the resultant nonlinear dynamics with the experimentally observed phenomenon of resonant collapse is discussed.
Effects of Reynolds number and Stokes number on particle-pair relative velocity in isotropic turbulence: a systematic experimental study
- Zhongwang Dou, Andrew D. Bragg, Adam L. Hammond, Zach Liang, Lance R. Collins, Hui Meng
-
- Published online by Cambridge University Press:
- 26 January 2018, pp. 271-292
-
- Article
- Export citation
-
The effects of Reynolds number ($R_{\unicode[STIX]{x1D706}}$) and Stokes number ($St$) on particle-pair relative velocity (RV) are investigated systematically using a recently developed planar four-frame particle tracking technique in a novel homogeneous and isotropic turbulence chamber. We compare the measured results with direct numerical simulation (DNS), verifying whether the conclusions of the DNS for simplified conditions and limited $R_{\unicode[STIX]{x1D706}}$ are still valid in reality. Two experiments are performed: varying $R_{\unicode[STIX]{x1D706}}$ between 246 and 357 at six $St$ values, and varying $St$ between 0.02 and 4.63 at five $R_{\unicode[STIX]{x1D706}}$ values. The measured mean inward particle-pair RV $\langle w_{r}^{-}\rangle$ as a function of separation distance $r$ is compared with the DNS under closely matched conditions. At all experimental conditions, an excellent agreement is achieved, except when the particle separation distance $r\lesssim 10\unicode[STIX]{x1D702}$ ($\unicode[STIX]{x1D702}$ is the Kolmogorov length scale), where the experimental $\langle w_{r}^{-}\rangle$ is consistently higher, possibly due to particle polydispersity and finite laser thickness in the experiments (Dou et al., arXiv:1712.07506, 2017). At any fixed $St,\langle w_{r}^{-}\rangle$ is essentially independent of $R_{\unicode[STIX]{x1D706}}$, echoing the DNS finding of Ireland et al. (J. Fluid Mech., vol. 796, 2016, pp. 617–658). At any fixed $R_{\unicode[STIX]{x1D706}}$, $\langle w_{r}^{-}\rangle$ increases with $St$ at small $r$, showing dominance of the path-history effect in the dissipation range when $St\gtrsim O(1)$, but decreases with $St$ at large $r$, indicating dominance of inertial filtering. We further compare the $\langle w_{r}^{-}\rangle$ and RV variance $\langle w_{r}^{2}\rangle$ from experiments with DNS and theoretical predictions by Pan & Padoan (J. Fluid Mech., vol. 661, 2010, pp. 73–107). For $St\lesssim 1$, experimental $\langle w_{r}^{-}\rangle$ and $\langle w_{r}^{2}\rangle$ match these values well at $r\gtrsim 10\unicode[STIX]{x1D702}$, but they are higher than both DNS and theory at $r\lesssim 10\unicode[STIX]{x1D702}$. For $St\gtrsim 1$, $\langle w_{r}^{-}\rangle$ from all three match well, except for $r\lesssim 10\unicode[STIX]{x1D702}$, for which experimental values are higher, while $\langle w_{r}^{2}\rangle$ from experiment and DNS are much higher than theoretical predictions. We discuss potential causes of these discrepancies. What this study shows is the first experimental validation of $R_{\unicode[STIX]{x1D706}}$ and $St$ effect on inertial particle-pair $\langle w_{r}^{-}\rangle$ in homogeneous and isotropic turbulence.
Water wave overwash of a step
- D. M. Skene, L. G. Bennetts, M. Wright, M. H. Meylan, K. J. Maki
-
- Published online by Cambridge University Press:
- 29 January 2018, pp. 293-312
-
- Article
- Export citation
-
Water wave overwash of a step by small steepness, regular incident waves is analysed using a computational fluid dynamics (CFD) model and a mathematical model, in two spatial dimensions. The CFD model is based on the two-phase, incompressible Navier–Stokes equations, and the mathematical model is based on the coupled potential-flow and nonlinear shallow-water theories. The CFD model is shown to predict vortices, breaking and overturning in the region where overwash is generated, and that the overwash develops into fast-travelling bores. The mathematical model is shown to predict bore heights and velocities that agree with the CFD model, despite neglecting the complicated dynamics where the overwash is generated. Evidence is provided to explain the agreement in terms of the underlying agreement of mass and energy fluxes.
Vortex dynamics and sound emission in excited high-speed jets
- Michael Crawley, Lior Gefen, Ching-Wen Kuo, Mo Samimy, Roberto Camussi
-
- Published online by Cambridge University Press:
- 29 January 2018, pp. 313-347
-
- Article
- Export citation
-
This work aims to study the dynamics of and noise generated by large-scale structures in a Mach 0.9 turbulent jet of Reynolds number $6.2\times 10^{5}$ using plasma-based excitation of shear layer instabilities. The excitation frequency is varied to produce individual or periodic coherent ring vortices in the shear layer. First, two-point cross-correlations are used between the acoustic near field and far field in order to identify the dominant noise source region. The large-scale structure interactions are then investigated by stochastically estimating time-resolved velocity fields using time-resolved near-field pressure traces and non-time-resolved planar velocity snapshots (obtained by particle image velocimetry) by means of an artificial neural network. The estimated time-resolved velocity fields show multiple mergings of large-scale structures in the shear layer, and indicate that disintegration of coherent ring vortices is the dominant aeroacoustic source mechanism for the jet studied here. However, the merging of vortices in the initial shear layer is also identified as a non-trivial noise source mechanism.
Dispersion coefficient in an electro-osmotic flow of a viscoelastic fluid through a microchannel with a slowly varying wall zeta potential
- J. C. Arcos, F. Méndez, E. G. Bautista, O. Bautista
-
- Published online by Cambridge University Press:
- 29 January 2018, pp. 348-386
-
- Article
- Export citation
-
The dispersion coefficient of a passive solute in a steady-state pure electro-osmotic flow (EOF) of a viscoelastic liquid, whose rheological behaviour follows the simplified Phan-Thien–Tanner (sPTT) model, along a parallel flat plate microchannel, is studied. The walls of the microchannel are assumed to have modulated and low $\unicode[STIX]{x1D701}$ potentials, which vary slowly in the axial direction in a sinusoidal manner. The flow field required to obtain the dispersion coefficient was solved using the lubrication approximation theory (LAT). The solution of the electric potential is based on the Debye–Hückel approximation for a symmetric $(z:z)$ electrolyte. The viscoelasticity of the fluid is observed to notably amplify the axial distribution of the effective dispersion coefficients due to the variation in the $\unicode[STIX]{x1D701}$ potentials of the walls. The problem was formulated for two cases: when the Debye layer thickness (EDL) was on the order of unity (thick EDL) and in the limit where the thickness of the EDL was very small compared with the height of the microchannel (thin EDL limit). Due to the coupling between the nonlinear governing equations and the sPTT fluid model, they were replaced by their approximate linearized forms and solved in the limit of $\unicode[STIX]{x1D700}\ll 1$ using the regular perturbation technique. Here $\unicode[STIX]{x1D700}$ is the amplitude of the sinusoidal function of the $\unicode[STIX]{x1D701}$ potentials. Additionally, the numerical solution of the simplified governing equations was also obtained for $\unicode[STIX]{x1D700}=O(1)$ and compared with the approximate solution, showing excellent agreement for $0\leqslant \unicode[STIX]{x1D700}\leqslant 0.3$. Note that the dispersion coefficient primarily depends on the Deborah number, on the ratio of the half-height of the microchannel to the Debye length, and on the assumed variation in the $\unicode[STIX]{x1D701}$ potentials of the walls.
Internal solitary wave generation by tidal flow over topography
- R. Grimshaw, K. R. Helfrich
-
- Published online by Cambridge University Press:
- 29 January 2018, pp. 387-407
-
- Article
- Export citation
-
Oceanic internal solitary waves are typically generated by barotropic tidal flow over localised topography. Wave generation can be characterised by the Froude number $F=U/c_{0}$, where $U$ is the tidal flow amplitude and $c_{0}$ is the intrinsic linear long wave phase speed, that is the speed in the absence of the tidal current. For steady tidal flow in the resonant regime, $\unicode[STIX]{x1D6E5}_{m}<F-1<\unicode[STIX]{x1D6E5}_{M}$, a theory based on the forced Korteweg–de Vries equation shows that upstream and downstream propagating undular bores are produced. The bandwidth limits $\unicode[STIX]{x1D6E5}_{m,M}$ depend on the height (or depth) of the topographic forcing term, which can be either positive or negative depending on whether the topography is equivalent to a hole or a sill. Here the wave generation process is studied numerically using a forced Korteweg–de Vries equation model with time-dependent Froude number, $F(t)$, representative of realistic tidal flow. The response depends on $\unicode[STIX]{x1D6E5}_{max}=F_{max}-1$, where $F_{max}$ is the maximum of $F(t)$ over half of a tidal cycle. When $\unicode[STIX]{x1D6E5}_{max}<\unicode[STIX]{x1D6E5}_{m}$ the flow is always subcritical and internal solitary waves appear after release of the downstream disturbance. When $\unicode[STIX]{x1D6E5}_{m}<\unicode[STIX]{x1D6E5}_{max}<\unicode[STIX]{x1D6E5}_{M}$ the flow reaches criticality at its peak, producing upstream and downstream undular bores that are released as the tide slackens. When $\unicode[STIX]{x1D6E5}_{max}>\unicode[STIX]{x1D6E5}_{M}$ the tidal flow goes through the resonant regime twice, producing undular bores with each passage. The numerical simulations are for both symmetrical topography, and for asymmetric topography representative of Stellwagen Bank and Knight Inlet.
Wave-induced mean flows in rotating shallow water with uniform potential vorticity
- Jim Thomas, Oliver Bühler, K. Shafer Smith
-
- Published online by Cambridge University Press:
- 01 February 2018, pp. 408-429
-
- Article
- Export citation
-
Theoretical and numerical computations of the wave-induced mean flow in rotating shallow water with uniform potential vorticity are presented, with an eye towards applications in small-scale oceanography where potential-vorticity anomalies are often weak compared to the waves. The asymptotic computations are based on small-amplitude expansions and time averaging over the fast wave scale to define the mean flow. Importantly, we do not assume that the mean flow is balanced, i.e. we compute the full mean-flow response at leading order. Particular attention is paid to the concept of modified diagnostic relations, which link the leading-order Lagrangian-mean velocity field to certain wave properties known from the linear solution. Both steady and unsteady wave fields are considered, with specific examples that include propagating wavepackets and monochromatic standing waves. Very good agreement between the theoretical predictions and direct numerical simulations of the nonlinear system is demonstrated. In particular, we extend previous studies by considering the impact of unsteady wave fields on the mean flow, and by considering the total kinetic energy of the mean flow as a function of the rotation rate. Notably, monochromatic standing waves provide an explicit counterexample to the often observed tendency of the mean flow to decrease monotonically with the background rotation rate.
A multifractal model for the velocity gradient dynamics in turbulent flows
- Rodrigo M. Pereira, Luca Moriconi, Laurent Chevillard
-
- Published online by Cambridge University Press:
- 01 February 2018, pp. 430-467
-
- Article
- Export citation
-
We develop a stochastic model for the velocity gradient dynamics along a Lagrangian trajectory in isotropic and homogeneous turbulent flows. Comparing with different attempts proposed in the literature, the present model, at the cost of introducing a free parameter known in turbulence phenomenology as the intermittency coefficient, gives a realistic picture of velocity gradient statistics at any Reynolds number. To achieve this level of accuracy, we use as a first modelling step a regularized self-stretching term in the framework of the recent fluid deformation (RFD) approximation that was shown to give a realistic picture of small-scale statistics of turbulence only up to moderate Reynolds numbers. As a second step, we constrain the dynamics, in the spirit of Girimaji & Pope (Phys. Fluids A, vol. 2, 1990, p. 242), in order to impose a peculiar statistical structure to the dissipation seen by the Lagrangian particle. This probabilistic closure uses as a building block a random field that fulfils the statistical description of the intermittency, i.e. multifractal, phenomenon. To do so, we define and generalize to a statistically stationary framework a proposition made by Schmitt (Eur. Phys. J. B, vol. 34, 2003, p. 85). These considerations lead us to propose a nonlinear and non-Markovian closed dynamics for the elements of the velocity gradient tensor. We numerically integrate this dynamics and observe that a stationary regime is indeed reached, in which (i) the gradient variance is proportional to the Reynolds number, (ii) gradients are typically correlated over the (small) Kolmogorov time scale and gradient norms over the (large) integral time scale, (iii) the joint probability distribution function of the two non-vanishing invariants $Q$ and $R$ reproduces the characteristic teardrop shape, (iv) vorticity becomes preferentially aligned with the intermediate eigendirection of the deformation tensor and (v) gradients are strongly non-Gaussian and intermittent, a behaviour that we quantify by appropriate high-order moments. Additionally, we examine the problem of rotation rate statistics of (axisymmetric) anisotropic particles as observed in direct numerical simulations. Although our realistic picture of velocity gradient fluctuations leads to better results when compared to the former RFD approximation, it is still unable to provide an accurate description for the rotation rate variance of oblate spheroids.
Numerical and theoretical analyses of the dynamics of droplets driven by electrowetting on dielectric in a Hele-Shaw cell
- Yasufumi Yamamoto, Takahiro Ito, Tatsuro Wakimoto, Kenji Katoh
-
- Published online by Cambridge University Press:
- 01 February 2018, pp. 468-488
-
- Article
- Export citation
-
Droplet movement by electrowetting on dielectric (EWOD) in a Hele-Shaw cell is analysed theoretically and numerically. We propose a simple theoretical model for the motion, which describes well the voltage dependency of droplet speed below the saturation voltage as measured experimentally. The simulation method for numerical analyses is constructed by using the Young–Lippmann equation to represent EWOD and the generalised Navier boundary condition to represent the moving contact line in the context of the front-tracking method. With an adjusted slip parameter, the present full three-dimensional numerical simulation reproduces well the shape evolution and movement speed of droplets as observed experimentally. We verify the proposed theoretical model in numerical experiments with various shapes and voltages. Furthermore, we analyse theoretically the behaviour of the contact line at the onset of droplet motion as observed in the simulation and experiment, and we are able to estimate very well the time scale on which the contact angle changes.
The kinematic genesis of vortex formation due to finite rotation of a plate in still fluid
- M. Jimreeves David, Manikandan Mathur, R. N. Govardhan, J. H. Arakeri
-
- Published online by Cambridge University Press:
- 02 February 2018, pp. 489-524
-
- Article
- Export citation
-
We present a combined experimental and numerical study of an idealized model of the propulsive stroke of the turning manoeuvre in fish. Specifically, we use the framework of Lagrangian coherent structures (LCSs) to describe the kinematics of the flow that results from a thin plate performing a large angle rotation about its tip in still fluid. Temporally and spatially well-resolved velocity fields are obtained using a two-dimensional, incompressible finite-volume solver, and are validated by comparisons with experimentally measured velocity fields and alternate numerical simulations. We then implement the recently proposed variational theory of LCSs to extract the hyperbolic and elliptic LCSs in the numerically generated velocity fields. Detailed LCS analysis is performed for a plate motion profile described by $\dot{\unicode[STIX]{x1D703}}(t)=\unicode[STIX]{x1D6FA}_{max}\sin ^{2}(\unicode[STIX]{x1D714}t)$ during $0\leqslant t\leqslant t_{o}$ and zero otherwise. The stopping time $t_{o}$ is given by $t_{o}=\unicode[STIX]{x03C0}/\unicode[STIX]{x1D714}=10~\text{s}$, the value of $\unicode[STIX]{x1D6FA}_{max}$ chosen to give a stopping angle of $\unicode[STIX]{x1D703}_{max}=90^{\circ }$, resulting in a Reynolds number $Re=c^{2}\unicode[STIX]{x1D6FA}_{max}/\unicode[STIX]{x1D708}=785.4$, where $c$ is the plate chord length and $\unicode[STIX]{x1D708}=10^{-6}~\text{m}^{2}~\text{s}^{-1}$ the kinematic viscosity of water. The flow comprises a starting and a stopping vortex, resulting in a pair of oppositely signed vortices of unequal strengths that move away from the plate in a direction closely aligned with the final plate orientation at $t/t_{o}\approx 2$. The hyperbolic LCSs are shown to encompass the fluid material that is advected away from the plate for $t>t_{o}$, henceforth referred to as the advected bulk. The starting and stopping vortices, identified using elliptic LCSs and hence more objective than Eulerian vortex detection methods, constitute only around two thirds of the advected bulk area. The advected bulk is traced back to $t=0$ to identify five distinct lobes of fluid that eventually form the advected bulk, and hence map the long-term fate of various regions in the fluid at $t=0$. The five different lobes of fluid are then shown to be delineated by repelling LCS boundaries at $t=0$. The linear momentum of the advected bulk region is shown to account for approximately half of the total impulse experienced by the plate in the direction of its final orientation, thus establishing its dynamical significance. We provide direct experimental evidence for the kinematic relevance of hyperbolic and elliptic LCSs using novel dye visualization experiments, and also show that attracting hyperbolic LCSs provide objective characterization of the spiral structures often observed in vortical flows. We conclude by showing that qualitatively similar LCSs persist for several other plate motion profiles and stopping angles as well.
Impinging jet flow and hydraulic jump on a rotating disk
- Yunpeng Wang, Roger E. Khayat
-
- Published online by Cambridge University Press:
- 02 February 2018, pp. 525-560
-
- Article
- Export citation
-
The free-surface flow formed by a circular jet impinging on a rotating disk is analysed theoretically. The study explores the effects of rotation and inertia on the thin-film flow. Both boundary-layer height and film thickness are found to diminish with rotation speed. A maximum film thickness develops in the supercritical region, which reflects the competition between the convective and centrifugal effects. Unlike the flow on a stationary disk, an increase in the wall shear stress along the radial direction is predicted, at a rate that strengthens with rotating speed. Our results corroborate well existing measurements. The location and height of the hydraulic jump are determined subject to the value of the thickness at the edge of the disk, which is established first for a stationary disk based on the capillary length, and then for a rotating disk using existing analyses and measurements in spin coating. The case of a stationary is revisited in an effort to predict the location and height of the jump uniquely. The formulated value of the height at the edge of the disk seems to give excellent results for a jet at moderately high flow rate (or low viscosity) where the jump structure is well identifiable in reality.
Droplet impact onto an elastic plate: a new mechanism for splashing
- Michael Pegg, Richard Purvis, Alexander Korobkin
-
- Published online by Cambridge University Press:
- 02 February 2018, pp. 561-593
-
- Article
- Export citation
-
During a droplet impact onto a substrate, splashing is known to be caused by the presence of surrounding gas or by surface roughness. Impact occurring in a vacuum onto a smooth rigid wall results in droplet spreading, rather than development of a corona or prompt splash. Here we present an analytical and numerical study of a third potential splashing mechanism, namely elastic deformation of the substrate. An axisymmetric Wagner-style model of droplet impact is formulated and solved using the method of normal modes, together with asymptotic analysis and numerical methods. We highlight the effect that a flexible substrate brings to the contact line velocity and jet behaviour, demonstrating that oscillation of the substrate can cause blow-up of the splash jet which is absent for a rigid substrate and indicate the onset of splashing.
Small-scale dynamics of settling, bidisperse particles in turbulence
- Rohit Dhariwal, Andrew D. Bragg
-
- Published online by Cambridge University Press:
- 02 February 2018, pp. 594-620
-
- Article
- Export citation
-
Mixing and collisions of inertial particles at the small scales of turbulence can be investigated by considering how pairs of particles move relative to each other. In real problems the two particles will have different sizes, i.e. they are bidisperse, and the effect of gravity on their motion is often important. However, how turbulence and gravity compete to control the motion of bidisperse inertial particles is poorly understood. Motivated by this, we use direct numerical simulations (DNS) to investigate the dynamics of settling, bidisperse particles in isotropic turbulence. In agreement with previous studies, we find that without gravity (i.e. $Fr=\infty$, where $Fr$ is the Froude number), bidispersity leads to an enhancement of the relative velocities, and a suppression of their spatial clustering. For $Fr<1$, the relative velocities in the direction of gravity are enhanced by the differential settling velocities of the bidisperse particles, as expected. However, we also find that gravity can strongly enhance the relative velocities in the ‘horizontal’ directions (i.e. in the plane normal to gravity). This non-trivial behaviour occurs because fast settling particles experience rapid fluctuations in the fluid velocity field along their trajectory, leading to enhanced particle accelerations and relative velocities. Indeed, the results show that even when $Fr\ll 1$, turbulence can still play an important role, not only on the horizontal motion, but also on the vertical motion of the particles. This is related to the fact that $Fr$ only characterizes the importance of gravity compared with some typical acceleration of the fluid, yet accelerations in turbulence are highly intermittent. As a consequence, there is a significant probability for particles to be in regions of the flow where the Froude number based on the local, instantaneous fluid acceleration is ${>}1$, even though the typically defined Froude number is $\ll 1$. This could imply, for example, that extreme events in the mixing of settling, bidisperse particles are only weakly affected by gravity even when $Fr\ll 1$. We also find that gravity drastically reduces the clustering of bidisperse particles. These results are strikingly different to the monodisperse case, for which recent results have shown that when $Fr<1$, gravity strongly suppresses the relative velocities in all directions, and can enhance clustering. Finally, we consider the implications of these results for the collision rates of settling, bidisperse particles in turbulence. We find that for $Fr=0.052$, the collision kernel is almost perfectly predicted by the collision kernel for bidisperse particles settling in quiescent flow, such that the effect of turbulence may be ignored. However, for $Fr=0.3$, turbulence plays an important role, and the collisions are only dominated by gravitational settling when the difference in the particle Stokes numbers is ${\geqslant}O(1)$.