Numerical Study on Spreading and Vaporization Process of Liquid Nitrogen Droplet Impinging on Heated Wall

: Micro-structured surfaces can affect heat transfer mechanisms because of enlarged speciﬁc surface areas. However, employing the Leidenfrost effect during liquid nitrogen (LN 2 ) droplet cooling of a heated micro-structured surface possessing a ﬁn with a spacing much smaller than the diameter of the droplet has not yet been explored. In the present work, a direct numerical simulation (DNS) is carried out to investigate heat transfer mechanisms of the LN 2 droplet, whose diameter is sufﬁciently larger than the structured spacing of ﬁn, impinging on a micro-structured surface with variable velocities. For a comparative study, a smooth surface is also employed in numerical simulations. The spreading mechanisms and vaporization behavior of the droplet along with liquid ﬁlm morphology at various conditions are investigated. Results show that a smaller ﬁn size inhibits LN 2 in entering into the grooves between the ﬁns and left the surface untouched by the droplet completely, and eventually, a thinner liquid ﬁlm is spread out in contrast to the smooth surface. Notably, at a low Weber number, the droplet can be shrunk or even rebounded away from the wall after impinging on the wall. The fastest vaporization behavior for both surfaces, namely smooth and micro-structured, is obtained at a Weber number of 180. Additionally, an effective heat transfer upon the micro-structured surface is observed at a low impinging velocity of the droplet.


Introduction
Spray cooling with liquid nitrogen (LN2) has substantial potential for application in the fields of high heat flow electronic equipment cooling [1,2], cryogenic wind tunnel [3], cryogenic surgery [4], and food transportation [5] due to its low evaporation temperature, accessibility, relatively low cost, high safety, and stability.The spray cooling of LN 2 atomizes the jet into tiny droplets that collide with the surface and remove heat through vaporization.Therefore, it is essential to elucidate the mechanisms of spreading and vaporization of the droplets impinging on the surface from a microcosmic view.
Most of the literature on spray cooling, droplet impingement, and vaporization has been focused on room temperature fluids as the working medium [6][7][8].Liu et al. [9][10][11] demonstrated that increasing pressure difference can increase the uniformity of droplet distribution and reduce particle size, though the far-field spray cone angle was decreased significantly with the increasing variation in pressure.Numerical simulations have been widely used to study jet atomization in order to obtain insight into the internal flow mechanism [12].Chen et al. [13] established a numerical model coupling the algebraic interfacial area density (AIAD) model and homogeneous relaxation model (HRM).They analyzed the internal mechanisms of LN 2 sprays such as flow pattern, pressure distribution, velocity distribution, and complete core length based on a two-fluid numerical model.Wang et al. [14] developed a numerical model to predict the heating, boiling, and flashing characteristics of LN 2 droplets and the spray patterns during the cooling process to optimize operating conditions.Ruan et al. [15] confirmed that reducing the diameter of the droplet can be a more effective approach than increasing the velocity of the droplet to improve vaporization.
An approach of passive improvement in heat transfer can be a favorable way to enhance heat transfer.The fluid dynamics change with respective surface morphology, and the corresponding heat transfer characteristics can be influenced as well.
Yao et al. [16] analyzed the static and dynamic wettability of the nanoscale hydrophobic surface, and found that the micronano surface can increase the droplet rebound height, which depends on the static contact angle.Yee et al. [17] studied the influence of groove direction, impact velocity, and surface temperature on the shrinkage dynamics of the droplet through high-speed imaging.The oscillating airflow under the droplet caused by grooves resulted in the temperature-induced anisotropic shrinkage.Qin et al. [18] reported that droplets are more likely to splash by hitting a wet surface.Zhan et al. [19] revealed a generation process of droplets when an upward spray impinging at the downward solid surface and the proportion of splashing droplets was increased with the increase in Weber number.The droplet size was proportional to the length of the most unstable Rayleigh-Taylor wave.
In general, most of the current research on LN 2 spray cooling has been based on the analysis of macroscopic flow characteristics of the spray field and the effect of the microstructured surface on spray cooling based on the experiments of ambient fluids.There are a few literature studies on the cooling effect of a single LN 2 droplet from a microcosmic view.
Among the reported simulation methods, direct numerical simulation has been applied to combustion, reaction flow, and spray processes by virtue of its high accuracy and the ability to provide in-depth flow information of each instantaneous flow field [20][21][22][23].In this work, the direct numerical simulation is used to analyze the spreading process of LN 2 droplets on smooth surface and micro-structured surface.The vaporization rate, residual radius, variation in temperature of LN 2 droplets at the different velocities, and respective surface morphology are quantitatively studied.The effects of Weber number and initial velocity on the collapse and vaporization of the droplets are discussed.Here, the reported mechanisms related to the microscopic flow and laws can ensure the credibility of the accurate establishment of the macroscopic multiphase flow model.

Numerical Model 2.1. Physical Model
The geometric model for the LN 2 droplet in free-falling state and impinging on a copper surface is shown in Figure 1.The LN 2 droplet with the initial diameter (d 0 ) of 2 mm width was considered at a distance (H 0 ) of 4 mm between the solid surface and center of the LN 2 droplet.The initial temperature (T 0 ) of LN 2 droplets was kept at 77 K, and the surface temperature (T w ) of copper remained constant at 300 K.The fins with a thickness of 0.4 mm, height of 0.4 mm, and spacing of 0.6 mm were designed.Some assumptions were considered for the vaporization mechanism of the LN2 droplet in simulation. 1.The effect of temperature on physical properties was disregarded.2. The interface temperature was assumed to be at saturation temperature.3. Heat generation induced by friction was not considered.4. The gas-liquid pair velocity was small, so it could be regarded as incompressible Some assumptions were considered for the vaporization mechanism of the LN 2 droplet in simulation.
The effect of temperature on physical properties was disregarded.2.
The interface temperature was assumed to be at saturation temperature.

3.
Heat generation induced by friction was not considered.4.
The gas-liquid pair velocity was small, so it could be regarded as incompressible flow.
The analysis was carried out using the dimensionless parameters, namely diameter (d*), velocity (v*), density (ρ*) and time (t*) are defined as follows: The Reynolds number and Weber number of the LN 2 droplet impacting the wall surface are as follows: The residual radius of droplet vaporization can be calculated from the residual liquid phase f res : where, f res represents the volume fraction of residual liquid.

Governing Equations
In the direct numerical simulation, the governing equations are solved directly by numerical calculation which avoids the limit of the selection of empirical parameters and the accuracy of the turbulence model.The droplet vaporization, a two-phase flow process, is described by a two-dimensional incompressible N-S equation.The momentum equation is as follows: where, u is a velocity vector, ρ is the fluid density, p is a pressure, µ is the fluid viscosity, and D is the deformation rate tensor.
The last term in Equation ( 8) determines the surface tension term through surface tension coefficient (σ), surface curvature (κ), and Dirac delta function (δ s ) so that the surface tension only acts on the phase interface.
The continuity equation yields: The gas-liquid continuity equation was established by using the homogeneous flow model and defined as: where, f represents the volume fraction of liquid, then the physical properties in the flow field can be obtained according to f : The subscripts l and v represent liquid and gas phase, respectively.The energy equation yields: where, T, c p , and λ are termed as the temperature field, heat capacity of the fluid, and thermal conductivity, respectively.In this paper, the numerical simulation software basilisk is used to solve the incompressible N-S equation by employing the geometric volume of fluid (VOF) convective method to capture the phase interface.The detailed numerical method can be seen in Popinet [24].

