Fluid-structure interaction with mean flow: Over-scattering and unstable resonance growth

It is known theoretically [1-3] that infinitely long fluid loaded plates in mean flow exhibit a range of unusual phenomena in the 'long time' limit. These include convective instability, absolute instability and negative energy waves which are destabilized by dissipation. However, structures are necessarily of finite length and may have discontinuities. Moreover, linear instability waves can only grow over a limited number of cycles before non-linear effects become dominant. We have undertaken an analytical and computational study to investigate the response of finite, discontinuous plates to ascertain if these unusual effects might be realized in practice. Analytically, we take a "wave scattering" [2,4] - as opposed to a "modal superposition" [5] - view of the fluttering plate problem. First, we solve for the scattering coefficients of localized plate discontinuities and identify a range of parameter space, well outside the convective instability regime, where over-scattering or amplified reflection/transmission occurs. These are scattering processes that draw energy from the mean flow into the plate. Next, we use the Wiener-Hopf technique to solve for the scattering coefficients from the leading and trailing edges of a baffled plate. Finally, we construct the response of a finite, baffled plate by a superposition of infinite plate propagating waves continuously scattering off the plate ends and solve for the unstable resonance frequencies and temporal growth rates for long plates. We present a comparison between our computational results and the infinite plate theory. In particular, the resonance response of a moderately sized plate is shown to be in excellent agreement with our long plate analytical predictions. Copyright © 2010 by ASME.

Eigenmodes of lined flow ducts with rigid splices

This paper presents an analytic expression for the acoustic eigenmodes of a cylindrical lined duct with rigid axially running splices in the presence of flow. The cylindrical duct is considered to be uniformly lined except for two symmetrically positioned axially running rigid liner splices. An exact analytic expression for the acoustic pressure eigenmodes is given in terms of an azimuthal Fourier sum, with the Fourier coefficients given by a recurrence relation. Since this expression is derived using a Greens function method, the completeness of the expansion is guaranteed. A numerical procedure is described for solving this recurrence relation, which is found to converge exponentially with respect to number of Fourier terms used and is in practice quick to compute; this is then used to give several numerical examples for both uniform and sheared mean flow. An asymptotic expression is derived to directly calculate the pressure eigenmodes for thin splices. This asymptotic expression is shown to be quantitatively accurate for ducts with very thin splices of less than 1 % unlined area and qualitatively helpful for thicker splices of the order of 6 % unlined area. A thin splice is in some cases shown to increase the damping of certain acoustic modes. The influences of thin splices and thin boundary layers are compared and found to be of comparable magnitude for the parameters considered. Trapped modes at the splices are also identified and investigated. © 2011 Cambridge University Press.

Slip flow and heat transfer in microbearings with fractal surface topographies

The effects of random surface roughness on slip flow and heat transfer in microbearings are investigated. A three-dimensional random surface roughness model characterized by fractal geometry is used to describe the multiscale self-affine roughness, which is represented by the modified two-variable Weierstrass- Mandelbrot (W-M) functions, at micro-scale. Based on this fractal characterization, the roles of rarefaction and roughness on the thermal and flow properties in microbearings are predicted and evaluated using numerical analyses and simulations. The results show that the boundary conditions of velocity slip and temperature jump depend not only on the Knudsen number but also on the surface roughness. It is found that the effects of the gas rarefaction and surface roughness on flow behavior and heat transfer in the microbearing are strongly coupled. The negative influence of roughness on heat transfer found to be the Nusselt number reduction. In addition, the effects of temperature difference and relative roughness on the heat transfer in the bearing are also analyzed and discussed. © 2012 Elsevier Ltd. All rights reserved.

CFD modeling of single-phase flow in a packed bed with MRI validation

Computational fluid dynamics (CFD) simulations are becoming increasingly widespread with the advent of more powerful computers and more sophisticated software. The aim of these developments is to facilitate more accurate reactor design and optimization methods compared to traditional lumped-parameter models. However, in order for CFD to be a trusted method, it must be validated using experimental data acquired at sufficiently high spatial resolution. This article validates an in-house CFD code by comparison with flow-field data obtained using magnetic resonance imaging (MRI) for a packed bed with a particle-to-column diameter ratio of 2. Flows characterized by inlet Reynolds numbers, based on particle diameter, of 27, 55, 111, and 216 are considered. The code used employs preconditioning to directly solve for pressure in low-velocity flow regimes. Excellent agreement was found between the MRI and CFD data with relative error between the experimentally determined and numerically predicted flow-fields being in the range of 3-9%. © 2012 American Institute of Chemical Engineers (AIChE).

Numerical computation of laminar flow in a heated rotating cavity with an axial throughflow of air

A heated rotating cavity with an axial throughflow of cooling air is used as a model for the flow in the cylindrical cavities between adjacent discs of a high-pressure gas-turbine compressor. In an engine the flow is expected to be turbulent, the limitations of this laminar study are fully realised but it is considered an essential step to understand the fundamental nature of the flow. The three-dimensional, time-dependent governing equations are solved using a code based on the finite volume technique and a multigrid algorithm. The computed flow structure shows that flow enters the cavity in one or more radial arms and then forms regions of cyclonic and anticyclonic circulation. This basic flow structure is consistent with existing experimental evidence obtained from flow visualization. The flow structure also undergoes cyclic changes with time. For example, a single radial arm, and pair of recirculation regions can commute to two radial arms and two pairs of recirculation regions and then revert back to one. The flow structure inside the cavity is found to be heavily influenced by the radial distribution of surface temperature imposed on the discs. As the radial location of the maximum disc temperature moves radially outward, this appears to increase the number of radial arms and pairs of recirculation regions (from one to three for the distributions considered here). If the peripheral shroud is also heated there appear to be many radial arms which exchange fluid with a strong cyclonic flow adjacent to the shroud. One surface temperature distribution is studied in detail and profiles of the relative tangential and radial velocities are presented. The disc heat transfer is also found to be influenced by the disc surface temperature distribution. It is also found that the computed Nusselt numbers are in reasonable accord over most of the disc surface with a correlation found from previous experimental measurements. © 1994, MCB UP Limited.

