담쾰러 수

유체역학 전산유체역학 마지막 수정: 2026-08-09 04:38:09

1. 개요[편집]

담쾰러 수 Damköhler number
기호	Da (Da_I, Da_II, …)
정의	유동(또는 확산) 시간 척도 / 화학 시간 척도
이름	게르하르트 담쾰러 (1936)
극한	Da → 0 동결 유동 · Da → ∞ 평형
대표 무대	반응기 설계, 연소, 재진입 화학, 촉매
수치적 함의	큰 Da = 강성 화학

다 섞이기 전에 타느냐, 다 타기 전에 섞이느냐. 그 승부의 배당률.

담쾰러 수(Damköhler number, $\mathrm{Da}$ )는 **유체가 어떤 영역에 머무르는 시간과 그 안에서 화학반응이 진행되는 데 걸리는 시간의 비를 재는 무차원수**다.

\mathrm{Da} = \frac{\tau_{\text{flow}}}{\tau_{\text{chem}}}

$\tau_{\text{flow}}$ 는 체류 시간·대류 시간·난류 회전 시간 같은 유동 쪽 척도이고, $\tau_{\text{chem}}$ 은 반응이 특성적으로 진행되는 시간이다(반응속도론). 이름은 독일 화학공학자 게르하르트 담쾰러(Gerhard Damköhler)가 1936년 반응로 성능을 무차원 그룹으로 정리한 연구에서 왔다.

$\mathrm{Da}$ 가 하는 일은 단순하다 — 화학과 수송 중 누가 율속인지 통보한다. 그런데 이 단순함이 실무에서 자주 배신하는데, 이유는 " $\mathrm{Da}$ "라고 불리는 무차원수가 하나가 아니기 때문이다.

2. 어느 Da를 말하는 건지부터 밝혀라[편집]

담쾰러 본인이 매긴 번호부터 여러 개였다. 1차 반응(속도상수 $k$ ), 대표 길이 $L$ , 속도 $u$ , 확산계수 $D$ 로 쓰면 가장 자주 쓰이는 둘은 이렇다.

기호	정의	재는 것
$\mathrm{Da_I}$	$kL/u$	반응 대 대류 수송
$\mathrm{Da_{II}}$	$kL^2/D$	반응 대 확산 수송
$\mathrm{Da_{III}}$	반응 발열 대 대류 열수송	열 버전
$\mathrm{Da_{IV}}$	반응 발열 대 전도 열수송	열 버전

$\mathrm{Da_I}$ 과 $\mathrm{Da_{II}}$ 는 독립이 아니다. 나눠 보면

\frac{\mathrm{Da_{II}}}{\mathrm{Da_I}} = \frac{uL}{D} = \mathrm{Pe}

즉 페클레 수가 둘을 잇는다. 대류가 지배적인 계에서는 두 값이 $\mathrm{Pe}$ 배, 즉 자릿수 단위로 벌어지므로 어느 쪽인지 안 밝히고 “Da가 크다”고 쓰면 문장이 아무 정보도 전달하지 않는다. 논문에서 $\mathrm{Da}$ 기호를 보면 정의식부터 찾는 것이 국룰인 이유다.¹

$\mathrm{Da_{II}}$ 는 촉매 공학에서 다른 이름으로도 산다. 다공성 촉매 입자 안에서 반응과 세공 확산이 겨루는 정도를 재는 틸레 계수(Thiele modulus) $\phi$ 가 $\phi^2 = kL^2/D_{\text{eff}}$ 로, 정확히 $\mathrm{Da_{II}}$ 다. 유효인자 $\eta$ 가 $\phi$ 의 함수로 떨어지고, $\phi \gg 1$ 이면 반응이 입자 표면 껍질에서만 일어나 안쪽 촉매는 놀고 있다는 결론이 나온다. 값비싼 백금을 입자 표면에만 담지하는 계란껍질(egg-shell) 촉매가 이 부등식의 산물이다.

3. 두 극한과 그 사이[편집]

극한에서는 계산이 오히려 쉬워진다.

$\mathrm{Da} \to 0$ — 동결(frozen) 유동. 반응이 너무 느려서 유체가 지나가는 동안 조성이 안 변한다. 화학 소스항을 지우고 대류-확산 방정식만 풀면 된다. 극초음속 노즐 팽창부에서 조성이 목 부근에서 “얼어붙는” 것이 이 극한이다.
$\mathrm{Da} \to \infty$ — 평형 / mixed-is-burnt. 반응이 무한히 빨라서 국소 조성이 즉시 화학평형에 도달한다. 이러면 종 수송 방정식 수십 개가 보존 스칼라 하나(혼합분율 $Z$ )와 룩업 테이블로 붕괴한다. 확산 화염의 고전 이론 전체가 이 극한 위에 서 있다.
$\mathrm{Da} \sim 1$ — 비평형. 둘 다 못 쓴다. 화학과 수송을 동시에 풀어야 하고, 착화·소염·오염물 생성처럼 실제로 알고 싶은 현상은 하필 전부 여기에 있다(극초음속 유동의 재진입 화학이 대표적이다).