Vaporization Model
The mass transfer rate is calculated through the heat balance at the interface by using the following relation.
In order to simulate the phase transition process, the surface vaporization flow (Stefan flow) caused by phase change was neglected, and the change in the volume of respective phases during the phase transition process was considered.The specific method to calculate the interface displacement velocity u I according to the mass transfer rate before employing VOF convection is to obtain the superposition velocity field with the original velocity field and restore the velocity field after completing VOF convection.
The temperature field was solved accurately before calculating the interface displacement velocity.The convection term of the energy equation is handled in the same way as that of the VOF convection equation.Since the phase transition occurs at the gas-liquid interface at the saturation temperature, the Dirichlet boundary condition at the interface is realized by adding a correction source term at the interface to solve the temperature diffusion equation.It is ensured that the interface temperature remains unchanged at the saturation temperature: where, τ c is the control time, and δ s is the Dirac function δ centering the interface.Furthermore, the temperature gradient at the interface is approximately estimated by calculating the temperature gradient on the surface of grid of the vapor and liquid phase adjacent to the interface (taking the two-dimensional case as an example), as shown in Figure 2: where, n x and n y are the component of the interface normal vector along the x-and ydirection, respectively, and ∆ is the length of the grid side.The source code for this method is visible in Magdelaine's Sandbox [25].
Energies 2022, 15, 8700 5 of 20 where,  and  are the component of the interface normal vector along the x-and ydirection, respectively, and Δ is the length of the grid side.The source code for this method is visible in Magdelaine's Sandbox [25].In the present model, the Peclet number, a ratio of the product of specific heat capacity and superheat to the latent heat of vaporization, is introduced to characterize the weight of convection rate and diffusion rate.

Boundary Condition
Initially, the gas is under the static state so the left and right sides of the calculation domain are at the symmetric boundary conditions.The wall surface is defined as the lower boundary at a constant temperature, and the upper boundary is a free flow condition.The domain was considered as 10 times the initial diameter of the droplet.The physical properties of the employed fluids are shown in Table 1, and the physical properties of the solid copper surface are shown in Table 2.The thermophysical properties of the fluids were derived from the software Refprop, which is a physical property database of National Institute of Standards and Technology, NIST.The thermophysical properties of copper are derived from the Thermophysical Properties of Matter Databases, TPMD.which is a manual published by the Center for Information and Numerical Data Analysis and Synthesis, CINDAS.In the present model, the Peclet number, a ratio of the product of specific heat capacity and superheat to the latent heat of vaporization, is introduced to characterize the weight of convection rate and diffusion rate.

Boundary Condition
Initially, the gas is under the static state so the left and right sides of the calculation domain are at the symmetric boundary conditions.The wall surface is defined as the lower boundary at a constant temperature, and the upper boundary is a free flow condition.The domain was considered as 10 times the initial diameter of the droplet.The physical properties of the employed fluids are shown in Table 1, and the physical properties of the solid copper surface are shown in Table 2.The thermophysical properties of the fluids were derived from the software Refprop, which is a physical property database of National Institute of Standards and Technology, NIST.The thermophysical properties of copper are derived from the Thermophysical Properties of Matter Databases, TPMD.which is a manual published by the Center for Information and Numerical Data Analysis and Synthesis, CINDAS.The open source code software basilisk is used to directly solve the above-mentioned N-S equation.The variable impact of the characteristic spatial scale of the droplet was considered, i.e., some regions remained stable with time and showed no obvious features, while other regions contained details of the impact such as the interaction of the droplet with the wall.In order to reduce the number of mesh, the adaptive mesh refinement (AMR) method was used to reduce computing resource consumption during DNS process.The software basilisk uses quadtree/octree to achieve effective AMR to accommodate details.For each time step, when the discretization error is greater than the set value, the grid is automatically encrypted.The multistage Poisson equation was used to solve the pressure, and the complex solid boundary is represented by the VOF method.

Grid Independence Verification
The process of droplets hitting the heated wall involves vaporization and spreading, and the grid independence of the phase transition model and motion model are verified respectively.
Three kinds of adaptive mesh were used to verify the grid independence.The levels of the mesh were coarse, fine, and refined, which means that the grid encryption levels are level 7, level 8, and level 9, respectively.The corresponding grid sizes are 2 × 10 −9 m, 2 × 10 −10 m, 2 × 10 −11 m, respectively.A schematic of the fine mesh is shown in Figure 3.During the calculation, the software Basilisk will also perform local encryption on the mesh, and the local mesh encryption level is 10, the corresponding mesh size is 2 × 10 −11 m.

Grid Independence Verification of Motion Model
The process of the oil droplet impinging on the copper surface at room temperature was simulated, The grid of different encryption grades is shown in Figure 5.When the encryption grade is coarse, the edge of the droplet is blurred.When the encryption grade is fine, the edge of the droplet is clearer.When the encryption grade is refined, the edge and inside of the droplet are well encrypted, which can accurately show the shape and contour of the droplet.Therefore, in order to ensure the accuracy of numerical simulation and save computing resources, the mesh encryption grade was considered as the fine one.

Grid Independence Verification of Phase Transition Model
The process of the LN2 droplet at 77 K impinging on the copper surface at 300 K was simulated.The remaining radius of the vaporized droplet at the different grid encryption levels is shown in Figure 4.All three kinds of mesh possess good convergence.When the encryption grade is coarse, the simulation results of the residual radius of the droplet are not accurate enough, while the results are basically consistent for the fine or refined mesh.

Grid Independence Verification of Motion Model
The process of the oil droplet impinging on the copper surface at room temperature was simulated, The grid of different encryption grades is shown in Figure 5.When the encryption grade is coarse, the edge of the droplet is blurred.When the encryption grade

Grid Independence Verification of Motion Model
The process of the oil droplet impinging on the copper surface at room temperature was simulated, The grid of different encryption grades is shown in Figure 5.When the encryption grade is coarse, the edge of the droplet is blurred.When the encryption grade is fine, the edge of the droplet is clearer.When the encryption grade is refined, the edge and inside of the droplet are well encrypted, which can accurately show the shape and contour of the droplet.Therefore, in order to ensure the accuracy of numerical simulation and save computing resources, the mesh encryption grade was considered as the fine one.

