24.12 유도 유동(Induced Flow)의 분포와 모델링

1. 유도 유동의 정의

유도 유동(induced flow)은 로터가 유체에 공력을 작용시켜 생성한 자유 흐름 이외의 추가적인 유체 운동이다. 로터 디스크 평면에서의 축방향 유도 속도 $v_i$ 와 접선 방향 유도 속도 $v_t$ 가 유도 유동의 주요 성분이다. 유도 유동의 분포는 로터의 공력 하중 분포, 블레이드 기하, 비행 상태에 의해 결정되며, 공력 해석과 성능 예측의 핵심 변수이다.

2. 균일 유도 유동 모형

가장 단순한 유도 유동 모형은 균일 유입(uniform inflow) 가정이다. 이 가정에서 로터 디스크 평면의 모든 지점에서 축방향 유도 속도가 동일하다고 간주한다. 운동량 이론의 결과 $v_i = \sqrt{T / (2 \rho A)}$ 가 적용된다. 이 모형은 초기 설계와 평균 성능 예측에 적합하지만, 실제 유도 속도의 반경·방위각 방향 변화를 재현하지 못한다.

3. 선형 유도 유동 모형

Glauert가 제시한 선형 유도 유동 모형은 전진 비행 상태의 비대칭성을 반영한다. 이 모형은 유도 속도가 반경과 방위각에 대해 다음과 같이 변화한다고 가정한다.

$v_i(r, \psi) = \bar{v}_i (1 + k_x \bar{r} \cos \psi + k_y \bar{r} \sin \psi)$

여기서 $\bar{v}_i$ 는 평균 유도 속도, $\bar{r} = r/R$ , $\psi$ 는 방위각, $k_x$ 와 $k_y$ 는 경사 유입에 의한 경사 계수이다. 이 모형은 전진 비행의 비균일성을 간단히 반영하며, 헬리콥터 해석에서 널리 사용된다.

4. 동적 유도 유동 모형

기체의 비정상 기동이나 로터 속도 변화에 따른 유도 유동의 시간 변화를 재현하기 위해 동적 유도 유동 모형이 사용된다. Pitt-Peters 모형은 3자유도 유도 유동을 1차 미분 방정식으로 표현한다.

$\mathbf{M} \dot{\boldsymbol{\lambda}} + \mathbf{L}^{-1} \boldsymbol{\lambda} = \mathbf{C}_T$

여기서 $\boldsymbol{\lambda} = (\lambda_0, \lambda_s, \lambda_c)^\top$ 는 균일, 종방향, 횡방향 유도 유동 성분, $\mathbf{C}_T$ 는 추력과 모멘트 계수, $\mathbf{M}$ 과 $\mathbf{L}$ 은 질량과 감쇠 행렬이다. 이 모형은 Pitt와 Peters의 Theoretical Prediction of Dynamic-Inflow Derivatives(Vertica, vol. 5, no. 1, 1981)에서 정식화되었다.

5. 블레이드 요소 운동량 이론에서의 유도 유동

블레이드 요소 운동량 이론(Blade Element Momentum Theory)에서 유도 유동은 반경 방향으로 변화하는 유도 인자 $a(r) = v_i / V$ 와 $a'(r) = v_t / (\omega r)$ 로 표현된다. 각 환형 요소의 유도 인자는 운동량 보존과 블레이드 요소 공력 평형의 반복 해법으로 결정된다. 이 접근은 실용적 정확도로 반경 방향 유도 분포를 제공한다.

6. 자유 와류 해석

가장 정확한 유도 유동 해석은 자유 와류 해석(free-wake analysis)을 통해 수행된다. 블레이드에서 방출된 와류 필라멘트가 상호 유도 속도에 의해 스스로 정합적 기하를 형성하도록 반복 계산한다. 각 계산 단계에서 Biot-Savart 법칙을 이용해 공간의 모든 지점에서 유도 속도를 산출한다. 이 접근은 Leishman의 Principles of Helicopter Aerodynamics(2nd ed., Cambridge University Press, 2006)에서 상세히 다루어진다.

7. 멀티로터에서의 유도 유동

멀티로터의 유도 유동은 다수 로터의 유도 속도장이 중첩되어 형성된다. 각 로터의 유도 유동은 다음과 같이 표현된다.

$\mathbf{v}_{\text{induced}}(\mathbf{x}) = \sum_{i=1}^{N} \mathbf{v}_i(\mathbf{x})$

여기서 $\mathbf{v}_i$ 는 $i$ 번째 로터의 유도 속도장이다. 이러한 중첩은 각 로터의 유입 조건을 변화시키고, 성능을 단일 로터의 단순 합과 다르게 만든다.

8. 유도 유동의 주파수 응답

유도 유동의 시간 응답 특성은 주파수 영역에서 1차 필터 특성으로 근사된다.

$\dfrac{v_i(s)}{T(s)} = \dfrac{K}{\tau_i s + 1}$

여기서 $\tau_i$ 는 유도 유동의 시정수, $K$ 는 정적 이득이다. 시정수는 일반적으로 블레이드 통과 주기와 같은 차수이며, 제어기 설계에서 유도 유동의 지연 특성을 반영하는 데 사용된다.

9. 유도 유동 분포의 측정

유도 유동 분포는 다음의 방법으로 실측된다. 첫째, 디스크 평면 근방에서의 피토-정압관 배열 측정. 둘째, PIV를 이용한 유동장 가시화. 셋째, 열선 유속계에 의한 고주파 성분 측정. 넷째, 디스크 하류의 후류 측정으로 간접 추정. 이러한 측정은 이론 해석과 수치 모형의 검증 자료를 제공한다.

10. 로봇공학적 적용

유도 유동의 모델링은 자율 비행 로봇의 다음 측면에 영향을 준다. 첫째, 비행 시뮬레이터의 공력 모형 정확도. 둘째, 급기동 시의 추력 응답 예측. 셋째, 로터 간 간섭에 의한 유도 유동 변화의 재현. 넷째, 제어기의 동적 응답 설계. 다섯째, 지면·벽면 근접 비행의 유도 유동 변화 예측. 이러한 적용에서 유도 유동 모형의 충실도는 자율 비행의 정밀도와 안정성을 결정한다.

11. 학습 기반 보완

최근 연구에서는 학습 기반 접근으로 유도 유동 모형을 보완한다. 첫째, 비행 데이터로부터 유도 유동 잔차를 학습. 둘째, 비선형 유동 변동을 신경망으로 모형화. 셋째, 환경 조건별 유도 유동 특성의 적응적 보정. 이러한 접근은 O’Connell 외의 Neural-Fly Enables Rapid Learning for Agile Flight in Strong Winds(Science Robotics, vol. 7, no. 66, 2022) 등에서 확장된 공력 모형의 예시로 이어진다.

12. 출처

Glauert, H. The Elements of Aerofoil and Airscrew Theory. Cambridge University Press, 1926.
Pitt, D. M., and Peters, D. A. “Theoretical Prediction of Dynamic-Inflow Derivatives.” Vertica, vol. 5, no. 1, 1981.
Leishman, J. G. Principles of Helicopter Aerodynamics, 2nd ed. Cambridge University Press, 2006.
Johnson, W. Helicopter Theory. Princeton University Press, 1980.
O’Connell, M., Shi, G., Shi, X., Azizzadenesheli, K., Anandkumar, A., Yue, Y., and Chung, S.-J. “Neural-Fly Enables Rapid Learning for Agile Flight in Strong Winds.” Science Robotics, vol. 7, no. 66, 2022.

13. 버전

v1.0 (2026-04-17)