Peng Robinson równanie stanu daje lepsze przybliżenie krytycznego współczynnika ściśliwości , gdzie dane eksperymentalne waha się między a ( Źródło ). Ponieważ obliczona wartość powinna wynosić dla dowolnej substancji ( Źródło ), oczekiwałbym, że otrzymam raczej bliską liczbę, stosując ponownie numerycznie Równanie Stanu, biorąc pod uwagę punkt na diagramie Clapeyrona.
Jednak gdy próbuję rozwiązać równanie stanu dla , daje ono ( Mathematica ), rozwiązanie o rozbieżności od wartości teoretycznej.
Również gdy używam wyraźnie innych parametrów dla tego EoS z Wikipedii , wynik numeryczny z tym samym kodem wynosi .
thermodynamics
chemical-engineering
Reza Afzalan
źródło
źródło
Zc
jest niezależnym parametrem. proszę sprawdzić moje rozwiązanieOdpowiedzi:
Oryginalny papier Peng Robinson (1976) ma równanie ułożone sześciennie w stosunku do Z.
W punkcie krytycznym A i B wynoszą odpowiednio 0,45724 i 0,07780.
Kiedy równanie sześcienne zostanie rozwiązane dla jego trzech pierwiastków w punkcie krytycznym, daje jeden rzeczywisty pierwiastek i dwa wyobrażone pierwiastki. Z c = 0,30041 - 0,01199 i Z c = 0,30041 + 0,01199 i
Średnia z tych trzech pierwiastków wynosi:
W punkcie krytycznym w systemie jest tylko jedna faza. Jednak równanie sześcienne ocenia się na podstawie tego, jak dobrze jest w znalezieniu trzech pierwiastków odpowiadających każdej fazie w układzie. Zakładam, że standardową praktyką przy porównywaniu sześciennych równań stanu jest przyjmowanie średniej z wszystkich trzech pierwiastków, niezależnie od tego, czy są one rzeczywiste, czy nie.
źródło
Kiedy sześcienny EOS zostanie rozwiązany w punkcie krytycznym, MUSI mieć 3 równe pierwiastki, ponieważ równania parametrów zostały wyprowadzone w celu spełnienia ograniczeń krytycznych (punkt przegięcia z zerowym nachyleniem w punkcie krytycznym). Np. Dla RK Eos Zc = 1/3 (Giorgio Soave)
źródło
Krytyczny współczynnik ściśliwości Zc dla równania stanu Peng Robinsona jest podawany jako 0,307 w Walas, Phase Equilibria in Engineering Engineering. Współczynniki wynosiły 0,45724 i 0,07780. Za pomocą Excela znalazłem Zc = 0,3213. Stosując 0,457235 i 0,077796, jak podano w pracy Stryjeka i Very, Zc = 0,3111. Widać, że rozwiązanie jest bardzo wrażliwe na używane wartości. Nadal jednak nie mogę uzyskać 0,307.
źródło
Ważne jest, aby zrozumieć, że punkt krytyczny czystej substancji charakteryzuje się pewnymi właściwościami. Na przykład, pierwsza i druga pochodna cząstkowa ciśnienia w stosunku do objętości molowej w stałej temperaturze znikają. Po drugie, ważne jest, aby zrozumieć, że punkt krytyczny podany przez równanie stanu (EoS) niekoniecznie jest tym samym punktem, co określony na podstawie pomiarów.
Aby znaleźć punkt krytyczny EoS Peng-Robinsona (PR), można obliczyć pierwszą i drugą pochodną cząstkową w odniesieniu do objętości molowej i poszukać danej pary objętość / temperatura, w której znikają obie pochodne cząstkowe. Robiąc to za pomocą solvera dla układów równań nieliniowych (np. Metoda Levenberga-Marquardta), znalazłem punkt krytyczny CO2 zbliżony do i , co daje ciśnienie krytyczne , co daje krytyczną ściśliwość z . 0,89 ρ c 1,0 p c Z c = 0,30741.0Tc 0.89ρc 1.0pc Zc=0.3074
Aby uzyskać dobre wyniki, zalecam:
źródło