Một số vấn đề thường gặp trong quá trình vận hành- 123docz.net

HNO 3 60% BẰNG PHẦN MỀM ASPEN HYSYS

3.6 Một số vấn đề thường gặp trong quá trình vận hành và cách xử l ý

Như phần tổng quan đã trình bày, một số sự cố thường gặp trong vận hành có thể mô phỏng bằng phần mềm là sự suy giảm hoạt tính xúc tác và sự thay đổi chất lượng nguyên liệu đầu vào, trong đó sự suy giảm hoạt tính xúc tác là khó xác định do đó trong luận văn chỉ nghiên cứu về ảnh hưởng của hàm lượng nước trong NH 3 nguyên liệu đến nhiệt độ lò đốt K101.

Do hiện tượng cân bằng pha, tại thiết bị hóa hơi NH 3 - E103, với mỗi giá trị nhiệt độ và áp suất có một giá trị xác định hơi nước bão hòa, tức là khi hàm lượng nước trong pha lỏng tăng thì hàm lượng nước trong pha hơi sẽ tăng nhưng đến giá trị bão hòa sẽ không tăng nữa.

Với điều kiện vận hành của dây chuyền, để đảm bảo hệ thống hoạt động ổn định, người ta khống chế nhiệt độ K101 tại giá trị 910±10oC, đồng thời giới hạn nhiệt độ lớn nhất tại E103 là 87oC (để tránh hiện tượng cuốn theo pha lỏng vào pha hơi làm nhiệt độ K101 bất ổn định ) [28].

Theo điều kiện trên, để khảo sát được tập chung, luận văn sẽ khảo sát sự phụ thuộc nhiệt độ K101 vào hàm lượng nước trong NH3 khi nhiệt độ tại E103 là 87oC.

Thiết lập sơ đồ mô phỏng như sơ đồ rút gọn nhưng chỉ đến thiết bị K101 đồng thời bổ xung một dòng vật liệu là nước vào dòng NH3 ban đầu. Kết quả khảo sát như sau: (xem phụ lục 4)

B ả n g 3 .2 0 Q u a n h ệ g iữ a n h iệ t đ ộ K 1 0 1 v à h à m lư ợ n g n ư ớ c tr o n g N H 3

Hàm lượng nước, % Nhiệt độ

K 1 O 1

( o C)

Ghi chú Dòng NH 3 trước

E103

Dòng NH3 sau E103

0 , 0 0 0 , 0 0 910,3

0,30 0,30 908,5

0,60 0,60 906,6

0,89 0,89 904,8

1,19 1,19 903,0

1,48 1,48 901,2

1,77 1,77 899,4

2,06 2,06 897,6

2,34 2,30 896,1

2,63 2,30 896,1

2,91 2,30 896,1

3,19 2,30 896,1

3,47 2,30 896,1

3,75 2,30 896,1

4,03 2,30 896,1

4,31 2,30 896,1

4,58 2,30 896,1

4,85 2,30 896,1

5,12 2,30 896,1

5,39 2,30 896,1

5,66 2,30 896,1

Mối quan hệ thể hiện dưới dạng đồ thị như hình 3.6.

Hàm lượng nước trong NH3 ban đầu, % khối lượng

H ìn h 3 .6 Q u a n h ệ n h iệ t đ ộ K 1 0 1 v à h à m lư ợ n g n ư ớ c tro n g N H 3

Kết quả trên cho thấy, để hạn chế ảnh hưởng của hàm lượng nước tới nhiệt độ K101 (theo giới hạn >900oC), yêu cầu hàm lượng nước không được lớn hơn 1,77%.

Qua đồ thị cũng cho thấy giới hạn Max của hàm lượng nước trong NH 3 ban đầu (theo điều kiện nhiệt độ tại E103 không quá 87oC) là 2,34% (khi vượt qua giới hạn này thì mức trong E103 sẽ luôn tăng).

KẾT LUẬN

Sau một thời gian nghiên cứu, khảo sát, tiến hành thực hiện, luận văn đã đạt mục đích, yêu cầu đề ra.

