26.41 지면 효과와 벽면 효과의 전산유체역학(CFD) 해석

1. 전산유체역학 해석의 필요성과 강점

전산유체역학(Computational Fluid Dynamics, CFD) 해석은 지면 효과와 벽면 효과와 같은 복합 공력 현상을 수치적으로 해석하는 강력한 도구이다. 실험적 풍동 시험이 비용과 시간 면에서 제약이 있는 반면, CFD 해석은 다양한 형상과 조건에 대한 광범위한 데이터를 상대적으로 낮은 비용으로 생성할 수 있다. 또한 실험에서 측정하기 어려운 공간 내부의 유동 구조, 압력 분포, 와류 발달 과정 등을 완전히 가시화할 수 있다는 장점이 있다. 본 절에서는 지면 효과 및 벽면 효과의 CFD 해석을 위한 지배 방정식, 난류 모델, 격자 구성, 경계 조건, 그리고 후처리 기법을 체계적으로 서술한다.

2. 지배 방정식

비압축성 점성 유동의 CFD 해석은 Navier-Stokes 방정식에 기반한다.

$\nabla \cdot \mathbf{u} = 0$

$\rho \frac{D\mathbf{u}}{Dt} = -\nabla p + \mu \nabla^2 \mathbf{u} + \mathbf{f}$

여기서 $\mathbf{u}$ 는 속도 벡터, $p$ 는 정압, $\rho$ 는 밀도, $\mu$ 는 동점성 계수, $\mathbf{f}$ 는 체적력이다. 회전 기준계에서 로터 해석을 수행하는 경우 원심력과 코리올리력이 체적력으로 추가된다. 이 방정식은 지면 및 벽면에서의 경계 조건과 함께 풀어져서 유동장이 결정된다.

3. 난류 모델링 기법

실제 회전익기 주변 유동은 고레이놀즈 수 난류이므로 난류 모델링이 필수적이다. CFD 해석에서 사용되는 주요 난류 모델은 다음과 같다. 첫째, Reynolds Averaged Navier-Stokes(RANS) 모델은 난류 변동을 시간 평균하여 풀이하며, 계산 비용이 낮다. 대표적 모델로는 Spalart-Allmaras, $k$ - $\epsilon$ , $k$ - $\omega$ SST가 있다. 둘째, Large Eddy Simulation(LES)는 큰 규모의 난류는 직접 해석하고 작은 규모는 모델링하며, RANS보다 정확하나 계산 비용이 크다. 셋째, Detached Eddy Simulation(DES)은 RANS와 LES의 혼합으로 벽면 근처는 RANS로, 자유 영역은 LES로 처리한다. 지면 근접 유동의 경우 후류와 경계층의 복합 현상이 포함되므로 DES 또는 LES가 선호된다.

4. 격자 구성

격자(mesh) 구성은 CFD 해석의 정확도에 결정적 영향을 미친다. 로터 해석에서는 회전 영역과 정지 영역을 구분하여 서로 다른 격자 시스템을 적용한다. 중첩 격자(overset mesh) 또는 미끄럼 격자(sliding mesh) 기법이 이 경우 사용된다. 지면 및 벽면 경계층을 정확히 해석하기 위해 무차원 벽거리 $y^+$ 를 1 이하로 유지하는 세밀한 경계층 격자가 요구된다. 전체 격자 수는 해석의 충실도에 따라 수백만에서 수십억 셀에 이를 수 있다. 격자 수렴성(mesh convergence) 검토는 해석 결과의 신뢰성 확보에 필수적이다.

5. 로터 모델링 접근

로터를 CFD에서 모델링하는 방법은 세 가지로 구분된다. 첫째, 완전 격자 해석(fully resolved analysis)은 블레이드 형상을 명시적으로 격자화하고 회전 운동을 모사한다. 가장 정확하나 계산 비용이 매우 크다. 둘째, 액추에이터 디스크 모델(actuator disk model)은 로터를 체적력 분포로 근사하여 계산 비용을 크게 감소시킨다. 셋째, 액추에이터 라인 모델(actuator line model)은 액추에이터 디스크의 중간 수준 정확도와 비용을 제공한다. 지면 효과와 벽면 효과 연구에서는 목적에 따라 적절한 수준이 선택된다.

6. 경계 조건의 설정

CFD 해석의 경계 조건은 물리적 상황을 수치적으로 표현한다. 주요 경계 조건은 다음과 같다. 첫째, 지면 및 벽면에서는 점착 조건(no-slip condition)과 비침투 조건(no-penetration condition)이 부과된다. 둘째, 이동 지면의 경우 벽면 접선 속도를 자유류 속도와 동일하게 설정한다. 셋째, 외부 경계에서는 자유류 조건(far-field condition) 또는 압력 출구 조건(pressure outlet)이 사용된다. 넷째, 회전 로터의 경우 블레이드 표면에서 점착 조건과 함께 회전 기준계 조건이 적용된다. 경계 조건의 오설정은 해석 결과를 크게 왜곡시킬 수 있다.

7. 시간 적분과 수렴성

정상 상태 해석에서는 시간 단계별 반복 계산을 통해 수렴된 해를 구한다. 비정상 해석에서는 시간 단계 크기 $\Delta t$ 가 중요하며, Courant-Friedrichs-Lewy(CFL) 조건을 만족해야 한다. 일반적으로 로터 해석에서는 시간 단계를 블레이드 1회전당 1000 단계 이상으로 설정한다. 수렴성은 잔차(residual)의 감소와 관심 변수(예: 추력 계수)의 수렴으로 판단된다.

8. 검증과 타당성 확인

