Tandem Cylinders: FW-H Against Literature

validation

Animation of the divergence of the velocity (the dilatation) on a plane through the two cylinders, which shows the sound radiated by the flow around and between them, with isosurfaces of the vortical structures coloured by the velocity magnitude.

Overview

Two circular cylinders of equal diameter, placed in tandem along the freestream direction, shed an unsteady turbulent wake whose surface pressure fluctuations radiate sound. This case solves the flow with detached-eddy simulation and predicts the radiated acoustic pressure at one observer using the solid-surface Ffowcs Williams-Hawkings analogy.

Problem Setup

The freestream Mach number is M∞=0.128 (U∞=43.4 m/s, T∞=288.15 K), with a cylinder diameter D=0.05715 m. The domain is periodic in the spanwise direction over 3D, and both cylinder surfaces are no-slip adiabatic walls. The turbulence closure is Spalart-Allmaras delayed detached-eddy simulation, with the FDS convective flux and its hybrid low-dissipation formulation, and second-order implicit time integration at Δt=1.695×10−5 s for 23,600 steps, on a block-structured mesh of 36.8 million hexahedral cells with 96 spanwise layers. A sponge layer that absorbs the outgoing sound starts 1 m from the mesh origin and is 2 m thick. The solid-surface Ffowcs Williams-Hawkings analogy integrates both cylinder walls and records the acoustic pressure at one observer in the mid-span plane, at (9.11D,32.49D) from the centre of the upstream cylinder. This is a microphone position in the Quiet Flow Facility experiment whose spectrum is shown for comparison.

The experimental cylinders span 16D and the simulated span is 3D, so the simulated spectrum is scaled to a 16D span before the comparison, following Seo and Moon (2007), using how coherent the simulated surface pressure is along the span. Away from the shedding tone the pressure is incoherent beyond a fraction of a diameter, and the scaling adds 7.1 dB. Near the tone the pressure stays coherent over the whole simulated span, and the scaling approaches the fully coherent 14.5 dB. The spectrum is averaged over the record from t=0.06 s on, with a frequency resolution of about 12 Hz.

Quantities of Interest

The time-averaged and RMS surface pressure coefficient around both cylinders, the surface pressure spectra at 135∘ on the upstream cylinder and 45∘ on the downstream cylinder, and the time-averaged streamwise velocity and two-dimensional turbulence kinetic energy along the centreline in the gap and in the wake, each compared with the measurements and the CFL3D simulation reported by Lockard (2011). The narrowband power spectral density of the acoustic pressure at the observer, compared with the spectrum measured in the Quiet Flow Facility, and the frequency and level of the spectral maximum between 150 and 220 Hz for both.

Sources

J. H. Seo and Y. J. Moon, “Aerodynamic noise prediction for long-span bodies,” Journal of Sound and Vibration 306, 564–579 (2007)., https://doi.org/10.1016/j.jsv.2007.05.042
D. P. Lockard, “Summary of the Tandem Cylinder Solutions from the Benchmark problems for Airframe Noise Computations-I Workshop,” AIAA Paper 2011-0353, 49th AIAA Aerospace Sciences Meeting, Orlando, FL, January 2011., https://ntrs.nasa.gov/citations/20110005482

Results

Mesh and domain

Mesh and domain
The block-structured mesh in the plane perpendicular to the cylinder axes (383,426 cells in the plane, extruded over 3D with 96 layers, 36.8 million hexahedra): (a) the full domain with the far field at 100D, (b) the inner block, (c) the two cylinders with their O-grid squares and (d) the wall boundary layer of the upstream cylinder, which has 50 layers, a first cell of 10−6 m and a growth ratio of 1.13. The cell size is D/50 in the core, D/25 in the wake and D/10 in the acoustic region.

Mean pressure coefficient on the cylinders

Mean pressure coefficient, upstream cylinder
Mean pressure coefficient, downstream cylinder
Time-averaged pressure coefficient around the upstream (left) and downstream (right) cylinder against the angle from the upstream stagnation point, averaged over three spanwise stations, with NASA's Quiet Flow Facility (QFF) and Basic Aerodynamic Research Tunnel (BART) measurements and NASA's CFL3D simulation, as reported by Lockard (2011). The upstream cylinder is tripped in every measurement; “trip” marks the measurements with the downstream cylinder tripped as well.

RMS pressure coefficient on the cylinders

RMS pressure coefficient, upstream cylinder
RMS pressure coefficient, downstream cylinder
Root-mean-square of the pressure coefficient fluctuation around the upstream (left) and downstream (right) cylinder, averaged over three spanwise stations, with NASA's Quiet Flow Facility (QFF) and Basic Aerodynamic Research Tunnel (BART) measurements and NASA's CFL3D simulation, as reported by Lockard (2011). The upstream cylinder is tripped in every measurement; “trip” marks the measurements with the downstream cylinder tripped as well.

