Evaluation of subcooled MQL in cBN hard turning of powder-based Cr-Mo-V tool steel using simulations and experiments

Metal cutting fluids for improved cooling and lubrication are an environmental risk and a health risk for workers. Minimizing water consumption in industry is also a goal for a more sustainable production. Therefore, metal cutting emulsions that contain hazardous additives and consume considerable amounts of water are being replaced with more sustainable metal cutting fluids and delivery systems, like vegetable oils that are delivered in small aerosol droplets, i.e., via minimum quantity lubrication (MQL). Since the volume of the cutting fluid in MQL is small, the cooling capacity of MQL is not optimal. In order to improve the cooling capacity of the MQL, the spray can be subcooled using liquid nitrogen. This paper investigates subcooled MQL with machining simulations and experiments. The simulations provide complementary information to the experiments, which would be otherwise difficult to obtain, e.g., thermal behavior in the tool-chip contact and residual strains on the workpiece surface. The cBN hard turning simulations and experiments are done for powder-based Cr-Mo-V tool steel, Uddeholm Vanadis 8 using MQL subcooled to −10 °C and regular MQL at room temperature. The cutting forces and tool wear are measured from the experiments that are used as the calibration factor for the simulations. After calibration, the simulations are used to evaluate the thermal effects of the subcooled MQL, and the surface residual strains on the workpiece. The simulations are in good agreement with the experiments in terms of chip morphology and cutting forces. The cutting experiments and simulations show that there is only a small difference between the subcooled MQL and regular MQL regarding the wear behavior, cutting forces, or process temperatures. The simulations predict substantial residual plastic strain on the workpiece surface after machining. The surface deformations are shown to have significant effect on the simulated cutting forces after the initial tool pass, an outcome that has major implications for inverse material modeling.


