19.26 전산 유체역학(CFD)의 기본 원리

전산 유체역학(Computational Fluid Dynamics, CFD)은 유체의 보존 법칙을 수치적 방법으로 해석하여 속도장, 압력장, 온도장, 농도장 등을 계산하는 학문 분야이다. 해석 대상 유동의 지배 방정식은 편미분 방정식 체계로 기술되며, 이를 이산화하여 대수 방정식 체계로 변환한 뒤 컴퓨터를 이용하여 해를 구한다. 본 절에서는 CFD의 이론적 기반, 이산화 기법, 해법 절차, 검증 과정 등 핵심 원리를 체계적으로 정리한다.

1. CFD의 정의와 역할

실험 유체역학과 이론 유체역학이 역사적으로 선행한 해석 수단이었던 반면, CFD는 1960년대 이후 컴퓨터 성능의 향상에 따라 독립적인 제3의 해석 수단으로 자리 잡았다. CFD는 복잡한 기하학적 형상과 경계 조건에 대해 비용 효율적 해석을 가능하게 하며, 실험으로는 측정이 어려운 내부 유동장 정보를 제공한다는 장점을 갖는다. 그러나 수치 오차, 모델 근사, 격자 의존성 등의 제약이 존재하므로 이론 및 실험 결과와의 교차 검증이 필수적이다.

2. 지배 방정식과 보존 형태

일반적인 CFD 해석에서는 질량, 운동량, 에너지 보존 법칙을 연립하여 해석한다. 비압축성 유동의 경우

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

$\dfrac{\partial \mathbf{u}}{\partial t} + (\mathbf{u} \cdot \nabla) \mathbf{u} = -\dfrac{1}{\rho}\nabla p + \nu \nabla^2 \mathbf{u}$

이, 압축성 유동의 경우 연속, 운동량, 에너지 방정식과 상태 방정식이 결합된 오일러 혹은 나비에-스토크스 방정식 체계가 지배 방정식이 된다. 수치 해법에서는 보존 법칙의 적분 형식을 기반으로 하는 보존 형태(conservative form)가 선호되며, 이는 수치해가 질량, 운동량, 에너지의 국소 보존을 엄밀히 만족하도록 보장한다.

3. 이산화 기법

편미분 방정식을 이산화하는 대표적 기법은 유한 차분법(Finite Difference Method, FDM), 유한 체적법(Finite Volume Method, FVM), 유한 요소법(Finite Element Method, FEM), 스펙트럴법(Spectral Method)이 있다.

유한 차분법은 미분 연산자를 테일러 전개에 기반한 차분 근사로 대체하는 가장 단순한 기법이며, 구조화 격자에 주로 적용된다. 유한 체적법은 제어 체적에 대한 적분 형태의 보존 법칙을 이산화하므로 보존성이 자연스럽게 보장되며, 복잡한 형상과 비구조 격자에 적용이 용이하여 상용 CFD 코드에서 가장 널리 사용된다. 유한 요소법은 기저 함수에 의한 약형(weak form) 전개를 기반으로 하며, 복잡한 기하학적 형상 처리와 경계 조건 적용이 유연하다. 스펙트럴법은 고차 정확도의 해를 제공하지만 간단한 형상에 제한적으로 적용된다.

4. 유한 체적법의 기본 원리

유한 체적법은 임의의 제어 체적 $V$ 에 대해 보존 법칙의 적분 형태를 적용한다. 일반적인 스칼라 수송 방정식은

$\dfrac{\partial}{\partial t}\int_V \rho \phi \, dV + \oint_{\partial V} \rho \phi \mathbf{u} \cdot d\mathbf{A} = \oint_{\partial V} \Gamma_\phi \nabla \phi \cdot d\mathbf{A} + \int_V S_\phi \, dV$

