Aeolian Tone of a Circular Cylinder: FW-H Against Direct CFD

verification

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 cylinder in a cross-flow sheds vortices alternately from its two sides; the resulting lift fluctuation radiates a tone, the Aeolian tone, whose frequency is set by the flow and whose angular pattern is that of a dipole aligned with the lift direction. This case computes the flow with an unsteady Reynolds-averaged simulation, predicts the radiated tone with the Ffowcs Williams-Hawkings analogy evaluated both on the cylinder wall and on a closed permeable surface in the same solution, and compares it with the pressure the flow solver computes directly at the same observer points.

Problem Setup

The flow is at M∞=0.25 (U∞=85.07 m/s, T∞=288.15 K, p∞=86985 Pa) over a cylinder of diameter D. The domain is periodic in the spanwise direction over 4D, and a sponge layer that absorbs the outgoing sound starts at 25D and is 12D thick. The flow is computed with an unsteady Reynolds-averaged simulation using the Spalart-Allmaras turbulence model, the fluid is an ideal gas with Sutherland viscosity, and time integration is second-order implicit with Δt=5.88×10−4 s. The permeable surface is a box spanning −2D≤x≤10D and −2D≤y≤2D over the full span, and both Ffowcs Williams-Hawkings integrals are repeated over 25 copies of the span to represent a long cylinder. There are 48 observers at mid-span: six microphones between 3D and 12D from the cylinder, and a ring of radius 12D with extra observers around the two quiet directions along the flow axis.

Quantities of Interest

The level of the shedding tone at the ring observers, from both Ffowcs Williams-Hawkings surfaces and from the directly simulated pressure, over the window 120≤tU∞/D≤150, and the pressure time histories at the six near-field microphones, from the three sources.

Results

Mesh and domain

Mesh and domain
The mesh in the plane perpendicular to the cylinder axis, with the permeable surface (green box), the start (r=25D) and the full strength (r=37D) of the sponge layer, the far-field boundary, the 25 observers on a ring of radius 12D and the reference microphones.

Tone level on the 12D observer ring

Tone level on the \( 12D \) observer ring
Level of the Aeolian tone at the observers on a ring of radius 12D, from the Ffowcs Williams-Hawkings analogy on the cylinder wall and on the permeable surface and from the pressure the solver samples at the same points. The angle is measured from the flow direction about the cylinder axis.

Microphone pressure at 3D, 6D and 12D

Microphone pressure at \( 3D \), \( 6D \) and \( 12D \)
Pressure at the three near-field microphones at 90∘ above the cylinder, at distances of 3D, 6D and 12D, over the analysis window, from the pressure the solver samples at the same points and from the cylinder-wall and permeable FW-H surfaces. Time is in convective units tU∞/D.

Microphone pressure at 45∘, 90∘ and 135∘

Microphone pressure at \( 45^\circ \), \( 90^\circ \) and \( 135^\circ \)
Pressure at the three near-field microphones at 5D from the cylinder axis, at 45∘, 90∘ and 135∘, over the analysis window, from the pressure the solver samples at the same points and from the cylinder-wall and permeable FW-H surfaces. Time is in convective units tU∞/D.

Discussion

The Aeolian tone is the hum a wire or a pole makes in the wind. Vortices shed alternately from the two sides of the cylinder push it up and down, so it radiates like a dipole: loud above and below the cylinder, quiet along the flow direction. The plot of the tone level on the 12D ring shows this pattern, and the Ffowcs Williams-Hawkings (FW-H) predictions from the cylinder wall and from the permeable surface lie on top of the pressure the solver samples at the same points.

Some difference between FW-H and direct CFD (computational fluid dynamics) is expected. The FW-H method takes the flow solution on a surface close to the cylinder and carries the sound from there to each observer with a formula that is exact for sound travelling through a uniform flow, so nothing is lost on the way. The direct CFD value is the wave after the solver itself has carried it across the grid, and every cell it passes through takes a little energy out of it. Because of this numerical dissipation, the directly computed level tends to sit at or slightly below the FW-H level, and the gap grows with distance.

Directly downstream, near 0∘, the ring observers sit in the wake. There the pressure the solver samples includes the passing vortices themselves and not only sound, and the direct CFD points lie above the FW-H curves.

The two microphone plots compare the pressure signals themselves over the analysis window, not just their levels. The tone rises and falls in strength over the window, and both FW-H predictions follow the solver's pressure through it, in shape and timing, at all six microphones. The permeable-surface prediction also matches its level at every microphone. The cylinder-wall prediction is shifted slightly at the two microphones closest to the shed vortices: 3D above the cylinder, and 5D away at 45∘ on the downstream side. This is expected. The cylinder-wall surface counts only the sources on the cylinder itself, while close to the cylinder and its wake the pressure the solver samples also includes the pressure of the vortices shed into the flow. The permeable surface encloses the cylinder and the near wake, so it accounts for them. Further out, or upstream, the pressure of the vortices has fallen off, since it decays with distance faster than sound does, and the two surfaces agree.

Solver configuration

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