多尺度热传导广义有限元方法的时间收敛性分析

《Computer Methods in Applied Mechanics and Engineering》：Temporal convergence analysis of the generalized finite element method for multi-scale heat transfer

　　

摘要：研究了具有时间相关富集函数（enrichment functions）的广义有限元方法（GFEM）在涉及局部、非平稳热载荷的多尺度瞬态热传导问题中的时间精度、收敛性和效率。GFEM能够通过分析或数值考虑构建针对特定解的富集函数，从而在粗固定有限元网格上同时解析空间和时间尺度。采用全局-局部广义有限元方法（GFEMgl）构建数值富集函数。使用广义梯形法（generalized trapezoidal method）作为时间积分器研究GFEM的时间行为，以评估在此框架内时间步进的可行性。针对具有尖锐时间梯度和高度局部化空间梯度的模型问题进行的数值实验表明，GFEMgl仅在使用瞬态局部问题时才能获得鲁棒收敛，这代表了相对于先前研究的进步。性能分析表明，与标准有限元法（FEM）公式相比，具有时间相关形函数的GFEM能够高效地提供高保真解。结果展示了GFEM在显著降低计算成本方面的巨大潜力。

论文解读

**研究背景与问题**

在现代工程中，许多关键问题涉及跨越不同空间和时间尺度的复杂多物理场行为，例如高速飞行器在极端环境下的气动热弹性响应、焊接力学中的局部加热、以及激光加热与成形过程。这些场景中存在强烈的、高度局部的、非平稳的热载荷，产生尖锐的温度梯度，对数值模拟提出了严峻挑战。传统有限元方法（FEM）在处理此类多尺度问题时，通常需要在整个计算域内进行广泛网格细化或随时间步重新划分网格，导致计算成本极高。因此，开发能够高效、准确模拟此类问题的数值方法成为迫切需求。

广义有限元方法（GFEM）通过将基于先验信息构建的富集函数（enrichment functions）引入近似空间，旨在缓解网格细化的需求。然而，现有研究主要集中于空间收敛行为，对时间收敛性的理解不足，特别是在三维瞬态热传导问题中。此外，先前的GFEM应用多局限于与计算网格对齐的静态热载荷或移动界面问题，忽略了具有任意空间取向的平移热载荷。全局-局部广义有限元方法（GFEMgl）通过求解局部问题来动态构建富集函数，但以往研究仅考虑稳态局部问题，其时间收敛行为在强时间梯度下尚不明确。因此，本研究旨在填补这一空白，系统表征GFEM在三维多尺度瞬态热传导问题中的时间收敛性。

**研究内容与结论**

研究人员对GFEM与时间相关富集函数（time-dependent enrichment functions）在高度局部化、强瞬态热载荷下的时间收敛性进行了系统分析。采用广义梯形法（generalized trapezoidal method）作为时间积分器，分别考察了后向欧拉法（BE）和克兰克-尼科尔森法（CN）的表现。通过两个模型问题开展数值实验：第一个问题采用解析解构造的指数富集函数，第二个问题模拟三维光束上施加高斯分布热流，使用GFEMgl构建数值富集函数。研究得出以下结论：GFEM结合时间相关、针对解定制的富集函数能够在粗网格上实现时间收敛，且收敛率与理论预期一致；GFEMgl仅在使用瞬态局部问题（TL-GFEMgl）时才能在强时间梯度下获得鲁棒的时间收敛行为，而使用稳态局部问题（SS-GFEMgl）会导致误差饱和和收敛率下降；与标准FEM及高阶多项式GFEM相比，GFEMgl在保持同等精度的前提下，可将自由度降低高达95%，显著减少计算成本。该论文发表在《Computer Methods in Applied Mechanics and Engineering》。

**主要关键技术方法**

研究人员采用了以下关键技术方法：
1. **广义有限元法（GFEM）**：通过引入基于先验信息或局部问题解构建的富集函数，在粗固定网格上解析多尺度特征。
2. **全局-局部广义有限元法（GFEMgl）**：通过求解局部边界值问题（local boundary value problems）动态构建数值富集函数。局部问题可采用稳态或瞬态热传导方程，使用高分辨率网格求解。
3. **时间离散化**：采用广义梯形法（后向欧拉法和克兰克-尼科尔森法）进行时间步进，并特别处理了时间相关形函数带来的离散化问题。
4. **误差分析**：使用基于内部能量的时间平均相对L2误差（L₂Error(U(t))）作为时间精度度量，确保空间误差最小化以聚焦于时间收敛行为。

