برنامه متلب برای رسم نمودار T-xy مخلوط های دو جزیی پیروی کننده از قانون رائولت

در این قسمت می خواهیم با استفاده از نرم افزار متلب کدی را بنویسیم که نمودار های T-xy و x-y را برای مخلوط های دو جزیی که از قانون رائولت پیروی می کنند را رسم نماید.

برای این کار ابتدا نرم افزار باید مقدار فشار را از کاربر دریافت کرده و در آن فشار نقطه جوش دو ماده ی خالص را محاسبه نماید و سپس در بازه ی بین این دو دما مقادیر x و y را محاسبه نماید.برای محاسبه ی نقطه جوش می توان از معادله ی آنتوان استفاده نمود.بنابراین کدی که خواهیم نوشت باید ضرائب A , B , C معادله ی آنتوان را نیز برای هر کدام از دو ماده از کاربر سوال کند.

Log10(Psat)=A-B/(T+C)

برای محاسبه ی x:

Pt=x1P1sat+(1-x1)P2sat

x=(Pt-P2sat)/(P1sat-P2sat)

دقت کنید که جز ۱ جز فرارتر در سیستم است.

برای محاسبه ی y:

yPt=xP1sat

y=xP1sat/Pt

کد نوشته شده برای برنامه را در زیر مشاهده می کنید.

۱٫ % This program plot the T-xy @ x-y diagrams for mixtures that obey from Raoult’s law.
۲٫ %———————————————————————————–
۳٫ clc
۴٫ clear all
۵٫
۶٫ A1=input(‘Enter value of A1: ‘);
۷٫ B1=input(‘Enter value of B1: ‘);
۸٫ C1=input(‘Enter value of C1: ‘);
۹٫
۱۰٫ A2=input(‘Enter value of A2: ‘);
۱۱٫ B2=input(‘Enter value of B2: ‘);
۱۲٫ C2=input(‘Enter value of C2: ‘);
۱۳٫
۱۴٫ P=input(‘Enter value of pressure (mmHg): ‘);
۱۵٫ T1=(B1/(A1-log10(P)))-C1;
۱۶٫ T2=(B2/(A2-log10(P)))-C2;
۱۷٫
۱۸٫ i=1;
۱۹٫ for T=T1:.1:T2
۲۰٫ TT(i)=T;
۲۱٫ x=(P-10^(A2-B2/(T+C2)))/(10^(A1-B1/(T+C1))-10^(A2-B2/(T+C2)));
۲۲٫ X(i)=x;
۲۳٫ Y(i)=(x*10^(A1-(B1/(T+C1))))/P;
۲۴٫ i=i+1;
۲۵٫ end
۲۶٫
۲۷٫ subplot(1,2,1)
۲۸٫ plot(X,TT)
۲۹٫ hold on
۳۰٫ plot(Y,TT)
۳۱٫ subplot(1,2,2)
۳۲٫ plot(X,X)
۳۳٫ hold on
۳۴٫ plot(X,Y)

نمودار را برای مخلوط بنزن و تولوئن رسم کردیم که نتیجه ی آن را در زیر می توانید مشاهده کنید.



منبع:انجمن مهندسی شیمی