직접수치모사

시뮬레이션 난류 마지막 수정: 2026-07-09 04:23:45

직접수치모사 Direct Numerical Simulation
약칭	DNS
분야	난류 모델링 × CFD
핵심	난류 모델 없이 모든 스케일을 직접 해상
격자 스케일링	$N \sim Re^{9/4}$

1. 개요[편집]

직접수치모사(Direct Numerical Simulation, DNS)는 난류 모델을 일절 쓰지 않고 나비에-스토크스 방정식을 있는 그대로 이산화해, 가장 큰 와류부터 가장 작은 소산 와류까지 유동의 모든 시공간 스케일을 직접 해상하는 난류 해석 기법이다. 난류 모델링의 스펙트럼에서 가장 오른쪽 끝, 즉 “모델링 0%, 계산 100%“의 극단에 위치한다.

RANS는 다 평균 내고, LES는 작은 놈만 모델링하지만, DNS는 아무것도 봐주지 않는다. 콜모고로프 미세 스케일에서 벌어지는 점성 소산까지 격자로 낱낱이 잡아낸다. 그래서 결과는 사실상 “수치 실험(numerical experiment)“에 가깝다. 대신 그 대가로 격자수와 계산량이 천문학적으로 불어난다. DNS는 정확함의 대명사이자 동시에 비용의 대명사다.¹

2. 모든 스케일을 해상한다는 것[편집]

난류는 다양한 크기의 와류가 겹겹이 쌓인 계층 구조다. 큰 와류가 에너지를 받아 점점 작은 와류로 쪼개지며 에너지를 넘겨주고(에너지 캐스케이드), 최종적으로 가장 작은 스케일에서 점성에 의해 열로 소산된다. DNS의 조건은 명확하다. 가장 큰 와류(적분 스케일 $L$ )부터 가장 작은 와류(콜모고로프 스케일 $\eta$ )까지를 하나의 격자와 하나의 시간 간격으로 전부 담아야 한다.

격자 간격은 $\eta$ 수준으로 촘촘해야 하고,
도메인은 $L$ 을 여러 개 담을 만큼 커야 하며,
시간 간격은 CFL 조건을 지키면서 가장 빠른 소산 과정을 좇아야 한다.

이 세 조건을 동시에 만족시키는 순간, 격자수는 걷잡을 수 없이 커진다.

3. 콜모고로프 스케일과 격자수 폭발[편집]

콜모고로프의 1941년 이론에 따르면, 가장 작은 소산 와류의 크기 $\eta$ 는 점성 $\nu$ 와 에너지 소산율 $\varepsilon$ 으로 결정된다.

\eta = \left( \frac{\nu^3}{\varepsilon} \right)^{1/4}

적분 스케일 $L$ 과 콜모고로프 스케일 $\eta$ 의 비율은 레이놀즈수로 이렇게 정리된다.

\frac{L}{\eta} \sim Re^{3/4}

즉 한 방향에 필요한 격자점 수가 $Re^{3/4}$ 에 비례한다. 3차원이니까 세제곱하면,

N \sim \left( \frac{L}{\eta} \right)^{3} \sim Re^{9/4}

그 유명한 $Re^{9/4}$ 스케일링이다. 시간 간격까지 고려하면 총 연산량은 대략 $Re^{3}$ 에 육박한다.² 레이놀즈수를 10배 올리면 격자는 약 178배, 연산량은 1000배 늘어난다는 뜻. 이래서 DNS는 레이놀즈수 앞에서 무릎을 꿇는다.

4. 왜 산업에서는 못 쓰는가[편집]

산업 유동의 레이놀즈수는 보통 $10^6$ 을 우습게 넘긴다. 자동차 외부 유동, 항공기 날개, 가스터빈 연소기 전부 그렇다. 이걸 $Re^{9/4}$ 에 대입하면 격자수가 $10^{13}$ 단위로 나온다. 지구상 어떤 슈퍼컴퓨터로도 이런 문제를 실무 일정 안에 풀 수 없다.³

그래서 DNS는 대개 다음 조건에서만 가능하다.

레이놀즈수가 낮거나 중간인 유동(채널 유동, 균질 등방성 난류, 제트 초기 영역 등)
단순한 형상(주기 경계가 걸린 정육면체, 평판 채널 같은)
학술 연구 목적, 그리고 넉넉한 슈퍼컴 할당량

산업 CFD가 RANS와 LES로 타협하는 이유가 바로 여기 있다. DNS는 “정답을 알지만 감당 못 하는” 방법이다.

5. 그럼에도 학술적 가치[편집]

돈도 시간도 잡아먹는 DNS를 왜 계속 하느냐. 답은 하나다. 다른 어떤 방법으로도 얻을 수 없는 완전한 데이터를 주기 때문이다.

난류 모델의 검증 기준(reference): RANS와 LES의 SGS 모델이 얼마나 맞는지 비교할 “정답지”가 필요한데, 그 정답지를 DNS가 제공한다. 실험으로는 측정 불가능한 압력 변동, 소산율, 응력 텐서 전 성분을 시공간 전 영역에서 뽑아낼 수 있다.
난류 물리의 근본 이해: 벽 난류의 구조, 에너지 캐스케이드, 간헐성(intermittency), 전이(transition) 메커니즘 같은 난제를 들여다보는 현미경 역할.
모델 없는 진실: 검증 및 확인의 관점에서 DNS는 모델 오차가 원리적으로 없다. 남는 건 이산화 오차와 통계 수렴 오차뿐이라, 격자만 충분하면 “진짜 나비에-스토크스의 해”에 무한히 가까워진다.

실제로 DNS는 스펙트럴 방법 같은 고차 정확도 기법과 궁합이 좋다. 인공 소산 없이 소산 스케일까지 깨끗하게 잡아야 하기 때문에, 정확도가 낮은 유한체적법보다 스펙트럴/고차 유한차분이 선호된다.

6. 관련 문서[편집]

최초의 본격적인 DNS는 Orszag & Patterson(1972)의 균질 등방성 난류 계산과, Kim, Moin & Moser(1987)의 채널 유동 DNS로 꼽힌다. 후자는 지금도 난류 모델 검증의 고전 데이터셋으로 인용된다. 무려 40년 넘게 우려먹히는 중. ↩
정확히는 공간 $Re^{9/4}$ 에 시간 스텝 수 $Re^{3/4}$ 을 곱해 $Re^{3}$ 이 된다는 계산. 상수와 정의에 따라 지수가 조금씩 오르내리지만, “레이놀즈수의 3승 언저리”라는 감각은 변하지 않는다. ↩
항공기 전기체 DNS는 2080년대에나 가능하리라는 Spalart(2000)의 추정이 유명하다. 무어의 법칙이 식어가는 지금은 그마저도 낙관이었다는 소리를 듣는다. DNS 하는 사람은 레이놀즈수를 낮춰서 논문을 쓰는 법을 먼저 배운다. ↩