一个求解Stefan自由边界问题的基于位势理论的直角网格方法