DOI QR코드

DOI QR Code

Viscosity of Binary Gas Mixture from the Calculation by Using the Brake Theory of Viscosity

Brake 점성이론으로 계산한 이성분기체의 점성

  • Kim, Won-Soo (Department of Metallurgical Engineering, Hongik University)
  • 김원수 (홍익대학교 금속공학과)
  • Published : 2004.06.20

Abstract

Brake theory of viscosity, which can sucessfully calculate the viscosity of real gases, dense gases and liquids, is extended to the binary gas mixture. Adjustable parameters are not involved, but the calculated results are good agreements with the experimental values at high pressure as well as low pressure. Corresponding state equation for viscosity can be obtained by using the Redlich-Kwong equation, so that we hope this equation may be useful for the supercritical fluid in engineering applications at high pressure around the critcal point.

Keywords

Binary Gas Mixture;Gas Viscosity;Corresponding State Equation

서 론

1860년 Maxwell이 낮은 압력하의 기체의 점성에 관한 점성식을 제시하였으나 이상기체에서만 적용될 뿐 실제기체에서는 그 실용성을 기대할 수 없었다. 분자간에 인력과 척력이 작용되는 실제기체에서는 Chapmann과 Enskog1,2가 최초로 분자간의 이중출돌과 탄성충돌등을 포함한 몇가지 가정을 하며 실용화를 위한 새로운 이론적 접근을 시도하였다. 그러나 실제로 기체의 점성을 계산하기 위해서는 분자간의 충돌직경과 collision integral값을 알지 아니하면 아니되었다. 분자간 충돌직경을 계산하기 위하여 여러 형태의 potential function이 검토되었고, collision integral 계산을 위한 수많은 연구3,4가 진행되었으나 만족할만한 결과가 발견되지 아니하였으며, 1972년 Neufeld 등5이 computer simulation을 통한 실험식의 형태로 충돌직경과 collison integral을 계산하기에 이르렀다. 이로부터 고압이 아닌 실제기체의 점성을 거의 완벽하게 계산 할 수 있었으며 낮은 압력하에서의 이성분기체의 점성에 대한 계산의 길이 열리게 된 것이다. 그러나 분자간의 다중충돌이 예상되는 dense gas의 경우는 계산이 더욱 복잡하였기에, Neufeld 이후 30년이 지난 현금에도 공업적이고 실용성이 있는 식들이 출현되고 있지 못하고 있는 실정이다.

Formal한 연구형태는 아니지만 각종 점성의 abnormal behavior6,7를 잘 설명하고 낮은 압력에서는 물론 dense gas 영역에서까지 점성의 계산이 가능했던 brake viscosity model8은 이성분기체의 경우에도 낮은 압력에서의 점성계산은 물론 dense gas영역에서도 통용될 수 있을 것으로 기대된다.

본 연구에서는 점성에 관한 brake viscosity model을 이용하여 이성분기체의 경우 낮은 압력에서 점성을 계산하고 실험값과 비교한 후 그 타당성이 인정되면 dense gas영역에서의 이성분기체의 점성도 계산코저 한다. dense gas 영역에서 그 계산값이 실험값과 잘 일치한다면 다양한 공업적 응용성을 기대할 수 있을 것이다. 실제로 대부분의 기체는 이성분이상의 혼합물로 이루어져 있으므로, 낮은 압력이건 높은 압력이건 다성분기체의 점성계산이 필요함은 물론이다. 초임계유체의 경우 밀도는 액체와 비슷한 값을 가지나 점성의 값은 훨씬 작아 기체의 점성의값에 가까우며 확산계수의 값은 매우 크므로 추출 기술의 중요한 변수가 됨은 물론 첨단분리기술의 하나로 기대될 수 있다. 이러한 의미에서 dense gas 영역에서의 이성분기체의 점성계산은 매우 바람직한 연구가 될 것이다.

 

이 론

유체의 단위면적당 작용하는 힘을 F, 이로부터 발생하는 z축으로의 유체의 속도(μ)의 기울기를 dμ/dz이라하면 유체의 점성 η는 다음과 같이 표현된다.

이상기체의 경우 유체에 작용하는 힘 F는 운동량의 변화에 단위시간에 충돌 수를 곱한 것과 같으므로 다음과 같이 표현된다.

여기서 N은 아보가드로의 수, υ는 분자 속도 m은 분자 질량에 해당한다. 그리고 λ는 단위거리 또는 분자의 평균행로에 해당된다. 따라서 η는

이것이 일반적으로 알려진 분자의 kinetic motion에 의한 이상기체의 점성이론이다.

실제기체나 dense gas의 경우는 intermolecular force의 영향을 받으므로 점성 η는 분자의 kinetic motion에의한 점성 ηK에 추가하여 intermolecular force에 의한 부분을 추가할 수 있다. 따라서 η는다음과 같이 표현된다.9

