LDF Approximation의 Mass Transfer Coefficient는 실제로 무엇을 나타내는가

흡착공정 모델을 보다 보면 복잡한 확산방정식 대신 하나의 간단한 식이 반복해서 등장합니다. 입자 내부의 평균 흡착량이 평형값을 향해 얼마나 빠르게 접근하는지를 하나의 계수로 표현하는 Linear Driving Force(LDF) 근사입니다.

계산은 간단하지만 여기서 자주 생기는 오해가 있습니다. LDF의 물질전달계수(mass transfer coefficient)를 흡착제 고유의 확산계수와 같은 값으로 받아들이는 것입니다.

둘은 같은 물리량이 아닙니다. LDF 계수는 분자가 기공 안에서 몇 m²/s의 확산계수를 갖는지를 직접 말해주는 값이라기보다, 주어진 입자와 실험조건에서 흡착량이 평형을 향해 완화되는 전체 속도를 하나의 시간척도로 압축한 값에 가깝습니다.

복잡한 입자 내부 농도분포를 하나의 평균값으로 바꾼다

실제 흡착입자를 생각하면 입자 내부의 모든 위치에서 흡착량이 동시에 변하지 않습니다.

외부 기체의 농도가 갑자기 바뀌면 입자 표면 근처가 먼저 반응하고, 시간이 지나면서 변화가 내부로 전달됩니다. 따라서 입자 중심과 표면 사이에는 일시적인 농도분포가 형성됩니다.

이를 엄밀하게 계산하려면 입자 내부의 위치와 시간에 따른 확산방정식을 풀어야 합니다.

LDF 근사는 이 공간분포를 직접 추적하는 대신 입자 전체의 평균 흡착량을 사용합니다. 가장 기본적인 형태는 다음과 같습니다.

dq̄/dt = k_LDF(q* − q̄)

여기서 q̄는 입자 내부의 평균 흡착량, q*는 현재 외부 조건과 평형을 이루는 흡착량입니다. 두 값의 차이인 q* − q̄가 평형에서 얼마나 떨어져 있는지를 나타내는 구동력(driving force)이고, k_LDF가 그 차이가 얼마나 빠르게 줄어드는지를 결정합니다.

즉 LDF는 실제 농도분포가 선형이라는 뜻이 아닙니다. 복잡한 입자 내부 확산의 결과를 평균 흡착량 하나의 동역학으로 바꾸는 근사입니다.

k_LDF의 단위가 말해주는 것이 있다

k_LDF의 단위는 보통 시간의 역수입니다.

예를 들어 s⁻¹로 표현된다면 이 값은 거리당 분자 이동속도를 직접 나타내는 것이 아니라 흡착입자가 평형을 향해 얼마나 빠르게 반응하는지를 나타내는 속도상수와 비슷한 역할을 합니다.

k_LDF가 크면 q̄와 q*의 차이가 빠르게 줄어들고, 작으면 입자가 새로운 평형상태에 접근하는 데 더 긴 시간이 필요합니다.

이 점만 봐도 m²/s 단위를 갖는 확산계수 D와 k_LDF가 동일한 물리량이 아니라는 것을 알 수 있습니다.

그렇다고 둘이 관계가 없는 것은 아닙니다. 구형 입자에서 내부 확산이 지배적이고 일정한 유효확산계수를 가정하면 LDF 계수를 다음과 같은 형태로 연결할 수 있습니다.

k_LDF ≈ 15D_eff / Rₚ²

여기서 D_eff는 모델에서 사용하는 유효확산계수(effective diffusivity), Rₚ는 입자 반지름입니다. 계수 15는 구형 입자의 확산 문제를 LDF 형태로 근사할 때 나타나는 대표적인 값이며, 사용하는 확산모델과 기하학적 정의에 따라 관계식은 달라질 수 있습니다.

중요한 것은 k_LDF가 D만으로 결정되지 않는다는 점입니다. 같은 확산계수를 가진 재료라도 입자 크기가 달라지면 물질전달의 시간척도가 달라집니다.

입자가 커지면 같은 재료도 느리게 보인다

위 관계에서 특히 눈에 띄는 것은 Rₚ²입니다.

동일한 흡착제를 사용하면서 입자 반지름만 두 배로 증가했다고 생각해보겠습니다. 내부 확산계수가 그대로라면 분자가 이동해야 하는 대표적인 거리가 길어지기 때문에 입자 전체가 평형에 접근하는 시간도 크게 증가합니다.