Verification of Phase Transition Model
Experimental data obtained by Volkov et al. [26] were used to verify the numerical model.The experiment was conducted by taking into account the vaporization of a suspended water droplet with an initial diameter of about 3 mm in a hollow transparent silicon glass cylinder with a radius of 0.1 m.A hot air fan was located below the cylinder that blew wind in an upward direction at 3~4.5 m/s.The initial diameter of the water drop was 3.06 mm at 30℃ temperature and the hot air at 400 ℃ and 800 ℃, respectively, blew at speed of 3 m/s is considered for numerical simulation verification.
The compared results are shown In Figure 6.It can be seen that the vaporization rate obtained by simulation is faster than that of the experiment, resulting in a sharp decrease in diameter of the droplet.However, the variation pattern of the diameter obtained from simulation is in good agreement with experiment results, i.e., an average error of 13.08% at 400 ℃ and 1.49% at 800 ℃.The ISM algorithm used in this paper does not consider the influence of disturbance speed on the movement speed of the phase interface.Compared with the working condition of 400 °C, the temperature gradient at 800 °C is larger, and the influence of the disturbance speed is small, so the numerical calculation results are more consistent with the test results.

Comparison with Experimental Results Verification of Phase Transition Model
Experimental data obtained by Volkov et al. [26] were used to verify the numerical model.The experiment was conducted by taking into account the vaporization of a suspended water droplet with an initial diameter of about 3 mm in a hollow transparent silicon glass cylinder with a radius of 0.1 m.A hot air fan was located below the cylinder that blew wind in an upward direction at 3~4.5 m/s.The initial diameter of the water drop was 3.06 mm at 30°C temperature and the hot air at 400 °C and 800 °C, respectively, blew at speed of 3 m/s is considered for numerical simulation verification.
The compared results are shown In Figure 6.It can be seen that the vaporization rate obtained by simulation is faster than that of the experiment, resulting in a sharp decrease in diameter of the droplet.However, the variation pattern of the diameter obtained from simulation is in good agreement with experiment results, i.e., an average error of 13.08% at 400 °C and 1.49% at 800 °C.The ISM algorithm used in this paper does not consider the influence of disturbance speed on the movement speed of the phase interface.Compared with the working condition of 400 • C, the temperature gradient at 800 • C is larger, and the influence of the disturbance speed is small, so the numerical calculation results are more consistent with the test results.

Verification of Motion Model
Figure 7 shows the morphological comparison of the simulation and experiment [27] when the water droplets hit the wall.In the environment of 0.1 mPa, the water droplets with an initial diameter of 2.6 mm hit a hydrophobic wall with a contact angle of 160 • vertically at an initial velocity of 0.5 m/s at a distance of 3 mm from the wall.The simulation results show that the water droplets have successively experienced spreading, retraction, rebound, and peristaltic rise after hitting the wall, which is consistent with the experimental results.
in diameter of the droplet.However, the variation pattern of the diameter obtained from simulation is in good agreement with experiment results, i.e., an average error of 13.08% at 400 ℃ and 1.49% at 800 ℃.The ISM algorithm used in this paper does not consider the influence of disturbance speed on the movement speed of the phase interface.Compared with the working condition of 400 °C, the temperature gradient at 800 °C is larger, and the influence of the disturbance speed is small, so the numerical calculation results are more consistent with the test results.

Verification of Motion Model
Figure 7 shows the morphological comparison of the simulation and experiment [27] when the water droplets hit the wall.In the environment of 0.1 mPa, the water droplets with an initial diameter of 2.6 mm hit a hydrophobic wall with a contact angle of 160° vertically at an initial velocity of 0.5 m/s at a distance of 3 mm from the wall.The simulation results show that the water droplets have successively experienced spreading, retraction, rebound, and peristaltic rise after hitting the wall, which is consistent with the experimental results.Figure 8 shows the spreading coefficient trend comparison of the simulation and experiment [28] when the water droplets hit the wall, where the abscissa t* represents dimensionless time, and the ordinate d/d0 is the ratio of the spreading diameter of the liquid film to the initial diameter of the droplet.In the experiment, a single droplet with a diameter of 0.5 mm hits the surface of the iron-nickel alloy at an initial velocity of 2.32 m/s, and the droplet diffusion diameter with a static contact angle of 57° is measured.The trend of the two curves is highly consistent, and the maximum deviation does not exceed 10%, which proves the reliability of the model.

Comparison with Analytical Results
The vaporization of an n-decane droplet in stationary air was simulated, and the nu- Figure 8 shows the spreading coefficient trend comparison of the simulation and experiment [28] when the water droplets hit the wall, where the abscissa t* represents dimensionless time, and the ordinate d/d0 is the ratio of the spreading diameter of the liquid film to the initial diameter of the droplet.In the experiment, a single droplet with a diameter of 0.5 mm hits the surface of the iron-nickel alloy at an initial velocity of 2.32 m/s, and the droplet diffusion diameter with a static contact angle of 57 • is measured.The trend of the two curves is highly consistent, and the maximum deviation does not exceed 10%, which proves the reliability of the model.

Verification of Motion Model
Figure 7 shows the morphological comparison of the simulation and experiment [27] when the water droplets hit the wall.In the environment of 0.1 mPa, the water droplets with an initial diameter of 2.6 mm hit a hydrophobic wall with a contact angle of 160° vertically at an initial velocity of 0.5 m/s at a distance of 3 mm from the wall.The simulation results show that the water droplets have successively experienced spreading, retraction, rebound, and peristaltic rise after hitting the wall, which is consistent with the experimental results.Figure 8 shows the spreading coefficient trend comparison of the simulation and experiment [28] when the water droplets hit the wall, where the abscissa t* represents dimensionless time, and the ordinate d/d0 is the ratio of the spreading diameter of the liquid film to the initial diameter of the droplet.In the experiment, a single droplet with a diameter of 0.5 mm hits the surface of the iron-nickel alloy at an initial velocity of 2.32 m/s, and the droplet diffusion diameter with a static contact angle of 57° is measured.The trend of the two curves is highly consistent, and the maximum deviation does not exceed 10%, which proves the reliability of the model.

Comparison with Analytical Results
The vaporization of an n-decane droplet in stationary air was simulated, and the numerical model was validated according to the reported  law [29].By considering the assumption of uniformity and quasi-steady state, the analytic solution for the behavior of

Comparison with Analytical Results
The vaporization of an n-decane droplet in stationary air was simulated, and the numerical model was validated according to the reported d 2 law [29].By considering the assumption of uniformity and quasi-steady state, the analytic solution for the behavior of droplets can be obtained by solving the conservation equations of mass and energy.The variation in the diameter can be expressed by d 2 law.
where K is the constant determined by the physical properties of gas and liquid phase.Figure 9 compares the numerical and analytical solutions for the d 2 value of the vaporized droplet.The numerical simulation calculates the diameter of the droplet based on liquid volume.It can be seen that the numerical results are quite consistent with the analytical results, and the vaporization rate is slightly higher than that of the quasi-steady state.At the end of the simulation, the d 2 value of droplet decreases from 3.98 to 1.4 mm 2 .
Energies 2022, 15, 8700 10 of 22 droplets can be obtained by solving the conservation equations of mass and energy.The variation in the diameter can be expressed by  law.
where K is the constant determined by the physical properties of gas and liquid phase.Figure 9 compares the numerical and analytical solutions for the  value of the vaporized droplet.The numerical simulation calculates the diameter of the droplet based on liquid volume.It can be seen that the numerical results are quite consistent with the analytical results, and the vaporization rate is slightly higher than that of the quasi-steady state.At the end of the simulation, the  value of droplet decreases from 3.98 to 1.4 mm 2 .

