Damkohler 수와 반응 침투

반응을 빠르게 하면 재료 내부까지 잘 사용하는가

다공성 촉매 표면에 반응물이 공급된다고 생각해보자. 분자는 내부로 확산하면서 반응으로 소모된다. 반응 속도를 키우면 전체 촉매를 더 잘 사용할 것 같지만, 너무 빠른 반응은 표면 근처에서 대부분의 분자를 소모하여 내부 공급을 부족하게 만들 수 있다. 국소 반응이 빨라진 배율만큼 전체 반응량이 늘어나는 것은 아니다.

비교할 것은 확산이 재료 내부를 보충하는 시간과 반응이 물질을 소모하는 시간이다. 이 두 시간의 비를 직접 유도하고 유한한 평판의 정확한 해를 구해보자. 이어서 실제 저장 질량과 소모량으로 수지를 점검한다. 과학 AI의 입력 변수, 매개변수 추정, 닮은 실험 사이의 전이를 생각할 수 있는 작은 모델이 된다.

기하와 물리 가정을 먼저 정한다

x=0은 외부에 노출된 표면이고 x=L은 밀폐면 또는 대칭면이다. 양면에서 같은 조건으로 공급되는 대칭 판이라면 L은 전체 두께의 절반이다. A는 이 반쪽 판의 실제 단면적이다. 표면 농도 c0를 지정하고 안쪽 면에서는 법선 방향 확산 플럭스를 0으로 둔다.

확산계수 D는 양의 상수다. 체적당 일차 소모율은 kc이며 k≥0인 상수다. 유체 이류·온도 변화·외부 전달 저항·비선형 포화·공극 구조 변화는 제외한다. 농도는 희석된 비음수 scalar이며 여러 화학종을 모두 추적하는 모델은 아니다.

$$Dc_{xx}-kc=0,\qquad c(0)=c_0,\quad c_x(L)=0,\qquad [D]=\mathrm{m^2/s},\quad[k]=\mathrm{s^{-1}}.$$

농도 단위는 kg/m³, 국소 소모 밀도는 kg/(m³ s), 체적 적분한 전체 소모율은 kg/s다. 무유속 경계는 농도 0을 뜻하지 않는다. 농도 그림이나 애니메이션을 볼 때도 기울기·농도·면 플럭스·체적 소모량을 구별해야 한다.

반응이 0이어도 유효한 척도로 무차원화한다

ξ=x/L로 놓고 확산 시간 L²/D를 사용한다. c0가 양수라면 농도를 표면 농도로 나눈다. 연쇄법칙을 적용하고 공통 계수를 나누면 정상 식에는 반응과 확산을 비교하는 수 하나가 남는다.

$$\xi=\frac{x}{L},\qquad C=\frac{c}{c_0},\qquad C_{\xi\xi}-\mathrm{Da}\,C=0,\qquad \mathrm{Da}=\frac{kL^2}{D}=\phi^2.$$

여기의 Da는 확산 시간을 사용하는 Damköhler 수이며 φ는 선언한 일차 반응·상수 물성·평판 관례의 Thiele 계수다. 다른 문제는 이류 시간이나 다른 반응식을 사용하여 Damköhler 수를 정의할 수 있다. 이름이 같다는 이유로 숫자만 비교하지 말고 정의와 척도를 확인해야 한다.

c0=0이면 물리 해는 영 농도장이다. 그러나 측정 비 c/c0는 미정이다. 표면 농도를 1로 선언한 기저 함수는 계산할 수 있지만, 이를 영 농도장에서 측정한 정규화 값이라고 표시하지 않는다.

시간과 침투 길이로 직관을 얻는다

양의 일차 반응 상수의 역수는 반응 시간이다. 이를 확산 시간과 비교하면 Da를 얻는다. 한 반응 시간 동안 확산이 도달하는 대표 거리는 침투 길이를 준다.

$$t_D=\frac{L^2}{D},\qquad t_R=\frac1k\quad(k>0),\qquad \mathrm{Da}=\frac{t_D}{t_R},\qquad \ell_p=\sqrt{\frac Dk},\qquad \phi=\frac L{\ell_p}.$$

