以脉宽为毫秒量级的单脉冲激光逆重力方向对金属的打孔实验为模拟对象，利用ANSYS软件并基于生死单元方法进行了数值仿真。根据靶材受热表面的气化压力反作用于熔融层的力学效应，建立了熔融液体克服表面张力运动的喷溅模型；分别用有限元法和有限差分法对二维轴对称的热传导方程和流体的一维欧拉运动方程进行耦合求解，得到了不同时刻模型的温度场和流场；应用生死单元法使发生喷溅的熔融液体单元在计算中失效，并得到了孔洞内壁和洞口表面的形状，其形貌与文献报道的实验结果相符合。数值结果显示出洞底附近的孔洞内壁具有不规则凹凸，而洞口则呈喇叭状，这是打孔过程中熔融液体在洞底附近的熔融层较薄处发生流动中断和位于洞口两侧的液体间张力作用所致。

The simulation of long pulsed laser drilling process is done by using the element birth and death method of ANSYS. According to the recoil force of evaporation on the heated surface, a model of ejection of melted materials is built. After coupled solving the two-dimensional axisymmetric differential equation of heat conduction and one-dimensional Eulerian equation of fluid motion by finite element and finite difference method respectively, the results of temperature field and flow field with temporal variation are obtained. The whole shape of the hole will appear by excluding the expelled element during computation, and the result accords with the experiment in literature. The simulation result explains that the irregular shape at the bottom of the hole and a trumpet shaped entrance appearing in the experiment are induced by the interruption of melted liquid within the thin melted depth and the tension between ejected liquid and melted materials adjacent to entrance， respectively.