Introduction
Machining operations are among the most common manufacturing processes used in industry. Therefore, even small changes in productivity, energy efficiency, and sustainability of machining operations have major influence in the global scale. Sustainability has become increasingly an important factor in manufacturing. One of the earliest mentions of sustainable manufacturing was in a paper by James Moore in 1977, where he envisioned the future goals of manufacturing organizations to be "less capital intensive, less energy consuming, more rewarding to workers, less demanding of raw materials and sustainable much further into the future" [1]. Since then, sustainability has become one of the focus areas of manufacturing research. Sustainability research builds upon Life Cycle Assessment (LCA), which sustainability has surpassed in publication volume at the beginning of the twenty-first century. Ong et al. (2001) [2] developed a Life Cycle Assessment tool for estimating the harmful emissions caused by material processing. Kumaran et al. (2001) [3] did similar research on environmental LCA for products. Duflou et al. (2012) [4] investigated the energy efficiency of production and their results show that there is not much that can be done in process mechanics level, but manufacturing system level optimization offers some possibilities for higher efficiency. Teece (2007) [5] considers the global manufacturing environment and concludes that companies aiming for sustainability require dynamic and flexible management. For more technical solutions for improved sustainability in machining, Jayal et al. (2010) investigate a sustainability index that considers the following criteria: (1) safety, health, and environment and (2) cost of machining, power consumption, and waste management. Their paper considers sustainable machining in three groups, i.e., dry machining, minimum quantity lubrication (MQL), and cryogenic machining, which all have their benefits and disadvantages [6]. Pušavec et al. (2010) analyzed LCA of different cooling strategies, flood, high-pressure cooling (HPC), and cryogenic cooling. Their results show that cryogenic cooling with liquid nitrogen has the highest potential for improved sustainability of machining superalloys [7]. Sivaiah et al. (2018) compared the machining of stainless steel with cryogenic, emulsion, and MQL cutting fluids against dry cutting conditions. Cryogenic coolant performed the best in all measured categories: cutting temperature, tool wear rate, and surface roughness. MQL was the second best option, third was flood cooling, and dry cutting was expectedly the worst [8]. A more fundamental investigation of the cooling capability of cryogenic coolants and flood emulsion is given in Pušavec et al. (2019) [9]. The tradeoff in optimization between sustainability and productivity requires more advanced predictive process models [6]. Regarding that goal, Avram et al. (2011) developed a multi-criterion optimization tool for sustainable machining. Their results show that MQL is the optimal solution for performance and sustainability [10].
In general, cutting fluids are used to improve tool life and surface quality during machining operations by reducing cutting temperature and friction. Typical cutting fluids are wateroil emulsions, with emulsifiers and additives to prevent corrosion and bacterial growth. The disposal of cutting fluids is expensive because of the toxic additives in the emulsions and the machine oils and heavy metals that are mixed with the fluids during the machining process. To improve the costefficiency and to reduce the environmental impact, the amount of cutting fluids applied need to be minimized and the toxic substances in the cutting fluids need to be removed. Both objectives can be achieved by using pure-vegetable oil with MQL delivery system. Vegetable oil is a sustainable, nontoxic, renewable, and biodegradable alternative to traditional cutting oils [11]. A comprehensive literature review by Sankaranarayanan et al. (2021) shows that bio-friendly cutting fluids are increasing their market share, which is motivated by the good lubricating properties of vegetable oils and growing environmental consciousness [12]. Khan et al. (2020) evaluated multiple cooling strategies on their case examples, regarding productivity and sustainability aspects. They used flood, MQL, and hybrid CryoMQL (HCM) coolants and calculated the energy consumption while machining titanium alloys, and the results show that the HCM method has the lowest sustainability impact and highest productivity [13]. Sudheer et al. (2014) investigated compressed refrigerated air, dry and Propylene Glycol MQL cooling mediums in a cutting of metal matrix composite (MMC) materials. Their results show that compressed refrigerated air improves the surface roughness slightly compared with dry cutting conditions and propylene glycol has significant effect by reducing the surface roughness by 50% compared with other mediums in the study [14]. It has also been shown that using MQL decreases the tool-chip contact length [15] and leads to reduced tool wear, cutting forces, and surface roughness, but the poor cooling effect is the major disadvantage of MQL [16]. On the other hand, in some cases, the lower cooling effect can be beneficial with regard to tool wear, for example in drilling, where the center of the tool has a higher wear rate with emulsion than with MQL [17]. The cooling effect of the MQL can be increased by using MQL combined with liquid N 2 or CO 2 , but these systems have a high initial cost. Khanna et al. (2021) review different cryogenic setups like hybrid cryo-MQL and high-pressure cutting fluids, etc. They conclude that the field is yet matured enough to give specific guidelines for selecting a cooling method for each case but investing in the development of cryogenic technologies is necessary for sustainable and economical machining [18]. A new and alternative approach is to use liquid N 2 in a heat exchanger to cool the MQL oil to subzero temperatures. This paper investigates a subcooled vegetable-oil MQL system patented by Accu-Svenska Ab [19] in hard machining of Vanadis 8 tool steel using 3D FEM simulations and experiments to compare the machining performance between MQL and subcooled MQL at −10°C. Cryogenic temperature has been defined to be below −150°C by National Institute of Standards and Technology [20]. Therefore, this paper refers to the Accu-Svenska system as subcooled MQL and other cryogenic MQL systems as CryoMQL. There have not been previous studies on machinability of Vanadis 8, and there are only very few papers investigating the Accu-Svenska CryoMQL system. In addition, the simulation methodology in this paper shows that using the approach proposed in this paper leads to better results than previously reported.

Literature review
The existing research in cryogenic coolants is focused on a machining of aerospace alloys like Ti-6Al-4V and Inconel 718, but some other work materials have also been investigated like tool grade aluminum and steel. Shokrani Shokrani et al. (2017), because their results show that LN 2 leads to highest tool wear and surface roughness but CryoMQL on the other hand leads to lowest tool wear and surface roughness. Adhesive wear was identified as the principal wear mechanism in all cooling conditions. The most pronounced chip serration was observed in cryogenic cutting conditions [22]. Interesting new approach was investigated by Olsson et al. (2021), where a LN 2 -cooling and induction heating was applied to a machining of pure tungsten and pure niobium. It was found that LN 2 and hybrid LN 2 with induction performed the best for machining the tungsten, but the conventional high-pressure coolant was optimal for machining niobium [23]. In this paper, the state-ofthe-art simulation practices of modeling cryogenic and subcooled cutting fluids are mapped in order to model the subcooled MQL investigated in this paper. In addition, publications investigating friction and wear models are reviewed to establish a method for modeling the lubricating effect of MQL in cutting conditions. Most of the simulations in the literature are done with Deform, AdvantEdge, or Abaqus using a high heat transfer coefficient between the tool and environment and workpiece and environment. In some cases, the cooling effect has also been simulated using CFD. All papers investigating tool wear used Usui's wear model. The modeling of friction in metal cutting is typically done in dry conditions, but in this paper, the focus is on MQL where the friction is in pure-oil environment. Many papers report an increasing friction coefficient with increasing cutting speed (sliding speed), which is opposite behavior from what is commonly observed in dry cutting conditions.

