21.32 전산 유체역학(CFD)의 공기역학 응용 (Computational Fluid Dynamics Applications in Aerodynamics)

21.32 전산 유체역학(CFD)의 공기역학 응용 (Computational Fluid Dynamics Applications in Aerodynamics)

1. 전산 유체역학의 개요

전산 유체역학(Computational Fluid Dynamics, CFD)은 유체의 운동을 지배하는 편미분 방정식(나비에-스토크스 방정식 등)을 컴퓨터를 이용하여 수치적으로 풀어 유동장의 속도, 압력, 온도 등의 분포를 구하는 학문 분야이다. CFD는 풍동 실험과 함께 현대 공기역학 해석의 양대 축을 구성하며, 풍동 실험이 어렵거나 비용이 과도한 경우에 강력한 대안을 제공한다.

CFD의 기본 절차는 다음과 같다:

  1. 전처리(preprocessing): 해석 대상의 기하학적 형상을 정의하고, 해석 영역을 이산화(discretization)하여 격자(mesh/grid)를 생성한다.
  2. 풀이(solving): 지배 방정식을 격자 위에서 수치적으로 이산화하고 반복 계산(iteration)하여 해를 구한다.
  3. 후처리(postprocessing): 계산 결과를 시각화하고, 공력 계수, 압력 분포, 유동 구조 등을 분석한다.

2. 지배 방정식

공기역학적 CFD 해석의 지배 방정식은 나비에-스토크스 방정식(Navier-Stokes equations)으로서, 질량 보존(연속 방정식), 운동량 보존, 에너지 보존의 법칙을 수학적으로 표현한 것이다. 비압축성 유동에서의 연속 방정식과 운동량 방정식은 다음과 같다:

\nabla \cdot \mathbf{V} = 0

\rho\left(\frac{\partial \mathbf{V}}{\partial t} + \mathbf{V} \cdot \nabla\mathbf{V}\right) = -\nabla p + \mu \nabla^2 \mathbf{V}

압축성 유동에서는 밀도의 시간 변화와 에너지 방정식이 추가적으로 고려된다. 이 방정식 체계는 비선형이고 결합되어 있어 일반적으로 해석적 풀이가 불가능하며, 수치적 방법에 의존한다.

3. 격자 생성

격자 생성(mesh generation)은 CFD 해석의 정확도와 효율에 결정적 영향을 미치는 핵심 전처리 단계이다. 주요 격자 유형은 다음과 같다:

구조 격자(structured grid): 격자의 위상(topology)이 규칙적인 배열을 따른다. 계산 효율이 높고 경계층 해상도의 정밀 제어가 용이하지만, 복잡한 형상에 대한 격자 생성이 어렵다. O형, C형, H형 격자가 익형 해석에 전통적으로 사용된다.

비구조 격자(unstructured grid): 삼각형(2D) 또는 사면체(3D) 요소의 비규칙적 배열로 구성된다. 복잡한 형상에 대한 격자 생성이 자동화되어 용이하지만, 계산 효율이 구조 격자보다 낮고 경계층 해상도의 제어가 상대적으로 어렵다.

혼합 격자(hybrid grid): 경계층 영역에는 프리즘(prism) 또는 육면체(hexahedral) 구조 격자를 적용하고, 나머지 영역에는 비구조 격자를 적용하는 방식이다. 경계층 해상도와 격자 생성 용이성을 동시에 확보할 수 있어 실용 해석에서 널리 사용된다.

경계층의 정확한 해석을 위해 물체 표면 부근의 격자는 충분히 조밀하여야 하며, 첫 번째 격자점의 벽면 법선 거리를 무차원 벽면 거리 y^+로 표현하여 관리한다. RANS 해석에서 벽면 함수를 사용하지 않는 경우 y^+ < 1이 요구된다.

4. 난류 모델링

대부분의 공기역학적 유동은 난류 상태이며, 난류의 모든 스케일을 직접 해석하는 것은 극도로 높은 계산 비용을 요구한다. 따라서 실용적 CFD 해석에서는 난류 모델(turbulence model)을 사용하여 난류의 효과를 근사적으로 모델링한다.

