Adjoint algorithms for the Navier-Stokes equations in the low Mach number limit

This paper describes a derivation of the adjoint low Mach number equations and their implementation and validation within a global mode solver. The advantage of using the low Mach number equations and their adjoints is that they are appropriate for flows with variable density, such as flames, but do not require resolution of acoustic waves. Two versions of the adjoint are implemented and assessed: a discrete-adjoint and a continuous-adjoint. The most unstable global mode calculated with the discrete-adjoint has exactly the same eigenvalue as the corresponding direct global mode but contains numerical artifacts near the inlet. The most unstable global mode calculated with the continuous-adjoint has no numerical artifacts but a slightly different eigenvalue. The eigenvalues converge, however, as the timestep reduces. Apart from the numerical artifacts, the mode shapes are very similar, which supports the expectation that they are otherwise equivalent. The continuous-adjoint requires less resolution and usually converges more quickly than the discrete-adjoint but is more challenging to implement. Finally, the direct and adjoint global modes are combined in order to calculate the wavemaker region of a low density jet. © 2011 Elsevier Inc.

Assessment of the contribution of surface roughness to airframe noise

The generation of sound by turbulent boundary-layer flow at low Mach number over a rough wall is investigated by applying a theoretical model that describes the scattering of the turbulence near field into sound by roughness elements. Attention is focused on the numerical method to approximately quantify the absolute level of far-field radiated roughness noise. Models for the source statistics are obtained by scaling smooth-wall data by the increased skin friction velocity and boundary-layer thickness for a rough surface. Numerical integration is performed to determine the roughness noise, and it reproduces the spectral characteristics of the available empirical formula and experimental data. Experiments are conducted to measure the radiated sound from two rough plates in an open jet The measured noise spectra of the rough plates are above that of a smooth plate in 1-2.5 kHz frequency and exhibit reasonable agreement with the predicted level. Estimates of the roughness noise for a Boeing 757 sized aircraft wing with idealized levels of surface roughness show that hi the high-frequency region the sound radiated from surface roughness may exceed that from the trailing edge, and higher overall sound pressure levels are observed for the roughness noise. The trailing edge noise is also enhanced by surface roughness somewhat A parametric study indicates that roughness height and roughness density significantly affect the roughness noise with roughness height having the dominant effect The roughness noise directivity varies with different levels of surface roughness. Copyright © 2007 by the American Institute of Aeronautics and Astronautics, Inc. All rights reserved.

Effects of surface roughness on airframe noise

The generation of sound by turbulent boundary layer flow at low Mach number over a rough wall is investigated by applying the theoretical model which describes the scattering of the turbulence near field into sound by roughness elements. Attention is focused on the numerical method to approximately quantify the absolute level of the roughness noise radiated to far field. Empirical models for the source statistics are obtained by scaling smooth-wall data through increased skin friction velocity and boundary layer thickness for the rough surface. Numerical integration is performed to determine the roughness noise, and it reproduces the spectral characteristics of the available empirical formula and experimental data. Experiments are conducted to measure the radiated sound from two rough plates in an open jet by four 1/2'' free-field condenser microphones. The measured noise spectra of the rough plates are above that of a smooth plate in 1-2.5 kHz frequency and exhibits encouraging agreement with the predicted spectra. Also, a phased microphone array is utilized to localize the sound source, and it confirms that the rough plates generate higher source strengthes in this frequency range. A parametric study illustrates that the roughness height and roughness density significantly affect the far-field radiated roughness noise with the roughness height having the dominant effect. The estimates of the roughness noise for a Boeing 757 sized aircraft wing show that in high frequency region the sound radiated from surface roughness may exceed that from the trailing edge, and higher overall sound pressure levels for the roughness noise are also observed.

Modal scattering at an impedance transition in a lined flow duct

An explicit Wiener-Hopf solution is derived to describe the scattering of duct modes at a hard-soft wall impedance transition in a circular duct with uniform mean flow. Specifically, we have a circular duct r = 1, - ∞ < x < ∞ with mean flow Mach number M > 0 and a hard wall along x < 0 and a wall of impedance Z along x > 0. A minimum edge condition at x = 0 requires a continuous wall streamline r = 1 + h(x, t), no more singular than h = Ο(x1/2) for x ↓ 0. A mode, incident from x < 0, scatters at x = 0 into a series of reflected modes and a series of transmitted modes. Of particular interest is the role of a possible instability along the lined wall in combination with the edge singularity. If one of the "upstream" running modes is to be interpreted as a downstream-running instability, we have an extra degree of freedom in the Wiener-Hopf analysis that can be resolved by application of some form of Kutta condition at x = 0, for example a more stringent edge condition where h = Ο(x3/2) at the downstream side. The question of the instability requires an investigation of the modes in the complex frequency plane and therefore depends on the chosen impedance model, since Z = Z (ω) is essentially frequency dependent. The usual causality condition by Briggs and Bers appears to be not applicable here because it requires a temporal growth rate bounded for all real axial wave numbers. The alternative Crighton-Leppington criterion, however, is applicable and confirms that the suspected mode is usually unstable. In general, the effect of this Kutta condition is significant, but it is particularly large for the plane wave at low frequencies and should therefore be easily measurable. For ω → 0, the modulus fends to |R001| → (1 + M)/(1 -M) without and to 1 with Kutta condition, while the end correction tends to ∞ without and to a finite value with Kutta condition. This is exactly the same behaviour as found for reflection at a pipe exit with flow, irrespective if this is uniform or jet flow.

