Dear mathias,
Thanks for your reply,
The problem has been solved by changing dynamics and reducing relaxation time.
May I aks you an another question, which is about boundary condition.
I want to computing phase change process in 3D structure with a cylinder heating surface. The boundary condition of “setRegularizedTemperatureBoundary” in 2D should be changed to which kind of boundary condition in 3D.