전산음향학

편집 역사 토론
시뮬레이션 전산유체역학 물리 마지막 수정: 2026-07-16 04:38:19

1. 개요[편집]

전산음향학
Computational Aeroacoustics
약칭CAA
이론적 기초라이트힐 음향 상사론 (1952)
표준 적분식FW-H 방정식 (1969)
핵심 난제에너지 스케일 차, 수치 분산
주 사용처항공기 소음, 팬·블로워, 자동차 풍절음

유동 에너지의 0.00001%가 소리로 새어 나간다. CFD 엔지니어에게는 반올림 오차, 소음 엔지니어에게는 밥줄.

전산음향학(Computational Aeroacoustics, CAA)은 유동이 만들어내는 소리의 발생과 전파를 수치적으로 예측하는 분야다. 정확히는 “공력음향학”의 계산 버전으로, 스피커 진동 같은 구조 기인 소음이 아니라 유동 자체가 음원인 문제를 다룬다. 제트 소음, 프로펠러 소음, 자동차 사이드미러 풍절음, 데이터센터 팬 소리 — 전부 여기 소관.

CAA가 그냥 CFD의 부록이 아니라 독립 분야가 된 이유는 스케일 때문이다. 아음속 유동에서 음향으로 변환되는 에너지는 유동 운동에너지의 M5M^5 정도 — 마하 0.2면 대략 10−410^{-4}의 0.3배 수준이다.1 압력장으로 보면, 수백 Pa의 유동 압력 변동 위에 얹힌 수 mPa짜리 음향 요동을 골라내야 한다. 5~7자릿수 아래의 신호를 잡아야 하는데, 보통 CFD 도식의 수치 소산이 그 신호보다 크다. 답이 오차보다 작은 상황 — 이것이 CAA의 출발점이자 영원한 고통이다.

2. 음향 파동 방정식[편집]

정지 매질에서 작은 요동에 대한 지배식은 익숙한 파동 방정식이다.

