연속 교반 탱크 반응기

계산화학 수치해석 마지막 수정: 2026-08-14 04:31:45

1. 개요[편집]

연속 교반 탱크 반응기 Continuous Stirred-Tank Reactor
약칭	CSTR
다른 이름	완전혼합 흐름 반응기 · WSR/PSR(연소 분야)
핵심 가정	완전혼합 — 탱크 안 = 유출 조성
대표 척도	체류시간 $\tau = V/q$
간판 현상	다중 정상상태 · 점화/소화 이력
분기	안장-마디 · 호프
대조군	관형 반응기(PFR)

넣는 순간 이미 다 섞였다고 치자. 이 뻔뻔한 가정 하나가 반응공학 절반을 만들었다.

연속 교반 탱크 반응기(continuous stirred-tank reactor, CSTR)는 원료가 연속으로 들어오고 생성물이 연속으로 나가되, 탱크 내부가 완벽히 균질하게 섞여 있어서 내부의 조성·온도가 곧 유출 스트림의 조성·온도와 같다고 보는 이상화된 반응기다. 공간 좌표가 사라지고 시간만 남으므로, 편미분방정식이 상미분방정식 몇 줄로 붕괴한다. 반응공학 교과서의 첫 번째 반응기이자, 화학 CFD가 셀 하나를 다룰 때 마음속에 그리는 모형이기도 하다.

CSTR이 학부 첫 장에 나오는 이유는 계산이 쉬워서만이 아니다. 이 단순한 모형이 다중 정상상태, 점화-소화 이력, 자기 지속 진동 같은 비선형 동역학의 대표 현상을 전부 품고 있다. 방정식이 두 줄인데 분기 이론 교과서 한 챕터가 나온다. 연소 분야에서는 같은 이상화를 완전혼합 반응기(well-stirred reactor, WSR) 또는 PSR이라 부르며, 이름만 다를 뿐 지배 방정식은 같다.

2. 완전혼합 가정과 물질수지[편집]

검사체적을 탱크 전체로 잡고 화학종 $j$ 에 대해 수지를 쓴다. 유입 유량 $q$ , 부피 $V$ , 유입 농도 $C_{j,\text{in}}$ 라 하면

V\frac{dC_j}{dt} = q\,(C_{j,\text{in}} - C_j) + V\,r_j

여기서 완전혼합 가정이 하는 일은 유출항에 $C_{j,\text{out}}$ 대신 $C_j$ 를 쓸 수 있게 해 주는 것 하나다. 그 한 줄 덕분에 공간 미분이 통째로 사라진다.

$\tau = V/q$ 를 체류시간(space time)이라 부른다. 1차 비가역 반응 $\mathrm{A}\to\mathrm{B}$ , $r_A = -kC_A$ 의 정상상태를 풀면

C_A = \frac{C_{A0}}{1+k\tau}, \qquad X = \frac{k\tau}{1+k\tau} = \frac{\mathrm{Da}}{1+\mathrm{Da}}

즉 전환율이 담쾰러 수 $\mathrm{Da}=k\tau$ 하나로 결정된다. 무차원수 하나로 반응기 설계가 끝나는, 화학공학의 대표적인 미담이다.

한편 체류시간은 어디까지나 평균이라는 점을 잊으면 안 된다. 이상 CSTR의 체류시간분포는 완전한 지수분포다.

E(t) = \frac{1}{\tau}e^{-t/\tau}

평균은 $\tau$ , 표준편차도 $\tau$ . 다시 말해 들어오자마자 나가는 분자와 영원히 머무는 분자가 공존한다. 이 지독한 분산이 뒤에 나올 “CSTR은 PFR보다 전환율이 낮다”의 물리적 정체다. 반대로 PFR의 체류시간분포는 폭이 0인 델타 함수다.

3. 열수지 — 발생열 곡선과 제거열 직선[편집]

