‘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.