Simulating cryogenic cutting fluids and MQL
Salame et al. (2019) simulated 2D orthogonal a cutting of Ti-6Al-4V using FEM model in Abaqus and use the cutting simulation temperature distributions as input to CFD model of cryogenic cooling (direct application of liquid CO 2 ) in Ansys. Their results show that cryogenic nozzle placement has significant effect on the cutting temperature and the position of the maximum temperature, somewhat counterintuitively so that farther nozzle position (4-mm vs. 8-mm distance) led to lower temperature and the maximum temperature moved from the secondary to the primary shear zone [24]. show that the heat transfer coefficient is strongly dependent on the temperature difference between the surface and gas, being 75 kW/m 2 K at 200°C temperature difference and reducing exponentially; ΔT = 160°C, h = 15 kW/m 2 K, and about 1.5 kW/m 2 K at temperature difference below 50°C [26]. Sun et al. (2015) investigated a machining of Ti-5553 alloy with MQL, flood, and cryogenic cooling using experiments and FEM simulations with Deform 2D. Cryogenic cooling was simulated using an additional heat exchange window under the tool flank, with environment temperature of −184°C and convection coefficient of 10000 kW/(m 2 K). The FEM simulations are in good agreement with experiments done using cryogenic cooling. Cryogenic cooling reduced cutting forces by 30% compared with flood and MQL [27].  followed the work of Sun et al. (2015) closely, but with more experimental data, more reliable thermal imaging, and more detailed descriptions of the simulations. The simulations are relatively accurate in comparison to experiments with 10-30% error with respect to cutting forces, temperature, and chip morphology. CO 2 coolant reduces cutting temperature by 15-52% compared with dry cutting [28]. Imbrogno et al. (2017) simulated Ti-6Al-4V cutting using Deform 3D and TANH model. LN 2 cryogenic cooling was simulated using local heat transfer windows with −196°C temperature and heat convection calculated by method implemented by Umbrello [29] and proposed by Astakhov [30]. The results show good agreement with experiments regarding cutting force and temperature and less so for feed force [31]. Davoudinejad et al. (2015) simulated cryogenic cooling of an orthogonal cutting using Third Wave Systems AdvantEdge FEM software. The heat transfer coefficient was calibrated to 2000 kW/m 2 K by testing with pre-selected values between 46.75 and 5000 kW/m 2 K. Simulations were first verified using dry machining conditions, which were in good comparison with the experiments having error less than 10% in all control variables: cutting forces, feed forces, and chip thickness. The simulations with cryogenic conditions also produced low error. The simulations and experiments show that the cutting forces and stresses are higher with cryogenic cooling than with dry cutting while chip thickness and serration frequency reduced with cryogenic cooling [32]. Mishra et al. (2019) simulated a machining of Ti-6Al-4V with a textured tool in cryogenic cooling conditions using AdvantEdge FEM software. The simulations show that realistic cutting temperature was achieved using 3000 kW/m 2 K heat convection coefficient at tool-chip interface. Overall accuracy of the simulations was not optimal [33]. Dix et al. (2014) modeled drilling assisted with LN 2 cryogenic coolant. Drill with internal coolant channels was simulated with Deform 3D. The heat transfer coefficients used in the simulations varied by work temperature 23.3 kW/m 2 K for −181°C and 46.8 kW/m 2 K for 650°C, while ambient temperature was set to −191°C. The simulations showed good agreement with the experiments [34]. Rotella et al. (2014) investigate a machining of aluminum AA7075-T651 with cryogenic cutting fluid using experiments and simulations. Dry conditions lead to higher cutting forces and lower workpiece hardness, which is credited to reduced thermal softening and grain growth during cryogenic cooling. The simulations are done with Deform 3D using 20 W/m 2 K heat convection coefficient between the tool and environment and workpiece and environment. The tool-chipinterface heat conduction was set to 55000 kW/m 2 K to enable reaching the thermal steady state quickly. A modified Johnson-Cook model was used for the flow stress modeling with room temperature set to −182°C. The Johnson-Cook model was modified by replacing the yield strength equivalent parameter with Hall-Petch relation, which is explained in Farrokh and Khan (2009) [35]. Model parameters were taken from Curle and Govender (2010) [36]. Cryogenic cooling effect was simulated using a local heat transfer window with the convection coefficient initially set to 0.2 kW/m 2 K, and the grain size was predicted using Zener-Hollomon equation. Simulations predicted cutting force components with decent accuracy and grain size and hardness with relatively good accuracy [37]. Yu Su (2015) investigated cryogenic MQL in a machining of AISI H13 steel using 3D FEM simulations. The MQL was assisted with refrigerated compressed air (RCA) at different temperatures [−10°C; −140°C]. The simulations were done with Deform 3D using the Johnson-Cook flow stress model and a local heat transfer window for the cryogenic effect. The heat transfer coefficient was 100 W/m 2 K. The results show that the optimal RCA temperature was dependent on cutting speed and was −30°C for 47 m/min and −50°C for 104 m/min. Cutting temperature and cutting forces were decreased significantly while using the RCA aided MQL. The accuracy of the simulations was not adequately evaluated [38]. Shi et al. (2019) used a hybrid CFD and FEM approach to simulate LN 2 cryogenic cooling in the machining of Ti6Al4V titanium. They use a Voce power law for flow stress modeling because the Johnson-Cook model gives invalid values below the reference temperature. They used a CFD to determine the convection coefficient between the LN 2 spray and the tool/workpiece, which was estimated in range from 24 to 50 kW/m 2 K. Their simulation results were in relatively good comparison with the experiments: 7% error regarding the chip temperature and 10-24% error in the cutting force components. They explain the high force error with a ploughing effect in the feed direction, friction model errors, and unaccounted tool wear [39].
Hard machining has some special characteristics that arise from the geometry and material properties of the workpiece: the cutting depth and feed are small; therefore, the effect of friction and work material residual strains has higher effect than in cutting with larger chip volume. The strain hardening behavior of hardened materials is steeper, and regular Johnson-Cook model is not optimal for capturing the shape of the strain hardening curve. The cutting force proportions are different from regular cutting, since the perpendicular or passive force becomes the largest cutting force component in hard machining, instead of the main cutting force. The temperature of the workpiece does not increase much since most of the heat is removed with the chip, and the heat generated at the tertiary deformation zone is small compared to the volume of the workpiece because of the small chip volume. In existing literature, there are only few publications on FEM modeling of hard machining. One such work is done by Umbrello et al. (2008), where they developed a flow stress model for AISI H13 based on the J-C model with hardness dependency. The model includes additional terms in the strain hardening equations, to take into account the hardness variations in the work material. The model gives relatively accurate (<10% error) simulation results regarding cutting forces and chip morphology in four cases of machining AISI H13 at different hardness . Only the feed forces were overestimated by 11-16% in all cases. The simulations were done with STFC Deform [40]. Much of the data and methods used in Umbrello et al. (2008) were inspired by Yan et al. (2007) [41]. Saez-de-Buruaga et al. (2019) present a flow stress model for hard machining of 42CrMoS4. Their work illustrates well the steep strain hardening and leveling of the flow stress in hardened steels, and why J-C model is not optimal for fitting such behavior. Their model is based on yield stress and saturation stress and filling the gap between with exponential relation. The model fits well to the presented data; however, the simulated cutting forces have 10% average error and contact lengths 29% average error [42]. To illustrate the poor applicability of Johnson-Cook model to strain hardening behavior of hardened steels, the Saez-de-Buruaga model and the best fit of Johnson-Cook model to the same data set are shown in Figure 1. To overcome this issue in this paper, the flow stress is modeled using tabular data.  Yadav et al. (2015) simulated MRR (material removal rate) and flank wear in a dry turning of Inconel 718, using Deform 3D. They used somewhat an unorthodox approach with an unusually long workpiece model with a coarse mesh, whereas it is more common to use exaggerated wear model parameters to simulate longer cutting distance. The simulations were built with the Johnson-Cook flow stress model, a Coulomb friction, and Usui's wear model. The simulations show good agreement (<7%) with experiment regarding MRR and flank wear, but the paper does not report the cutting forces, so the simulation accuracy might be specific to MRR and flank wear in the local cutting conditions used [43]. Binder et al. (2015Binder et al. ( , 2017 used Deform 3D and Usui's wear model to simulate an in situ tool wear by updating the tool geometry while machining AISI 1045 with coated and uncoated carbide tools. The results show good agreement with the experimental tool wear, but the paper does not compare the chip morphology or the cutting forces [44,45]. Schulze and Zanger (2011) simulated tool wear in 2D orthogonal a machining of Ti-6Al-4V using Abaqus and Usui's wear model. The WC tool was modeled rigid and the contact between the tool and workpiece was modeled with Coulomb friction. Ti-6Al-4V was modeled with Johnson-Cook flow stress model J-C failure model. The results show good comparison with experimental tool wear and reasonably good comparison with specific cutting force and chip temperature, both being within ±15% [46]. Yen et al. (2004) developed a wear simulation with Deform 2D and Konti-Cut user script that updates the tool geometry based on the Usui's wear model. Work material is AISI 1045 and tool is uncoated carbide. Their results showed that the adopted wear model was not able to capture wear at the tool edge, i.e., flank wear, but the wear was concentrated on the rake face near the end of tool-chip contact length [47]: "Besides, it is known that high wear rates are normally seen at the initial stage of a cut for a sharp edge due to the weaker edge strength and possible micro-chipping." Laakso et al. (2017) showed that the plastic deformation and creep at the tool edge contribute to the high initial wear rate [48] and the initial wear was shown to contribute a major proportion of the feed force in [49], and the same conclusion was made later in Peng et al. (2019) [50]. In order to circumvent the issue with flank wear not progressing as fast as the wear at rake face, Yen et al. set a 60-μm initial flank wear land with 3°inclination and repeated the simulation. This simulation showed defined flank wear progression over time, the V B growing to 120 μm in 6 min. The paper did not report any validation of cutting forces or chip morphology. Tamizharasan et al. (2012) optimized tool geometry regarding tool wear in a machining of AISI 1045 using Deform 3D with Usui's wear model and analysis of variance to investigate the effect of geometry parameters. Their results did not reflect the reality of experiments with average cutting force error of 82%, and the wear depth from simulations was 10 −4 mm and 0.537 mm in experiments [51]. Szczotkarz et al. (2021) investigate tool wear in machining of 316L stainless steel with MQL and dry cutting conditions. Their results show that MQL reduced tool wear by 18% compared to dry cutting. Their results also show that MQL can very well reduce adhesive wear, but for example notch wear is not significantly affected by MQL compared to dry conditions [52].