1. Đã xây dựng được mô hình mô phỏng tổng thể dây chuyền sản xuất axit HNO 3 60% loại đơn áp. Xây dựng mô hình động học quá trình sản xuất axit HNO 3 bằng công nghệ oxi hóa NH 3 với 5 bước phản ứng chính gồm:

- Phản ứng oxi hóa NH 3 thành NO với lưới xúc tác 90% Pt và 10% Rh - Phản ứng oxi hóa NO thành NO 2

- Phản ứng dime hóa NO 2 thành N 2 O 4

- Phản ứng và cân bằng hấp thụ NO 2 vào H 2 O thành HNO 3 và NO - Phản ứng trong tháp tẩy trắng.

2. Đã khảo sát và lựa chọn mô hình chuyển hóa NH 3 thành NO phù hợp với công nghệ của công ty hóa chất 95, kết quả mô phỏng là hợp lý.

3. Phương pháp nghiên cứu hiện đại, ứng dụng phần mềm mô phỏng Aspen Hysys mô phỏng tĩnh phân xưởng sản xuất HNO 3 60%, thiết lập các phản ứng, thiết lập thành phần các cấu tử cũng như dòng nguyên liệu vào hệ thống, thiết lập các thông số thiết bị.

4. So sánh kết quả mô phỏng với số liệu thiết kế với các chỉ tiêu: tỷ lệ thành phần dòng nguyên liệu, nhiệt độ, tỷ lệ thành phần nguyên liệu và độ chuyển hóa, tính toán cân bằng vật chất và năng lượng, tối ưu hóa. Kết quả mô phỏng và số liệu thiết kế là tương đương nhau.

5. Đã đánh giá việc sử dụng năng lượng hiện tại và hướng cải tiến theo công nghệ tối ưu. Đã khảo sát yếu tố gây sự cố mất ổn định nhiệt độ K101 và đưa ra yêu cầu về hàm lượng nước trong NH 3 nguyên liệu <1,77% (giá trị giới hạn là 2,34).

TÀI LIỆU THAM KHẢO

1. Gutierrez-Canas, C, Arias, PL, and Legarreta, JA (1989), "Industrial nitrogen oxides absorption simulation", C o m p u t e r s & c h e m ic a l e n g in e e r in g . 13(9), pp.

985-1002.

2. Lee, Hyo C and Farrauto, Robert J (1989), "Catalyst deactivation due to transient behavior in nitric acid production", In d u s t r ia l & e n g in e e r in g c h e m is t r y r e s e a r c h .

28(1), pp. 1-5.

3. Liu, Ji Dong, et al. (2013), Simulation and Optimization of Absorption Column of the Dual-Pressure Nitric Acid Production Process, A d v a n c e d M a t e r ia ls R e s e a r c h , Trans Tech Publ, pp. 3830-3834.

4. Mummey, MJ and Schmidt, LD (1981), "Decomposition of NO on clean Pt near atmospheric pressures: I. Steady state kinetics", S u rfa c e S c ie n c e . 109(1), pp. 29-42.

5. Sperner, Franz and Hohmann, Wolfgang (1976), "Rhodium-platinum gauzes for ammonia oxidation " , P la t in u m M e t a ls R e v ie w . 20(1), pp. 12-20.

6. Takoudis, CG and Schmidt, LD (1983), "Reaction between nitric oxide and ammonia on polycrystalline platinum. 1. Steady-state kinetics", T h e J o u r n a l o f P h y s ic a l C h e m is t r y . 87(6), pp. 958-963.

7. Vernikovskaya, NV, Pinaeva, LG, and Isupova, LA (2014), "Oxidation of ammonia to NOx in a two bed reactor (Pt gauzes+ oxide monolytic layer):

Experimental studies and mathematical modelling " , C h e m ic a l E n g in e e r in g J o u r n a l. 238, pp. 140-147.