발열반응이 들어오면 이야기가 완전히 달라진다. 반응열 $(-\Delta H_r)$ , 밀도 $\rho$ , 비열 $c_p$ , 냉각 자켓의 총괄 열전달 $UA$ , 냉매 온도 $T_c$ 라 하면 에너지 수지는

\rho c_p V \frac{dT}{dt} = \rho c_p q\,(T_{\text{in}}-T) + (-\Delta H_r)\,V k(T)\,C_A - UA\,(T-T_c)

정상상태에서는 물질수지가 $C_A$ 를 온도에 종속시키므로( $C_A = C_{A0}/(1+k(T)\tau)$ ), 열 문제가 온도 하나만의 방정식으로 줄어든다. 이때 우변을 발생열과 제거열로 갈라 놓는 것이 고전적인 도해법이다.

Q_g(T) = (-\Delta H_r)\,q\,C_{A0}\,\frac{k(T)\tau}{1+k(T)\tau}, \qquad Q_r(T) = (\rho c_p q + UA)\,(T - T^{*})

여기서 $T^{*} = (\rho c_p q\,T_{\text{in}} + UA\,T_c)/(\rho c_p q + UA)$ 다. 두 함수의 생김새가 극적으로 다르다.

$Q_g$ 는 S자 곡선이다. 속도상수 $k$ 가 아레니우스 방정식 $k = A e^{-E_a/RT}$ 를 따르므로 저온에서는 거의 0, 중간 온도에서 지수적으로 폭발, 고온에서는 $X\to1$ 로 포화한다. 포화의 원인이 화학이 아니라 반응물이 다 떨어져서라는 점이 중요하다. 아무리 뜨거워도 들어온 것보다 많이 태울 수는 없다.
$Q_r$ 은 직선이다. 기울기 $\rho c_p q + UA$ , 가로축 절편 $T^{*}$ . 유량을 키우거나 냉각을 강화하면 기울기가 서고, 냉매 온도를 올리면 직선이 오른쪽으로 평행이동한다.

정상상태는 이 S자와 직선의 교점이다. 그리고 S자와 직선은 교점이 하나일 수도, 셋일 수도 있다.

4. 다중 정상상태와 점화-소화 이력[편집]

직선이 S자의 완만한 아랫부분만 자르면 저온 교점 하나, 윗부분만 자르면 고온 교점 하나. 그런데 직선의 기울기가 적당하고 위치가 알맞으면 교점이 세 개가 된다. 같은 반응기, 같은 유량, 같은 공급 조건인데 정상상태가 셋이라는 뜻이다.

발열 CSTR 의 정상상태를 온도 y 하나의 스칼라 방정식으로 정리해, 발생열 S자와 제거열 직선(기울기 1+β)의 교점을 이분법으로 풀고 야코비안 2×2 의 대각합·행렬식으로 안정성을 판정한다 — 가운데 교점은 det<0 안장이다. 오른쪽은 냉각수 온도 y_c 를 올렸다 내리는 준정적 스윕을 RK4 로 적분한 것으로, γ=20·Da=0.02·B=11·β=0.4 에서 점화 y_c=4.626, 소화 y_c=2.599, 폭 2.027 의 히스테리시스 고리가 닫힌다. 냉각계수 β 를 0.81 이상으로 올리면 접점이 사라져 고리도 사라진다.

세 해의 성격은 다르다. 아래·위 두 해는 안정하고, 가운데 해는 불안정해서 실제로는 관측되지 않는다. 여기서 나오는 것이 반응공학에서 가장 유명한 그림, 점화-소화 이력 곡선이다.

저온 가지에서 냉매 온도(또는 체류시간)를 서서히 올리면, 어느 지점에서 직선이 S자에 접하고 저온 교점과 중간 교점이 충돌해 소멸한다. 그 순간 반응기는 남은 유일한 안정해 — 고온 가지 — 로 점프한다. 이것이 점화(ignition)다.
고온 가지에서 반대로 냉매 온도를 낮추면, 다른 접점에서 고온 교점이 소멸하며 저온 가지로 떨어진다. 이것이 소화(extinction)다.
두 사건이 서로 다른 조건에서 일어나므로 위·아래 경로가 갈라진다. 한 번 불붙은 반응기는 붙일 때보다 훨씬 나쁜 조건까지 버틴다. 가스터빈의 재점화 고도가 소염 고도보다 낮은 것과 같은 구조다.

