Saya mencoba menghitung tekanan titik gelembung metanol dalam campuran metanol-isobutana pada 110 degC menggunakan parameter PC-SAFT. Kode tersebut memberikan kesalahan yang sama atau nilai tekanan yang berbeda secara acak. Apakah saya salah mengatur parameter k_ij atau apakah masalahnya terletak pada hal lain?

Kode:

import numpy as np
from pcsaft import flashTQ # documentation: https://pcsaft.readthedocs.io/en/latest/

'''
Data for methanol: https://pubs.acs.org/doi/full/10.1021/ie010954d
Data for isobutane: https://pubs.acs.org/doi/full/10.1021/ie0003887
'''
t = 110 + 273.15 # K
q = 0.2 # fraction, that is evaporated
x0 = 0.75 # mole fraction of methanol in mixture
x = np.asarray([x0, 1-x0])
pcsaft_params = {'m': np.asarray([1.5255, 2.2616]), # [methanol, isobutane]
                 's': np.asarray([3.2300, 3.7574]),
                 'e': np.asarray([188.90, 216.53]),
                 'e_assoc': np.asarray([0.035176]),
                 'vol_a': np.asarray([2899.5])}
pcsaft_params['k_ij'] = np.asarray([[0, 0],
                                    [0, 0]])

pvap, xl, xv = flashTQ(t, q, x, pcsaft_params)
print('Bubble point pressure: {} Pa'.format(pvap))
1
Sofi 10 Mei 2021, 03:31

1 menjawab

Jawaban Terbaik

Ada dua masalah:

  1. Anda perlu menambahkan nol untuk isobutana ke array parameter e_assoc dan vol_a
  2. Parameter vol_a dan e_assoc telah dibalik dan dicocokkan dengan kunci yang salah

Berikut adalah bagaimana parameter dict akan terlihat:

pcsaft_params = {'m': np.asarray([1.5255, 2.2616]), # [methanol, isobutane]
                 's': np.asarray([3.2300, 3.7574]),
                 'e': np.asarray([188.90, 216.53]),
                 'vol_a': np.asarray([0.035176, 0]),
                 'e_assoc': np.asarray([2899.5, 0])}

Dan izinkan saya menjelaskan mengapa nol ekstra itu diperlukan untuk isobutana. Jika salah satu dari parameter tambahan tersebut digunakan, maka nilai perlu diberikan untuk setiap komponen dalam sistem. Jika tidak, modul psaft tidak tahu komponen mana yang harus ditetapkan oleh parameter asosiasi. Untuk senyawa tanpa suku tambahan, gunakan saja nol.

0
Zach 12 Mei 2021, 05:22