CFD 해석 결과는 실험 데이터와의 비교를 통해 검증된다. 일반적인 검증 대상은 추력 계수, 토크 계수, 그리고 특정 위치에서의 속도 및 압력이다. 검증 기준은 응용 목적에 따라 다르며, 공학적 설계에서는 5% 이내의 오차, 학술 연구에서는 1% 이내의 오차가 요구되는 경우가 있다. 검증 통과한 해석 결과는 실험으로 접근하기 어려운 조건에 대한 예측에 사용된다.

9. 지면 효과 CFD 연구 사례

지면 효과의 CFD 연구 사례로는 다음이 있다. Griffiths, Ananthan, and Leishman의 “Predictions of Rotor Performance in Ground Effect Using a Free-Vortex Wake Model”(Journal of the American Helicopter Society, 2005)은 자유 와류 방법을 사용한 지면 효과 해석을 제시한다. Nathan과 Green의 “Measurements of a Rotor Flow in Ground Effect and Visualisation of the Brownout Phenomenon”(American Helicopter Society 64th Annual Forum, 2008)은 지면 효과와 먼지 비산 현상의 해석을 다룬다. 이러한 연구들은 CFD 해석의 다양한 적용 사례를 제공한다.

10. 벽면 효과 CFD 연구 사례

벽면 효과의 CFD 연구는 최근 멀티로터 드론의 실내 및 구조물 근접 비행 증가와 함께 확대되고 있다. Kan, Teo, Lai, and Ng의 “Computational Investigation of Rotor Aerodynamics in Proximity to Vertical Walls”(AIAA SciTech Forum, 2018) 및 관련 연구들은 다양한 벽면 거리에서의 공력 특성을 수치적으로 해석하였다. 이러한 연구들은 벽면 효과에 대한 해석 모델의 검증과 개선에 기여한다.

11. 복합 환경의 CFD 해석

실내 비행, 협소 공간, 도심 캐니언 등 복합 환경의 CFD 해석은 대규모 계산 자원을 요구한다. 일반적으로 수천만에서 수억 격자 셀이 필요하며, 고성능 컴퓨팅(HPC) 자원이 활용된다. 해석 비용을 감소시키기 위한 접근으로는 다단계 격자(multi-level mesh), 적응 격자 세분화(adaptive mesh refinement), 그리고 병렬 계산이 있다. 또한 저차 대리 모델(reduced-order surrogate model)을 고차 CFD 결과로부터 학습하여 실시간 예측에 사용하는 접근도 활발히 연구되고 있다.

12. 상용 및 오픈소스 CFD 소프트웨어

지면 효과 및 벽면 효과 해석에 사용되는 주요 CFD 소프트웨어는 다음과 같다. 상용 소프트웨어로는 ANSYS Fluent, STAR-CCM+, Siemens Simcenter가 있으며, 오픈소스로는 OpenFOAM, SU2, CFDSim이 있다. 각 소프트웨어는 강점과 약점이 있으므로 응용 목적에 따라 선택된다. 로터 해석에 특화된 소프트웨어로는 Helios(미국 육군 항공성 개발), ROTCFD 등이 있다.

13. 실용적 설계 지침

CFD 해석을 성공적으로 수행하기 위한 실용적 지침은 다음과 같다. 첫째, 단순 사례에 대한 검증을 먼저 수행하여 설정의 타당성을 확인한다. 둘째, 격자 수렴성 연구를 통해 격자 크기가 결과에 미치는 영향을 정량화한다. 셋째, 적절한 난류 모델을 선택하고 그 한계를 인지한다. 넷째, 경계 조건을 물리적으로 타당하게 설정한다. 다섯째, 시간 단계와 반복 횟수가 충분한지 수렴성을 확인한다. 여섯째, 결과를 실험 데이터 또는 이론 예측과 비교하여 타당성을 검증한다.

14. 미래 연구 방향

지면 효과 및 벽면 효과의 CFD 해석은 다음 방향으로 발전하고 있다. 첫째, 고정확도 LES 및 DNS(Direct Numerical Simulation) 해석의 확대이다. 둘째, 기계 학습을 활용한 CFD 가속화와 대리 모델 개발이다. 셋째, 멀티 피지크스 해석(공력-구조-열 결합)의 확대이다. 넷째, 실시간 CFD 기반 제어(실시간 보상) 연구이다. 다섯째, 클라우드 컴퓨팅 기반 CFD 서비스의 확산이다. 이러한 발전은 지면 효과와 벽면 효과의 더욱 정확한 이해와 드론 설계 및 제어의 향상을 지원한다.

15. 출처

Griffiths, D. A., Ananthan, S., and Leishman, J. G., “Predictions of Rotor Performance in Ground Effect Using a Free-Vortex Wake Model,” Journal of the American Helicopter Society, Vol. 50, No. 4, 2005.
Nathan, N. D., and Green, R. B., “Measurements of a Rotor Flow in Ground Effect and Visualisation of the Brownout Phenomenon,” American Helicopter Society 64th Annual Forum, 2008.
Kan, X., Teo, Y. H., Lai, A. C. H., and Ng, B. F., “Computational Investigation of Rotor Aerodynamics in Proximity to Vertical Walls,” AIAA SciTech Forum, 2018.
Leishman, J. G., “Principles of Helicopter Aerodynamics,” 2nd ed., Cambridge University Press, 2006.
Ferziger, J. H., Perić, M., and Street, R. L., “Computational Methods for Fluid Dynamics,” 4th ed., Springer, 2020.
Spalart, P. R., “Detached-Eddy Simulation,” Annual Review of Fluid Mechanics, Vol. 41, 2009.

16. 버전

문서 버전: v1.0
작성 기준일: 2026-04-17