unicfdlab / libAcoustics

libAcoustics - OpenFOAM library for far-field noise computation

Repository from Github https://github.comunicfdlab/libAcousticsRepository from Github https://github.comunicfdlab/libAcoustics

On the choice of the surface for FWH Analogy

windseacn opened this issue · comments

Hello, I'm interested in CAA and find this library is very useful for beginners. I tested the vortexshed tutorial with different surfaces, the cylinder surface (default in the tutorial), a square box with the side length equal to the diameter of the cylinder and one larger cylindrical surface with its radius equal to the cylinder diameter, these three surfaces are shown below:
surface

The sound pressure levels for these three cases are shown below.

untitled

From my limited knowledge on this subject, the results from these three surface configurations should be very close as the Mach number is very low. However, the result with the square surface is quite different from the others, there exists a second peak in the spectrum. While the result from larger cylindrical surface is close that from rigid body surface. Does this mean the results are sensitive to the shape of the surface chosen?

Hello,

I think that, the problem with square is related to the fact that the square surface intersects the cylinder patch which caused the problems with data interpolation. Try to use square with side length equal to 4R.
The difference between patch surface and 2R circle caused by the influence of the vortex wake.

Regards.

Thanks for your reply. The difference between patch surface and 2R circle is minor, and could be expected.

I have tried square with side length equal to 3R, which is not good either as below.
untitled1

I will try whether the result for square with side length equal to 4R and report later.

For four different surfaces as below
surface2

The SPLs for these four configurations are
untitled2

It's interesting to note that the frequency of the maximum SPL changes if a larger square surface is used.

Could you share your case with square surfaces?

The case is the same as the vortexShed tutorial from the library, with more surfaces added.

vortexShedm2d.tar.gz

I have tested more surfaces, some low resolution cylinder surfaces with different segments and radius 2R
surface4

The SPL results are shown below
untitled4
It can be seen that for these surfaces, there exist additional SPL peaks around double and quadruple base frequencies. As the number of segments increases, i.e., the cylinder becomes smoother, the contributions from higher frequencies become lower. I wonder whether this is physical result, or induced by some numerical issues?

Cai et al, have tested different square boxes for flow around cylinder with FWH analogy, from which this phenomenon is not observed.
surface3


A numerical study of transient flow around a cylinder and aerodynamic sound radiation
https://link.springer.com/article/10.1134/S0869864318030022

The FW-H stl surfaces have poor resolution (2 cells per square side), try to create and use square surface with more cells per side.

Regards

Thanks for your suggestion, which indeed solves the problem.