레이놀즈 평균 나비에-스토크스(RANS) 방법: 속도와 압력을 시간 평균값과 변동 성분으로 분해하고, 변동 성분에 의한 레이놀즈 응력(Reynolds stress)을 난류 모델로 닫는(close) 방법이다. 대표적 RANS 난류 모델은 다음과 같다:

  • Spalart-Allmaras(SA) 모델: 1방정식 모델로서 외부 공기역학 해석에 널리 사용된다. 부착 유동에서 우수한 예측 정확도를 제공한다.
  • k-\omega SST 모델: Menter(1994)가 개발한 2방정식 모델로서, 벽면 부근에서 k-\omega 모델의 장점과 외부 유동에서 k-\epsilon 모델의 장점을 결합한다. 역압력 구배와 박리 유동에서 SA 모델보다 우수한 예측을 제공하며, 공기역학 분야에서 표준적으로 사용된다.

대와류 시뮬레이션(Large Eddy Simulation, LES): 난류의 대규모 와류 구조를 직접 해석하고, 격자 해상도보다 작은 소규모 와류만을 모델링하는 방법이다. RANS보다 정확한 비정상 유동 해석이 가능하지만, 계산 비용이 RANS의 수십~수백 배에 달한다.

직접 수치 시뮬레이션(Direct Numerical Simulation, DNS): 난류의 모든 스케일을 격자로 직접 해상하는 방법으로, 모델링 오차가 없지만 계산 비용이 극도로 높다(Re^3에 비례). 저레이놀즈 수 유동의 벤치마크 연구에 제한적으로 사용된다.

5. 천이 모델링

저레이놀즈 수 비행 로봇의 공력 해석에서 경계층 천이(laminar-turbulent transition)의 정확한 예측은 특히 중요하다. 표준 RANS 난류 모델은 경계층 천이를 예측하지 못하므로, 천이 모델을 추가적으로 결합하여야 한다.

\gamma-Re_\theta 천이 모델(Langtry & Menter, 2009): 간헐도(intermittency) \gamma와 천이 운동량 두께 레이놀즈 수 Re_\theta에 대한 두 개의 수송 방정식을 추가하여 자연 천이, 바이패스 천이, 박리 유도 천이를 예측한다. 상용 CFD 코드에 널리 구현되어 있다.

e^N 방법과의 결합: RANS 해석에서 경계층 안정성 해석을 후처리적으로 수행하여 천이 위치를 예측하고, 이를 다시 RANS 해석에 반영하는 반복적 방법이다.

6. 수치 방법과 이산화

지배 방정식을 격자 위에서 이산화하는 주요 수치 방법은 다음과 같다:

유한 체적법(Finite Volume Method, FVM): 해석 영역을 제어 체적(control volume)으로 분할하고, 각 제어 체적에 대해 보존 법칙을 적분 형태로 적용한다. 보존 성질(conservation property)이 본질적으로 보장되어 공기역학 CFD에서 가장 널리 사용된다.

유한 요소법(Finite Element Method, FEM): 해석 영역을 요소(element)로 분할하고, 시험 함수(test function)와 가중 잔차(weighted residual) 원리를 이용하여 약형식(weak form)으로 이산화한다. 복잡한 기하학과 비구조 격자에 유연하다.

시간 이산화에서 정상 유동 해석에는 의사 시간 진행(pseudo time-stepping)이, 비정상 유동 해석에는 이중 시간 진행(dual time-stepping) 또는 직접 시간 적분이 사용된다.

7. 익형과 날개의 CFD 해석

2차원 익형의 CFD 해석은 격자 의존성 연구(grid convergence study)와 난류/천이 모델의 선택이 결과의 정확도를 좌우한다. 일반적인 2차원 익형 RANS 해석의 절차는 다음과 같다:

  1. C형 또는 O형 구조 격자를 생성하며, 원방 경계(far-field boundary)는 시위 길이의 20~50배 이상으로 설정한다.
  2. 벽면 법선 방향 첫 번째 격자 간격을 y^+ < 1로 설정한다.
  3. 적절한 난류/천이 모델을 선택한다.
  4. 격자 수렴(mesh convergence) 연구를 수행하여 격자 독립적 결과를 확보한다.

