Response Surface Methodology Routed Optimization of Performance of Hydroxy Gas Enriched Diesel Fuel in Compression Ignition Engines

: In this study, the response surface methodology (RSM) optimization technique was employed for investigating the impact of hydroxy gas (HHO) enriched diesel on performance, acoustics, smoke and exhaust gas emissions of the compression ignition (CI) engine. The engine was operated within the HHO ﬂow rate range of 0–10 L/min and engine loads of 15%, 30%, 45%, 60% and 75%. The results disclosed that HHO concentration and engine load had a substantial inﬂuence on the response variables. Analysis of variance (ANOVA) results of developed quadratic models indicated the appropriate ﬁt for all models. Moreover, the optimization of the user-deﬁned historical design of an experiment identiﬁed an optimum HHO ﬂow rate of 8 L/min and 41% engine load, with composite desirability of 0.733. The responses corresponding to optimal study factors were 25.44%, 0.315 kg/kWh, 117.73 ppm, 140.87 ppm, 99.37 dB, and 1.97% for brake thermal efﬁciency (BTE), brake speciﬁc fuel consumption (BSFC), CO, HC, noise, and smoke, respectively. The absolute percentage errors (APEs) of RSM were predicted and experimental results were below 5%, which vouched for the reliable use of RSM for the prediction and optimization of acoustics and smoke and exhaust emission characteristics along with the performance of a CI engine.


Introduction
The oil reserves are depleting rapidly and are only sufficient to meet the drastically increasing power demand for the next fifty years [1]. Energy demand is soaring at an unprecedented pace and the available sources are too meagre to satisfy the needs [2,3]. In this scenario, the consumption of diesel as a transportation fuel has also increased by about 40% over the last decade [4]. The agriculture sector makes a major contribution to diesel consumption [5,6] as heavy machinery uses diesel as a fuel [7]. Moreover, automotive diesel engines share 26% of total greenhouse gas emissions into the environment, which is an unignorable threat to the stability of the Earth [8][9][10]. This has motivated researchers to investigate alternative fuels, such as hydroxy gas (HHO) for the versatile dual-fuel

Materials and Methods
This study used a 30 kW Perkins (AD 3.152), 91.4 mm bore length, four stroke diesel engine for experimentation. The specifications of the test engine are shown in Table 1. The experiments were performed for different flow rates of HHO, ranging from 0 L/min to 10 L/min while loads were varied from 0-75%, with an equal increment of 15% for each strategic test run. The physicochemical properties of diesel and hydrogen are presented in Table 2. Diesel fuel was directly supplied to the engine through the fuel injectors. However, HHO gas was supplied to the test engine intake manifold at varying flow rates for the diesel HHO mixture. The schematic of the experimental setup, which includes HHO generator, noise measuring meter, smoke meter, emission analyzer and electric heaters, is displayed in Figure 1. Brake thermal efficiency and brake specific fuel consumption were calculated numerically by utilizing calorific value (CV), brake power (BP), and fuel consumption (FC). The brake power was measured from the integrated control panel with heaters, which indicates the value of voltage and current at different load conditions, varied through electrical switches to turn on/off the heaters. In the experimental setup, heaters acted as a resistance load. The test engine was equipped with three phase AC generator having five breakers. Each breaker was equivalent to a load of 15%, which was applied to the engine through the generator. Simply, all the breakers turned on means the engine is operating on a load of 75%. Fuel consumption was determined by measuring time for the consumption of 100 mL of liquid fuel indicated by a gauged cylinder fixed adjacent to diesel containing tank while calorific value was obtained from Pakistan State Oil (PSO). HHO flow rate was ascertained from the rotameter connected at the output of the HHO generator. An emission analyzer (TESTO 350) was employed as a CO and HC emission content recorder, with a measuring sensitivity of 1 ppm CO, and 10 ppm HC. A smoke (opacity) meter (Wager 6500 manufactured by GasTech), which was in full compliance with the requirements of the SAE J1667 test criteria, was used to notice the smoke within a range of 0.0-100.0%. The engine noise level was measured with a sound level meter (UNI-T UT353), which can accurately measure 30-130 dB sound at a frequency response of 31.5 Hz to 8KHz. HHO was generated by electrolysis of water, which is the most commonly used method. Alkaline hydroxides, for example, KOH and NaOH and so forth, were used for speeding up the reaction [36]. The HHO generation system included AC supply, load controller, transformer, rectifier, reactor, and bubbler. The maximum production capacity Processes 2021, 9,1355 4 of 18 of the unit was 10 L/min. Potassium hydroxide (KOH) was used as a catalyst owing to its higher solubility and affinity for water [22,37]. The flow rate of the produced gas was controlled using a potentiometer which was directly dependent on the current passing through the cell. HHO was generated by electrolysis of water, which is the most commonly used method. Alkaline hydroxides, for example, KOH and NaOH and so forth, were used for speeding up the reaction [36]. The HHO generation system included AC supply, load controller, transformer, rectifier, reactor, and bubbler. The maximum production capacity of the unit was 10 L/min. Potassium hydroxide (KOH) was used as a catalyst owing to its higher solubility and affinity for water [22,37]. The flow rate of the produced gas was controlled using a potentiometer which was directly dependent on the current passing through the cell.