**研究结果**

**6.1 一维方向尖锐瞬态梯度的制造解问题**
- **空间离散化影响**：线性FEM和二次GFEM在粗网格上无法解析热尖峰，而二次+指数GFEM（即使用时间相关解析富集函数）能够准确捕捉并预测热尖峰的运动，实现时间收敛。
- **时间收敛性分析**：对于二次+指数GFEM，后向欧拉法达到线性收敛率（~0.998），克兰克-尼科尔森法达到近似二次收敛率（~2.42），与理论预期吻合。网格尺寸（h_α）和体积热容（ρc）对收敛率影响较小，但较小时间步长需先解析热尖峰运动。

**6.2 三维光束上施加尖锐热流的数值实验**
- **参考解确定**：通过hp-GFEM（高阶多项式GFEM）进行网格和时间步收敛性研究，确定空间和时间参考解，精度约O(10^-4)。
- **通量角度影响**：对于稳态热流，通量角度（控制y方向梯度）对时间误差影响极小，GFEMgl与hp-GFEM的收敛率与理论值一致（BE≈1，CN≥2），但SS-GFEMgl在CN下出现误差饱和。
- **通量速度影响**：平移热流速度增加时，后向欧拉法的误差增大，需更小时间步才能达到渐近收敛；克兰克-尼科尔森法对速度不敏感，但SS-GFEMgl在CN下出现误差饱和。GFEMgl的收敛率不受速度影响，仅影响达到渐近收敛所需的时间步。
- **时间上升系数γ影响**：γ增大（热流增长更快，时间梯度更强）导致相同时间步下误差增大。SS-GFEMgl在CN下收敛率显著下降（如γ=∞时收敛率约0.6），而hp-GFEM仍接近理论率。通过使用瞬态局部问题（TL-GFEMgl），所有情况下均恢复最优或接近最优收敛率，误差饱和水平降低一个数量级。
- **计算成本分析**：GFEMgl全局自由度相比hp-GFEM减少约3-4个数量级，局部问题自由度也远少于全局细化所需。TL-GFEMgl在精度上与hp-GFEM相当，而SS-GFEMgl在CN下精度较低。GFEMgl通过并行求解局部问题进一步降低计算时间。

**讨论与结论**

研究结论部分翻译如下：本文研究了使用一般时间积分器的GFEM在多尺度热传导问题中的时间收敛性。结果证明，时间积分器的行为取决于所考虑的参数空间和热方程的形式。一旦解析了底层特征的特征时间尺度，后向欧拉法和克兰克-尼科尔森法均表现出鲁棒行为，达到预期收敛准则。本研究考虑的GFEMgl策略能有效捕捉广泛细尺度特征，但使用稳态局部问题构建全局-局部富集函数（SS-GFEMgl）时，在更强烈的时间梯度下会导致误差增大和收敛率下降。通过使用瞬态局部问题（TL-GFEMgl）改善了该行为，使得尖锐时间梯度能够被捕捉。因此，只要选择适当的富集函数来近似局部特征，具有针对解定制、时间相关富集函数的GFEMgl能够生成鲁棒且精确的解。总体而言，具有针对解定制、时间相关富集函数的GFEM结合后向欧拉法和克兰克-尼科尔森法，在研究的算例中相比标准有限元方法取得了显著的计算优势。此类公式缓解了多项式近似所需的对随时间和空间演化的局部现象进行空间解析的广泛细尺度网格需求。与多项式近似相比，针对解定制、时间相关富集函数在获得时间解析解时实现了高达95%的自由度减少。GFEM的一个固有问题是需要先验信息来构建合适的富集函数以捕捉细尺度现象。GFEMgl通过使用局部问题的解原位计算富集函数来解决这一问题。该方法的一个关键优势是每个定义的局部问题可以独立并行求解，为持续改进分析类似问题的求解能力提供了灵活性。GFEM被视作一种高阶且灵活的数值方法，与之前可用于此类问题的标准有限元方法相比，能够在显著降低计算成本的同时提供精确解。