Biodiesel Production
Current work focuses on maximizing biodiesel yield using Taguchi’s approach and response surface methodology (RSM). In Table 2, the rank of all the parameters is computed along with the S/N ratios at different levels. Here, the delta value is used for determining the ranks of performance parameters, the highest delta value denotes rank 1, and so on. Hence, the molar ratio with rank 1 is identified as the most effective performance parameter in biodiesel production.
Table 2. Response table of operating parameters
Parameters
|
Level
|
Delta
|
Rank
|
1
|
2
|
3
|
Molar Ratio (oil/methanol)
|
38.17
|
38.43
|
38.79
|
0.60
|
1
|
Reaction Time (minutes)
|
38.38
|
38.55
|
38.48
|
0.23
|
3
|
Catalyst Concentration (wt.%)
|
38.59
|
38.48
|
38.31
|
0.13
|
2
|
Table 3 depicts the biodiesel yield (in percentage) obtained for each run (or experiment) in each day session. Data shows less than 1% deviation in biodiesel yield; thus, the mean value of three sessions is considered for further calculations. The mean biodiesel yield (in percentage) at various experiments (or run) and the S/N ratios are also given in Table 3. The mean biodiesel yield ranges from 77.45% (minimum) to 89.14% (maximum), and similar results are obtained by Kumar et al. (2017) [42] and Kumar et al. (2018) [43]. The highest mean yield of 89.14% is obtained in run 26 by MR 9:1, RT 120 mins. and CC 0.75 wt.%. Also, the optimum condition, which is measured by the highest value of the S/N ratio (38.94), is achieved by two experiments, i.e., run 22 and run 26, respectively.
Table 3. Biodiesel yield and S/N ratio of L27 experiments performed
Run
|
Molar Ratio
|
Reaction Time
(min.)
|
Catalyst Conc.
(wt. %)
|
Biodiesel Yield (%)
|
Mean Yield (%)
|
S/N ratio
|
Morning
|
Afternoon
|
Post-Afternoon
|
1
|
6:1
|
90
|
0.5
|
80.23
|
80.32
|
80.29
|
80.28
|
38.18
|
2
|
6:1
|
90
|
0.75
|
81.22
|
81.29
|
81.21
|
81.24
|
38.15
|
3
|
6:1
|
90
|
1
|
77.43
|
77.44
|
77.48
|
77.45
|
37.74
|
4
|
6:1
|
105
|
0.5
|
84.59
|
84.67
|
84.57
|
84.61
|
38.56
|
5
|
6:1
|
105
|
0.75
|
82.11
|
82.23
|
82.17
|
82.17
|
38.24
|
6
|
6:1
|
105
|
1
|
78.73
|
78.97
|
78.91
|
78.87
|
37.97
|
7
|
6:1
|
120
|
0.5
|
84.76
|
85.18
|
85.06
|
85.00
|
38.49
|
8
|
6:1
|
120
|
0.75
|
81.67
|
81.71
|
81.72
|
81.70
|
38.34
|
9
|
6:1
|
120
|
1
|
77.81
|
78.00
|
77.95
|
77.92
|
37.84
|
10
|
7.5:1
|
90
|
0.5
|
83.97
|
84.11
|
84.13
|
84.07
|
38.42
|
11
|
7.5:1
|
90
|
0.75
|
81.59
|
81.76
|
81.78
|
81.71
|
38.30
|
12
|
7.5:1
|
90
|
1
|
83.08
|
83.11
|
83.17
|
83.12
|
38.42
|
13
|
7.5:1
|
105
|
0.5
|
86.34
|
86.43
|
86.49
|
86.42
|
38.75
|
14
|
7.5:1
|
105
|
0.75
|
82.57
|
82.67
|
82.71
|
82.65
|
38.33
|
15
|
7.5:1
|
105
|
1
|
84.89
|
84.91
|
84.96
|
84.92
|
38.58
|
16
|
7.5:1
|
120
|
0.5
|
83.58
|
83.57
|
83.62
|
83.59
|
38.50
|
17
|
7.5:1
|
120
|
0.75
|
82.12
|
82.16
|
82.23
|
82.17
|
38.26
|
18
|
7.5:1
|
120
|
1
|
82.19
|
82.27
|
82.32
|
82.26
|
38.28
|
19
|
9:1
|
90
|
0.5
|
85.99
|
86.12
|
86.19
|
86.10
|
38.69
|
20
|
9:1
|
90
|
0.75
|
88.31
|
88.34
|
88.34
|
88.33
|
38.92
|
21
|
9:1
|
90
|
1
|
85.29
|
85.33
|
85.40
|
85.34
|
38.64
|
22
|
9:1
|
105
|
0.5
|
88.73
|
88.79
|
88.79
|
88.77
|
38.94
|
23
|
9:1
|
105
|
0.75
|
87.14
|
87.12
|
87.16
|
87.14
|
38.87
|
24
|
9:1
|
105
|
1
|
86.83
|
86.78
|
86.76
|
86.79
|
38.73
|
25
|
9:1
|
120
|
0.5
|
86.94
|
87.00
|
86.94
|
86.96
|
38.83
|
26
|
9:1
|
120
|
0.75
|
89.11
|
89.16
|
89.15
|
89.14
|
38.94
|
27
|
9:1
|
120
|
1
|
84.51
|
84.55
|
84.50
|
84.52
|
38.56
|
Initially, the effects of operating parameters on the response are analyzed using Taguchi’s approach, and Fig. 5 shows the impact of different operating parameters on the biodiesel yield at all three levels. The peak value defines the optimum value of all parameters for different S/N ratio plots. Consequently, MR 9:1 (level 3), RT 105 min. (level 2), and CC 0.5 wt.% (level 1) achieves the highest biodiesel yield. It indicates that the above process parameters are required to obtain the maximum biodiesel yield at the optimized conditions.
Regression analysis is performed by the data acquired in biodiesel yield from eq. 2, and the percentage of the coefficient of regression (97.72%) is given by eq. 3: -
Yield = 5.9 + 1.72 MR + 1.379 RT + 1.79 CC + 0.584 (MR)2 - 0.00588 (RT)2 - 0.306 (CC)2 - 0.0267 MR*RT + 0.873 MR*CC - 0.0351 RT*CC (2)
R2 = 97.72% (3)
Experimental and predicted values of biodiesel yield and the error percentage in different experiments of Taguchi’s L27 approach are given in Table 4. The result shows a very close relationship between the experimental and predicted values of biodiesel yield. The experimental yield values (%age) show a minimal deviation from predicted yield values (%age), and the difference is measured as the mean error percentage for each experimentation. The highest mean error percentage is about 1.17%, which indicates that the experimental yield values are not more than 1.2% deviated from the predicted yield values.
Table 4. Experimental and predicted values for Taguchi's L27 approach
Experiment (Run)
|
Molar Ratio
|
Reaction Time
(min.)
|
Catalyst Conc.
(wt. %)
|
Yield (%)
|
Error
(%)
|
Experimental
|
Predicted
|
1
|
6:1
|
90
|
0.5
|
80.28
|
81.07
|
0.9884
|
2
|
6:1
|
90
|
0.75
|
81.24
|
80.81
|
-0.5186
|
3
|
6:1
|
90
|
1
|
77.45
|
77.07
|
-0.4807
|
4
|
6:1
|
105
|
0.5
|
84.61
|
84.79
|
0.2045
|
5
|
6:1
|
105
|
0.75
|
82.17
|
81.64
|
-0.6481
|
6
|
6:1
|
105
|
1
|
78.87
|
79.23
|
0.4558
|
7
|
6:1
|
120
|
0.5
|
85.00
|
84.04
|
-1.1371
|
8
|
6:1
|
120
|
0.75
|
81.70
|
82.65
|
1.1675
|
9
|
6:1
|
120
|
1
|
77.92
|
77.93
|
0.0163
|
10
|
7.5:1
|
90
|
0.5
|
84.07
|
83.37
|
-0.8268
|
11
|
7.5:1
|
90
|
0.75
|
81.71
|
82.19
|
0.5908
|
12
|
7.5:1
|
90
|
1
|
83.12
|
83.33
|
0.2556
|
13
|
7.5:1
|
105
|
0.5
|
86.42
|
86.55
|
0.1503
|
14
|
7.5:1
|
105
|
0.75
|
82.65
|
82.48
|
-0.1973
|
15
|
7.5:1
|
105
|
1
|
84.92
|
84.95
|
0.0390
|
16
|
7.5:1
|
120
|
0.5
|
83.59
|
84.16
|
0.6761
|
17
|
7.5:1
|
120
|
0.75
|
82.17
|
81.85
|
-0.3890
|
18
|
7.5:1
|
120
|
1
|
82.26
|
82.01
|
-0.2984
|
19
|
9:1
|
90
|
0.5
|
86.10
|
86.00
|
-0.1143
|
20
|
9:1
|
90
|
0.75
|
88.33
|
88.27
|
-0.0696
|
21
|
9:1
|
90
|
1
|
85.34
|
85.50
|
0.1873
|
22
|
9:1
|
105
|
0.5
|
88.77
|
88.47
|
-0.3413
|
23
|
9:1
|
105
|
0.75
|
87.14
|
87.84
|
0.7983
|
24
|
9:1
|
105
|
1
|
86.79
|
86.40
|
-0.4524
|
25
|
9:1
|
120
|
0.5
|
86.96
|
87.36
|
0.4616
|
26
|
9:1
|
120
|
0.75
|
89.14
|
88.50
|
-0.7114
|
27
|
9:1
|
120
|
1
|
84.52
|
84.75
|
0.2754
|
Analysis of variance (ANOVA)
ANOVA determines the significance of various parameters in the process. The ANOVA analysis (as represented in Table 5) defines the importance of process parameters for biodiesel yield. The significance of a model is determined by its p-value (<0.0500), where a higher p-value states the process parameter is insignificant. In this case, MR and CC have a p-value of 0.000, making them highly significant model terms, while RT has a p-value of 0.011; thus, it is also a significant model term but less effective than MR and CC.
Table 5. ANOVA table for biodiesel yield
Source
|
Degree of Freedom
|
Seq. SS
|
Adj. SS
|
Adj. MS
|
F-value
|
P-value
|
Contribution
(%)
|
MR
|
2
|
163.115
|
163.115
|
81.5573
|
110.46
|
0.000
|
75.67
|
RT
|
2
|
12.262
|
12.262
|
6.1308
|
8.30
|
0.011
|
5.69
|
CC
|
2
|
34.278
|
34.278
|
17.1391
|
23.21
|
0.000
|
15.90
|
MR*RT
|
4
|
4.834
|
4.834
|
1.2084
|
1.64
|
0.256
|
---
|
MR*CC
|
4
|
29.056
|
29.056
|
7.2641
|
9.84
|
0.004
|
---
|
RT*CC
|
4
|
9.320
|
9.320
|
2.3300
|
3.16
|
0.078
|
---
|
Residual Error
|
8
|
5.907
|
5.907
|
0.7383
|
|
|
2.74
|
Total
|
26
|
258.771
|
|
|
|
|
100
|
Percentage contribution of MR, RT, and CC is also computed using the ANOVA table (Table 5), where Seq. SS terms are considered for calculating the percentage contribution of a parameter in the process. The highest contributing parameter is the molar ratio (75.67%), followed by catalyst concentration (15.9%) and reaction time (5.69%), while a residual error of 2.74% is also observed. It indicates that a slight variation in the molar ratio will significantly affect biodiesel yield.
Response surface methodology (RSM)
RSM is based upon the different mathematical and statistical practices established through the experimental design of adequate empirical models. It shows the relationship between the observed and theoretically obtained values from different empirical models. Fig. 6 depicts a connection between predicted values and actual values for biodiesel yield in percentage; it also shows the closeness of the two values.
Response surface methodology (RSM) also gives the interactive effects between different factors in surface-contour plots, which helps better understand the effects of two factors on response [44]. The surface-contour plots for the above experiments are shown in Fig. 7, 8, and 9. In Fig. 7, the effects of MR and CC are depicted, while the RT is taken as 105 min. (actual factor). The value of RT is taken from Taguchi’s approach, which gives optimized biodiesel yield at 105 min. The figure shows that biodiesel production enhances when MR (alcohol to oil) is higher that is similar to the results obtained by Fan et al. (2011) [44]. Consequently, the catalyst concentration between 0.6 and 0.9 wt.% shows a slight decline, going against the findings of Yadav et al. (2018) [45]. However, the plot shows the potential for higher biodiesel yield with the rise in molar ratio beyond 9:1.
The interactive effects between CC and RT are represented in Fig. 8. Here, the curve's flatness shows a lower percentage contribution of these parameters, as described in the ANOVA table (table 5). CC and RT's percentage contributions are less than MR; thus, putting molar ratio to 9:1 (actual factor) will show less fluctuation in the surface counter curve. The transesterification process in solar-assisted biodiesel production takes more time than conventional methods; as a result, the optimized RT is observed between 95 and 110 mins. However, after 110 mins, a decline in biodiesel yield is observed due to reverse reaction. Similar observations are also noticed by Mihankhah et al. (2016) [12].
In Fig. 9, the interaction of RT with MR is shown for biodiesel yield (as a response). As discussed, the biodiesel yield rises with the molar ratio; thus, a similar trend is observed in the figure. Higher molar ratio with RT more than 95 min. gives more than 83% of biodiesel yield. From different published works, it was observed that the transesterification process requires a 3:1 M ratio (alcohol to triglyceride) for alkyl ester (biodiesel) conversion [46, 47]; however, to get a biodiesel yield greater than 80%, an MR of 9:1 or more is favorable. In contrast, the higher MR results in extra alcohol content, creating difficulty in glycerol separation. Thus, a relevant MR between 9 and 12 is acceptable for biodiesel production, while a higher MR of 15:1 reduces the biodiesel yield due to difficulty in glycerol recovery [33, 45]. RSM analysis concluded that the optimized parameters for desirable biodiesel yield (91.1%) are molar ratio 8.92:1, reaction time 108.97 minutes, and catalyst concentration 0.61 wt.%, respectively. Five consecutive experiments are performed to confirm the values obtained through the RSM model, and the average of the five experiments is found to be 91.9% biodiesel yield.
Fuel properties
For considering any fuel for internal combustion (IC) engines, the fuel properties must match the engine input parameters as per ASTM standards. Thus, the desirable fuel properties are measured for biodiesel (LSME) and raw linseed oil, and compared with diesel. Table 6 compares desirable fuel properties of LSME and diesel as per ASTM testing methods.
Table 6. Desirable fuel properties for CI engine
Properties
|
ASTM Testing method
|
SI Unit
|
Diesel
|
Raw
linseed oil
|
Linseed methyl ester (LSME)
|
This work
|
[48]
|
[49]
|
Density at 15℃
|
D-1298
|
kg/m3
|
809.6
|
914.0
|
869.8
|
860
|
872
|
Kinematic viscosity at 40℃
|
D-445
|
×10-6 m2/s
|
2.815
|
27.317
|
5.453
|
5
|
8.2
|
Calorific value
|
D-240
|
×103 kJ/kg
|
45.5
|
35.2
|
36.7
|
35.6
|
37.5
|
Flash point
|
D-93
|
℃
|
72
|
>210
|
162
|
187
|
161
|
Pour point
|
D-97
|
℃
|
-9
|
-15
|
-11
|
---
|
---
|
Cetane index
|
D-4737
|
---
|
48.1
|
40.2
|
42.7
|
---
|
55
|
The density of the conventional diesel is 809.6 kg/m3, whereas raw linseed oil has a higher density of 914.0 kg/m3. After the transesterification of linseed oil, the linseed methyl ester (LSME) is produced [50], and the density is reduced to 869.8 kg/m3, which is under the ASTM standards [1]. In contrast, the kinematic viscosity of LSME is found to be 5.453 × 10-6 m2/s which is twice of diesel (2.815 × 10-6 m2/s) and is close to the upper limit of the ASTM standards [1, 13].
The calorific value (CV) of a fuel is the heating capacity of the fuel or the energy released when the fuel is burnt. The CV of diesel is measured maximum (45.5 × 103 kJ/kg), and LSME is measured about 36.7 × 103 kJ/kg. Flash point and pour point of LSME are 162℃ and -11℃, while the cetane index is 42.7, respectively. Hence, it is concluded that the desirable properties of biodiesel are comparable to diesel (under ASTM standards) [1, 13]. As a result, a biodiesel blend with diesel could be used in the existing CI engines with no modifications, or pure biodiesel could directly be used with few modifications in the CI engine.
Fatty acid methyl ester (FAME) composition
The FAME composition of the oil is an essential parameter. It gives the composition of all the FAME content present in the biodiesel. Ultima make series 2100 gas chromatography (specifications in Table 7) is used to measure the FAME composition of biodiesel (LSME).
Table 7. Specifications of gas chromatography
Make
|
Ultima
|
Model
|
Series 2100
|
Detector
|
Flame ionization detector (FID)
|
Carrier gas
|
Nitrogen
|
Other gases
|
Hydrogen and zero air
|
Oven temperature
|
240℃
|
Rate of temperature rise
|
10℃/min
|
The gas chromatogram for LSME is shown in Fig. 10. The curve shows the different FAME compositions in percentage, while Table 8 presents the FAME composition of linseed methyl ester (LSME).
Table 8. FAME composition in wt.%
Fatty acid
|
Chemical formula
|
Degree of Saturation
|
Linseed methyl ester
|
This work*
|
[51]*
|
[52]*
|
Palmitic
|
C16H32O2
|
16:0
|
5.41
|
6.58
|
5.69
|
Stearic
|
C18H36O2
|
18:0
|
4.94
|
4.43
|
5.58
|
Oleic
|
C18H34O2
|
18:1
|
18.73
|
18.51
|
20.59
|
Linoleic
|
C18H32O2
|
18:2
|
15.95
|
17.25
|
15.80
|
Linolenic
|
C18H30O2
|
18:3
|
54.12
|
53.21
|
51.38
|
Total saturated
|
-
|
-
|
10.45
|
11.01
|
11.90
|
Total unsaturated
|
-
|
-
|
89.35
|
88.99
|
88.10
|
*All values are in wt.%
The linolenic acid content is highest (54.12%) in LSME, whereas other higher fatty acids are linoleic acid (15.95%) and oleic acid (18.73%). However, a small quantity of palmitic acid (5.41%) and stearic acid (4.94%) is also present in the linseed methyl ester (as depicted in Table 8). The experimentally obtained FAME values are comparable to earlier work, and it is found that the values obtained (Table 8) do not differ significantly from other published data [36, 51–53].
Comparative analysis of results
A detailed comparison of various biodiesel production methods with feed oil is represented in Table 9. It is concluded that solar-assisted biodiesel production is well competitive with other techniques in biodiesel yield. However, the zero-power requirement in the solar-assisted process makes it superior to other methods.
Table 9. Comparative analysis of present work and other published works
Reference
|
Feed oil
|
Methods used
|
Power
|
Biodiesel yield
|
Predicted
|
Linseed oil
|
---
|
---
|
91.1%
|
Experimental (RSM model)
|
Linseed oil
|
Solar-assisted
|
0
|
91.9%
|
Experimental (ANOVA)
|
Linseed oil
|
Solar-assisted
|
0
|
91.6%
|
Mohite et al. [13]
|
Karanja-Linseed oil mixture
|
Mechanical starring and conventional heating
|
180 W
|
78.9%
|
Singh et al. [15]
|
Microalgae oil
|
Ultrasonic
|
50 W
|
98%
|
Uddin et al. [16]
|
Waste cooking oil
|
Mechanical starring and conventional heating
|
210 W
|
79%
|
Hsiao et al. [20]
|
Soybean oil
|
Microwave-assisted
|
300 W
|
96.6%
|
Kumar et al. [21]
|
Pongamia pinnata seed oil
|
Microwave-assisted
|
300 W
|
96%
|
Mohan et al. [40]
|
Semal oil
|
Ultrasonic cavitation and mechanical stirring
|
780 W
|
90.7%
|
Kumar et al. [43]
|
Pongamia oil
|
Mechanical starring and conventional heating
|
180 W
|
90.2%
|
*Bold is used to represent present work results