반응기 공학에서는 이 비가 아예 전환율 공식으로 나온다. 1차 반응·정상 상태에서 완전혼합 흐름 반응기(CSTR)와 관형 반응기(PFR)의 전환율은

X_{\mathrm{CSTR}} = \frac{\mathrm{Da}}{1+\mathrm{Da}}, \qquad X_{\mathrm{PFR}} = 1 - e^{-\mathrm{Da}}

로, 여기서 $\mathrm{Da} = k\tau$ ( $\tau$ 는 체류 시간)다. 같은 $\mathrm{Da}=2$ 에서 CSTR은 67%, PFR은 86%. 반응기 형식 선택이 무차원수 하나로 정량화되는, 화학공학 교과서의 첫 장면이다.

4. 점화-소염 S자 곡선 — 접힘이 두 개[편집]

담쾰러 수가 단순한 비율 이상의 물건이 되는 순간은 발열이 반응속도를 되먹일 때다. 완전혼합 반응기(WSR/PSR)에서 체류 시간 $\tau$ (곧 $\mathrm{Da}$ )를 천천히 늘리며 정상 상태 온도를 그려 보자. 정상 상태 조건은 발열률과 유출에 의한 제열률이 같아지는 지점이고, 발열률은 아레니우스 방정식 때문에 온도의 지수 함수다.

결과는 유명한 S자 곡선이다. $\mathrm{Da}$ 가 작을 때는 저온 비반응 가지 하나뿐이다가, 어느 값을 넘으면 해가 세 개가 되고, 더 키우면 다시 고온 연소 가지 하나만 남는다. 곡선이 되접히는 두 지점이 각각 점화점과 소염점이다.

저온 가지를 따라 $\mathrm{Da}$ 를 키우다 점화점에 닿는 순간, 그 아래 정상 상태가 소멸해 계가 고온 가지로 튀어 오른다. 이것이 **자착화**다.
고온 가지에서 $\mathrm{Da}$ 를 줄이다 소염점에 닿으면 반대로 저온 가지로 떨어진다. 이것이 소염(blow-off)이다.
두 사건이 다른 $\mathrm{Da}$ 에서 일어나므로 이력(hysteresis) 이 생긴다. 한번 붙은 불이 붙일 때보다 낮은 유량까지 버티는 이유이고, 가스터빈 재점화 고도가 소염 고도보다 훨씬 낮은 이유다.

이 두 접힘은 비유가 아니라 문자 그대로 안장-마디 분기 한 쌍이다. 중간 가지는 안장형 불안정 해라 물리적으로 실현되지 않고, 곡선 위 어느 점에서도 자코비안이 특이해지는 곳이 정확히 두 접힘점이다. 그래서 이 곡선을 수치로 그릴 때는 $\mathrm{Da}$ 를 독립변수로 잡고 뉴턴법을 돌리면 접힘점에서 반드시 발산하고, 수치 연속법의 호장법으로 곡선을 매개변수화해야 중간 가지까지 통과할 수 있다.² 연소 이론의 S자 곡선과 분기 이론의 이력 고리가 같은 그림이라는 점은, 이 무차원수가 왜 단순 비율 이상인지를 잘 보여준다.

5. 난류 연소 체계 분류에서의 위치[편집]

난류 예혼합 화염에서는 $\tau_{\text{flow}}$ 를 난류 적분 척도의 회전 시간 $l_t/u'$ 로, $\tau_{\text{chem}}$ 을 층류 화염의 통과 시간 $\delta_L/S_L$ 로 잡는다.

