Multi-Objective Optimization of Energy Consumption and Surface Quality in Nanofluid SQCL Assisted Face Milling

Considering the significance of improving the energy efficiency, surface quality and material removal quantity of machining processes, the present study is conducted in the form of an experimental investigation and a multi-objective optimization. The experiments were conducted by face milling AISI 1045 steel on a Computer Numerical Controlled (CNC) milling machine using a carbide cutting tool. The Cu-nano-fluid, dispersed in distilled water, was impinged in small quantity cooling lubrication (SQCL) spray applied to the cutting zone. The data of surface roughness and active cutting energy were measured while the material removal rate was calculated. A multi-objective optimization was performed by the integration of the Taguchi method, Grey Relational Analysis (GRA), and the Non-Dominated Sorting Genetic Algorithm (NSGA-II). The optimum results calculated were a cutting speed of 1200 rev/min, a feed rate of 320 mm/min, a depth of cut of 0.5 mm, and a width of cut of 15 mm. It was also endowed with a 20.7% reduction in energy consumption. Furthermore, the use of SQCL promoted sustainable manufacturing. The novelty of the work is in reducing energy consumption under nano fluid assisted machining while paying adequate attention to material removal quantity and the product’s surface quality.


Introduction
Recently, immoderate energy consumption has become one of the most severe problems in the manufacturing industry, and conserving energy has become a necessity for the industry.In manufacturing processes, machine tools are the primary energy consuming devices [1].Gutowski et al. [2] have concluded that 14.8% of total energy was utilized during the material removal process.According to Newman, 6-40% of energy can be saved by choosing optimum cutting parameters, cutting tools and cutting paths [3].Additionally, the energy consumption of machine tools is also associated with an increased environmental cost.Therefore, there is a dire need for the optimization of cutting parameters based on minimum energy criteria.Material removal rate, cutting force, tool life, surface roughness (SR), power and energy consumption are the responses measured during the milling process.However, the optimization of a single objective has many limitations.The conclusion of the above discussion is that this problem necessitates the need for a framework to optimize conflicting objectives simultaneously.
The manufacturing industry is facing a new challenge to reduce energy consumption due to the higher demand forenergy and increasing concerns regarding fossil fuel depletion [4].However, more difficulties arise when reducing energy consumption without compromising on the surface quality of the product.In the machining process, the typical interactions between the cutting tool-work surface and the cutting tool-chip cause severe friction and intensive material deformation.Even, the harder cutting tool becomes worn due to an excessive mechanical load on the tool edges.Moreover, the thermal adversities caused by friction at the adjacent surfaces instigate unfavourability in machining.This unfavourability demeans the machining productivity and consumes higher amounts of energy.A significant portion of consumed power is lost due to friction during the cutting process [5].
Thermal upsetting at machining interfaces can be pacified by using cooling agents.Likewise, the effects of friction can be limited by employing lubricating agents.Owing to this fact, the use of cooling-lubricating agents in machining is widespread to date.However, the growing concerns for the environment and human health-related problems have forced manufacturers to look for an alternative as well as better cooling-lubrication modes [6][7][8].
Considering those aspects, for decades, the alternative mode of cooling-lubrication in machining has been the use of small droplets of agents under high air pressure so that the impinged sprays penetrate in between the interfaces.These droplets serve two functions: (i) they reduce friction by creating a thin film in between the material bodies [9], and (ii) they enhance heat dissipation by absorbing and removing heat from that zone [10].This approach is reported in the literature under a number of different terms such as: small quantity cooling lubrication (SQCL), minimum quantity cooling lubrication (MQCL), minimum quantity lubrication (MQL), and near dry machining (NDM) etc., although the different terms exhibit slight differences in the way the cooling is delivered.Many research works show the positive sides of SQCL [11,12].However, more advantages are gained when SQCL was applied using nano-fluids in suspension.The improved tribological properties of the nano-fluids themselves result in several other benefits in machining, such as a reduced coefficient of friction, lower cutting forces, improved surface finish, reduced cutting temperature, and lower tool wear associated with prolonged tool life.
The use of nanoparticles such as Cu, CuO, Al 2 O 3 , MoS 2 , SiO 2, and diamond/graphite nanoparticles in the base fluids in the machining process under MQL conditions has become a new trend.However, in depth analysis of the application of water-based nano-fluids in the machining process is still a common research area [13].Mintsa et al. [14] measured the thermal conductivity of alumina and copper oxide mixed in water (water-based nano-fluids).They showed that the effective thermal conductivity increased when increasing the particle volume fraction and when reducing the size of nanoparticles.In another study, Zhang et al. [15] also measured the effective thermal conductivity of distilled water based nanofluids.The results revealed that effective thermal conductivity increases with the increase of nanoparticle concentration.Lvet al. [16] utilized graphene oxide/silicon dioxide hybrid nanoparticle in water-based nano-fluids in evaluating the machining characteristics of MQL technology.Results revealed a significant reduction in the worn scar diameter (WSD) and the coefficient of friction (COF) using hybrid GO/SiO 2 water-based MQL technology in the milling process.
Modern manufacturing aims at reducing energy consumption as well as achieving high quality within small timeframes [17].The International Organization for Standardization is currently working on an ISO/DIS 14955 and ISO/DIS 50001 series which would assess the energy efficiency of machine tools and their environmental effects.Kara and Li [18] developed an empirical model to characterize the relationship between energy consumption and machining parameters and reached the apparent conclusion that specific energy consumption decreases when increasing the material removal rate.Abhang and Hameedullah [19] investigated the productivity and power consumption of machine Energies 2019, 12, 710 3 of 22 tools in the turning of EN-31 steel.Results showed that when the tool nose radius, cutting speed, feed rate and depth of cut were at lower levels, less power was consumed by the machine tool.Bhattacharya et al. [20] performed experiments to study surface quality and power consumption during the high-speed machining of AISI 1045 steel.The cutting speed was found to be the most significant cutting parameter for both surface roughness and power consumption.However, all experiments were performed under dry cutting conditions.Camposeco-Negrete [21] optimized the cutting parameters for minimum energy consumption in the turning of AISI 6061.It was found that feed rate has the most significant effect on the surface roughness and energy consumption; lower values of feed rate resulted in better surface finish, but with higher energy consumption.Pusavec et al. [22] conducted many experiments in the dry, MQL, cryogenic and hybrid cooling assisted machining of nickel-based Inconel alloy 718.Cutting forces, power, tool life, material removal rate and cooling and lubrication performance were studied.The cooling and lubrication techniques were found useful for reducing cutting forces and improving tool life.
Velchevet al. [23] developed an approach for the optimization of cutting parameters for minimum energy in the turning process.A mathematical model was developed that considered the cutting parameters to express the direct energy consumption of machine tool.Kumar et al. [24] performed a multi-objective optimization by assigning different weights to energy consumption, surface quality, and productivity.The authors have concluded that the depth of cut and the use of radius were vital factors for analytical hierarchy process (AHP) and entropy weight methods respectively.Maruda et al. [25] studied the influence of cooling conditions on-chip formation and tool wear.The authors have concluded that the application of MQCL reduced cutting tool wear.
To conserve resources while maintaining improved machinability, there is a need for multiple objectives optimization using different algorithms.Grey relation algorithm (GRA) is one of the most efficient algorithms used for multi-objective optimization [26].The GRA is based on a simple calculation, and it can provide an accurate analysis of small sample size.Moreover, it can assign variable weights to each response.It reduces complex multi-objective objective problems into a single object effectively and efficiently.Mia et al. [27] successfully employed a GRA method to optimize the cutting force and surface roughness in end milling under MQL conditions.On the other side, Deb et al. [28] proposed a fast and elitist heuristic method know as Non-dominated Sorting Genetic Algorithm NSGA-II; researchers have been successfully using the NSGA II to optimize machining parameters.Rao et al. [29] employed NSGA II for a multi-objective optimization for multiple machining processes.Then, Li et al. [30] reported the use of NSGA-II for the machining of Ti-6Al-4V alloy machining.
It is evident that the chip-tool interface endures friction, which requires lubrication and cooling.At the same time, energy reduction is a prime motto of the sustainable production system.For that reason, the whole machining system needs to be optimized considering multiple objectives.On the one hand, GRA cannot provide us with a Pareto front solution to find trade-off relationships.However, on the other hand, NSGA-II cannot provide information on a complex problem as a single characteristic.So, to overcome the drawbacks of each algorithm and to present a solution to a complex problem in detail, the multiple objectives are simultaneously optimized in this study using both a conventional method (GRA) as well as using a heuristic one (NSGA-II).Besides these aspects, other novelties include the use of Cu nano-fluid in distilled water and its use in the form of SQCL, reduction of energy consumption, and improvement of material removal rate.

