numerical simulations of nano-structured surface boiling can be attributed to difficulties in descriptions of the micro-layer beneath the growing bubble and the nano-structure on the boiling surface. For the description of the micro-layer, its evaporation has not been included in most previous simulations [19; 20; 21]. The rapid evaporation of the micro-layer is considered to have a significant contribution to bubble growth and boiling heat transfer [22]. The micro-layer evaporation accounts for approximately from 40% to 60% of the total heat and mass transfer during the boiling process [23], meaning the micro-layer's existence cannot be ignored in numerical simulations. The establishment of the micro-layer model relies on the accurate experimental measurement of the rapid growth of boiling bubbles. [PERSON] and [PERSON] [24] first proposed a micro-layer model based on their experimental results, and the model was then adopted in [PERSON] and [PERSON]'s simulations [25]. In addition, [PERSON] et al. [26] also deduced a micro-layer model according to the velocity of the vapor-liquid interface. However, many speculations of the micro-layer thickness in those models directly influence the model's accuracy. With the development of the measuring technique, the micro-layer thickness distribution was reported by [PERSON] et al. [27; 28] using the laser extinction method. [PERSON]'s up-to-date micro-layer thickness equations have been widely adopted by recent numerical simulations [29; 30; 31; 32], and good robustness of the micro-layer thickness equations has been revealed in those studies. Modeling the nano-structures on a boiling surface is more challenging. The wickability of the super hydrophilic nano-structure will enhance the rewetting performance on the boiling surface, and the evaporation of the wicked fluid inside the nano-structure near the triple contact line also contributes to the bubble growth. However, a CFD method cannot establish fine modeling of the centimeter scale bubble and the nano-scale surface structure at the same time. Here, we provide a brief overview of numerical simulations of micro- and nano-structured surfaces boiling. [PERSON] et al. [33] used the level-set method to capture the gas-liquid interface and numerically simulated the nuclear boiling process on a micro-structured surface. They analyzed the bubble dynamics and heat transfer characteristics in each case. [PERSON] et al. [34] used the Lattice Boltzmann Method (LBM) to simulate the nucleate boiling process on the surface of two-dimensional micro-column structures. [PERSON] et al. [35] simulated boiling heat transfer on hydrophilic-hydrophobic mixed surfaces with 3D LBM; however, in their study, the effect of nano-structures only manifested in the influence on the wettability of the heating surface, while simulation and analysis were absent for the additional heat transfer caused by the nano-structures. The molecular dynamics (MD) simulation is an ideal tool to investigate the boiling process in nano-scale. [PERSON] et al. [36] studied the sub-cooled boiling on the different nano-structured surfaces at different degrees of wall superheat. [PERSON] et al. [37] analyzed the nucleation process on a platinum surface with/without the triangular nano-groove. [PERSON] et al. [38] carried out MD simulations for the boiling on a concave hemispherical nanostructure surface. It should be acknowledged that molecular dynamics research plays an important role in revealing the micro-mechanism of bubble nucleation and CHF formation on nano-structured surfaces. However, the biggest problem with MD methods is that the simulated time and spatial scales are extremely small, limiting their applicability to the macroscopic phenomenon of bubble boiling. Recently, several experimental investigations about boiling on superhydrophilic nano-structured surfaces have been carried out, and creative models about the wicked fluid inside the nano-structure have been proposed. [PERSON] et al.