∂2p′∂t2−c02∇2p′=0\frac{\partial^2 p'}{\partial t^2} - c_0^2 \nabla^2 p' = 0

이건 그냥 선형화된 오일러 방정식에서 유도된다. 즉 음향은 압축성 유동의 부분집합이다. 비압축 CFD가 소리를 못 내는 이유가 여기 있다 — 음속을 무한대로 놓는 순간 파동이 사라진다.

시간 조화 가정(p′=p^e−iωtp' = \hat{p}e^{-i\omega t})을 넣으면 헬름홀츠 방정식이 된다.

∇2p^+k2p^=0,k=ω/c0\nabla^2 \hat{p} + k^2 \hat{p} = 0, \qquad k = \omega / c_0

이 형태는 맥스웰 방정식 기반 전자기 해석과 수학적으로 쌍둥이다. 그래서 경계요소법이나 완전정합층 같은 도구가 두 분야를 그대로 오간다. 소음 엔지니어와 안테나 엔지니어가 같은 논문을 읽는 이유.

문제는 유동이 있으면 이 깔끔함이 무너진다는 것. 대류, 굴절, 산란이 끼어들고, 비선형 항이 음원 역할을 한다. 그 지저분함을 정리한 것이 라이트힐이다.

3. 라이트힐 음향 상사론[편집]

1952년 제임스 라이트힐(M. J. Lighthill)은 제트 엔진 소음 문제를 붙들고 있다가 기막힌 재배열을 해냈다. 나비에-스토크스 방정식과 연속 방정식을 조합해 정확한 항등식을 얻은 것이다.

∂2ρ′∂t2−c02∇2ρ′=∂2Tij∂xi∂xj\frac{\partial^2 \rho'}{\partial t^2} - c_0^2 \nabla^2 \rho' = \frac{\partial^2 T_{ij}}{\partial x_i \partial x_j}

좌변은 정지 매질의 파동 연산자, 우변이 음원이다. 라이트힐 응력 텐서는 다음과 같다.

Tij=ρuiuj+(p′−c02ρ′)δij−τijT_{ij} = \rho u_i u_j + (p' - c_0^2 \rho')\delta_{ij} - \tau_{ij}

여기서 강조할 게 있다. 이건 근사가 아니라 항등식이다. 나비에-스토크스를 다시 쓴 것에 불과하다. 근사가 되는 순간은 우변 TijT_{ij}를 유동 해석에서 계산해 넣고, 그 음원이 되돌아와 유동에 영향을 주지 않는다고 가정할 때다. 그래서 “상사론(analogy)” — 실제 유동을 정지 매질 속 등가 음원의 분포로 비유한다는 뜻이다.2

이 이론의 최대 성과는 M8M^8 법칙이다. TijT_{ij}의 이중 발산은 사중극자(quadrupole) 음원이고, 차원 해석을 하면 아음속 난류 제트의 음향 파워가 제트 속도의 8제곱에 비례한다는 결론이 나온다.

Pac∼ρ0U8D2c05P_{\text{ac}} \sim \rho_0 \frac{U^8 D^2}{c_0^5}

속도를 반으로 줄이면 소음 파워가 256분의 1. 초창기 제트기가 귀를 찢던 것과, 현대 고바이패스 터보팬이 배기 속도를 낮춰 조용해진 것이 이 지수 하나로 설명된다. 이론 하나가 항공기 엔진 형상을 바꿔버린 드문 사례.3

이후 커얼(Curle)이 고정 벽면 효과(쌍극자)를, Ffowcs Williams와 Hawkings가 움직이는 벽면을 확장했다. 음원 세기 순서는 단극자 > 쌍극자 > 사중극자 — 그래서 저마하 유동에서 벽이 있으면 벽면 쌍극자가 지배한다. 사이드미러 소음이 대표적.

4. FW-H 방정식[편집]

Ffowcs Williams–Hawkings 방정식(1969)은 CAA 실무의 진짜 주인공이다. 임의로 움직이는 표면이 있는 유동에 대해 라이트힐 상사론을 일반화한 것으로, 일반화 함수를 써서 다음 형태로 정리된다.

□2p′=∂∂t[Qδ(f)]−∂∂xi[Liδ(f)]+∂2∂xi∂xj[TijH(f)]\Box^2 p' = \frac{\partial}{\partial t}\left[ Q \delta(f) \right] - \frac{\partial}{\partial x_i}\left[ L_i \delta(f) \right] + \frac{\partial^2}{\partial x_i \partial x_j}\left[ T_{ij} H(f) \right]

우변 세 항이 각각 두께 소음(단극자 — 물체가 공기를 밀어내며 생기는), 하중 소음(쌍극자 — 표면 압력 변동), 사중극자 소음(체적 난류)이다. 헬리콥터 로터에서 두께 소음과 하중 소음이 왜 따로 논의되는지가 여기서 나온다.

실무에서 FW-H가 사랑받는 이유는 명확하다. 적분 해로 풀리기 때문이다. 유동 해석은 물체 근처 좁은 영역만 하고, 그 표면(투과성 FW-H 면을 쓰면 물체를 감싼 가상의 면도 가능)에서 압력·속도 이력을 저장한 뒤, 지연 시간(retarded time)을 고려한 표면 적분으로 수백 미터 밖 관측점 소음을 계산한다. 격자를 관측점까지 깔 필요가 없다. ANSYS Fluent, STAR-CCM+, OpenFOAM 모두 FW-H 모듈을 갖고 있다.

투과성 면 위치 선정은 CAA 실무의 대표적 지뢰밭이다. 면이 너무 안쪽이면 사중극자를 놓치고, 너무 바깥이면 격자가 조대해져 그 면을 지나는 음파가 이미 뭉개져 있다. 게다가 면이 후류(wake)를 가로지르면 와류가 통과하며 가짜 음원을 만든다. 그래서 여러 개의 면을 놓고 결과를 평균하거나, 후류 방향 끝면을 열어두는 편법이 관행이다.

5. 직접 계산 vs 하이브리드[편집]

CAA 접근법은 크게 두 갈래다.

방식격자 범위비용주 사용처
직접 계산 (DNC)음원 + 원거리 전파 전체살인적연구, 초음속 제트
하이브리드음원 영역만 + 적분/전파현실적산업계 대부분

직접 계산은 압축성 나비에-스토크스를 음원부터 관측점까지 통째로 푼다. 원리적으로 완벽하지만 문제가 셋이다. (1) 파장이 물체 크기보다 훨씬 길어서 도메인이 터무니없이 커야 한다. (2) 시간 스텝이 음속 CFL 조건에 묶인다. (3) 앞서 말한 진폭 스케일 차. 항공기 전기체 직접 계산은 DNS와 비슷한 팔자다.

하이브리드는 음원 생성과 전파를 분리한다.

  1. 음원 영역만 대와류 모사나 분리 와류 모사로 푼다. RANS는 시간 평균이라 원리적으로 광대역 소음을 못 낸다.4
  2. 표면/체적 음원 데이터를 저장한다.
  3. FW-H 적분이나 선형화 오일러 방정식(LEE), 음향 유한요소법/경계요소법으로 원거리까지 전파시킨다.

전파 단계가 선형이라 비용이 극적으로 떨어진다. 대신 단방향 가정 — 음향이 다시 유동을 흔드는 피드백이 있으면(휘슬, 헬름홀츠 공명, 캐비티 톤) 하이브리드가 놓친다. 이럴 땐 직접 계산 외에 답이 없다.

6. 수치 요구 조건[편집]

CAA 도식이 일반 CFD 도식과 다른 이유는 파동을 수백 파장 거리까지 형태를 유지한 채 날려야 하기 때문이다.

  • 저분산·저소산 고차 도식: 2차 상류차분은 CAA에서 사형이다. 소산으로 음파를 몇 파장 만에 지워버린다. 대신 6차 이상 콤팩트 차분, DRP(Dispersion-Relation-Preserving) 도식, 스펙트럴 방법 계열, 고차 유한체적법/DG를 쓴다. CFD의 국룰인 “2차면 충분”이 CAA에서는 통하지 않는다.
  • 파장당 격자점(PPW): 2차 도식은 파장당 2030점을 요구하지만, 고차 도식은 68점이면 된다. 3차원에서 이 차이는 격자 수 수십 배다. 고차 도식이 비싸 보여도 결국 싸게 먹히는 이유.
  • 시간 적분: 룽게-쿠타법 계열의 저장 최적화 LDDRK가 표준. 위상 오차가 작아야 한다.
  • 무반사 경계: 도메인 끝에서 음파가 튕겨 돌아오면 그 반사파가 신호와 같은 크기다. 완전정합층(PML), 특성 기반 비반사 경계(NSCBC), 격자를 일부러 조대하게 만들어 음파를 죽이는 스펀지 존(sponge zone)이 쓰인다. 스펀지 존은 “수치 소산을 무기로 쓴다”는 점에서 우아하지는 않지만 잘 먹힌다.

7. 활용 사례[편집]

  • 항공기: 이착륙 소음 인증(ICAO Chapter 14)이 법적 규제라 예측이 곧 돈이다. 팬 소음, 제트 소음, 랜딩기어·플랩의 기체 소음(airframe noise). 요즘은 UAM/드론의 프로펠러 소음이 뜨거운 주제다.
  • 팬·블로워: 데이터센터, 냉장고, 노트북. 블레이드 통과 주파수(BPF) 톤과 광대역 소음을 나눠 본다. 톤 소음은 로터-스테이터 간섭에서 오므로 터보기계 해석과 붙어 있다.
  • 자동차: 시속 100 km에서 사이드미러·A필러 풍절음이 엔진 소음을 이긴다. 전기차 시대에 이게 더 중요해졌다 — 엔진 소리가 사라지자 그동안 가려져 있던 것들이 다 들리기 시작했다.
  • HVAC: 덕트 소음, 팽창밸브 소음.

8. 관련 문서[편집]

9. Footnotes[편집]

  1. 라이트힐의 M8M^8 파워 법칙을 운동에너지 플럭스(∼ρU3D2\sim \rho U^3 D^2)로 나누면 음향 효율이 M5M^5에 비례한다는 결과가 나온다. 마하 0.2면 0.25=3.2×10−40.2^5 = 3.2 \times 10^{-4} 스케일. 계수까지 따지면 실제로는 이보다 더 작다. 제트기 옆에서 고막이 찢어지는데 그게 에너지의 0.01%도 안 된다는 게 음향의 잔인함이다. ↩

  2. “상사론”이라는 이름 때문에 근사 이론으로 오해받는데, 방정식 자체는 완전히 엄밀하다. 근사는 오직 TijT_{ij}를 어떻게 조달하느냐에서 들어온다. 라이트힐 본인도 “나는 아무것도 근사하지 않았다, 다만 모르는 것을 우변에 몰아넣었을 뿐”이라는 취지의 말을 남겼다. 미지의 항을 우변으로 옮기는 것 — 물리학의 유서 깊은 전통이다. ↩

  3. 초기 터보젯의 배기 속도는 600 m/s 수준이었고, 현대 고바이패스 터보팬은 훨씬 느린 대량의 공기를 민다. 같은 추력을 내면서 속도를 낮춘 것인데, U8U^8이라 소음이 극적으로 줄었다. 물론 대신 팬이 커지면서 팬 소음이 새 주인공이 되었다. 소음 공학은 두더지 잡기다. ↩

  4. 정확히는 RANS로도 통계적 소음 모델(SNGR, 확률적 음원 재구성)을 붙이면 광대역 예측 시늉은 가능하다. 다만 “시늉”이라는 단어를 강조하고 싶다. 진짜 필요한 건 시간 해상 난류이고, 그래서 CAA 예산의 대부분은 LES 돌리는 데 들어간다. ↩