“…Physically, this continuum hydrodynamic model is closely related to the so-called model "H" which was originally devised to describe the critical dynamics of thermal fluctuations [41,42]. Supplemented with the hydrodynamic boundary conditions derived in our previous works [22,25], this model can be used to fully take into account the various physical processes involved in the contact line dynamics, including phase transition (evaporation or condensation) and capillary flow in the bulk fluid region [26,27], boundary slip of fluid [25,39,40,43], temperature slip across the fluid-solid interface [44,45], and mechanical-thermal cross coupling at the fluid-solid interface [46]. Due to the use of diffuse-interface method, the stress and thermal singularities are resolved automatically.…”