다성분 Maxwell-Stefan Cross-Term을 무시하면 분리 예측이 달라지는 이유

두 성분 A와 B가 같은 제올라이트를 통과한다고 생각해보겠습니다. 순수성분 실험에서는 A가 B보다 훨씬 빠르고, 흡착평형까지 고려했을 때도 A가 먼저 빠져나올 것으로 예상됩니다.

그런데 실제 혼합물에서는 예상한 차이가 줄어들거나 특정 조건에서 한 성분의 플럭스가 크게 달라질 수 있습니다. 순수성분 확산계수를 잘못 측정한 것도 아니고 기공 구조가 갑자기 변한 것도 아닙니다.

문제는 두 성분을 각각 독립적으로 움직이는 분자로 취급했다는 데 있을 수 있습니다. 같은 미세기공을 공유하는 순간 A의 이동은 B와 완전히 분리된 문제가 아니게 됩니다.

A의 플럭스에 왜 B의 구동력이 들어올까

가장 단순한 확산 모델에서는 A의 플럭스를 A의 농도구배와 연결하고, B의 플럭스는 B의 농도구배와 따로 연결합니다.

다성분 수송(multicomponent transport)에서는 이런 분리가 항상 성립하지 않습니다. 일반화된 형태로 플럭스와 구동력의 관계를 나타내면 다음처럼 생각할 수 있습니다.

J_A = −D_AA∇c_A − D_AB∇c_B

여기서 첫 번째 항은 A 자신의 농도 변화가 A의 플럭스에 미치는 영향이고, 두 번째 항은 B의 농도 변화가 A의 플럭스에 미치는 영향을 나타냅니다.

D_AB처럼 다른 성분의 구배와 연결되는 항이 교차항(cross-term)입니다.

실제 Maxwell–Stefan 모델에서는 단순한 농도구배보다 화학퍼텐셜 구배(chemical potential gradient)에 의한 구동력과 성분 사이의 상대운동을 이용해 문제를 기술하지만, 교차항이 말해주는 물리적 의미는 분명합니다. A의 이동을 계산할 때 A만 보고 끝낼 수 없다는 것입니다.

교차항은 단순한 계산상의 보정값이 아니다

Maxwell–Stefan 관점에서는 각 성분이 받는 구동력이 여러 종류의 마찰과 균형을 이룹니다.

A는 제올라이트에 대해 움직이면서 저항을 받고, 동시에 B와 서로 다른 속도로 움직이려 할 때도 추가적인 저항을 받습니다. 이때 A와 B 사이의 상대운동에 해당하는 부분에는 두 성분 사이의 Maxwell–Stefan 교환확산계수(exchange diffusivity) Đ_AB가 관여합니다.

A와 B의 속도가 비슷하면 상대운동은 작습니다. 반대로 A가 빠르게 이동하려 하고 B가 거의 움직이지 않는다면 두 성분의 상대속도는 커집니다. 이때 성분 간 마찰을 무시하면 A가 실제보다 자유롭게 이동하는 것으로 계산될 수 있습니다.

따라서 교차항은 결과를 맞추기 위해 나중에 추가하는 경험적인 보정값이라기보다, 혼합물 안에서 성분들의 이동이 서로 연결된다는 사실을 수송식에 반영하는 부분입니다.

이동성이 크게 다를수록 오차가 커질 수 있다

모든 혼합물에서 교차항이 똑같이 중요하지는 않습니다.

두 성분의 이동성이 비슷하고 서로의 상대운동에 따른 저항이 작다면 교차항을 생략해도 결과가 크게 달라지지 않는 조건이 있을 수 있습니다.

하지만 빠른 A와 느린 B가 좁은 기공을 함께 이용한다면 상황이 달라집니다. A만 놓고 계산하면 빠르게 통과할 것으로 예상되지만, 실제 혼합물에서는 느린 B와의 상대운동 때문에 A의 이동도 영향을 받을 수 있습니다.

특히 흡착량이 증가해 두 성분이 같은 이동 경로를 자주 공유하거나 서로 지나치기 어려운 구조에서는 이런 운동학적 연결(kinetic coupling)이 더 중요해질 수 있습니다.