Normal shock boundary layer control with various vortex generator geometries

Various vortex generators which include ramp, split-ramp and a new hybrid concept "ramped-vane" are investigated under normal shock conditions with a diffuser at Mach number of 1.3. The dimensions of the computational domain were designed using Reynolds Average Navier-Stokes studies to be representative of the flow in an external-compression supersonic inlet. Using this flow geometry, various vortex generator concepts were studied with Implicit Large Eddy Simulation. In general, the ramped-vane provided increased vorticity compared to the other devices and reduced the separation length downstream of the device centerline. In addition, the size, edge gap and streamwise position respect to the shock were studied for the ramped-vane and it was found that a height of about half the boundary thickness and a large trailing edge gap yielded a fully attached flow downstream of the device. This ramped-vane also provided the largest reduction in the turbulent kinetic energy and pressure fluctuations. Additional benefits include negligible drag while the reductions in boundary layer displacement thickness and shape factor were seen compared to other devices. © 2010 by Sang Lee.

Normal shock boundary layer control with various vortex generator geometries

Various vortex generators which include ramp, split-ramp and a new hybrid concept "ramped-vane" are investigated under normal shock conditions with a diffuser at Mach number of 1.3. The dimensions of the computational domain were designed using Reynolds Average Navier-Stokes studies to be representative of the flow in an external-compression supersonic inlet. Using this flow geometry, various vortex generator concepts were studied with Implicit Large Eddy Simulation. In general, the ramped-vane provided increased vorticity compared to the other devices and reduced the separation length downstream of the device centerline. In addition, the size, edge gap and streamwise position respect to the shock were studied for the ramped-vane and it was found that a height of about half the boundary thickness and a large trailing edge gap yielded a fully attached flow downstream of the device. This ramped-vane also provided the largest reduction in the turbulent kinetic energy and pressure fluctuations. Additional benefits include negligible drag while the reductions in boundary layer displacement thickness and shape factor were seen compared to other devices. © 2011 Elsevier Ltd.

Simulation of three-dimensional viscous flow in turbomachinery geometries using a solution-adaptive unstructured mesh methodology

This paper presents a numerical method for the simulation of flow in turbomachinery blade rows using a solution-adaptive mesh methodology. The fully three-dimensional, compressible, Reynolds-averaged Navier-Stokes equations with k-ε turbulence modeling (and low Reynolds number damping terms) are solved on an unstructured mesh formed from tetrahedral finite volumes. At stages in the solution, mesh refinement is carried out based on flagging cell faces with either a fractional variation of a chosen variable (like Mach number) greater than a given threshold or with a mean value of the chosen variable within a given range. Several solutions are presented, including that for the highly three-dimensional flow associated with the corner stall and secondary flow in a transonic compressor cascade, to demonstrate the potential of the new method.

The development of turbine exit flow in a swan-necked inter-stage diffuser

This paper describes both the migration and dissipation of flow phenomena downstream of a transonic high-pressure turbine stage. The geometry of the HP stage exit duct considered is a swan-necked diffuser similar to those likely to be used in future engine designs. The paper contains results both from an experimental programme in a turbine test facility and from numerical predictions. Experimental data was acquired using three fast-response aerodynamic probes capable of measuring Mach number, whirl angle, pitch angle, total pressure and static pressure. The probes were used to make time-resolved area traverses at two axial locations downstream of the rotor trailing edge. A 3D time-unsteady viscous Navier-Stokes solver was used for the numerical predictions. The unsteady exit flow from a turbine stage is formed from rotordependent phenomena (such as the rotor wake, the rotor trailing edge recompression shock, the tip-leakage flow and the hub secondary flow) and vane-rotor interaction dependant phenomena. This paper describes the time-resolved behaviour and three-dimensional migration paths of both of these phenomena as they convect downstream. It is shown that the inlet flow to a downstream vane is dominated by two corotating vortices, the first caused by the rotor tip-leakage flow and the second by the rotor hub secondary flow. At the inlet plane of the downstream vane the wake is extremely weak and the radial pressure gradient is shown to have caused the majority of the high loss wake fluid to be located between the mid-height of the passage and the casing wall. The structure of the flow indicates that between a high pressure stage and a downstream vane simple two-dimensional blade row interaction does not occur. The results presented in this paper indicate that the presence of an upstream stage is likely to significantly alter the structure of the secondary flow within a downstream vane. The paper also shows that vane-rotor interaction within the upstream stage causes a 10° circumferential variation in the inlet flow angle of the 2nd stage vane.

