Laminar Circular Cylinder at Re = 100: FW-H Against Literature

validation

Animation of the divergence of the velocity (the dilatation) on a plane through the cylinder, which shows the sound radiated by the vortex shedding, with isosurfaces of the vortical structures in the wake.

Overview

A circular cylinder in a laminar cross-flow sheds an alternating vortex street whose surface-pressure fluctuations radiate a tone at the shedding frequency. This case predicts the far-field sound of a Reynolds-number-100, Mach-0.15 cylinder with the Ffowcs Williams-Hawkings analogy, evaluated on the cylinder wall and on a closed permeable surface in the same flow solution, and compares the angular distribution of the shedding tone with a published reference.

Problem Setup

The flow is laminar with M∞=0.15 (U∞=51.04 m/s, T∞=288.15 K, p∞=86985 Pa) and a constant dynamic viscosity chosen so that ReD=100 for a cylinder of diameter D. The O-grid is one diameter deep and periodic in the spanwise direction, and a sponge layer that absorbs the outgoing sound starts at 16D and is 18D thick. Time integration is second-order implicit with Δt=2.5×10−4 s.

The Ffowcs Williams-Hawkings analogy is evaluated on the cylinder wall and on a permeable box spanning −8D≤x≤15D and −5D≤y≤5D, open on its spanwise faces. Both integrals are repeated over 50 copies of the span to represent a long cylinder. The acoustic pressure is recorded at 36 observers spaced 10∘ apart on an arc of radius 12.9D about x=1.86D, at mid-span.

Quantities of Interest

The frequency of the shedding tone, and its level around the observer arc from both integration surfaces, normalized by its maximum between 25∘ and 175∘ and compared with the reference curve, with the RMS difference between the two over that arc.

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 mesh in the plane perpendicular to the cylinder axis, every cell edge drawn, with the permeable surface (green box), the start and the full strength of the sponge layer, and the 36 observers on a ring of radius 12.9D; the light blue arc is the scored angle range, 25∘ to 175∘.

Far-field directivity of the shedding tone

Far-field directivity of the shedding tone
Level of the shedding tone around the observer arc, relative to its maximum over 25∘ to 175∘, from the cylinder-wall and permeable FW-H surfaces, with the reference curve.

Reference: 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). (digitized far-field directivity of the laminar cylinder), https://doi.org/10.2514/1.42583

Directivity: difference from the reference curve

FW-H surfaceShedding tone (Hz)RMS difference from the reference curve (dB)
Solid 8.19460.9304
Permeable 8.3290.98396

Shedding-tone frequency and the RMS difference between the normalized directivity levels and the reference curve over the same arc.

Reference: 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). (digitized far-field directivity of the laminar cylinder), https://doi.org/10.2514/1.42583

Discussion

At this low Reynolds number the wake is a regular street of vortices that leave the top and the bottom of the cylinder in turn. Each one pushes the cylinder sideways, so the cylinder radiates a single clean tone, here at 8.2 Hz from the cylinder-wall surface and 8.3 Hz from the permeable surface. The sound has a dipole pattern: loud above and below the cylinder, quiet along the flow direction. In the directivity plot that shows up as two lobes, one up and one down.

The plot compares the shape of this pattern, with each curve divided by its own maximum between 25∘ and 175∘. Both Ffowcs Williams-Hawkings (FW-H) surfaces reproduce the reference: the two lobes, their width, and the quiet directions along the flow axis. Over the scored arc the root-mean-square difference from the reference is 0.93 dB for the cylinder wall and 0.98 dB for the permeable surface, about 11% in pressure amplitude.

Solver configuration

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