Submission of Term Work': Subject - Cpms

‘SUBMISSION OF TERM WORK’
SUBJECT - CPMS
Submitted by
PRAKHAR VERMA
PRN – 1614110020
SEAT NO. – 2014390209
Department of Chemical Engineering

BHARATI VIDYAPEETH DEEMED (TO BE) UNIVERSITY
COLLEGE OF ENGINEERING
PUNE-411043
2019-2020
BHARATI VIDYAPEETH ( DEEMED TO BE) UNIVERSITY
COLLEGE OF ENGINEERING
PUNE-4110
CERTIFICATE
This is to certify that, Prakhar Verma of class B-Tech Chemical Sem -8 , roll no. 09, has
completed all the Practical Work, satisfactorily in the subject of Chemical Process
Modelling and Simulation, in the department of Chemical Engineering as prescribed by
Bharati Vidyapeeth (Deemed to be University), College of Engineering in the academic year
2019 - 2020.
Date- 15/9/20
DR. PV CHAVAN DR. S. N. SHINDE EXTERNAL

(HOD) (subject teacher)
INDEX
S.N. Experiment Name Pg.No

1 1-2
SIMULATION OF II ORDER
BATCH REACTION
2 3-11
GRAVITY FLOW TANK
3 12-17
BUBBLE POINT
CALCULATION
4 18-22
SIMULATION OF THREE
CSTRs IN SERIES
5 23-24
NON- ISOTHERMAL CSTR
EXPERIMENT 1
SIMULATION OF II ORDER BATCH REACTION (MATLAB)
t= C=
0 10.0000
0.5000 9.7590
1.0000 9.5350
1.5000 9.3261
2.0000 9.1311
2.5000 8.9485
3.0000 8.7774
3.5000 8.6167
4.0000 8.4655
4.5000 8.3230
5.0000 8.1887
5.5000 8.0617
6.0000 7.9417
6.5000 7.8280
7.0000 7.7201
7.5000 7.6178
8.0000 7.5206
8.5000 7.4282
9.0000 7.3401
9.5000 7.2563
10.0000 7.1763
10.5000 7.1000
11.0000 7.0272
11.5000 6.9575
12.0000 6.8909
12.5000 6.8272
13.0000 6.7661
13.5000 6.7076
14.0000 6.6515
14.5000 6.5977
15.0000 6.5461
15.5000 6.4965
16.0000 6.4488
16.5000 6.4030
17.0000 6.3590
17.5000 6.3166
18.0000 6.2758
18.5000 6.2365
19.0000 6.1986
19.5000 6.1622
20.0000 6.1270
ans =0
>>
.M FILE CODE
functiondC=second(t,C)%defining ODE function

%C = concentration, kmol/m^3
dC=zeros(1,1); % A column vector for ODE function
F=@(t,C)[(-.01*C^2)+ (0.05*C)];
[t,C]=ode45(F,[0 20],[10])
plot(t,C(:,1))
xlabel('time,s')
ylabel('concentration,kmol/m^3')
title('second order batch reaction')
GRAPH
EXPERIMENT -2
GRAVITY FLOW TANK (MATLAB)
>>.m file code
function dy=gravity(t,y)
dy=zeros(2,1);
dy(1)=(0.0107*y(2))-(0.00205*(y(1))^2);
dy(2)=0.311-(0.0624*y(1));
>>.m file code
[T,Y]=ode45(@gravity,[0 1200],[3.4 2.05])