Da가 작으면 확산이 소모보다 빠르게 공급을 보충하여 농도가 거의 균일하다. Da가 크면 판 전체를 보충하기 전에 소모가 진행되어 반응이 표면 가까이에 집중된다. 침투 길이는 대표 척도이지 그 깊이에서 농도가 갑자기 0이 되는 날카로운 경계가 아니다.

k=0이면 유한한 반응 시간이나 반응 침투 길이를 만들지 않는다. 물리 해와 확산 시간은 여전히 유효하다. 임의의 큰 유한값을 넣으면 반응 기작이 없다는 상태를 숨기게 된다.

유한 판의 정확한 해를 유도한다

특성방정식의 두 근은 ±√(k/D)다. 두 지수 함수를 밀폐면을 중심으로 하는 cosh로 묶으면 그 면의 미분 조건이 바로 만족된다. 남은 계수는 표면 농도로 정한다.

$$c(x)=c_0\frac{\cosh[\phi(1-\xi)]}{\cosh\phi},\qquad c_x=-\frac{c_0\phi}{L}\frac{\sinh[\phi(1-\xi)]}{\cosh\phi}.$$

표면에서는 cosh의 비가 1이다. 안쪽 면에서는 기울기가 0이지만 농도는 c0/coshφ이므로 양의 c0와 유한한 φ에서 양수다. 내부로 갈수록 농도는 감소하며 반응과 농도가 양수인 곳의 곡률은 양수다. 이는 D cxx=kc에 필요한 모양이지 내부에서 반응물을 생성한다는 뜻이 아니다.

k가 정확히 0이면 원래 식과 같은 경계조건에서 농도가 일정하고 기울기와 소모가 0인 해를 얻는다. 쌍곡함수 해의 극한도 이 결과와 일치한다.

실제 질량 수지를 복원한다

양의 x 방향 확산 플럭스는 −Dcx다. 이를 미분하면 국소 소모의 음수가 된다. 실제 길이 구간에서 적분하고 실제 단면적을 곱한다.

$$j=-Dc_x,\qquad j_x=-kc,\qquad A j(0)=A\int_0^L kc(x)\,dx.$$

안쪽 면의 플럭스가 0이므로 유입된 물질은 그 면을 통해 빠져나가기 전에 모두 소모된다. 정상 상태의 유입률과 소모율은 같지만 내부에는 일정한 비영 질량이 저장될 수 있다. 반응이 없던 이류·확산 단원과 달리 여기의 플럭스는 공간적으로 일정하지 않다.

실제 공간 적분에서 저장 질량과 전체 소모율을 구하면 다음과 같다. 두 수의 숫자가 같을 때도 차원은 다르다.

$$M=A\int_0^L c(x)\,dx =Ac_0L\frac{\tanh\phi}{\phi},\qquad R=kM=Ac_0\sqrt{Dk}\,\tanh\phi\quad(k>0).$$

양면에 같은 조건이 주어진 대칭 판 전체의 저장량과 소모율은 이 반쪽 결과의 2배다. φ 안의 L을 전체 두께로 바꾸면서 적분도 다시 2배 하면 기하를 중복하여 잘못 계산한다.

효과도와 침투 길이는 다른 관측이다

판 내부 전체가 표면 농도로 유지된다는 이상적인 기준의 소모율은 kc0AL이다. 이 기준이 비영일 때 실제 소모율을 나누면 효과도를 얻는다.

$$\eta_{\rm eff}=\frac{R}{kc_0AL} =\frac{\tanh\phi}{\phi},\qquad \lim_{\phi\to0}\eta_{\rm eff}=1.$$

효과도는 반응률의 비인 무차원 값이고 침투 길이는 길이 차원을 가진다. φ가 크면 효과도는 대략 1/φ다. 이 모델에서 k를 키우면 전체 소모율은 계속 증가하지만 큰 k에서는 k에 비례하기보다 √k에 비례하는 선도 거동을 보인다.