Response Surface Methodology
RSM is a statistical technique used for the estimation of relationships between input and response variables. It adopts linear, quadratic, or higher-order polynomial functions to investigate the statistical significance of the study factors and their interactions. Moreover, the regression concept is used for the prediction and optimization of responses. Over the years, the use of RSM in engineering fields has shown welcome results in the prediction of complex systems.
In the current study, the examined parameters were engine load and HHO concentration. Design-Expert version 11 was used for defining the multi-level historical design. The candidate set was created using user-defined discrete levels. Engine load and HHO concentration were assigned to four and six levels, respectively. The response variables measured were BTE, BSFC, CO, HC, noise, and smoke. The best fit model for each response was selected and analysis of variance (ANOVA) was applied for a better understanding of model attributes. In ANOVA, F is a probability distribution in different samplings, Df is degrees of freedom and the p-value is a statistical measure of variations in samples of a particular property. The decision rule for significance was benchmarked as a p-value less than 0.05. The percentage contribution (PC%) of each model term was calculated, which is a ratio of an aggregate of squared deviations to an individual sum of squares (SOS). PC% is a tool that provides a rough idea about the relative importance of study factors and the interactions.

ANOVA Results
The ANOVA results and fit statistics for BTE are presented in Table 3. The F-value of 1980.51 and p-value less than 0.0001 show that the model for BTE is significant. Moreover, the R 2 value of 0.9976 (refer to Table 4) is close to positive unity and there is sufficient

Response Surface Methodology
RSM is a statistical technique used for the estimation of relationships between input and response variables. It adopts linear, quadratic, or higher-order polynomial functions to investigate the statistical significance of the study factors and their interactions. Moreover, the regression concept is used for the prediction and optimization of responses. Over the years, the use of RSM in engineering fields has shown welcome results in the prediction of complex systems.
In the current study, the examined parameters were engine load and HHO concentration. Design-Expert version 11 was used for defining the multi-level historical design. The candidate set was created using user-defined discrete levels. Engine load and HHO concentration were assigned to four and six levels, respectively. The response variables measured were BTE, BSFC, CO, HC, noise, and smoke. The best fit model for each response was selected and analysis of variance (ANOVA) was applied for a better understanding of model attributes. In ANOVA, F is a probability distribution in different samplings, Df is degrees of freedom and the p-value is a statistical measure of variations in samples of a particular property. The decision rule for significance was benchmarked as a p-value less than 0.05. The percentage contribution (PC%) of each model term was calculated, which is a ratio of an aggregate of squared deviations to an individual sum of squares (SOS). PC% is a tool that provides a rough idea about the relative importance of study factors and the interactions.