A novel MPT noise prediction methodology for highly-integrated propulsion systems with inlet flow distortion

Embedded propulsion systems, such as for example used in advanced hybrid-wing body aircraft, can potentially offer major fuel burn and noise reduction benefits but introduce challenges in the aerodynamic and acoustic integration of the high-bypass ratio fan system. A novel approach is proposed to quantify the effects of non-uniform flow on the generation and propagation of multiple pure tone noise (MPTs). The new method is validated on a conventional inlet geometry first. The ultimate goal is to conduct a parametric study of S-duct inlets in order to quantify the effects of inlet design parameters on the acoustic signature. The key challenge is that the mechanism underlying the distortion transfer, noise source generation and propagation through the non-uniform flow field are inherently coupled such that a simultaneous computation of the aerodynamics and acoustics is required. The technical approach is based on a body force description of the fan blade row that is able to capture the distortion transfer and the MPT noise generation mechanisms while greatly reducing computational cost. A single, 3-D full-wheel unsteady CFD simulation, in which the Euler equations are solved to second-order spatial and temporal accuracy, simultaneously computes the MPT noise generation and its propagation in distorted mean flow. Several numerical tools were developed to enable the implementation of this new approach. Parametric studies were conducted to determine appropriate grid and time step sizes for the propagation of acoustic waves. The Ffowcs-Williams and Hawkings integral method is used to propagate the noise to far field receivers. Non-reflecting boundary conditions are implemented through the use of acoustic buffer zones. The body force modeling approach is validated and proof-of-concept studies demonstrate the generation of disturbances at both blade-passing and shaft-order frequencies using the perturbed body force method. The full methodology is currently being validated using NASA's Source Diagnostic Test (SDT) fan and inlet geometry. Copyright © 2009 by Jeff Defoe, Alex Narkaj & Zoltan Spakovszky.

Scaling of transitionally-rough flow structures in DNSes of channels with riblets up to Reτ=550

The linear, drag-reducing effect of vanishingly small riblets breaks down once their size is in the transitionally-rough regime. We have previously reported that this breakdown is caused by the additional Reynolds stresses produced by the appearance of elongated spanwise rollers just above the riblet surface. These rollers are related with the Kelvin--Helmholtz instability of free shear layers, and to similar structures appearing over other rough and porous surfaces. However, because of the limited Reτ=180 in our previous DNSes, it could not be determined whether those structures scaled in inner or outer units. Furthermore, it is questionable if results in the transitionally-rough regime at Reτ=180 can be extrapolated to configurations of practical interest. At such small Reynolds numbers, roughness of transitional size can perturb a large portion of the boundary layer, which is not the case in most industrial and atmospheric applications. To clarify these issues we have conducted a set of DNSes at Reτ=550. Our results indicate that the spanwise rollers scale in wall units, and support the validity of the extrapolation to configurations of practical interest.

The acoustic analogy in an annular duct with swirling mean flow

This paper is concerned with modelling the effects of swirling flow on turbomachinery noise. We develop an acoustic analogy to predict sound generation in a swirling and sheared base flow in an annular duct, including the presence of moving solid surfaces to account for blade rows. In so doing we have extended a number of classical earlier results, including Ffowcs Williams & Hawkings' equation in a medium at rest with moving surfaces, and Lilley's equation for a sheared but non-swirling jet. By rearranging the Navier-Stokes equations we find a single equation, in the form of a sixth-order differential operator acting on the fluctuating pressure field on the left-hand side and a series of volume and surface source terms on the right-hand side; the form of these source terms depends strongly on the presence of swirl and radial shear. The integral form of this equation is then derived, using the Green's function tailored to the base flow in the (rigid) duct. As is often the case in duct acoustics, it is then convenient to move into temporal, axial and azimuthal Fourier space, where the Green's function is computed numerically. This formulation can then be applied to a number of turbomachinery noise sources. For definiteness here we consider the noise produced downstream when a steady distortion flow is incident on the fan from upstream, and compare our results with those obtained using a simplistic but commonly used Doppler correction method. We show that in all but the simplest case the full inclusion of swirl within an acoustic analogy, as described in this paper, is required. © 2013 Cambridge University Press.

Aeroacoustic source measurement methods for characterizing the sound generated by ducted flow devices with higher-order modes.

The International Organization for Standardization (ISO) method 5136 is widely used in industry and academia to determine the sound power radiated into a duct by fans and other flow devices. The method involves placing the device at the center of a long cylindrical duct with anechoic terminations at each end to eliminate reflections. A single off-axis microphone is used on the inlet and outlet sides that can theoretically capture the plane-wave mode amplitudes but this does not provide enough information to fully account for higher-order modes. In this study, the "two-port" source model is formulated to include higher-order modes and applied for the first three modes. This requires six independent surface pressure measurements on each side or "port." The resulting experimental set-up is much shorter than the ISO rig and does not require anechoic terminations. An array of six external loudspeaker sources is used to characterize the passive part of the two-port model and the set-up provides a framework to account for transmission of higher-order modes through a fan. The relative importance of the higher-order modes has been considered and their effect on inaccuracies when using the ISO method to find source sound power has been analyzed.