8. Abdel Salam, Mohamed (2016), "Modelling and Simulation of Gauze Reactor of Ammonia Oxidation " , A m e r ic a n J o u r n a l o f C h e m ic a l E n g in e e r in g . 4(1), p. 16.

9. Belghaieb, J. , et al. (2010), "Energy optimization of a network of exchangers- reactors in a nitric acid production plant " , C h e m ic a l E n g in e e r in g T r a n s a c t io n s

21, pp. 271-276.

10. Dalaouti, Natassa and Seferlis, Panos (2005), "Design sensitivity of reactive absorption units for improved dynamic performance and cleaner production: the NOx removal process", J o u r n a l o f C l e a n e r P r o d u c t io n . 13(15), pp. 1461-1470.

11. Emig, G, Wohlfahrt, K, and Hoffmann, U (1979), "Absorption with simultaneous complex reactions in both phases, demonstrated by the modeling and calculation of a counter-current flow column for the production of nitric acid", C o m p u t e r s &

C h e m ic a l E n g in e e r in g . 3(1-4), pp. 143-150.

12. Fila, Vlastimil and Bernauer, Bohumil (1994), "A mathematical model of a gauze reactor for the ammonia oxidation", C o lle c t io n o f C z e c h o s lo v a k c h e m ic a l c o m m u n ic a t io n s . 59(4), pp. 855-874.

13. Gutierrez-Caras, C. , Arias, P. L., and Legarreta, J. A. (1989), "Industrial nitrogen oxides absorption simulation gutierrezcaas", C o m p u t e r s rh em . E n g n g . 13(9) pp.

985-1002.

14. Hamid, Mohd Kamaruddin Abd (2007), "HYSYS®: An Introduction to Chemical Engineering Simulation", A p o s t ila d e H a m id .

15. Hickman, Daniel A. and Schmidt, Lanny D. (1990), "Modeling Catalytic Gauze Reactors: Ammonia Oxidation", In d . E n g . C h e m . R e s . 30, pp. 50-55.

16. Pignet, T. and Schmidt, L. D. (1974), "Selectivity of nh3 oxidation on platinum",

C h e m ic a l E n g in e e r in g S c ie n c e . 29 pp. 1123-1131.

17. Pignet, T. and Schmidt, L. D. (1975), "Kinetics of NH3 Oxidation on Pt, Rh, and Pd", J o u r n a l o f c a t a ly s is . 40, pp. 212-225.

18. Savenkov, A. S., et al. (2016), "Modeling of ammonia oxidation on a platinoid catalyst, taking into account the N2O formation " , R u s s ia n J o u r n a l o f A p p l ie d C h e m is t r y . 88(10), pp. 1570-1575.

19. Suchak, N. J. and Joshi, J. B. (1994), "Simulation and Optimization of NO, Absorption System in Nitric Acid Manufacture ", A l C h E J o u r n a l. 40(6).

20. Sweeney, A.J. and Liu, Y. A. (2011), "Use of Simulation To Optimize NOx Abatement by Absorption and Selective Catalytic Reduction " , In d . E n g . C h e m . R e s . 40, pp. 2618-2627.

21. Thiemann, Michael, Scheibler, Erich, and Wiegand, Karl Wilhelm (2000),

"Nitric acid, nitrous acid, and nitrogen oxides " , U llm a n n 's e n c y c lo p e d ia o f in d u s t r ia l c h e m is t r y .

22. Van Nieuwenhuyse, Ave E (2000), "Production of Nitric acid", E u r o p e a n F e r t iliz e r M a n u fa c t u r e r s ' A s s o c ia t io n ( E F M A ) , B r u s s e ls B e lg iu m . Booklet No. 2 of 8.

23. Verheijen, Dr.Ir. P.l.T. and De Leeuw, Ir. V. (1992), M o d e lin g a n d s t o c h a s t ic s im u la t io n o f the r e a c t io n s a n d a b s o r p t io n o f n o x in w a te r, Michiel Spoor Laveibos Zoetenneer The Netherlands.