plot(T,Y(:,1),'-o',T,Y(:,2),'-*')
xlabel('Process time, s')
ylabel('height ,ft')
title('Gravity Flow tank')
>> gravityj
T=
1.0e+03 *
0.0010
0.0021
0.0031
0.0042
0.0094
0.0146
0.0198
0.0250
0.0336
0.0421
0.0507
0.0593
0.0696
0.0800
0.0904
0.1007
0.1107
0.1208
0.1308
0.1409
0.1511
0.1614
0.1717
0.1819
0.1946
0.2073
0.2200
0.2327
0.2458
0.2588
0.2719
0.2850
0.2975
0.3101
0.3227
0.3352
0.3522
0.3691
0.3861
0.4030
0.4194
0.4357
0.4521
0.4684
0.4877
0.5070
0.5262
0.5455
0.5671
0.5887
0.6103
0.6319
0.6579
0.6840
0.7101
0.7362
0.7649
0.7937
0.8224
0.8511
0.8811
0.9111
0.9411
0.9711
1.0011
1.0311
1.0611
1.0911
1.1184
1.1456
1.1728
1.2000
Y=
3.4000 2.0500
3.3987 2.1530
3.3987 2.2561
3.3997 2.3592
3.4019 2.4621
3.4290 2.9728
3.4813 3.4705
3.5561 3.9474
3.6502 4.3967
3.8389 5.0616
4.0560 5.6173
4.2862 6.0518
4.5159 6.3626
4.7760 6.5787
5.0020 6.6377
5.1834 6.5657
5.3159 6.3916
5.3985 6.1544
5.4399 5.8790
5.4458 5.5908
5.4229 5.3084
5.3766 5.0404
5.3152 4.8072
5.2458 4.6173
5.1738 4.4733
5.0873 4.3577
5.0113 4.3065
4.9514 4.3096
4.9082 4.3541
4.8806 4.4294
4.8697 4.5200
4.8728 4.6119
4.8862 4.6968
4.9054 4.7665
4.9276 4.8199
4.9497 4.8551
4.9700 4.8733
4.9926 4.8751
5.0071 4.8576
5.0129 4.8298
5.0121 4.7992
5.0077 4.7717
5.0010 4.7502
4.9938 4.7371
4.9874 4.7312
4.9811 4.7313
4.9775 4.7371
4.9768 4.7454
4.9779 4.7535
4.9800 4.7611
4.9824 4.7655
4.9841 4.7658
4.9850 4.7640
4.9855 4.7613
4.9852 4.7586
4.9843 4.7577
4.9837 4.7581
4.9837 4.7589
4.9838 4.7595
4.9840 4.7595
4.9841 4.7591
4.9840 4.7588
4.9839 4.7587
4.9839 4.7590
4.9840 4.7592
4.9840 4.7592
4.9840 4.7590
4.9840 4.7590
4.9839 4.7590
4.9840 4.7591
4.9840 4.7591
4.9840 4.7591
4.9840 4.7591
>>
GRAPH
EXPERIMENT - 3
BUBBLEPOINT CALCULATION (MATLAB)
>>.m file code
%Ques Find the Bubble point Temperature of VLE given that pressure is 760 mm of Hg(1atm) and liquid
phase composition of
%binary mixture is 0.5 and antoine constant A1=-3848.08,B1=17.53,A2=-4328.11, B2=17.591,
P=760,Xa=0.5;Xb=0.5,T=373,h=0.1.
A1=-3848.08;
B1=17.53;
A2=-4328.11;
B2=17.591;
%Pressure= P mm of Hg
P=760;
Xa=0.5;
Xb=0.5;
%Temperature= T K
T=373;
h=0.5;
%Saturated Pressure= Psat1 mm of Hg
%Saturated Pressure= Psat2 mm of Hg
Psat1=2.718^((A1/T)+B1);
Psat2=2.718^((A2/T)+B2);
%Vapour Pressure= Pvap mm of hg
Pvap=Xa*Psat1+Xb*Psat2;
while(abs(Pvap-P))>=1
if(Pvap>P)
T=T-h
elseif(Pvap<P)
T=T+h
end
Psat1=2.718^((A1/T)+B1);
Psat2=2.718^((A2/T)+B2);
Pvap=Xa*Psat1+Xb*Psat2;
end
Psat1
Psat2
Pvap
T
>>Bubblepoint
T=
372.5000
T=
372
T=
371.5000
T=
371
T=
370.5000
T=
370
T=
369.5000
T=
369
T=
368.5000
T=
368
Psat1 =
1.1791e+03
Psat2 =
340.0957
Pvap =
759.6101
T=
368
>>
EXPERIMENT - 4
SIMULATION OF THREE CSTRs IN SERIES (MATLAB)
>>.m file code
function dy = CSTR_f(t,y);
dy = zeros(3,1);
t=2;
k=0.5;
CA0=1.8;
dy(1)=((CA0-y(1))/t)-(k*y(1));
dy(2)=((y(1)-y(2))/t)-(k*y(2));
dy(3)=((y(2)-y(3))/t)-(k*y(3));
>>.m file code
[T,Y]=ode45(@CSTR_f,[0 2.9],[0.4 0.2 0.1])
plot(T,Y(:,1),'*',T,Y(:,2),'-o',T,Y(:,3),'.-')
xlabel('Time, s')
ylabel('Concentration, kmol A/m^3')
title('3 CSTR in Series')
>> seriescstr
T=
0
0.0402
0.0804
0.1206
0.1608
0.2333
0.3058
0.3783
0.4508
0.5233
0.5958
0.6683
0.7408
0.8133
0.8858
0.9583
1.0308
1.1033
1.1758
1.2483
1.3208
1.3933
1.4658
1.5383
1.6108
1.6833
1.7558
1.8283
1.9008
1.9733
2.0458
2.1183
2.1908
2.2633
2.3358
2.4083
2.4808
2.5533
2.6258
2.6983
2.7708
2.8031
2.8354
2.8677
2.9000
Y=
0.4000 0.2000 0.1000