Experimental Procedure
The face milling experiments were performed on a CNC machine tool (Carver 400M_RT) with a spindle power of 5.6 KW and a maximum rotation speed of 6000 rpm.AISI-1045 steel was used as work material.A three fluted carbide cutting tool with a diameter of 24 mm was used for machining.Figure 1 shows the schematic, experimental setup, and equipment used in measurements.For nano-fluid based SQCL, a solution of copper nanoparticles suspended in distilled water was employed.Diionized water (DI water) was used as a coolant due to the fact that it has no harmful effects on workers' health, and it is of low cost compared to other lubricants.Moreover, the addition of copper nanoparticles in distilled water increases its thermal stability, as reported by Mintsa et al. [14].Copper nanoparticles with a 60 nm diameter were mixed with distilled water to prepare the solution (encouraged by reference [31]).The dispersion stability of nanoparticles is essential to investigate their performance.Thereby, the sodium dodecyl sulfate (SDS) was used as the surfactant For nano-fluid based SQCL, a solution of copper nanoparticles suspended in distilled water was employed.Diionized water (DI water) was used as a coolant due to the fact that it has no harmful effects on workers' health, and it is of low cost compared to other lubricants.Moreover, the addition of copper nanoparticles in distilled water increases its thermal stability, as reported by Mintsa et al. [14].
Copper nanoparticles with a 60 nm diameter were mixed with distilled water to prepare the solution (encouraged by reference [31]).The dispersion stability of nanoparticles is essential to investigate their Energies 2019, 12, 710 5 of 22 performance.Thereby, the sodium dodecyl sulfate (SDS) was used as the surfactant in the mixture to achieve good stability.The stability mechanism of nano-fluid was measured in terms of zeta potential Z. Zeta potential values of more than 30 mV showed the high stability of the nano-fluid and its value was found to be positive 45 mV, depicting that the nano-fluid was highly stable.Figure 2 shows the preparation procedure of the Cu based nanofluid.in the mixture to achieve good stability.The stability mechanism of nano-fluid was measured in terms of zeta potential Z. Zeta potential values of more than 30 mV showed the high stability of the nanofluid and its value was found to be positive 45 mV, depicting that the nano-fluid was highly stable.
Figure 2 shows the preparation procedure of the Cu based nanofluid.The focus of the study is to investigate the effect of cutting parameters on the responses above.Therefore, the technological parameters related to the SQCL system were kept constant by doing preliminary experiments and following the existing literature [11,12].The nano-fluid SQCL was sprayed using a single nozzle at an air pressure of 4 bars, a lubricant flow rate of 300 ml/h; and a nozzle angle of 45°.The maximum material removal volume (MRV) was set to 2400 mm 3 .Experiments were designed according to the Taguchi L27 orthogonal array.Cutting parameters and their levels were selected by making preliminary trial runs and also from the published literature [32].The focus of the study is to investigate the effect of cutting parameters on the responses above.Therefore, the technological parameters related to the SQCL system were kept constant by doing preliminary experiments and following the existing literature [11,12].The nano-fluid SQCL was sprayed using a single nozzle at an air pressure of 4 bars, a lubricant flow rate of 300 mL/h; and a nozzle angle of 45 • .The maximum material removal volume (MRV) was set to 2400 mm 3 .Experiments were designed according to the Taguchi L27 orthogonal array.Cutting parameters and their levels were selected by making preliminary trial runs and also from the published literature [32].
Figure 3 shows the cause-effect diagram for the efficiency of the machining process.In the current study, machining parameters (spindle speed, feed rate, depth of cut, and width of cut) were selected.Table 1 contains the levels of cutting parameters.Figure 3 shows the cause-effect diagram for the efficiency of the machining process.In the current study, machining parameters (spindle speed, feed rate, depth of cut, and width of cut) were selected.Table 1 contains the levels of cutting parameters.

