Numerical homogenization of a heat transfer problem with phase transitions in heterogeneous frozen soils

在广袤而寒冷的北方，大地不仅仅是泥土；它是一个结构性的基础，支撑着世界上规模最大的城市、管道和道路。几个世纪以来，这片多年冻土一直保持着坚固状态，其冰冻状态为人类建筑提供了稳固的岩床。然而，随着气候变暖，这片冻土正开始融化。当土壤中的冰融化时，地面会失去强度，导致不均匀的沉降和建筑物的开裂。精确预测这种融化会达到多深，关系到安全与工程设计，但土壤本身却提出了一个巨大的复杂性问题。真实的土壤并非一块均匀的冰，它是由不同材料组成的层、透镜体和囊状物构成的混乱混合物，其属性随距离的变化可能小至几厘米。要完美地模拟如此复杂的地面融化，计算机需要追踪每一个细微的层面，这项任务所需的计算能力超出了目前进行大规模、长期预测的能力。

为了解决这个问题，研究人员转向了一种称为“数值均质化”（numerical homogenization）的技术。想象一下你在预测水流过海绵的过程。你不需要模拟每一个微小的孔洞或纤维，而是将海绵视为一种具有“平均”渗透能力的单一材料。这就是均质化的本质：用一个更简单、经过平均处理的模型来取代混乱且细节丰富的现实，从而使模型运行得足够快，能够由计算机执行。然而，当地面正在融化时，这种“平均化”变成了一个陷阱。冻土与融化土壤之间的边界——即“前缘”（front）——并不会整齐地停在计算机计算块的边缘，而是会直接切穿其中。这产生了一个困境：如果一个土块一半是冻结的，一半是融化的，计算机应该如何决定它的属性？它应该表现得更像冰，更像水，还是介于两者之间？这种选择将决定模型的预测结果是建筑物会下沉还是屹立不动。

在最近的一项研究中，来自北东联邦大学（North-Eastern Federal University）的 Timur Zakharov 针对处理这些“混合”土块的具体困境，测试了三种不同的处理方式。研究重点关注了一座坐落在冻土上的加热建筑，在这种情景下，来自结构的温暖会缓慢融化下方的土壤。团队创建了这个环境的“数字孪生”，使用一种真实的、随机分布的土壤导热模式来模拟自然界中地面的混沌状态。随后，他们在粗糙的网格上运行了模拟，在这些网格中，每个土块的尺寸都远大于实际的土层和冰层，这迫使相变前缘（phase-change front）切穿了这些土块的中间。他们对比了三种用于决定这些被切分的土块内部情况的不同方法。

第一种方法被称为“表格法”（table approach），它预先计算出每种可能的土壤状态下的属性列表，并根据土块的平均温度来查找答案。第二种方法是“节点方案”（nodal scheme），它通过确定土块顶角的土壤状态并对这些决策进行平均来做出判断。第三种方法称为“h-平滑法”（h-smoothing），它根据线性温度梯度计算土块融化的精确比例，本质上是在测量土块中已经变成水的切片部分。研究人员在二维和三维空间中都运行了这些模拟，并将结果与一个追踪了每一层土壤的极其详细的“参考”解进行了对比。

研究结果揭示了一个令人惊讶的权衡关系，这对工程师如何选择工具提出了挑战。当目标是预测土壤温度时，h-平ч滑法被证明是最准确的。它以最高的精度追踪热流，尤其是在其他方法难以应对的粗糙网格上。然而，当研究人员观察融化前缘（即冻土与融化土壤之间的实际边界）的位置时，情况却完全改变了。h-平滑法虽然在温度方面表现出色，但它始终高估了融化的深度。它预测的地面融化程度比实际情况更深。

这种差异的原因在于能量守恒定律。h-平滑法在努力精确描述融化区域形状的过程中，不经意间在每个时间步长中向系统添加了少量的额外热量。在长达一年的时间内，这种微小的盈余不断累积，导致模拟的地面比真实的地面融化得更多。相比之下，另外两种方法，特别是表格法和节点方案，是“保守型”的。它们既不增加也不减少能量，因此即使它们的温度读数稍欠精确，也能更准确地预测融化深度。研究表明，如果研究人员需要知道土壤的确切温度，h-平滑法是更优的选择；但如果目标是了解由于融化导致的地面沉降深度，那么保守型方法则更为可靠。

研究还探讨了当土壤熔点在不同层之间不一致时的情形。在这种更复杂的情景下，使用详细属性查询的表格法在预测融化前缘位置方面，表现优于较简单的双张量（two-tensor）方法。这表明，对于高度多变的地面条件，需要一个更详细的预计算土壤状态图。最终，这项研究为工程师和科学家提供了一个清晰的指南：不存在单一的“最佳”方法。选择取决于所提出的问题。如果优先考虑热力学精度，一种方法胜出；如果优先考虑融化的物理边界，另一种方法则占据领先地位。通过澄清这些差异，这项工作有助于确保用于保护北极地区建筑的模型不仅运行迅速，而且值得信赖。