22.41 양력과 항력의 전산유체역학(CFD) 해석

1. CFD 해석의 수학적 기반

전산유체역학(Computational Fluid Dynamics, CFD)은 Navier-Stokes 방정식을 이산화하여 수치적으로 해석하는 기법이다. 압축성·점성 유동의 지배 방정식은 다음과 같이 보존 형식으로 표현한다.

$\frac{\partial \mathbf{U}}{\partial t} + \nabla \cdot \mathbf{F}(\mathbf{U}) - \nabla \cdot \mathbf{F}_v(\mathbf{U}, \nabla \mathbf{U}) = 0$

여기서 $\mathbf{U} = (\rho, \rho \mathbf{u}, \rho E)^\top$ 는 보존 변수, $\mathbf{F}$ 는 비점성 플럭스, $\mathbf{F}_v$ 는 점성·열전도 플럭스이다. 비압축성 흐름에서는 연속 방정식 $\nabla \cdot \mathbf{u} = 0$ 과 운동량 방정식의 결합으로 표현되며, 압력-속도 결합 알고리즘(SIMPLE, PISO, PIMPLE)이 적용된다. 이와 같은 지배 방정식의 수치 해석은 Ferziger와 Perić이 Computational Methods for Fluid Dynamics(3rd ed., Springer, 2002)에서 체계적으로 정리한 바 있다.

2. 공간 이산화 기법

CFD 코드에서는 유한 차분법(finite difference method), 유한 체적법(finite volume method), 유한 요소법(finite element method), 스펙트럴 방법(spectral method)이 사용된다. 실용 공력 해석에서는 보존성 요구 때문에 유한 체적법이 지배적이다. 대류 항의 이산화에는 풍상 스킴(upwind scheme), MUSCL(Monotonic Upstream-centered Scheme for Conservation Laws), WENO(Weighted Essentially Non-Oscillatory) 등이 선택되며, 점성 항은 중심 차분으로 이산화된다. 시간 적분에는 암시적(Implicit) 방법과 명시적(Explicit) Runge-Kutta 방법이 모두 활용되며, 정상 해석에는 국소 시간 스텝 기법(local time stepping)과 다중 격자법(multigrid method)이 수렴 가속을 위해 사용된다.

3. 격자 생성과 경계층 해상

격자(mesh, grid)는 유한 체적을 구성하는 이산 단위로, 구조 격자(structured), 비구조 격자(unstructured), 혼합(hybrid), 중첩 격자(overset), 카테시안 잘라내기 격자(Cartesian cut-cell) 등이 있다. 양력과 항력 해석에서는 경계층 해상이 결정적 요인이다. 경계층 첫 격자점의 무차원 거리 $y^+$ 는 난류 모형에 따라 권장 범위가 다르며, 벽 해결형(wall-resolved) 접근에서는 $y^+ \approx 1$ 을 유지한다. 벽 함수(wall function) 접근에서는 $30 \lesssim y^+ \lesssim 300$ 범위의 격자 구성이 사용된다. 이러한 격자 전략은 NASA Turbulence Modeling Resource의 권장 기준 및 다수의 AIAA Drag Prediction Workshop 보고서에서 공통적으로 적용된다.

4. 난류 모형의 선택

난류 유동은 Navier-Stokes 방정식의 Reynolds 평균(Reynolds-Averaged Navier-Stokes, RANS) 또는 Large Eddy Simulation(LES) 등으로 처리한다. 양력과 항력 해석에 자주 사용되는 난류 모형은 다음과 같다.

난류 모형	특징	대표 적용
Spalart-Allmaras	단일 이송 방정식, 벽면 유동에 우수	순항 조건 고속 항공기 항력 예측
$k$ - $\omega$ SST	벽면·자유 전단류 양자에 균형	박리·실속 영역 일반 해석
$k$ - $\varepsilon$	고 $Re$ 자유 전단류에 적합	실내 환기, 환경 흐름
Reynolds Stress Model	이방성 난류 변동 재현	강한 선회·3D 박리 유동
Detached Eddy Simulation	RANS + LES 혼합	대규모 박리를 포함하는 비정상 흐름
Large Eddy Simulation	대와류 직접 해석	저 $Re$ 박리, BVI, 와류 발산

Spalart와 Allmaras가 A One-Equation Turbulence Model for Aerodynamic Flows(AIAA Paper 92-0439, 1992), Menter가 Two-Equation Eddy-Viscosity Turbulence Models for Engineering Applications(AIAA Journal, vol. 32, no. 8, 1994)에서 각 모형의 정식화를 제시한 이래, 다수의 공력 해석에서 표준 모형으로 자리 잡았다.

5. 층류-난류 천이 모형

저 $Re$ 및 고고도 비행체의 양력과 항력 해석에서는 층류-난류 천이 예측이 필수이다. Langtry와 Menter가 A Correlation-Based Transition Model Using Local Variables—Part I: Model Formulation(Journal of Turbomachinery, vol. 128, no. 3, 2006)에서 제시한 $\gamma$ - $\tilde{Re}_{\theta t}$ 모형은 지역 변수를 이용하는 전이 모형으로 널리 사용된다. 또한 $e^N$ 방법은 경계층 선형 안정성 해석을 통해 전이 위치를 예측하는 고전적 기법으로, Drela의 XFOIL이 이를 2차원 익형 해석에 적용한 대표 예이다.