Response Measurements
The responses were measured for three successive runs, and finally, the average of the measured data was adopted as the mean value for further computation.

Production Rate
The production rate was determined by the material removal rate (mm 3 /min) which was computed using Equation (1) [32].
where , , are the number of tool flutes, depth of cut and width of cut, respectively.

Surface Quality
The quality of the produced machined surface was defined by the average of the roughness profile that the machining process was generating.For that purpose, the average surface roughness parameter (SR) was selected as the representative index of surface quality.It was measured using a

Response Measurements
The responses were measured for three successive runs, and finally, the average of the measured data was adopted as the mean value for further computation.

Production Rate
The production rate was determined by the material removal rate (mm 3 /min) which was computed using Equation (1) [32].
where z, a p , a e are the number of tool flutes, depth of cut and width of cut, respectively.

Surface Quality
The quality of the produced machined surface was defined by the average of the roughness profile that the machining process was generating.For that purpose, the average surface roughness parameter (SR) was selected as the representative index of surface quality.It was measured using a portable stylus profilometer Mitutoyo Surftest SJ-410 at three equally separated points, and the average values were taken to mitigate errors [33].

Active Cutting Energy
The Portable Power Cell model PPC-3 (Figure 4) was used to measure power demand during the cutting process.To measure the active energy consumption, a power sensor was attached to the main bus of the electrical cabinet of the CNC machine at a sampling frequency of 10Hz.Equation ( 2) can be used to estimate active energy consumption [32].
Energies 2019, 12, 710 where ACE; active cutting energy (kJ), P i ; active power measured in watts (W), t i ; time, and k; data samples during the process.In the experimental analysis, the lubricant delivery system used a separate power resource that was assumed to be constant.The measurement system for surface roughness and power consumption is shown in Figure 4.
The Portable Power Cell model PPC-3(Figure 4) was used to measure power demand during the cutting process.To measure the active energy consumption, a power sensor was attached to the main bus of the electrical cabinet of the CNC machine at a sampling frequency of 10Hz.Equation ( 2) can be used to estimate active energy consumption [32].
where ACE; active cutting energy (kJ), ; active power measured in watts (W), ; time, and ; data samples during the process.In the experimental analysis, the lubricant delivery system used a separate power resource that was assumed to be constant.The measurement system for surface roughness and power consumption is shown in Figure4.