ANOVA Results
The ANOVA results and fit statistics for BTE are presented in Table 3. The F-value of 1980.51 and p-value less than 0.0001 show that the model for BTE is significant. Moreover, the R 2 value of 0.9976 (refer to Table 4) is close to positive unity and there is sufficient agreement between predicted and adjusted R 2 . The p values from the ANOVA table show that both load and concentration of HHO are significant. However, the load is significantly contributing to aggregated variations with a PC% of 84.4 compared to fuel concentration (3.5%). The best fitted quadratic model from the fit summary was selected owing to the poor fit and aliased nature of linear and cubic models, respectively. The actual regression equation for BTE is given by Equation (1). (1) The contour plot (see Figure 2a) reveals the impact of load and fuel addition on BTE variation. The red color of the contour region advocates high engine BTE at high load and high HHO concentration. The color gradually shifted to red with an increase in HHO amount. The more explicit variation of response (BTE) is noticeable in Figure 2b. The 3D surface plot shows the rising curve of BTE with positive moments along the load and fuel axes. The maximum thermal efficiency is observable at an engine load of 75% and a 10 L/min flow rate of HHO. The improvement in BTE with HHO enrichment is due to the complete combustion of diesel in the presence of hydroxy gas which resulted from higher mean effective pressure near TDC, owing to the faster flame travel in the case of hydrogen. The dark and light circles above and below the surface represent the experimental and predicted values, respectively. Furthermore, the accuracy of the given models could be assessed using certain diagnostics tests and graphs. In general, small deviations between experimental and predicted results are desirable for efficient models. Figure 3 shows a comparison of actual and predicted BTE. The minimal deviations of predicted values from actual data sets are testimony to a good fit of the quadratic regression model.
ANOVA results for the second response variable, BSFC, are shown in Table 5. The model is significant owing to an F value of 169.80, a p-value less than the designated range, and R 2 (0.9725) close to 1, as indicated in Table 6. The ANOVA findings show the significant effect of both load and HHO on fuel consumption. However, compared on a comparative scale, the load variations were found to have a greater impact on an engine than HHO concentration, as evidenced by PCs of 72.9% and 1.4%, respectively. The quadratic regression equation for BSFC on an actual scale is shown by Equation (2).
(2)  ANOVA results for the second response variable, BSFC, are shown in Table 5. The model is significant owing to an F value of 169.80, a p-value less than the designated range, and R 2 (0.9725) close to 1, as indicated in Table 6. The ANOVA findings show the significant effect of both load and HHO on fuel consumption. However, compared on a comparative scale, the load variations were found to have a greater impact on an engine than HHO concentration, as evidenced by PCs of 72.9% and 1.4%, respectively. The quadratic regression equation for BSFC on an actual scale is shown by Equation (2)  ANOVA results for the second response variable, BSFC, are shown in Table 5. The model is significant owing to an F value of 169.80, a p-value less than the designated range, and R 2 (0.9725) close to 1, as indicated in Table 6. The ANOVA findings show the significant effect of both load and HHO on fuel consumption. However, compared on a comparative scale, the load variations were found to have a greater impact on an engine than HHO concentration, as evidenced by PCs of 72.9% and 1.4%, respectively. The quadratic regression equation for BSFC on an actual scale is shown by Equation (2) (2)   Table 6. Coefficient of determination for BSFC.