로 표현되며, 좌변의 시간 변화항, 대류항과 우변의 확산항, 생성항으로 구성된다. 각 항의 면적분과 체적분을 이산화하기 위해 면 값 보간 방식(예: 상류 차분, 중앙 차분, QUICK, TVD)과 수치 플럭스 정의(예: Roe, HLLC, AUSM)가 사용된다. 이산화된 방정식은 이웃 셀의 값과 선형 관계를 형성하며, 결과적으로 연립 대수 방정식 시스템을 얻는다.

19.26.5 시간 적분 기법

시간 종속 문제에서 시간 도함수의 이산화는 명시적(explicit)과 음해적(implicit) 기법으로 구분된다. 명시적 기법은 각 시간 단계에서 이전 시점의 값만을 이용하여 현재 값을 직접 계산하며, 구현이 간단하나 시간 단계가 Courant-Friedrichs-Lewy(CFL) 조건에 의해 제한된다. 음해적 기법은 현재 시점의 값을 포함하는 방정식을 반복 해법을 통해 푸는 방식으로, 무조건 안정성을 확보할 수 있으나 계산 비용이 증가한다. 대표적 시간 적분법으로는 전진 오일러, 후진 오일러, Crank-Nicolson, Runge-Kutta, Backward Differentiation Formula(BDF) 등이 있다.

CFL 조건은 명시적 기법의 안정성 판정 기준으로서, 1차원 대류 문제의 경우

$C = \dfrac{u \Delta t}{\Delta x} \le C_{\max}$

로 정의된다. $C_{\max}$ 는 이산화 방식과 기법에 따라 결정되는 값이다.

5. 압력-속도 결합 알고리즘

비압축성 유동에서는 연속 방정식이 시간 미분을 포함하지 않아 압력과 속도의 직접 결합 해법이 곤란하다. 이를 해결하기 위해 압력 보정 알고리즘이 개발되었다. Patankar와 Spalding이 제안한 SIMPLE(Semi-Implicit Method for Pressure-Linked Equations) 알고리즘은 추정 압력장으로부터 임시 속도를 계산하고, 연속 방정식을 이용하여 압력 보정 방정식을 유도한 후 속도와 압력을 갱신한다. 이후 SIMPLEC, PISO, Coupled solver 등 개선된 알고리즘들이 개발되었으며, 각각 반복 수렴 속도와 안정성 측면에서 고유한 장점을 가진다.

6. 난류 모델링

대부분의 공학적 유동은 난류 영역에 속하며, 나비에-스토크스 방정식을 직접 수치 해석(Direct Numerical Simulation, DNS)으로 푸는 것은 레이놀즈 수가 높은 문제에서 계산 비용이 과도하여 제한적이다. 이에 따라 시간 평균 또는 공간 필터링을 통해 축소된 방정식 체계를 해석하는 기법이 개발되어 있다.

Reynolds-Averaged Navier-Stokes(RANS) 접근법은 유동 변수를 평균 성분과 변동 성분의 합으로 분해하고, 평균 운동량 방정식에 나타나는 Reynolds 응력 텐서를 난류 모델로 닫는다. 대표적 모델로는 $k$ - $\varepsilon$ 모델, $k$ - $\omega$ 모델, $k$ - $\omega$ SST 모델, Spalart-Allmaras 모델 등이 있다. Large Eddy Simulation(LES)은 큰 스케일의 난류 구조를 직접 해석하고 작은 스케일은 서브그리드 모델로 근사하는 기법이며, Detached Eddy Simulation(DES)과 같은 하이브리드 기법도 널리 사용된다.

7. 경계 조건과 초기 조건