3차원 날개 및 전기체(complete vehicle)의 CFD 해석에서는 격자 수가 수백만~수천만 셀에 달하며, 병렬 컴퓨팅(parallel computing)이 필수적이다. 오버셋 격자(overset/chimera grid) 기법은 복잡한 다물체 구성(로터+기체, 날개+동체 등)의 해석에 유용하다.

8. 프로펠러와 로터의 CFD 해석

프로펠러와 로터의 CFD 해석에서는 블레이드의 회전을 처리하기 위한 특수한 기법이 필요하다:

다중 기준 좌표계(Multiple Reference Frame, MRF): 정상 유동 해석에서 회전 영역에 회전 기준 좌표계를 적용하여 블레이드가 정지한 것처럼 해석한다. 계산 비용이 낮지만, 블레이드-와류 상호 작용과 같은 비정상 효과를 포착할 수 없다.

슬라이딩 메시(sliding mesh): 회전 영역과 고정 영역 사이의 계면에서 격자가 시간에 따라 미끄러지며, 비정상 유동을 직접 해석한다. 로터 간 간섭, 후류 구조의 시간적 발전 등을 포착할 수 있지만, 계산 비용이 높다.

액추에이터 디스크/액추에이터 라인 모델: 블레이드를 직접 격자로 해석하지 않고, 체적력(body force)으로 모델링하는 저충실도 방법이다. 다중 로터 시스템의 간섭 해석이나 전기체 해석에서 계산 효율을 크게 향상시킨다.

9. CFD 해석의 검증과 확인

CFD 결과의 신뢰성을 확보하기 위한 검증(verification)과 확인(validation) 절차는 필수적이다:

검증(verification): 수치적 방법이 지배 방정식을 올바르게 풀고 있는지 확인하는 과정이다. 격자 수렴 연구, 반복 수렴(iteration convergence) 확인, 보존량의 잔차(residual) 검사가 포함된다. Richardson 외삽법(Richardson extrapolation)에 의한 격자 수렴 지수(Grid Convergence Index, GCI)의 산출이 권장된다.

확인(validation): 수치적 결과가 물리적 현실을 정확히 반영하는지 실험 데이터와 비교하여 확인하는 과정이다. 압력 분포, 양력·항력 계수, 유동 구조 등을 풍동 실험 데이터와 체계적으로 비교한다.

10. 로봇 공학에서의 CFD 활용

비행 로봇의 설계와 해석에서 CFD는 다음과 같이 활용된다:

  • 익형 및 블레이드 설계: 다양한 익형 후보에 대해 저레이놀즈 수 공력 특성을 예측하고 최적 익형을 선정한다.
  • 다중 로터 간섭 해석: 멀티로터 시스템에서 로터 간의 공기역학적 간섭을 정량적으로 평가한다.
  • 지면/벽면 효과 해석: 구조물 근접 비행에서의 공력 변화를 예측한다.
  • 전기체 공력 해석: 날개, 동체, 프로펠러, 미익 등 전체 구성 요소를 포함한 통합 공력 해석을 수행한다.
  • 비행 시뮬레이션 모델 생성: CFD에서 산출된 공력 데이터베이스를 비행 역학 시뮬레이션에 공급한다.

오픈소스 CFD 코드인 OpenFOAM과 상용 코드인 ANSYS Fluent, STAR-CCM+ 등이 비행 로봇의 공기역학 해석에 널리 활용되고 있다 (Anderson, 2017).

참고 문헌

  • Anderson, J. D. (2017). Fundamentals of Aerodynamics (6th ed.). McGraw-Hill Education.
  • Langtry, R. B., & Menter, F. R. (2009). Correlation-based transition modeling for unstructured parallelized computational fluid dynamics codes. AIAA Journal, 47(12), 2894–2906.
  • Menter, F. R. (1994). Two-equation eddy-viscosity turbulence models for engineering applications. AIAA Journal, 32(8), 1598–1605.
  • Versteeg, H. K., & Malalasekera, W. (2007). An Introduction to Computational Fluid Dynamics: The Finite Volume Method (2nd ed.). Pearson Education.

버전: 2026-04-16-v1.0