k=0 또는 c0=0이면 실제 반응률의 비는 0/0이다. 기저 효과도 함수에는 φ=0에서 1이라는 유용한 연속 연장이 있지만, 이를 실제 측정 반응률 비와 구분해 표시해야 한다. 여기서는 A를 양수로 제한하므로 면적 0에 따른 별도의 분모 예외는 없다.

실제 단위로 계산해본다

L=1m, D=0.25m²/s, k=1/s, c0=2kg/m³, A=1.5m²를 택하자. 확산 시간은 4s, 반응 시간은 1s, Da=4, φ=2다. 침투 길이는 0.5m다.

밀폐면 농도는 약 0.53160446kg/m³, 효과도는 약 0.48201379다. 표면 플럭스는 약 0.96402758kg/(m² s)이므로 전체 소모율은 약 1.44604137kg/s다. 저장 질량은 약 1.44604137kg이다. 마지막 두 숫자가 같은 것은 k=1/s라는 선택 때문이며 같은 물리량이라는 뜻은 아니다.

이상적인 균일 공급 기준의 소모율은 3kg/s다. 실제 소모율을 이 값으로 나누면 같은 효과도를 얻는다. 안쪽의 비영 농도도, 시간에 따라 일정한 저장량도 무유속 경계나 정상 질량 수지를 위반하지 않는다.

균일하게 반올림되어도 고갈은 남을 수 있다

Da가 작을 때 1에서 거의 같은 cosh 비를 직접 빼면 실제 고갈을 잃는다. 쌍곡함수 곱의 항등식을 사용하면 안정적으로 계산할 수 있다.

$$1-C= \frac{2\sinh[\phi(1-\xi/2)]\sinh(\phi\xi/2)}{\cosh\phi} =\frac{\mathrm{Da}}2\xi(2-\xi)+O(\mathrm{Da}^2).$$

안쪽 면의 선도 고갈은 Da/2, 중앙점은 3Da/8, 표면은 정확히 0이다. 농도 자체가 c0로 반올림되어도 고갈을 따로 계산하면 양의 값을 보존할 수 있다. 효과도도 1로 반올림되는 동안 양의 효과도 부족분은 남을 수 있으며 선도항은 Da/3이다.

kL²/D를 일찍 일반 숫자로 곱하면 작은 양의 값이 언더플로한다. 로그에서 계산하면 비영 그룹을 보존하고 Da가 일반 숫자로 표현되지 않아도 φ는 표현할 수 있다. 실제 0, 일반 숫자, 로그로만 남는 비영 상태를 구별해야 한다. 표현 불가를 반응이 없다는 뜻으로 해석하면 안 된다.

큰 φ에서 밀폐면의 꼬리를 보존한다

직접 계산한 cosh 두 개가 모두 넘칠 수 있으므로 큰 φ에서는 로그 비를 사용한다. 감쇠와 밀폐면의 반사 기여를 분리하면 다음 표현을 얻는다.

$$\log C=-\phi\xi+ \log\!\left(1+e^{-2\phi(1-\xi)}\right) -\log\!\left(1+e^{-2\phi}\right),\qquad C(1)\sim2e^{-\phi}.$$

밀폐면의 계수 2가 중요하다. 반무한 영역의 단일 지수는 다른 먼 경계조건을 가지므로 유한 판 끝의 정확한 농도를 대신할 수 없다. 아주 작은 양의 농도가 일반 숫자로 표시되지 않아도 로그는 유한할 수 있다. 이 사실은 벽의 기울기가 정확히 0이라는 조건을 바꾸지 않는다.

닮은 판과 다른 무차원 수를 연결한다

양의 a,b에 대해 L′=aL, D′=bD, k′=(b/a²)k로 바꾸면 Da·φ와 같은 ξ에서의 정규화 곡선이 보존된다. 확산 시간과 반응 시간은 모두 a²/b배다. 표면 농도를 유지하고 면적도 기하학적으로 닮게 바꾸면 저장 질량은 a³배, 전체 소모율은 ab배다.

