[论文解读] One-phase Stefan problem with a latent heat depending on the position of the free boundary and its rate of change

Julieta Bollati, Domingo A. Tarzia|arXiv (Cornell University)|Mar 21, 2017

Lattice Boltzmann Simulation Studies参考文献 14被引用 7

一句话总结

本文提出了一类一相Stefan问题的解析解，其中潜热依赖于移动自由边界的位移和速度，推广了经典及先前的变潜热公式。通过相似变量和合流超几何函数（Kummer函数），作者推导出一个显式解，当特定参数取值时可恢复已知结果，为具有位置与速率相关潜热的相变动力学提供了统一的数值验证基准与物理解读。

ABSTRACT

Fil: Bollati, Julieta. Universidad Austral. Facultad de Ciencias Empresariales. Departamento de Matematicas; Argentina. Consejo Nacional de Investigaciones Cientificas y Tecnicas. Centro Cientifico Tecnologico Conicet - Rosario; Argentina

研究动机与目标

• 通过将潜热表示为自由边界位移与变化率的函数，扩展现有的一相Stefan问题模型。
• 利用相似性方法与特殊函数，提供一个显式解析解，确保数学上的严谨性与普适性。
• 统一并推广文献中已有的结果，包括常数、线性及与速度倒数相关的潜热公式。
• 通过计算示例与使用牛顿法的数值参数估计，验证解的正确性。
• 为具有变潜热的Stefan问题数值模拟提供基准解。

提出的方法

• 建立一个一相Stefan问题，其中自由边界由热传导方程与Stefan条件控制，潜热为 L(s(t), ṡ(t)) = γs^β(t)ṡ^δ(t)。
• 应用相似性变换，将偏微分方程组转化为关于相似变量 η 的常微分方程。
• 利用Kummer的合流超几何函数求解所得的二阶常微分方程，通解以 M(a,b,z) 和 ηM(a+1,b+1,z) 的形式表示。
• 在 η = 0 和 η → ∞ 处施加边界条件，以确定常数及自由边界参数 ξ。
• 利用Kummer函数理论与互补误差函数重复积分的性质，确保解在整数与非整数参数下均有效。
• 采用数值求根方法（牛顿法）估计给定 β、δ 和 T₀ 值下的自由边界参数 ξ。

实验结果

研究问题

• RQ1 在引入位置与速率相关潜热后，一相Stefan问题中自由边界的动力学行为如何变化？
• RQ2 能否推导出一个统一的解析解，使其在 β 与 δ 变化时可作为已知文献解的特例？
• RQ3 当 β、δ 及初始温度系数 T₀ 变化时，自由边界参数 ξ 的行为如何？
• RQ4 该解如何保持与物理直觉的一致性，例如初始温度越高，相变越快？
• RQ5 该解析解在多大程度上可作为具有变潜热的Stefan问题数值方法的基准？

主要发现

• 推导出一相Stefan问题中潜热 L(s(t), ṡ(t)) = γs^β(t)ṡ^δ(t) 的显式解析解，适用于任意实数 β 与 δ。
• 当 β = 0 且 δ = 0 时，解退化为经典Stefan问题，恢复标准相似性解。
• 该解推广了先前结果：当 β = 1、δ = 0 时，可恢复线性潜热情形；当 β = 0、δ = -1 时，可恢复文献[15]与[16]中研究的逆速度依赖情形。
• 当 δ 固定时，ξ 随 β 增大而增加；反之亦然，表明 β 或 δ 值越大，相变越快。
• ξ 随 T₀ 增大而增加，证实了初始温度越高，相变过程越快的物理直觉。
• 通过牛顿法进行的数值验证表明，解与已知结果一致，且该解析解可作为具有变潜热的Stefan问题数值方案的可靠基准。