23.10 블레이드 요소-운동량 이론(BEMT)의 결합 해석

1. BEMT의 기본 구조

블레이드 요소-운동량 이론(Blade Element Momentum Theory, BEMT)은 반경 방향 환형(annular) 제어 체적에 대한 운동량 보존 해석과 블레이드 요소의 공력 해석을 결합하여, 프로펠러의 유도 속도와 공력 하중을 정량 산출하는 해석 기법이다. 각 환형 구역에서 공력 평형과 운동량 평형이 동시에 성립하도록 유도 인자를 반복적으로 결정한다. 이는 Glauert가 The Elements of Aerofoil and Airscrew Theory(Cambridge University Press, 1926)에서 정식화한 이래, 프로펠러와 풍력 터빈, 헬리콥터 로터 해석의 표준 방법으로 자리 잡았다.

2. 환형 요소의 운동량 평형

반경 $r$ , 두께 $dr$ 의 환형 요소에 대해 통과 유량은 $d\dot{m} = \rho \cdot 2 \pi r \, dr \cdot (V + v_a)$ 이며, 여기서 $V$ 는 자유 흐름 속도, $v_a$ 는 축방향 유도 속도이다. 운동량 보존으로부터 환형 요소에 작용하는 미소 추력은 다음과 같이 표현된다.

$dT_{\text{mom}} = d\dot{m} \cdot 2 v_a = 4 \pi \rho (V + v_a) v_a r \, dr$

또한 회전 성분에 대한 운동량 보존은 다음과 같다.

$dQ_{\text{mom}} = 4 \pi \rho (V + v_a) v_t r^3 \omega \, dr$

여기서 $v_t$ 는 접선 방향 유도 속도이다. 유도 인자 $a = v_a / V$ 와 $a' = v_t / (\omega r)$ 를 도입하면 위 식은 다음과 같이 정리된다.

$dT_{\text{mom}} = 4 \pi \rho V^2 a (1 + a) r \, dr, \quad dQ_{\text{mom}} = 4 \pi \rho V \omega r^3 a' (1 + a) \, dr$

3. 블레이드 요소의 공력 표현

동일 환형 구역에 대해 블레이드 요소 이론으로부터 다음의 미소 추력과 미소 토크가 얻어진다.

$dT_{\text{BE}} = \frac{1}{2} \rho U^2 B c \, (C_l \cos \phi - C_d \sin \phi) \, dr$

$dQ_{\text{BE}} = \frac{1}{2} \rho U^2 B c r \, (C_l \sin \phi + C_d \cos \phi) \, dr$

여기서 $U$ 는 국부 상대 속도, $\phi$ 는 국부 유입각, $B$ 는 블레이드 수이다. 국부 받음각 $\alpha = \beta - \phi$ 에서 $C_l$ 과 $C_d$ 는 블레이드 단면 익형의 공력 극선으로부터 결정된다.

4. 결합 해석과 유도 인자 결정

BEMT의 핵심은 각 환형 구역에서 운동량 이론과 블레이드 요소 이론에서 산출된 미소 추력 및 토크를 서로 일치시키는 것이다. 즉, 다음의 조건이 성립해야 한다.

$dT_{\text{mom}} = dT_{\text{BE}}, \quad dQ_{\text{mom}} = dQ_{\text{BE}}$

이 조건은 유도 인자 $a$ , $a'$ 에 대한 비선형 연립 방정식을 제공하며, 반복 계산으로 해를 구한다. 일반적인 반복 절차는 다음과 같다. 첫째, 초기 추정값 $a_0 = 0$ , $a'_0 = 0$ 을 설정한다. 둘째, 유입각 $\phi$ 와 받음각 $\alpha$ 를 계산한다. 셋째, 블레이드 단면 자료로부터 $C_l$ , $C_d$ 를 결정한다. 넷째, 운동량 방정식과 블레이드 요소 방정식으로부터 새로운 $a$ , $a'$ 를 계산한다. 다섯째, 수렴 조건( $|a_{k+1} - a_k| < \epsilon$ )이 만족될 때까지 반복한다.

5. 팁 손실과 허브 손실의 보정

Prandtl 팁 손실 인자 $F_t$ 와 허브 손실 인자 $F_h$ 는 운동량 이론의 유도 인자에 직접 적용된다. 보정된 운동량 미소 추력과 미소 토크는 다음과 같이 표현된다.

$dT_{\text{mom}}^{F} = F \cdot 4 \pi \rho V^2 a (1 + a) r \, dr, \quad F = F_t F_h$

