초록
외부 유동소음 문제를 다루는 대부분의 산업현장에서 FW-H(Ffowcs Williams and Hawkings) 방정식을 이용한 복합전산공력음향 기법이 수치적인 효율성으로 인하여 널리 사용되고 있다. 그러나 사중극자항을 무시한 면적분 형태의 FW-H 방정식을 사용할 때 경우에 따라 무시할 수 없는 비물리적인 소음이 발생한다고 알려져 있다. 특히, 수중 프로펠러와 같이 날개 끝 와류 공동(tip vortex cavitation)이 하류방향으로 길게 형성되는 유동에 대해서는 적절하게 모델링하지 않으면 소음 예측의 정확도가 떨어지게 된다. 따라서 본 연구에서는 사중극자 보정항을 추가하여 적분면에서 FW-H 방정식으로부터 발생하는 비물리적인 음향을 저감시키고자 하였다. 먼저 FW-H 방정식에 기초하여 개발한 내부 예측코드의 정확성을 확인하기 위하여 에어컨 실외기에 사용되는 축류팬을 대상으로 검증을 수행하였으며, ANSYS Fluent의 결과와 비교하여 잘 일치하는 것을 확인하였다. 사중극자 보정항의 효과를 확인하기 위하여 등엔트로피 와류 전파에 대한 소음 해석을 수행하였으며, 사중극자 보정항에 의한 오차의 저감 효과가 발생하는 것을 확인하였다. 마지막으로 Clark-Y 수중익형에서 발생하는 공동 유동장을 대상으로 소음 해석을 수행하였으며, 공동이적분면을 통과할 때 발생하는 오차를 사중극자 보정항을 이용하여 저감할 수 있다는 것을 확인하였다.
In most industry fields concerning external flow noise problems, the hybrid computational aeroacoustic techniques based on the FW-H (Ffowcs Williams and Hawkings) equation are widely used for its numerical efficiency. However, when the surface integral form of FW-H equation is used without volume quadrupole sources, it is known to generate significant non-physical noise in a certain case. Especially, in the case of a flow in which the tip vortex cavitation is formed in the distant downstream direction such as flow driven by an underwater propeller, the accuracy in noise prediction becomes poor unless it is not properly modelled. Therefore, in this study, the nonphysical acoustic waves caused by the surface integral form of FW-H equation is reduced by adding the quadrupole correction term. First, to verify the accuracy of the in-house code of FW-H equation, the noise by an axial fan used in the outdoor unit of air conditioner was calculated and compared with the results of ANSYS Fluent. In order to verify the effects of the quadrupole correction term, the noise prediction for isentropic vortex convection is performed and it is confirmed that the error is reduced by the quadrupole correction term. Finally, the noise prediction is performed for the flow field generated by the Clark-Y hydrofoil in underwater. It is confirmed that the error caused by the cavitation passing through the integral surface can be reduced by the quadrupole correction term.