Results and Discussion
In this section, the multi-objective optimization is performed using Grey relational analysis (GRA) and NSGA-II methods to deal with difficult decision-making situations.Furthermore, the influence of the milling parameter on the surface quality, material removal rate and the cutting energy consumption is discussed.Finally, the mechanism of nano fluid assisted machining favorability is presented.The framework of the multi-objective optimization problem is shown in Figure 5.The framework explains the importance of energy consumption as a response.A comparison has been performed between the conventional method when only two responses (surface quality and productivity) are considered and a proposed method where energy is a sustainable parameter has been added.In the conventional method, equal weights are given to both responses to give them equal importance.However, in the proposed method an entropy method has been employed to the weights of each response according to the sensitivity of the response.In the proposed method, optimization has been performed using two algorithms: GRA and NSGA-II.The reason for employing both optimization algorisms is to get benefits from both methods to understand the complex multi-objective problem fully.The application of GRA is straight forward and quick.It converts the multi-objective optimization problem into a single objective and provides Grey

Results and Discussion
In this section, the multi-objective optimization is performed using Grey relational analysis (GRA) and NSGA-II methods to deal with difficult decision-making situations.Furthermore, the influence of the milling parameter on the surface quality, material removal rate and the cutting energy consumption is discussed.Finally, the mechanism of nano fluid assisted machining favorability is presented.The framework of the multi-objective optimization problem is shown in Figure 5.The framework explains the importance of energy consumption as a response.A comparison has been performed between the conventional method when only two responses (surface quality and productivity) are considered and a proposed method where energy is a sustainable parameter has been added.In the conventional method, equal weights are given to both responses to give them equal importance.However, in the proposed method an entropy method has been employed to the weights of each response according to the sensitivity of the response.In the proposed method, optimization has been performed using two algorithms: GRA and NSGA-II.The reason for employing both optimization algorisms is to get benefits from both methods to understand the complex multi-objective problem fully.The application of GRA is straight forward and quick.It converts the multi-objective optimization problem into a single objective and provides Grey relational grades (GRG).The GRG is a performance index of the GRA and provides information as a single characteristic in optimization problems.For example, a GRG value of 0.8678 highlights the best experiment in the L27 design.However, the heuristics are mostly recommended for the optimization instead of classic approaches, due to their capacity for a more precise global search.It is nearly impossible to get the optimal values of all responses simultaneously.So, NSGA-IIhas been used to obtain Pareto fronts to find the tradeoff relationships between a pair of two objects.GRA provides us the local optimal results.However, NSGA-II provides global optimal results.single characteristic in optimization problems.For example, a GRG value of 0.8678 highlights the best experiment in the L27 design.However, the heuristics are mostly recommended for the optimization instead of classic approaches, due to their capacity for a more precise global search.It is nearly impossible to get the optimal values of all responses simultaneously.So, NSGA-IIhas been used to obtain Pareto fronts to find the tradeoff relationships between a pair of two objects.GRA provides us the local optimal results.However, NSGA-II provides global optimal results.Experimental design L27 along with process parameters and obtained results are presented in Table 2.  Experimental design L27 along with process parameters and obtained results are presented in Table 2.