24. Wiegand, Karl Wilhelm, Scheibler, Erich, and Thiemann, Michael (1990),

"Computation of plate columns for NOx absorption by a new stage-to-stage method", C h e m ic a l E n g in e e r in g & T e c h n o lo g y . 13(1), pp. 289-297.

25. Aipom£HKO, B.H and KaprHH, C.H (1962), TexH on oeuH Á 3 0 m H 0 ũ K u c n o m u ,

MocKOBa.

26. BopõbeB, HHKonan HBaHOBHH (2011), TexH on oeuH c8H 3aH Hoeo a 3 o m a u

a.3om Hbix y ò o ố p e H U ũ , EenopyccKHH rocygapcTBeHHbiH TexHonornnecKHH

yHHBepcmeT.

27. ronoBHa, EB (2008), O K u e n e H u e a M M u a K a H a m a m u H O u Ò H b ix c e m K a x u ốnoHHOM OKcuỜHOM K a m a n m a m o p e c o m o e o ũ c m p y K m y p u , M.: p x y HM.

MeHgeaeeBa. 2009.

28. Tổng cục công nghiệp quốc phòng -Z195 (2003), "Quy trình công nghệ sản xuất axit Nitric 60% ", p. 44.

29. Nguyễn Thị Minh Hiền (2014), M ô p h ỏ n g c ô n g n g h ệ h ó a h ọ c , Nhóm Mô phỏng Công nghệ Hóa học và Dầu khí Trường Đại học Bách khoa Hà Nội.

30. Tổng cục công nghiệp quốc phòng-Z195 (2012), "TCCS 01:2012/CT95 - Axit Nitric - Yêu cầu kỹ thuật", p. 16.

31. Nguyễn Toàn Quyền (2014), C ơ s ở lý th u yế t v à ứ n g d ụ n g c ủ a p h ư ơ n g p h á p p h â n tíc h E x e r g y đ ể đ á n h g iá h iệ u q u ả c ủ a c á c h ệ th ố n g n h iệ t -lạ n h , Trường Đại học Bách Khoa Hà Nội.

32. Nguyễn Hoa Toàn và Lê Thị Mai Hương (2005), C ô n g n g h ệ c á c h ợ p c h ấ t v ô c ơ c ủ a N itơ , Nhà xuất bản khoa học và kỹ thuật.

PHỤ LỤC 1

1. Phản ứng ôxi hóa NH 3 Biểu thức mô hình 1:

