二维二相Stefan问题的数值模拟

本文利用线方法和变区域的有限元法,成功地实现了二维二相Stefan问题的数值模拟.本方法的优点是运动边界上的两个条件被同时满足,因此,它具有较高的精度,即使在粗网格情况.本文把一个正方形冰块的熔解问题作为典型的计算例子.在各种不同参数(热扩散系数、熔解潜热等)和两类边界条件下,给出了温度分布和分界面位置的数值结果.在退化的单相情况下,与Meyer的结果进行了比较,符合是很好的.

NUMERICAL SIMULATION OF TWO-DIMENSIONAL,TWO-PHASE STEFAN PROBLEMS

Using the method of line and the variable domain finite element method, the numerical simulation of 2-D., two-phase Stefan problems has been carried out successfully. The advantage of the method is that two boundary conditions on the moving boundary are simultaneously satisfied. Honce, it has a higher precision even for the coarse meshes.A melting problem of a square ice will be as a typical calculating example. For the different parameters (the thermal diffusivity, the latent heat of melting etc.and two types of the boundary condition, the numerical results of the temperature distribution and the position of the interface have been obtained. For the degenerated one-phase case, our results are compared with Meyer's results7]. The agreement is quite good.