## 1 Introduction  
One may outline two issues when modeling heat transfer at a real scale. The first is associated with theoretical modeling that requires empirical data in terms of empirical correlations for calculation of the heat transfer coefficient, and the second is the application of numerical simulation software that is, when mass transfer is considered, very expensive in terms of the needed computational resources and time consumption when a fundamental (basic) approach is utilized in obtaining the numerical solution. To mitigate this, we used a lumped heat conduction model in the estimation of the heat transfer coefficient based on temperature distribution in a material and the corresponding heat conduction time, together with the novel mass transfer model that is used in conjunction with the frozen turbulence modeling approach. Within the scope of fundamental issues related to the boiling phenomena and correspondingly the associated heat and mass transfer, a single bubble dynamic has been experimentally studied in [PERSON] [1]. The downward-facing heating element was used to generate the bubble in a subcooled Perfluorhexan liquid (FC-72), and the bubble growth characteristics were examined among other features of the flow in the vicinity of the bubble. Shortly thereafter, the ISO 9950:1995 [2] standard probe made of Inconel 600 nickel alloy was reported in [PERSON] [3]. The oil and water quenchants at different temperature levels were examined using inverse heat transfer analysis (IHTA). Furthermore, at a fundamental research scale, the saturated and subcooled film boiling around a cylinder was studied both theoretically and experimentally in [PERSON] et al. [4]. The heat transfer surface was partitioned in the bottom, vertical, and top, and by application of different correlations for the calculation of the heat transfer coefficient, the heat flow rate for the entire cylinder was obeyed. The boiling issues related to the immersion quenching process were studied experimentally and numerically (finite element method, FEM) in conjunction with the inverse heat transfer analysis (IHTA) in [PERSON] [5]. The author used different geometries and different immersion velocities in studying rewetting zones and multipoint sampling of temperature data in the estimation of the heat transfer coefficients via IHTA. The enhanced heat transfer due to the existence of bubbles in the vicinity of the vertical heated surface has been reported within the experimental work described in [PERSON] et al. [6]. The noticeable increase in the heat flux was noticed in the presence of the bubbles in contrast to the pure, single-phase, bulk fluid. The heat transfer coefficient during the quenching process was determined by [PERSON] et al. [7]. The authors used a Lee-based mass transfer model in a mixture model formulation, with the dispersed phase being modeled using a bubble crowding approach. The irrelevance of imposing the heat transfer coefficient, as is mandatory in the structural thermal analysis of quenching using IHTA-generated heat transfer coefficient data, is explicitly stated in [PERSON] [8]. The variational multiscale (VMS) approach, which is claimed to be a large eddy simulation (LES) analog within FEM modeling and further explained in [PERSON] [9], with the well-established level-set method (LS), which is a direct numerical simulation (DNS) of an interface motion (not of turbulence, as clearly noted in [PERSON] [10]) method known for efficient interface curvature estimation, has been applied in the numerical modeling of the quenching process. Going back to the application of the inverse heat conduction approach, the heat transfer coefficients, together with the cooling rates, were estimated by [PERSON] et al. [11]. The immersion quenching process using water-liquid and polymer-based solutions at different temperature levels has been considered by the authors.