Friction and lubrication
Metal cutting is typically modeled in dry cutting conditions for simplicity and much of the research work is focused on modeling the friction in dry conditions. In dry cutting conditions, friction is shown to reduce with increased cutting speeds, shown in, e.g., Mane et al. (2020) [53]. This paper investigates cutting with MQL and subcooled MQL, and therefore, the friction is modeled in lubricated environment, in which the behavior of friction is significantly different from dry condition. Claudin et al. (2010) investigate the effect of oil on friction at tool-workpiece interface. They provide insightful review on friction and lubrication in metal cutting, bringing the controversial issues to discussion: (1) if cutting lubricant can penetrate the tool-chip interface because of the high contact pressure, (2) are tribological tests like pin-on-disc method representative to metal cutting conditions, and (3) different approaches for friction testing developed specifically for metal cutting. The authors present a tribometer that they use to test friction in dry and oil-lubricated machining conditions for AISI 4140 workpiece and TiN-coated carbide tool. The results of the tribometer test show an exponentially decreasing friction coefficient μ∈[0.6;0.2] with increasing sliding velocity for dry condition and the static friction μ≈0.1 in lubricated conditions. Their results also show that oil does indeed penetrate the contact interface even at very high~2.5 GPa contact pressures, objecting the notion, e.g., in [54]. The oil needs to be continuously reapplied since it stays in the interface only for a short time (~1 s) varying with sliding speed. Adding oil to a dry condition starts to affect the contact after 0.2 to 0.6 s [55]. Pottirayil et al. (2010) experimented with paraffinic, naphthenic, and aromatic oil in metal cutting, using a tribometer with the cutting tool and a friction pin submerged in oil. Their results show a steady friction coefficient of 0.1 and peak friction increasing between 0.2 and 0.6 with increasing sliding speed (0.15-0.38 m/s). They conclude that friction testing for metal cutting must be done by in situ experiment; otherwise, the environmental contamination and oxidation will affect the results [56]. Cassin and Boothroyd investigated lubrication in a machining of copper, where they observed a clear increase in the friction angle with increasing cutting speed while using various cutting fluids and constant friction angle with dry cutting [57]. Campen et al. (2012) dedicated their research efforts on investigating the increase of lubricated friction with increasing sliding speed. Their paper scope is limited to organic lubricants in AISI 512100 ball-disc contact. Their results show that the friction coefficient is small <0.18 in all cases, and a typical trend is increasing friction with increasing sliding speed ∈[10 −7 to 10 −2 m/s]. Their conclusions suggest that the small friction at low sliding speed that increases with speed is caused by molecular behavior of the oils: sliding between ordered close-packed monolayers [58]. Faverjon [60]. In addition to the sliding speed, the surface finish and the contact temperature are found to affect the friction coefficient [61].

