The rhythm of fountains: The length and time scales of rise height fluctuations at low and high Froude numbers

The magnitude and frequency of vertical fluctuations of the top of an axisymmetric miscible Boussinesq fountain forms the focus of this work. We present measurements of these quantities for saline-aqueous fountains in uniform quiescent surroundings. Our results span source Froude numbers 0.3 ≤ Fr 0 ≤ 40 and, thereby, encompass very weak, weak, intermediate and forced classes of fountain. We identify distinct scalings, based on known quantities at the fountain source, for the frequency of fountain height fluctuations which collapse our data within bands of Fr0. Notably, our scalings reveal that the (dimensionless) frequency takes a constant value within each band. These results highlight characteristic time scales for the fluctuations which we decompose into a single, physically apparent, length scale and velocity scale within each band. Moreover, within one particular band, spanning source Froude numbers towards the lower end of the full range considered, we identify unexpectedly long-period fluctuations indicating a near balance of inertia and (opposing) buoyancy at the source. Our analysis identifies four distinct classes of fluctuation behaviour (four bands of Fr 0) and this classification matches well with existing classifications of fountains based on rise heights. As such, we show that an analysis of the behaviour of the fountain top alone, rather than the entire fountain, provides an alternative approach to classifying fountains. The similarity of classifications based on the two different methods confirms that the boundaries between classes mark tangible changes in the physics of fountains. For high Fr0 we show that the dominant fluctuations occur at the scale of the largest eddies which can be contained within the fountain near its top. Extending this, we develop a Strouhal number, Strtop, based on experimental measures of the fountain top, defined such that Strtop = 1 would suggest the dominant fluctuations are caused by a continual cycle of eddies forming and collapsing at this largest physical scale. For high- Fr 0 fountains we find Strtop ≈ 0. 9. © 2013 Cambridge University Press.

Fluid-structure interaction simulations of a tension-cone inflatable aerodynamic decelerator

A series of fluid-structure interaction simulations of an aerodynamic tension-cone supersonic decelerator prototype intended for large mass payload deployment in planetary explorations are discussed. The fluid-structure interaction computations combine large deformation analysis of thin shells with large-eddy simulation of compressible turbulent flows using a loosely coupled approach to enable quantification of the dynamics of the vehicle. The simulation results are compared with experiments carried out at the NASA Glenn Research Center. Reasonably good agreement between the simulations and the experiment is observed throughout a deflation cycle. The simulations help to illuminate the details of the dynamic progressive buckling of the tension-cone decelerator that ultimately results in the collapse of the structure as the inflation pressure is decreased. Furthermore, the tension-cone decelerator exhibits a transient oscillatory behavior under impulsive loading that ultimately dies out. The frequency of these oscillations was determined to be related to the acoustic time scale in the compressed subsonic region between the bow shock and the structure. As shown, when the natural frequency of the structure and the frequency of the compressed subsonic region approximately match, the decelerator exhibits relatively large nonaxisymetric oscillations. The observed response appears to be a fluid-structure interaction resonance resulting from an acoustic chamber (pistonlike) mode exciting the structure. Copyright © 2013 by Christopher Porter, R. Mark Rennie, Eric J. Jumper.

Aerodynamic noise from a poroelastic trailing edge with implications for the silent flight of owls

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 scattered aerodynamic noise. The scattering problem is solved using the Wiener-Hopf technique for constant plate properties to identify their scaling dependence on the resulting aerodynamic noise, including the dependence on flight velocity, where special attention is paid to the limiting cases of rigid, porous and elastic, impermeable plate conditions. Results from these analyses attempt to address how trailing edge noise may be mitigated by porosity and seek to deepen the understanding of how owls hunt in acoustic stealth. © 2012 by Justin W. Jaworski and Nigel Peake. Published by the American Institute of Aeronautics and Astronautics, Inc.

Kelvin-Helmholtz-like instabilities in turbulent flows over complex surfaces

