وبلاگ

اکنون معادلات (۶-۱۲۷) و (۶-۱۲۸) را ترکیب می کنیم و با کمی عملیات ریاضی کار را ادامه می دهیم.

مرحله ی آخر استفاده از مساوی بودن جمع کسرهای مولی هر فاز با یک است.

که این معادله شبیه به معادله ی Rachford-Rice در حل مساله جداسازی آنی هم دما است.
با داشتن تخمین مناسبی از K_LL، می توان از معادله ی (۶-۱۳۰) با روش حدس و خطای نیوتن ابتدا مقدار L_f و بعد x_i^aq. و x_i^HC را بدست آورد. برای تخمین اولیه ی K_LL از مقادیر در نظر گرفته شده برای کسر مولی C₁₄ و میزان تخمینی فاز هیدرو کربنی استفاده می کنیم.

این مقادیر انتخابی از آن جهت لازم است که می خواهیم دامنه ای از کسر مولی کلی (z_i) را بگونه ای در نظر بگیریم که اثرات میزان نسبی آب به گلیکول، مقدار نسبی C₁₄ به C₁₆ و اندازه ی فاز هیدروکربنی همه در قالب سه حلقه بدست آید.
حال که مقادیر بهتری برای کسر مولی ها بدست آمد به سراغ معادله ی NRTL^[74] می رویم و با توجه به پارامتر های گرفته شده از نرم افزار Aspen برای این معادله ی ضریب فعالیت، مقدار ضریب فعالیت (γ_i) را چه برای فاز آبی و چه برای فاز آلی می یابیم. همچنین، بدلیل اینکه تا اینجا فاز مایع هیدروکربنی تا حدودی اصلاح شده است، معادله ی تعادلی آن را با فاز جامد مورد توجه قرار می دهیم. برای این کار از مدل وان اصلاح شده^[۷۵] استفاده می شود. این معادله به ضریب فعالیت اجزا در فاز مایع و جامد نیازمند است که اولی از NRTL و دومی از فرم خاصی از UNIQUAC نوشته شده برای فاز جامد و اجزای هیدروکربنی محاسبه می شود.
در این فرم خاص UNIQUAC ثابت های r و q برای اجزای هیدروکربنی از دو معادله ی تجربی زیر برآورد می شود.

برای پارامتر τ_ij نیاز به گرمای نهان ویژه ی تصعید و عدد کووردیناسیون (Z) توصیف کننده ی ساختار بلوری فاز جامد داریم که اولی از ترکیب چند معادله ی تجربی که تابع دما، دمای بحرانی اجزا، عدد کربنی دو جزء (CN)و ضریب ω است بدست می آید. همچنین، با در نظر گرفتن ساختار بلوری اورتورومبیک^[۷۶] عدد Z مساوی ۶ اختیار می شود و برای کسر مولی دو جزء در فاز جامد نیز از  برای C₁₆و C₁₄استفاده می شود.
پس با تفاسیر بالا، ضریب فعالیت فاز جامد بدست می آید که همراه با ضریب فعالیت فاز مایع آلی در معادله ی مربوط به مدل پدرسن اصلاح شده قرار می گیرد که در زیر آمده است.

که اعداد A₄ و A₅ به ترتیب  و  هستند. همچنین، با محاسبه شدن سمت راست معادله ی (۶-۱۰۸) و اصلاح شدن کسر مولی فاز آلی با معادله ی (۶-۱۲۹) مقدار اصلاح شده ی کسر مولی جامد هم محاسبه می گردد.
تا این مرحله یک بار حلقه ی حدس و خطا را اجرا کرده ایم. تنها کافی است که با مقادیر کنونی x_i^HC و x_i^aq. مقدار جدید K_LL را از نسبت ضرایب فعالیت بیابیم و حلقه ی حدس و خطا را تکرار نماییم (یعنی از معادله ی ۱۳۰).