자착화

계산화학 전산유체역학 시뮬레이션 마지막 수정: 2026-08-24 04:46:18

1. 개요[편집]

자착화 Autoignition
정의	외부 점화원 없이 혼합기가 스스로 발화
핵심 지표	착화지연시간 τ_ig
표준 계산	0차원 단열 반응기 + 강성 ODE 적분
표준 실측	충격파관 · 급속압축기
특이 현상	음의 온도계수(NTC) · 냉염 · 2단 착화
응용	디젤 · HCCI · 노킹 · 화재 안전

불을 붙이지 않았는데 붙었다. 디젤 엔진에서는 설계이고 가솔린 엔진에서는 사고다.

자착화(autoignition, 자발 착화)는 불꽃·열선 같은 외부 점화원 없이, 연료-산화제 혼합기가 스스로의 발열 화학반응이 만든 열과 라디칼의 되먹임만으로 발화하는 현상이다. 압축이나 가열로 혼합기가 어떤 온도·압력에 도달하면, 처음에는 눈에 띄지 않을 만큼 느리게 진행되던 반응이 지수적으로 가속해 유한한 시간 뒤 폭발적으로 완결된다.

이 “유한한 시간”이 자착화 연구의 거의 전부다. 착화지연시간(ignition delay time) $\tau_{ig}$ 라 부르며, 디젤 엔진의 연소 위상, 가솔린 엔진의 노킹 한계, HCCI의 제어 가능성, 가스터빈 예혼합기의 안전 여유, 화학공정의 폭발 위험도가 전부 이 한 값으로 환원된다.

담쾰러 수 문서에 나온 점화-소염 S자 곡선의 아래쪽 접힘점이 바로 자착화이고, 그 접힘은 비유가 아니라 안장-마디 분기다. 여기서는 그 사건을 시간 영역에서 어떻게 정의하고, 계산하고, 측정하고, 검증하는지를 다룬다.

2. 착화지연시간을 어떻게 정의할 것인가[편집]

의외로 여기서부터 합의가 없다. 문헌에 쓰이는 정의만 나열해도 이렇다.

압력 기준 — $dp/dt$ 가 최대인 시각. 정적 반응기와 급속압축기의 표준.
온도 기준 — $dT/dt$ 가 최대인 시각. 계산에서 가장 흔하다.
OH* 기준 — 여기 상태 OH 라디칼의 화학발광이 최대인 시각. 충격파관 광학 계측의 표준.
외삽 기준 — 신호의 최대 기울기 접선을 그어 기저선과 만나는 점. 신호가 완만할 때 재현성이 좋다.
임계값 기준 — 압력 상승이 초기값의 몇 %를 넘는 시각.

고온·단발 착화에서는 이 정의들이 수 마이크로초 안에서 일치하므로 아무거나 써도 된다. 문제는 저온 2단 착화 영역이다. 1단 착화(냉염)에서 이미 압력이 꿈틀하고 CH₂O 발광이 뜨는데, “압력 상승 임계값” 정의를 쓰면 그 꿈틀거림을 착화로 잡아 버려 $\tau_{ig}$ 가 절반 이하로 튄다. 같은 연료에 대해 실험실마다 배로 다른 값을 보고하는 사태의 상당수가 여기서 나온다. 논문에서 $\tau_{ig}$ 를 보면 정의부터 찾는 것이 국룰이고, 이건 담쾰러 수의 정의를 먼저 찾는 습관과 똑같은 이유다.¹

3. 표준 계산 절차 — 0차원 반응기[편집]

계산 쪽 정의는 훨씬 깔끔하다. 공간을 완전히 지우고 균질한 단열 반응기를 상정한다. 정적(constant volume) 반응기라면

\frac{dY_k}{dt} = \frac{\dot\omega_k W_k}{\rho}, \qquad \frac{dT}{dt} = -\frac{1}{\rho\, c_v}\sum_k u_k\,\dot\omega_k W_k

정압(constant pressure) 반응기라면 $c_v \to c_p$ , $u_k \to h_k$ 로 바꾼다. $\dot\omega_k$ 는 아레니우스 방정식 기반 상세 메커니즘이 주는 생성률이다(반응속도론). 초기 조성·온도·압력을 넣고 적분하면 온도가 완만하다가 어느 순간 수직으로 솟는 곡선이 나오고, 그 무릎이 $\tau_{ig}$ 다.