Riblets are small surface protrusions aligned with the flow direction, which confer an anisotropic roughness to the surface [6]. We have recently reported that the transitional-roughness effect in riblets, which limits their performance, is due to a Kelvin–Helmholtz-like instability of the overlying mean flow [7]. According to our DNSs, the instability sets on as the Reynolds number based on the roughness size of the riblets increases, and coherent, elongated spanwise vortices begin to develop immediately above the riblet tips, causing the degradation of the drag-reduction effect. This is a very novel concept, since prior studies had proposed that the degradation was due to the interaction of riblets with the flow as independent units, either to the lodging of quasi-streamwise vortices in the surface grooves [2] or to the shedding of secondary streamwise vorticity at the riblet peaks [9]. We have proposed an approximate inviscid analysis for the instability, in which the presence of riblets is modelled through an average boundary condition for an overlying, spanwise-independent mean flow. This simplification lacks the accuracy of an exact analysis [4], but in turn applies to riblet surfaces in general. Our analysis succeeds in predicting the riblet size for the onset of the instability, while qualitatively reproducing the wavelengths and shapes of the spanwise structures observed in the DNSs. The analysis also connects the observations with the Kelvin–Helmholtz instability of mixing layers. The fundamental riblet length scale for the onset of the instability is a ‘penetration length,’ which reflects how easily the perturbation flow moves through the riblet grooves. This result is in excellent agreement with the available experimental evidence, and has enabled the identification of the key geometric parameters to delay the breakdown. Although the appearance of elongated spanwise vortices was unexpected in the case of riblets, similar phenomena had already been observed over other rough [3], porous [1] and permeable [11] surfaces, as well as over plant [5,14] and urban [12] canopies, both in the transitional and in the fully-rough regimes. However, the theoretical analyses that support the connection of these observations with the Kelvin–Helmholtz instability are somewhat scarce [7, 11, 13]. It has been recently proposed that Kelvin–Helmholtz-like instabilities are a dominant feature common to “obstructed” shear flows [8]. It is interesting that the instability does not require an inflection point to develop, as is often claimed in the literature. The Kelvin-Helmholtz rollers are rather triggered by the apparent wall-normal-transpiration ability of the flow at the plane immediately above the obstructing elements [7,11]. Although both conditions are generally complementary, if wall-normal transpiration is not present the spanwise vortices may not develop, even if an inflection point exists within the roughness [10]. REFERENCES [1] Breugem, W. P., Boersma, B. J. & Uittenbogaard, R. E. 2006 J. Fluid Mech. 562, 35–72. [2] Choi, H., Moin, P. & Kim, J. 1993 J. Fluid Mech. 255, 503–539. [3] Coceal, O., Dobre, A., Thomas, T. G. & Belcher, S. E. 2007 J. Fluid Mech. 589, 375–409. [4] Ehrenstein, U. 2009 Phys. Fluids 8, 3194–3196. [5] Finnigan, J. 2000 Ann. Rev. Fluid Mech. 32, 519–571. [6] Garcia-Mayoral, R. & Jimenez, J. 2011 Phil. Trans. R. Soc. A 369, 1412–1427. [7] Garcia-Mayoral, R. & Jimenez, J. 2011 J. Fluid Mech. doi: 10.1017/jfm.2011.114. [8] Ghisalberti, M. 2009 J. Fluid Mech. 641, 51–61. [9] Goldstein, D. B. & Tuan, T. C. 1998 J. Fluid Mech. 363, 115–151. [10] Hahn, S., Je, J. & Choi, H. 2002 J. Fluid Mech. 450, 259–285. [11] Jimenez, J., Uhlman, M., Pinelli, A. & G., K. 2001 J. Fluid Mech. 442, 89–117. [12] Letzel, M. O., Krane, M. & Raasch, S. 2008 Atmos. Environ. 42, 8770–8784. [13] Py, C., de Langre, E. & Moulia, B. 2006 J. Fluid Mech. 568, 425–449. [14] Raupach, M. R., Finnigan, J. & Brunet, Y. 1996 Boundary-Layer Meteorol. 78, 351–382.

The interaction of riblets with wall-bounded turbulence