Effect of the Micro-Structured Surface
After impinging on the wall, the droplet's behavior is mainly affected by inertia force, surface tension, viscous force, inner wall pressure, and boundary layer resistance.The liquid viscosity and surface tension can be the critical parameter to impact the collision behavior [29].At an appropriate initial velocity, the droplets can impinge the wall and experience the reciprocating movement of the spreading and contraction under the combination of the forces mentioned above.Under the consumption of viscous force and boundary layer resistance, the amplitude of spreading radius as well as velocity is decreased gradually and can reach the equilibrium by the combined effect of the internal pressure and surface intension.
Firstly, the process of droplet impinging on the smooth surface was simulated.The initial impact velocity of LN2 was 1 m/s.The low surface tension and high impact velocity resulted in the fine deformation of the droplet during falling down.After hitting the bottom wall, a strong deformation occurred at the bottom surface of the droplet, while the upper part of the droplet remained in perfect shape.The liquid splashed at both ends of the droplet, and droplets started to flatten, break up, and spread out over a film that is much thinner than the initial radius of the droplet.The simulated Weber number and Reynolds number were 180 and 9914, respectively.According to the morphological determination formula given by Sommerfeld [30], i.e.,  = √, the calculated value of K was 133.87 that is consistent with the mode of crushing and splashing.After impinging on the wall, the droplet's behavior is mainly affected by inertia force, surface tension, viscous force, inner wall pressure, and boundary layer resistance.The liquid viscosity and surface tension can be the critical parameter to impact the collision behavior [29].At an appropriate initial velocity, the droplets can impinge the wall and experience the reciprocating movement of the spreading and contraction under the combination of the forces mentioned above.Under the consumption of viscous force and boundary layer resistance, the amplitude of spreading radius as well as velocity is decreased gradually and can reach the equilibrium by the combined effect of the internal pressure and surface intension.
Firstly, the process of droplet impinging on the smooth surface was simulated.The initial impact velocity of LN 2 was 1 m/s.The low surface tension and high impact velocity resulted in the fine deformation of the droplet during falling down.After hitting the bottom wall, a strong deformation occurred at the bottom surface of the droplet, while the upper part of the droplet remained in perfect shape.The liquid splashed at both ends of the droplet, and droplets started to flatten, break up, and spread out over a film that is much thinner than the initial radius of the droplet.The simulated Weber number and Reynolds number were 180 and 9914, respectively.According to the morphological determination formula given by Sommerfeld [30], i.e., K = We √ Re, the calculated value of K was 133.87 that is consistent with the mode of crushing and splashing.
The spreading morphology of a droplet with a dimensionless diameter of 1 impacted the micro-structured surface with the grooves' spacing of 0.2 mm, the height of 0.4 mm, and a thickness of 0.4 mm, as shown in Figure 10a.The small spacing between fins caused the gas to become trapped between the adjacent fins when the droplet hits the wall, and trapped gas pushed the droplet upward at the micro-structured surface.The droplets cannot fully spread out smoothly between the adjacent structures, conforming to the Cassie-Baxter model, as shown in Figure 10b.The solid area and liquid phase showed a very small contact area and failed to achieve the large contact area.Figure 11 shows the comparison of the spreading morphology between LN 2 droplets impinging on the smooth surface and the micro-structured surface, respectively.At the dimensionless time, t* = 5, both ends of the droplet retracted and curled up on the micro-structured surface, impeding the continuous spreading of the droplet.However, at t* = 9, the droplet dispersed into several large liquid films at the micro-structured surface, while the liquid films on the smooth surface remained nearly uniform and intact.The diameter of the liquid film on the smooth and the micro-structured surface varies with time, as shown in Figure 12.Before t* = 5, the dimensionless diameter of the liquid film on smooth surface and the micro-structured surface is almost the same, and the spreading diameter on the smooth surface is slightly larger than that on the micro-structured surface.Beyond t* = 5, the spreading speed of the droplet on the micro-structured surface is significantly reduced due to the contraction and breakup of the droplets into multiple liquid films, as shown in Figure 11.The spreading morphology of a droplet with a dimensionless diameter of 1 impacted the micro-structured surface with the grooves' spacing of 0.2 mm, the height of 0.4 mm, and a thickness of 0.4 mm, as shown in Figure 10a.The small spacing between fins caused the gas to become trapped between the adjacent fins when the droplet hits the wall, and trapped gas pushed the droplet upward at the micro-structured surface.The droplets cannot fully spread out smoothly between the adjacent structures, conforming to the Cassie-Baxter model, as shown in Figure 10b.The solid area and liquid phase showed a very small contact area and failed to achieve the large contact area.Figure 11 shows the comparison of the spreading morphology between LN2 droplets impinging on the smooth surface and the micro-structured surface, respectively.At the dimensionless time, t* = 5, both ends of the droplet retracted and curled up on the micro-structured surface, impeding the continuous spreading of the droplet.However, at t* = 9, the droplet dispersed into several large liquid films at the micro-structured surface, while the liquid films on the smooth surface remained nearly uniform and intact.The diameter of the liquid film on the smooth and the micro-structured surface varies with time, as shown in Figure 12.Before t* = 5, the dimensionless diameter of the liquid film on smooth surface and the microstructured surface is almost the same, and the spreading diameter on the smooth surface is slightly larger than that on the micro-structured surface.Beyond t* = 5, the spreading speed of the droplet on the micro-structured surface is significantly reduced due to the contraction and breakup of the droplets into multiple liquid films, as shown in Figure 11.The spreading morphology of a droplet with a dimensionless diameter of 1 impacted the micro-structured surface with the grooves' spacing of 0.2 mm, the height of 0.4 mm, and a thickness of 0.4 mm, as shown in Figure 10a.The small spacing between fins caused the gas to become trapped between the adjacent fins when the droplet hits the wall, and trapped gas pushed the droplet upward at the micro-structured surface.The droplets cannot fully spread out smoothly between the adjacent structures, conforming to the Cassie-Baxter model, as shown in Figure 10b.The solid area and liquid phase showed a very small contact area and failed to achieve the large contact area.Figure 11 shows the comparison of the spreading morphology between LN2 droplets impinging on the smooth surface and the micro-structured surface, respectively.At the dimensionless time, t* = 5, both ends of the droplet retracted and curled up on the micro-structured surface, impeding the continuous spreading of the droplet.However, at t* = 9, the droplet dispersed into several large liquid films at the micro-structured surface, while the liquid films on the smooth surface remained nearly uniform and intact.The diameter of the liquid film on the smooth and the micro-structured surface varies with time, as shown in Figure 12.Before t* = 5, the dimensionless diameter of the liquid film on smooth surface and the microstructured surface is almost the same, and the spreading diameter on the smooth surface is slightly larger than that on the micro-structured surface.Beyond t* = 5, the spreading speed of the droplet on the micro-structured surface is significantly reduced due to the contraction and breakup of the droplets into multiple liquid films, as shown in Figure 11.When the droplet impinges the liquid film downward, the wall will exert a reaction force on the droplet, forming a certain dynamic pressure inside the droplet.The raised microstructure surface brings greater diffusion resistance, increases kinetic energy loss, and when the droplet impinges the microstructure surface, the diffusion speed of the liquid film is affected.From the analysis of the micro-structured surface, the maximum diameter and spreading speed of the LN2 droplets is smaller than that of a smooth surface.In addition, the existence of air gaps between fins causes the LN2 droplet to wet the surface incompletely.

