23.47 격자 볼츠만 방법(LBM)을 이용한 프로펠러 해석

23.47 격자 볼츠만 방법(LBM)을 이용한 프로펠러 해석

1. LBM의 개념

격자 볼츠만 방법(Lattice Boltzmann Method, LBM)은 유체의 거시적 거동을 통계역학 기반의 볼츠만 방정식을 이산화하여 해결하는 수치 기법이다. 전통적 Navier-Stokes 방정식 기반 CFD가 유체의 거시적 보존 법칙을 직접 해결하는 것과 달리, LBM은 입자 분포 함수(particle distribution function)의 시간 발전을 격자 상에서 추적하여 유동장을 계산한다. Chen과 Doolen이 Lattice Boltzmann Method for Fluid Flows(Annual Review of Fluid Mechanics, vol. 30, 1998)에서 체계적으로 정리한 바 있다.

2. LBM의 기본 방정식

LBM의 지배 방정식은 다음과 같은 이산 볼츠만 방정식이다.

f_i(\mathbf{x} + \mathbf{e}_i \Delta t, t + \Delta t) = f_i(\mathbf{x}, t) + \Omega_i(f)

여기서 f_i는 속도 방향 \mathbf{e}_i의 입자 분포 함수, \Omega_i는 충돌 연산자이다. 가장 널리 사용되는 충돌 연산자는 Bhatnagar-Gross-Krook(BGK) 충돌 모형이다.

\Omega_i^{\text{BGK}} = -\dfrac{1}{\tau} (f_i - f_i^{\text{eq}})

여기서 \tau는 이완 시간, f_i^{\text{eq}}는 평형 분포 함수이다. 거시적 밀도 \rho와 속도 \mathbf{u}f_i의 0차와 1차 모멘트로 계산된다.

\rho = \sum_i f_i, \quad \rho \mathbf{u} = \sum_i \mathbf{e}_i f_i

적절한 Chapman-Enskog 전개에 따르면, LBM은 저 마하 수 영역에서 Navier-Stokes 방정식을 재현한다.

3. LBM의 장점

LBM은 다음과 같은 장점을 가진다. 첫째, 격자의 국지적 계산으로 인해 병렬화가 용이하며, GPU 기반 가속이 효과적이다. 둘째, 복잡한 경계 형상을 격자 외부에 있는 노드로 처리하는 간단한 방법(bounce-back method)이 가능하다. 셋째, 비압축성 유동에 대해서도 명시적 시간 적분으로 해결되므로 안정성이 높다. 넷째, 다상 유체, 다성분 유체, 복잡한 경계 조건의 확장이 비교적 간단하다. 다섯째, 명시적 격자 이동 없이 복잡한 경계 운동을 처리할 수 있다.

4. LBM의 한계

LBM은 다음과 같은 한계도 가진다. 첫째, 전통적 D3Q19, D3Q27 등의 격자 구조는 마하 수가 약 0.3 이상인 영역에서 정확도가 감소한다. 둘째, 압축성 유동과 충격파 문제의 직접 해석에는 특수 확장이 필요하다. 셋째, 메모리 요구량이 Navier-Stokes 기반 CFD에 비해 상대적으로 높다. 넷째, 벽면 근처 경계층 해상을 위해서는 적절한 격자 밀집화와 벽 모형이 필요하다. 이러한 한계에도 LBM은 저·중속 비압축성 유동, 복잡한 경계, 비정상 유동의 해석에 우수한 성능을 보인다.

5. 프로펠러 해석에의 적용

프로펠러 CFD 해석에 LBM을 적용하는 접근은 다음과 같다. 첫째, 회전하는 블레이드를 이동 경계(moving boundary)로 처리하는 bounce-back 방법 또는 imbedded boundary method를 활용한다. 둘째, 격자를 블레이드 주변에 세밀하게 배치하고, 후류 영역은 점진적으로 확장한다. 셋째, 난류 해석에는 Large Eddy Simulation(LES) 기반 LBM이 사용된다. 넷째, 소음 예측을 위해 LBM-FW-H 결합 접근이 활용된다.

