Stefan Problem
斯特凡问题
在铸造模拟中听到了"斯特凡问题"这个词，它和普通热传导有什么区别吗？
斯特凡问题的理论基础
斯特凡问题的定义
我在CAE论文中看到了"斯特凡问题"这个术语，这只是指包含相变的热传导问题吗？
严格地说，这是指伴有移动边界的相变问题的经典数学模型。例如，冰融化成水时，固相和液相的边界（相界面）会随时间移动。这个移动边界的位置本身是未知数，在该位置施加潜热吸收·释放条件（斯特凡条件）是其特点。
"斯特凡条件"具体用什么公式表示？它和普通热传导方程有什么区别？
在每一相内，支配方程是普通的非稳态热传导方程，例如
相变发生的温度是恒定的吗（例如0°C）？现实中的材料不是纯物质吧。
很好的指出。经典斯特凡问题假设纯物质的等温相变。但在实际应用中，大部分是合金凝固问题，它们有"凝固温度范围"。此时采用"焓法"或"温度恢复法"等近似求解方法。例如，铝合金A356的凝固范围约为575°C到615°C。在这个温度范围内，固相和液相共存的"糊状区"中，边界不再尖锐。
斯特凡问题的数值计算方法
移动边界的处理方式
在数值计算中追踪移动的边界，是否需要改变网格本身？
主要有两种方法。一是"界面追踪法"，改变和重新构造网格以适应界面。二是"界面捕捉法"，在固定网格上用辅助变量（如水平集函数或VOF函数）表示界面位置。界面追踪法精度高，但大变形时网格质量容易恶化。铸造模拟通常采用后者的方法。
用固定网格捕捉界面时，潜热如何处理？当界面穿过一个网格单元时，该单元的热容量怎么定义？
这时经常使用"焓法"。用包含温度和相变的焓
焓法的计算中，数值上不会振荡吗？
会的。这是个大问题。由于温度-焓关系接近不连续，迭代计算中会发生"超调"，导致单个时间步中单元过度融化·凝固而产生发振。为防止这种情况，开发了"焓热源项法"和"温度恢复法"等技术。在实现时需要以"隐式"方式处理潜热，并仔细监控迭代收敛。许多商用求解器默认内置了这种稳定化处理。
斯特凡问题的实际应用
模拟设置的要点
实际进行铸造或熔融模拟时，材料物性中特别需要注意什么？
有三点很重要。第一是"潜热"的值。要从JIS或ASM手册等信得过的数据库获取。第二是"固相率-温度关系"。除非是纯物质，否则不是线性的，需要基于Scheil公式或杠杆规则的数据。第三是"凝固·熔融温度范围内的热导率和比热"。这些通常随固相率变化，使用常数会导致误差。例如，镁合金AZ91的凝固范围约为470°C到595°C，在这个范围内物性变化很大。
网格的细度需要多少？只在界面附近细化就可以吗？
因为事先不知道界面位置，整个可能发生相变的区域都需要一定的细度。根据经验，相变区域厚度方向至少需要5～10个单元。例如，如果预期凝固层10mm，网格尺寸应在1～2mm以下。此外，温度梯度陡峭的金属模附近或冷却水管周围还需要进一步细化。网格依赖性验证是必须的。粗网格中潜热释放速度偏快，凝固完成时间往往比实际早。
初始条件设为整体都是熔融温度时，金属模接触面的传热系数怎样确定？
这是不确定性最大的参数之一。在铸造中，金属模涂层有无、金属模温度、压力都会使其大幅变化。例如，用压铸法铸造铝合金时，金属模接触面的传热系数初期非常高（10,000 W/m²K以上），但凝固收缩产生空隙后急剧下降（1,000 W/m²K以下）。没有实测数据时，使用文献值或类似工艺的推测值，并进行灵敏度分析来评估对结果的影响是现实的做法。日本铸造工程学会（JFS）规格等也有参考数据。
斯特凡问题的软件比较
各求解器的方法
用Ansys Fluent或COMSOL求解斯特凡问题时，用户是否需要直接设置界面条件？
不需要。通常使用更高层次的功能。Ansys Fluent中只需启用"凝固/融化"模型，背后就会自动运行焓法。用户只需为材料定义潜热和固液相线温度即可。COMSOL Multiphysics中向"热传导模块"添加"相变"节点，类似地设置物性即可。两者都没有直接输入经典斯特凡条件的界面，而是实现了实用中的近似求解方法。
专用铸造模拟软件（如FLOW-3D CAST或ProCAST）和通用CFD软件的凝固融化模型根本上有什么区别？
专用软件具有针对铸造工艺特化的物理模型和完善的前后处理。例如，ProCAST（ESI集团）标准配备了预测缩孔所需的"二相区达西流模型"和"CAFE法"微观组织预测。还擅长考虑金属模加热循环和铸液注入工艺的连铸模拟。而通用软件如Fluent能计算凝固融化的基本过程，但要重现铸造特有现象需要用户自己构建追加模型，难度很高。
用Abaqus做热应力耦合分析时，相变引起的体积变化怎样处理？
在Abaqus/Standard中，一般先在"热电"步计算温度场，再在"结构"步施加热应变进行弱耦合分析。相变引起的体积变化（如铸铁石墨化膨胀）可用用户子程序"UEXPAN"或"UMAT"定义为温度和相（固相率）的函数，从而表示。也可用Abaqus的"模型改变"功能，动态改变凝固单元的物性（热膨胀系数）。无论如何，都需要超越纯粹斯特凡问题的材料工程知识。
斯特凡问题的故障排除
计算发散·振荡时
凝固模拟中，计算有时会发散或温度围绕融点剧烈振荡无法收敛。原因是什么？
最常见的原因是"时间步过大"。潜热的释放·吸收是巨大的热源项，显式处理会导致不稳定。对策首先是尝试减小时间步（如从0.01秒改为0.001秒）。其次检查求解器的"松弛系数"。将温度或焓的松弛系数从默认值1.0降低到0.5或0.7有时能稳定计算。Ansys Fluent中还可通过调整"凝固/融化参数"的"糊状区参数"（将默认值10^5增加到10^6～10^7），提高凝固区域的伪粘性，防止发散。
计算能收敛，但凝固完成时间明显比测量值早（或晚）。应该检查哪里？
第一，重新确认"潜热值"和"固液相线温度"。数据表中的值可能是99.99%高纯度物质的，与实际合金不同。第二，"界面传热系数"。太小则热散失困难，凝固延迟；太大则凝固加速。第三，"网格尺寸"。如前所述，过粗会使热扩散加快，凝固提前完成。需要进行网格依赖性验证，通过逆向计算来调整界面传热系数使其与测量值吻合（标定）。这在实务中常见。
"糊状区"比预期宽得多。这是物性设置错误吗？
很可能是"固相率-温度关系"的设置不恰当。纯物质设置（固相率在熔点瞬间从0跳到1）时，数值扩散影响会使糊状区看起来不自然地扩大。合金情况下，必须输入基于平衡状态图的准确固相率曲线。另外，如果没有将热导率设为固相率的函数，糊状区内的传热会不准确，导致区域不自然扩大。查看软件手册，确认固相率和物性联动设置的方法。
参考价值
更详细
错误