Coefficient of Determination Value
The effect of HHO addition and load on the fuel consumption trend of an engine is shown in Figure 4a,b. The contour plot in Figure 4a shows that fuel economy improved with successive addition of HHO to diesel at high loads. Moreover, it is also evident that, for the load range of 15-45%, there are more abrupt variations in BSFC compared to high loads, as indicated by a multi-color region. The response surface curve in Figure 4b shows the decreasing increasing trend of BSFC with load and fuel concentration. The sudden lift in the curve at the culmination is due to increased fuel demand at a high load owing to ample friction resistance. The improved fuel economy with HHO enrichment is primarily because of the higher calorific value of HHO and efficient combustion due to the lean diesel-HHO-air mixture [38,39]. The comparison of predicted and actual BSFC, as given in Figure 5, shows a bit of disorder data near the regression line. The disorderliness is due to the manual use of equipment in calculating BSFC. However, the deviations are not so large and therefore the model is acceptable. The effect of HHO addition and load on the fuel consumption trend of an engine is shown in Figures 4a,b. The contour plot in Figure 4a shows that fuel economy improved with successive addition of HHO to diesel at high loads. Moreover, it is also evident that, for the load range of 15%-45%, there are more abrupt variations in BSFC compared to high loads, as indicated by a multi-color region. The response surface curve in Figure 4b shows the decreasing increasing trend of BSFC with load and fuel concentration. The sudden lift in the curve at the culmination is due to increased fuel demand at a high load owing to ample friction resistance. The improved fuel economy with HHO enrichment is primarily because of the higher calorific value of HHO and efficient combustion due to the lean diesel-HHO-air mixture [38,39]. The comparison of predicted and actual BSFC, as given in Figure 5, shows a bit of disorder data near the regression line. The disorderliness is due to the manual use of equipment in calculating BSFC. However, the deviations are not so large and therefore the model is acceptable.   Similar to performance, quadratic models for emissions were also analyzed using ANOVA. The defined model for CO emissions was significant as shown in Table 7. The results revealed that both factors were significant; however, the percentage contribution of load to overall variations was greater compared to fuel concentration. Moreover, the R 2 value was 0.9819 (refer to Table 8) and there was a reasonable agreement between adjusted and predicted R 2 . In an attempt to see the accuracy of the selected model, the actual versus predicted description in Figure 6  Similar to performance, quadratic models for emissions were also analyzed using ANOVA. The defined model for CO emissions was significant as shown in Table 7. The results revealed that both factors were significant; however, the percentage contribution of load to overall variations was greater compared to fuel concentration. Moreover, the R 2 value was 0.9819 (refer to Table 8) and there was a reasonable agreement between adjusted and predicted R 2 . In an attempt to see the accuracy of the selected model, the actual versus predicted description in Figure 6  (3) Figure 6. Comparison of actual and predicted CO.    Similar to performance, quadratic models for emissions were also analyzed using ANOVA. The defined model for CO emissions was significant as shown in Table 7. The results revealed that both factors were significant; however, the percentage contribution of load to overall variations was greater compared to fuel concentration. Moreover, the R 2 value was 0.9819 (refer to Table 8) and there was a reasonable agreement between adjusted and predicted R 2 . In an attempt to see the accuracy of the selected model, the actual versus predicted description in Figure 6  (3) Figure 6. Comparison of actual and predicted CO. CO emission pattern of the engine subjected to various loads. The emissions are shown with the multi-color scheme, where blue stands for the minimum and red for the maximum. The response surface in Figure 7b depicts the CO variations with load and HHO. The main root of carbon monoxide generation is the partial burning of fuel inside the engine. The addition of hydroxy gas not only reduces the carbon content but also facilitates complete combustion which consequently reduces the emissions [40]. Therefore, a curve is seen to be following a decreasing trend in the presence of HHO.