The purpose of this thesis is to give answer to the question: why do riblets stop working for a certain size? Riblets are small surface grooves aligned in the mean direction of an overlying turbulent flow, designed specifically to reduce the friction between the flow and the surface. They were inspired by biological surfaces, like the oriented denticles in the skin of fastswimming sharks, and were the focus of a significant amount of research in the late eighties and nineties. Although it was found that the drag reduction depends on the riblet size scaled in wall units, the physical mechanisms implicated have not been completely understood up to now. It has been explained how riblets of vanishing size interact with the turbulent flow, producing a change in the drag proportional to their size, but that is not the regime of practical interest. The optimum performance is achieved for larger sizes, once that linear behavior has broken down, but before riblets begin adopting the character of regular roughness and increasing drag. This regime, which is the most relevant from a technological perspective, was precisely the less understood, so we have focused on it. Our efforts have followed three basic directions. First, we have re-assessed the available experimental data, seeking to identify common characteristics in the optimum regime across the different existing riblet geometries. This study has led to the proposal of a new length scale, the square root of the groove crosssection, to substitute the traditional peak-to-peak spacing. Scaling the riblet dimension with this length, the size of breakdown of the linear behavior becomes roughly universal. This suggests that the onset of the breakdown is related to a certain, fixed value of the cross-section of the groove. Second, we have conducted a set of direct numerical simulations of the turbulent flow over riblets, for sizes spanning the full drag reduction range. We have thus been able to reproduce the gradual transition between the different regimes. The spectral analysis of the flows has proven particularly fruitful, since it has made possible to identify spanwise rollers immediately above the riblets, which begin to appear when the riblet size is close to the optimum. This is a quite surprising feature of the flow, not because of the uniqueness of the phenomenon, which had been reported before for other types of complex and porous surfaces, but because most previous studies had focused on the detail of the flow above each riblet as a unit. Our novel approach has provided the adequate tools to capture coherent structures with an extended spanwise support, which interact with the riblets not individually, but collectively. We have also proven that those spanwise structures are responsible for the increase in drag past the viscous breakdown. Finally, we have analyzed the stability of the flow with a simplified model that connects the appearance of rollers to a Kelvin–Helmholtz-like instability, as is the case also for the flow over plant canopies and porous surfaces. In spite of the model emulating the presence of riblets only in an averaged, general fashion, it succeeds to capture the essential attributes of the breakdown, and provides a theoretical justification for the scaling with the groove cross-section.

An investigation of higher-order multi-objective optimisation for 3D aerodynamic shape design

We investigate the performance of different variants of a suitably tailored Tabu Search optimisation algorithm on a higher-order design problem. We consider four objective func- tions to describe the performance of a compressor stator row, subject to a number of equality and inequality constraints. The same design problem has been previously in- vestigated through single-, bi- and three-objective optimisation studies. However, in this study we explore the capabilities of enhanced variants of our Multi-objective Tabu Search (MOTS) optimisation algorithm in the context of detailed 3D aerodynamic shape design. It is shown that with these enhancements to the local search of the MOTS algorithm we can achieve a rapid exploration of complicated design spaces, but there is a trade-off be- tween speed and the quality of the trade-off surface found. Rapidly explored design spaces reveal the extremes of the objective functions, but the compromise optimum areas are not very well explored. However, there are ways to adapt the behaviour of the optimiser and maintain both a very efficient rate of progress towards the global optimum Pareto front and a healthy number of design configurations lying on the trade-off surface and exploring the compromise optimum regions. These compromise solutions almost always represent the best qualitative balance between the objectives under consideration. Such enhancements to the effectiveness of design space exploration make engineering design optimisation with multiple objectives and robustness criteria ever more practicable and attractive for modern advanced engineering design. Finally, new research questions are addressed that highlight the trade-offs between intelligence in optimisation algorithms and acquisition of qualita- tive information through computational engineering design processes that reveal patterns and relations between design parameters and objective functions, but also speed versus optimum quality. © 2012 AIAA.

Application of the multi-objective Alliance algorithm to a benchmark aerodynamic optimization problem

This paper introduces a new version of the multiobjective Alliance Algorithm (MOAA) applied to the optimization of the NACA 0012 airfoil section, for minimization of drag and maximization of lift coefficients, based on eight section shape parameters. Two software packages are used: XFoil which evaluates each new candidate airfoil section in terms of its aerodynamic efficiency, and a Free-Form Deformation tool to manage the section geometry modifications. Two versions of the problem are formulated with different design variable bounds. The performance of this approach is compared, using two indicators and a statistical test, with that obtained using NSGA-II and multi-objective Tabu Search (MOTS) to guide the optimization. The results show that the MOAA outperforms MOTS and obtains comparable results with NSGA-II on the first problem, while in the other case NSGA-II is not able to find feasible solutions and the MOAA is able to outperform MOTS. © 2013 IEEE.

Experimental study of surface roughness noise