Effect of the Thermophysical Properties
The viscosity and surface tension of the fluid have a great influence on the movement of the liquid film, and the particularity of the cryogenic fluid causes its spread to be different from that of the normal temperature fluid.The LN2 droplet and oil droplet impinging the smooth wall under the same initial conditions are simulated.In the numerical simulation, the initial diameter is 2 mm, the initial speed is 1 m/s, and the drop height is 1 mm.Figures 13 and 14 show the spreading process of oil droplets and liquid nitrogen droplets impinging the smooth wall, respectively.When the droplet impinges the liquid film downward, the wall will exert a reaction force on the droplet, forming a certain dynamic pressure inside the droplet.The raised microstructure surface brings greater diffusion resistance, increases kinetic energy loss, and when the droplet impinges the microstructure surface, the diffusion speed of the liquid film is affected.From the analysis of the micro-structured surface, the maximum diameter and spreading speed of the LN 2 droplets is smaller than that of a smooth surface.In addition, the existence of air gaps between fins causes the LN 2 droplet to wet the surface incompletely.

Effect of the Thermophysical Properties
The viscosity and surface tension of the fluid have a great influence on the movement of the liquid film, and the particularity of the cryogenic fluid causes its spread to be different from that of the normal temperature fluid.The LN2 droplet and oil droplet impinging the smooth wall under the same initial conditions are simulated.In the numerical simulation, the initial diameter is 2 mm, the initial speed is 1 m/s, and the drop height is 1 mm.Figures 13 and 14 show the spreading process of oil droplets and liquid nitrogen droplets impinging the smooth wall, respectively.Different from the LN2 droplet, the oil droplet maintains a good shape during the impact of the wall.At the moment t* = 1.8, the oil droplet reaches its maximum spreading diameter, after which its edge warps and begins to contract towards the center.Until t* = 4.2, it returns to its initial shape and then rebounds off the wall.Under the same initial conditions, the surface tension and viscosity of the LN2 droplet are much smaller than that of the oil droplet.After the initial impact on the wall, the LN2 droplet continues to spread.At the moment of t* = 1.8, the LN2 droplet shows a trend of high center and low on both sides.The spreading speed is significantly higher than that of the oil droplet, and the rebound phenomenon does not occur.
Figure 15 is the comparison of the spreading coefficient of the nitrogen droplet and the oil droplet.For the oil droplet, it first spread and then retracted after impinging the wall, and the spreading diameter shows a trend of first increasing and then decreasing, reaching the maximum diffusion value at the time of t* = 1.8.For the LN2 droplet, it continues to spread after hitting the wall, the diameter of the liquid film continues to expand, the liquid film becomes thinner and thinner, and the diffusion rate gradually slows down.Viscosity and surface tension together determine the morphological change process of the droplet after impinging the wall, of which, the LN2 droplets are much smaller than the oil Different from the LN2 droplet, the oil droplet maintains a good shape during the impact of the wall.At the moment t* = 1.8, the oil droplet reaches its maximum spreading diameter, after which its edge warps and begins to contract towards the center.Until t* = 4.2, it returns to its initial shape and then rebounds off the wall.Under the same initial conditions, the surface tension and viscosity of the LN2 droplet are much smaller than that of the oil droplet.After the initial impact on the wall, the LN2 droplet continues to spread.At the moment of t* = 1.8, the LN2 droplet shows a trend of high center and low on both sides.The spreading speed is significantly higher than that of the oil droplet, and the rebound phenomenon does not occur.
Figure 15 is the comparison of the spreading coefficient of the nitrogen droplet and the oil droplet.For the oil droplet, it first spread and then retracted after impinging the wall, and the spreading diameter shows a trend of first increasing and then decreasing, reaching the maximum diffusion value at the time of t* = 1.8.For the LN2 droplet, it continues to spread after hitting the wall, the diameter of the liquid film continues to expand, the liquid film becomes thinner and thinner, and the diffusion rate gradually slows down.Viscosity and surface tension together determine the morphological change process of the droplet after impinging the wall, of which, the LN2 droplets are much smaller than the oil droplets.Under the same initial conditions, the surface energy of the oil droplet is sufficient to overcome its kinetic energy, and the wall has a greater resistance to it.Therefore, after the impact of the oil droplet, the liquid film first spreads, the spreading speed is gradually smaller, and then begins to contract under the action of surface tension.In contrast, after the impact of the LN2, no contraction occurs, the kinetic energy has been dominant, and the liquid film continues to spread, but the spreading speed decreases under the action of the surface tension.Under the same initial conditions, the surface energy of the oil droplet is sufficient to overcome its kinetic energy, and the wall has a greater resistance to it.Therefore, after the impact of the oil droplet, the liquid film first spreads, the spreading speed is gradually smaller, and then begins to contract under the action of surface tension.In contrast, after the impact of the LN2, no contraction occurs, the kinetic energy has been dominant, and the liquid film continues to spread, but the spreading speed decreases under the action of the surface tension.The vaporization of the LN 2 droplet with an initial diameter of 2 mm and a velocity of 1 m/s on the smooth and micro-structured hot wall, respectively, was simulated by considering influences of gravity during falling and impinging.Initially, the dimensionless diameter and velocity are kept at 1, and the initial residual r* is 0.5.The vaporization process of the LN 2 droplet is studied within t*=10.Figure 16 shows the morphology of the droplet during the vaporization on smooth surface and micro-structured surface at t* = 2.The temperature of the solid wall surface is much higher than the boiling temperature of LN 2 that causes an obvious Leidenfrost phenomenon.The vaporization of the LN2 droplet with an initial diameter of 2 mm and a velocity of 1 m/s on the smooth and micro-structured hot wall, respectively, was simulated by considering influences of gravity during falling and impinging.Initially, the dimensionless diameter and velocity are kept at 1, and the initial residual r* is 0.5.The vaporization process of the LN2 droplet is studied within t*=10.Figure 16 shows the morphology of the droplet during the vaporization on smooth surface and micro-structured surface at t* = 2.The temperature of the solid wall surface is much higher than the boiling temperature of LN2 that causes an obvious Leidenfrost phenomenon.The presence of the fins and trapped gas inhibits the droplet to get in contact with the surface, resulted in the detached droplet that doesn't encroach into grooves.
As the residual radius decrease with time being (Figure 17), it can be seen that the droplet vaporizes more slowly when adding fins to the wall.The pressure field marked by contour lines at t* = 2.5 is shown in Figure 18.At the micro-structured surface, there is The presence of the fins and trapped gas inhibits the droplet to get in contact with the surface, resulted in the detached droplet that doesn't encroach into grooves.
As the residual radius decrease with time being (Figure 17), it can be seen that the droplet vaporizes more slowly when adding fins to the wall.The pressure field marked by contour lines at t* = 2.5 is shown in Figure 18.At the micro-structured surface, there is a negative pressure area under the droplet, which makes the droplet lifetime longer during vaporization.At t* = 5, the temperature of the smooth surface is significantly higher than that of the micro-structured surface, and the temperature curve shows violent fluctuation as shown in Figure 19.The gas produced during droplet vaporization is responsible for fluctuation due to a lower thermal conductivity than that of LN2.The temperature of the area covered by gas shows a high level while the temperature of the wall that is exposed to LN2 shows a significant drop.The overall temperature curve shows a trend of high on both sides and low in the middle, confirming the fact that cooling capacity is transferred to both sides from the area where the droplet falls and hits.At t* = 5, the temperature of the smooth surface is significantly higher than that of the micro-structured surface, and the temperature curve shows violent fluctuation as shown in Figure 19.The gas produced during droplet vaporization is responsible for fluctuation due to a lower thermal conductivity than that of LN 2 .The temperature of the area covered by gas shows a high level while the temperature of the wall that is exposed to LN 2 shows a significant drop.The overall temperature curve shows a trend of high on both sides and low in the middle, confirming the fact that cooling capacity is transferred to both sides from the area where the droplet falls and hits.Specifically, the droplet boils upon the film adjacent to the wall at high temperature, and a clear vapor layer can be seen between the droplet and hot wall.In the case of the micro-structured surface, when the droplet is about to impact the wall surface, the entrapped gas further hinders the droplet to come into contact with the wall surface.Thus, a negative pressure is formed between the fins, and the vaporization rate is lower than that of the smooth surface.