여기서 $F_t$ 는 블레이드 팁 근방에서 $F_t \to 0$ 으로 수렴하여 유도 인자의 발산을 억제하고, $F_h$ 는 허브 근방에서 유도 인자를 감소시킨다. 이러한 보정은 실측과의 일치도를 크게 향상시키며, 표준 BEMT 절차의 일부로 포함된다.

6. BEMT 결과의 무차원 성능 계수

BEMT의 반경 방향 적분 결과로 얻어지는 총 추력 $T$ 와 총 토크 $Q$ 는 다음과 같이 무차원 성능 계수로 표현된다.

$C_T = \frac{T}{\rho n^2 D^4}, \quad C_Q = \frac{Q}{\rho n^2 D^5}, \quad C_P = \frac{P}{\rho n^3 D^5} = 2 \pi C_Q$

$J = \frac{V}{n D}, \quad \eta = \frac{C_T}{C_P} J$

여기서 $n$ 은 회전 속도(rev/s), $D = 2R$ 는 직경, $P = 2 \pi n Q$ 는 축동력이다. 이 정의는 McCormick이 Aerodynamics, Aeronautics, and Flight Mechanics(2nd ed., Wiley, 1995)에서 표준 형식으로 제시하였다.

7. 해석의 한계와 확장

BEMT는 축대칭 유입과 시간 평균 유동을 전제하므로, 경사 유입(yaw inflow)이나 비정상 유동 조건에서는 정확도가 감소한다. 이러한 경우 Glauert 경사 유입 모형, 동적 유입(dynamic inflow) 모형, Pitt-Peters 모형 등의 확장이 적용된다. 또한 블레이드-블레이드 상호작용과 3차원 유동 효과가 지배적인 영역에서는 자유 와류 해석(free-wake analysis)이나 CFD 해석으로 보완한다. Leishman이 Principles of Helicopter Aerodynamics(2nd ed., Cambridge University Press, 2006)에서 이러한 확장 기법을 정리하였다.

8. 해석 절차의 도해

단계	처리 내용
1	반경을 $N$ 개의 환형 요소로 분할
2	각 요소에서 피치각 $\beta$ , 시위 $c$ 및 익형 극선 입력
3	유도 인자 초기값 $a = a' = 0$ 설정
4	유입각 $\phi$ 및 받음각 $\alpha$ 계산
5	익형 자료로 $C_l$ , $C_d$ 결정
6	운동량-블레이드 요소 평형 조건에서 새 $a$ , $a'$ 계산
7	수렴 시 다음 요소로 이동, 불수렴 시 4로 회귀
8	모든 요소 합으로 $T$ , $Q$ , $P$ , $\eta$ 산출

이 절차는 각종 공학 소프트웨어(QPROP, XROTOR, JavaProp 등) 및 자체 개발 BEMT 코드에서 공통적으로 사용된다. Drela의 QPROP Formulation(MIT Technical Note, 2006)은 BEMT의 수치 절차를 상세히 다룬 대표적 문서이다.

9. 로봇공학적 적용

BEMT는 멀티로터, 소형 고정익 무인기, 도심항공교통 기체의 프로펠러 및 로터 해석에서 광범위하게 사용된다. 특히 다음과 같은 목적에 적합하다. 첫째, 임의 운용 조건에서의 추력·토크·효율 예측. 둘째, 블레이드 기하의 설계 최적화. 셋째, 프로펠러-모터 매칭. 넷째, 시뮬레이션 공력 모형의 테이블 생성. 다섯째, 실시간 추력 추정기에 사용되는 파라미터 식별. 소형 드론의 공력 모델링 정확도 향상을 위해 측정 기반 블레이드 형상 자료와 BEMT를 결합하는 접근이 최근 활발하다.

10. 출처

Glauert, H. The Elements of Aerofoil and Airscrew Theory. Cambridge University Press, 1926.
Leishman, J. G. Principles of Helicopter Aerodynamics, 2nd ed. Cambridge University Press, 2006.
McCormick, B. W. Aerodynamics, Aeronautics, and Flight Mechanics, 2nd ed. Wiley, 1995.
Drela, M. “QPROP Formulation.” MIT Technical Note, 2006.
Brandt, J. B., and Selig, M. S. “Propeller Performance Data at Low Reynolds Numbers.” 49th AIAA Aerospace Sciences Meeting, AIAA Paper 2011-1255, 2011.

11. 버전

v1.0 (2026-04-17)