여기서 ηK는 kinetic motion에 의한 점성의 값으로써 kinetic pressure PK에 비례하며, ηI는 각각 intermolecular force에 의한 점성의 값으로써 internal pressure PI에 비례한다. 이상기체의 경우 식 (3)으로 표현된 점성 η는 ηK에 해당된다.

이 때 실제 기체의 점성 η는 다음과 같다.10

식 (5)에서의 점성의 크기는 절대압력(=kinetic pressure + internal pressure) 과 분자 충돌시간의 곱으로 표현되어 있다. 자동차의 brake의 제동력은 실린다의 마찰력과 제동시간의 곱에 비례하며이 점에서 식 (5)으로 표현된 점성의 크기는 자동차의 brake 세기와 유사하다하여 점성에 관한 brake model로 이름지은바 있다.

분자의 평균행로 λ와 분자속도 υ는 각각 다음과 같이 표현할 수 있으며

여기서 V, σ, k, T.P는 각각 몰부피, 경질구의 직경, 볼츠만의 상수, 절대온도 그리고 압력을 나타낸다. 한편 실제기체의 경우 이성분기체의 상태방정식으로서 다음과 같은 van der Waals 상태식을 이용할 수 있다.

여기서 R은 기체상수, Pm은 이성분기체의 압력, Vm은 이성분혼합기체의부피, a(x) b(x)는 몰분률 x의 변화에 따르는 혼합물의 van der Waals 상수이다. a(x), b(x)는 다음과 같이 표현된다.11

여기서

ζ는 상호작용인력에 관련되는 상수이다.

한편

으로 나타낸다.12

Aaa, Abb는 섞이기 전의 각성분의 인력을 나타내는 상수이며 Baa, Bbb는 섞이기 전의 순수한 상태의 부피를 나타내는 두 물질의 상수이다. (8) 식에서

Pm이 1기압일 때 Vm은 대략 다음과 같다.

식 (11-14)을 식 (15-16)에 대입하고 식 (15-17)을 식 (5)에 대입하면 이 성분기체의 점성을 계산할 수 있다.

 

계산결과 및 고찰

Sutherland8에 의하면 충돌직경

Table 1.Comparison of calculated and observed viscosities for nonpolar mixtures at 20 ℃ (1 atm)

여기서 Tb는 끓는 온도이다. 여기서 σo= (αb(x))1/3로 놓으면 된다. α값을 0.65, ζ값을 1로 놓고 비극성 기체 혼합물의 점성을 계산하여 실험치와 비교한 것이 Table 1에, 계산하기 비교적 어려운 극성혼합물의 점성값을 계산하여 Table 2에 표시하였다. 극성물질이건 비극성 물질이건 계산결과가 비교적 실험값과 잘 일치하였으므로 brake model은 낮은 압력하에서 이성분기체에 점성계산에 타탕한 모델임을 알 수 있으며 이로부터 dense gas 영역에서의 활용의 가능성이 있다고 판단된다. dense gas 영역에서 점성을 계산하는데는 van der Waals상태 방정식은 적합하지 않으므로 다음과 같이 dense gas영역에서 일반적으로 통용되는 Redlich-Kwong상태방정식을 사용하였다. 순수한 물질에서의 Redlich-Kwong상태 방정식은

Table 2.Comparison of calculated and observed viscosities for polar gas mixture at 1 atm. (T c; NH3 406 K, (C2H5)2O 467 K, H2S 373 K)

여기서 a', b'은 상수로서 다음과 같이 표현된다.14

식 (19-20)에서 Tc, Pc는 각각 임계온도, 임계압력을 나타낸다. 한편 혼합물의 상태방정식은

로 표현되며 여기서 Pm, Vm는 혼합물의 압력 및 부피이며, a(x), b(x)는 혼합물의 상수로서 다음과 같이 표현된다.

여기서 yi는 혼합물 각성분의 몰분률이며, Tm, Pm은 혼합물의 온도 및압력, Tr, Pr은 혼합물의 환산온도 및 환산압력으로 T/Tm, P/Pm으로 표현된다. 한편 혼합물의 임계온도 Tcm과 Pcm은 다음과 같다.14

기체 혼합물의 kinetic pressure Pk와 Internal pressure Pi는 각각 다음과 같으므로

마찬가지 방법으로 기체의 평균속도와 평균행로를 구하여 혼합기체의 점성 ηm을 나타내면 다음과 같다.

여기서 K는 상수, σm는 혼합물의 collision diameter, Mm은 혼합물의 질량이다. Sutherland correction을 사용하여 표시하면 σm=σmo(1 + 0.6 Tcm/T)1/2가 된다. α를 기하학적 상수라 할 때

로 표시된다. α에 0.536을 대입하고 식 (29)의 점성값의 단위를 μ-poise로 나타낸다면