여기서 유동은 하나도 없다. 그래서 “0차원 반응기 계산”이라 부르고, 유일한 난점은 순수하게 수치적이다.

강성. 라디칼(H, O, OH)의 특성 시간은 나노초, 유도기간은 밀리초다. 자코비안 고유값의 비가 $10^6$ $10^9$ 을 넘고, 이건 정의 그대로 강성 방정식이다. 명시적 룽게-쿠타법을 쓰면 안정성이 가장 빠른 라디칼에 묶여 밀리초 하나를 적분하는 데 스텝 수십억 개가 필요하다. 그래서 실무는 예외 없이 후진 미분 공식(BDF) 계열 암시적 적분기를 쓴다 — SUNDIALS의 CVODE, DASSL, Radau5. 매 스텝 자코비안 행렬을 만들고 뉴턴-랩슨법으로 푸는 구조라 종 수 $K$ 에 대해 비용이 $K^2$ $K^3$ 으로 붙는다.

정확도 설정의 함정. 착화는 지수적 성장이라 유도기간의 미세한 오차가 $\tau_{ig}$ 에 그대로 곱해진다. 상대허용오차를 $10^{-6}$ 정도로 느슨하게 잡고 돌리면 $\tau_{ig}$ 가 수 % 흔들리는 일이 흔하고, 메커니즘 비교 논문에서는 보통 $10^{-8}$ 이하로 조인다. 해석적 자코비안을 쓰는지 수치 미분으로 때우는지도 여기서 갈린다.

Cantera나 Chemkin 같은 도구가 이 절차를 한 줄로 감싸 놓았기 때문에 오늘날 $\tau_{ig}$ 계산 자체는 노트북에서 초 단위로 끝난다. 진짜 비용은 이 계산을 3차원 CFD의 셀 100만 개에서 스텝마다 반복해야 할 때 발생하며, 그때 연산자 분리와 테이블링이 등장한다(연소 시뮬레이션).

참고로 유입·유출이 있는 완전혼합 반응기(WSR/PSR)는 같은 화학에 체류 시간을 붙인 이상화이고, 그 정상 상태 해석은 연속 교반 탱크 반응기 쪽에서 다룬다.

4. 아레니우스 스케일링[편집]

계산 결과를 압축해 공학적으로 쓰려면 상관식이 필요하다. 널리 쓰이는 꼴은

\tau_{ig} = A\,\phi^{a}\,p^{-n}\,\exp\!\left(\frac{E_a}{R T}\right)

이다. $\phi$ 는 당량비, $E_a$ 는 총괄 활성화 에너지(개별 소반응의 것이 아니라 전체 착화 과정을 하나의 아레니우스로 근사했을 때의 유효값)다. 지수 앞의 부호에 주의 — 반응 속도가 $\exp(-E_a/RT)$ 인 것과 반대로, 그 역수인 지연시간은 $\exp(+E_a/RT)$ 다. 온도가 오르면 지연이 짧아진다.

압력 지수 $n$ 은 대개 1 근처다. 반응이 대체로 2차이므로 농도가 압력에 비례하고, 착화 시간이 농도에 반비례하는 것이 1차 근사다. 이 스케일 때문에 압축비를 높이면 착화가 빨라지고, 그래서 디젤 엔진이 성립하며, 그래서 가솔린 엔진의 말단 가스가 노킹한다.

$\ln \tau_{ig}$ 를 $1/T$ 에 대해 그린 그림(착화 지연 아레니우스 선도)이 이 분야의 표준 도표다. 고온에서는 기울기 $E_a/R$ 의 직선이 잘 나오는데, 온도를 낮추면 이 직선이 반드시 배신한다.

5. NTC와 냉염 — 온도를 올렸는데 착화가 늦어진다[편집]

큰 탄화수소(n-헵탄, iso-옥탄, 디젤 유사 연료)의 아레니우스 선도에는 기울기가 뒤집히는 구간이 있다. 대략 650~900 K, 수십 bar 영역에서 온도를 올렸는데 $\tau_{ig}$ 가 오히려 길어진다. 이것이 음의 온도계수(Negative Temperature Coefficient, NTC) 영역이다.

원인은 저온 산화 경로의 경쟁이다.