Coefficient of Determination
Value R 2 0.9819 Adjusted R 2 0.9781 Predicted R 2 0.9684 The variations in emissions of carbon monoxide with load and HHO concentrations are shown in Figure 7a,b. The contour plot (Figure 7a) provides a general illustration of the CO emission pattern of the engine subjected to various loads. The emissions are shown with the multi-color scheme, where blue stands for the minimum and red for the maximum. The response surface in Figure 7b depicts the CO variations with load and HHO. The main root of carbon monoxide generation is the partial burning of fuel inside the engine. The addition of hydroxy gas not only reduces the carbon content but also facilitates complete combustion which consequently reduces the emissions [40]. Therefore, a curve is seen to be following a decreasing trend in the presence of HHO. Similarly, Table 9 presents the ANOVA results of HC emission. The model selected and input variables are significant because of p values less than 0.005. The coefficient of determination, the R 2 value, however, is shown in Table 10. Engine load and HHO concentration had percentage contributions of 82.4% and 10.3% respectively. The comparison Similarly, Table 9 presents the ANOVA results of HC emission. The model selected and input variables are significant because of p values less than 0.005. The coefficient of determination, the R 2 value, however, is shown in Table 10. Engine load and HHO concentration had percentage contributions of 82.4% and 10.3% respectively. The comparison of actual and predicted HC emissions in Figure 8 shows that the selected model is accurate. Equation (4)  (4)  The detailed effect of varying factors on hydrocarbon emissions is shown in Figure  9a,b. The addition of HHO reduced HC emissions for all concentrations and the minimum emissions were found to be for 10 L/min, as shown in Figure 9a. Similarly, the response surface shows the emission variations of each fuel combination and is seen following a decreasing trend. The presence of hydroxy gas reduces HC, while carbon present in lubricating oil and primary diesel fuel is oxidized by excessive oxygen and high combustion temperatures inside the cylinder. Moreover, a relatively short quenching distance and a wider flammability range in the case of gaseous fuel have improved the engine performance in this regard [41]. The detailed effect of varying factors on hydrocarbon emissions is shown in Figure 9a,b. The addition of HHO reduced HC emissions for all concentrations and the minimum emissions were found to be for 10 L/min, as shown in Figure 9a. Similarly, the response surface shows the emission variations of each fuel combination and is seen following a decreasing trend. The presence of hydroxy gas reduces HC, while carbon present in lubricating oil and primary diesel fuel is oxidized by excessive oxygen and high combustion temperatures inside the cylinder. Moreover, a relatively short quenching distance and a wider flammability range in the case of gaseous fuel have improved the engine performance in this regard [41]. In addition to the performance and emission of an engine, factors of noise and smoke have also been considered. When the piston oscillates in the cylinder, it creates vibrations which consequently cause high noise levels. Moreover, when sudden ignition of fuel occurs inside the combustion chamber, it generates pressure waves that increase the intensity of the vibrations [42]. The smoke is produced as the result of a rich air-fuel mixture and lubricant burning in the combustion chamber [12]. Tables 11 and 12 present the ANOVA results for noise and smoke. The quadratic models and study factors for both In addition to the performance and emission of an engine, factors of noise and smoke have also been considered. When the piston oscillates in the cylinder, it creates vibrations which consequently cause high noise levels. Moreover, when sudden ignition of fuel occurs inside the combustion chamber, it generates pressure waves that increase the intensity of the vibrations [42]. The smoke is produced as the result of a rich air-fuel mixture and lubricant burning in the combustion chamber [12]. Tables 11 and 12 present the ANOVA results for noise and smoke. The quadratic models and study factors for both responses were significant. Variations in noise would be more due to HHO concentration rather than load, as shown by percentage contributions of 13.19% and 74.30%. Similarly, the smoke model unveils that both load and fuel amount have a significant impact on smoke produced. Moreover, a model R 2 value close to one (based on Tables 13 and 14) and actual versus predicted diagnostic descriptions (Figures 10 and 11) evidenced the accuracy of the selected models. Equations (5) and (6)      The effect of load and HHO on noise could be studied using the contour plots and response surface presented in Figure 12a,b. The red-colored region at the right top corner of Figure 12a indicates that, with the addition of HHO, the noise level increased and is at a maximum for 10 L/min HHO. The same trend could be seen more explicitly in Figure  12b where the response surface shows the gradual increase in noise level. The increased noise level with the addition of hydroxy gas could be apprehended by improved thermal efficiency and excessive combustion at high pressures inside the chamber [42,43]. The opacity is seen following a decreasing trend with a rise in fuel enrichment and load, as shown by Figure 13a,b. The contour plot and 3D response surface show that the least smoke is found for a blend of diesel with 10 L/min of HHO. The improved performance of an engine in terms of smoke emissions could be attributed to reduced HC emissions, high flame propagation, and high flame temperature of hydrogen [44]. The effect of load and HHO on noise could be studied using the contour plots and response surface presented in Figure 12a,b. The red-colored region at the right top corner of Figure 12a indicates that, with the addition of HHO, the noise level increased and is at a maximum for 10 L/min HHO. The same trend could be seen more explicitly in Figure 12b where the response surface shows the gradual increase in noise level. The increased noise level with the addition of hydroxy gas could be apprehended by improved thermal efficiency and excessive combustion at high pressures inside the chamber [42,43]. The opacity is seen following a decreasing trend with a rise in fuel enrichment and load, as shown by Figure 13a,b. The contour plot and 3D response surface show that the least smoke is found for a blend of diesel with 10 L/min of HHO. The improved performance of an engine in terms of smoke emissions could be attributed to reduced HC emissions, high flame propagation, and high flame temperature of hydrogen [44].
The effect of load and HHO on noise could be studied using the contour plots and response surface presented in Figure 12a,b. The red-colored region at the right top corner of Figure 12a indicates that, with the addition of HHO, the noise level increased and is at a maximum for 10 L/min HHO. The same trend could be seen more explicitly in Figure  12b where the response surface shows the gradual increase in noise level. The increased noise level with the addition of hydroxy gas could be apprehended by improved thermal efficiency and excessive combustion at high pressures inside the chamber [42,43]. The opacity is seen following a decreasing trend with a rise in fuel enrichment and load, as shown by Figure 13a,b. The contour plot and 3D response surface show that the least smoke is found for a blend of diesel with 10 L/min of HHO. The improved performance of an engine in terms of smoke emissions could be attributed to reduced HC emissions, high flame propagation, and high flame temperature of hydrogen [44].

