The methods of local mesh refined,dynamic time step and lumped heat capacity matrix were introduced to resolve the numerical oscillation problem.

• 为了解决数值振荡问题,采用了局部网格细化、动态时间步长和集中热容矩阵等方法。