تعیین روابط لازم برای پیش بینی فشار بخار مواد خالص و ترکیب های دوتایی آنها با استفاده از معادله حالت MMM

---

شرح پروژه

معادله حالت در مقاله پیوست ارائه شده. برای پیش بینی خواص تعادلی مثل فشار بخار و دانسیته فاز مایع و سایر نتایج محاسبه شده در جداول مقاله به روابط فوگاسیته نیاز است.

پیش بینی خواص تعادلی مثل فشار بخار رو با معادلات حالت PR و SRK پیاده کردم که شامل یک روش سعی و خطا است. برای معادله ارائه شده در این مقاله به نتیجه نرسیدم.

میخوام با معادله حالت این مقاله همون کاری رو انجام بدم که مثلاً با PR برای پیش بینی فشار بخار انجام میشه.
مراحل این طوری هست:
1- حدس فشار بخار از رابطه Lee kesler
2- محاسبه ضریب تراکم پذیری برای فاز مایع و بخار از معادله حالت (PR یا هر کدوم دیگه)
3- محاسبه فوگاسیته فاز مایع و بخار از رابطه متناسب با معادله حالت مرحله قبل
4- مقایسه فوگاسیته دو فاز و در صورتی که اختلاف معناداری داشته باشند اصلاح فشار بخار حدس در مرحله اول
5- تکرار محاسبات از مرحله دوم تا جایی که اختلاف فوگاسیته دو فاز اندک باشد و فشار بخار ثابت شود

برای این کار به چند رابطه نیاز هست:
1- فرم کلی معادله حالت این مقاله بر حسب ضریب تراکم پذیری.
که مثلاً برای پنگ رابینسون این هست:
z^3-(1-B)Z^2+(A-3B^2-2B)Z-(AB-B^2-B^3)=0
2- فرم کلی رابطه فوگاسیته بر حسب ضریب تراکم پذیری برای معادله حالت این مقاله. که برای پنگ رابینسون همچین چیزی هست:
ln(f/p)= Z-1-ln(Z-B)-(A/2.82B)*ln.......

دوتا رابطه بدست اوردم ولی توی متلب میزنم جواب نمیده انگار یه جایی اشتباهه
اگه رابطه رو بدست اوردید من کد متلبش رو برای SRK و PR نوشتم فقط کافیه به جای اون روابط جایگزین بشن تا ببینیم جواب میده یا نه

این پروژه شامل 1 فایل مهم است، لطفا قبل از ارسال پیشنهاد حتما نسبت به بررسی این فایل اقدام فرمایید.

مهارت ها و تخصص های مورد نیاز

---

---