Papers
Two-dimensional numerical study of vortex shedding regimes of oscillatory flow past two circular cylinders in side-by-side and tandem arrangements at low Reynolds numbers
- Ming Zhao, Liang Cheng
-
- Published online by Cambridge University Press:
- 16 June 2014, pp. 1-37
-
- Article
- Export citation
-
Oscillatory flow past two circular cylinders in side-by-side and tandem arrangements at low Reynolds numbers is simulated numerically by solving the two-dimensional Navier–Stokes (NS) equations using a finite-element method (FEM). The aim of this study is to identify the flow regimes of the two-cylinder system at different gap arrangements and Keulegan–Carpenter numbers (KC). Simulations are conducted at seven gap ratios $\def \xmlpi #1{}\def \mathsfbi #1{\boldsymbol {\mathsf {#1}}}\let \le =\leqslant \let \leq =\leqslant \let \ge =\geqslant \let \geq =\geqslant \def \Pr {\mathit {Pr}}\def \Fr {\mathit {Fr}}\def \Rey {\mathit {Re}}G$ ($G=L/D$ where $L$ is the cylinder-to-cylinder gap and $D$ the diameter of a cylinder) of 0.5, 1, 1.5, 2, 3, 4 and 5 and KC ranging from 1 to 12 with an interval of 0.25. The flow regimes that have been identified for oscillatory flow around a single cylinder are also observed in the two-cylinder system but with different flow patterns due to the interactions between the two cylinders. In the side-by-side arrangement, the vortex shedding from the gap between the two cylinders dominates when the gap ratio is small, resulting in the gap vortex shedding (GVS) regime, which is different from any of the flow regimes identified for a single cylinder. For intermediate gap ratios of 1.5 and 2 in the side-by-side arrangement, the vortex shedding mode from one side of each cylinder is not necessarily the same as that from the other side, forming a so-called combined flow regime. When the gap ratio between the two cylinders is sufficiently large, the vortex shedding from each cylinder is similar to that of a single cylinder. In the tandem arrangement, when the gap between the two cylinders is very small, the flow regimes are similar to that of a single cylinder. For large gap ratios in the tandem arrangement, the vortex shedding flows from the gap side of the two cylinders interact and those from the outer sides of the cylinders are less affected by the existence of the other cylinder and similar to that of a single cylinder. Strong interaction between the vortex shedding flows from the two cylinders makes the flow very irregular at large KC values for both side-by-side and tandem arrangements.
On the structure and origin of pressure fluctuations in wall turbulence: predictions based on the resolvent analysis
- M. Luhar, A. S. Sharma, B. J. McKeon
-
- Published online by Cambridge University Press:
- 16 June 2014, pp. 38-70
-
- Article
- Export citation
-
We generate predictions for the fluctuating pressure field in turbulent pipe flow by reformulating the resolvent analysis of McKeon and Sharma (J. Fluid Mech., vol. 658, 2010, pp. 336–382) in terms of the so-called primitive variables. Under this analysis, the nonlinear convective terms in the Fourier-transformed Navier–Stokes equations (NSE) are treated as a forcing that is mapped to a velocity and pressure response by the resolvent of the linearized Navier–Stokes operator. At each wavenumber–frequency combination, the turbulent velocity and pressure field are represented by the most-amplified (rank-1) response modes, identified via a singular value decomposition of the resolvent. We show that these rank-1 response modes reconcile many of the key relationships among the velocity field, coherent structure (i.e. hairpin vortices), and the high-amplitude wall-pressure events observed in previous experiments and direct numerical simulations (DNS). A Green’s function representation shows that the pressure fields obtained under this analysis correspond primarily to the fast pressure contribution arising from the linear interaction between the mean shear and the turbulent wall-normal velocity. Recovering the slow pressure requires an explicit treatment of the nonlinear interactions between the Fourier response modes. By considering the velocity and pressure fields associated with the triadically consistent mode combination studied by Sharma and McKeon (J. Fluid Mech., vol. 728, 2013, pp. 196–238), we identify the possibility of an apparent amplitude modulation effect in the pressure field, similar to that observed for the streamwise velocity field. However, unlike the streamwise velocity, for which the large scales of the flow are in phase with the envelope of the small-scale activity close to the wall, we expect there to be a $\def \xmlpi #1{}\def \mathsfbi #1{\boldsymbol {\mathsf {#1}}}\let \le =\leqslant \let \leq =\leqslant \let \ge =\geqslant \let \geq =\geqslant \def \Pr {\mathit {Pr}}\def \Fr {\mathit {Fr}}\def \Rey {\mathit {Re}}\pi /2$ phase difference between the large-scale wall-pressure and the envelope of the small-scale activity. Finally, we generate spectral predictions based on a rank-1 model assuming broadband forcing across all wavenumber–frequency combinations. Despite the significant simplifying assumptions, this approach reproduces trends observed in previous DNS for the wavenumber spectra of velocity and pressure, and for the scale-dependence of wall-pressure propagation speed.
The role of advance ratio and aspect ratio in determining leading-edge vortex stability for flapping flight
- R. R. Harbig, J. Sheridan, M. C. Thompson
-
- Published online by Cambridge University Press:
- 16 June 2014, pp. 71-105
-
- Article
- Export citation
-
The effects of advance ratio and the wing’s aspect ratio on the structure of the leading-edge vortex (LEV) that forms on flapping and rotating wings under insect-like flight conditions are not well understood. However, recent studies have indicated that they could play a role in determining the stable attachment of the LEV. In this study, a numerical model of a flapping wing at insect Reynolds numbers is used to explore the effects of these parameters on the characteristics and stability of the LEV. The word ‘stability’ is used here to describe whether the LEV was attached throughout the stroke or if it was shed. It is demonstrated that increasing the advance ratio enhances vorticity production at the leading edge during the downstroke, and this results in more rapid growth of the LEV for non-zero advance ratios. Increasing the wing aspect ratio was found to have the effect of shortening the wing’s chord length relative to the LEV’s size. These two effects combined determine the stability of the LEV. For high advance ratios and large aspect ratios, the LEV was observed to quickly grow to envelop the entire wing during the early stages of the downstroke. Continued rotation of the wing resulted in the LEV being eventually shed as part of a vortex loop that peels away from the wing’s tip. The shedding of the LEV for high-aspect-ratio wings at non-zero advance ratios leads to reduced aerodynamic performance of these wings, which helps to explain why a number of insect species have evolved to have low-aspect-ratio wings.
Electrohydrodynamics of particle-covered drops
- Malika Ouriemi, Petia M. Vlahovska
-
- Published online by Cambridge University Press:
- 16 June 2014, pp. 106-120
-
- Article
- Export citation
-
We experimentally investigate the effect of surface-absorbed colloidal particles on the dynamics of a leaky dielectric drop in a uniform DC electric field. Depending on the particle polarizabilty, coverage and the electrical field intensity, particles assemble into various patterns such as an equatorial belt, pole-to-pole chains or a band of dynamic vortices. The particle structuring changes droplet electrohydrodynamics: under the same conditions where a particle-free drop would be a steady oblate spheroid, the belt can give rise to unsteady behaviours such as sustained drop wobbling or tumbling. Moreover, particle chaining can be accompanied by prolate drop deformation and tip-streaming.
The Burnett equations in cylindrical coordinates and their solution for flow in a microtube
- Narendra Singh, Amit Agrawal
-
- Published online by Cambridge University Press:
- 16 June 2014, pp. 121-141
-
- Article
- Export citation
-
The Burnett equations constitute a set of higher-order continuum equations. These equations are obtained from the Chapman–Enskog series solution of the Boltzmann equation while retaining second-order-accurate terms in the Knudsen number $\def \xmlpi #1{}\def \mathsfbi #1{\boldsymbol {\mathsf {#1}}}\let \le =\leqslant \let \leq =\leqslant \let \ge =\geqslant \let \geq =\geqslant \def \Pr {\mathit {Pr}}\def \Fr {\mathit {Fr}}\def \Rey {\mathit {Re}}\mathit{Kn}$. The set of higher-order continuum models is expected to be applicable to flows in the slip and transition regimes where the Navier–Stokes equations perform poorly. However, obtaining analytical or numerical solutions of these equations has been noted to be particularly difficult. In the first part of this work, we present the full set of Burnett equations in cylindrical coordinates in three-dimensional form. The equations are reported in a generalized way for gas molecules that are assumed to be Maxwellian molecules or hard spheres. In the second part, a closed-form solution of these equations for isothermal Poiseuille flow in a microtube is derived. The solution of the equations is shown to satisfy the full Burnett equations up to $\mathit{Kn} \leq 1.3$ within an error norm of ${\pm }1.0\, \%$. The mass flow rate obtained analytically is shown to compare well with available experimental and numerical results. Comparison of the stress terms in the Burnett and Navier–Stokes equations is presented. The significance of the Burnett normal stress and its role in diffusion of momentum is brought out by the analysis. An order-of-magnitude analysis of various terms in the equations is presented, based on which a reduced model of the Burnett equations is provided for flow in a microtube. The Burnett equations in full three-dimensional form in cylindrical coordinates and their solution are not previously available.
Effect of large bulk viscosity on large-Reynolds-number flows
- M. S. Cramer, F. Bahmani
-
- Published online by Cambridge University Press:
- 17 June 2014, pp. 142-163
-
- Article
- Export citation
-
We examine the inviscid and boundary-layer approximations in fluids having bulk viscosities which are large compared with their shear viscosities for three-dimensional steady flows over rigid bodies. We examine the first-order corrections to the classical lowest-order inviscid and laminar boundary-layer flows using the method of matched asymptotic expansions. It is shown that the effects of large bulk viscosity are non-negligible when the ratio of bulk to shear viscosity is of the order of the square root of the Reynolds number. The first-order outer flow is seen to be rotational, non-isentropic and viscous but nevertheless slips at the inner boundary. First-order corrections to the boundary-layer flow include a variation of the thermodynamic pressure across the boundary layer and terms interpreted as heat sources in the energy equation. The latter results are a generalization and verification of the predictions of Emanuel (Phys. Fluids A, vol. 4, 1992, pp. 491–495).
How flexibility affects the wake symmetry properties of a self-propelled plunging foil
- Xiaojue Zhu, Guowei He, Xing Zhang
-
- Published online by Cambridge University Press:
- 18 June 2014, pp. 164-183
-
- Article
- Export citation
-
The wake symmetry properties of a flapping-foil system are closely associated with its propulsive performance. In the present work, the effect of the foil flexibility on the wake symmetry properties of a self-propelled plunging foil is studied numerically. We compare the wakes of a flexible foil and a rigid foil at a low flapping Reynolds number of 200. The two foils are of the same dimensions, flapping frequency, leading-edge amplitude and cruising velocity but different bending rigidities. The results indicate that flexibility can either inhibit or trigger the symmetry breaking of the wake. We find that there exists a threshold value of vortex circulation above which symmetry breaking occurs. The modification of vortex circulation is found to be the pivotal factor in the influence of the foil flexibility on the wake symmetry properties. An increase in flexibility can result in a reduction in the vorticity production at the leading edge because of the decrease in the effective angle of attack, but it also enhances vorticity production at the trailing edge because of the increase in the trailing-edge flapping velocity. The competition between these two opposing effects eventually determines the strength of vortex circulation, which, in turn, governs the wake symmetry properties. Further investigation indicates that the former effect is related to the streamlined shape of the deformed foil while the latter effect is associated with structural resonance. The results of this work provide new insights into the functional role of passive flexibility in flapping-based biolocomotion.
Drops of power-law fluids falling on a coated vertical fibre
- Liyan Yu, John Hinch
-
- Published online by Cambridge University Press:
- 19 June 2014, pp. 184-215
-
- Article
- Export citation
-
We study the solitary wave solutions in a thin film of a power-law fluid coating a vertical fibre. Different behaviours are observed for shear-thickening and shear-thinning fluids. For shear-thickening fluids, the solitary waves are larger and faster when the reduced Bond number is smaller. For shear-thinning fluids, two branches of solutions exist for a certain range of the Bond number, where the solitary waves are larger and faster on one and smaller and slower on the other as the Bond number decreases. We carry out an asymptotic analysis for the large and fast-travelling solitary waves to explain how their speeds and amplitudes change with the Bond number. The analysis is then extended to examine the stability of the two branches of solutions for the shear-thinning fluids.
Quasi-geostrophic approximation of anelastic convection
- Friedrich H. Busse, Radostin D. Simitev
-
- Published online by Cambridge University Press:
- 18 June 2014, pp. 216-227
-
- Article
- Export citation
-
The onset of convection in a rotating cylindrical annulus with parallel ends filled with a compressible fluid is studied in the anelastic approximation. Thermal Rossby waves propagating in the azimuthal direction are found as solutions. The analogy to the case of Boussinesq convection in the presence of conical end surfaces of the annular region is emphasised. As in the latter case, the results can be applied as an approximation for the description of the onset of anelastic convection in rotating spherical fluid shells. Reasonable agreement with three-dimensional numerical results published by Jones, Kuzanyan & Mitchell (J. Fluid Mech., vol. 634, 2009, pp. 291–319) for the latter problem is found. As in those results, the location of the onset of convection shifts outwards from the tangent cylinder with increasing number $\def \xmlpi #1{}\def \mathsfbi #1{\boldsymbol {\mathsf {#1}}}\let \le =\leqslant \let \leq =\leqslant \let \ge =\geqslant \let \geq =\geqslant \def \Pr {\mathit {Pr}}\def \Fr {\mathit {Fr}}\def \Rey {\mathit {Re}}N_{\rho }$ of density scale heights until it reaches the equatorial boundary. A new result is that at a much higher number $N_{\rho }$ the onset location returns to the interior of the fluid shell.
The quiescent core of turbulent channel flow
- Y. S. Kwon, J. Philip, C. M. de Silva, N. Hutchins, J. P. Monty
-
- Published online by Cambridge University Press:
- 18 June 2014, pp. 228-254
-
- Article
- Export citation
-
The identification of uniform momentum zones in wall-turbulence, introduced by Adrian, Meinhart & Tomkins (J. Fluid Mech., vol. 422, 2000, pp. 1–54) has been applied to turbulent channel flow, revealing a large ‘core’ region having high and uniform velocity magnitude. Examination of the core reveals that it is a region of relatively weak turbulence levels. For channel flow in the range $\def \xmlpi #1{}\def \mathsfbi #1{\boldsymbol {\mathsf {#1}}}\let \le =\leqslant \let \leq =\leqslant \let \ge =\geqslant \let \geq =\geqslant \def \Pr {\mathit {Pr}}\def \Fr {\mathit {Fr}}\def \Rey {\mathit {Re}}Re_{\tau } = 1000\text {--}4000$, it was found that the ‘core’ is identifiable by regions bounded by the continuous isocontour lines of the streamwise velocity at $0.95U_{CL}$ (95 % of the centreline velocity). A detailed investigation into the properties of the core has revealed it has a large-scale oscillation which is predominantly anti-symmetric with respect to the channel centreline as it moves through the channel, and there is a distinct jump in turbulence statistics as the core boundary is crossed. It is concluded that the edge of the core demarcates a shear layer of relatively intense vorticity such that the interior of the core contains weakly varying, very low-level turbulence (relative to the flow closer to the wall). Although channel flows are generally referred to as ‘fully turbulent’, these findings suggest there exists a relatively large and ‘quiescent’ core region with a boundary qualitatively similar to the turbulent/non-turbulent interface of boundary layers, jets and wakes.
Inertial wave excitation and focusing in a liquid bounded by a frustum and a cylinder
- Marten Klein, Torsten Seelig, Michael V. Kurgansky, Abouzar Ghasemi V., Ion Dan Borcia, Andreas Will, Eberhard Schaller, Christoph Egbers, Uwe Harlander
-
- Published online by Cambridge University Press:
- 18 June 2014, pp. 255-297
-
- Article
- Export citation
-
The mechanism of localized inertial wave excitation and its efficiency is investigated for an annular cavity rotating with $\def \xmlpi #1{}\def \mathsfbi #1{\boldsymbol {\mathsf {#1}}}\let \le =\leqslant \let \leq =\leqslant \let \ge =\geqslant \let \geq =\geqslant \def \Pr {\mathit {Pr}}\def \Fr {\mathit {Fr}}\def \Rey {\mathit {Re}}\Omega _0$. Meridional symmetry is broken by replacing the inner cylinder with a truncated cone (frustum). Waves are excited by individual longitudinal libration of the walls. The geometry is non-separable and exhibits wave focusing and wave attractors. We investigated laboratory and numerical results for the Ekman number $E\approx 10^{-6}$, inclination $\alpha =5.71^\circ $ and libration amplitudes $\varepsilon \leq 0.2$ within the inertial wave band $0 < \omega < 2\Omega _0$. Under the assumption that the inertial waves do not essentially affect the boundary-layer structure, we use classical boundary-layer analysis to study oscillating Ekman layers over a librating wall that is at an angle $\alpha \neq 0$ to the axis of rotation. The Ekman layer erupts at frequency $\omega =f_{*}$, where $f_{*}\equiv 2 \Omega _0 \sin \alpha $ is the effective Coriolis parameter in a plane tangential to the wall. For the selected inclination this eruption occurs for the forcing frequency $\omega /\Omega _0=0.2$. For the librating lids eruption occurs at $\omega /\Omega _0=2$. The study reveals that the frequency dependence of the total kinetic energy $K_{\omega }$ of the excited wave field is strongly connected to the square of the Ekman pumping velocity $w_{{E}}(\omega )$ that, in the linear limit, becomes singular when the boundary layer erupts. This explains the frequency dependence of non-resonantly excited waves. By the localization of the forcing, the two configurations investigated, (i) frustum libration and (ii) lids together with outer cylinder in libration, can be clearly distinguished by their response spectra. Good agreement was found for the spatial structure of low-order wave attractors and periodic orbits (both characterized by a small number of reflections) in the frequency windows predicted by geometric ray tracing. For ‘resonant’ frequencies a significantly increased total bulk energy was found, while the energy in the boundary layer remained nearly constant. Inertial wave energy enters the bulk flow via corner beams, which are parallel to the characteristics of the underlying Poincaré problem. Numerical simulations revealed a mismatch between the wall-parallel mass fluxes near the corners. This leads to boundary-layer eruption and the generation of inertial waves in the corners.
Scaling of the turbulent/non-turbulent interface in boundary layers
- Kapil Chauhan, Jimmy Philip, Ivan Marusic
-
- Published online by Cambridge University Press:
- 19 June 2014, pp. 298-328
-
- Article
- Export citation
-
Scaling of the interface that demarcates a turbulent boundary layer from the non-turbulent free stream is sought using theoretical reasoning and experimental evidence in a zero-pressure-gradient boundary layer. The data-analysis, utilising particle image velocimetry (PIV) measurements at four different Reynolds numbers ($\def \xmlpi #1{}\def \mathsfbi #1{\boldsymbol {\mathsf {#1}}}\let \le =\leqslant \let \leq =\leqslant \let \ge =\geqslant \let \geq =\geqslant \def \Pr {\mathit {Pr}}\def \Fr {\mathit {Fr}}\def \Rey {\mathit {Re}}\delta u_{\tau }/\nu =1200\mbox{--}14\, 500$), indicates the presence of a viscosity dominated interface at all Reynolds numbers. It is found that the mean normal velocity across the interface and the tangential velocity jump scale with the skin-friction velocity $u_{\tau }$ and are approximately $u_{\tau }/10$ and $u_{\tau }$, respectively. The width of the superlayer is characterised by the local vorticity thickness $\delta _{\omega }$ and scales with the viscous length scale $\nu /u_{\tau }$. An order of magnitude analysis of the tangential momentum balance within the superlayer suggests that the turbulent motions also scale with inner velocity and length scales $u_{\tau }$ and $\nu /u_{\tau }$, respectively. The influence of the wall on the dynamics in the superlayer is considered via Townsend’s similarity hypothesis, which can be extended to account for the viscous influence at the turbulent/non-turbulent interface. Similar to a turbulent far-wake the turbulent motions in the superlayer are of the same order as the mean velocity deficit, which lends to a physical explanation for the emergence of the wake profile in the outer part of the boundary layer.
Reconnection of skewed vortices
- Y. Kimura, H. K. Moffatt
-
- Published online by Cambridge University Press:
- 20 June 2014, pp. 329-345
-
- Article
- Export citation
-
Based on experimental evidence that vortex reconnection commences with the approach of nearly antiparallel segments of vorticity, a linearised model is developed in which two Burgers-type vortices are driven together and stretched by an ambient irrotational strain field induced by more remote vorticity. When these Burgers vortices are exactly antiparallel, they are annihilated on the strain time-scale, independent of kinematic viscosity $\def \xmlpi #1{}\def \mathsfbi #1{\boldsymbol {\mathsf {#1}}}\let \le =\leqslant \let \leq =\leqslant \let \ge =\geqslant \let \geq =\geqslant \def \Pr {\mathit {Pr}}\def \Fr {\mathit {Fr}}\def \Rey {\mathit {Re}}\nu $ in the limit $\nu \rightarrow 0$. When the vortices are skew to each other, they are annihilated under this action over a local extent that increases exponentially in the stretching direction, with clear evidence of reconnection on the same strain time-scale. The initial helicity associated with the skewed geometry is eliminated during the process of reconnection. The model applies equally to the reconnection of weak magnetic flux tubes under the action of a strain field, when Lorentz forces are negligible.
Analysis of a model for foam improved oil recovery
- P. Grassia, E. Mas-Hernández, N. Shokri, S. J. Cox, G. Mishuris, W. R. Rossen
-
- Published online by Cambridge University Press:
- 20 June 2014, pp. 346-405
-
- Article
- Export citation
-
During improved oil recovery (IOR), gas may be introduced into a porous reservoir filled with surfactant solution in order to form foam. A model for the evolution of the resulting foam front known as ‘pressure-driven growth’ is analysed. An asymptotic solution of this model for long times is derived that shows that foam can propagate indefinitely into the reservoir without gravity override. Moreover, ‘pressure-driven growth’ is shown to correspond to a special case of the more general ‘viscous froth’ model. In particular, it is a singular limit of the viscous froth, corresponding to the elimination of a surface tension term, permitting sharp corners and kinks in the predicted shape of the front. Sharp corners tend to develop from concave regions of the front. The principal solution of interest has a convex front, however, so that although this solution itself has no sharp corners (except for some kinks that develop spuriously owing to errors in a numerical scheme), it is found nevertheless to exhibit milder singularities in front curvature, as the long-time asymptotic analytical solution makes clear. Numerical schemes for the evolving front shape which perform robustly (avoiding the development of spurious kinks) are also developed. Generalisations of this solution to geologically heterogeneous reservoirs should exhibit concavities and/or sharp corner singularities as an inherent part of their evolution: propagation of fronts containing such ‘inherent’ singularities can be readily incorporated into these numerical schemes.
Long-wave dynamics of an inextensible planar membrane in an electric field
- Y.-N. Young, Shravan Veerapaneni, Michael J. Miksis
-
- Published online by Cambridge University Press:
- 20 June 2014, pp. 406-431
-
- Article
- Export citation
-
In this paper the dynamics of an inextensible capacitive elastic membrane under an electric field is investigated in the long-wave (lubrication) leaky dielectric framework, where a sixth-order nonlinear differential equation with an integral constraint is derived. Steady equilibrium profiles for a non-conducting membrane in a direct current (DC) field are found to depend only on the membrane excess area and the volume under the membrane. Linear stability analysis on a tensionless flat membrane in a DC field gives the growth rate in terms of membrane conductance and electric properties in the bulk. Numerical simulations of a capacitive conducting membrane under an alternating current (AC) field elucidate how variation of the membrane tension correlates with the nonlinear membrane dynamics. Different membrane dynamics, such as undulation and flip-flop, are found at different electric field strength and membrane area. In particular a travelling wave on the membrane is found as a response to a periodic AC field in the perpendicular direction.
Perturbation theory and numerical modelling of weakly and moderately nonlinear dynamics of the incompressible Richtmyer–Meshkov instability
- A. L. Velikovich, M. Herrmann, S. I. Abarzhi
-
- Published online by Cambridge University Press:
- 23 June 2014, pp. 432-479
-
- Article
- Export citation
-
A study of incompressible two-dimensional (2D) Richtmyer–Meshkov instability (RMI) by means of high-order perturbation theory and numerical simulations is reported. Nonlinear corrections to Richtmyer’s impulsive formula for the RMI bubble and spike growth rates have been calculated for arbitrary Atwood number and an explicit formula has been obtained for it in the Boussinesq limit. Conditions for early-time acceleration and deceleration of the bubble and the spike have been elucidated. Theoretical time histories of the interface curvature at the bubble and spike tip and the profiles of vertical and horizontal velocities have been calculated and favourably compared to simulation results. In our simulations we have solved 2D unsteady Navier–Stokes equations for immiscible incompressible fluids using the finite volume fractional step flow solver NGA developed by Desjardins et al. (J. Comput. Phys., vol. 227, 2008, pp. 7125–7159) coupled to the level set based interface solver LIT (Herrmann, J. Comput. Phys., vol. 227, 2008, pp. 2674–2706). We study the impact of small amounts of viscosity on the flow dynamics and compare simulation results to theory to discuss the influence of the theory’s ideal inviscid flow assumption.
The coalescence of liquid drops in a viscous fluid: interface formation model
- James E. Sprittles, Yulii D. Shikhmurzaev
-
- Published online by Cambridge University Press:
- 24 June 2014, pp. 480-499
-
- Article
- Export citation
-
The interface formation model is applied to describe the initial stages of the coalescence of two liquid drops in the presence of a viscous ambient fluid whose dynamics is fully accounted for. Our focus is on understanding (a) how this model’s predictions differ from those of the conventionally used one, (b) what influence the ambient fluid has on the evolution of the shape of the coalescing drops and (c) the coupling of the intrinsic dynamics of coalescence and that of the ambient fluid. The key feature of the interface formation model in its application to the coalescence phenomenon is that it removes the singularity inherent in the conventional model at the onset of coalescence and describes the part of the free surface ‘trapped’ between the coalescing volumes as they are pressed against each other as a rapidly disappearing ‘internal interface’. Considering the simplest possible formulation of this model, we find experimentally verifiable differences with the predictions of the conventional model showing, in particular, the effect of drop size on the coalescence process. According to the new model, for small drops a non-monotonic time dependence of the bridge expansion speed is a feature that could be looked for in further experimental studies. Finally, the results of both models are compared to recently available experimental data on the evolution of the liquid bridge connecting coalescing drops, and the interface formation model is seen to give a better agreement with the data.
Discrete-vortex method with novel shedding criterion for unsteady aerofoil flows with intermittent leading-edge vortex shedding
- Kiran Ramesh, Ashok Gopalarathnam, Kenneth Granlund, Michael V. Ol, Jack R. Edwards
-
- Published online by Cambridge University Press:
- 23 June 2014, pp. 500-538
-
- Article
- Export citation
-
Unsteady aerofoil flows are often characterized by leading-edge vortex (LEV) shedding. While experiments and high-order computations have contributed to our understanding of these flows, fast low-order methods are needed for engineering tasks. Classical unsteady aerofoil theories are limited to small amplitudes and attached leading-edge flows. Discrete-vortex methods that model vortex shedding from leading edges assume continuous shedding, valid only for sharp leading edges, or shedding governed by ad-hoc criteria such as a critical angle of attack, valid only for a restricted set of kinematics. We present a criterion for intermittent vortex shedding from rounded leading edges that is governed by a maximum allowable leading-edge suction. We show that, when using unsteady thin aerofoil theory, this leading-edge suction parameter (LESP) is related to the $\def \xmlpi #1{}\def \mathsfbi #1{\boldsymbol {\mathsf {#1}}}\let \le =\leqslant \let \leq =\leqslant \let \ge =\geqslant \let \geq =\geqslant \def \Pr {\mathit {Pr}}\def \Fr {\mathit {Fr}}\def \Rey {\mathit {Re}}A_0$ term in the Fourier series representing the chordwise variation of bound vorticity. Furthermore, for any aerofoil and Reynolds number, there is a critical value of the LESP, which is independent of the motion kinematics. When the instantaneous LESP value exceeds the critical value, vortex shedding occurs at the leading edge. We have augmented a discrete-time, arbitrary-motion, unsteady thin aerofoil theory with discrete-vortex shedding from the leading edge governed by the instantaneous LESP. Thus, the use of a single empirical parameter, the critical-LESP value, allows us to determine the onset, growth, and termination of LEVs. We show, by comparison with experimental and computational results for several aerofoils, motions and Reynolds numbers, that this computationally inexpensive method is successful in predicting the complex flows and forces resulting from intermittent LEV shedding, thus validating the LESP concept.
Critical layer and radiative instabilities in shallow-water shear flows
- Xavier Riedinger, Andrew D. Gilbert
-
- Published online by Cambridge University Press:
- 23 June 2014, pp. 539-569
-
- Article
- Export citation
-
In this study a linear stability analysis of shallow-water flows is undertaken for a representative Froude number $\def \xmlpi #1{}\def \mathsfbi #1{\boldsymbol {\mathsf {#1}}}\let \le =\leqslant \let \leq =\leqslant \let \ge =\geqslant \let \geq =\geqslant \def \Pr {\mathit {Pr}}\def \Fr {\mathit {Fr}}\def \Rey {\mathit {Re}}F=3.5$. The focus is on monotonic base flow profiles $U$ without an inflection point, in order to study critical layer instability (CLI) and its interaction with radiative instability (RI). First the dispersion relation is presented for the piecewise linear profile studied numerically by Satomura (J. Meterol. Soc. Japan, vol. 59, 1981, pp. 148–167) and using WKBJ analysis an interpretation given of mode branches, resonances and radiative instability. In particular surface gravity (SG) waves can resonate with a limit mode (LM) (or Rayleigh wave), localised near the discontinuity in shear in the flow; in this piecewise profile there is no critical layer. The piecewise linear profile is then continuously modified in a family of nonlinear profiles, to show the effect of the vorticity gradient $Q^{\prime } = - U^{\prime \prime }$ on the nature of the modes. Some modes remain as modes and others turn into quasi-modes (QM), linked to Landau damping of disturbances to the flow, depending on the sign of the vorticity gradient at the critical point. Thus an interpretation of critical layer instability for continuous profiles is given, as the remnant of the resonance with the LM. Numerical results and WKBJ analysis of critical layer instability and radiative instability for more general smooth profiles are provided. A link is made between growth rate formulae obtained by considering wave momentum and those found via the WKBJ approximation. Finally the competition between the stabilising effect of vorticity gradients in a critical layer and the destabilising effect of radiation (radiative instability) is studied.
Low-Reynolds-number wakes of elliptical cylinders: from the circular cylinder to the normal flat plate
- Mark C. Thompson, Alexander Radi, Anirudh Rao, John Sheridan, Kerry Hourigan
-
- Published online by Cambridge University Press:
- 24 June 2014, pp. 570-600
-
- Article
- Export citation
-
While the wake of a circular cylinder and, to a lesser extent, the normal flat plate have been studied in considerable detail, the wakes of elliptic cylinders have not received similar attention. However, the wakes from the first two bodies have considerably different characteristics, in terms of three-dimensional transition modes, and near- and far-wake structure. This paper focuses on elliptic cylinders, which span these two disparate cases. The Strouhal number and drag coefficient variations with Reynolds number are documented for the two-dimensional shedding regime. There are considerable differences from the standard circular cylinder curve. The different three-dimensional transition modes are also examined using Floquet stability analysis based on computed two-dimensional periodic base flows. As the cylinder aspect ratio (major to minor axis) is decreased, mode A is no longer unstable for aspect ratios below 0.25, as the wake deviates further from the standard Bénard–von Kármán state. For still smaller aspect ratios, another three-dimensional quasi-periodic mode becomes unstable, leading to a different transition scenario. Interestingly, for the 0.25 aspect ratio case, mode A restabilises above a Reynolds number of approximately 125, allowing the wake to return to a two-dimensional state, at least in the near wake. For the flat plate, three-dimensional simulations show that the shift in the Strouhal number from the two-dimensional value is gradual with Reynolds number, unlike the situation for the circular cylinder wake once mode A shedding develops. Dynamic mode decomposition is used to characterise the spatially evolving character of the wake as it undergoes transition from the primary Bénard–von Kármán-like near wake into a two-layered wake, through to a secondary Bénard–von Kármán-like wake further downstream, which in turn develops an even longer wavelength unsteadiness. It is also used to examine the differences in the two- and three-dimensional near-wake state, showing the increasing distortion of the two-dimensional rollers as the Reynolds number is increased.