Effect of Weber Number
The Weber number, a decisive parameter of droplet breakage, is determined by the ratio of inertial force to surface tension.Therefore, the spreading of the droplet after impinging on the wall was simulated with various Weber numbers for the smooth wall and micro-structured wall.The simulated Weber numbers were 30, 180, 360, and 450.
The spreading morphology of droplet on the smooth surface at Weber number of 30 is shown in Figure 20.At the Weber number of 30, the spreading morphology of LN2 droplets on the smooth and micro-structured surface was quite different from that of other cases of different Weber numbers.At the smooth surface, the droplet did not splash at the Weber number of 30, and the droplet's motor pattern changed from the vertical motion to horizontal motion due to the contact with the solid wall.The symmetrical shape of liquid crowns appeared at both ends of the droplet with inward contraction, and the droplet spread out from the middle to form a thin film.The thickness of the film decreased continuously.After reaching the maximum spreading coefficient (d/d0) of 4.39, the droplet split from the middle and contracted with the shape of liquid crowns on both sides, and Specifically, the droplet boils upon the film adjacent to the wall at high temperature, and a clear vapor layer can be seen between the droplet and hot wall.In the case of the micro-structured surface, when the droplet is about to impact the wall surface, the entrapped gas further hinders the droplet to come into contact with the wall surface.Thus, a negative pressure is formed between the fins, and the vaporization rate is lower than that of the smooth surface.

Effect of Weber Number
The Weber number, a decisive parameter of droplet breakage, is determined by the ratio of inertial force to surface tension.Therefore, the spreading of the droplet after impinging on the wall was simulated with various Weber numbers for the smooth wall and micro-structured wall.The simulated Weber numbers were 30, 180, 360, and 450.
The spreading morphology of droplet on the smooth surface at Weber number of 30 is shown in Figure 20.At the Weber number of 30, the spreading morphology of LN 2 droplets on the smooth and micro-structured surface was quite different from that of other cases of different Weber numbers.At the smooth surface, the droplet did not splash at the Weber number of 30, and the droplet's motor pattern changed from the vertical motion to horizontal motion due to the contact with the solid wall.The symmetrical shape of liquid crowns appeared at both ends of the droplet with inward contraction, and the droplet spread out from the middle to form a thin film.The thickness of the film decreased continuously.After reaching the maximum spreading coefficient (d/d 0 ) of 4.39, the droplet split from the middle and contracted with the shape of liquid crowns on both sides, and the dimensionless height (h/d 0 ) of 1.55 was obtained for the liquid crowns after the contraction.As the Weber number became greater than 180, the dominant effect of the surface tension weakened, the droplet struck the surface of the wall, broke, and splashed, and the droplet spread out to a maximum radius in comparison to the initial value, obtaining a larger contact area.At the different Weber numbers, the width of the liquid film decreased on the microstructured surface, because a few droplets entered the grooves, which reduces certain droplet volume fractions that can spread out on the surface.In addition, the available fins promote the breaking up of the droplet.At t* = 6, a large number of tiny droplets break away from the initial droplet and get spattered.The existed fins can change the contact angle of the droplet.In addition, simulation results show that spreading on the smooth surface provides a relatively stable liquid film surface while spreading on the micro-structured surface causes a certain degree of fluctuation in the liquid film surface.This may be attributed to the stronger pressure wave inside the transmission of the liquid.The collision of the droplet with the solid surface formed a pressure wave that interacts with the surface tension, resulting in fluctuations, and eventually affecting the flow and motion form of the droplet.
Furthermore, the vaporization of the droplet impacting the smooth surface and micro-structured surface under the different Weber number was simulated, and the residual radius of the droplet was quantified.
The effect of Weber number on droplet vaporization on the smooth and micro-structured surface is shown in Figures 21 and 22, respectively.In the case of smooth surface, the variation in Weber number showed no obvious effect on the vaporization process of the droplet.However, at the micro-structured surface, the vaporization rate is the slowest at 30, and the highest vaporization rate is observed at the Weber number of 180.However, beyond the 180, i.e., Weber number of 360 and 450, the vaporization rate showed no significant increment.The smooth surface offered a better vaporization effect than that of the micro-structured surface.At the different Weber numbers, the width of the liquid film decreased on the microstructured surface, because a few droplets entered the grooves, which reduces certain droplet volume fractions that can spread out on the surface.In addition, the available fins promote the breaking up of the droplet.At t* = 6, a large number of tiny droplets break away from the initial droplet and get spattered.The existed fins can change the contact angle of the droplet.In addition, simulation results show that spreading on the smooth surface provides a relatively stable liquid film surface while spreading on the micro-structured surface causes a certain degree of fluctuation in the liquid film surface.This may be attributed to the stronger pressure wave inside the transmission of the liquid.The collision of the droplet with the solid surface formed a pressure wave that interacts with the surface tension, resulting in fluctuations, and eventually affecting the flow and motion form of the droplet.
Furthermore, the vaporization of the droplet impacting the smooth surface and microstructured surface under the different Weber number was simulated, and the residual radius of the droplet was quantified.
The effect of Weber number on droplet vaporization on the smooth and microstructured surface is shown in Figures 21 and 22, respectively.In the case of smooth surface, the variation in Weber number showed no obvious effect on the vaporization process of the droplet.However, at the micro-structured surface, the vaporization rate is the slowest at 30, and the highest vaporization rate is observed at the Weber number of 180.However, beyond the 180, i.e., Weber number of 360 and 450, the vaporization rate showed no significant increment.The smooth surface offered a better vaporization effect than that of the micro-structured surface.

