COMSOL中的固体瞬态加热建模
适用于任何涉及瞬态加热的情况
将一个空间上均匀分布的热载荷施加在一个具有均匀初始温度的圆柱体材料顶面的圆形区域内。最开始载荷很高,但在一段时间后会逐渐下降。除了施加热载荷外,还添加了一个边界条件来模拟整个顶面的热辐射,它使零件重新冷却。假设材料属性(热导率、密度和比热)和表面辐射率在预期温度范围内保持不变,并且假设没有其他作用的物理场。我们的建模目标是用它来计算圆柱体材料内随时间变化的温度分布。图1为建模场景。
在图1中,可以看到几何体和载荷是围绕中心线轴向对称的,所以我们可以合理地推断,解也将是轴向对称的。因此,我们可以将模型简化为二维轴对称建模平面。在中间的圆形区域内,热通量是均匀的。最简单的建模方法是通过在二维域的边界上引入一个点来修改几何形状。这个点将边界划分为受热和未受热的部分。在几何形状上增加这个点,可以确保所产生的网格与热通量的变化完全一致。考虑到这些,我们可以创建一个等效于三维模型的二维轴对称计算模型(图2)。
此外,我们还考虑了施加的热通量大小的瞬时变化的情况;在 t=0.25s 时,它的值变得较低。我们可能也想知道求解器采取的时间步长,这可以通过修改求解器的设置,按求解器的步长输出结果,然后就可以绘制零件顶部中心点的温度,如图3所示。
接下来,我们用不同的求解器相对容差值重新运行该模型,并在图中进行比较(图4)。这类图表明,像预期的那样,随着公差变小,解迅速向同一个值收敛。
另一个可以计算的量是进入该域的总能量。我们可以对通过边界的总热通量的表达式 ht.nteflux 进行积分,使用 timeint() 算子对时间进行积分,得到总能量。积分的结果在下面的表格中列出,用于增加时间步长的相对容差。
从数据中我们可以观察到,进入系统的总能量实际上几乎与时间步长容差无关。乍一看,这似乎是对我们模型的一个奇妙的验证。然而,需要指出的是,我们在这里观察到的是有限元法(FEM)的基本数学特性。简单说,就是总能量总是会很好地平衡。这并不意味着模型中没有错误,错误只是出现在不同的地方……接下来,我们就去寻找错误。
从哪些角度去寻找错误呢,答案是:输入错误、离散化错误、时间步长误差、插值错误等
到目前为止,我们已经观察了模型中某一点的解,并观察到随着我们完善瞬态求解器的相对容差,解似乎收敛得很好,现在,我们来看看空间温度分布。我们将从沿中心线的温度开始,对于最宽松的容差 1e-2,看初始时间的解以及求解器采取的第一个时间步长,如下图所示。
从初始值图中,我们可以看到,沿中心线的温度与规定的初始温度不一致——有些地方甚至低于初始值。这是由于COMSOL使用了所谓的一致初始化。
在这种情况下,考虑一致初始化的方式是,调整温度场使之与施加的载荷和边界条件相一致。
简单来说,为了使这个解的离散化误差最小化,我们需要在场变化剧烈的位置划分更精细的网格。根据我们的经验,在非常接近表面和边界的法线方向上,场的变化非常大,但在内部则变得更加平滑。这正是需要边界层网格划分的情况,如图6所示,它在边界的法线上创建薄的单元。
现在,我们可以重新运行模拟,并绘制出初始时间和下一个时间步长的解。
在图7中,我们可以观察到,在初始值时,温度的下调在空间上更为局部。因此,通过这种细化的网格,我们减少了空间离散和时间步长误差。我们还可以看一下沿建模域顶部边界的结果,代表暴露表面的温度分布。图8中显示的是使用 1e-2 的容差绘制的初始时间和第一个时间步长。在这些图中,我们可以观察到空间中相当剧烈的震荡场。这是空间离散化的一个表现。
解与之前类似,但现在我们必须在过渡的位置细化网格。对于这个问题,可以对划线点应用更加精细的大小设置,从而得到如下所示的网格。
从图10的温度结果中,我们可以看到,现在解中的振荡已经减少了,在空间和时间上的传播也没有那么多。即使求解器的相对容差为 1e-2,求解的结果也已经有了很大的改善。
在这个例子中,施加在边界上的热载荷并没有及时移动,所以划分边界的方法是合理的。如果热载荷分布要移动,那么受热面的整个网格就需要更加精细。
我们假设材料属性随温度的变化而保持不变,并且不依赖于任何其他物理场。这是一个重大的简化,因为所有的材料属性都会随温度变化。材料甚至可以经历相变,如融化。
我们也可以很容易地预见到这是一个多物理场问题,如涉及热固化的方程式,甚至是材料非线性电磁加热问题。在这种情况下,我们不仅需要监测温度场的收敛性,还需要监测所有其他正在求解的场变量,甚至可能是它们的时间和空间导数。这些情况可能都需要在建模域的所有地方采用非常精细的网格。
综上,这是一个典型的传热建模问题。我们注意到,在空间和时间上,在载荷突然变化的情况下,解会出现某些错误。我们应该对这些错误类型有所了解,并知道它们是有限元方法的固有结果,就像所有的数值方法一样,只是对现实的一种近似。尽管这些误差看起来很大,但由于瞬态传热方程的扩散性,它们的大小在空间和时间上都会衰减。
COMSOL的学习也需要多看文献,可以根据文献中的方法去做一些尝试,另外还有一些不错的博客和公众号,像“COMSOL仿真交流”这个公众号我觉得挺不错的,上面分享了很多comsol的案例,可以多学习借鉴别人的方法。