23.46 프로펠러의 전산유체역학(CFD) 해석
1. 프로펠러 CFD 해석의 의의
전산유체역학(Computational Fluid Dynamics, CFD) 해석은 프로펠러 주변 유동의 Navier-Stokes 방정식을 수치적으로 해결하여, 공력 성능, 후류 구조, 소음 원천, 비정상 공력 응답 등 고차원 정보를 제공하는 해석 기법이다. 블레이드 요소 운동량 이론(BEMT)과 같은 저차원 해석이 공급하지 못하는 세부 유동 구조를 재현할 수 있어, 설계의 최종 검증과 복잡한 조건의 해석에 사용된다. CFD 해석의 표준 절차는 Ferziger와 Perić의 Computational Methods for Fluid Dynamics(3rd ed., Springer, 2002)에서 체계적으로 기술되어 있다.
2. 격자 생성
프로펠러 CFD 해석의 격자 생성은 다음의 요소를 고려한다. 첫째, 블레이드 표면의 경계층을 해상하기 위해 블레이드에 수직한 방향으로 세밀한 격자가 배치되며, 첫 격자점의 무차원 거리 y^+는 1 이하로 유지한다. 둘째, 블레이드 전연부와 팁 주변은 곡률이 크고 박리 가능성이 있으므로 격자를 밀집시킨다. 셋째, 블레이드 후류가 전파되는 영역은 중간 해상도로 격자화한다. 넷째, 블레이드 회전을 처리하기 위해 회전 영역과 정지 영역을 분리한 격자 구성이 사용된다.
3. 회전 처리 기법
| 기법 | 설명 |
|---|---|
| 회전 참조계 (MRF) | 정지 해석에서 회전계 효과를 포함하는 평균화된 접근 |
| 슬라이딩 격자 (Sliding Mesh) | 회전 영역과 정지 영역의 경계에서 격자가 미끄러지는 방식 |
| 오버셋 격자 (Overset, Chimera) | 독립적으로 구축된 격자가 서로 중첩되어 보간으로 결합 |
| 적응 격자 세분화 (AMR) | 해상이 필요한 영역에 격자를 자동 세분화 |
| 액추에이터 디스크 모형 | 블레이드를 분포 체적력으로 대체하여 격자 간소화 |
이 표는 프로펠러 CFD 해석의 회전 처리 기법을 요약한 것이다. 정상 상태 해석에는 MRF가, 비정상 해석에는 슬라이딩 격자 또는 오버셋 격자가 주로 사용된다.
4. 난류 모형
프로펠러의 CFD 해석에서 사용되는 난류 모형은 다음과 같다. 첫째, Spalart-Allmaras 1방정식 모형은 부착 유동에서 계산 효율이 높고 안정적이다. 둘째, k-\omega SST 모형은 벽면과 자유 전단류 모두에 우수하며, 박리를 포함한 일반적 해석에 널리 사용된다. 셋째, 전이 모형(Langtry-Menter \gamma-\tilde{Re}_{\theta t} 모형)은 층류-난류 천이를 예측하여 저 Re 프로펠러 해석에 필수적이다. 넷째, Detached Eddy Simulation(DES)과 Large Eddy Simulation(LES)은 대규모 박리와 BVI를 포함하는 비정상 해석에 사용된다.
5. 경계 조건 설정
프로펠러 CFD 해석의 경계 조건은 다음과 같이 설정된다. 첫째, 자유 흐름 경계는 멀리 떨어진 영역에 배치되며, 전진 속도와 정적 압력을 지정한다. 둘째, 블레이드 표면은 비점착 조건(no-slip)이다. 셋째, 회전 축 주위에 축대칭 경계 조건을 적용하여 계산 영역을 감소시킬 수 있다. 넷째, 블레이드 간 주기 경계 조건(periodic boundary)은 단일 블레이드 영역만으로 전체 해석을 수행할 수 있게 한다.
6. 해석 절차
프로펠러 CFD 해석의 일반적 절차는 다음과 같다. 첫째, 블레이드의 3차원 CAD 모형을 구축하고 격자를 생성한다. 둘째, 유체 영역을 정의하고 적절한 회전 처리 기법을 선택한다. 셋째, 초기 조건으로 자유 흐름 상태를 설정한다. 넷째, 정상 상태 해석을 먼저 수행하여 수렴된 해를 획득한다. 다섯째, 필요에 따라 비정상 해석으로 진행한다. 여섯째, 블레이드 표면 압력과 점성 응력을 적분하여 추력과 토크를 계산한다. 일곱째, 후류 구조, 유도 속도장, 와류 분포를 분석한다.
7. 수렴성과 정확도 검증
CFD 해석의 수렴성은 잔차(residual), 질량 보존 오차, 공력 계수의 반복 변화를 관찰하여 평가한다. 정확도 검증은 격자 의존성 시험(grid convergence study)으로 수행된다. Richardson 외삽을 이용한 GCI(Grid Convergence Index) 방법은 격자 수렴성을 정량화하는 표준 절차이다. Roache가 Verification and Validation in Computational Science and Engineering(Hermosa Publishers, 1998)에서 이 방법을 정식화하였다.
8. 실측 자료와의 검증
CFD 해석 결과는 풍동 시험과 비행 시험의 실측 자료와 비교하여 검증된다. 주요 검증 대상은 다음과 같다. 첫째, 추력과 동력 계수의 J 의존성. 둘째, 후류 내 속도 분포와 와류 위치. 셋째, 블레이드 표면 압력 분포. 넷째, 소음 스펙트럼. NASA와 AIAA의 국제 공력 해석 워크숍(Drag Prediction Workshop, Rotorcraft Hover Prediction Workshop 등)은 CFD 해석의 체계적 검증을 수행하는 대표적 협력 활동이다.
9. 대표 CFD 코드
프로펠러 CFD 해석에 사용되는 대표적 상용 및 공개 코드는 다음과 같다.
| 코드 | 특성 |
|---|---|
| NASA OVERFLOW | 오버셋 격자 기반 회전익 전용 |
| ANSYS Fluent | 상용 범용 CFD 코드 |
| Star-CCM+ | 상용 범용 CFD 코드 |
| OpenFOAM | 공개 소스 C++ 기반 |
| SU2 | Stanford 개발 공개 소스 |
| CFX | ANSYS 회전 기계 전용 |
이 표는 프로펠러 해석에 사용되는 대표 CFD 코드를 예시한 것이다. 각 코드는 특정 해석 유형에 강점을 가지며, 적용 목적에 따라 선택된다.
10. 소음 예측과의 결합
프로펠러 CFD 해석은 공력 음향 예측과 결합되어 소음 분석에 활용된다. 일반적 접근은 다음과 같다. 첫째, CFD에서 비정상 블레이드 표면 압력을 계산한다. 둘째, 계산된 압력을 Ffowcs Williams-Hawkings(FW-H) 방정식에 입력하여 원격 음장의 음압을 산출한다. 셋째, 광대역 소음은 경계층 난류 통계를 바탕으로 별도 모형으로 예측한다. 이러한 결합은 CFD-CAA(Computational AeroAcoustics) 접근이라 한다.
11. 로봇공학적 적용
자율 비행 로봇 및 도심항공교통 기체의 설계에서 CFD 해석은 다음과 같이 활용된다. 첫째, 복잡한 기체-프로펠러 상호작용의 정량 해석. 둘째, 근접 비행 조건의 로터 간섭 평가. 셋째, 지면 근접 재순환 유동의 재현. 넷째, 저소음 설계의 물리적 원리 이해. 다섯째, 학습 기반 공력 모형의 훈련 자료 생성. 이러한 해석은 단순 이론으로 포착 불가능한 세부 유동 구조를 제공하여, 자율 비행 시스템의 정밀 설계에 기여한다.
12. 출처
- Ferziger, J. H., and Perić, M. Computational Methods for Fluid Dynamics, 3rd ed. Springer, 2002.
- Roache, P. J. Verification and Validation in Computational Science and Engineering. Hermosa Publishers, 1998.
- Spalart, P. R., and Allmaras, S. R. “A One-Equation Turbulence Model for Aerodynamic Flows.” AIAA Paper 92-0439, 1992.
- Menter, F. R. “Two-Equation Eddy-Viscosity Turbulence Models for Engineering Applications.” AIAA Journal, vol. 32, no. 8, 1994.
- Ventura Diaz, P., and Yoon, S. “High-Fidelity Computational Aerodynamics of Multi-Rotor Unmanned Aerial Vehicles.” 56th AIAA Aerospace Sciences Meeting, AIAA Paper 2018-1266, 2018.
13. 버전
v1.0 (2026-04-17)