Grey Relational Analysis
Grey relational analysis was used to solve the complicated multi-objective optimization problem.There are three simple steps for this algorithm.In the first step, to avoid variability and different units, each response value was normalized.If the required objective is 'the larger the better', then the results were normalized using Equation (3).If the objective is 'the smaller the better', then the results were normalized using Equation (4) [26], The grey relation coefficients (GRC) for each factor were calculated by Equation (5).Finally, Grey relational grades (GRG) were computed using Equation (6) [26].
∆0 i (k) is known as the deviation sequence and it is the absolute value.Minimum and maximum values of ∆0 i (k) are denoted by ∆ min and ∆ max , respectively γ i represents the GRG of each experiment provided in Appendix A. The values of ξ lies between 0 and 1 and it is called a distinguishing coefficient, normally it is taken as 0.5.The w k represents the normalized weightage of factor k.
Traditionally, equal weights are assigned to objectives in simple problems.However, in complex problems, it is not easy to assume equal weights for each criterion.In complex multi-attribute decision making (MADM) problems, the weight of each criterion is calculatedby prioritization, expert opinion, and evaluation [34].Finding an appropriate and effective way to assign the weights is a dire need of the MADM.Subjective and objective are the two well-known weight assignment methods.The examples of subjective weight assignment methods are the AHP, Delphi method, and weight least square method; in these methods, weights are assigned according to the preference of the decision makers.
Moreover, subjective weight methods are less accurate and have some drawbacks [35].However, in the objective based weight assignment method, weights are assigned on the inherent information of indexes to determine weights of indexes.Entropy-based weight assignment is an objective weight assignment method which not only avoids human-made disturbances and assumptions but also provides results according to the sensitivity of the objectives.In that way, more realistic and accurate Energies 2019, 12, 710 10 of 22 weights are calculated [35].For these reasons, in this paper, the weight of each response was calculated in a systematic manner using entropy weighting.
In common practice, the concept of entropy is originated from thermo-science.This method finds the weight factor based on the influence of various input parameters on responses [24].According to which, a higher weight should be assigned to the objective with higher variation in its values.In the first step, the mean values of each coefficient for each response were calculated based on Table 3.For example, to calculate the mean values of the Grey relation coefficients of the width of cut at cutting level 1 (5 mm), all experiments were selected with a width of cut at 5 mm and then the mean was taken.Similarly, the mean was taken at level 2 (10 mm) and level 3 (15 mm), respectively.The range was calculated from the difference between the maximum and minimum values of the mean of GRC.The Ranges (R) for each factor were computed using Equation (7).

Development of RSM Models and Analysis of Variance
The second order Response Surface Methodology (RSM) models for the three objective functions: material removal rate, surface roughness, and active cutting energy were formulated using the measured data presented in Table 2. Process parameters in coded form affected the response more evenly, thereby accelerating the convergence speed of the algorithm and improving the precision of the algorithm.Coded variables were sometimes plotted on the x-axis and can easily be calculated by Equations ( 10)- (13).
A = (N − 1700)/500 Coded Variables = 1, 0, −1 It is necessary to understand the relationship between cutting parameters and Grey relational grades forthe optimization of cutting parameters.Response surface methodology was employed to generate the second-order RSM model for MRR, SR, and ACE.Mia et al. [36] worked on the investigation of machinability characteristics and concluded that well-fitting mathematical models could be obtained using RSM.Equations ( 14)-( 16) describe a second-order empirical model for MRR, SR, ACE respectively.Moreover, Equation (17) shows the mathematical relation of Grey relational grade regarding input parameters.

Non-Dominated Sorting Genetic Algorithm NSGA-II
The genetic algorithm is a heuristic algorithm to optimize NP-hard problems inspired by a natural selection process in bio-science systems.Srinivas and Deb [37] proposed NSGA, and later Table 5 presents the analysis of variance (ANOVA) results for the regression model.ANOVA results illustrated that the coefficient of determination R-square value was 99.02%, and R-(adj) was 98.03%, which indicated the model's excellent compatibility with experimental results.Experiment No. 9 (in bold in Table 4) was identified as the optimal trade-off having combination (N = 1200 rev min , f = 320 mm min , a p = 0.5 mm, a e = 15 mm) and can be verified from Figure 6 as having the highest GRG value.

