Session: Modeling, Simulation and Optimization
Session: Modeling, Simulation and Optimization
Session: Modeling, Simulation and Optimization
By
The subject of least square analyses also plays central role in the
application of the theory of statistics, which treats problems involving
random.
The subject of random variables and statistics is beyond the scope of our
presentation, we will therefore use least squares in our case
studies.
General equation:
y = a0 + a1x (2.1)
where,
a0 = (∑yi / n) - a1 (∑xi/n) (2.2)
General equation:
(2.5)
On solving the above matrix, we can get the values of a0, a1 & a2.
Case Studies
Three case studies has been
analyzed, they are:
Fitting Gilliland-Sherwood
data.
Case I
Cp = a0 + a1T
150 1. 426
160 1. 447
170 1. 469
180 1. 492
190 1. 516
200 1. 541
210 1.567
220 1. 596
230 1. 627
240 1. 661
250 1. 696
260 1. 732
270 1. 770
280 1. 808
290 1. 848
300 1. 888
On solving these data’s by Least Square method, we got the
values of a0 and a1 and it was found to be 0.96 and 0.00297
respectively.
On substituting the values of a0 and a1 in the equation 2.1, we got
the predicted value, which seems to be near when compared to
the experimental values.
y = 0.96 + 0.00297x
PredictedValues Experimentalvalues
1.405 1.426
1.434 1.447
1.457 1.469
1.489 1.492
1.509 1.516
1.538 1.541
1.558 1.567
1.587 1.596
1.611 1.627
1.658 1.661
1.684 1.696
1.725 1.732
1.768 1.770
1.795 1.808
1.821 1.848
1.867 1.888
Goodness of fit
2
1.8
1.6
1.4
1.2
1 Fig I (c)
EXPERIMENT
AL
0.8 PREDICTED
0.6
0.4
0.2
0
The above graphEXPshows the
EXPEXPEXPEXPGoodness ofEXP
EXPEXPEXPEXP FitEXP
for
EXPthe Heat
EXPEXP EXP Capacity Data
of Methylcyclohexane.
1 2 3 The 4 5experimental
6 7 8 9 10 and11 predicted
12 13 14 15values showed
above shows the minimum percentage of error. The error percentage
would approximately lies between 2-5%.
Case II
Vapour Liquid Equilibrium to
Wilson Equation
Vapour Liquid equilibrium data were taken from Heptane-
Toluene binary system at 1 atm pressure.
Here we fitted activity coefficient data to Wilson Equation
As it required Non-Linear regression, we have used
polynomial method and we have got the predicted values
which are nearer to the experimental.
xi yi
1.000 0.0000
0.790 0.1259
0.596 0.1509
0.480 0.1392
0.390 0.1250
0.293 0.1111
0.220 0.0950
0.150 0.0707
0.065 0.0290
0.000 0.0000
On solving these data's by Polynomial method, we got the values
of a0, a1 and a2 and it was found to be -0.00425, 0.575, -0.5568
respectively.
On substituting these values on the equation 2.4, we got the
predicted values which were very nearer to the experimental
values.
Predicted Experimental
0.0012 0.000
0.1154 0.1259
0.1495 0.1509
0.1435 0.1392
0.1343 0.1250
0.1102 0.1111
0.0925 0.0950
0.0695 0.0707
0.0278 0.0290
0.0000 0.0000
Goodness of fit
0.16
0.14
0.12
0.1
0.08 EXPERIMENTAL
0.06
PREDICTED Fig II (c)
0.04
0.02
0
EXP 1 EXP 2 EXP 3 EXP 4 EXP 5 EXP 6 EXP 7 EXP 8 EXP 9
The above graph shows the Goodness of Fit for the Vapour Liquid
Equilibrium data. The experimental and predicted lines showed above
shows the minimum percentage of error. The error percentage was
approximately lies between 6 – 7%.
Case III
Mass Transfer Data of Gilliland-
Sherwood Equation
The Mass Transfer Data's were taken and analyzed by Gilliland-
Sherwood Equation
As it required Non-Linear regression, we used the polynomial method
and we have got the predicted values which are nearer to the
experimental values.
xi yi
43.7 0.60
24.2 1.80
51.6 1.87
32.3 1.86
26.1 2.16
92.8 2.17
On solving these data’s by polynomial method, we got the values
of a0, a1 and a2 and it was found to be 16.11, -0.7588 and 0.0053
respectively.
On substituting these values on the equation 2.4, we got the
predicted values which were very nearer to the experimental
values.
PREDICTED EXPERIMENTAL
0.48 0.60
1.62 1.80
1.69 1.87
1.70 1.86
1.92 2.16
1.95 2.17
Goodness of fit
Fig III (c)
2.5
1.5
EXPERIMENT
AL
1 PREDICTED
0.5
0
EXP 1 EXP 2 EXP 3 EXP 4 EXP 5 EXP 6
The above graph shows the Goodness of Fit for the Mass Transfer
Data from Gilliland-Sherwood equation. The experimental and
predicted lines showed above shows the minimum percentage of
error. The error percentage was approximately lies between 20 –
25%. The error was little high since the data was Non-Linear.
Results and Discussion
The three case studies analyzed using least squares and
polynomial regression yield good results.
From the graph we have concluded that the second case study
has come out well with minimum error %. As we see the graph
we can see the two lines very close indicating that the
regression was successful with very less error.
Discussion - Case III
In the third case study, the table III (a) shows the data taken for
Performing regression.
The table III (b) is the table containing the experimental data and predicted
data. Using those values we have plotted a graph III (c) which shows the
goodness of fit.
From the graph we have concluded that the second case study has come out
well with minimum error %. As we see the graph we can see the two lines
very close indicating that the regression was successful with very less error.
The error % is high when compared to the first two cases because the data
is non-linear.
Algorithm
General algorithm for all the three cases:
Step 1: The data’s were taken from the case
studies.
Step 2: Regression was applied to the data using
the formulas which are stated in the beginning.
Step 3: The necessary values of Ao and A1 was
determined.
Step 4: These values are substituted in the general
equations (2.1 for case study I and 2.4 for
Case study II and III).
Step 5: Using that we have determined the
predicted values.
Step 6: A graph was drawn between the
experimental and the predicted data.
Step 7: The error percentage was calculated.
Step 8: The graph plotted is the goodness of fit.
Conclusion
The three sets of data's were analyzed and good results have
been obtained in all three case studies.
The first 2 case studies have come perfectly with minimum error
and the goodness of fit graph is plotted and was found to be
good.
For the third case study the error was little high when compared
to the other two, this is because the data is too non-linear.
Goodness of fit graph was plotted for the third case study and
has come out well.