대류-확산 방정식

전산유체역학 수치해석 마지막 수정: 2026-07-13 04:35:58

1. 개요[편집]

대류-확산 방정식(Convection-Diffusion equation)은 어떤 스칼라량 $\phi$ (온도·농도·화학종 질량분율 등)가 유동에 실려 이동(대류)하는 동시에 분자 확산으로 퍼져 나가는 과정을 기술하는 스칼라 수송 방정식이다. 전산유체역학의 거의 모든 스칼라 수송 계산의 기본 골격이며, 형태는 다음과 같다.

\frac{\partial \phi}{\partial t} + \mathbf{u}\cdot\nabla\phi = \nabla\cdot(\Gamma\nabla\phi) + S

좌변 둘째 항 $\mathbf{u}\cdot\nabla\phi$ 가 대류(convection), 우변 첫째 항 $\nabla\cdot(\Gamma\nabla\phi)$ 가 확산(diffusion), $S$ 가 생성/소멸 항이다. 겉보기엔 순해 보이지만, 이 방정식은 CFD를 처음 배우는 사람에게 “수치 진동”과 “거짓 확산”이라는 두 개의 통과의례를 안겨주는 것으로 악명 높다.¹

2. 대류와 확산의 성격 차이[편집]

이 방정식이 까다로운 근본 이유는 두 항의 수학적 성질이 정반대라는 데 있다.

확산 항은 등방적(isotropic)이다. 정보가 사방으로 대칭하게 퍼지며, 급격한 기울기를 완만하게 문질러 준다. 수학적으로 타원형/포물형 성질을 띠어 수치적으로 온순하다. 중심차분(차분 도식) 같은 대칭 이산화와 궁합이 좋다.
대류 항은 방향성(directional)을 가진다. 정보가 유동 방향(상류→하류)으로만 흐르는 쌍곡형 성질이라, 상류에서 무슨 일이 벌어졌는지가 하류를 결정한다. 이 방향성을 무시하고 대칭 이산화를 쓰면 반드시 탈이 난다.

이 둘의 상대적 세기를 재는 무차원수가 페클레 수 $Pe = UL/\Gamma$ 다. $Pe$ 가 작으면 확산 지배, 크면 대류 지배 문제가 된다. 그리고 CFD의 온갖 골칫거리는 대부분 $Pe$ 가 클 때 튀어나온다.

3. 셀 페클레 수와 수치 진동[편집]

이산화 후에는 격자 셀 하나 크기 $\Delta x$ 를 기준으로 한 셀 페클레 수(cell Péclet number)가 관건이 된다.

Pe_{\text{cell}} = \frac{U\,\Delta x}{\Gamma}

1차원 정상 대류-확산 문제를 중심차분으로 이산화하면, $Pe_{\text{cell}} > 2$ 일 때 계수 행렬이 대각 우세를 잃으면서 해에 **비물리적 진동(wiggle)**이 나타난다. 매끄러워야 할 온도장이 톱니처럼 튀거나, 음의 농도 같은 물리적으로 불가능한 값이 등장한다.² 이는 폰 노이만 안정성 해석이나 계수 부호 분석으로 예측 가능한 결정론적 현상이지, 버그가 아니다.

해결책은 두 갈래다. 격자를 충분히 조밀하게 만들어 $Pe_{\text{cell}} < 2$ 를 강제하거나(고 $Pe$ 유동에서는 격자 수가 폭발한다), 대류 항 이산화 자체를 방향성을 존중하는 방식으로 바꾸는 것.

4. 상류차분과 거짓 확산[편집]

가장 단순한 처방이 **상류차분(upwind scheme)**이다. 셀 경계의 $\phi$ 값을 하류가 아닌 상류 쪽 값으로 취한다. 물리적 정보 흐름 방향을 존중하므로 셀 페클레 수와 무관하게 진동 없는(유계, bounded) 해를 준다. 안정성 하나는 확실하다.

문제는 공짜가 아니라는 것. 1차 상류차분은 절단 오차의 지배항이 확산과 똑같은 $\partial^2\phi/\partial x^2$ 형태라서, 방정식에 존재하지도 않는 **가짜 확산(false diffusion / numerical diffusion)**을 슬쩍 끼워 넣는다.³ 유동이 격자선과 비스듬히 흐를 때 특히 심해서, 날카로워야 할 경계층이나 화염면이 실제보다 뭉개져 나온다. “진동은 없앴지만 해상도를 팔아먹은” 셈이다.

5. 고차 도식: QUICK와 TVD[편집]

정확도와 안정성을 동시에 잡으려는 시도가 고차 대류 도식이다.

QUICK(Quadratic Upstream Interpolation for Convective Kinematics) — 상류 쪽 두 점과 하류 한 점을 2차 포물선으로 보간해 셀 경계값을 구한다. 3차 정확도라 거짓 확산이 크게 줄지만, 급격한 기울기 근처에서는 여전히 약간의 언더슈트/오버슈트가 생길 수 있다.
TVD(Total Variation Diminishing) 계열 — 해의 전변동(total variation)이 시간이 지나도 증가하지 않도록 강제하는 부류. 유량 제한자(flux limiter, 예: minmod·van Leer·superbee)를 써서 매끄러운 영역에서는 고차 정확도를, 불연속 근처에서는 상류차분으로 매끄럽게 전환한다. 진동을 원천 봉쇄하면서 해상도를 지키는, 압축성 유동의 충격파 포착에서도 표준이 된 접근이다.

실무에서는 유한체적법 기반 상용/오픈소스 코드(OpenFOAM, ANSYS Fluent 등)가 이 도식들을 옵션으로 제공한다. “일단 상류차분으로 수렴시키고, 안정되면 고차 도식으로 갈아타 정확도를 확보한다”가 흔한 국룰 전략이다.

6. 활용과 확장[편집]

대류-확산 방정식은 CFD의 스칼라 수송 전반을 관장한다. 열전달 해석의 에너지 방정식(온도 수송), 다상유동의 상 분율 이동, 연소·오염물 확산의 화학종 수송, 심지어 난류 모델링의 $k$ · $\epsilon$ · $\omega$ 수송 방정식까지 전부 이 형태다. 나비에-스토크스 방정식의 운동량 방정식도 $\phi=\mathbf{u}$ 인 벡터 버전의 대류-확산으로 읽을 수 있어서, 여기서 배운 이산화 교훈이 유동 자체의 해석에 그대로 이어진다. 시간 항이 있으면 크랭크-니콜슨법 같은 시간 적분과 결합해 비정상 문제로 확장된다.

7. 관련 문서[편집]

대류-확산 방정식은 CFD 교과서(파탄카, 페릭 등)에서 나비에-스토크스보다 먼저 등장한다. 이 스칼라 방정식 하나에서 이산화의 거의 모든 이슈(보존성·유계성·정확도·안정성)가 응축되어 나타나기 때문. 여기서 좌절하면 본편은 더하다. ↩
흔히 “위글(wiggle)“이라 부른다. 처음 본 사람은 대개 코드 버그를 의심하며 며칠을 태우다가, 교과서에서 “이건 정상입니다”라는 문장을 발견하고 허탈해한다. 통과의례. ↩
거짓 확산은 존재하지 않는 점성을 공짜로 얻는 것과 같다. 해가 안정적으로 보여서 초심자는 오히려 좋아하지만, 실제로는 물리를 뭉개고 있는 것이다. “수렴 잘 되는 결과일수록 의심하라”는 격언의 대표 사례. ↩