Non-Dominated Sorting Genetic Algorithm NSGA-II
The genetic algorithm is a heuristic algorithm to optimize NP-hard problems inspired by a natural selection process in bio-science systems.Srinivas and Deb [37] proposed NSGA, and later Deb et al. [38] proposed NSGA-II, which elevates all the short-comings of NSGA, NSGA-II can obtain true Pareto optimal solutions based on the rank of non-domination.In NSGA-II, minimization of fitness is assumed, and each solution is given a rank which is equal to its non-domination level, e.g., one is considered the best level.Therefore, it obtains solutions not just within the experimental runs, but global space is used for searching for the optimal solution.A flowchart of GRA coupled with NSGA-II has been shown in Figure 7. assignment of weights does not consider the priority of responses.Therefore, the optimal solution is not guaranteed.In random weight assignment methods, weights are assigned to each objective function and then combined them to make a single objective.However, although these methods are simple to use, weight selection is still a challenging problem.In this paper, the weight of each response was calculated in a systematic manner using entropy weighting of GRA as discussed in Section 4.1.To convert multi-objective optimization into a single objective optimization problem, each objective function is multiplied with its respective weight.Improvement in the response highly depends on the weight assigned.In weight-based GA (WBGA), weights are random; this random assignment of weights does not consider the priority of responses.Therefore, the optimal solution is not guaranteed.In random weight assignment methods, weights are assigned to each objective function and then combined them to make a single objective.However, although these methods are simple to use, weight selection is still a challenging problem.In this paper, the weight of each response was calculated in a systematic manner using entropy weighting of GRA as discussed in Section 4.1.

Optimization Variables and Objective Function
Equations ( 14)-( 16) serve as the objective function of NSGA-II.In this method, NSGA-II was the main technique, and GRA entropy was utilized to generate the weight of each response.NSGA-II was used to find four optimized milling parameters (spindle speed, feed rate, depth of cut, and width of cut) so that the best performance measures (MRR, SR, and ACE) were achieved without violating any of the constraints considered for this problem.For the simultaneous optimization of goals, the three objective functions were given as follows: Objective 1 = material removal rate; Objective 2 = surface roughness; Objective 3 = active cutting energy.The measured response must meet the constraints in the CNC milling process.Following the practical situation, the following constraints were considered.
Machine tool: CNC milling conditions must be within an allowable range of machine tool parameter values; Equation (18) explains this constraint for processing.
x min < x i < x max (i = 1, 2, 3, 4) ( 18) x i denotes milling parameters, F c represents the maximum cutting force, P max is nominal motor power, and η themachine efficiency.In the optimization procedure using NSGA-II, the following parameters were listed based on the study to get optimal solutions with low computational effort and have been set by trial and error through many repeats of simulations and also determined from the literature [39]: crossover probability: 0.95; distribution index for the crossover: 10; mutation probability: 0.01; distribution index formulation: 10; population size: 100.However, ACE increased as the MRR increased.On the other hand, surface roughness values also increased as the MRR increased.As seen in Figure 8a, region A consisted of a Pareto-front solution, which was good for reducing ACE and at the same time not causing an increase in R a values.The optimized values of cutting parameters and response values were: (N:1622.38;f :300.394;ap:0.4999;ae:13.7012;SR:2.022;ACE:40.2928) at point A. In addition, the Pareto optimal solutions in region B were good for getting maximum productivity.However, the ACE was higher in this region.The optimized milling parameters at point B were: (N:1200; f :320; ap:0.5;ae:15) and the corresponding response values were: (MRR:2387.5;ACE:64.60).Similarly, point C was a tradeoff between SR and MRR with optimized milling parameters at point: (N:1217; f :320.001;ap:0.3;ae:15) and the response values were: (MRR:2135.73;SR:0.66).Besides, other solutions can be selected after considering the desired criteria.A wide range of Pareto-optimum solutions demonstrated the applicability of this method in machining optimization that considers conflicting objectives.