Materials and methods
The subcooled MQL is evaluated using cutting experiments and simulations. The work material is Vanadis 8, the tool is a cBN insert, and the cutting experiments are done with a lathe. The simulations are done with Deform 3D FEM software.

Cutting experiments
The cutting experiments were done using an EMCO Turn 365 CNC lathe with Kistler 9275 force measurement system and Accu-Svenska's prototype subcooled MQL system. The cooling and lubrication are delivered only at the rake face, in front of the tool-chip interface. The coolant nozzle placement is shown in Figure 2. Tool wear was measured in 50-mm cutting length intervals with Zeiss Discovery V20 optical microscope and with FEI/Philipps XL-30 SEM. The work material is a powder-based tool steel from Uddeholm, commercially named as Vanadis 8. The alloying elements of Vanadis 8 are presented in Table 1. The material is machined in hardened state at~62-64 HRC. Sandvik cBN 7125 grade insert was used for all experiments. The cutting parameters were selected based on Sandvik recommendations and initial cutting tests. The cutting parameters are cutting speed v c ∈ [50,75,100 m/min], cutting feed f = 0.15 mm/rev, and cutting depth a p = 0.15 mm. The experiments were done with MQL and subcooled MQL that are compared with simulated results. The experimental work done parallel to this paper shows that the −10 degrees is the optimal temperature for this Accu-Svenska Cryo MQL setup. Going below −10 degrees does not improve the tool wear. Tool wear is visualized using Logit model (Eq. 1) that is capable of capturing the wear progression over time, instead of only the flank wear length V B at the end of tool life. Logit model is described in detail by Laakso and Johansson (2019) [62]. In this paper, the t max is solved from Taylor's tool life Eq. 2 [63].
Logit model, where the parameters a, b, and m are for calibrating the wear progression shape Taylor's model, where n is the wear rate exponent and C is the wear coefficient

