분리 와류 모사

1. 개요[편집]

분리 와류 모사 Detached Eddy Simulation
약칭	DES
분류	RANS–LES 하이브리드 난류 해석
제안	Spalart 외 (1997)
벽 근처	RANS로 처리
박리 후류	LES로 처리
주요 변종	DDES, IDDES

분리 와류 모사(Detached Eddy Simulation, DES)는 벽 근처의 얇은 경계층은 RANS로 싸게 처리하고, 벽에서 떨어져 나와 큰 와류가 판치는 유동박리 후류 영역은 대와류 모사(LES)로 해상하는 하이브리드 난류 해석 기법이다. 1997년 필립 스팔라트(Philippe Spalart)가 “항공기 전체를 LES 하는 건 2080년에나 가능하니, LES가 정말 필요한 곳에만 LES를 쓰자”는 지극히 현실적인 타협안으로 제안했다.¹

핵심 발상은 이렇다. 경계층에서 LES가 비싼 이유는 벽에 붙은 자잘한 난류 스트리크를 해상하느라 격자가 폭발하기 때문이다. 그런데 이 부분은 RANS가 수십 년간 잘 튜닝해온 영역이다. 반대로 박리된 후류의 거대한 비정상 와류는 RANS가 젬병이고 LES가 잘한다. 그러면 각자 잘하는 데를 맡기면 되지 않나? 그 발칙한 아이디어가 DES다.

2. RANS와 LES 사이 어디쯤[편집]

전산유체역학의 난류 해석은 계산비용과 정확도의 스펙트럼 위에 늘어서 있다. 왼쪽 끝에 만능이지만 천문학적으로 비싼 직접수치모사(DNS), 오른쪽 끝에 싸지만 시간평균만 주는 RANS가 있고, 그 사이 대와류 모사(LES)가 있다.

기법	해상 범위	격자 규모	벽 근처 비용
DNS	모든 스케일	$Re^{9/4}$ 급	살인적
LES	큰 와류 해상, 작은 건 모델	크다	매우 큼
DES	벽=RANS, 박리후류=LES	중간	감당 가능
RANS	전부 모델링	작다	저렴

DES는 이 표에서 정확히 LES와 RANS 사이 틈을 메운다. 벽 근처의 콜모고로프 스케일급 미세 난류를 해상하는 짐을 RANS에 떠넘김으로써, 순수 LES라면 벽 해상에 쏟아야 할 격자를 후류 해상에 몰아줄 수 있다. 항공기 착륙장치, 후방 계단, 실속 날개처럼 대규모 박리가 지배하는 문제에서 가성비가 폭발한다.

3. 작동 원리 — 길이 척도 바꿔치기[편집]

원조 DES(DES97)는 스팔라트-알마라스 모델(S-A) 위에서 정의된다. S-A 모델의 소산항에는 벽까지의 거리 $d$가 길이 척도로 들어가는데, DES는 이 $d$를 다음처럼 격자 크기와 경쟁시킨다.

$$\tilde{d} = \min\left(d,\; C_{DES}\,\Delta\right)$$

여기서 $\Delta$는 국소 격자 크기(보통 세 방향 중 최대), $C_{DES}\approx 0.65$는 상수다. 논리는 간단하다.

• 벽 근처에서는 $d$가 작아서 $\tilde{d}=d$ → 모델이 RANS 모드로 작동.
• 벽에서 멀고 격자가 촘촘한 후류에서는 $C_{DES}\Delta < d$ → $\tilde{d}=C_{DES}\Delta$ → 소산이 격자 크기에 묶이며 사실상 LES의 스마고린스키형 아공간 모델처럼 작동.

즉 하나의 방정식이 격자와 벽거리의 비교만으로 자기가 RANS인지 LES인지를 스스로 판단한다. 별도의 도메인 분할 없이 매끄럽게 모드가 전환된다는 점이 DES의 우아함이자, 동시에 온갖 사고의 원천이다.

4. DDES와 IDDES — 사고 수습의 역사[편집]