연료 $\mathrm{RH}$ 가 라디칼에 수소를 빼앗겨 알킬 라디칼 $\mathrm{R}$ 이 된다.
$\mathrm{R} + \mathrm{O_2} \rightleftharpoons \mathrm{RO_2}$ — 산소 첨가. 가역이며, 온도가 오르면 평형이 왼쪽으로 밀린다.
$\mathrm{RO_2}$ 가 분자 내 수소 이동으로 $\mathrm{QOOH}$ 가 되고, 여기에 산소가 한 번 더 붙어 $\mathrm{O_2QOOH}$ → 케토하이드로퍼옥사이드로 간다.
케토하이드로퍼옥사이드가 깨지면 라디칼이 두 개 이상 나온다. 즉 연쇄분지. 이게 저온 착화의 엔진이다.

온도가 올라가면 2단계의 평형이 왼쪽으로 밀려 $\mathrm{RO_2}$ 가 살아남지 못하고, 대신 $\mathrm{R} + \mathrm{O_2} \to$ 알켄 $+\ \mathrm{HO_2}$ 같은 연쇄전파(라디칼 수가 안 늘어나는) 경로가 이긴다. 분지가 꺼지므로 전체 반응이 느려진다. 더 온도를 올리면 $\mathrm{H_2O_2} \to 2\,\mathrm{OH}$ 와 $\mathrm{H} + \mathrm{O_2} \to \mathrm{O} + \mathrm{OH}$ 라는 고온 분지 경로가 켜지면서 다시 정상적인 아레니우스 거동으로 복귀한다. NTC는 두 분지 메커니즘 사이의 골짜기인 셈이다.

이 저온 분지가 실제로 만드는 관측 가능한 사건이 냉염(cool flame)이다. 온도가 50~200 K 정도만 오르고 멈추는 1단 착화로, 포름알데히드(CH₂O) 여기종이 내는 희미한 푸른빛 때문에 그런 이름이 붙었다. 온도가 오르면 NTC 때문에 반응이 스스로 감속하므로 폭주가 중단되고, 계는 잠시 뜸을 들이다가 고온 화학이 켜지면서 2단 착화로 완결된다. 그래서 저온 자착화의 압력 이력에는 계단이 두 개 찍힌다.

NTC는 학술적 진기명기가 아니라 엔진 제어의 물리적 근거다. HCCI가 부하 영역을 넓히기 어려운 것도, 디젤 연소가 분사 시기에 예민한 것도, 옥탄가라는 지표가 단일 숫자로 안 끝나고 RON/MON 두 개로 갈리는 것도 전부 이 골짜기의 존재 때문이다.²

6. 열폭주 판정 — 세미노프와 프랑크-카메네츠키[편집]

$\tau_{ig}$ 를 상세 화학으로 적분하기 이전에, “애초에 폭주가 일어나는가”를 판정하는 고전 이론이 둘 있다. 둘 다 발열과 방열의 경쟁이라는 같은 그림에서 출발하지만 온도의 공간 분포를 어떻게 다루느냐가 다르다.

세미노프(Semenov) 이론 — 내부 온도 균일. 비오 수가 작아 계 내부 온도가 하나라고 보고, 발열률 $q V A e^{-E/RT}$ 와 뉴턴 냉각 $\chi S (T-T_a)$ 를 비교한다. 발열은 $T$ 의 지수 함수, 방열은 선형 함수이므로 두 곡선은 교차하지 않거나, 두 점에서 만나거나, 접한다. 접하는 순간이 임계이고, 그때의 온도 상승은

\Delta T_{cr} = T_{cr} - T_a = \frac{R T_a^2}{E_a}

라는 놀랍도록 단순한 결과가 나온다. 상온 근처·전형적 활성화 에너지에서 이 값은 수십 K 수준이다 — 겨우 그만큼 데워지고 나면 되돌릴 수 없다는 뜻이며, 화재·저장 안전 설계의 출발점이다.

프랑크-카메네츠키(Frank-Kamenetskii) 이론 — 내부 전도 지배. 반대로 열이 전도로만 빠져나가는 고체·정체 유체라면 온도 분포를 풀어야 한다. 무차원 매개변수

\delta = \frac{Q E_a \rho A L^2}{\lambda R T_a^2}\exp\!\left(-\frac{E_a}{R T_a}\right)