Simulations
FEM simulation have been used in metal cutting research intensively for nearly half a century. Alternative methods have been proposed like BEM, PFEM, or analytical models like one proposed by Lazoglu and Islam (2012) [64]. The simulations in this paper are done with Deform 3D that is a FEM solver based on Lagrange formulation and implicit time integration. The simulation setup is presented in Fig. 3. Time step for the simulations is 1×10 −6 s, the element size is shown in Figure 4 for the workpiece and in Figure 5 for the tool, and the workpiece is meshed with an adaptive remeshing to avoid excessive element distortions. The tool is modeled as elastic and workpiece is elastic plastic. The simulations run for two consecutive cutting revolutions to capture the effect of the previously deformed workpiece layer on the cutting mechanics and possible temperature buildup in the workpiece. The effect of the previously deformed workpiece surface has been shown to affect the flow stress of the work material significantly in Laakso (2020) [65].

Tool geometry
The tool was modeled with Creo Parametric 6.0.2.0 using the insert geometry values provided by Sandvik and microscope measurements. The chosen strategy of modeling the tool wear in advance was inspired by the existing research on the progressive modeling of the tool wear, which has been shown to be not ideally accurate and unproductive regarding the CPU hours. By modeling the tool wear in advance, based on the wear geometry obtained from the experiments, the error of inaccurate tool wear prediction in simulations is minimized and the simulation solving time is faster. The modeled tool geometry and mesh are shown in Figure 5. The tool was modeled with wear geometry obtained from experiments at 80-μm flank wear and corresponding crater wear, shown in Figure 6.

Friction
Friction in the simulations was used between the tool and the workpiece. The friction model is shear friction, presented in Eq. 3, where k is the shear yield stress and m is the shear friction coefficient. The shear friction coefficients at different cutting speeds were determined using inverse analysis (presented in Section 3.3) and initial assumption based on the literature review suggesting that friction increases with increasing sliding speed in lubricated contact conditions. The determined shear friction coefficients were found to concede with the velocity dependent behavior found in literature review, as shown in Figure 7. Shear friction model, where τ f is the shear friction