RSM Based Optimization
Optimization is the study of maximizing output. An RSM-based optimization is a method to identify optimized conditions by maximizing or minimizing the study factors. In the current work, the emission and performance parameters of the engine are optimized using the numerical optimization feature of the Design Expert. In the optimization setup shown in Table 15, the goal of maximum was assigned to BTE only, while for smoke, noise, BSFC, CO, and HC the minimum criteria were selected. Moreover, the default in the range criterion for study factors was selected.
The engine operating conditions identified by optimization were 41% engine load and blend of diesel with 8 L/min HHO, both rounded to the nearest whole number. The response variables, corresponding to optimized operating conditions, were 25.44% BTE, 0.315 kg/kWh BSFC, 117.7 3 ppm of CO, 140.86 ppm of HC, 99.4 dB of noise, and smoke of 1.97%. The optimum gained values of study factors and response variables are shown by the red and blue dots in Figure 14. The experimentation and RSM models in the previous sections advocated the use of 10 L/min blended diesel for boosted performance and

RSM Based Optimization
Optimization is the study of maximizing output. An RSM-based optimization is a method to identify optimized conditions by maximizing or minimizing the study factors. In the current work, the emission and performance parameters of the engine are optimized using the numerical optimization feature of the Design Expert. In the optimization setup shown in Table 15, the goal of maximum was assigned to BTE only, while for smoke, noise, BSFC, CO, and HC the minimum criteria were selected. Moreover, the default in the range criterion for study factors was selected. The engine operating conditions identified by optimization were 41% engine load and blend of diesel with 8 L/min HHO, both rounded to the nearest whole number. The response variables, corresponding to optimized operating conditions, were 25.44% BTE, 0.315 kg/kWh BSFC, 117.7 3 ppm of CO, 140.86 ppm of HC, 99.4 dB of noise, and smoke of 1.97%. The optimum gained values of study factors and response variables are shown by the red and blue dots in Figure 14. The experimentation and RSM models in the previous sections advocated the use of 10 L/min blended diesel for boosted performance and reduced emissions. However, at the same time, all the blend percentages were unfavorable for noise and therefore an optimum concentration of 8 L/min sounds reasonable.

