《Communications in Nonlinear Science and Numerical Simulation》：Numerical investigation of the 2D unsteady natural convection heat transfer equation with tempered fractional constitutive relationship
编辑推荐：
本文针对复杂流体（如纳米流体、多孔介质）中经典傅里叶定律难以准确描述异常热扩散与记忆效应的问题，提出并发展了一种求解二维非定常自然对流传热方程的高阶全离散间断伽辽金(DG)方法。该方法通过缓变分数阶导数对傅里叶定律进行推广，以更精准地刻画非局部和记忆性传热现象。研究人员构建了高效的缓变分数阶导数有限差分格式，在此基础上建立了稳定且收敛的全离散DG方案，并采用分级时间网格以精确捕捉初始时刻的弱奇异性。数值实验验证了该方法在求解复杂流体系统传热问题中的高效性与高精度，为模拟具有复杂微观结构和历史效应的热输运过程提供了强有力的数值工具。
在自然界和工程领域，从大气的环流到微纳尺度散热器的设计，流体因温度梯度驱动的密度变化而产生的自然对流传热无处不在。然而，在描述如纳米流体、异质多孔介质等复杂材料中的热传导时，经典的傅里叶定律（即热流与瞬时温度梯度成正比）常常“失灵”。因为这些系统中的热输运不仅关乎“现在”的温度差，还受到“过去”热历史的深刻影响，并且热量能在空间中跨越较远距离进行非局部的传递，表现出“异常扩散”的特性。为了捕捉这些复杂的记忆效应和非局部相互作用，研究者们引入了分数阶微积分，它像为传统微积分增添了“历史滤镜”，能更好地描绘这种具有记忆性的物理过程。但标准的分数阶模型有时会“记住”太多太久的历史，导致物理上不切实际的无限长记忆或方差，这促使了缓变分数阶导数的诞生——它在分数阶导数的基础上加入了一个指数阻尼项，既能保留关键的短到中期记忆，又能避免物理上不合理的极端“长尾”效应，从而更贴合真实世界中热输运的衰减记忆行为。本文发表于《Communications in Nonlinear Science and Numerical Simulation》，由Leilei Wei、Libo Feng、Ian Turner、Zhi Mao和Fawang Liu合作完成，他们聚焦于发展一套高精度的数值“工具箱”，来求解基于这种更真实的缓变分数阶本构关系的二维非定常自然对流传热方程。
为了攻克这一数值模拟难题，研究人员主要运用了三个关键技术。首先，他们设计了一种高效、高精度的有限差分格式来逼近关键的缓变分数阶时间导数。其次，他们构建了一套基于分段多项式空间的全离散间断伽辽金（DG）方法，该方法在空间上对控制方程进行离散，其核心优势在于能灵活处理复杂几何形状、适应非均匀网格，并能通过局部高阶多项式近似实现高精度。最后，为了精准捕捉方程在初始时刻（t=0）附近解可能出现的弱奇异性，他们采用了非均匀的、在初始点附近加密的分级时间网格，并相应地提出了改进的全离散DG格式，从而确保了数值解在全时间域的准确性。
稳定性与收敛性分析：研究人员首先从理论上对所提出的全离散DG方法进行了严格的数学证明，证实了该数值格式是无条件稳定的，并且其数值解会收敛到真实解。这为方法的可靠性和鲁棒性奠定了坚实的数学基础，意味着无论时间步长和空间网格如何选取，计算过程都不会因误差积累而失控，且结果能随着网格加密而逼近真解。
数值实验验证：为了展示方法的有效性，论文进行了系统的数值实验。通过构造具有精确解析解的问题，他们计算了数值解与精确解之间的误差，并观察了误差随着网格加密（h-精化）和多项式阶数升高（p-精化）而减小的规律。结果表明，所提出的DG方法能够达到理论预期的高阶收敛精度，验证了格式的正确性。这些基准测试为方法应用于更复杂的实际问题提供了可信度。
参数影响研究：研究深入探讨了缓变分数阶模型中的两个关键参数——分数阶α和缓变参数γ对传热行为的影响。数值模拟揭示，较小的α值增强了系统的“历史记忆”，导致热量扩散得更慢，温度分布变化更为平缓；而较大的γ值则意味着记忆效应衰减得更快，相当于“遗忘”了过去更早的热历史，从而使热响应更接近经典（无记忆）扩散行为。这种参数敏感性分析不仅验证了模型捕捉物理现象的能力，也为通过调整α和γ来匹配特定材料（如特定纳米流体或多孔介质）的实验观测数据提供了可能。
研究的结论与讨论部分强调了这项工作的核心贡献与重要意义。首先，论文成功构建了一个用于求解二维缓变分数阶自然对流传热方程的高阶全离散DG框架。这个框架是“高配”的：它基于分段Pr多项式，意味着可以在空间上实现任意高阶精度；它采用了全离散格式，将时间和空间离散紧密耦合，确保了整体计算的协调性。其次，研究并非“黑箱”操作，而是提供了严格的稳定性与收敛性证明，这像为数值方法颁发了“质量合格证书”，保证了其在科学计算中的可靠性与数学严谨性。再者，所提出的方法具有高度的灵活性与扩展潜力。尽管本文为简化分析聚焦于线性模型，但所建立的DG框架完全有能力处理更普遍的非线性自然对流传热问题。文中提到，非线性项可以自然地嵌入当前离散结构中，并可通过高效的迭代算法（如雅可比自由牛顿-克里洛夫法）求解，这为模拟真实世界中温度依赖热导率等非线性效应铺平了道路。
这项研究的最终意义在于，它为解决复杂流体和材料中的非经典传热模拟提供了一套强大、可靠且高精度的数值工具。传统的数值方法在处理此类具有非局部性和记忆效应的分数阶模型时，往往面临精度不足、理论分析缺失或难以处理多维问题等局限。本文的工作不仅突破了这些限制，其采用的DG方法还能轻松应对复杂几何和非均匀网格，非常适合于实际工程应用中常见的非规则计算域。通过缓变分数阶导数，模型能够更真实地描述纳米流体内部粒子-流体相互作用、多孔介质中因结构约束导致的亚扩散到正常扩散的转变等现象。因此，该研究对于深入理解先进能源系统、环境工程、生物医学工程等领域中涉及异常热输运的过程，以及优化相关设备的设计与性能，具有重要的理论价值和应用前景。