\mathrm{Da} = \frac{l_t/u'}{\delta_L/S_L}

이 $\mathrm{Da}$ 와 카를로비츠 수 $\mathrm{Ka} = \tau_{\text{chem}}/\tau_\eta$ (콜모고로프 스케일의 시간)를 축으로 그린 것이 보르기 선도(Borghi diagram)다. 보통 $u'/S_L$ 대 $l_t/\delta_L$ 의 log-log 평면에 그리는데, 그러면 $\mathrm{Da}=1$ 은 기울기 1의 대각선이 되고 $\mathrm{Ka}=1$ (클리모프-윌리엄스 기준)이 또 하나의 선을 긋는다. 세 무차원수는 독립이 아니라

\mathrm{Re}_t = \mathrm{Da}^2\,\mathrm{Ka}^2

로 묶여 있어, 평면 위 한 점을 찍으면 나머지가 전부 결정된다.

구획은 대략 이렇게 나뉜다.

영역	조건	물리
주름진 화염편	$u'/S_L < 1$	난류가 화염면을 살짝 구길 뿐
굴곡진 화염편	$\mathrm{Ka} < 1$ , $u'/S_L>1$	화염면이 접히고 주머니가 생김
얇은 반응층	$\mathrm{Ka}>1$ , $\mathrm{Da}>1$	소용돌이가 예열층은 뚫지만 반응층은 못 뚫음
완전교반 반응기	$\mathrm{Da}<1$	난류가 화염 구조 자체를 찢음

실무적 의미는 모델 선택이다. $\mathrm{Da} \gg 1$ 인 화염편 영역에서는 화염면 내부 구조가 층류 화염과 같다고 볼 수 있으므로 flamelet/FGM 테이블링이 정당화된다. $\mathrm{Da} \lesssim 1$ 로 내려가면 그 전제가 깨져서 수송 PDF나 상세 화학 기반 모델이 필요해진다(연소 시뮬레이션). MILD/무염 연소처럼 일부러 $\mathrm{Da}$ 를 1 근처로 끌고 가는 최신 연소 방식이 모델링에서 유독 까다로운 것도 같은 이유다.³

6. 큰 Da가 수치해석에 청구하는 것[편집]

$\mathrm{Da}$ 가 크다는 말은 화학 시간 척도가 유동 시간 척도보다 몇 자릿수 짧다는 뜻이고, 그건 곧 강성 방정식이라는 말이다. 자코비안 고유값의 비가 $10^6$ ~ $10^9$ 까지 벌어지면 명시적 룽게-쿠타법의 시간 간격이 가장 빠른 라디칼에 묶여 계산이 사실상 불가능해진다. 그래서 반응성 유동 코드는 예외 없이 다음 세 가지 중 하나 이상을 쓴다.

연산자 분리. 대류·확산과 반응을 한 스텝 안에서 번갈아 푼다(스트랑 분할이면 2차 정확도). 유동 항은 명시적으로, 반응 항은 셀마다 독립 ODE로 암시적으로 — 물리가 다른 두 문제를 각자 잘 맞는 적분기에 맡기는 것이다. 대가는 분리 오차이며, 그 오차가 커지는 지점이 하필 $\mathrm{Da}\sim1$ 인 화염면 근처다.
암시적 강성 적분기. BDF 계열(CVODE/DVODE), Radau5 등. 매 스텝 자코비안 행렬을 만들고 뉴턴 반복을 돌리므로 종 수의 제곱~세제곱으로 비용이 커진다.
테이블링·축소. 애초에 온라인으로 안 푼다. 화염편 테이블(FGM/FPI), ISAT 같은 적응 표작성, QSSA·DRG 기반 메커니즘 축소가 여기 속한다. $\mathrm{Da} \to \infty$ 극한이 참일수록 테이블이 잘 맞는다는 점에서, 테이블링은 큰 $\mathrm{Da}$ 에 대한 대응인 동시에 큰 $\mathrm{Da}$ 를 전제로 하는 방법이다.

격자 설계에도 직접 영향을 준다. $\mathrm{Da}$ 가 크면 반응층이 얇아지므로( $\delta \sim \sqrt{D\tau_{\text{chem}}}$ ), 그 층을 해상하려면 격자 간격이 그에 맞춰 줄어야 한다. 반응 소스항을 격자에 못 담으면 화염 위치와 소염 판정이 통째로 틀어지고, 격자를 줄이면 답이 좋아진다는 CFD의 상식이 배신당하는 몇 안 되는 현장이 된다. 국소적으로만 얇은 층이 생긴다는 성질 자체는 적응 격자 세분화와 궁합이 아주 좋다.

7. 관련 문서[편집]

심지어 역수로 정의하는 문헌도 있다. 어떤 대기화학 논문은 화학 시간을 수송 시간으로 나눈 값을 $\mathrm{Da}$ 라 부르는데, 그러면 “ $\mathrm{Da}$ 가 크다”의 물리적 의미가 정반대가 된다. 무차원수 이름 하나로 논문 두 편의 결론이 충돌하는 광경은 생각보다 자주 벌어진다. ↩
소염점 근처에서 반응기 코드가 “수렴 실패”를 뱉는 것은 코드가 못나서가 아니라, 그 지점에서 실제로 정상 상태가 소멸하기 때문이다. 완화 계수를 아무리 만져도 없는 해는 안 나온다. 수치가 아니라 물리가 거절하고 있다는 걸 알아채는 데 보통 며칠이 걸린다. ↩
“ $\mathrm{Da}$ 를 키우려면 반응을 빠르게 하거나 유동을 느리게 하면 된다”는 말은 맞지만, 연소기 설계에서 유동을 느리게 하는 것은 대개 출력을 줄이는 것과 같다. 그래서 실제 손잡이는 온도 하나뿐이고, 온도는 아레니우스 지수 안에 앉아 있다. 모든 길은 활성화 에너지로 통한다. ↩