Effects of Impinging Velocity
The spreading and vaporization process of the droplet impinging on both smooth and the micro-structured surfaces were simulated at a falling velocity of 0.5 m/s, 1 m/s, 2 m/s, and 3 m/s, respectively.A higher velocity upon a smooth surface provides more residual energy of the initial kinetic energy under certain surface tension and viscous force as well as higher the spreading velocity resulted in a short time required to reach the maximum spreading diameter.At the Weber number of 180, the surface tension is not enough to maintain the complete shape of the droplet so that the droplet hit the wall, broke up, and spattered into some small secondary droplet.Figure 23a shows the variation curve of the residual radius of LN2 droplet vaporizing on the smooth surface.Since the initial velocity of the droplet is different, the beginning of the vaporization is also different so the time origin to the moment when the lower end of the droplet touches the wall is translated, and the translated variation curves of the residual radius are shown in Figure 23b.It can be seen that the initial vaporization rate of the droplet at the different velocities is consistent after the droplet contacts the wall surface.At t* > 3, the droplet with the initial velocity of 0.5 m/s maintains a higher vaporization rate and obtains the minimum residual radius.With the increasing impinging velocity, the vaporization rate gradually decreases, and the value of the residual radius is positively correlated with the value of the initial velocity.
On the other hand, for the micro-structured surface, the variation curve of the residual radius is shown in Figure 24b, which is translated at the time origin to the moment when the lower end of the droplet touches the wall (refer to Figure 24a).The effect of

Effects of Impinging Velocity
The spreading and vaporization process of the droplet impinging on both smooth and the micro-structured surfaces were simulated at a falling velocity of 0.5 m/s, 1 m/s, 2 m/s, and 3 m/s, respectively.A higher velocity upon a smooth surface provides more residual energy of the initial kinetic energy under certain surface tension and viscous force as well as higher the spreading velocity resulted in a short time required to reach the maximum spreading diameter.At the Weber number of 180, the surface tension is not enough to maintain the complete shape of the droplet so that the droplet hit the wall, broke up, and spattered into some small secondary droplet.Figure 23a shows the variation curve of the residual radius of LN2 droplet vaporizing on the smooth surface.Since the initial velocity of the droplet is different, the beginning of the vaporization is also different so the time origin to the moment when the lower end of the droplet touches the wall is translated, and the translated variation curves of the residual radius are shown in Figure 23b.It can be seen that the initial vaporization rate of the droplet at the different velocities is consistent after the droplet contacts the wall surface.At t* > 3, the droplet with the initial velocity of 0.5 m/s maintains a higher vaporization rate and obtains the minimum residual radius.With the increasing impinging velocity, the vaporization rate gradually decreases, and the value of the residual radius is positively correlated with the value of the initial velocity.
On the other hand, for the micro-structured surface, the variation curve of the residual radius is shown in Figure 24b, which is translated at the time origin to the moment when the lower end of the droplet touches the wall (refer to Figure 24a).The effect of

Effects of Impinging Velocity
The spreading and vaporization process of the droplet impinging on both smooth and the micro-structured surfaces were simulated at a falling velocity of 0.5 m/s, 1 m/s, 2 m/s, and 3 m/s, respectively.A higher velocity upon a smooth surface provides more residual energy of the initial kinetic energy under certain surface tension and viscous force as well as higher the spreading velocity resulted in a short time required to reach the maximum spreading diameter.At the Weber number of 180, the surface tension is not enough to maintain the complete shape of the droplet so that the droplet hit the wall, broke up, and spattered into some small secondary droplet.Figure 23a shows the variation curve of the residual radius of LN 2 droplet vaporizing on the smooth surface.Since the initial velocity of the droplet is different, the beginning of the vaporization is also different so the time origin to the moment when the lower end of the droplet touches the wall is translated, and the translated variation curves of the residual radius are shown in Figure 23b.It can be seen that the initial vaporization rate of the droplet at the different velocities is consistent after the droplet contacts the wall surface.At t* > 3, the droplet with the initial velocity of 0.5 m/s maintains a higher vaporization rate and obtains the minimum residual radius.With the increasing impinging velocity, the vaporization rate gradually decreases, and the value of the residual radius is positively correlated with the value of the initial velocity.became larger.At the falling velocity of 0.5 m/s, the droplet maintained a complete shape before touching the wall surface.After slowly contacting the wall surface, it spread out into a liquid film about twice the size of the initial diameter and sufficiently transferred heat to the micro-structured surface, resulting in a fast vaporization effect.With the increasing velocity, the droplet quickly broke up and spattered after impinging on the surface with fins.The small secondary droplets produced rolled onto the micro-structured surface and failed to transfer heat.At t* = 10, no vaporization was observed in the droplet.

Conclusions
In summary, a direct numerical simulation was carried out to investigate the heat transfer performance of the LN2 droplet impinging a micro-structured surface composed of fins that were designed with a spacing much smaller than the diameter of the droplet.For a comparable study, numerical simulations with a smooth surface were also conducted.The structural parameters of the micro-structured surface affect the heat transfer.On the other hand, for the micro-structured surface, the variation curve of the residual radius is shown in Figure 24b, which is translated at the time origin to the moment when the lower end of the droplet touches the wall (refer to Figure 24a).The effect of impinging velocity and the droplet vaporization on the micro-structured surface is similar to the above-mentioned smooth surface.At the velocity of 0.5 m/s, the droplet showed the highest vaporization rate after contacting the surface.With the increasing velocity, the droplet maintained a high vaporization rate for a shorter time, and the remaining radius became larger.At the falling velocity of 0.5 m/s, the droplet maintained a complete shape before touching the wall surface.After slowly contacting the wall surface, it spread out into a liquid film about twice the size of the initial diameter and sufficiently transferred heat to the micro-structured surface, resulting in a fast vaporization effect.With the increasing velocity, the droplet quickly broke up and spattered after impinging on the surface with fins.The small secondary droplets produced rolled onto the micro-structured surface and failed to transfer heat.At t* = 10, no vaporization was observed in the droplet.impinging velocity and the droplet vaporization on the micro-structured surface is similar to the above-mentioned smooth surface.At the velocity of 0.5 m/s, the droplet showed the highest vaporization rate after contacting the surface.With the increasing velocity, the droplet maintained a high vaporization rate for a shorter time, and the remaining radius became larger.At the falling velocity of 0.5 m/s, the droplet maintained a complete shape before touching the wall surface.After slowly contacting the wall surface, it spread out into a liquid film about twice the size of the initial diameter and sufficiently transferred heat to the micro-structured surface, resulting in a fast vaporization effect.With the increasing velocity, the droplet quickly broke up and spattered after impinging on the surface with fins.The small secondary droplets produced rolled onto the micro-structured surface and failed to transfer heat.At t* = 10, no vaporization was observed in the droplet.

Conclusions
In summary, a direct numerical simulation was carried out to investigate the heat transfer performance of the LN2 droplet impinging a micro-structured surface composed of fins that were designed with a spacing much smaller than the diameter of the droplet.For a comparable study, numerical simulations with a smooth surface were also conducted.The structural parameters of the micro-structured surface affect the heat transfer.