접점이 곧 **안장-마디 분기**다. 안정해와 불안정해가 만나서 함께 사라지고, 그 순간 자코비안의 행렬식이 0이 된다. 좌굴 문제의 하중 극한점(스냅스루)과 수학적으로 같은 사건이라는 점은 여러 번 강조할 만하다. 구조물이 “튀는” 것과 반응기가 “붙는” 것이 같은 정규형에서 나온다.¹

실무적 함의는 무섭다. 다중 정상상태 영역에서 운전하는 반응기는 일시적인 교란 하나로 다른 상태로 넘어가서 돌아오지 않는다. 공급 온도가 잠깐 튀거나 교반기가 몇 초 멈춘 것으로 저온 가지에 있던 반응기가 점화되면, 그다음은 열폭주 시나리오다. 반응기 안전 설계가 정상상태 계산이 아니라 분기 도표 계산인 이유다.

5. 안정성, 그리고 호프 분기[편집]

“발생열 곡선이 제거열 직선보다 완만하면 안정하다”는 기울기 조건( $dQ_g/dT < dQ_r/dT$ )이 오래된 경험칙이다. 물리적으로는 자명하다 — 온도가 살짝 올랐을 때 제거열이 발생열보다 더 많이 늘어야 되돌아온다.

문제는 이것이 필요조건일 뿐 충분조건이 아니라는 점이다. CSTR은 $(C_A, T)$ 두 상태변수를 가진 2차원 동역학계이고, 제대로 된 판정은 자코비안 행렬의 대각합과 행렬식으로 한다.

\det J > 0 \ \ \text{그리고} \ \ \operatorname{tr} J < 0 \quad \Longleftrightarrow \quad \text{점근 안정}

기울기 조건은 사실상 $\det J>0$ 을 말한 것이라, $\operatorname{tr} J$ 는 놓친다. $\det J>0$ 을 유지한 채 $\operatorname{tr} J$ 가 0을 가로지르면 **호프 분기**가 일어나고, 정상상태가 불안정해지면서 그 주위에 극한 순환이 태어난다. 즉 반응기가 정상상태에 앉지 않고 스스로 진동한다.

메커니즘은 직관적이다. 온도가 오른다 → 반응이 가속해 반응물을 태워 없앤다 → 연료가 떨어져 발생열이 급감한다 → 냉각이 이겨 온도가 떨어진다 → 반응물이 다시 채워진다 → 재점화. 농도와 온도 두 시간척도가 어긋나 있으면 이 순환이 닫히고, 이완 진동에 가까운 톱니 파형이 나온다. 실험실 반응기에서 온도가 주기적으로 출렁이는 것을 보고 “온도계가 고장 났다”고 결론 내리기 전에 의심해 볼 것.²

6. PFR과의 비교, 그리고 탱크 직렬[편집]

같은 1차 반응, 같은 $\mathrm{Da}=k\tau$ 에서 두 이상 반응기의 전환율은

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

$\mathrm{Da}=2$ 에서 67% 대 86%, $\mathrm{Da}=5$ 에서 83% 대 99.3%. 전환율이 높아질수록 격차가 벌어져서, 99% 전환을 목표로 하면 CSTR은 PFR의 20배 넘는 부피를 요구한다. 이유는 앞서 본 체류시간분포에 있다. CSTR은 언제나 유출 농도에서 반응한다 — 즉 반응기 전체가 가장 묽은 조건에서 돌아간다. PFR은 입구의 진한 농도에서 출발해 점점 묽어지므로 평균 반응속도가 높다.