로 나타낼 수 있다. 식 (31)을 이용하면 온도와 압력에 따르는 이성분기체를 비롯한 혼합기체 dense gas 의 점성값을 계산할수 있다. 높은 압력에서 H2-N2의 혼합기체에 대한 점성에 대해 계산한 결과 Table 3과 같다.

Table 3.Comparison of calculated and observed viscosities at high pressure for H2-N2 gas mixture at 20 ℃(Pc; H2 12.8 atm, N2 33.5 atm)

동일한 온도와 동일한 압력하에서 같은 몰수의 기체는 동일한 부피를 차지함으로 기체의 상태방정식이 성립될 수 있었다. 기체의 상태방정식의 특성으로 기체의 대응상태방정식을 구성할 수 있는 것처럼 기체의 점성이기체의 상태방정식내의 열역학적 성질로 구성되는 본점성이론은 모든 혼합기체의 임계점에서의 compressibility factor PcmVcm/(RTcm)값을 대략 0.292로 정하면 다음과 같이 점성에 관한 대응상태 방정식을 구성할 수 있다. 여기서 Vcm은 임계점에서의 부피를 의미한다.

여기서 ηc는 임계점에서의 혼합유체의 점성이며, ηr은 환산점성(Reduced viscosity)으로서 환산온도 Tr, 환산부피 Vr(=V/Vm), 환산압력 Pr로 표시되어 있다. 대응상태의 방정식은 정확한 정량적 계산은 할 수 없다 하더라도 정성적 판단과 상당한 정량적 분석을 통하여 여러 온도 및 압력 범위내에서 다양한 공업적 활용성을 기대할 수 있을 것이다. 원하는 온도 및 압력하에서 최저의 점성값이 무엇인가를 계산하면 초임계유체의 용질 추출에 활용이 기대된다.

 

결 론

액체 및 기체와 임계점에서의 점성의 계산 및 그 특성을 잘 설명했던 brake 점성이론을 이성분계의 기체에 적용한 결과 계산결과는 실험치와 잘 일치하였으며, 다음 특징처럼 식이 간단하고 계산이 용이하여 공업적 활용성이 뛰어남을 발견하였다.

(1) adjustable parameter가 없으므로 물질마다 parameter를 정하는 불편이 없다.

(2) 임계점에서의 점성값과 점성의 대응상태방정식을 구성할 수 있으므로 다양한 공업적 응용성이 기대된다.

(3) 낮은 압력은 물론 고압의 dense gas 점성 계산에 두루 활용할 수 있으며, 압력과 온도를 조정함으로써 밀도와 점성을 변화시킬 있는 혼합 초임계유체의 용매의 사용이 가능함으로서 용질추출의 기술를 제고시킬 수 있다.

본 연구는 2003년도 훙익대학교 교내연구비에 의해서수행되었으며, 연구에 도움을 주신 학교당국과 brake 점성이론을 개발하여 많은 논문을 쓰도록 해주신 고려대학교 화학과 최동식 교수님께 감사드립니다.

References

  1. Baker, J. A.; Fork, W.; Smith, F. Phys. Fluids 1964, 7,897. https://doi.org/10.1063/1.1711301
  2. Klein, M.; Smith, F. J. J. Res. Bur. Stand 1968, 72A,359. https://doi.org/10.6028/jres.072A.033
  3. Kim, W. S.; Chair, T. S.; Pak, H. Bull. Chem. Soc.1988, 9, 213.
  4. Kim, W. S.; Chair, T. S.; Pak, H. Bull. Chem. Soc.1989, 10, 372.
  5. Kim, W. S.; Chair, T. S. Korean J. Chem. Eng. 1992,10, 124. https://doi.org/10.1007/BF02697404
  6. Kim, W. S.; Kim, J. Y.; Chair, T. S. J. Kor. Chem. Soc.1992, 36, 376.
  7. Kestin, J.; Yata, J. J. Chem. Phys. 1968, 49, 4780. https://doi.org/10.1063/1.1669960
  8. Arun, K. P.; Dhattacharyya, P. K. J. Chem. Phys. 1969,51, 928.
  9. Hirshfielder, J. O.; Curtiss, C. F.; Bird, R. B. MolecularTheory of Gases and Liquids; Wiley New York; 1954.
  10. Akhlaq A. Pakistan J. Scientific Reserch 1980, 32, 44.
  11. Chapman, S.; Cowling T.G. The Mathematical Theoryof Nonuniform Gases, Cambridge University Press;New York; 1969.
  12. Kim, W.; Chair, T. S. J. Kor. Chem. Soc. 1995, 39, 891.
  13. Reid, R. C.; Prausnitz, J. M.; Sherwood, T. K. TheProperties of Gases and Liquids; McGraw-Hill BookCo; 3rd ed.; 197; 38.
  14. Neufeld, P. D.; Janzen, A. R.; Aziz, R. A. J. Chem.Phys. 1972, 57, 1100. https://doi.org/10.1063/1.1678363