Influence of Milling Parameters on MRR, ACE, and SR
It is clear from Figure 9 that for active cutting energy (ACE) and MRR, the width of the cut had the most significant effect, while for spindle speed, this effect was the least significant.Furthermore, GRC of ACE and MRR were in a linear relationship with the width of cut, so the highest value of the width of cut produced the lowest values of ACE and highest MRR.Similarly, for SR, the width of the cut also showed the most significant effect.From Table 2, the maximum SR values were found at a 5 mm width of cut, and the lowest SR value was found at 15 mm.
Extensive dispersion of Cu nanoparticles in the distilled water facilitated by cool high-pressure mist in workpiece and tool interface showed an excellent performance in reducing surface roughness.Furthermore, GRC of ACE and MRR were in a linear relationship with the width of cut, so the highest value of the width of cut produced the lowest values of ACE and highest MRR.Similarly, for SR, the width of the cut also showed the most significant effect.From Table 2, the maximum SR values were found at a 5 mm width of cut, and the lowest SR value was found at 15 mm.
Extensive dispersion of Cu nanoparticles in the distilled water facilitated by cool high-pressure mist in workpiece and tool interface showed an excellent performance in reducing surface roughness.The mist from nano-fluid assisted SQCL systems exhibited more efficient feeding into the machining zone than the conventional cooling method [40].This was because, in conventional wet cooling, lubricant splashes off by the rotating tool, but, mist sticks to the tool-workpiece interface and thus assists in improving cutting performance.In addition to this, the Cu nanoparticle in the SQCL provided a polishing effect, since an interacting force was then induced between the tool-particle interface, which is shown in Figure 9a [41].This induced interacting force travelled upon the cutting tool interacting surface and produced power on the surface.Hence, the surface energy of the workpiece was significantly reduced due to a high interacting force between the particle workpiece interface [42].
On the one hand, the breaking of nanoparticles required a certain energy for the removal of asperities from the workpiece surface, which generated a new surface with better finish.On the other hand, the process of making Cu nanoparticles transferred potential energy from a cutting tool, which later became kinetic energy of workpiece surface atoms.Therefore, a 300 ml/h flow rate provided a number of Cu nanoparticles, which transferred more kinetic energy to surface atoms and enhanced heat dissipation from the cutting zone.Figure 9b explains how nanoparticles behave as a spacer and effectively minimize the frictional effects between interacting surfaces.In the high speed-milling process, heat generated in the cutting zone altered elasto-hydrodynamic lubrication into boundary lubrication, and thus, this may create a possibility of rolling effect (Figure 9b) which could result in a reduction in coefficient of friction [43].
Moreover, the increase in the concentration of Cu particles in distilled water strengthened the thin protective film on the workpiece surface.Thiswas due to the strong chemical interaction between the newly generated surface and nanoparticles.This phenomenon increased the surface quality of the workpiece [42].The decrease in surface roughness was mainly due to the nanoparticles between the workpiece-tool interface.A good surface finish was mainly due to the excellent lubrication that slides the chips away from tool surface; better chip quality resulted in less surface roughness.Furthermore, a mist formation that entered the machining zone reduced the friction due to a ball bearing effect due to billions of rolling nanoparticles, and high-pressure spray removed the chips, ultimately leading to a better surface finish [44].A higher concentration of solid fine Cu nanoparticles of a large size provided effective lubrication under aggravated sliding conditions.Because these particles easily penetrated the surface between rubbing face, the original sliding friction divided into sliding as well as rolling friction, and that reduced the coefficient of friction [45].

Comparative Analysis and Validation
Traditionally, specialized handbooks and worker experience are used to determine the desired cutting parameters, but these sources provide starting values of cutting parameters, and these values are not necessarily suitable for optimal conditions.Initially, the levels of milling parameters were set at (N2, f2, a p 1, a e 1), and this combination is shown in Table 2 at experiment number 13.In this study, Figure 2 describes conventional multi-objective optimization, different from our purposed optimization method, only surface roughness, and material removal rate were selected as the optimization objectives in conventional optimization.Milling parameters (N = 2200 rev/min, f = 320 mm/min, a p = 0.5 mm and a e = 15 mm) were identified as the optimal setting for conventional multi-objective optimization.In order to perform the comparison between conventional and proposed multi-objective optimization method, we have also obtained optimal milling parameters from conventional optimization with two objectives.Table 6 shows the comparison of the experimental results of two optimization methods.Initially, the GRA based optimization results compared with the initial settings of milling parameters.For the comparison of two methods, material removal volume (MRV) was fixed to the maximum level, i.e., 2400 mm 3 .As a result, SR and ACE decreased from 2.85 µm to 2.29 µm and 269.6 kJ to 53.98 kJ respectively.Parameters settings based on the proposed optimization method resulted in a 79.9% reduction in energy consumption as compared to the initial setting.It is clear that the proposed optimization method can improve multi-response in the cutting process.The feasibility of the proposed method is indicated by the promising results obtained.
As shown in Table 6, the GRA based optimization resulted in a 20.7% reduction in ACE as compared with conventional objectives optimization.Moreover, along with traditional objectives, the proposed method also added a sustainable objective, i.e., active energy consumption, which promotes environmental benefits and saves energy.

