Papers
Experimental study of parametric subharmonic instability for internal plane waves
- Baptiste Bourget, Thierry Dauxois, Sylvain Joubaud, Philippe Odier
-
- Published online by Cambridge University Press:
- 16 April 2013, pp. 1-20
-
- Article
- Export citation
-
Internal waves are believed to be of primary importance as they affect ocean mixing and energy transport. Several processes can lead to the breaking of internal waves and they usually involve nonlinear interactions between waves. In this work, we study experimentally the parametric subharmonic instability (PSI), which provides an efficient mechanism to transfer energy from large to smaller scales. It corresponds to the destabilization of a primary plane wave and the spontaneous emission of two secondary waves, of lower frequencies and different wave vectors. Using a time-frequency analysis, we observe the time evolution of the secondary waves, thus measuring the growth rate of the instability. In addition, a Hilbert transform method allows the measurement of the different wave vectors. We compare these measurements with theoretical predictions, and study the dependence of the instability with primary wave frequency and amplitude, revealing a possible effect of the confinement due to the finite size of the beam, on the selection of the unstable mode.
Low-Reynolds-number flow through two-dimensional shunts
- A. Setchi, A. J. Mestel, K. H. Parker, J. H. Siggers
-
- Published online by Cambridge University Press:
- 16 April 2013, pp. 21-39
-
- Article
- Export citation
-
Motivated by numerous biological and industrial applications relating to bypasses, mixing and leakage, we consider low-Reynolds-number flow through a shunt between two channels. An analytical solution for the streamfunction is found by matching biorthogonal expansions of Papkovich–Fadle eigenfunctions in rectangular subregions. The general solution can be adapted to model a variety of interesting problems of flow through two-dimensional shunts by imposing different inlet and outlet flux distributions. We present several such flow profiles but the majority of results relate to the particular problem of a side-to-side anastomosis in the small intestine. We consider different flux fractions through the shunt with particular emphasis on the pressure and recirculating regions, which are important factors in estimating health risks pertaining to this surgical procedure.
Quasi-geostrophic shallow-water doubly-connected vortex equilibria and their stability
- Hanna Płotka, David G. Dritschel
-
- Published online by Cambridge University Press:
- 16 April 2013, pp. 40-68
-
- Article
- Export citation
-
We examine the form, properties, stability and evolution of doubly-connected (two-vortex) relative equilibria in the single-layer $f$-plane quasi-geostrophic shallow-water model of geophysical fluid dynamics. Three parameters completely describe families of equilibria in this system: the ratio $\gamma = L/ {L}_{D} $ between the horizontal size of the vortices and the Rossby deformation length; the area ratio $\alpha $ of the smaller to the larger vortex; and the minimum distance $\delta $ between the two vortices. We vary $0\lt \gamma \leq 10$ and $0. 1\leq \alpha \leq 1. 0$, determining the boundary of stability $\delta = {\delta }_{c} (\gamma , \alpha )$. We also examine the nonlinear development of the instabilities and the transitions to other near-equilibrium configurations. Two modes of instability occur when $\delta \lt {\delta }_{c} $: a small-$\gamma $ asymmetric (wave 3) mode, which is absent for $\alpha \gtrsim 0. 6$; and a large-$\gamma $ mode. In general, major structural changes take place during the nonlinear evolution of the vortices, which near ${\delta }_{c} $ may be classified as follows: (i) vacillations about equilibrium for $\gamma \gtrsim 2. 5$; (ii) partial straining out, associated with the small-$\gamma $ mode, where either one or both of the vortices get smaller for $\gamma \lesssim 2. 5$ and $\alpha \lesssim 0. 6$; (iii) partial merger, occurring at the transition region between the two modes of instability, where one of the vortices gets bigger, and (iv) complete merger, associated with the large-$\gamma $ mode. We also find that although conservative inviscid transitions to equilibria with the same energy, angular momentum and circulation are possible, they are not the preferred evolutionary path.
Coarsening and solidification via solvent-annealing in thin liquid films
- Tony S. Yu, Vladimir Bulović, A. E. Hosoi
-
- Published online by Cambridge University Press:
- 16 April 2013, pp. 69-90
-
- Article
- Export citation
-
We examine solidification in thin liquid films produced by annealing amorphous ${\mathrm{Alq} }_{3} $ (tris-(8-hydroxyquinoline) aluminium) in methanol vapour. Micrographs acquired during annealing capture the evolution of the film: the initially-uniform film breaks up into drops that coarsen, and single crystals of ${\mathrm{Alq} }_{3} $ nucleate randomly on the substrate and grow as slender ‘needles’. The growth of these needles appears to follow power-law behaviour, where the growth exponent, $\gamma $, depends on the thickness of the deposited ${\mathrm{Alq} }_{3} $ film. The evolution of the thin film is modelled by a lubrication equation, and an advection–diffusion equation captures the transport of ${\mathrm{Alq} }_{3} $ and methanol within the film. We define a dimensionless transport parameter, $\alpha $, which is analogous to an inverse Sherwood number and quantifies the relative effects of diffusion- and coarsening-driven advection. For large $\alpha $-values, the model recovers the theory of one-dimensional, diffusion-driven solidification, such that $\gamma \rightarrow 1/ 2$. For low $\alpha $-values, the collapse of drops, i.e. coarsening, drives flow and regulates the growth of needles. Within this regime, we identify two relevant limits: needles that are small compared to the typical drop size, and those that are large. Both scaling analysis and simulations of the full model reveal that $\gamma \rightarrow 2/ 5$ for small needles and $\gamma \rightarrow 0. 29$ for large needles.
An explicit algebraic Reynolds-stress and scalar-flux model for stably stratified flows
- W. M. J. Lazeroms, G. Brethouwer, S. Wallin, A. V. Johansson
-
- Published online by Cambridge University Press:
- 16 April 2013, pp. 91-125
-
- Article
- Export citation
-
This work describes the derivation of an algebraic model for the Reynolds stresses and turbulent heat flux in stably stratified turbulent flows, which are mutually coupled for this type of flow. For general two-dimensional mean flows, we present a correct way of expressing the Reynolds-stress anisotropy and the (normalized) turbulent heat flux as tensorial combinations of the mean strain rate, the mean rotation rate, the mean temperature gradient and gravity. A system of linear equations is derived for the coefficients in these expansions, which can easily be solved with computer algebra software for a specific choice of the model constants. The general model is simplified in the case of parallel mean shear flows where the temperature gradient is aligned with gravity. For this case, fully explicit and coupled expressions for the Reynolds-stress tensor and heat-flux vector are given. A self-consistent derivation of this model would, however, require finding a root of a polynomial equation of sixth-order, for which no simple analytical expression exists. Therefore, the nonlinear part of the algebraic equations is modelled through an approximation that is close to the consistent formulation. By using the framework of a $K\text{{\ndash}} \omega $ model (where $K$ is turbulent kinetic energy and $\omega $ an inverse time scale) and, where needed, near-wall corrections, the model is applied to homogeneous shear flow and turbulent channel flow, both with stable stratification. For the case of homogeneous shear flow, the model predicts a critical Richardson number of 0.25 above which the turbulent kinetic energy decays to zero. The channel-flow results agree well with DNS data. Furthermore, the model is shown to be robust and approximately self-consistent. It also fulfils the requirements of realizability.
Some insights for the prediction of near-wall turbulence
- Farid Karimpour, Subhas K. Venayagamoorthy
-
- Published online by Cambridge University Press:
- 16 April 2013, pp. 126-139
-
- Article
- Export citation
-
In this paper, we revisit the eddy viscosity formulation to highlight a number of important issues that have direct implications for the prediction of near-wall turbulence. For steady wall-bounded turbulent flows, we make the equilibrium assumption between rates of production ($P$) and dissipation ($\epsilon $) of turbulent kinetic energy ($k$) in the near-wall region to propose that the eddy viscosity should be given by ${\nu }_{t} \approx \epsilon / {S}^{2} $, where $S$ is the mean shear rate. We then argue that the appropriate velocity scale is given by $\mathop{(S{T}_{L} )}\nolimits ^{- 1/ 2} {k}^{1/ 2} $ where ${T}_{L} = k/ \epsilon $ is the turbulence (decay) time scale. The difference between this velocity scale and the commonly assumed velocity scale of ${k}^{1/ 2} $ is subtle but the consequences are significant for near-wall effects. We then extend our discussion to show that the fundamental length and time scales that capture the near-wall behaviour in wall-bounded shear flows are the shear mixing length scale ${L}_{S} = \mathop{(\epsilon / {S}^{3} )}\nolimits ^{1/ 2} $ and the mean shear time scale $1/ S$, respectively. With these appropriate length and time scales (or equivalently velocity and time scales), the eddy viscosity can be rewritten in the familiar form of the $k$–$\epsilon $ model as ${\nu }_{t} = \mathop{(1/ S{T}_{L} )}\nolimits ^{2} {k}^{2} / \epsilon $. We use the direct numerical simulation (DNS) data of turbulent channel flow of Hoyas & Jiménez (Phys. Fluids, vol. 18, 2006, 011702) and the turbulent boundary layer flow of Jiménez et al. (J. Fluid Mech. vol. 657, 2010, pp. 335–360) to perform ‘a priori’ tests to check the validity of the revised eddy viscosity formulation. The comparisons with the exact computations from the DNS data are remarkable and highlight how well the equilibrium assumption holds in the near-wall region. These findings could prove to be useful in near-wall modelling of turbulent flows.
Lagrangian approach to laminar–turbulent interfaces in transitional pipe flow
- Markus Holzner, Baofang Song, Marc Avila, Björn Hof
-
- Published online by Cambridge University Press:
- 16 April 2013, pp. 140-162
-
- Article
- Export citation
-
Transition in shear flows is characterized by localized turbulent regions embedded in the surrounding laminar flow. These so-called turbulent spots or puffs are observed in a variety of shear flows and in certain Reynolds-number regimes, and they are advected by the flow while keeping their characteristic length. We show here for the case of pipe flow that this seemingly passive advection of turbulent puffs involves continuous entrainment and relaminarization of laminar and turbulent fluid across strongly convoluted interfaces. Surprisingly, interface areas are almost two orders of magnitude larger than the pipe cross-section, while local entrainment velocities are much smaller than the mean speed. Even though these velocities were shown to be small and proportional to the Kolmogorov velocity scale (in agreement with a prediction by Corrsin) in a flow without mean shear before, we find that, in pipe flow, local entrainment velocities are about an order of magnitude smaller than this scale. The Lagrangian method used to study the dynamics of the laminar–turbulent interfaces allows accurate determination of the leading and trailing edge speeds. However, to resolve the highly complex interface dynamics requires much higher numerical resolutions than for ordinary turbulent flows. This method also reveals that the volume flux across the leading edge has the same radial dependence but the opposite sign as that across the trailing edge, and it is this symmetry that is responsible for the puff shape remaining constant.
A spherical squirming swimmer in unsteady Stokes flow
- Kenta Ishimoto
-
- Published online by Cambridge University Press:
- 16 April 2013, pp. 163-189
-
- Article
- Export citation
-
The motion of a spherical squirmer in unsteady Stokes flow is investigated for a deeper understanding of unsteady inertial effects on swimming micro-organisms and differences of swimming strokes between a wave pattern and a flapping motion. An asymptotic analysis with respect to the small amplitude and the small inertia is performed, and the average swimming velocity after a long period of time under an assumption of a time-periodic stroke is obtained. This analysis shows that the scallop theorem still holds in a long-time asymptotic sense for tangential deformation, but that the time variation of the shape generates a net velocity even for a reciprocal swimmer. It is also found that the inertial effects on the swimming velocity are significant for a flapping swimmer, as contrasted with little influence on that of a swimmer with a wave pattern. The inertial effect is also illustrated with a simple squirmer, so that a reciprocal motion can be almost an optimal stroke under a constraint on energy consumption.
Solution of the quasi-one-dimensional linearized Euler equations using flow invariants and the Magnus expansion
- Ignacio Duran, Stephane Moreau
-
- Published online by Cambridge University Press:
- 16 April 2013, pp. 190-231
-
- Article
- Export citation
-
The acoustic and entropy transfer functions of quasi-one-dimensional nozzles are studied analytically for both subsonic and choked flows with and without shock waves. The present analytical study extends both the compact nozzle solution obtained by Marble & Candel (J. Sound Vib., vol. 55, 1977, pp. 225–243) and the effective nozzle length proposed by Stow, Dowling & Hynes (J. Fluid Mech., vol. 467, 2002, pp. 215–239) and by Goh & Morgans (J. Sound Vib., vol. 330, 2011, pp. 5184–5198) to non-zero frequencies for both modulus and phase through an asymptotic expansion of the linearized Euler equations. It also extends the piecewise-linear approximation of the velocity profile in the nozzle proposed by Moase, Brear & Manzie (J. Fluid Mech., vol. 585, 2007, pp. 281–304) to any arbitrary profile or equivalently any nozzle geometry. The equations are written as a function of three variables, namely the dimensionless mass, total temperature and entropy fluctuations, yielding a first-order linear system of differential equations with varying coefficients, which is solved using the Magnus expansion. The solution shows that both the modulus and the phase of the transfer functions of the nozzle have a strong dependence on the frequency. This holds for both choked flows and subsonic converging–diverging nozzles. The method is used to compare two different nozzle geometries with the same inlet and outlet Mach numbers, showing that, even if the compact solution predicts no differences between the transfer functions of the two nozzles, significant differences are found at non-zero frequencies. A parametric study is finally performed to calculate the indirect to direct noise ratio for a model combustor, showing that this ratio decreases at higher frequencies.
Worst-case amplification of disturbances in inertialess Couette flow of viscoelastic fluids
- Binh K. Lieu, Mihailo R. Jovanović, Satish Kumar
-
- Published online by Cambridge University Press:
- 16 April 2013, pp. 232-263
-
- Article
- Export citation
-
Amplification of deterministic disturbances in inertialess shear-driven channel flows of viscoelastic fluids is examined by analysing the frequency responses from spatio-temporal body forces to the velocity and polymer stress fluctuations. In strongly elastic flows, we show that disturbances with large streamwise length scales may be significantly amplified even in the absence of inertia. For fluctuations without streamwise variations, we derive explicit analytical expressions for the dependence of the worst-case amplification (from different forcing to different velocity and polymer stress components) on the Weissenberg number ($\mathit{We}$), the maximum extensibility of the polymer chains ($L$), the viscosity ratio and the spanwise wavenumber. For the Oldroyd-B model, the amplification of the most energetic components of velocity and polymer stress fields scales as ${\mathit{We}}^{2} $ and ${\mathit{We}}^{4} $. On the other hand, the finite extensibility of polymer molecules limits the largest achievable amplification even in flows with infinitely large Weissenberg numbers: in the presence of wall-normal and spanwise forces, the amplification of the streamwise velocity and polymer stress fluctuations is bounded by quadratic and quartic functions of $L$. This high amplification signals low robustness to modelling imperfections of inertialess channel flows of viscoelastic fluids. The underlying physical mechanism involves interactions of polymer stress fluctuations with a base shear, and it represents a close analogue of the lift-up mechanism that initiates a bypass transition in inertial flows of Newtonian fluids.
Near-wall turbulent fluctuations in the absence of wide outer motions
- Yongyun Hwang
-
- Published online by Cambridge University Press:
- 16 April 2013, pp. 264-288
-
- Article
- Export citation
-
Numerical experiments that remove turbulent motions wider than ${ \lambda }_{z}^{+ } \simeq 100$ are carried out up to ${\mathit{Re}}_{\tau } = 660$ in a turbulent channel. The artificial removal of the wide outer turbulence is conducted with spanwise minimal computational domains and an explicit filter that effectively removes spanwise uniform eddies. The mean velocity profile of the remaining motions shows very good agreement with that of the full simulation below ${y}^{+ } \simeq 40$, and the near-wall peaks of the streamwise velocity fluctuation scale very well in the inner units and remain almost constant at all the Reynolds numbers considered. The self-sustaining motions narrower than ${ \lambda }_{z}^{+ } \simeq 100$ generate smaller turbulent skin friction than full turbulent motions, and their contribution to turbulent skin friction gradually decays with the Reynolds number. This finding suggests that the role of the removed outer structures becomes increasingly important with the Reynolds number; thus one should aim to control the large scales for turbulent drag reduction at high Reynolds numbers. In the near-wall region, the streamwise and spanwise velocity fluctuations of the motions of ${ \lambda }_{z}^{+ } \leq 100$ reveal significant lack of energy at long streamwise lengths compared to those of the full simulation. In contrast, the losses of the wall-normal velocity and the Reynolds stress are not as large as those of these two variables. This implies that the streamwise and spanwise velocities of the removed motions penetrate deep into the near-wall region, while the wall-normal velocity and the Reynolds stress do not.
Multilayer shallow water equations with complete Coriolis force. Part 3. Hyperbolicity and stability under shear
- Andrew L. Stewart, Paul J. Dellar
-
- Published online by Cambridge University Press:
- 16 April 2013, pp. 289-317
-
- Article
- Export citation
-
We analyse the hyperbolicity of our multilayer shallow water equations that include the complete Coriolis force due to the Earth’s rotation. Shallow water theory represents flows in which the vertical shear is concentrated into vortex sheets between layers of uniform velocity. Such configurations are subject to Kelvin–Helmholtz instabilities, with arbitrarily large growth rates for sufficiently short-wavelength disturbances. These instabilities manifest themselves through a loss of hyperbolicity in the shallow water equations, rendering them ill-posed for the solution of initial value problems. We show that, in the limit of vanishingly small density difference between the two layers, our two-layer shallow water equations remain hyperbolic when the velocity difference remains below the same threshold that also ensures the hyperbolicity of the standard shallow water equations. Direct calculation of the domain of hyperbolicity becomes much less tractable for three or more layers, so we demonstrate numerically that the threshold for the velocity differences, below which the three-layer equations remain hyperbolic, is also unchanged by the inclusion of the complete Coriolis force. In all cases, the shape of the domain of hyperbolicity, which extends outside the threshold, changes considerably. The standard shallow water equations only lose hyperbolicity due to shear parallel to the direction of wave propagation, but the complete Coriolis force introduces another mechanism for loss of hyperbolicity due to shear in the perpendicular direction. We demonstrate that this additional mechanism corresponds to the onset of a transverse shear instability driven by the non-traditional components of the Coriolis force in a three-dimensional continuously stratified fluid.
Turbulent boundary-layer noise: direct radiation at Mach number 0.5
- Xavier Gloerfelt, Julien Berland
-
- Published online by Cambridge University Press:
- 16 April 2013, pp. 318-351
-
- Article
- Export citation
-
Boundary layers constitute a fundamental source of aerodynamic noise. A turbulent boundary layer over a plane wall can provide an indirect contribution to the noise by exciting the structure and a direct noise contribution. The latter part can play a significant role even if its intensity is very low, explaining why it is difficult to measure. In the present study, the aerodynamic noise generated by a spatially developing turbulent boundary layer is computed directly by solving the compressible Navier–Stokes equations. This numerical experiment aims at giving some insight into the noise radiation characteristics. The acoustic wavefronts have a large wavelength and are oriented in the direction opposite to the flow. Their amplitude is only 0.7 % of the aerodynamic pressure for a flat-plate flow at Mach 0.5. The particular directivity is mainly explained by convection effects by the mean flow, giving an indication about the compactness of the sources. These vortical events correspond to low frequencies and thus have a large lifetime. They cannot be directly associated with the main structures populating the boundary layer such as hairpin or horseshoe vortices. The analysis of the wall pressure can provide a picture of the flow in the wavenumber–frequency space. The main features of wall pressure beneath a turbulent boundary layer as described in the literature are well reproduced. The acoustic domain, corresponding to supersonic wavenumbers, is detectable but can hardly be separated from the convective ridge at this relatively high speed. This is also due to the low frequencies of sound emission as noted previously.
Diffusion in hydrogel-supported phospholipid bilayer membranes
- Chih-Ying Wang, Reghan J. Hill
-
- Published online by Cambridge University Press:
- 16 April 2013, pp. 352-373
-
- Article
- Export citation
-
We model a cylindrical inclusion (lipid or membrane protein) translating with velocity $U$ in a thin planar membrane (phospholipid bilayer) that is supported above and below by Brinkman media (hydrogels). The total force $F$, membrane velocity, and solvent velocity are calculated as functions of three independent dimensionless parameters: $\Lambda = \eta a/ ({\eta }_{m} h)$, ${\ell }_{1} / a$ and ${\ell }_{2} / a$. Here, $\eta $ and ${\eta }_{m} $ are the solvent and membrane shear viscosities, $a$ is the particle radius, $h$ is the membrane thickness, and ${ \ell }_{1}^{2} $ and ${ \ell }_{2}^{2} $ are the upper and lower hydrogel permeabilities. As expected, the dimensionless mobility $4\mathrm{\pi} \eta aU/ F= 4\mathrm{\pi} \eta aD/ ({k}_{B} T)$ (proportional to the self-diffusion coefficient, $D$) decreases with decreasing gel permeabilities (increasing gel concentrations), furnishing a quantitative interpretation of how porous, gel-like supports hinder membrane dynamics. The model also provides a means of inferring hydrogel permeability and, perhaps, surface morphology from tracer diffusion measurements.
The quasi-geostrophic theory of the thermal shallow water equations
- Emma S. Warneford, Paul J. Dellar
-
- Published online by Cambridge University Press:
- 16 April 2013, pp. 374-403
-
- Article
- Export citation
-
The thermal shallow water equations provide a depth-averaged description of motions in a fluid layer that permits horizontal variations in material properties. They typically arise through an equivalent barotropic approximation of a two-layer system, with a spatially varying density contrast due to an evolving temperature field in the active layer. We formalize a previous derivation of the quasi-geostrophic (QG) theory of these equations, by performing a direct asymptotic expansion for small Rossby number. We then present a second derivation as the small Rossby number limit of a balanced model that projects out high-frequency dynamics due to inertia-gravity waves. This latter derivation has wider validity, not being restricted to mid-latitude $\beta $-planes. We also derive their local energy conservation equation from the QG limit of a thermal shallow water pseudo-energy conservation equation. This derivation involves the ageostrophic correction to the leading-order geostrophic velocity that is eliminated in the usual derivation of a closed evolution equation for the QG potential vorticity. Finally, we derive the non-canonical Hamiltonian structure of the thermal QG equations from a decomposition in Rossby number of a pseudo-energy and Poisson bracket for the thermal shallow water equations.
Large-eddy simulation of passive scalar dispersion in an urban-like canopy
- D. A. Philips, R. Rossi, G. Iaccarino
-
- Published online by Cambridge University Press:
- 16 April 2013, pp. 404-428
-
- Article
- Export citation
-
Results from large-eddy simulations of short-range dispersion of a passive scalar from a point source release in an urban-like canopy are presented. The computational domain is that of a variable height array of buildings immersed in a pressure-driven, turbulent flow with a roughness Reynolds number ${\mathit{Re}}_{\tau } = 433$. A comparative study of several cases shows the changes in plume behaviour for different mean flow directions and source locations. The analysis of the results focuses on utilizing the high-fidelity datasets to examine the three-dimensional flow field and scalar plume structure. The detailed solution of the flow and scalar fields within the canopy allows for a direct assessment of the impact of local features of the building array geometry. The staggered, skewed and aligned arrangements of the buildings with respect to the oncoming flow were shown to affect plume development. Additional post-processing quantified this development through parameters fundamental to reduced-order Gaussian dispersion models. The parameters include measures of concentration decay with distance from the source as well as plume trajectory and spread. The horizontal plume trajectory and width were found to be more sensitive to source location variations, and hence local geometric features, than vertical plume parameters.
Wall turbulence without walls
- Yoshinori Mizuno, Javier Jiménez
-
- Published online by Cambridge University Press:
- 16 April 2013, pp. 429-455
-
- Article
- Export citation
-
We perform direct numerical simulations of turbulent channels whose inner layer is replaced by an off-wall boundary condition synthesized from a rescaled interior flow plane. The boundary condition is applied within the logarithmic layer, and mimics the linear dependence of the length scales of the velocity fluctuations with respect to the distance to the wall. The logarithmic profile of the mean streamwise velocity is recovered, but only if the virtual wall is shifted to a position different from the location assumed by the boundary condition. In those shifted coordinates, most flow properties are within 5–10 % of full simulations, including the Kármán constant, the fluctuation intensities, the energy budgets and the velocity spectra and correlations. On the other hand, buffer-layer structures do not form, including the near-wall energy maximum, and the velocity fluctuation profiles are logarithmic, strongly suggesting that the logarithmic layer is essentially independent of the near-wall dynamics. The same agreement holds when the technique is applied to large-eddy simulations. The different errors are analysed, especially the reasons for the shifted origin, and remedies are proposed. It is also shown that the length rescaling is required for a stationary logarithmic-like layer. Otherwise, the flow evolves into a state resembling uniformly sheared turbulence.
Aerodynamic noise from a poroelastic edge with implications for the silent flight of owls
- Justin W. Jaworski, N. Peake
-
- Published online by Cambridge University Press:
- 16 April 2013, pp. 456-479
-
- Article
- Export citation
-
The interaction of a turbulent eddy with a semi-infinite poroelastic edge is examined with respect to the effects of both elasticity and porosity on the efficiency of aerodynamic noise generation. The scattering problem is solved using the Wiener–Hopf technique to identify the scaling dependence of the resulting aerodynamic noise on plate and flow properties, including the dependence on a characteristic flow velocity $U$. Special attention is paid to the limiting cases of porous-rigid and impermeable–elastic plate conditions. Asymptotic analysis of these special cases reveals parametric limits where the far-field acoustic power scales like ${U}^{6} $ for a porous edge, and a new finite range of ${U}^{7} $ behaviour is found for an elastic edge, to be compared with the well-known ${U}^{5} $ dependence for a rigid impermeable edge. Further numerical results attempt to address how trailing-edge noise may be mitigated by porosity and flexibility and seek to deepen the understanding of how owls hunt in acoustic stealth.
Receptivity of the boundary layer to vibrations of the wing surface
- A. I. Ruban, T. Bernots, D. Pryce
-
- Published online by Cambridge University Press:
- 16 April 2013, pp. 480-528
-
- Article
- Export citation
-
In this paper we study the generation of Tollmien–Schlichting waves in the boundary layer due to elastic vibrations of the wing surface. The subsonic flow regime is considered with the Mach number outside the boundary layer $M= O(1)$. The flow is investigated based on the asymptotic analysis of the Navier–Stokes equations at large values of the Reynolds number, $\mathit{Re}= {\rho }_{\infty } {V}_{\infty } L/ {\mu }_{\infty } $. Here $L$ denotes the wing section chord; and ${V}_{\infty } $, ${\rho }_{\infty } $ and ${\mu }_{\infty } $ are the free stream velocity, air density and dynamic viscosity, respectively. We assume that in the spectrum of the wing vibrations there is a harmonic that comes in to resonance with the Tollmien–Schlichting wave on the lower branch of the stability curve; this happens when the frequency of the harmonic is a quantity of the order of $({V}_{\infty } / L){\mathit{Re}}^{1/ 4} $. The wavelength, $\ell $, of the elastic vibrations of the wing is assumed to be $\ell \sim L{\mathit{Re}}^{- 1/ 8} $, which has been found to represent a ‘distinguished limit’ in the theory. Still, the results of the analysis are applicable for $\ell \gg L{\mathit{Re}}^{- 1/ 8} $ and $\ell \ll L{\mathit{Re}}^{- 1/ 8} $; the former includes an important case when $\ell = O(L)$. We found that the vibrations of the wing surface produce pressure perturbations in the flow outside the boundary layer, which can be calculated with the help of the ‘piston theory’, which remains valid provided that the Mach number, $M$, is large as compared to ${\mathit{Re}}^{- 1/ 4} $. As the pressure perturbations penetrate into the boundary layer, a Stokes layer forms on the wing surface; its thickness is estimated as a quantity of the order of ${\mathit{Re}}^{- 5/ 8} $. When $\ell = O({\mathit{Re}}^{- 1/ 8} )$ or $\ell \gg {\mathit{Re}}^{- 1/ 8} $, the solution in the Stokes layer appears to be influenced significantly by the compressibility of the flow. The Stokes layer on its own is incapable of producing the Tollmien–Schlichting waves. The reason is that the characteristic wavelength of the perturbation field in the Stokes layer is much larger than that of the Tollmien–Schlichting wave. However, the situation changes when the Stokes layer encounters a wall roughness, which are plentiful in real aerodynamic flows. If the longitudinal extent of the roughness is a quantity of the order of ${\mathit{Re}}^{- 3/ 8} $, then efficient generation of the Tollmien–Schlichting waves becomes possible. In this paper we restrict our attention to the case when the Stokes layer interacts with an isolated roughness. The flow near the roughness is described by the triple-deck theory. The solution of the triple-deck problem can be found in an analytic form. Our main concern is with the flow behaviour downstream of the roughness and, in particular, with the amplitude of the generated Tollmien–Schlichting waves.
On the problem of large-scale magnetic field generation in rotating compressible convection
- B. Favier, P. J. Bushby
-
- Published online by Cambridge University Press:
- 16 April 2013, pp. 529-555
-
- Article
- Export citation
-
Mean-field dynamo theory suggests that turbulent convection in a rotating layer of electrically conducting fluid produces a significant $\alpha $-effect, which is one of the key ingredients in any mean-field dynamo model. Provided that this $\alpha $-effect operates more efficiently than (turbulent) magnetic diffusion, such a system should be capable of sustaining a large-scale dynamo. However, in the Boussinesq model that was considered by Cattaneo & Hughes (J. Fluid Mech., vol. 553, 2006, pp. 401–418) the dynamo produced small-scale, intermittent magnetic fields with no significant large-scale component. In this paper, we consider the compressible analogue of the rotating convective layer that was considered by Cattaneo & Hughes (2006). Varying the horizontal scale of the computational domain, we investigate the dependence of the dynamo upon the rotation rate. Our simulations indicate that these turbulent compressible flows can drive a small-scale dynamo but, even when the layer is rotating very rapidly (with a mid-layer Taylor number of $Ta= 1{0}^{8} $), we find no evidence for the generation of a significant large-scale component of the magnetic field on a dynamical time scale. Like Cattaneo & Hughes (2006), we measure a negligible (time-averaged) $\alpha $-effect when a uniform horizontal magnetic field is imposed across the computational domain. Although the total horizontal magnetic flux is a conserved quantity in these simulations, the (depth-dependent) horizontally averaged magnetic field always exhibits strong fluctuations. If these fluctuations are artificially suppressed within the code, we measure a significant mean electromotive force that is comparable to that found in related calculations in which the $\alpha $-effect is measured using the test-field method, even though we observe no large-scale dynamo action.