unicfdlab / libAcoustics

libAcoustics - OpenFOAM library for far-field noise computation
140 stars 57 forks source link

On the choice of the surface for FWH Analogy #78

Closed windseacn closed 2 years ago

windseacn commented 2 years ago

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?

unicfdlab commented 2 years ago

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.

windseacn commented 2 years ago

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.

windseacn commented 2 years ago

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.

unicfdlab commented 2 years ago

Could you share your case with square surfaces?

windseacn commented 2 years ago

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

vortexShedm2d.tar.gz

windseacn commented 2 years ago

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?

windseacn commented 2 years ago

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

unicfdlab commented 2 years ago

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

windseacn commented 2 years ago

Thanks for your suggestion, which indeed solves the problem.