Vapor Liquid Equilibrium (VLE) Binary System Formic acid-Acetic acid Using UNIQUAC (MATLAB)

#UNIQUAC #thermodynamics #VLE #matlab #activity
Vapor Liquid Equilibrium (VLE) Binary System Formic acid-Acetic acid Using UNIQUAC (MATLAB).
Focus on Activity Coefficient and Bubble Temperature Calculation at 1 atm.
Anderson, T.F., Prausnitz, J.M., 1978. Application of the UNIQUAC Equation to Calculation of Multicomponent Phase Equilibria. 1. Vapor-Liquid Equilibria. Ind. Eng. Chem. Proc. Des. Dev. 17, 552-561. doi.org/10.102...
Conti, J.J., Othmer, D.F., Gilmont, R., 1960. Composition of Vapors from Boiling Binary Solutions. Systems Containing Formic Acid, Acetic Acid, Water, and Chloroform. J. Chem. Eng. Data 5, 301-307. doi.org/10.102...
clc;clear;close all;
%Data of mole fraction of formic acid
x1d=[1.000;0.926;0.869;0.792;0.709;0.614;0.521;0.447;0.374;0.279;0.218;0.175;0.127;0.000];
y1d=[1.000;0.942;0.894;0.832;0.762;0.672;0.599;0.522;0.442;0.352;0.280;0.230;0.180;0.000];
%Data temp from degc to K
Td=[100.6;101.2;101.9;102.9;104;105.3;106.6;107.9;109.4;110.8;112.3;113.3;114.2;118.1]+273.15;
%input x1
n=19; %discreet
x1=linspace(1e-5,0.9999,n);
x1=x1';
x2=1-x1;
Pt=760; %mmHg Total Pressure in System
function P=Psat(T)
%Function for vapor press calculation
A1=[27.928 -2597.6 -7.2489 6.4110E-10 3.9421E-06];
A2=[28.376 -2973.4 -7.032 -1.5051E-09 2.1806E-06];
P1=10.^(A1(1)+A1(2)./T+A1(3).*log10(T)+A1(4).*T+A1(5).*T.^2);
P2=10.^(A2(1)+A2(2)./T+A2(3).*log10(T)+A2(4).*T+A2(5).*T.^2);
P=[P1 P2];
end

[@essidabouhligha3914]

Good job , keep going

[@rahmatsunarya]

Thank you Sir 👍😊

[@essidabouhligha3914]

@@rahmatsunarya even me I have some Matlab codes about thermodynamics modeling. I like thermodynamics modeling.

[@rahmatsunarya]

Perfect. Can you tell me your works in thermodynamic model?

[@essidabouhligha3914]

@@rahmatsunarya Activity Coefficients Cubic Equation of State. Especially in Cubic Equations of State & mixing rules.

[@essidabouhligha3914]

@@rahmatsunarya Two phase flash calculations using equations of state

[@infanatasyawalaraihan362]

Assalamualaikum, permisi bapak. Izin bertanya, apakah saya bisa mendapatkan script tersebut untuk keperluan studi saya dalam Termodinamika

[@rahmatsunarya]

Ok. Sudah dikirimkan via email

[@user-hx6oz3ws8m]

If we were interested in reading then why we come on u tube to listen lecture

[@rahmatsunarya]

I am so sorry, If you need more assistance, please let me know. thank you

[@velascomartinezeverdavid3660]

Sorry, could you share me your script, please?
I had problems with the funtion @fres

[@rahmatsunarya]

Surely, please get me in gmail 🙏🏻

[@fitrihakimisai2048]

What is z?

[@rahmatsunarya]

Hello. z is UNIQUAC coordination number

[@slimebox347]

Excuse me sir, perhaps can you upload this matlab script in some drive or web, cause i want to use it for some experienting with multiple component on vle

[@rahmatsunarya]

Please sent me email at [e-mail] . Thanks

[@slimebox347]

Thank you sir, i already send email to you