金属凝固传热仿真的核心难点与潜热处理算法解析

我做金属凝固过程的传热学仿真最早是从一个很简单的想法开始的把“金属凝固”当成一个带冷却的瞬态导热问题去解。真正把模型跑起来之后结果却让我怀疑人生——按一个普通热传导问题算出来的中心温度在十几秒内就掉到了固相线以下而现场的工艺数据告诉我这个位置的凝固至少还要持续几分钟。问题出在哪出在我把凝固过程的难点想简单了。金属凝固数值模拟和普通导热分析最大的区别就是它必须在能量守恒里处理“相变潜热”而潜热在数学上会带来一个非常陡的非线性项处理不好算出来的温度场完全是错的。这篇文章我就从传热学仿真的角度把金属凝固过程数值模拟里最核心的物理模型、潜热处理算法、完整的算例搭建过程以及我在调试中踩过的坑全部讲一遍。无论你是刚开始接触凝固模拟的学生还是已经在用商业软件做铸件、焊件热仿真的工程师这篇文章都能够帮你把底层逻辑理顺。1. 凝固传热和普通导热模拟的差距藏在三个“低温细节”里很多人在刚接触时都会问凝固不就发生在液相线温度和固相线温度之间吗我只要在这个温度区间多放一点“热容量”不就行了思路没有错但在具体执行时有三个细节会让仿真结果从“差不多能用”变成“差得离谱”。1.1 潜热释放比过热度显热大一个数量级必须当成能量储存来处理我先用一个简化但很能说明问题的例子来算一笔账。以某铝合金为例比热容约900 J/(kg·K)凝固潜热按390 kJ/kg计。假如浇注温度比液相线高40 K那么过热段能释放的显热是900×4036000 J/kg也就是36 kJ/kg而凝固过程中潜热释放量是390 kJ/kg是过热显热的十倍以上。更夸张的是金属从液相线继续冷却到室温时全部显热加起来也不一定比潜热多多少。这意味着什么呢意味着如果你只把潜热简单当成一个“在冷却过程某一段释放的热量源项”去处理数值上是很容易出问题的。潜热释放不是一个固定的集中热源它取决于局部固相分数的变化速率也就是ρ ∂H/∂t ∇·(k∇T)这里的焓H包含显热和潜热两部分。潜热的多少由固相分数fs决定而fs是温度的函数。比如在一个二元合金里液相线到固相线之间温度降低固相分数从0逐渐增加到1这个耦合关系非常强。如果不把焓作为一个独立的状态变量求解而只盯着温度那么数值格式在凝固区间里就很容易失去守恒性。我见过不少初学传热学仿真的同学在商业软件里把“潜热”填在材料表里就以为万事大吉结果中心冷却曲线严重偏快位置越靠中心的点偏得越厉害。核心原因就是他们没有意识到凝固过程中的焓不是一个只跟“当前节点有没有凝固完”相关的量它还和凝固时间跨度、空间分辨率和数值耗散全部耦合在一起。1.2 固液两相的物性变化让问题从一开始就是强非线性的凝固仿真的第二个难点是物性。液相的导热系数和固相不同密度也不同比热容在相变区间附近也不是常数。以钢来说固态时的导热系数大约在30~50 W/(m·K)液态时虽然本身的导热系数不高但流动性带来的自然对流使得等效换热能力远高于单纯导热。这种非线性的物性变化会让“温度场”和“固相分数场”相互纠缠。温度决定固相分数固相分数反过来决定热容、导热系数和界面热阻。如果我们抱着“常温导热模拟”的惯性把材料参数设置成常数或者只是粗略地做个线性插值那算出来的凝固时间、温度梯度、糊状区分布基本上是没有什么定量参考价值的。对做铸件凝固工艺的人来说我们最关心的往往不是某一个点的温度而是最后凝固的液相孤立区在哪、糊状区的形态怎么演化。这些东西对物性模型的敏感程度远高于对边界条件中辐射率的敏感程度。所以在做金属凝固过程数值模拟之前要把材料物性随温度的变化整理清楚。不需要强求实验室级别的精确度但至少固相、液相、固液两相区的三个典型值必须给出来并设置成随固相分数平滑过渡。2. 做凝固仿真前先定下这三件事相变区间、液相流动、界面换热很多新手一上来就找潜热处理算法或者急着调时间步长。我的建议是反过来先把下面三个物理问题定下来。这三个问题不定清楚后面无论用什么算法都是空中楼阁。2.1 相变区间模型固相分数和温度的关系不是随便画条线就行合金的凝固过程通常在一个温度区间里完成这个区间由液相线温度和固相线温度确定。但在这两条线之间固相分数随温度的演化规律并不是唯一的常见处理有三种。第一种是线性模型认为固相分数从液相线到固相线均匀增加。这个模型最简单计算成本低适合做基本工艺判断。第二种是杠杆定律模型它假设固相和液相内部完全扩散均匀固相分数随温度以非线性的方式变化。这种模型适合分析凝固速度较慢、扩散时间充分的场合但多数快速凝固工况其实不满足假设。第三种是Scheil模型假设固相内没有扩散而液相内完全混合这个模型在铸件凝固模拟中很常用对于枝晶偏析严重的合金更接近实际情况。在传热学仿真里固相分数曲线确定之后潜热释放的速率分布也就确定了。同样的凝固区间采用不同固相分数模型虽然总潜热一样但释放的时机和峰值速度不同最终的温度场和冷却曲线也会有差别。如果项目目标是做宏观温度场评估用线性模型通常够用如果后续要做微观组织预测或者缩孔分析建议至少采用Scheil模型。2.2 液相对流到底算不算算不清楚就先用有效导热系数金属液在凝固过程中并不是静止的。冷壁附近的金属液温度低、密度大会向下流动形成自然对流如果有电磁搅拌或者浇注带来的初始流动那就是强制对流。液体的对流会显著改变凝固前沿的温度分布尤其是在铸件上表面和侧壁附近对流强的地方壳层不容易长厚。在一个完整的物理过程中凝固传热仿真可以细分为传热、流动、传质三部分。但如果项目目标是偏宏观的传热学仿真通常不会直接去求解液相区的Navier-Stokes方程因为液固两相区的流动、糊状区渗透率、结晶间流动模型都会让计算变得极其复杂稳定性也更难控制。一个折衷方案是采用有效导热系数法也就是把液相自然对流对热量输运的贡献近似等效成一个比纯导热大得多的导热系数。具体放大多少倍要根据熔体的Pr数、Grashof数和容器几何形状来估计没有统一常数多数工程实践里取3到10倍。做参数敏感性分析时可以把这个倍数作为一个不确定变量先看它对凝固时间的影响有多大。如果影响很大说明你的工况不能省掉流动模拟如果影响很小那就可以安心用传热模型。2.3 铸件—铸型界面的换热系数比合金物性更值得花时间做金属凝固传热学仿真最容易栽跟头的其实不是材料热物性而是铸件和铸型之间的界面换热系数。这个数值的波动范围极大。举个具体的例子铝合金砂型铸造时铸件和砂型之间的界面换热系数可能只有几百W/(m²·K)而入金属模或水冷铜模时这个换热系数可能上千甚至数千W/(m²·K)。同一个铸件换热系数差一倍凝固时间可能差出30%以上。很多模拟者喜欢把这个量设成一个常数然后用“反算标定”的方式去拟合实测热电偶曲线。如果能拟合得好那当然没问题。但要注意随着凝固收缩金属外壳会和铸型脱开形成气隙气隙越大界面换热系数下降越明显。所以界面换热系数其实是随凝固时间变化的。在比较精细的仿真里需要把气隙厚度作为反馈量联动修改界面换热系数或者至少用分段常数近似。我的建议是如果项目还没有可靠的界面换热系数数据不要花太多时间去折腾潜热算法。先用一组典型参数把模型跑通再结合现场测温反算这个换热系数这样循环几次模型才能真正有预测价值。3. 潜热处理的三类算法实测等效热容、焓法、温度回升法的真实取舍物理模型确定之后棘手的问题才真正到来怎么在数值计算里处理潜热我在实际项目中用过三种主流方案下面把它们各自的实现思路、优缺点和应用场景都摊开来说。3.1 等效热容法能跑通但不要为了收敛去拉宽凝固区间等效热容法的逻辑非常直观。在固液两相区内温度每变化1 K材料除了显热比热之外还会释放一部分潜热于是可以把凝固潜热折算成比热容的一部分cp_eff cp L / (TL - TS)这个公式理解起来毫无门槛。实现方式就是在材料属性表里把固液两相区的比热容设为一个峰值更大的值。很多商业软件自带这个功能用户只需要输入潜热量和凝固区间就行。但等效热容法有一个非常隐蔽的坑如果凝固区间很窄比如纯金属或共晶成分合金只有1~2 K的凝固温度范围那么这个等效比热可能比实际比热大上百倍。在数值计算里网格节点温度刚刚跨过液相线时它会瞬间“吃掉”大量热量导致温度被压低到固相线附近凝固平台的温度位置被严重扭曲。为了缓解这个问题有的人会人为把凝固区间拉宽比如把20 K的真实区间改成200 K。这样计算是稳定了但糊状区的空间范围被数值误差人为放大后面所有基于糊状区位置的分析都会失真。所以使用等效热容法时我的经验是凝固区间内部至少要覆盖3到5个网格时间步长要保证单个步长内温度变化不超过凝固区间宽度的1/10同时尽量不要为了