외부 전달 저항, 온도 의존 반응률, 다른 기하나 다른 반응 차수가 중요해지면 확산 기반 Da 하나를 맞추는 것으로 충분하지 않다. 추가 매개변수와 방정식이 필요할 수 있다. AI 학습에서 Da에 따라 자료를 겹쳐 보는 것은 선언한 가족 안에서 유용하지만 가족 밖의 기작에 대한 증거까지 지워서는 안 된다.

앞선 점탄성과 연결하면 Weissenberg 수는 완화 시간과 변형률 속도 척도의 곱으로 정의할 수 있으며 대표 속도 구배 U/L을 사용할 수 있다. Deborah 수는 완화 시간과 선언한 관측 시간 T의 비다. T가 주기라면 ωτ와 τ/T에는 2π 차이가 있다. Mach·Froude 수는 또 다른 항을 비교한다. 각각의 기준 척도와 가정을 밝혀야 하며, 두 scalar 예제가 이 모든 계를 시뮬레이션하는 것은 아니다.

풀이가 있는 여섯 문제

문제 1 — 전체 두께

전체 두께가 2m이고 양면에 같은 농도가 공급되는 대칭 판에서 이 단원의 L은 얼마인가? 전체 소모율은 어떻게 구하는가?

풀이: 한쪽 표면에서 대칭면까지인 L=1m를 사용한다. 반쪽 판의 소모율을 구해 2배 한다. Da에 2m를 그대로 넣으면 선언한 관례보다 4배 큰 값을 얻는다.

문제 2 — 밀폐 경계

안쪽 면의 플럭스가 0이면 그 면의 농도도 0인가?

풀이: 아니다. 양의 D에서 무유속은 기울기 0을 뜻한다. 농도는 c0/coshφ이며 유한한 φ와 양의 c0에서 양수다.

문제 3 — 반응 제거

양의 표면 농도를 유지하고 k=0으로 바꾸면 저장량·소모율·효과도는 어떻게 되는가?

풀이: 농도가 균일하여 저장량은 Ac0L이고 소모율은 정확히 0이다. 실제·이상 소모가 모두 0이므로 측정 효과도는 미정이다. 해석적 연장은 1이며 유한한 반응 시간은 없다.

문제 4 — 작은 고갈

Da=10⁻²⁰에서 중앙의 정규화 고갈과 효과도 부족분을 추정하라.

풀이: 선도 고갈은 3.75×10⁻²¹, 효과도 부족분은 약 3.3333×10⁻²¹이다. C와 효과도가 1로 반올림되어도 이 보정이 정확히 0인 것은 아니다.

문제 5 — 빨라진 반응

이미 φ가 매우 클 때 다른 입력을 유지하고 k를 4배 하면 대략 어떻게 되는가?

풀이: φ는 2배, 침투 길이는 절반, 효과도는 대략 절반, 전체 소모율은 대략 2배다. 확산이 제한하는 상태에서 국소 반응 속도 4배가 전체 소모 4배를 뜻하지 않는다.

문제 6 — 닮음과 질량

a=2, b=1로 놓고 면적도 닮게 바꾸면 Da가 같은 동안 무엇이 달라지는가?

풀이: L은 2배, k는 1/4배, 면적은 4배다. 두 시간은 4배, 저장량은 8배, 소모율은 2배다. 곡선이 겹친다는 이유로 실제 재고와 반응률이 같아지는 것은 아니다.

자료와 적용 한계

MIT 3.185 Transport Phenomena의 Chapter 2는 대칭 판의 정상 일차 반응·확산 해를 유도한다. 원문의 전체 두께 관례를 여기의 반두께로 변환해야 한다. Finlayson의 저자 제공 반응확산 교재는 Thiele 시간 비와 기하·경계 전달의 역할을 설명한다. 이 단원은 선언한 반쪽 판의 수지와 닮음을 직접 유도했다. 외부 전달 저항·반응열·비선형 다성분 반응은 모델 범위 밖이다.