그래서 순수성분 확산계수(pure-component diffusivity) 두 개를 혼합물 모델에 그대로 넣는 것만으로는 충분하지 않은 경우가 생깁니다. 순수한 A에는 A-B 상대운동이 존재하지 않기 때문에 혼합된 뒤 새롭게 나타나는 성분 간 저항을 그 데이터만으로 직접 알 수 없기 때문입니다.

열역학적 교차효과와 마찰에 의한 연결은 구분해야 한다

여기서 19번에서 다룬 역농도 확산과 연결되는 부분이 있습니다.

다성분 흡착에서는 B의 조성 변화가 A의 화학퍼텐셜을 바꿀 수 있습니다. 이것은 혼합흡착 평형에서 비롯되는 열역학적 연결(thermodynamic coupling)입니다.

반면 Maxwell–Stefan의 성분 간 마찰은 A와 B가 서로 다른 속도로 움직일 때 발생하는 상대운동의 문제입니다. 이것은 운동학적 연결에 해당합니다.

두 효과는 최종 플럭스 식에서는 함께 나타날 수 있지만 물리적인 출발점은 다릅니다.

따라서 계산에서 A의 플럭스가 예상과 다르다고 해서 모든 차이를 Đ_AB 하나로 설명해서는 안 됩니다. 먼저 혼합흡착 열역학이 구동력을 얼마나 바꿨는지 확인하고, 그다음 성분 간 상대운동에 의한 연결이 얼마나 추가되는지를 살펴봐야 합니다.

이 둘을 구분하지 않으면 모델은 실험값과 비슷하게 맞더라도 왜 그런 플럭스가 나왔는지는 잘못 해석할 수 있습니다.

플럭스 차이가 결국 분리성능 차이로 이어진다

다성분 수송에서 중요한 것은 확산계수 자체보다 그 결과로 각 성분의 플럭스가 어떻게 달라지는가입니다.

A와 B의 플럭스 비가 바뀌면 막이나 결정에서 나타나는 동역학적 선택성(kinetic selectivity)도 달라질 수 있습니다. 고정층에서도 성분별 이동속도가 달라지면 농도 전선의 진행과 파과 시점에 영향을 줄 수 있습니다.

교차항을 무시한 모델이 빠른 성분 A의 이동을 지나치게 크게 계산한다면 실제보다 높은 분리성능을 예측할 수 있습니다. 반대로 성분 사이의 연결이 특정 조건에서 이동을 돕는다면 교차항을 제거한 모델이 플럭스를 낮게 계산하는 경우도 가능합니다.

즉 교차항의 효과는 항상 ‘확산을 느리게 만든다’는 한 방향의 보정으로 이해하면 안 됩니다. 조성과 흡착량, 각 성분의 이동성, 성분 간 상대운동에 따라 결과가 달라질 수 있습니다.

교차항을 무시해도 되는지는 계산 전에 정해져 있지 않다

다성분 모델을 단순하게 만들기 위해 교차항을 생략하는 것은 가능하지만, 그것이 항상 안전한 가정은 아닙니다.

두 성분의 이동성이 얼마나 다른지, 흡착량이 어느 정도인지, 같은 기공 경로를 얼마나 강하게 공유하는지, 혼합물 조성이 변할 때 각 성분의 플럭스가 얼마나 민감하게 반응하는지를 먼저 살펴봐야 합니다.

무엇보다 순수성분 확산 데이터만 잘 맞는다는 사실은 다성분 모델도 정확하다는 증거가 아닙니다. 혼합물이 된 뒤에만 나타나는 상대운동과 성분 간 연결은 순수성분 실험에서는 검증할 수 없기 때문입니다.

Maxwell–Stefan 교차항의 의미도 여기에 있습니다. A가 움직일 때 B를 함께 계산해야 한다는 것은 단순히 식을 복잡하게 만드는 문제가 아닙니다. 같은 미세기공을 공유하는 순간 두 성분의 이동이 서로의 수송 조건이 된다는 뜻입니다.

따라서 교차항을 생략하기 전에 물어야 할 것은 ‘계산이 얼마나 간단해지는가’보다 ‘이 혼합물에서 A와 B를 정말 독립적으로 움직이는 성분으로 취급해도 되는가’입니다.

참고자료

Taylor & Krishna, Multicomponent Mass Transfer (1993)

Krishna, Chemical Society Reviews (2012)

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