Chapter 14

COMSOL 시뮬레이션 응용

CFD 시뮬레이션 워크플로, Lid-driven cavity 벤치마크, 검증 및 확인(V&V) 방법론을 다룹니다.

COMSOL 벤치마크 검증 Ghia et al.

CFD 시뮬레이션 워크플로

CFD 시뮬레이션의 표준 워크플로는 (1) 기하 모델링, (2) 격자 생성, (3) 물리 설정(지배 방정식, 경계 조건, 물성치), (4) 수치 해석, (5) 후처리 및 분석의 5단계로 구성됩니다. COMSOL Multiphysics는 이 전 과정을 하나의 GUI 환경에서 수행할 수 있는 상용 FEM 기반 다물리 시뮬레이션 소프트웨어입니다.

격자 생성 시 중요한 고려사항으로는 경계층 격자(boundary layer mesh), 격자 품질(skewness, aspect ratio), 적응적 격자 세밀화(adaptive mesh refinement) 등이 있습니다. 물리 설정에서는 난류 모델 선택($k$-$\epsilon$, $k$-$\omega$ SST 등), 벽 함수 사용 여부, 정상/비정상 해석 선택이 결과에 큰 영향을 미칩니다.

Lid-Driven Cavity 벤치마크

Lid-driven cavity (뚜껑 구동 공동)는 CFD의 가장 표준적인 벤치마크 문제입니다. 정사각형 공동의 상단 벽이 일정 속도 $U$로 움직이고, 나머지 세 벽은 고정(no-slip)입니다. Ghia et al. (1982)의 고해상도 수치 데이터가 표준 검증 자료로 널리 사용되며, Re = 100에서 1000까지의 속도 프로파일이 제공됩니다.

낮은 Re에서는 공동 중앙에 하나의 주 와류(primary vortex)만 형성됩니다. Re가 증가하면 주 와류의 중심이 공동 중앙 쪽으로 이동하고, 하단 모서리에 2차 와류(secondary vortex)가 나타납니다. Re = 1000 이상에서는 3차 와류도 관찰됩니다.

$$Re = \frac{U L}{\nu}$$
Figure 14.1 Lid-Driven Cavity 속도 프로파일 (Ghia et al. 1982)
Reynolds number Re = 100
Re 100
주 와류 위치
umin
Q. Lid-driven cavity에서 Re가 증가하면 주 와류의 중심은 어디로 이동하는가?
낮은 Re에서 주 와류 중심은 뚜껑(이동 벽) 가까이에 위치하지만, Re가 증가하면 관성 효과가 강해져 와류가 공동의 기하학적 중심 쪽으로 이동합니다. Ghia et al. 데이터에서 이 경향을 정량적으로 확인할 수 있습니다.

검증과 확인 (V&V)

검증(Verification)과 확인(Validation)은 CFD 결과의 신뢰성을 확보하기 위한 체계적 절차입니다. 검증은 "방정식을 올바르게 풀고 있는가?"(코드 검증, 해 검증)에 대한 답이고, 확인은 "올바른 방정식을 풀고 있는가?"(모델이 현실을 얼마나 잘 반영하는가)에 대한 답입니다.

코드 검증은 해석해가 있는 문제(MMS: Method of Manufactured Solutions 포함)로 수행하며, 이론적 수렴 차수를 확인합니다. 해 검증은 격자 수렴 검증(GCI)으로 이산화 오차를 정량화합니다. 모델 확인은 실험 데이터와의 비교로, 실험 불확실성과 모델 불확실성을 모두 고려해야 합니다.

$$GCI = \frac{F_s \, |e|}{r^p - 1}$$
Q. "올바른 방정식을 풀고 있는가?"에 답하는 절차는?
검증(Verification)은 "방정식을 올바르게 풀고 있는가?"(수치적 정확성)에, 확인(Validation)은 "올바른 방정식을 풀고 있는가?"(물리적 적합성)에 답합니다. 격자 수렴 연구는 검증의 일부입니다.

핵심 정리

01
워크플로
기하 모델링 -> 격자 생성 -> 물리 설정 -> 해석 -> 후처리. 각 단계에서의 적절한 선택이 결과 품질을 결정합니다.
02
벤치마크
Lid-driven cavity는 Re에 따른 와류 구조 변화를 보여주는 CFD의 표준 검증 문제입니다. Ghia et al. 데이터가 기준입니다.
03
V&V
Verification은 수치 정확성, Validation은 물리적 적합성. 둘 다 갖춰야 시뮬레이션 결과를 신뢰할 수 있습니다.