지배 방정식의 수학적 성격(타원형, 포물형, 쌍곡형)에 따라 적절한 경계 조건이 요구된다. 대표적 경계 조건으로는 속도 입구(velocity inlet), 압력 출구(pressure outlet), 벽면 조건(no-slip wall, slip wall), 대칭 조건(symmetry), 주기 조건(periodic), 원거리 조건(far-field) 등이 있다. 난류 경계층을 해석할 때에는 벽 함수(wall function)를 사용하여 점성 하부층(viscous sublayer) 해석을 우회하거나, 충분히 조밀한 격자를 통해 직접 해석(low-Reynolds wall treatment)을 수행한다. 첫 격자점의 무차원 거리 $y^+$ 는 난류 해석의 품질을 가늠하는 핵심 지표이다.

8. 수렴성, 일관성, 안정성

Lax 정리는 선형 문제에 대해 일관성(consistency)과 안정성(stability)이 성립하면 수렴성(convergence)이 보장됨을 명시한다. 일관성은 이산화 방정식이 격자 간격과 시간 간격이 0으로 수렴함에 따라 원래의 편미분 방정식으로 환원됨을, 안정성은 수치 해가 유한 시간 구간에서 유계 상태를 유지함을, 수렴성은 이산 해가 정확한 해로 수렴함을 각각 의미한다. 비선형 문제에서는 이 정리의 직접적 적용에 한계가 있으나, 이들 세 속성은 여전히 CFD 해석의 기본 지침으로 작용한다.

9. 수치 오차와 검증

CFD 결과에는 이산화 오차, 반복 수렴 오차, 모델링 오차, 프로그래밍 오차 등 다양한 오차가 포함된다. 이를 정량적으로 평가하기 위해 미국 기계학회(ASME) 등의 가이드라인은 Verification and Validation(V&V) 절차를 제시하고 있다. Verification은 코드가 지배 방정식을 정확히 해석하는지를 검사하며, Grid Convergence Index(GCI) 기반의 Richardson 외삽이 대표적 기법이다. Validation은 수치 결과가 실제 물리 현상을 얼마나 잘 반영하는지를 실험 결과와의 비교를 통해 평가한다.

10. 병렬화와 고성능 계산

CFD는 막대한 계산 자원을 요구하므로 분산 메모리 병렬화(MPI), 공유 메모리 병렬화(OpenMP), 그래픽 처리 장치(GPU) 가속 등의 고성능 계산 기법이 필수적으로 도입된다. 도메인 분할(domain decomposition)을 기반으로 한 병렬화가 가장 일반적이며, 할로 교환(halo exchange)을 통해 이웃 도메인 간 경계 정보를 주고받는다.

11. 로봇 공학에서의 응용과 의의

로봇 공학에서 CFD는 다음과 같은 응용에 기여한다. 무인 항공기의 동체 및 날개 공력 해석을 통해 양력, 항력, 모멘트 계수가 산출되며, 이는 비행 제어기 설계의 입력이 된다. 수중 로봇의 선체 항력 최적화, 프로펠러 성능 예측, 캐비테이션 해석에도 CFD가 활용된다. 소프트 로봇의 내부 유동 해석, 생체 모사 로봇의 추진 메커니즘 해석, 열관리 및 냉각 유동 해석 또한 주요 응용 분야이다. 본 절에서 기술한 CFD의 원리는 후속되는 격자 생성, 수치 해법, 시각화, 검증 절차의 이론적 토대를 제공한다.

12. 출처

Ferziger, J. H., Perić, M., and Street, R. L., Computational Methods for Fluid Dynamics, 4th ed., Springer, 2020.
Versteeg, H. K., and Malalasekera, W., An Introduction to Computational Fluid Dynamics: The Finite Volume Method, 2nd ed., Pearson, 2007.
Anderson, J. D., Computational Fluid Dynamics: The Basics with Applications, McGraw-Hill, 1995.
Hirsch, C., Numerical Computation of Internal and External Flows, 2nd ed., Butterworth-Heinemann, 2007.
Patankar, S. V., Numerical Heat Transfer and Fluid Flow, Hemisphere, 1980.
Wilcox, D. C., Turbulence Modeling for CFD, 3rd ed., DCW Industries, 2006.

13. 버전

v1.0