Cylinder at Re = 10⁴: FW-H Against Literature

verification

Animation of the divergence of the velocity (the dilatation) on a plane through the cylinder, which shows the sound radiated by the turbulent wake, with isosurfaces of the vortical structures coloured by the velocity magnitude in m/s.

Overview

A circular cylinder at a Reynolds number of 104 sheds a turbulent wake that radiates sound at the shedding frequency. This case simulates the flow with a large-eddy simulation and predicts the far-field sound with the Ffowcs Williams-Hawkings analogy, evaluated on the cylinder wall and on a permeable surface around the cylinder in the same flow solution. The wake statistics, the far-field spectra and the directivity are compared with Khalighi et al. (2010).

Problem Setup

The flow is at M∞=0.2 (U∞=68.06 m/s, T∞=288.15 K, p∞=86985 Pa) with a constant dynamic viscosity chosen so that ReD=104 for a cylinder of diameter D. The O-grid is πD deep and periodic in the spanwise direction, and a sponge layer that absorbs the outgoing sound starts at 26D and is 22D thick. The large-eddy simulation uses the WALE subgrid-scale model and the FDS convective flux with its hybrid low-dissipation formulation, with second-order implicit time integration and Δt=6.06×10−5 s.

Starting from a fully developed flow, the statistics and the acoustic record are collected over 86,368 time steps. The acoustic pressure is recorded at 38 observers at mid-span: a ring of 36, 10∘ apart, of radius 16.8D about x=3.26D, and the two microphones of the reference at (−1.2D,16.2D) and (−13.5D,0). The permeable surface is the box −6D≤x≤25D, −6D≤y≤6D over the full span, open on its spanwise faces and on its downstream face. Both Ffowcs Williams-Hawkings integrals are repeated over 25 copies of the span to represent a long cylinder.

Quantities of Interest

The mean and RMS streamwise velocity and the Lighthill-stress spectra in the wake, the far-field pressure spectra at the two microphones, and the level of the sound at the shedding frequency around the ring, from both integration surfaces, each compared with Khalighi et al. (2010).

Sources

Y. Khalighi, A. Mani, F. Ham and P. Moin, “Prediction of Sound Generated by Complex Flows at Low Mach Numbers,” AIAA Journal 48(2), 306–316 (2010)., https://doi.org/10.2514/1.42583

Results

Mesh and domain

Mesh and domain
The O-grid mesh in the plane perpendicular to the cylinder axis, with the permeable surface (green box), the start (r=26D) and the full strength (r=48D) of the sponge layer, the 36 observers of the directivity arc at r=16.8D and the two microphones of the spectra. Mesh lines are drawn one in four; the solver mesh is finer.

Velocity profiles in the wake

Velocity profiles in the wake
Mean and RMS streamwise velocity across the wake at x/D=2 and 5, with Khalighi et al. (2010, Fig. 7).

Lighthill stress spectra

Lighthill stress spectra
Spectra of the Lighthill stress components u1′u1′ and u2′u2′ at x/D=5 and 15, with Khalighi et al. (2010, Fig. 8).

Far-field pressure spectra

Far-field pressure spectra
Far-field pressure spectra at the two microphones of Khalighi et al. (2010, Fig. 9), from the FW-H analogy on the cylinder wall and on the permeable surface, with the reference curves, plotted against fD/(StU∞).

Far-field directivity at the shedding frequency

Far-field directivity at the shedding frequency
Level of the far-field pressure at the shedding frequency around the observer ring on a half-polar plot, one panel per integration surface, normalized by its maximum, with Khalighi et al. (2010, Fig. 10).

Discussion

At a Reynolds number of 104 the wake is turbulent, so the sound is a peak at the shedding frequency on top of a broad hiss from the turbulence. The sound depends on the flow that makes it, so the page first compares the wake. The mean and RMS (root-mean-square, a measure of how strongly the velocity fluctuates) velocity profiles behind the cylinder follow Khalighi et al. at both stations, and the Lighthill-stress spectra, which describe the turbulent fluctuations that act as sound sources, peak at the same frequency as the reference with the same level at low frequency.

The far-field spectra are plotted against frequency divided by the shedding frequency, so the shedding tone is at 1, and each spectrum is scaled to its own tone. Above the cylinder (105∘ from downstream) the strongest peak is at 1. This is the lift, which swings up and down once per shedding cycle. Directly upstream (180∘) the strongest peak is at 2. This is the drag, which pulses twice per cycle, once for each vortex. Both Ffowcs Williams-Hawkings (FW-H) surfaces reproduce these peaks, the smaller harmonics, and the fall of the spectrum at higher frequencies. The largest differences are at the lowest frequencies, well below the tone.

The directivity plot shows how the strength of the tone at the shedding frequency varies with direction, each curve scaled to its own maximum. Both surfaces give the single broad lobe of the reference, loud above the cylinder and quiet along the flow axis, with the computed lobe leaning further upstream.

Solver configuration

features exercised
Verification3DDES/LESIdeal gasUnsteadyEnergy equationFW-H (solid / permeable surface)