Name: Rutuja Pathak

Roll No.: BT21CME017

EXPERIMENT 6
2nd Order Complex Reaction Scheme

AIM: To develop and simulate mathematical model equation describing batch

THEORY: All the reactions are 1st order irreversible reactions. The system is
assumed to have a constant holdup and batch operation. Heat is not
accounted since we are only concerned about the concentrations of the
components and the variations due to rate coefficients.
PROBLEM SYNTAX:

#include<stdio.h>

void main()
{
float ca=0.15, cb=0.1, cc=0,cd=0,ce=0.0,cf=0,cg=0;
float k1=0.16, k2=0.05,k3=0.15,k4=0.14,k5=0.03,k6=0.01,k7=0.05,k8=0.05;
float cadot,cbdot,ccdot,cddot,cedot,cfdot,cgdot,tprint=0, time=0,delta=1;
int i;
printf("\nTime\tca\tcb\tcc\tcd\tce\tcf\tcg”);
for(i=0; i<=100; i++)
 {
cadot = - k1*ca*cb -k2*ca*cc - k7*ca*cg - k8*ca*cf;
cbdot = - k1*ca*cb - k4*cb*cd - k6*cb*cc;
ccdot = k1*ca*cb - k2*ca*cc - k3*cc*cd- k5*cc*cc - k6*cc*cb;
cddot = k2*ca*cc - k3*cc*cd - k4*cb*cd;
cedot = k3*cc*cd + k8*ca*cf;
cfdot = k4*cb*cd + k5*cc*cc + k7*ca*cg - k8*ca*cf;
cgdot = k6*cb*cc - k7*ca*cg;
printf("\n%f\t%f\t%f\t%f\t%f\t%f\t%f\t%f",time, ca, cb, cc, cd, ce, cf, cg);
if(time<=tprint){
ca = ca + cadot*delta;
cb = cb + cbdot*delta;
cc = cc + ccdot*delta;
cd = cd + cddot*delta;
ce = ce + cedot*delta;
cf = cf + cfdot*delta;
cg = cg + cgdot*delta;
time = time + delta;
}
tprint = tprint + delta;
}

RESULTS:
Time ca cb cc cd ce cf cg
0 0.15 0.1 0 0 0 0 0
1 0.1476 0.0976 0.0024 0 0 0 0
2 0.145277 0.095293 0.004685 0.000018 0 0 0.000002
3 0.143028 0.093073 0.006861 0.000051 0 0.000001 0.000007
4 0.140849 0.090936 0.008934 0.0001 0 0.000003 0.000013
5 0.138737 0.088877 0.010909 0.000161 0 0.000007 0.000021
6 0.136688 0.086893 0.012793 0.000235 0.000001 0.000013 0.000031
7 0.1347 0.084978 0.014589 0.000319 0.000001 0.000021 0.000042
8 0.13277 0.083131 0.016303 0.000413 0.000002 0.000031 0.000054
9 0.130895 0.081346 0.017938 0.000515 0.000003 0.000044 0.000067
10 0.129073 0.079622 0.019499 0.000625 0.000005 0.000059 0.000081
11 0.127302 0.077955 0.020989 0.000742 0.000007 0.000078 0.000096
12 0.12558 0.076343 0.022411 0.000865 0.00001 0.000099 0.000112
13 0.123904 0.074783 0.023769 0.000994 0.000013 0.000124 0.000128
14 0.122272 0.073272 0.025066 0.001127 0.000018 0.000151 0.000145
15 0.120684 0.071809 0.026305 0.001265 0.000023 0.000182 0.000163
16 0.119137 0.070391 0.027488 0.001406 0.000029 0.000215 0.000181
17 0.117629 0.069016 0.028618 0.00155 0.000036 0.000251 0.000199
18 0.116159 0.067682 0.029698 0.001697 0.000044 0.00029 0.000217
19 0.114725 0.066388 0.030729 0.001845 0.000053 0.000333 0.000236
20 0.113327 0.065132 0.031714 0.001996 0.000064 0.000377 0.000255
21 0.111963 0.063912 0.032655 0.002148 0.000075 0.000425 0.000275
22 0.110631 0.062727 0.033554 0.002301 0.000088 0.000476 0.000294
23 0.109331 0.061575 0.034412 0.002455 0.000103 0.000529 0.000313
24 0.108061 0.060456 0.035232 0.002609 0.000118 0.000584 0.000333
25 0.106821 0.059367 0.036015 0.002764 0.000135 0.000642 0.000352
26 0.105608 0.058308 0.036762 0.002918 0.000153 0.000702 0.000372
27 0.104423 0.057278 0.037475 0.003072 0.000173 0.000765 0.000391
28 0.103265 0.056275 0.038155 0.003226 0.000194 0.00083 0.000411
29 0.102131 0.055298 0.038804 0.003379 0.000217 0.000897 0.00043
30 0.101023 0.054347 0.039424 0.003532 0.000241 0.000966 0.000449
31 0.099938 0.05342 0.040014 0.003683 0.000267 0.001037 0.000468
32 0.098877 0.052517 0.040577 0.003833 0.000295 0.001109 0.000487
33 0.097837 0.051636 0.041113 0.003982 0.000323 0.001184 0.000506
34 0.09682 0.050778 0.041624 0.00413 0.000354 0.00126 0.000525
35 0.095823 0.049941 0.04211 0.004276 0.000386 0.001338 0.000544
36 0.094846 0.049124 0.042572 0.004421 0.000419 0.001417 0.000562
37 0.09389 0.048328 0.043013 0.004565 0.000454 0.001498 0.00058
38 0.092952 0.04755 0.043431 0.004706 0.00049 0.00158 0.000598
39 0.092033 0.046791 0.043828 0.004846 0.000528 0.001663 0.000616
40 0.091132 0.04605 0.044206 0.004984 0.000568 0.001748 0.000634
41 0.090248 0.045326 0.044564 0.00512 0.000609 0.001833 0.000651
42 0.089381 0.044618 0.044903 0.005255 0.000651 0.00192 0.000669
43 0.088531 0.043927 0.045225 0.005387 0.000695 0.002008 0.000686
44 0.087696 0.043252 0.045529 0.005518 0.000741 0.002097 0.000703
45 0.086878 0.042592 0.045816 0.005646 0.000788 0.002186 0.000719
46 0.086074 0.041947 0.046088 0.005773 0.000836 0.002276 0.000736
47 0.085285 0.041316 0.046345 0.005897 0.000886 0.002367 0.000752
48 0.08451 0.040699 0.046586 0.00602 0.000937 0.002459 0.000768
49 0.083749 0.040095 0.046814 0.00614 0.000989 0.002551 0.000783
50 0.083002 0.039505 0.047027 0.006259 0.001043 0.002644 0.000799
51 0.082268 0.038927 0.047228 0.006375 0.001098 0.002737 0.000814
52 0.081547 0.038362 0.047415 0.006489 0.001155 0.002831 0.000829
53 0.080838 0.037808 0.047591 0.006602 0.001212 0.002925 0.000844
54 0.080141 0.037266 0.047754 0.006712 0.001271 0.00302 0.000859
55 0.079457 0.036735 0.047906 0.00682 0.001331 0.003114 0.000873
56 0.078783 0.036216 0.048048 0.006927 0.001393 0.003209 0.000887
57 0.078122 0.035707 0.048178 0.007031 0.001455 0.003305 0.000901
58 0.077471 0.035208 0.048299 0.007133 0.001519 0.0034 0.000915
59 0.07683 0.034719 0.048409 0.007233 0.001584 0.003496 0.000928
60 0.076201 0.034241 0.048511 0.007332 0.00165 0.003591 0.000941
61 0.075581 0.033771 0.048603 0.007428 0.001717 0.003687 0.000954
62 0.074971 0.033312 0.048686 0.007522 0.001785 0.003782 0.000967
63 0.074372 0.032861 0.048761 0.007615 0.001854 0.003878 0.00098
64 0.073781 0.032419 0.048827 0.007705 0.001924 0.003974 0.000992
65 0.0732 0.031985 0.048886 0.007794 0.001995 0.004069 0.001004
66 0.072628 0.03156 0.048937 0.007881 0.002067 0.004165 0.001016
67 0.072065 0.031143 0.048981 0.007966 0.00214 0.00426 0.001028
68 0.07151 0.030734 0.049018 0.008049 0.002214 0.004355 0.00104
69 0.070964 0.030332 0.049048 0.008131 0.002289 0.00445 0.001051
70 0.070426 0.029939 0.049072 0.00821 0.002365 0.004544 0.001062
71 0.069896 0.029552 0.049089 0.008288 0.002441 0.004639 0.001073
72 0.069374 0.029173 0.0491 0.008365 0.002518 0.004733 0.001084
73 0.06886 0.028801 0.049105 0.008439 0.002596 0.004827 0.001094
74 0.068353 0.028435 0.049105 0.008512 0.002675 0.00492 0.001105
75 0.067853 0.028076 0.049099 0.008583 0.002755 0.005013 0.001115
76 0.067361 0.027724 0.049088 0.008653 0.002835 0.005106 0.001125
77 0.066876 0.027378 0.049072 0.008721 0.002916 0.005199 0.001135
78 0.066398 0.027038 0.049051 0.008787 0.002997 0.005291 0.001144
79 0.065927 0.026704 0.049025 0.008852 0.00308 0.005382 0.001154
80 0.065462 0.026377 0.048995 0.008916 0.003162 0.005474 0.001163
81 0.065003 0.026054 0.04896 0.008978 0.003246 0.005565 0.001172
82 0.064551 0.025738 0.048922 0.009038 0.00333 0.005655 0.001181
83 0.064106 0.025427 0.048879 0.009097 0.003414 0.005745 0.00119
84 0.063666 0.025121 0.048832 0.009155 0.0035 0.005834 0.001199
85 0.063232 0.024821 0.048782 0.009211 0.003585 0.005923 0.001207
86 0.062804 0.024526 0.048728 0.009266 0.003671 0.006012 0.001215
87 0.062382 0.024236 0.04867 0.009319 0.003758 0.0061 0.001224
88 0.061966 0.02395 0.048609 0.009371 0.003845 0.006187 0.001231
 89 0.061554 0.02367 0.048545 0.009422 0.003932 0.006274 0.001239
90 0.061149 0.023394 0.048478 0.009472 0.00402 0.006361 0.001247
91 0.060748 0.023123 0.048408 0.00952 0.004109 0.006447 0.001255
92 0.060353 0.022856 0.048335 0.009567 0.004197 0.006532 0.001262
93 0.059963 0.022594 0.04826 0.009613 0.004286 0.006617 0.001269
94 0.059578 0.022335 0.048181 0.009658 0.004376 0.006701 0.001276
95 0.059198 0.022082 0.048101 0.009701 0.004466 0.006785 0.001283
96 0.058822 0.021832 0.048017 0.009744 0.004556 0.006868 0.00129
97 0.058452 0.021586 0.047932 0.009785 0.004646 0.00695 0.001297
98 0.058086 0.021344 0.047844 0.009825 0.004737 0.007032 0.001303
99 0.057724 0.021106 0.047754 0.009864 0.004828 0.007114 0.00131
100 0.057367 0.020872 0.047662 0.009902 0.004919 0.007194 0.001316

GRAPH:

ca, cb, cc vs time

0.08
0.06
0.04
0.02
0
0
4
8
12
16
20
24
28
32
36
40
44
48
52
56
60
64
68
72
76
80
84
88
92
96
100

Time

ca cb cc
 concentraition concentration

0.002
0.004

0.006

0.008

0.012
0.001
0.002
0.003
0.004
0.005
0.006
0.007
0.008

0.01
0

0
0 0
4 4
8 8
12 12
16 16
20 20
24 24
28 28

cd, ce vs time
cf, cg vs time
32 32
36 36
40 40

cd
cf
44 44
48 48

time

time
cg 52 52

ce
56 56
60 60
64 64
68 68
72 72
76 76
80 80
84 84
88 88
92 92
96 96
100 100