가 임계값을 넘으면 정상 온도 분포가 존재하지 않는다. 임계값은 형상으로 결정되어 무한 평판(반두께 $L$ ) $\delta_{cr}=0.878$ , 무한 원기둥 $2.00$ , 구 $3.32$ 다. $\delta \propto L^2$ 이므로 같은 물질이라도 크게 쌓으면 터진다. 퇴비 더미, 석탄 야적장, 도료 걸레 뭉치, 리튬 배터리 팩의 자연발화가 전부 이 부등식의 응용이며, 위험물 저장 규정의 “적재 높이 제한”이 여기서 나온다.

두 이론이 주는 것은 $\tau_{ig}$ 가 아니라 폭주 여부의 임계 조건이다. 상세 화학 계산이 값싸진 오늘날에도 이 판정식이 살아 있는 이유는, 매개변수 세 개로 안전 여유를 손으로 계산할 수 있기 때문이다.

7. 어떻게 측정하고 무엇을 검증하나[편집]

$\tau_{ig}$ 는 화학반응 메커니즘 검증의 1순위 표적이다. 층류 화염 속도가 전파 화학을, 화학종 농도 프로파일이 경로를 검증한다면, $\tau_{ig}$ 는 분지 화학 전체를 하나의 스칼라로 압축해 시험한다. 주력 장비는 둘이다.

장비	온도 범위(대략)	시험 시간	주된 계통 오차
충격파관	1000 K 이상	수 ms	비이상 압력 상승, 저온에서의 불균일 착화
급속압축기(RCM)	600~1100 K	~100 ms	압축 후 열손실, 와류에 의한 혼합

충격파관은 반사 충격파로 시험 기체를 순간 가열해 균일한 고온·고압 상태를 만든다. 문제는 경계층 성장 때문에 시험부 압력이 가만히 있지 않고 $dp/dt$ 가 수 %/ms로 슬금슬금 오른다는 것. 지연이 짧을 때는 무시할 수 있지만 밀리초를 넘어가면 이 압력 이력을 모형에 넣어 줘야(구속 체적 대신 지정 체적 이력으로) 계산과 실험이 만난다. 1000 K 아래로 내려가면 벽면·불순물이 만드는 국소 착화 커널이 먼저 터지는 비균질 착화가 나타나 0차원 가정 자체가 깨진다.

급속압축기는 피스톤으로 수십 ms 만에 압축해 저온·고압 상태를 오래 유지한다. NTC 영역을 직접 볼 수 있는 거의 유일한 장비. 대신 압축이 끝난 뒤에도 벽으로 열이 계속 빠지므로, 같은 조건에서 비반응 기체(연료를 질소로 치환)로 한 번 더 실험해 압력 이력을 얻고, 그로부터 “유효 체적 이력”을 역산해 0차원 계산에 먹인다. 이 절차를 안 거치면 계산이 실험보다 항상 빨리 착화한다. 피스톤 머리에 크레비스(crevice)를 파 코너 와류를 가두는 설계도 균일성 확보를 위한 표준이 되었다.

이렇게 얻은 데이터로 검증하는 화학반응 메커니즘은 큰 연료일수록 거대해진다. n-헵탄 상세 메커니즘이 종 수백 개·반응 수천 개, iso-옥탄은 그 이상이다. 실제 CFD에 넣으려면 축소가 불가피하고, 축소 과정에서 무엇을 보존할 것인가의 답이 대개 “ $\tau_{ig}$ 를 관심 영역에서 몇 % 이내로”다. 민감도 해석으로 $\tau_{ig}$ 에 영향이 큰 소반응을 골라내는 작업이 그래서 이 분야의 일상이다.

8. 엔진과 난류 속에서[편집]

디젤. 분사부터 착화까지의 지연은 두 가지가 직렬로 붙어 있다. 물리적 지연(분무 미립화·증발·공기와의 혼합)과 화학적 지연( $\tau_{ig}$ ). 세탄가가 높다는 것은 화학적 지연이 짧다는 뜻이고, 그래야 예혼합 연소 비율이 줄어 초기 급격 연소와 소음이 억제된다. 자착화가 설계 사양인 유일한 엔진.