Surface pressure spectra on the cylinders

Surface pressure spectrum, upstream cylinder
Surface pressure spectrum, downstream cylinder
Power spectral density of the surface pressure on the upstream cylinder at θ=135∘ (left) and the downstream cylinder at θ=45∘ (right), in dB re (20μPa)2/Hz, averaged over three spanwise stations, with NASA's Quiet Flow Facility (QFF) and Basic Aerodynamic Research Tunnel (BART) measurements and NASA's CFL3D simulation, as reported by Lockard (2011).

Mean streamwise velocity on the centreline

Mean streamwise velocity in the gap
Mean streamwise velocity in the wake
Time-averaged streamwise velocity on the centreline y=0 in the gap (left) and the wake (right), normalized by the free-stream velocity, averaged over three spanwise stations, with the BART particle image velocimetry with a large field of view (BART LFOV) and with a smaller one (BART), and NASA's CFL3D simulation, as reported by Lockard (2011). The upstream cylinder is tripped in every measurement; “trip” marks the measurements with the downstream cylinder tripped as well.

Two-dimensional turbulence kinetic energy on the centreline

Two-dimensional turbulence kinetic energy in the gap
Two-dimensional turbulence kinetic energy in the wake
Two-dimensional turbulence kinetic energy (u′u′―+v′v′―)/2 on the centreline y=0 in the gap (left) and the wake (right), normalized by the squared free-stream velocity, averaged over three spanwise stations, with the BART particle image velocimetry with a large field of view (BART LFOV) and with a smaller one (BART), and NASA's CFL3D simulation, as reported by Lockard (2011). The upstream cylinder is tripped in every measurement; “trip” marks the measurements with the downstream cylinder tripped as well.

Narrowband spectrum at (9.11D,32.49D)

Narrowband spectrum at \( (9.11D, 32.49D) \)
Power spectral density of the acoustic pressure at the observer, corrected from the simulated spanwise extent of 3D to the experiment's 16D with the frequency-dependent spanwise coherence length of the simulation (Seo and Moon 2007), with the spectrum measured in NASA's Quiet Flow Facility (QFF) and NASA's CFL3D simulation, as reported by Lockard (2011).

Reference: D. P. Lockard, “Summary of the Tandem Cylinder Solutions from the Benchmark problems for Airframe Noise Computations-I Workshop,” AIAA Paper 2011-0353, 49th AIAA Aerospace Sciences Meeting, Orlando, FL, January 2011. (digitized from Fig. 17), https://ntrs.nasa.gov/citations/20110005482

Tonal peak

SourceTonal peak frequency (Hz)Tonal peak level (dB/Hz)
Luminary 173.7393.289
QFF (exp.) 174.6892.522

Frequency and span-corrected level of the spectral maximum between 150 and 220 Hz.

Reference: D. P. Lockard, “Summary of the Tandem Cylinder Solutions from the Benchmark problems for Airframe Noise Computations-I Workshop,” AIAA Paper 2011-0353, 49th AIAA Aerospace Sciences Meeting, Orlando, FL, January 2011. (digitized from Fig. 17), https://ntrs.nasa.gov/citations/20110005482

Discussion

The upstream cylinder sheds vortices that hit the downstream cylinder. The surface pressure on the downstream cylinder fluctuates several times more strongly than on the upstream one (compare the two RMS pressure plots), which makes it the louder of the two. The sound is a strong tone at the shedding frequency with a broad hiss on either side.

Most of the page compares the flow, because the sound depends on the flow that makes it. On the upstream cylinder the mean and RMS (root-mean-square, a measure of how strongly the pressure fluctuates) surface pressure follow the Quiet Flow Facility (QFF) measurements. On the downstream cylinder the pressure fluctuations peak at the same angles as measured, and the surface pressure spectrum follows the measurement within a few decibels from 200 Hz to 2 kHz. In the gap between the cylinders the turbulence kinetic energy follows the measurement, and behind the downstream cylinder the flow follows the Basic Aerodynamic Research Tunnel (BART) measurements made with both cylinders tripped.

The narrowband spectrum compares the predicted sound at the microphone with the QFF measurement. The simulation covers 3 diameters of span and the experiment 16, and a longer cylinder radiates more sound, so the simulated spectrum is scaled up to the experiment's span before the comparison, as described in the problem setup. The tone is at 173.7 Hz against the measured 174.7 Hz, a difference smaller than the 12 Hz frequency resolution of the spectrum, and its level is within 0.8 dB of the measurement. From about 300 Hz to 1.5 kHz the computed spectrum is within about 4 dB of the measured broadband level, including the smaller peak near 530 Hz, and it continues to follow the measured decline up to 3 kHz, up to 10 dB lower.

Solver configuration

features exercised
Validation3DDES/LESSpalart–AllmarasIdeal gasUnsteadyEnergy equationFW-H (solid / permeable surface)