6. 대표 연구

연구내용
Romani, Casalino프로펠러 및 로터의 LBM-LES 해석
Avallone 등프로펠러 소음 예측의 LBM 적용
Casalino 외eVTOL 기체의 LBM 공력 음향 해석
Casalino, Lockard회전익 기체의 LBM 비정상 해석

이 표는 프로펠러와 로터의 LBM 해석에 대한 대표적 최근 연구를 예시한 것이다. LBM은 특히 공력 음향(aeroacoustics) 해석에서 Navier-Stokes 기반 CFD 대비 효율적 접근으로 각광받고 있다.

7. 공력 음향 해석에의 활용

LBM은 다음의 이유로 공력 음향 해석에 적합하다. 첫째, 명시적 시간 적분으로 충분한 시간 해상도를 제공하여, 소리 전파에 필요한 정밀 압력 변동을 재현할 수 있다. 둘째, 낮은 수치 확산으로 원거리 음장의 음압 예측 정확도가 높다. 셋째, 복잡한 경계를 포함한 기체 구조의 회절과 산란을 직접 재현할 수 있다. 넷째, GPU 가속을 통해 고해상도 해석이 실용 시간 내에 가능하다. 이러한 특성은 도심항공교통 기체의 저소음 설계 연구에서 LBM의 활용이 급증하는 이유이다.

8. 상용 및 공개 LBM 코드

LBM 기반 상용 코드로는 PowerFLOW(Simulia, 현재 Dassault Systèmes)가 대표적이며, 자동차 산업과 항공 공력 음향 해석에 널리 활용되어 왔다. 공개 소스 코드로는 OpenLB, Palabos, waLBerla 등이 존재한다. 이들 코드는 각기 다른 격자 구조, 병렬화 전략, 물리 모형을 제공한다.

9. 해석 사례와 검증

LBM 기반 프로펠러 및 로터 해석의 검증은 풍동 실측 자료와의 비교, 전통적 CFD 해석과의 비교, 비행 실측 자료와의 비교를 통해 수행된다. 검증 항목은 다음을 포함한다. 첫째, 정지 추력과 효율. 둘째, 후류 속도 분포와 팁 와류 경로. 셋째, 블레이드 표면 압력 분포. 넷째, 원거리 음장의 소음 스펙트럼. LBM 해석 결과는 기존 Navier-Stokes 기반 CFD와 우수한 일치를 보이며, 특히 공력 음향 영역에서 장점이 두드러진다.

10. 로봇공학적 함의

자율 비행 로봇 분야에서 LBM은 다음의 응용에 기여한다. 첫째, 도심항공교통 기체의 공력 음향 설계. 둘째, 복잡한 기체-프로펠러 상호작용의 정량 해석. 셋째, 저소음 드론 개발의 설계 도구. 넷째, 실내 자율 비행 로봇의 환경 상호작용 해석. 다섯째, 고해상도 학습 기반 공력 모형의 훈련 자료 생성. LBM의 발전은 자율 비행 로봇의 공력 설계 정밀도와 설계 효율을 크게 향상시키고 있다.

11. 출처

  • Chen, S., and Doolen, G. D. “Lattice Boltzmann Method for Fluid Flows.” Annual Review of Fluid Mechanics, vol. 30, 1998.
  • Succi, S. The Lattice Boltzmann Equation: For Fluid Dynamics and Beyond. Oxford University Press, 2001.
  • Krüger, T., Kusumaatmaja, H., Kuzmin, A., Shardt, O., Silva, G., and Viggen, E. M. The Lattice Boltzmann Method: Principles and Practice. Springer, 2017.
  • Romani, G., and Casalino, D. “Rotorcraft Blade-Vortex Interaction Noise Prediction Using the Lattice Boltzmann Method.” Journal of Sound and Vibration, vol. 442, 2019.
  • Casalino, D., Diozzi, F., Sannino, R., and Paonessa, A. “Aircraft Noise Reduction Technologies: A Bibliographic Review.” Aerospace Science and Technology, vol. 12, 2008.

12. 버전

v1.0 (2026-04-17)