노킹. 가솔린 엔진에서는 화염면이 도달하기 전에 말단 가스(end gas)가 압축·가열되어 스스로 터진다. 여러 지점에서 거의 동시에 자착화하면 급격한 압력파가 발생해 실린더가 울린다. 실무에서 널리 쓰이는 예측이 리븐굿-우 적분이다.

\int_0^{t_{\text{knock}}} \frac{dt}{\tau_{ig}\big(T(t),\,p(t)\big)} = 1

매 순간의 조건에서 계산한 $\tau_{ig}$ 의 역수를 누적해 1이 되는 시점을 노킹 개시로 본다. 유도기간을 “쌓이는 양”으로 취급하는 근사이며, 놀랍도록 잘 맞는다.

HCCI. 예혼합기를 통째로 압축해 자착화시키는 방식. 화염 전파가 없으므로 국소 고온 영역이 없고 NOx가 극적으로 낮으며 효율도 높다. 문제는 착화 시점을 직접 조종할 손잡이가 없다는 것 — 점화 플러그도 분사 시기도 아닌, 오직 화학이 결정한다. 흡기 온도·EGR·압축비로 간접 제어할 수밖에 없고, 부하가 바뀌면 그 손잡이들이 동시에 흔들린다. 자착화 화학을 실시간으로 예측하지 못하면 상용화가 안 되는 구조이며, 이 분야가 $\tau_{ig}$ 정밀도에 가장 목마른 곳이다.

난류-화학 상호작용. 실제 연소기에서 자착화는 균질하게 일어나지 않는다. 온도·조성 요동 때문에 $\tau_{ig}$ 가 공간적으로 분포하고, 가장 짧은 지점(“착화 커널”)에서 먼저 터진 뒤 그 사건이 이웃으로 번진다. 번지는 속도는 화염 전파 속도가 아니라 반응성 구배로 결정되는 자발 착화 전파 속도

u_{sp} = \left|\frac{\partial \tau_{ig}}{\partial x}\right|^{-1}

이고(젤도비치의 구배 메커니즘), 이 속도가 음속에 가까워지면 압력파와 반응이 공명해 데토네이션으로 발전할 수 있다. 초강력 노킹(superknock)의 유력한 시나리오이며, 반대로 $u_{sp}$ 가 충분히 느리면 얌전한 순차 착화로 끝난다.

이 경쟁을 무차원으로 재는 것이 결국 담쾰러 수다. 난류가 화학보다 빠르면( $\mathrm{Da}\ll1$ ) 요동이 평준화되어 균질 착화에 가까워지고, 화학이 빠르면( $\mathrm{Da}\gg1$ ) 국소 커널이 그대로 살아남아 불균질 착화가 된다. 그리고 하필 실용적으로 중요한 조건들이 $\mathrm{Da}\sim1$ 근처에 몰려 있다는 것도, 이 바닥의 익숙한 심술이다.³

9. 관련 문서[편집]

계산 쪽도 깨끗하지는 않다. $dT/dt$ 최대 시각을 쓰면 적분기의 출력 간격이 성기게 잡혀 있을 때 첨두를 놓쳐 $\tau_{ig}$ 가 계단처럼 찍힌다. “메커니즘 A가 B보다 3% 빠르다”는 주장이 사실은 출력 간격 차이였던 사례가 있다. 후처리 스크립트가 물리보다 결론에 크게 기여하는 순간. ↩
RON(리서치 옥탄가)과 MON(모터 옥탄가)은 같은 연료를 다른 회전수·흡기온도로 시험한 값이고, 둘의 차이(감도, sensitivity)가 큰 연료일수록 NTC 거동이 약하다. 에탄올이 대표적. 즉 옥탄가 두 개는 “이 연료의 아레니우스 선도가 어떻게 휘었는가”를 두 점으로 요약한 값에 가깝다. 주유소 간판에 화학반응 메커니즘의 그림자가 걸려 있는 셈. ↩
자착화 예측이 얼마나 어려운지를 보여 주는 통계 하나. 같은 조건에 대해 널리 쓰이는 상세 메커니즘 여럿을 돌리면 $\tau_{ig}$ 가 배 이상 갈리는 일이 드물지 않고, 그 편차가 가장 큰 곳이 하필 NTC 영역이다. 엔진 설계자가 “그래서 어느 걸 쓰면 되냐”고 물으면 연소학자는 조용히 실험 조건표를 꺼낸다. ↩

자착화