Heat transfer
Heat transfer was modeled between the workpiece and the tool, between the tool and the environment, between the workpiece and the environment, between the tool and the MQL spray, and between the workpiece and the MQL spray. The heat transfer between the workpiece and the tool is modeled as a constant heat transfer coefficient of 40 kW/m 2 K, which is commonly used in FEM studies on tool-chip interface. The convection coefficient from tool and workpiece to environment is 20 W/m 2 K. The convection coefficient from the tool and the workpiece to a constant −10°C subcooled MQL spray is 10 kW/m 2 K, which is in line with Hribersek et al. (2017) and Shi et al. (2019) considering that the coolant in this research is not LN 2 but subcooled vegetable oil [26,39]. The subcooled MQL spray was modeled using an environmental window shown in Figure 8 that follows the tool. The subcooled MQL spray for the workpiece was similarly defined with an environmental window, shown in Figure 9.

Mechanical properties of Vanadis 8
Mechanical properties of Vanadis 8 hardened to 62.9 HRC are modeled using tabular flow stress data that was compiled from tensile testing data provided by Uddeholm, shown in Table 2, a thermal softening curve determined from temperaturedependent modulus of elasticity data calibrated with inverse methods and a rate hardening multiplier from the Johnson-Cook model that was calibrated with an inverse method. The flow stress strain dependency is shown in Figure 10, thermal softening is shown in Figure 11, and the Johnson-Cook rate hardening in Eq. 4. The rate hardening parameters are 2.39×10 −4 s −1 for the reference strain rateε ref and the rate hardening constant C is 2.2392×10 −2 . Rate Hardening Multiplier

Inverse calibration of material model parameters
The inverse calibration of the material model parameters was done with 2nd degree polynomial response surface method (Eq. 5) using a cutting speed, a friction coefficient, the maximum flow stress, a thermal softening modifier, and the rate hardening coefficient as variables, and the output criteria were the force components. Crosstalk between the parameters was included by multiplying the individual polynomials like in Eq. 6. A total of 25 calibration simulations were performed to find parameters that gave acceptable error in the final simulations. The response surface itself was optimized after every additional calibration simulation and after 25 simulations, the response surface error was 5%. The final optimal output of the response surface for the maximum flow stress, the rate hardening coefficient, and the thermal softening modifier and friction are given in the above material data sets. The response surface parameters that gave the optimal values are shown in Table 3. Polynomial response surface for "x" where "x" denotes any of the input variables and "i" denotes the force component Total crosstalk effect of all variables for force component "i"

Results
The results of the experiments and simulations are presented for tool wear, cutting forces, thermal behavior, chip morphology, and residual plastic strain on the workpiece surface.

Tool wear
Tool wear data is presented as data points and corresponding Logit-model fit (presented with line) shown in Figure 12. The difference between subcooled MQL and regular MQL is small, but one might see a slight tool life benefit in favor of the subcooled MQL. The Logit-model parameters and Taylor's model parameters are shown in Table 4.

Cutting forces
Cutting forces from subcooled MQL experiments and simulations are shown in Figure 13, Figure 14, and Figure 15. The experimental forces are extracted from the measurement date from the time interval where the tool wear was equivalent of V B 80 μm. Total average error of the simulated cutting forces on the second revolution compared with experimentally determined values is 2.6% and maximum error is 5.2% at 100 m/ min regarding the passive force. The same errors for the first revolution are 28% and 42%, which clearly indicates that the previously deformed layer has significant effect on the flow stress properties of the work material, and thus, the cutting forces. The results for the regular MQL experiments and simulations are shown in Figure 16, Figure 17, and Figure 18. The obtained forces do not have any significant difference compared with the subcooled MQL forces. The simulations are also in good agreement with the experiments, average error being 3% and maximum individual error 5.1%.

Chip morphology
The chip shapes from simulations are presented in Figure 19. The chip formation is continuous, as observed in the experiments as well. The chip has a distinct curl that also was noticed in the experiments, as shown in Figure 20. The simulated chip width at 100 m/min cutting speed with subcooled MQL is 0.58 mm while the corresponding chip width from the experiment is 0.537 mm, so there is an acceptable 8% error. The chip width also increased in simulations after the first cutting

