화학반응 메커니즘

계산화학 시뮬레이션 수치해석 마지막 수정: 2026-08-27 04:37:52

1. 개요[편집]

교과서의 총괄반응식은 회계장부다. 실제로 일어나는 일이 아니다.

화학반응 메커니즘(chemical reaction mechanism)은 하나의 총괄반응을 실제로 구성하는 소반응(elementary reaction)들의 집합과, 각 소반응의 속도계수·열역학 데이터를 함께 묶은 것이다. 시뮬레이션 하는 사람 입장에서 더 정확한 정의는 이쪽이다 — 메커니즘이란 화학종 목록, 반응 목록, 속도계수, 열역학 다항식, (필요하면) 수송 계수를 담은 하나의 데이터 파일이며, 그 파일이 있어야 비로소 반응 유동 코드가 돌아간다. 사실상의 표준 포맷은 Chemkin 입력 형식과 그 후계인 Cantera의 YAML이다.

왜 총괄반응으로는 안 되는가. 메테인 연소를 이렇게 쓴다고 하자.

\mathrm{CH_4} + 2\,\mathrm{O_2} \;\longrightarrow\; \mathrm{CO_2} + 2\,\mathrm{H_2O}

이건 원자 수를 맞춘 회계일 뿐 물리적 사건이 아니다. 메테인 분자 하나와 산소 분자 둘이 동시에 충돌해 정확히 이 조합으로 재배열될 확률은 사실상 0이다. 실제로는 라디칼 수십 종이 등장해 수백 번의 2분자 충돌을 릴레이로 이어가고, 그 릴레이의 구조가 착화 시점·화염 속도·오염물질 생성량을 전부 결정한다. 총괄식으로는 노킹도 NOx도 냉염도 설명할 수 없다.

속도식 일반론과 반응 차수는 반응속도론이, 온도 의존성의 물리는 아레니우스 방정식이 다룬다. 여기서는 메커니즘이라는 물건 자체 — 어떻게 생겼고, 어떻게 만들어지고, 왜 그렇게 비싸고, 어떻게 줄이는가 — 를 본다.

2. 소반응이라는 최소 단위[편집]

소반응은 중간체 없이 한 번의 사건으로 일어나는 반응이고, 등장하는 분자 수를 분자도(molecularity)라 한다.

분자도	예	성격
단분자	$\mathrm{H_2O_2} \to 2\,\mathrm{OH}$	결합 해리·이성질화. 압력 의존이 심하다
2분자	$\mathrm{H} + \mathrm{O_2} \to \mathrm{OH} + \mathrm{O}$	절대다수. 충돌 한 번에 원자 하나 옮기기
3분자	$\mathrm{H} + \mathrm{O_2} + \mathrm{M} \to \mathrm{HO_2} + \mathrm{M}$	재결합. 여분의 에너지를 3체 $\mathrm M$ 이 가져간다

4분자 이상은 사실상 없다고 본다 — 네 입자가 동시에 만날 확률이 무시할 만하기 때문이다. 소반응에서만 반응 차수가 화학량론 계수와 일치한다. 그래서 소반응 단위로 쪼개 놓으면 질량작용 법칙이 그대로 성립하고, 전체 계는 종 개수만큼의 상미분방정식이 된다.