6. 양력 및 항력 계수의 추출

CFD 해석 결과로부터 양력 $L$ 과 항력 $D$ 는 모형 표면의 압력과 점성 응력 적분으로 산출한다. 표면적 $S$ 상의 외향 단위 법선 벡터 $\hat{\mathbf{n}}$ 과 표면 응력 $\boldsymbol{\tau}$ 에 대해 공력 $\mathbf{F}$ 는 다음과 같다.

$\mathbf{F} = \oint_S \left( -p \hat{\mathbf{n}} + \boldsymbol{\tau} \cdot \hat{\mathbf{n}} \right) dS$

자유 흐름 방향 단위 벡터 $\hat{\mathbf{e}}_\infty$ 와 양력 방향 단위 벡터 $\hat{\mathbf{e}}_L$ 을 이용하면 $D = \mathbf{F} \cdot \hat{\mathbf{e}}_\infty$ , $L = \mathbf{F} \cdot \hat{\mathbf{e}}_L$ 이며, 무차원화하여 $C_L$ , $C_D$ 계수가 얻어진다. 항력 분해가 필요한 경우 파동 항력, 유도 항력, 점성 항력 성분을 원방 운동량 보존을 이용하여 분리하는 far-field 해석이 사용된다.

7. 수렴성, 검증, 확인

CFD 자료의 신뢰성을 확보하기 위해 다음 세 단계를 적용한다. 첫째, 수렴성(convergence) 확인은 잔차(residual)와 공력 계수의 반복 수렴을 모두 점검한다. 둘째, 검증(verification)은 수치 해가 지배 방정식의 정확한 해에 접근함을 확인하는 절차이며, 격자 수렴성 연구(Grid Convergence Study)와 Richardson 외삽이 사용된다. Roache가 Verification and Validation in Computational Science and Engineering(Hermosa Publishers, 1998)에서 GCI(Grid Convergence Index) 절차를 정식화하였다. 셋째, 확인(validation)은 수치 해가 실험 또는 해석 해와 물리적으로 일치함을 검토하며, 동일 조건의 풍동 시험 자료나 벤치마크 해를 사용한다. AIAA G-077-1998 Guide for the Verification and Validation of Computational Fluid Dynamics Simulations는 이 절차를 공식화한 지침이다.

8. 특수 해석 기법

양력과 항력 해석에서는 다음과 같은 특수 기법이 추가적으로 활용된다. 첫째, 오버셋 격자(overset grid, Chimera) 기법은 움직이는 물체와 정지 물체를 각기 구축한 격자에서 해석하고 상호 보간하는 방법으로, 로터와 기체 상호작용, 멀티로터 해석에 적용된다. 둘째, 액추에이터 디스크(actuator disk) 또는 액추에이터 라인(actuator line) 모형은 로터 또는 프로펠러를 분포형 힘원으로 근사하여 회전 블레이드를 전체 해상하지 않고도 후류와 주변 기체에 미치는 영향을 계산한다. 셋째, 적응 격자 세분화(adaptive mesh refinement, AMR)는 와류, 박리 영역 등 해상이 필요한 구역에서 국지적으로 격자를 세밀화한다. 이러한 기법은 Johnson의 Helicopter Theory(Princeton University Press, 1980)에서 다룬 회전익 해석의 전통과 NASA OVERFLOW, Stanford SU2, OpenFOAM 등 공개 및 상용 CFD 코드의 최신 기능을 통해 실현되고 있다.

9. 드론 및 로봇 응용

소형 무인기, 멀티로터, 협동 비행 로봇의 CFD 해석에서는 다음의 고려가 필요하다. 첫째, 저 $Re$ 조건에서의 전이 예측을 위해 전이 모형과 세밀한 경계층 격자가 필수이다. 둘째, 프로펠러와 기체의 상호작용을 정확히 재현하기 위해 오버셋 격자 또는 slide mesh 기법이 사용된다. 셋째, 지면 효과, 벽 효과, 실내 비행 환경의 경계 조건을 올바르게 반영해야 한다. 이러한 응용 사례와 수치 접근의 세부 비교는 Ventura Diaz와 Yoon의 High-Fidelity Computational Aerodynamics of Multi-Rotor Unmanned Aerial Vehicles(56th AIAA Aerospace Sciences Meeting, AIAA Paper 2018-1266, 2018) 등에서 확인할 수 있다.

10. 출처

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.
Langtry, R. B., and Menter, F. R. “A Correlation-Based Transition Model Using Local Variables—Part I: Model Formulation.” Journal of Turbomachinery, vol. 128, no. 3, 2006.
AIAA Standards. Guide for the Verification and Validation of Computational Fluid Dynamics Simulations. AIAA G-077-1998.
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.

11. 버전

v1.0 (2026-04-17)