원조 DES97에는 치명적 함정이 있었다. 격자를 “잘못” 촘촘하게 깔면 두꺼운 경계층 안에서 $C_{DES}\Delta < d$ 조건이 걸려버려, RANS로 처리돼야 할 벽 근처가 성급하게 LES 모드로 전환된다. 그런데 그 격자는 LES를 지탱할 만큼은 촘촘하지 않아서, 모델링된 난류 점성은 줄었는데 해상된 난류는 아직 안 생긴 어정쩡한 상태가 된다. 결과는 모델 응력 소진(Modeled Stress Depletion, MSD) → 벽 전단응력이 비현실적으로 낮아지고, 심하면 격자만 바꿨는데 멀쩡하던 경계층이 유동박리하는 **격자 유발 박리(Grid-Induced Separation, GIS)**까지 벌어진다.²

이를 막으려고 2006년 등장한 것이 **DDES(Delayed DES)**다. 경계층 안에 있는지를 감지하는 실드 함수 $f_d$를 도입해, 벽 근처에서는 격자가 아무리 촘촘해도 RANS 모드를 강제로 붙잡아 둔다.

$$\tilde{d} = d - f_d\,\max\!\left(0,\; d - C_{DES}\Delta\right)$$

한 발 더 나아간 **IDDES(Improved DDES)**는 여기에 벽 모델 LES(WMLES) 기능을 얹어, 유입 난류가 이미 존재할 때 벽 근처에서도 LES가 작동하도록 확장했다. 오늘날 상용 코드에서 “DES 돌린다”고 하면 사실상 DDES나 IDDES를 쓰는 경우가 대부분이다. 원조 DES97을 날것으로 쓰는 건 이제 지뢰밭을 맨발로 걷는 격.

5. Grey Area 문제[편집]

DES 계열의 영원한 골칫거리가 **그레이 에어리어(grey area)**다. 유동이 RANS 영역(모델링된 난류만 있음)에서 LES 영역(해상된 난류 필요)으로 넘어가는 전이 구간에서, 해상된 난류 변동이 즉시 발달하지 못하고 한참 뜸을 들이는 문제다. 얇은 전단층이 박리 직후 아직 큰 와류로 말려들기 전 구간이 특히 취약하다. 이 구간에서 모델은 “나는 이제 LES야”라고 선언했지만 정작 격자 위엔 해상할 난류 구조가 아직 없어서, 난류 응력이 실종된 회색 지대가 생긴다.

완화책으로는 아공간 난류 점성을 국소 와도로 낮춰 전단층의 켈빈-헬름홀츠 불안정을 빨리 촉발시키는 저소산(low-dissipation) 도식, $\Delta$를 세 방향 최대가 아니라 와류 방향을 고려해 재정의하는 $\tilde{\Delta}_\omega$ 척도, 그리고 아예 유입면에 합성 난류를 주입하는 방법 등이 쓰인다. 그럼에도 “박리 위치가 명확한 문제엔 DES가 잘 듣고, 박리가 완만하거나 전이가 지배적인 문제엔 여전히 조심하라”는 게 현장의 국룰이다.

6. 관련 문서[편집]

• 난류 모델링 · 대와류 모사 · 직접수치모사
• 스팔라트-알마라스 모델 · k-오메가 SST 모델
• 유동박리 · 경계층 · 콜모고로프 스케일
• 레이놀즈수 · 벽함수
• 전산유체역학

7. Footnotes[편집]

1. Spalart, P. R. et al. (1997). 원 논문의 문제의식이 노골적이다. 전기체 LES는 수십 년 뒤에나 가능하니, 지금 당장 쓸 수 있는 물건을 만들자는 실용주의의 승리. 스팔라트 본인이 스팔라트-알마라스 모델의 아버지이기도 하다. ↩

2. 격자를 더 좋게 만들었는데 결과가 더 나빠지는, 수치해석자의 멘탈을 가장 효과적으로 파괴하는 현상. “격자 독립성 검사 했더니 물리가 바뀌었어요”라는 문의의 상당수가 이 GIS다. ↩