따라서 서로 다른 크기의 펠릿에서 얻은 k_LDF를 단순 비교해 한쪽 재료의 분자 확산이 더 느리다고 결론내리면 안 됩니다.

분말 결정과 성형 펠릿을 비교할 때는 문제가 더 복잡해집니다. 결정 내부의 미세기공 확산뿐 아니라 결정 사이의 공극, 바인더, 거대기공 수송 같은 저항이 추가될 수 있기 때문입니다.

이 경우 피팅으로 얻은 하나의 LDF 계수는 어느 하나의 미세한 이동 과정만 가리키기보다 여러 저항이 합쳐진 입자 규모의 응답을 나타낼 수 있습니다.

실험에서 얻은 k_LDF는 무엇을 포함하고 있을까

실험 데이터를 LDF 모델에 맞춰 k_LDF를 얻었다고 해도 그 값을 바로 결정내 확산계수로 환산하기 전에 어떤 과정이 측정응답을 지배하는지 확인해야 합니다.

외부 막 물질전달(external film mass transfer)이 충분히 빠르지 않다면 기체 벌크에서 입자 표면으로 이동하는 과정이 전체 속도에 영향을 줍니다. 펠릿 내부에 거대기공과 미세기공이 함께 존재한다면 서로 다른 길이 스케일의 저항도 겹칠 수 있습니다.

강한 흡착열이 발생하는 조건에서는 온도 변화까지 흡착속도에 영향을 줄 수 있습니다. 이때 실험곡선 하나를 단순 LDF 모델에 잘 맞췄다고 해서 얻어진 계수가 순수한 미세기공 확산만을 나타낸다고 보기는 어렵습니다.

그래서 k_LDF는 어떻게 얻었는지가 중요합니다. 알려진 D와 입자 크기에서 계산한 계수인지, 특정 실험의 응답을 피팅해 얻은 계수인지에 따라 물리적인 해석의 범위가 달라집니다.

하나의 k_LDF가 모든 압력에서 유지될까

LDF 식이 간단하다는 이유로 하나의 k_LDF를 전체 흡착량 범위에 사용하는 경우가 있습니다. 하지만 실제 흡착계에서는 이것이 항상 성립하지 않습니다.

앞에서 살펴본 것처럼 미세기공의 분자 이동성은 흡착량에 따라 변할 수 있습니다. 흡착등온선의 비선형성이 강하면 외부 농도 변화와 입자 내부 흡착량 변화의 관계도 일정하지 않습니다.

따라서 넓은 압력 범위에서는 일정한 k_LDF 하나가 모든 동역학을 정확하게 표현하지 못할 수 있습니다. 필요하다면 흡착량이나 압력에 따라 달라지는 LDF 계수를 사용하거나, 입자 내부 확산을 직접 계산하는 모델과 비교해야 합니다.

이 때문에 LDF 모델의 정확성은 식이 얼마나 단순한지가 아니라 해당 조건에서 하나의 대표 시간척도로 실제 수송을 얼마나 잘 압축할 수 있는지에 달려 있습니다.

LDF는 확산을 없애는 모델이 아니라 압축하는 모델이다

LDF 근사를 사용하면 입자 내부의 상세한 농도분포는 보이지 않습니다. 그렇다고 내부 확산의 물리가 사라지는 것은 아닙니다.

입자 크기, 유효확산계수, 기공 구조와 여러 물질전달 저항이 만들어낸 결과가 k_LDF라는 하나의 계수에 압축되어 들어갑니다.

그래서 LDF 계수를 볼 때 중요한 질문은 ‘이 값이 얼마인가’만이 아닙니다. 이 계수가 어떤 입자 크기와 압력, 온도, 흡착량에서 얻어졌는지, 그리고 어떤 수송 저항을 대표하도록 정의됐는지를 함께 확인해야 합니다.

이 구분이 명확하면 LDF는 복잡한 입자 내부 확산을 PSA나 고정층 같은 공정 모델에 연결하는 매우 효율적인 근사가 됩니다. 반대로 k_LDF를 재료 고유의 확산계수처럼 사용하면 서로 다른 실험이나 입자에서 얻은 값을 잘못 비교할 수 있습니다.

참고자료

Glueckauf & Coates, Journal of the Chemical Society (1947)

Ruthven, Principles of Adsorption and Adsorption Processes (1984)

Kärger, Ruthven & Theodorou, Diffusion in Nanoporous Materials (2012)

Sircar & Hufton, Adsorption (2000)