The development of aerodynamic performance measurements in a transient test facility

Transient test facilities offer the potential for the simultaneous study of turbine aerodynamic performance, unsteady flow phenomena and the heat transfer characteristics of a turbine stage. This paper describes the development of aerodynamic performance measurement techniques in the Oxford Rotor Facility (ORF). The solutions to the technological issues involved with transient testing presented in this paper are expected to achieve levels of precision uncertainty comparable with traditional steady flow test rigs. The theoretical background to the measurement of aerodynamic performance is presented together with a comprehensive pre-test uncertainty analysis. The instrumentation scheme for the measurement of stage mass flow rate is discussed in detail, the measurements of shaft power, total inlet enthalpy, and stage pressure ratio are also outlined. The current working section features a 62% scale, 1-1/2 stage, high-pressure shroudless transonic turbine. The required inlet flow conditions are provided by an Isentropic Light Piston Tunnel (ILPT) with a quasi-steady state run time of approximately 70ms. The testing is conducted at engine representative specific speed, pressure ratio, gas-to-wall temperature ratio, Mach number and Reynolds number.

The effect of work processes on the casing heat transfer of a transonic turbine

This paper considers the effect of the rotor tip on the casing heat load of a transonic axial flow turbine. The aim of the research is to understand the dominant causes of casing heat-transfer. Experimental measurements were conducted at engine-representative Mach number, Reynolds number and stage inlet to casing wall temperature ratio. Time-resolved heat-transfer coefficient and gas recovery temperature on the casing were measured using an array of heat-transfer gauges. Time-resolved static pressure on the casing wall was measured using Kulite pressure transducers. Time-resolved numerical simulations were undertaken to aid understanding of the mechanism responsible for casing heat load. The results show that between 35% and 60% axial chord the rotor tip-leakage flow is responsible for more than 50% of casing heat transfer. The effects of both gas recovery temperature and heat transfer coefficient were investigated separately and it is shown that an increased stagnation temperature in the rotor tip gap dominates casing heat-transfer. In the tip gap the stagnation temperature is shown to rise above that found at stage inlet (combustor exit) by as much as 35% of stage total temperature drop. The rise in stagnation temperature is caused by an isentropic work input to the tip-leakage fluid by the rotor. The size of this mechanism is investigated by computationally tracking fluid path-lines through the rotor tip gap to understand the unsteady work processes that occur. Copyright © 2005 by ASME.

Unsteady shock motion on a NACA0012 aerofoil at low reduced frequencies

Using transonic blowdown windtunnel experiments, the 2D unsteady shock motion on a NACA0012 aerofoil is examined at various frequencies typical for helicopter blades in forward flight. The aerofoil is subjected to freestream velocities oscillating periodically between M = 0.66 and M = 0.77. Unsteady pressure traces and schlieren images are analyzed over a range of low reduced frequencies to provide information on shock location and strength throughout the cycle. Unsteady effects were noticeable even at very low reduced frequencies (down to O(0.01). However, through the range of frequencies investigated, and within experimental error, the unsteady shock location showed no discernible lag compared to the quasi-steady behaviour. On the other hand, significant variations were observed in shock strengths with the upstream running part of the cycle (decreasing Mach number) displaying considerably stronger shocks than during the accelerating part of the cycle. It could be shown that this variation in shock strength is primarily caused by the shock motion modifying the relative shock Mach number. As a result is was possible to use the quasi-steady results to predict the unsteady shock behaviour at the frequencies investigated here (below 0(0.1)).

Unsteady shock motion on a NACA0012 aerofoil at low reduced frequencies

Using transonic blowdown windtunnel experiments, the 2D unsteady shock motion on a NACA0012 aerofoil is examined at various frequencies typical for helicopter blades in forward flight. The aerofoil is subjected to freestream velocities oscillating periodically between M = 0.66 and M = 0.77. Unsteady pressure traces and schlieren images are analyzed over a range of low reduced frequencies to provide information on shock location and strength throughout the cycle. Unsteady effects were noticeable even at very low reduced frequencies (down to O(0.01). However, through the range of frequencies investigated, and within experimental error, the unsteady shock location showed no discernible lag compared to the quasi-steady behaviour. On the other hand, significant variations were observed in shock strengths with the upstream running part of the cycle (decreasing Mach number) displaying considerably stronger shocks than during the accelerating part of the cycle. It could be shown that this variation in shock strength is primarily caused by the shock motion modifying the relative shock Mach number. As a result is was possible to use the quasi-steady results to predict the unsteady shock behaviour at the frequencies investigated here (below 0(0.1)).