Conclusions
In summary, a direct numerical simulation was carried out to investigate the heat transfer performance of the LN 2 droplet impinging a micro-structured surface composed of fins that were designed with a spacing much smaller than the diameter of the droplet.For a comparable study, numerical simulations with a smooth surface were also conducted.The structural parameters of the micro-structured surface affect the heat transfer.Accordingly, the droplet's behavior by means of the motion and liquid film morphology under the various conditions were investigated.The major findings are as follows.
(1) The Leidenfrost phenomenon occurs on both surfaces, namely micro-structured and smooth surfaces owing to a large temperature difference.The produced vapor layer hindered the heat transfer between the hot surfaces and the LN 2 droplet.The droplet cannot completely wet the micro-structured surface, resulting in a small spreading diameter in the counterpart of the droplet at the smooth surface.The slow vaporization process is observed at the micro-structured surface rather than that of the smooth surface.(2) The Weber number affects the heat transfer of liquid nitrogen droplets on both smooth and micro-structured surfaces.On the smooth surface, the Weber number greater than 180 showed a minor effect on the droplet's spreading morphology, and no obvious effect was observed at the vaporization rate.On the other hand, the micro-structured surface exhibited the highest vaporization rate at the Weber number of 180.(3) Impinging velocity showed a noticeable impact on the heat transfer behavior.At the micro-structured surface, a large impinging velocity can lead to severe crushing and splashing of the droplet.The secondary droplet rolls over on the micro-structured surface, and the heat transfer performance becomes worse than that of the smooth surface.

Figure 1 .
Figure 1.Geometric model of the LN 2 droplet impinging on a copper surface.

Figure 2 .
Figure 2. The calculation of temperature gradient.(The red surface is to be the calculated surface, and the gradient values of the green and yellow surfaces are used for approximation).

Figure 2 .
Figure 2. The calculation of temperature gradient.(The red surface is to be the calculated surface, and the gradient values of the green and yellow surfaces are used for approximation).

Figure 4 .
Figure 4. Residual radius of droplet vaporization at different grid encryption levels.

Figure 4 .
Figure 4. Residual radius of droplet vaporization at different grid encryption levels.

Figure 4 .
Figure 4. Residual radius of droplet vaporization at different grid encryption levels.

Figure 5 .
Figure 5.Comparison of grids with different encryption grades.

Figure 6 .
Figure 6.Comparison between numerical simulation and experimental results of Phase Transition Model.Figure 6.Comparison between numerical simulation and experimental results of Phase Transition Model.

Figure 6 .
Figure 6.Comparison between numerical simulation and experimental results of Phase Transition Model.Figure 6.Comparison between numerical simulation and experimental results of Phase Transition Model.

Figure 7 .
Figure 7.Comparison between numerical simulation and experimental results of Motion Model.

Figure 8 .
Figure 8.The spreading coefficient trend comparison of the simulation and experiment.

Figure 7 .
Figure 7.Comparison between numerical simulation and experimental results of Motion Model.

Figure 7 .
Figure 7.Comparison between numerical simulation and experimental results of Motion Model.

Figure 8 .
Figure 8.The spreading coefficient trend comparison of the simulation and experiment.

Figure 8 .
Figure 8.The spreading coefficient trend comparison of the simulation and experiment.

Figure 9 .
Figure 9.Comparison between numerical and analytical results in terms of d2 value.

Figure 9 .
Figure 9.Comparison between numerical and analytical results in terms of d2 value.

Figure 11 .Figure 12 .
Figure 11.Comparison of liquid films formed by droplet impinging on smooth and micro-structured surface, respectively.

Figure 12 .
Figure 12.Comparison of spreading diameter between smooth surface and micro-structured surface.

Energies 2022, 15 , 8700 13 of 22 Figure 13 .
Figure 13.The spreading process of the oil droplet impinging the smooth wall.Figure 13.The spreading process of the oil droplet impinging the smooth wall.

Figure 13 .
Figure 13.The spreading process of the oil droplet impinging the smooth wall.Figure 13.The spreading process of the oil droplet impinging the smooth wall.

Figure 13 .
Figure 13.The spreading process of the oil droplet impinging the smooth wall.

Figure 14 .
Figure 14.The spreading process of the LN2 droplet impinging the smooth wall.

Figure 14 .
Figure 14.The spreading process of the LN2 droplet impinging the smooth wall.

Figure 15 .
Figure 15.Comparison of spreading coefficient between the oil droplet and the LN2 droplet.

4. 2 .
LN2 Droplet Vaporization Process 4.2.1.Effect of the Micro-Structured SurfaceThe vaporization of the LN2 droplet with an initial diameter of 2 mm and a velocity of 1 m/s on the smooth and micro-structured hot wall, respectively, was simulated by

Figure 15 .
Figure 15.Comparison of spreading coefficient between the oil droplet and the LN2 droplet.

Figure 15 .
Figure 15.Comparison of spreading coefficient between the oil droplet and the LN2 droplet.

Energies 2022 ,Figure 17 .
Figure 17.Variation in residual radius of vaporized droplet on the smooth and micro-structured surfaces, respectively.

Figure 17 .
Figure 17.Variation in residual radius of vaporized droplet on the smooth and micro-structured surfaces, respectively.

Figure 17 .
Figure 17.Variation in residual radius of vaporized droplet on the smooth and micro-structured surfaces, respectively.

Energies 2022, 15 , 8700 17 of 22 Figure 20 .
Figure 20.The spreading morphology of droplet on the smooth surface at Weber number of 30.

Figure 20 .
Figure 20.The spreading morphology of droplet on the smooth surface at Weber number of 30.
The spreading diameter of droplet at the different Weber numbers.(b) The residual radius at the different Weber numbers.

Figure 21 .
Figure 21.Effect of the Weber number on the vaporization of droplet on the smooth surface.
The spreading diameter of droplet at the different Weber numbers.(b) The residual radius at the different Weber numbers.

Figure 22 .
Figure 22.Effect of the Weber number on droplet vaporization on the micro-structured surface.

Figure 21 .
Figure 21.Effect of the Weber number on the vaporization of droplet on the smooth surface.

Figure 21 .
Figure 21.Effect of the Weber number on the vaporization of droplet on the smooth surface.
The spreading diameter of droplet at the different Weber numbers.(b) The residual radius at the different Weber numbers.

Figure 22 .
Figure 22.Effect of the Weber number on droplet vaporization on the micro-structured surface.

Figure 22 .
Figure 22.Effect of the Weber number on droplet vaporization on the micro-structured surface.
Original curve of residual radius.( b) Translated curve of residual radius.

Figure 23 .
Figure 23.Droplet vaporization on the smooth surface with the various impinging velocities.

Figure 24 .
Figure 24.Droplet vaporization on the micro-structured surface with various impinging velocities.

Figure 23 .
Figure 23.Droplet vaporization on the smooth surface with the various impinging velocities.
Original curve of residual radius.( b) Translated curve of residual radius.

Figure 23 .
Figure 23.Droplet vaporization on the smooth surface with the various impinging velocities.

Figure 24 .
Figure 24.Droplet vaporization on the micro-structured surface with various impinging velocities.

Figure 24 .
Figure 24.Droplet vaporization on the micro-structured surface with various impinging velocities.

Table 2 .
Physical properties of copper.