Conclusions
With an objective to simultaneously reduce the energy consumption and improve the surface quality, this study presented an experimental and optimization approach to reduce energy consumption and improve surface quality without compromising the material removal quantity under the use of nano-fluid assisted SQCL.The optimization was performed as a 'multi-objective' using GRA and NSGA-II methods.The following concluding remarks are found:

•
The GRA based optimization method has revealed that a cutting speed of 1200 rev/min, a feed rate of 320 mm/min, a depth of cut of 0.5 mm, and a width of cut of 15 mm were the optimal cutting levels.It also found that the width of the cut hold had the most significant effect on response measurements.

•
Comparison of experimental results revealed that a lower spindle speed and higher width of cut were suitable for energy efficiency in the practiced milling process.The reduction of energy consumption was 20.7% as compared to the conventional optimization method.

•
In NSGA-II, three conflicting objective functions, i.e., MRR, SR, and ACE, were considered for the optimization problem.The optimal solution for milling process was found from Pareto-fronts.
The GRA converted the complex multi-objective problem into a single objective, and a maximum value of GRG represented the initial optimal cutting condition.

•
The results showed that NSGA-II gave an efficient result and similarly, GRA predicted results accurately with less than 5% error.

•
Comparison of experimental results revealed that lower spindle speed and higher levels of the width of cut were suitable for energy efficiency of machine tools.Hence, the proposed algorithm in this paper provides broader solutions for engineers to choose cutting parameters according to optimal results; therefore, the proposed method demonstrated applicability for manufacturing engineers working in the metal processing industries.
In conclusion, the implementation of the current study in the industry can reduce a significant amount of energy consumption.

Figure 1 .
Figure 1.Experimental setup for a small quantity cooling lubrication (SQCL) assisted milling process (a) schematic diagram (b) Experimental setup

Figure 1 .
Figure 1.Experimental setup for a small quantity cooling lubrication (SQCL) assisted milling process (a) schematic diagram (b) Experimental setup.

Figure 2 .
Figure 2.The preparation of copper based nanofluid.

Figure 2 .
Figure 2. The preparation of copper based nanofluid.

Figure 3 .
Figure 3.The effect of process parameters on process efficiency.

Figure 3 .
Figure 3.The effect of process parameters on process efficiency.

Figure 4 .
Figure 4. Experimental equipment used, showing: (a) power measurement using a Portable Power Cell PPC-3 (b) surface roughness measurement.

Figure 4 .
Figure 4. Experimental equipment used, showing: (a) power measurement using a Portable Power Cell PPC-3 (b) surface roughness measurement.

Figure 5 .
Figure 5.A framework for the energy reduction problem in the milling process.

Figure 5 .
Figure 5.A framework for the energy reduction problem in the milling process.

+ 1 . 17 )Figure 6
Figure6presents the relationship between the measured and predicted values of GRG.Measured values were close to predicted values, and the average percentage deviation was found to be 4.69, which indicated that experimental and predicted results were in good agreement.

Figure 6 .
Figure 6.Comparison of measured GRG and predicted GRG values.

Figure 6 .
Figure 6.Comparison of measured GRG and predicted GRG values.

Figure 7 .
Figure 7. Flow chart of the milling process optimization algorithm with GRA Coupled with theNon-Dominated Sorting Genetic Algorithm (NSGA-II).

Figure 7 .
Figure 7. Flow chart of the milling process optimization algorithm with GRA Coupled with theNon-Dominated Sorting Genetic Algorithm (NSGA-II).

2 .
Constraints and Parameters for the NSGA-II Algorithm

Table 1 .
Cutting parameters and their levels.

Table 1 .
Cutting parameters and their levels.

Table 2 .
Experimental design with measured values of responses.

Table 2 .
Experimental design with measured values of responses.

Table 3 .
Weight calculation for each response using Grey entropy weight.

Table 4 .
Grey relation coefficients (GRC) and Grey relation grade (GRG) values for each response.

Table 5 .
ANOVA table for the Grey relational grade (GRG).

Table 5 .
ANOVA table for the Grey relational grade (GRG).