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

[chemblog]

[h=2]برنامه متلب برای رسم نمودار T-xy مخلوط های دو جزیی پیروی کننده از قانون رائولت[/h]

 

[P O U R I A]

در این قسمت می خواهیم با استفاده از نرم افزار متلب کدی را بنویسیم که نمودار های T-xy و x-y را برای مخلوط های دو جزیی که از قانون رائولت پیروی می کنند را رسم نماید.
برای این کار ابتدا نرم افزار باید مقدار فشار را از کاربر دریافت کرده و در آن فشار نقطه جوش دو ماده ی خالص را محاسبه نماید و سپس در بازه ی بین این دو دما مقادیر x و y را محاسبه نماید.برای محاسبه ی نقطه جوش می توان از معادله ی آنتوان استفاده نمود.بنابراین کدی که خواهیم نوشت باید ضرائب A , B , C معادله ی آنتوان را نیز برای هر کدام از دو ماده از کاربر سوال کند.
Log[SUB]10[/SUB](P[SUP]sat[/SUP])=A-B/(T+C)
برای محاسبه ی x:
P[SUB]t[/SUB]=x[SUB]1[/SUB]P[SUB]1[/SUB][SUP]sat[/SUP]+(1-x[SUB]1[/SUB])P[SUB]2[/SUB][SUP]sat[/SUP]
x=(P[SUB]t[/SUB]-P[SUB]2[/SUB][SUP]sat[/SUP])/(P[SUB]1[/SUB][SUP]sat[/SUP]-P[SUB]2[/SUB][SUP]sat[/SUP])
دقت کنید که جز ۱ جز فرارتر در سیستم است.
برای محاسبه ی y:
yP[SUB]t[/SUB]=xP[SUB]1[/SUB][SUP]sat[/SUP]
y=xP[SUB]1[/SUB][SUP]sat[/SUP]/P[SUB]t[/SUB]
کد نوشته شده برای برنامه را در زیر مشاهده می کنید.

۱٫   % This program plot the T-xy @ x-y diagrams for mixtures that obey from Raoult’s law.
2.   %———————————————————————————–
۳٫   clc
4.   clear all
5.
6.   A1=input(‘Enter value of A1: ‘);
7.   B1=input(‘Enter value of B1: ‘);
8.   C1=input(‘Enter value of C1: ‘);
9.
10.   A2=input(‘Enter value of A2: ‘);
11.   B2=input(‘Enter value of B2: ‘);
12.   C2=input(‘Enter value of C2: ‘);
13.
14.   P=input(‘Enter value of pressure (mmHg): ‘);
15.   T1=(B1/(A1-log10(P)))-C1;
16.   T2=(B2/(A2-log10(P)))-C2;
17.
18.   i=1;
19.    for T=T1:.1:T2
20.        TT(i)=T;
21.        x=(P-10^(A2-B2/(T+C2)))/(10^(A1-B1/(T+C1))-10^(A2-B2/(T+C2)));
22.        X(i)=x;
23.        Y(i)=(x*10^(A1-(B1/(T+C1))))/P;
24.        i=i+1;
25.    end
26.
27.    subplot(1,2,1)
28.    plot(X,TT)
29.    hold on
30.    plot(Y,TT)
31.    subplot(1,2,2)
32.    plot(X,X)
33.    hold on
34.    plot(X,Y)

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

​