A turbulent boundary-layer flow over a rough wall generates a dipole sound field as the near-field hydrodynamic disturbances in the turbulent boundary-layer scatter into radiated sound at small surface irregularities. In this paper, phased microphone arrays are applied to the experimental study of surface roughness noise. The radiated sound from two rough plates and one smooth plate in an open jet is measured at three streamwise locations, and the beamforming source maps demonstrate the dipole directivity. Higher source strengths can be observed in the rough plates than the smooth plate, and the rough plates also enhance the trailing-edge noise. A prediction scheme in previous theoretical work is used to describe the strength of a distribution of incoherent dipoles over the rigid plate and to simulate the sound detected by the microphone array. Source maps of measurement and simulation exhibit encouraging similarities in both source pattern and source strength, which confirms the dipole nature and the predicted magnitude of roughness noise. The simulations underestimate the streamwise gradient of the source strengths and overestimate the source strengths at the highest frequency. © 2007 by Yu Liu and Ann P. Dowling.

Aerodynamic design of high end wall angle turbine stages-part I: Methodology development

A design methodology is presented for turbines in an annulus with high end wall angles. Such stages occur where large radial offsets between the stage inlet and stage outlet are required, for example in the first stage of modern low pressure turbines, and are becoming more prevalent as bypass ratios increase. The turbine vanes operate within s-shaped ducts which result in meridional curvature being of a similar magnitude to the bladeto-blade curvature. Through a systematic series of idealized computational cases, the importance of two aspects of vane design are shown. First, the region of peak end wall meridional curvature is best located within the vane row. Second, the vane should be leant so as to minimize spanwise variations in surface pressure-this condition is termed "ideal lean." This design philosophy is applied to the first stage of a low pressure turbine with high end wall angles. © 2014 by ASME.

Aerodynamic design of high end wall angle turbine stages- part II: Experimental verification

An experimental investigation of a turbine stage featuring very high end wall angles is presented. The initial turbine design did not achieve a satisfactory performance and the difference between the design predictions and the test results was traced to a large separated region on the rear suction-surface. To improve the agreement between computational fluid dynamics (CFD) and experiment, it was found necessary to modify the turbulence modeling employed. The modified CFD code was then used to redesign the vane, and the changes made are described. When tested, the performance of the redesigned vane was found to have much closer agreement with the predictions than the initial vane. Finally, the flowfield and performance of the redesigned stage are compared to a similar turbine, designed to perform the same duty, which lies in an annulus of moderate end wall angles. A reduction in stage efficiency of at least 2.4% was estimated for the very high end wall angle design. © 2014 by ASME.

Influence of film cooling hole angles and geometries on aerodynamic loss and net heat flux reduction

Turbine design engineers have to ensure that film cooling can provide sufficient protection to turbine blades from the hot mainstream gas, while keeping the losses low. Film cooling hole design parameters include inclination angle (a), compound angle (b), hole inlet geometry, and hole exit geometry. The influence of these parameters on aerodynamic loss and net heat flux reduction is investigated, with loss being the primary focus. Low-speed flat plate experiments have been conducted at momentum flux ratios of IR=0.16, 0.64, and 1.44. The film cooling aerodynamic mixing loss, generated by the mixing of mainstream and coolant, can be quantified using a three-dimensional analytical model that has been previously reported by the authors. The model suggests that for the same flow conditions, the aerodynamic mixing loss is the same for holes with different a and b but with the same angle between the mainstream and coolant flow directions (angle k). This relationship is assessed through experiments by testing two sets of cylindrical holes with different a and b: one set with k=35 deg, and another set with k=60 deg. The data confirm the stated relationship between α, β, k and the aerodynamic mixing loss. The results show that the designer should minimize k to obtain the lowest loss, but maximize b to achieve the best heat transfer performance. A suggestion on improving the loss model is also given. Five different hole geometries (α=35.0 deg, β=0 deg) were also tested: cylindrical hole, trenched hole, fan-shaped hole, D-Fan, and SD-Fan. The D-Fan and the SD-Fan have similar hole exits to the fan-shaped hole but their hole inlets are laterally expanded. The external mixing loss and the loss generated inside the hole are compared. It was found that the D-Fan and the SD-Fan have the lowest loss. This is attributed to their laterally expanded hole inlets, which lead to significant reduction in the loss generated inside the holes. As a result, the loss of these geometries is≈50% of the loss of the fan-shaped hole at IR=0.64 and 1.44. © 2013 by ASME.