당연히 예외도 있다. 반응차수가 음이거나 자기촉매 반응이면 “묽은 곳에서 계속 돈다”가 오히려 유리해져 CSTR이 이긴다. 그리고 발열반응의 고온 가지에 눌러앉아 안정적으로 운전하는 능력, 열 제거 면적을 확보하기 쉬운 구조, 슬러리·촉매 현탁액 처리 능력은 CSTR만의 장점이다.

절충안이 탱크 직렬(tanks-in-series)이다. 같은 총 부피를 $N$ 개의 CSTR로 쪼개면

X_N = 1 - \frac{1}{(1+\mathrm{Da}/N)^{N}} \;\xrightarrow[N\to\infty]{}\; 1-e^{-\mathrm{Da}}

$N\to\infty$ 에서 PFR로 수렴한다. 이 관계는 설계뿐 아니라 모형화 도구로 더 자주 쓰인다. 실제 반응기의 체류시간분포를 측정한 뒤 그것을 재현하는 $N$ 을 역산해 축방향 혼합 정도를 요약하는 것이다. $N=1$ 이면 완전혼합, $N$ 이 크면 플러그 흐름에 가깝다.

7. 수치해석에서의 CSTR[편집]

정상상태 곡선 그리기. 온도를 미지수로 뉴턴법을 돌리면 접점(안장-마디)에서 반드시 발산한다. 없는 해로는 수렴하지 않으니 완화계수를 만져도 소용없다. 정답은 매개변수를 미지수로 승격시키는 수치 연속법이며, 구조해석의 호장법과 같은 아이디어다. S자 전체(불안정 중간 가지 포함)를 완주하려면 이것 말고는 없다.
과도 응답 적분. 상세 반응 기구를 붙이면 라디칼의 시간척도가 유동 시간척도보다 몇 자릿수 짧아 전형적인 강성 방정식이 된다. 명시적 룽게-쿠타법은 사실상 불가능하고 BDF 계열 암시적 적분기가 국룰이다.
완전혼합 가정의 검증. 실제 교반 탱크는 데드존·단락류·불완전 혼합을 갖는다. 추적자 실험으로 체류시간분포를 재거나 전산유체역학으로 유동장을 풀어 확인하고, 어긋나면 이상 반응기 몇 개를 이어 붙인 반응기 네트워크 모형(compartment model)으로 보정한다. 반응 시간척도가 난류 혼합 시간척도보다 짧으면 거시 혼합이 아무리 좋아도 분자 수준 혼합(micromixing)이 율속이 되어, 완전혼합 가정 자체가 틀린다.
CFD 셀과의 관계. 반응성 유동 해석에서 연산자 분리로 반응 항만 떼어 셀마다 독립 ODE로 푸는 순간, 각 셀은 사실상 작은 0차원 반응기가 된다. 그래서 CSTR/PSR 해석에서 얻은 강성 대응 노하우가 그대로 연소 시뮬레이션 코드에 재활용된다.³

8. 관련 문서[편집]

그래서 구조해석 하던 사람이 반응공학 세미나에 들어가면 슬라이드의 S자 곡선을 보고 “저거 우리 하중-변위 곡선인데요”라고 중얼거리게 된다. 실제로 같은 것이 맞고, 두 분야가 각자 이름을 붙였을 뿐이다. 극한점, 접힘점, 점화점, 스냅 하중 — 전부 한 식구. ↩
진동하는 반응기를 처음 본 대학원생의 반응은 대체로 (1) 열전대 접촉 불량 의심, (2) 교반기 RPM 의심, (3) 유량 펌프 맥동 의심 순이고, 호프 분기를 의심하는 것은 대개 네 번째다. 참고로 진짜 맥동 펌프인 경우도 많아서 (3)까지는 합리적이다. ↩
반대 방향의 오해도 흔하다. “CSTR은 0차원이니까 싸다”고 생각했다가, 화학종 500개짜리 상세 기구를 넣는 순간 자코비안이 500×500이 되고 한 점 계산에 수 분이 걸리는 광경을 보게 된다. 차원이 낮다고 싼 게 아니라 미지수가 적어야 싼 것이다. ↩