0.4197 0.2002 0.1000
0.4386 0.2008 0.1000
0.4568 0.2017 0.1000
0.4743 0.2029 0.1001
0.5040 0.2058 0.1002
0.5317 0.2096 0.1005
0.5575 0.2140 0.1009
0.5814 0.2189 0.1014
0.6037 0.2243 0.1020
0.6244 0.2301 0.1028
0.6437 0.2362 0.1038
0.6616 0.2425 0.1049
0.6783 0.2490 0.1062
0.6938 0.2556 0.1076
0.7082 0.2622 0.1091
0.7216 0.2689 0.1108
0.7341 0.2755 0.1125
0.7457 0.2821 0.1144
0.7565 0.2887 0.1164
0.7665 0.2951 0.1185
0.7759 0.3015 0.1206
0.7845 0.3077 0.1228
0.7926 0.3137 0.1251
0.8001 0.3196 0.1274
0.8071 0.3254 0.1298
0.8136 0.3310 0.1322
0.8197 0.3364 0.1346
0.8253 0.3416 0.1371
0.8305 0.3467 0.1395
0.8354 0.3516 0.1420
0.8399 0.3563 0.1444
0.8441 0.3608 0.1469
0.8480 0.3651 0.1493
0.8516 0.3693 0.1517
0.8550 0.3733 0.1541
0.8582 0.3772 0.1564
0.8611 0.3809 0.1587
0.8638 0.3844 0.1610
0.8663 0.3878 0.1632
0.8687 0.3910 0.1654
0.8697 0.3924 0.1664
0.8707 0.3937 0.1674
0.8716 0.3950 0.1683
0.8725 0.3964 0.1693
>>
EXPERIMENT - 5
Non Isothermal CSTR
>>.m file code
function dy = non_f(t,y);
dy = zeros(4,1);
f0=40;
v=48;
ca0=0.5;
ca=0.245;
t=600;
tj=594.6;
fj=49.9;
t0=530;
vj=3.85;
e=30000;
r=1.99;
u=150;
ah=250;
tj0=530;
cp=0.75;
cj=1;
kc=4;
tset=600;
f=40-(10*(48-y(1)));
fj=49.9-(kc*(600-y(3)));
k=7.08*(exp(-(e)/(r*y(3))));
dy(1)=f0-f;
dy(2)=((f0*ca0)-(f*y(2))-(v*k*y(2)))/y(1);
dy(3)=((f0*t0)-(f*y(3)))-((-30000*v*k*y(2))/(50*cp))-(((u*ah)*(y(3)-y(4)))/(50*cp))/y(1);
dy(4)=((fj*(tj0-y(4)))/vj)+(((u*ah)*(y(3)-y(4)))/(62.3*cj*vj));
>>.m file code
[T,Y]=ode45(@non_f,
[0,0.2,0.41,0.61,0.81,1.01,1.21,1.61,1.81,2.21,2.41,2.61,2.81,3.01,3.21,3.41,3.61,3.81,4.010],[48 0.245 600
594.6]);
plot(T,Y(:,1),'-.',T,Y(:,2),'-+',T,Y(:,3),'-o',T,Y(:,4),'-*')
xlabel('Time, h')
ylabel('V,Ca,T,Tj')
title('Nonisothermal CSTR')
legend('V,Ca,T,Tj')
CONCLUSION
From the experiments written in the practical, we can conclude that we have studied modelling and
simulation for a
1. II ORDER BATCH REACTION

2. GRAVITY FLOW TANK
3. THREE CSTRs IN SERIES and
4. A NON- ISOTHERMAL CSTR
We also studied the method to find the bubble point temperature of VLE for a given condition.

Submission of Term Work': Subject - Cpms

Uploaded by

Copyright:

Available Formats

Submission of Term Work': Subject - Cpms

Uploaded by

Document Information

Original Title

Copyright

Available Formats

Share this document

Share or Embed Document

Sharing Options

Did you find this document useful?

Is this content inappropriate?

Copyright:

Available Formats

Submission of Term Work': Subject - Cpms

Uploaded by

Copyright:

Available Formats

‘SUBMISSION OF TERM WORK’

SEAT NO. – 2014390209

Department of Chemical Engineering

DR. PV CHAVAN DR. S. N. SHINDE EXTERNAL

S.N. Experiment Name Pg.No

functiondC=second(t,C)%defining ODE function

>>.m file code

>>.m file code

[T,Y]=ode45(@gravity,[0 1200],[3.4 2.05])

>>.m file code

SIMULATION OF THREE CSTRs IN SERIES (MATLAB)

>>.m file code

[T,Y]=ode45(@CSTR_f,[0 2.9],[0.4 0.2 0.1])

0.4000 0.2000 0.1000

Non Isothermal CSTR

>>.m file code

1. II ORDER BATCH REACTION

You might also like