Thermal effects
The simulated workpiece surface temperature after one revolution with subcooled MQL at 100 m/min is 38.1°C and with MQL 41.3°C. The maximum temperature at the chip was 715°C for both cases. The tool temperature after the tool-chip engagement ended decreased slightly faster with the subcooled MQL compared with regular MQL, and the maximum tool temperature was 176°C with subcooled MQL and 178°C with regular MQL. Figure 21 shows the measurement points where the residual plastic strain was extracted from the simulations. Figure 22 and Figure 23 show the extracted data from simulations with  [71]. The difference between the subcooled MQL and regular MQL is relatively small. The strain increases by about 50% even after the second tool pass. One additional set of simulations was done to evaluate the residual strain formation with regard to the cooling-lubrication effect. In these simulations, the subcooled MQL windows were extended to the flank side of the tool, thus simulating two MQL nozzles, one on the rake face and other on the flank face. The simulation results are shown in Figure 24, where the residual stains from the subcooled MQL, regular MQL, and the dual-nozzle subcooled MQL at 100 m/min are compared. It can be seen

Discussion and conclusions
This paper investigated the subcooled MQL compared with regular MQL with simulations and experiments.
The following conclusions were made: Compared with regular MQL, subcooling the MQL to −10°C has a small beneficial effect on the cutting force, tool wear, and cutting temperature, observed with experiments and simulations. One of the reasons why there is no more significant difference is the high abrasivity of Vanadis 8, which leads to high abrasive wear of the tool, which is not sensitive to temperature. Optimizing the delivery of the subcooled MQL is expected to improve the tool wear, since the direct delivery of the lubricant at the flank face of the tool should improve the contact conditions. Also, the surface integrity of the machined surface is expected to improve with better lubrication and cooling at the flank face that was already predicted with the trial simulations done in this paper. The lubricative effect of the MQL was investigated with friction coefficient. The simulations suggest that the friction coefficient increases with increasing cutting speed. The literature review supports this behavior and the possible reasons for this are the change in the behavior of the oil in high pressure and temperature and inadequate volume of the lubrication oil flow, which both have more significant effect when the cutting speed is higher. The second nozzle is expected to improve the lubrication even at higher cutting speeds.
Simulations show that the second tool pass over the workpiece increases cutting forces because the surface layer of the workpiece is deformed during the initial tool pass, which can be well observed from the residual strain diagrams above. This has major implications for inverse determination of material properties: if the inverse method is used with simulations that have only one cutting pass, and the simulated workpiece material is not previously deformed, the inverse model gives overestimated values for flow stress. The researchers in the field often disregard this notion by stating that a pristine workpiece was used in the experiments, but do not realize that in every cutting operation, the surface is not pristine after the first revolution of the tool or the workpiece. Only broaching is capable of a long cutting pass with a pristine material, but even with this operation, the workpiece surface is deformed during the second pass, and the same error exists for any consecutive tool passes. Another notable fact is that if the simulations are modeled with real tool wear geometry and Figure 22 Residual plastic strain after 1st and 2nd tool pass with subcooled MQL Figure 23 Residual plastic strain after 1st and 2nd tool pass with regular MQL multiple tool passes, the ever-existing feed force error disappears. This paper provided further evidence of the notions about the importance of tool geometry, the initial wear of the tool, and surface properties of the workpiece made in Laakso et al. (2018Laakso et al. ( , 2019Laakso et al. ( , 2020 [48,49,65]. The simulation model needs to be further improved by calibrating the heat transfer coefficients using experimental setup with IR camera.
Author contribution Sampsa Laakso collaborated in the design of experiments, developed the simulation models and inverse analysis routines, did the literature review, wrote the manuscript, and analyzed the results.
Dinesh Mallipeddi did the design of experiments, conducted the experiments and microscopy, and participated in the writing process and analysis of the results.
Peter Krajnik acquired the funding, provided the facilities for the research, and participated in the writing process and quality assurance.
Funding Open access funding provided by Chalmers University of Technology. This work was performed under the sponsorship of the "Strategic vehicle research and innovation (FFI)" program (project SUSTAIN-CRYO, 2017-03057)-a partnership program run jointly by the Swedish state and Swedish automotive industry.
Data Availability All data and materials are presented in the paper Code availability The FEM solver used in this work is Scientific Forming Technologies Deform, which is commercially available software.

Declarations
Ethics approval This paper does not include any ethical considerations.
Consent to participate Not applicable as no humans or animals were the object of the research.
Consent for publication All the authors agree on publishing this work in the International Journal of Advanced Manufacturing Technology.

Competing interests The authors declare no competing interests.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.

Figure 24
Comparison of residual plastic strain after 1st and 2nd tool pass at 100 m/min