양액 고체 반응 모델의 양자화 해법

초록

본 논문은 부피반응, 입자모델, 무작위공극모델, 핵생성모델 등 대표적인 유체‑고체 반응 모델에 대해 최근 제안된 양자화 방법(Quantized Method)을 적용하고, 전통적인 수치 해법과의 정확도·계산 효율을 비교한다. 두 개의 비선형 결합 편미분 방정식을 시간 구간별로 상수화하여 선형화함으로써 빠른 근사 해를 얻으며, 다양한 매개변수 조건에서 검증 결과를 제시한다.

상세 분석

유체‑고체 반응은 석유 정제, 금속 제련, 촉매 재생 등 산업 현장에서 핵심적인 역할을 한다. 이러한 반응을 기술하기 위해서는 반응물의 확산, 표면 흡착, 내부 반응 속도 등을 동시에 고려해야 하므로, 일반적으로 부피반응 모델(Volume Reaction Model), 입자모델(Grain Model), 무작위공극 모델(Random Pore Model), 핵생성 모델(Nucleation Model)과 같은 다중 물리·화학 현상을 포함하는 모델이 사용된다. 각 모델은 반응 진행에 따라 고체의 구조적 변화와 반응면적 감소를 다르게 가정하지만, 공통적으로 두 개의 비선형 결합 편미분 방정식(CPDE)—하나는 물질 확산을, 다른 하나는 반응 진행을 기술—을 포함한다. 이러한 방정식은 비선형성, 경계 조건의 복합성, 그리고 시간·공간 변수의 강한 결합 때문에 직접 해석이 불가능하고, 전통적인 유한 차분법(FDM), 유한 요소법(FEM) 혹은 전산 유체역학(CFD) 기반 수치 해법에 의존한다. 그러나 고해상도 격자와 작은 시간 스텝을 요구하므로 계산 비용이 급격히 증가하고, 수렴성 문제가 발생하기도 한다.

양자화 방법은 이러한 문제점을 완화하기 위해 ‘시간을 아주 작은 구간으로 양자화(quantize)’하고, 각 구간 내에서는 물성치와 반응계수를 고정값으로 간주한다. 구체적으로는 (1) 전체 시간 축을 Δt 단위의 미소 구간으로 분할하고, (2) 각 구간에서 비선형 항을 이전 구간의 해값을 이용해 선형화한다. 이렇게 하면 원래의 비선형 CPDE가 매 구간마다 선형 상미분·편미분 방정식 체계로 변환되어, 해석적 혹은 반해석적 해법(예: 라플라스 변환, 변분 원리)으로 빠르게 풀 수 있다. 양자화 단계가 충분히 작으면 실제 연속 해와의 차이가 통계적으로 무시될 정도로 작아지며, 동시에 계산량은 전통적 수치법 대비 1~2 차수 감소한다.

논문에서는 네 가지 모델 각각에 대해 양자화 해법을 적용한 절차를 상세히 제시한다. 부피반응 모델에서는 반응률 상수와 확산 계수를 각각 고정하고, 라플라스 변환을 이용해 농도 프로파일을 구한 뒤, 반응 진행도(α)를 시간에 따라 업데이트한다. 입자모델에서는 입자 내부의 구역을 ‘핵‑껍질’ 구조로 가정하고, 핵 성장 속도를 양자화된 α에 의존하도록 설정한다. 무작위공극 모델은 공극 부피와 표면적 감소를 경험적 식으로 표현하고, 이를 양자화 단계마다 재계산한다. 핵생성 모델은 핵 생성 속도와 성장 속도를 별도 식으로 분리하고, 각각을 시간 구간별 상수화한다.

각 모델에 대해 실험 데이터와 기존 수치 해법(전형적인 Crank‑Nicolson 스킴)과의 비교를 수행하였다. 결과는 양자화 해법이 평균 상대 오차 25% 이내로 높은 정확도를 유지하면서도, 동일한 시간·공간 해상도에서 계산 시간이 3070% 정도 감소함을 보여준다. 특히 고반응성(큰 Damköhler 수) 상황에서 전통적 방법은 시간 스텝을 극히 작게 잡아야 하는 반면, 양자화 방법은 안정적인 결과를 얻을 수 있다. 다만, 매우 급격한 구조 변화(예: 급격한 공극 붕괴)가 발생하는 경우에는 양자화 구간을 더 세분화해야 하며, 이때는 계산 효율이 다소 감소한다는 한계점도 논의된다.

전반적으로 양자화 방법은 복잡한 비선형 CPDE를 다루는 유체‑고체 반응 모델에 대해 ‘빠르고 충분히 정확한’ 근사 해를 제공하는 유망한 접근법이며, 향후 다중 반응·다중 물질 시스템, 3‑차원 형상, 그리고 실시간 공정 제어와 같은 분야에 확장될 가능성이 크다.