\frac{d[X_k]}{dt} \;=\; \sum_{r} \nu_{k,r}\left( k_{f,r}\prod_j [X_j]^{\nu'_{j,r}} - k_{b,r}\prod_j [X_j]^{\nu''_{j,r}} \right)

역방향 계수 $k_b$ 를 따로 주는 일은 드물다. 미시적 가역성 때문에 $k_b = k_f / K_c(T)$ 이고, $K_c$ 는 열역학 데이터(대개 NASA 7항 또는 9항 다항식)에서 나온다. 이 규약 덕에 메커니즘은 자동으로 평형과 모순되지 않는다. 정방향과 역방향 계수를 각자 피팅해 넣으면 고온에서 조용히 열역학 제2법칙을 어기는 메커니즘이 완성된다.

3. 연쇄반응 — 개시·전파·분지·종결[편집]

연소 메커니즘의 뼈대는 라디칼 연쇄다. 역할별로 네 종류가 있다.

개시(initiation) — 안정한 분자에서 라디칼을 처음 만든다. $\mathrm{CH_4} + \mathrm{O_2} \to \mathrm{CH_3} + \mathrm{HO_2}$ 같은 것으로, 활성화 에너지가 커서 매우 느리다. 라디칼 농도가 0에서 출발한다는 사실이 유도기간의 존재 이유다.
전파(propagation) — 라디칼 하나가 들어가 라디칼 하나가 나온다. 라디칼 풀 크기는 그대로, 연료만 소모된다.
분지(branching) — 라디칼 하나가 들어가 둘 이상이 나온다. 여기가 폭발의 심장이다. 종결보다 분지가 빠르면 라디칼 수가 지수적으로 증가하고, 어느 순간 연료가 한꺼번에 소진되며 열이 풀린다.
종결(termination) — 라디칼이 서로 만나 안정한 분자가 되거나 벽에서 죽는다.

착화란 결국 분지가 종결을 이기는 순간이다. 세묘노프가 1930년대에 이 그림을 정립해 노벨 화학상을 받았고, 열폭주 관점의 판정 기준은 자착화 문서가 다룬다.

4. H₂/O₂ — 메커니즘 교과서의 표준 예제[편집]

수소-산소계는 종이 10개 남짓, 반응이 20여 개뿐인데도 메커니즘의 모든 개념이 들어 있어 늘 첫 예제로 나온다. 주인공은 두 반응이다.

\mathrm{H} + \mathrm{O_2} \;\rightarrow\; \mathrm{OH} + \mathrm{O} \quad(\text{분지}), \qquad \mathrm{H} + \mathrm{O_2} + \mathrm{M} \;\rightarrow\; \mathrm{HO_2} + \mathrm{M} \quad(\text{종결})

같은 반응물이 온도와 압력에 따라 정반대 역할을 한다는 것이 요점이다. 위쪽은 활성화 에너지가 커서 고온에서 이기고, 아래쪽은 3체 반응이라 압력(즉 $[\mathrm M]$ )에 비례해 커진다. 두 속도가 같아지는 교차 조건

2\,k_{\text{분지}}(T) \;\approx\; k_{\text{종결}}(T)\,[\mathrm M]

이 $p$ – $T$ 평면에 그리는 선이 그 유명한 제2 폭발한계다. 계수 2는 분지 한 번에 라디칼이 순증 하나가 아니라 둘 늘어나는 사슬 구조에서 온다.

전체 그림은 더 재미있다. 온도를 고정하고 압력을 올리면 비폭발 → 폭발 → 비폭발 → 폭발이 차례로 나타나 한계가 셋이 되고, $p$ – $T$ 평면에서 폭발 영역이 반도(peninsula) 모양으로 삐져나온다.

제1 한계 — 압력이 너무 낮으면 라디칼이 벽까지 확산해 죽는다. 벽 종결이 지배. 용기 크기와 재질에 값이 의존한다는, 화학치고는 굴욕적인 성질이 있다.
제2 한계 — 위의 3체 종결이 기상에서 분지를 이긴다.
제3 한계 — 압력이 더 오르면 죽은 줄 알았던 $\mathrm{HO_2}$ 가 되살아난다. $\mathrm{HO_2} + \mathrm{H_2} \to \mathrm{H_2O_2} + \mathrm{H}$ 로 이어지고 $\mathrm{H_2O_2} + \mathrm M \to 2\,\mathrm{OH}$ 가 새로운 분지 경로를 연다. 느린 분지라 열 축적과 결합해 폭발한다.

이 구조를 이해하면 저온·고압 연소의 절반은 이해한 것이다. 실제 엔진 조건(20~~60 bar, 700~~1000 K)은 정확히 3한계 영역이고, 그래서 $\mathrm{H_2O_2}$ 분해가 자착화 시점을 좌우한다.

5. 속도계수를 어떻게 적어 넣는가[편집]

수정 아레니우스. 표준 형식은 세 파라미터다.

k(T) \;=\; A\,T^{\,b}\exp\!\left(-\frac{E_a}{R T}\right)

$b$ 는 충돌 빈도의 온도 의존과 전이상태의 분배함수 효과를 흡수하는 자리로, 전이상태 이론에서 자연스럽게 나온다. 넓은 온도 구간을 하나의 아레니우스 식으로 못 맞추는 반응은 두 세트를 겹쳐 쓰기도 한다(이 경우 개별 $E_a$ 에 물리적 의미를 부여하면 안 된다).

3체 반응. $[\mathrm M]$ 은 모든 분자의 농도 합인데, 에너지를 받아가는 능력이 분자마다 달라 효율 계수를 곱한다. $[\mathrm M] = \sum_k \alpha_k [X_k]$ 이고 $\mathrm{H_2O}$ 는 대개 $\alpha \approx 10\text{--}20$ 이다. 연소 생성물인 물이 스스로 3체로서 반응을 바꾼다는 뜻이라, 습도나 EGR 비율이 착화에 미치는 영향의 상당 부분이 여기서 나온다.

압력 의존(falloff). 단분자 분해와 재결합은 압력에 따라 반응 차수 자체가 변한다. 저압에서는 활성화 충돌이 율속이라 2차, 고압에서는 활성화된 분자의 분해가 율속이라 1차다. 린데만-힌셜우드가 이 전이를 처음 모형화했고, 실제 데이터는 그 곡선보다 완만해서 보정 인자가 필요하다.

형식	성격
Lindemann	$k = k_\infty \dfrac{P_r}{1+P_r}$ , $P_r = k_0[\mathrm M]/k_\infty$ . 전이 구간이 실제보다 급하다
Troe	위 식에 폭 인자 $F_{\text{cent}}$ 를 곱해 전이를 완만하게. 파라미터 3~4개. 사실상 표준
SRI	Troe의 대안. 형태만 다르고 목적은 같다
PLOG	여러 압력에서 각각 아레니우스 3인자를 주고, $\log p$ 로 선형 보간. 이론 모형 없이 계산 결과를 그대로 담을 때
Chebyshev	$(\log T,\log p)$ 2차원 다항식 피팅. 가장 유연하고 가장 외삽에 취약

PLOG와 Chebyshev가 2000년대 이후 급격히 늘어난 이유는 명확하다. 속도계수를 실험이 아니라 양자화학 + 마스터 방정식 계산으로 얻는 것이 표준이 되면서, 계산 결과가 Troe 형식에 안 들어맞는 경우가 흔해졌기 때문이다. 코드가 이 형식들을 지원하지 않으면 최신 메커니즘을 아예 읽을 수 없다.

6. 실제 메커니즘은 얼마나 큰가[편집]

메커니즘	종	반응	용도
H₂/O₂	~10	~25	개념 학습, 수소 연소
GRI-Mech 3.0	53	325	천연가스 연소 + NOx. 가장 많이 쓰인 메커니즘
iso-옥테인 상세(LLNL, 2002)	857	3606	옥탄가 기준연료의 저온~고온 산화
가솔린·디젤 대체연료	1000~5000+	수만	실연료 모사

GRI-Mech 3.0은 특이한 물건이다. 개별 속도계수를 문헌에서 그대로 가져다 쓴 것이 아니라, 착화지연·화염속도·종 프로파일 같은 목표 실험 수십 건에 대해 불확실 파라미터를 전역 최적화해 맞춘 메커니즘이다. 그래서 목표 조건 안에서는 놀랄 만큼 잘 맞고, 밖으로 나가면 급격히 나빠진다.¹ 메커니즘을 쓰기 전에 개발자가 검증한 조건 범위(온도·압력·당량비·연료)를 확인하는 것이 반드시 선행돼야 하는 이유다.

저온 산화의 RO₂ 화학이 크기 폭발의 주범이다. 알케인 $\mathrm{RH}$ 가 700 K 근처에서 겪는 경로는 대략 이렇다.

\mathrm R + \mathrm{O_2} \rightleftharpoons \mathrm{RO_2} \;\to\; \mathrm{QOOH} \;\xrightarrow{+\mathrm{O_2}}\; \mathrm{O_2QOOH} \;\to\; \text{케토하이드로퍼옥사이드} + \mathrm{OH} \;\to\; \text{분해} + 2\,\mathrm{OH}

산소를 두 번 붙여야 라디칼이 순증하는 저온 분지가 완성된다. 문제는 $\mathrm{RO_2}$ 의 내부 수소 이동이 어느 탄소에서 일어나느냐에 따라 $\mathrm{QOOH}$ 이성질체가 여러 개 생기고, 그 각각이 또 여러 갈래로 갈라진다는 것이다. 탄소 수가 늘면 이성질체 수가 조합적으로 폭발한다. iso-옥테인 하나에 857종이 필요한 이유가 이것이고, 왜 이 화학이 냉염과 NTC 거동, 나아가 옥탄가의 물리적 근거인지도 여기서 나온다.

7. 왜 이렇게 비싼가 — 강성[편집]

메커니즘을 CFD에 붙이면 계산비가 폭발한다. 이유는 두 가지다.

첫째, 강성. 라디칼의 특성 시간은 나노초 수준인데 관심 있는 현상(착화, 화염 전파)은 밀리초 이상이다. 시간척도 비가 $10^6$ ~ $10^9$ 에 이르는 전형적인 강성 방정식이고, 명시적 룽게-쿠타법을 쓰면 안정성 때문에 시간 간격이 가장 빠른 라디칼에 묶여 버린다. 그래서 화학 적분은 거의 예외 없이 후진 미분 공식(BDF) 기반 암시적 적분기(DVODE, CVODE 계열)를 쓴다. 암시적이라는 것은 매 스텝 뉴턴 반복을 돈다는 뜻이고, 뉴턴 반복은 자코비안 행렬과 그 LU 분해를 요구한다.

둘째, 규모. 종이 $N_s$ 개면 야코비안은 $N_s \times N_s$ 이고 조밀 LU는 $O(N_s^3)$ 이다. 종 수가 50에서 500으로 늘면 비용은 1000배가 된다. 여기에 CFD 격자 셀 하나하나가 독립적인 ODE 계라는 사실을 곱하면, 실제 연소 해석에서 화학이 전체 실행 시간의 대부분(흔히 80% 이상)을 먹는다는 결론이 나온다. 연소 시뮬레이션의 비용 구조가 다른 CFD와 근본적으로 다른 이유다.

완화 수단은 다음과 같다.

연산자 분리. 수송과 반응을 분리해 각각 적합한 적분기로 푼다(스트랑 분할). 표준 관행이지만 분리 오차가 강성 화학에서 무시 못 할 수 있다.
야코비안 희소화·해석적 야코비안. 대부분의 종은 서로 직접 반응하지 않아 야코비안이 희소하다. 희소행렬 분해로 $O(N_s^3)$ 을 크게 낮춘다. 수치 미분 대신 해석적 야코비안을 코드 생성으로 만드는 것도 표준.
ISAT / 표 작성. 이미 계산한 조성 공간의 점 근처를 선형 근사로 재사용한다. 조성 공간을 저차원 다양체로 미리 표로 만들어 두는 flamelet/FGM 계열도 같은 발상.
동적 적응 화학(DAC). 셀마다 그 순간 필요한 종만 남기고 나머지는 끈다.

8. 메커니즘 축소[편집]

가장 확실한 절약은 메커니즘 자체를 줄이는 것이다. 접근은 두 갈래다.

골격 축소(skeletal reduction) — 종과 반응을 지운다.

DRG(Directed Relation Graph, 루·로 2005). 종을 노드로 두고, 종 $A$ 를 지웠을 때 $B$ 의 생성률이 얼마나 망가지는지를 간선 가중치로 삼아 유향 그래프를 만든다. 관심 종(연료, 산소, 열방출 지표)에서 출발해 임계값 이상의 간선을 따라 도달 가능한 종만 남긴다. 그래프 탐색 한 번이라 값이 거의 공짜이고, 임계값을 올려가며 종 수-오차 곡선을 그릴 수 있다는 것이 결정적 장점이다.
DRGEP(with Error Propagation, 페피오-데자르댕·피치 2008). DRG는 경로를 따라가며 오차가 계속 누적된다고 보는데, 실제로는 경로가 길어질수록 영향이 곱해져 줄어든다. DRGEP는 간선 가중치를 곱으로 전파해 이를 반영하고, 같은 오차 허용치에서 DRG보다 더 작은 메커니즘을 낸다.
민감도 해석. 위 방법으로 후보를 줄인 뒤 남은 종에 대해 직접 민감도를 계산해 마무리한다. 비싸지만 정확하다.

시간척도 축소 — 빠른 자유도를 대수식으로 바꾼다.

QSSA(준정상상태 근사). 생성과 소멸이 모두 매우 빠른 라디칼은 $d[X]/dt \approx 0$ 으로 두고 대수방정식으로 푼다. 미분방정식 개수가 줄고 강성이 직접 제거되는 것이 핵심 이득이다. 대신 QSSA 종끼리 얽히면 비선형 연립 대수방정식이 생겨 그것대로 골치다.
부분평형 근사(PEA). 정·역방향이 모두 빠른 반응 쌍을 평형으로 놓는다.
ILDM(마스·포프 1992). 화학 자코비안의 고윳값을 빠른 것과 느린 것으로 나누고, 빠른 방향이 이미 이완된 저차원 끌개 다양체를 조성 공간에서 미리 계산해 표로 만든다. 자유도를 2~3개로 줄여 버린다.
CSP(Computational Singular Perturbation, 램·구시스). 같은 아이디어를 특이섭동 이론으로 체계화해, 어떤 종이 QSSA 대상이고 어떤 반응이 부분평형인지를 사람의 직관이 아니라 계산으로 판정한다.

실무 파이프라인은 대개 “상세 메커니즘 → DRGEP로 골격 축소 → QSSA로 강성 제거 → 목표 조건에서 착화지연·화염속도 재검증”으로 굳어져 있다. 줄인 메커니즘은 반드시 원본과 대조 검증한다 — 착화지연 오차 5% 안쪽으로 맞춘 축소 메커니즘이 NOx는 3배 틀리는 일이 흔하다. 축소는 목표를 정해 놓고 하는 작업이지 범용 압축이 아니다.²

9. 어떻게 믿을 것인가[편집]

메커니즘은 이론적 구조물이 아니라 검증 대상 모델이다. 표준 검증 실험은 정해져 있다.

착화지연 — 충격파관(고온, 1000 K 이상)과 급속압축기(저온·고압, 엔진 조건). 저온 영역의 2단 착화가 RO₂ 화학을 직접 겨눈다.
층류 화염 속도 — 열방출과 수송의 결합을 본다. 화염 전파 예측력의 1차 지표.
종 농도 프로파일 — 완전혼합 반응기(JSR)나 층류 유동 반응기에서 중간체를 직접 잰다. 착화지연만 맞춘 메커니즘은 여기서 자주 무너진다.
소염 스트레인율 — 대향류 화염에서 잰다. 난류 연소 모형이 쓰는 양이다.

여기에 각 속도계수의 불확실성을 확률분포로 두고 예측의 불확실성을 전파하는 불확실성 정량화가 최근 표준으로 자리잡았다. 결국 메커니즘 개발은 역문제다 — 관측(착화지연, 화염속도)에서 파라미터(속도계수 수백 개)를 추정하는 문제이고, 파라미터가 관측보다 훨씬 많으므로 답이 유일하지 않다. 서로 다른 두 메커니즘이 같은 검증 데이터를 똑같이 잘 맞추면서 새 조건에서는 다른 예측을 내놓는 일이 실제로 벌어진다.³

10. 관련 문서[편집]

GRI-Mech 3.0의 공식 문서는 “메테인·천연가스 연소용이며 그 밖의 연료에 쓰지 말라”고 명시하고 있다. 그럼에도 이 메커니즘이 프로페인·에탄올·심지어 등유 해석 논문에 등장하는 광경은 20년째 반복된다. 파일이 공짜고 작아서 그렇다. ↩
메커니즘 축소 도구를 쓸 때 가장 흔한 사고는 목표 조건을 좁게 잡는 것이다. 당량비 1.0, 1 atm에서만 축소해 놓고 EGR 40%짜리 엔진 해석에 넣으면, 지운 종들이 정확히 그 조건에서 중요했다는 사실을 결과가 다 나온 뒤에 알게 된다. ↩
그래서 연소 학계는 검증 데이터를 표준화하는 데 상당한 공을 들여 왔다. 데이터셋과 메커니즘이 함께 공개되지 않은 예측 결과는 재현할 방법이 없다 — “우리 사내 메커니즘으로 계산했습니다”는 V&V 관점에서 아무 말도 하지 않은 것과 같다. ↩