- Biểu thức tốc độ phản ứng theo mô hình 1 có đơn vị (molcule.cm '2 .s -1 ): số phân tử phản ứng trên 1 đơn vị diện tích bề mặt xúc tác trong 1 giây, cần phải chuyển thành dạng đơn vị kmol.m "3 .s -1 .

- Theo [15], cần nhân tử số của biểu thức tốc độ (ở mô hình 1, chương 1) với A/(V.N o ):

A/V là số diện tích bề mặt có trong 1 đơn vị thể tích xúc tác; N o là số Avogadro=6,022E+23;

- Thông số lưới súc tác 90%Pt, 10%Rh:

+ Độ tinh khiết của Platin và Rhodi:

+ Khối lượng riêng của súc tác:

+ Khối lượng trên 1 m2 tiết diện lưới (g/m2):

+ Diện tích bề mặt trên đơn vị khối lượng:

+ Đường kính sợi Platin d w (mm):

+ Đường kính ngoài của lưới (mm):

+ Số lưới sử dụng (n L ):

+ Đường kính hữu dụng của lưới:

+ Tiết diện ngang lưới (S =1,85 2 x3,14/4):

> 99,95 %

19450 kg/m3 (19,45g/cm3) 630 ± 5 %

2,65 m2/kg 0,076 ± 0,003

1910 ± 5 12

1850mm 2,688 m2

+ Khối lượng hữu dụng của súc tác (m=S.n L .0,630): 20,3215 kg + Diện tích bề mặt của chất xúc tác trên một đơn vị thể tích khí:

A/V = A /(£V k ) = 4,26.104 m2/m3 = 4,26.108 cm2/m3 + Chiều dầy 12 lớp lưới: L = 12.d w :

- Theo [12], A, V xác định như sau:

+ A: Diện tích bề mặt trên 1 m 2 tiết diện lưới

9,12.10"4 m

A = 2n.n2w.dw.lw = 1,57 m2 + V: Thể tích lưới trên 1m2 tiết diện lưới V=£V k

V k = n.n2w.d2w.lw = 5,98.10"5m3 l w = [d 2w + (1/n w ) 2 ] 1/2

+ n w : số lỗ trên 1 đơn vị dài của lưới, n w = 80 lỗ/in = 32 lỗ/cm = 3200 lỗ/m + d w : đường kính lưới xúc tác = 0,076.10 "3 m

+ Thể tích tự do: f = 1-0,5nd w (n w ) 0,5 =61,8%

Như vậy, biểu thức tốc độ phản ứng theo mô hình 1 là:

Phản ứng Tốc độ, (kmol/m 3 .s);

P, torr; T, K NH3 —> 1/2N2 + 3/2H2 6,98.10 3.exp(-46^^).PNH3

1+9,2.10- 5 .exp(l 2682 ).P WH3 NH 3 + 5 / 4 O 2 —>NO+3/2H 2 O 1,48.101.exp (10Ệ 0).PN h 3.P o /22

1 + 1 , 6.10 3 .exp( 2^3 ).PW tf 3 NH 3 + 3 / 2 NO—>5/4N 2 +3/2H 2 O 1,05.102.exp(4^ ).P wo.P^ /23

[1+5.1O- S.exp( 8 0 02 ).Pwo + ũ,O145.exp(^^9,5).Pí^/í23]2 NO —> 1 / 2 N 2 + 1 / 2 O 2 3,91.101.exp ( - 2625).PN0

1 + 6 , 95 . 10 - 4 .exp( 34 F 7 ).Pwo + 1,56.exp(4775).Po2 H 2 + 1 / 2 O 2 --- > H 2 O 1,06.104 Po 2 .P „2

Biểu thức mô hình 2:

- Đơn vị hằng số A của biểu thức hằng số tốc độ là Nm3/(kg.s.atm): thể tích quy về điều kiện tiêu chuẩn của khí phản ứng trên đơn vị khối lượng xúc tác trên giây và áp suất. Để chuyển về đơn vị kmol/(m3.s.atm), tử số chuyển về số kmol (1kmol = 22,4Nm3, mẫu số cần chuyển đơn vị khối lượng xúc tác thành thể tích, hay phải nhân thêm giá trị mK/V =mK/(Vk.f) =1/(5,98.10-5.0,6 1 8.2,65) = 1,0211.104 (kg/m3)

(mK: khối lượng xúc tác trên 1 đơn vị diện tích bề mặt =1/2,65 kg/m2) Biểu thức hằng số k theo mô hình 2 là:

Phản ứng Biểu thức tốc độ A E

(kJ/kmole) NH3+1,25O2 ^ NO + 1,5H2O r 1 = k1.PNH3 1,66.10 5

kmol/(m 3 .s.atm) 73,6 NH3+0,75O2 ^0,5N2+1,5H2O r 2 = k2.PNH3 3,87.10 3

kmol/(m 3 .s.atm) -28,4

NO ^ 0,5N2 + 0,5O2 r3 = k3.P2NO 8,21

kmol/(m 3 .s.atm 2 ) 3,76 2. Phản ứng oxi hóa NO

2NO + O2 ằ 2 NO2 7/ = ệ f . V no . Vo2, [kmo1 m-3 s-l]

rr = . K . p ị o 2, [kmol m-3 s-l]

log 10 ^ = 652,1/T - 1,0366; [atm-2 s -1],

P n o .P o 2 5749 lgK = lg ^ 2 ^

P n o 2 T + l,7 5 lg T — 0,0005T + 2,839, a tm Chuyển sang dang biểu thức r f = - . p ị o . P o ? , [kmol m-3 s -1]

Tính lại 2 giá trị hằng số k -k /R tại 2 nhiệt độ khác nhau, sau đó xác định lại k o , E Ta được biểu thức sau (R= 0,08205 lit.atm.mol -1 .K -1 ; E đơn vị kJ/kmol):

_ _ 1,1203.exp ( - _

rf = T — . p 2N 0 .P o 2 , [kmol m-3 s -1]

r r =

7,89.107 .exp( 10 r 2 t 681) 2

— P n o 2 , [kmol m-3 s-l]

3. Phản ứng đime hóa

2NO 2

N 2 O 4

r = — . ( P 2N 0 2 — ^ Ị 2r ± )

k= 5,7.10 5 atm -1 .s -1

lgK = lg(p N2O4 /p 2NO2 ) = 2993/T - 9,226; atm-1

Tính lại k tương tự như phản ứng ôxi hóa NO, ta được các biểu thức:

6,95.106 2

T f = . P n o 2 , [kmol m-3 s-l]

,16 —57297}

r r = _ 1,17.10 16.exp {-57297)

T P n 204 , [kmol m-3 s-l]

Hoặc có thể sử dụng biểu thức K:

P2 no 2692 ,

LgK = lg[ ! ^ ) = — + 1,75lgT + 0,00484.T — 7,144.10- 6 . T 2 + 3,062

,P n 20. 2U4 ' T

Chuyển sang dạng Ln:

L n K -1 = i J Pn2°4 P2 NŨ2

6198,6

T — 1,75lnT — 0 ,0 1 1 1 4 .T + 1,645.10-5 . T 2 — 7,051

4. Phản ứng hấp thụ NO 2 hình thành axit nitric

3 NO2 + H 2O ^ 2 H N O 3 + N O ; K = [p2 HNO3 .pNo ].[p3 NO2 .pH2o ]"1 lg K = 7 0 7 / T - 4,71; (atm-1)

Hoặc K2 = 1,307.10'9.exp(4376/T) Chuyển sang dạng Ln :

( P h n o ? . P n o \ 1628

L n K = l n ( 2. 0 --- ) = --- --- 10,85, a t m 1

\ p ' no 2 - P h 2 o J T

5. Phản ứng pha lỏng trong tháp tẩy trắng Áp dụng mô hình

R ea ction S to ich iom etry R a te ex p ression

R5.1 2 NO 2 + H 2 O ^ HNO 2 + HNO 3 r5.1f = k 5.1. [NŨ2]2

R5.2 N 2 O 4 + H 2 O ^ HNO 3 + HNO 2 r5.2f = k 5 . 2 . [^ 2 0 ị ]

R5.3 3 HNO 2 ^ HNO 3 + H 2 O + 2NO , [HNO 2 ] 4

r5.3f = k 5.3. 2

______ E m ____

Giá trị hằng số k cho theo bảng sau:

k 5.1 20 4,67209 m3.kmol-1.s-1

k 5.2 jq ((-4139/T)+16,3415) s-1

k 5.3 Ị q ((-6200/T)+20,1979) atm2.m9.kmol-3.s-1

Tương tự trên, chuyển biểu thức k về dạng phương trình Arrenius, riêng phản ứng 5.3, trước khi chuyển ta đưa biểu thức tốc độ có áp suất NO về dạng nồng độ bằng việc thay [NO] = H no .P no . ; hằng số Henry H no = 1,84x10"5.e(1500((1/T) "1/298)) (Kmol.m"3.kPa"1) Ta được các biểu thức:

rs. 1 f = 46999. [ N0 2]2 T 5 . 2 f = 2 , 2 . 1 0 16 . exp ( ( 19235\

.[N 2 O 4 ] K RT )

r 5 . 3 f = 2,33.1016. exp ( ' 93145\ [HNO 2 ] . RT ) . [NO ] 2

PHỤ LỤC 2

1 2 3 4 5

L E G E N D S B u r lin g to n , M A U S A

C a s e N a m e : S O D O P F D D A Y D U . H S C

U n it S e t: SI

D a te /T im e : W e d A p r 1 8 0 9 : 2 4 : 3 8 2 0 1 8

6 7 8

F lu id P a c k a g e : S R K

Material Stream: 1.3

P r o p e r t y P a c k a g e : S R K - T w u

9 10

C O N D I T I O N S

11 O v e r a ll V a p o u r P h a s e

12 V a p o u r / P h a s e F r a c tio n 1 .0 0 0 0 1 .0 0 0 0

13 T e m p e r a t u r e : (C ) 1 2 2 .0 1 2 2 .0

14 P r e s s u r e : (k P a ) 7 5 0 . 0 7 5 0 . 0

15 M o la r F lo w (k g m o le /h ) 3 3 1 . 5 3 3 1 . 5

16 M a s s F lo w (k g /h ) 9 5 2 6 9 5 2 6

17 S td Id e a l L iq V o l F lo w (m 3 /h ) 1 0 .9 5 1 0 .9 5

18 M o la r E n t h a lp y ( k J /k g m o le ) - 1 7 5 6 - 1 7 5 6

19 M o la r E n t r o p y ( k J /k g m o le - C ) 1 4 4 .6 1 4 4 .6

2 0 H e a t F lo w (k J /h ) - 5 . 8 2 0 e + 0 0 5 - 5 . 8 2 0 e + 0 0 5

21 L iq V o l F lo w @ S t d C o n d (m 3 /h ) 7 8 3 5 7 8 3 5

2 2 2 3

C O M P O S I T I O N

2 4 2 5

O v e r a l l P h a s e V a p o u r F r a c tio n 1 .0 0 0 0

2 6 2 7

C O M P O N E N T S M O L A R F L O W

(k g m o le /h )

M O L E F R A C T I O N M A S S F L O W (k g /h )

M A S S F R A C T I O N L IQ U ID V O L U M E F L O W (m 3 /h )

L IQ U I D V O L U M E F R A C T I O N

2 8 A m m o n ia 0 .0 0 0 0 0 . 0 0 0 0 0 .0 0 0 0 0 . 0 0 0 0 0 .0 0 0 0 0 . 0 0 0 0

2 9 N itr o g e n 2 5 4 . 7 4 6 9 0 . 7 6 8 5 7 1 3 6 . 2 2 6 2 0 .7 4 9 1 8 .8 4 9 8 0 . 8 0 8 4

3 0 O x y g e n 6 7 . 5 7 7 9 0 . 2 0 3 9 2 1 6 2 . 4 9 2 8 0 . 2 2 7 0 1 .9 0 0 8 0 . 1 7 3 6

31 A r g o n 2 .8 6 1 6 0 . 0 0 8 6 1 1 4 . 3 1 6 8 0 . 0 1 2 0 0 .0 8 3 4 0 . 0 0 7 6

3 2 H 2 O 6 .2 9 2 7 0 . 0 1 9 0 1 1 3 . 3 6 4 2 0 . 0 1 1 9 0 .1 1 3 6 0 . 0 1 0 4

3 3 H y d r o g e n 0 .0 0 0 0 0 . 0 0 0 0 0 .0 0 0 0 0 . 0 0 0 0 0 .0 0 0 0 0 . 0 0 0 0

3 4 N O 0 .0 0 0 0 0 . 0 0 0 0 0 .0 0 0 0 0 . 0 0 0 0 0 .0 0 0 0 0 . 0 0 0 0

3 5 N O 2 0 .0 0 0 0 0 . 0 0 0 0 0 .0 0 0 0 0 . 0 0 0 0 0 .0 0 0 0 0 . 0 0 0 0

3 6 N 2 O 4 0 .0 0 0 0 0 . 0 0 0 0 0 .0 0 0 0 0 . 0 0 0 0 0 .0 0 0 0 0 . 0 0 0 0

3 7 H N O 3 0 .0 0 0 0 0 . 0 0 0 0 0 .0 0 0 0 0 . 0 0 0 0 0 .0 0 0 0 0 . 0 0 0 0

3 8 N 2 O 3 * 0 .0 0 0 0 0 . 0 0 0 0 0 .0 0 0 0 0 . 0 0 0 0 0 .0 0 0 0 0 . 0 0 0 0

3 9 H N O 2 * 0 .0 0 0 0 0 . 0 0 0 0 0 .0 0 0 0 0 . 0 0 0 0 0 .0 0 0 0 0 . 0 0 0 0

4 0 T o ta l 3 3 1 . 4 7 9 2 1 .0 0 0 0 9 5 2 6 . 4 0 0 0 1 .0 0 0 0 1 0 .9 4 7 6 1 .0 0 0 0

41 4 2

V a p o u r P h a s e P h a s e F r a c tio n 1 .0 0 0

4 3 4 4

C O M P O N E N T S M O L A R F L O W

(k g m o le /h )

M O L E F R A C T I O N M A S S F L O W (k g /h )

M A S S F R A C T I O N L IQ U ID V O L U M E F L O W (m 3 /h )

L IQ U I D V O L U M E F R A C T I O N

4 5 A m m o n ia 0 .0 0 0 0 0 . 0 0 0 0 0 .0 0 0 0 0 . 0 0 0 0 0 .0 0 0 0 0 . 0 0 0 0

4 6 N itr o g e n 2 5 4 . 7 4 6 9 0 . 7 6 8 5 7 1 3 6 . 2 2 6 2 0 .7 4 9 1 8 .8 4 9 8 0 . 8 0 8 4

4 7 O x y g e n 6 7 . 5 7 7 9 0 . 2 0 3 9 2 1 6 2 . 4 9 2 8 0 . 2 2 7 0 1 .9 0 0 8 0 . 1 7 3 6

4 8 A r g o n 2 .8 6 1 6 0 . 0 0 8 6 1 1 4 . 3 1 6 8 0 . 0 1 2 0 0 .0 8 3 4 0 . 0 0 7 6

4 9 H 2 O 6 .2 9 2 7 0 . 0 1 9 0 1 1 3 . 3 6 4 2 0 . 0 1 1 9 0 .1 1 3 6 0 . 0 1 0 4

5 0 H y d r o g e n 0 .0 0 0 0 0 . 0 0 0 0 0 .0 0 0 0 0 . 0 0 0 0 0 .0 0 0 0 0 . 0 0 0 0

51 N O 0 .0 0 0 0 0 . 0 0 0 0 0 .0 0 0 0 0 . 0 0 0 0 0 .0 0 0 0 0 . 0 0 0 0

5 2 N O 2 0 .0 0 0 0 0 . 0 0 0 0 0 .0 0 0 0 0 . 0 0 0 0 0 .0 0 0 0 0 . 0 0 0 0

5 3 N 2 O 4 0 .0 0 0 0 0 . 0 0 0 0 0 .0 0 0 0 0 . 0 0 0 0 0 .0 0 0 0 0 . 0 0 0 0

5 4 H N O 3 0 .0 0 0 0 0 . 0 0 0 0 0 .0 0 0 0 0 . 0 0 0 0 0 .0 0 0 0 0 . 0 0 0 0

5 5 N 2 O 3 * 0 .0 0 0 0 0 . 0 0 0 0 0 .0 0 0 0 0 . 0 0 0 0 0 .0 0 0 0 0 . 0 0 0 0

5 6 H N O 2 * 0 .0 0 0 0 0 . 0 0 0 0 0 .0 0 0 0 0 . 0 0 0 0 0 .0 0 0 0 0 . 0 0 0 0

5 7 T o ta l 3 3 1 . 4 7 9 2 1 .0 0 0 0 9 5 2 6 . 4 0 0 0 1 .0 0 0 0 1 0 .9 4 7 6 1 .0 0 0 0

5 8 5 9 6 0 61 6 2 6 3 6 4 6 5 6 6