펭로빈슨상태방정식으로 Pi를계산하려하는대요

1.계산의최종목적

탑,바텀,플래쉬단에서 Pi 액체,Pi 기체를계산하여 Ki값 계산후 기포점계산 sigma XiKi=1을만족할때까지

2.막힌부분

제가 탑최상단에서 액체를기준으로 Pi를계산할때, 최상단액체조성을써서

A,B값계산후 Z값의근을계산하려니

근3개중 1개는 1<Z<0밖,

1개는1<Z<0 안

1개는 Z' <0인음수가나옵니다..

왜음수가나오는지 AB값에문제가있는거같아서 한번 파일을검사해주셨으면합니다.

8시간째 책보고 반복계산을하였는데,

혼합물의 펭로빈슨계산법은 예제가찾기어려운데 어디물어볼곳이없어염치불구하고 남겨봅니다.

5성분증류.xls는 데이터시트이며,잘못된곳이있는지한번확인해주시면 감사하겠습니다.

저는맞다고생각하지만 개인적으로 ai*알파i 곱하는부분이먼가문제가있는게아닐까생각해봅니다

일반적으로 Z 3차식에서 A>B인거같은대 저는 A<B값으로나왔습니다.

물성및 저의계산결과는 PR만든것.xls파일에 있습니다.

계산과정에서 잘못된부분이 있다면 지적해주시면감사드리겠습니다.

그리고제가궁금한거 한가지가더있는대용

UNDERWOOD method에서1-q를 만족하는 파이값을찾는대

그과정에서 1-q, Q는포화액체=1이면 1-q=0이고

안토인으로계산시 너무작은최소단수가나왔습니다

문제에서는 플래쉬되는분율을 0.489 , Feed액체의상태는 포화액체로 주었는대요

이럴땐 언더우드법사용시 1-q=f=0.489를사용할지 아니면 1-q=0을 사용할지 질문드립니다.