Frequency Domain Calculations Of Acoustic Propagation

Two complex geometry problems are solved using the linearized Euler equations. The impedance mismatch method1 is used to impose the solid surfaces without the need to use a body-fitted grid. The problem is solved in the frequency domain to avoid long run times. Although the harmonic assumption eliminates all time dependence, a pseudo-time term is added to allow conventional iterative methods to be employed. A Jameson type, Runge-Kutta scheme is used to advance the solution in pseudo time. The spatial operator is based on a seven-point, sixth-order finite difference. Constant coefficient, sixth-derivative artificial dissipation is used throughout the domain. A buffer zone technique employing a complex frequency to damp all waves near the boundaries is used to minimize reflections. The results show that the method is capable of capturing the salient features of the scattering, but an excessive number of grid points are required to resolve the phenomena in the vicinity of the solid bodies because the wavelength of the acoustics is relatively short compared with the size of the bodies. Smoothly transitioning into the immersed boundary condition alleviates the difficulties, but a fine mesh is still required.