RSM Based Optimization
Optimization is the study of maximizing output. An RSM-based optimizatio method to identify optimized conditions by maximizing or minimizing the study f In the current work, the emission and performance parameters of the engine are mized using the numerical optimization feature of the Design Expert. In the optimi setup shown in Table 15, the goal of maximum was assigned to BTE only, while for s noise, BSFC, CO, and HC the minimum criteria were selected. Moreover, the def the range criterion for study factors was selected.
The engine operating conditions identified by optimization were 41% engin and blend of diesel with 8 L/min HHO, both rounded to the nearest whole numbe response variables, corresponding to optimized operating conditions, were 25.44% 0.315 kg/kWh BSFC, 117.7 3 ppm of CO, 140.86 ppm of HC, 99.4 dB of noise, and of 1.97%. The optimum gained values of study factors and response variables are s by the red and blue dots in Figure 14. The experimentation and RSM models in the ous sections advocated the use of 10 L/min blended diesel for boosted performan reduced emissions. However, at the same time, all the blend percentages were unf ble for noise and therefore an optimum concentration of 8 L/min sounds reasonabl  The statistical identification of how optimization involved the overall responses could be studied through composite desirability (D). It is a unitless value in the range of 0-1, with 1 for the best and 0 for the worst case. In the current study, composite desirability is 0.733, which is a clear indication that the optimization settings have achieved favorable outcomes for all responses. The contour plot of desirability is shown in Figure 15. Moreover, the impact of individual responses on the overall setting could be assessed through the individual desirability (d) of each response, as shown by the bar graph in Figure 16. It is evident that d is largest for CO (0.877) and lowest for noise (0.364). The numerical values show that minimizing carbon monoxide emissions would have the greatest impact on the overall settings compared, and minimizing noise would impart the least impact to the setting as a whole. ver, the impact of individual responses on the overall setting could be assessed throug the individual desirability (d) of each response, as shown by the bar graph in Figure 16. I is evident that d is largest for CO (0.877) and lowest for noise (0.364). The numerical value show that minimizing carbon monoxide emissions would have the greatest impact on th overall settings compared, and minimizing noise would impart the least impact to th setting as a whole.

Validation of RSM Results
The obtained RSM multi-optimization results were validated using experimentation The engine was operated on the optimal values of load and HHO concentration and th responses were recorded. The absolute percentage error (APE) between the RSM pre dicted and experimentally obtained results was calculated as shown in Table 16.

Validation of RSM Results
The obtained RSM multi-optimization results were validated using experimentation. The engine was operated on the optimal values of load and HHO concentration and the responses were recorded. The absolute percentage error (APE) between the RSM predicted and experimentally obtained results was calculated as shown in Table 16. The APE shows that the developed RSM models and optimization results are accurate. The predicted results showed a reasonable agreement with the experimental results, with APE of all responses being below 5%. However, the maximum APE of 4.76% was evaluated for BSFC, which may be due to inefficient desirability resulting from manual recording during experimentation. Collectively, the predicted results of the developed models were efficient, which promised the simplification of complex performances with the least investment of time, effort and capital.

Conclusions
The purpose of the current investigation was to examine the impact of the blends of diesel with HHO on performance, noise, smoke and tailpipe emissions. Engine load and blend percentage of HHO were the varying factors. The following conclusions could be obtained from the research.

•
ANOVA analysis of all the developed quadratic models indicated suitable fits.

•
The 10 L/min HHO blended diesel proved valuable for improving performance, smoke, and for containing emissions.

•
The noise increased for all the blended fuels and was maximum for 10 L/min HHO.

•
The optimum blend flow rate among 0-10 L/min was 8 L/min for an engine load of 41%.

•
Optimization revealed a composite desirability value of 0.733 with 25.44%, 0.315 kg/kWh, 117.73 ppm, 140.87 ppm, 99.37 dB, and 1.97% for BTE, BSFC, CO, HC, noise, and smoke respectively. • In the optimization model, the most and least significant factors affecting desirability (D) were CO and noise, respectively. • APE predicted that experimental results were below 5%.
The ANOVA and optimization results indicated the potential of hydroxy gas to be used as an alternative fuel in a CI engine. Thus, the use of HHO in blend percentages with diesel will help to save the stability of the Earth from deteriorating due to significantly reduced exhaust emissions compared to pure diesel. Moreover, the use of the RSM technique is beneficial and could save time and capital.