EFGM求解传热边界形状辨识和相变问题及自适应EFGM计算

【摘要】：传热边界形状辨识和相变传热问题的数值模拟具有重要的工程应用背景与理论探讨价值。由于这两个问题都涉及到边界变化的问题,因此,如何避免网格重构造成的计算困难成为相关数值模拟中所关注的一个重要问题。无网格伽辽金法(Element-free Galerkin method, EFGM)基于点的近似,不需要网格的生成和重构,具有良好的精度、稳定性和收敛性。因此,本文以EFGM为基础,探讨了传热边界形状辨识和相变问题的数值求解技术。 自适应EFGM只对计算精度不满足要求的区域进行模型改进,和均匀加密相比,可在提高精度的同时,显著减少所需的节点数量,降低计算量,具有重要的工程应用价值。因此,本文还对EFGM的自适应计算进行了研究。 主要的研究工作包括： 1.将EFG和水平集方法(Level set method, LSM)相结合,提出了一个求解传热边界形状辨识问题的数值模型。EFG求解边界变化的的温度场,不需网格重构；LSM可在固定网格下方便和准确地描述形状演化。通过一个曲线边界的识别问题对所提模型进行了数值验证,考虑了初值、EFG节点/LSM有限差分网格和测量误差的影响,计算结果令人满意。 2.将EFG与函数光滑化技术相结合,提出了一个求解相变传热问题的光滑等效热容数值模型。借助Sigmoid函数得到光滑连续的等效热容曲线,有效避免了等效热容阶跃变化引起的计算误差,不需使用传统的平均近似,对集中/非集中质量热容阵均适用。采用EFGM/有限差分法进行空间/时域离散。通过平板/角域凝固的数值算例对所提模型进行了数值验证,考虑了EFG节点、光滑函数参数的影响,并与平均近似方法对比,计算结果令人满意。 3.建立了一个基于节点误差估计的自适应EFG数值模型,并应用到弹性力学问题中。基于节点的误差估计可减小估计误差在节点间的虚假振荡。和Chung-Belytschko误差估计的自适应相比,可在引入较少节点数量的条件下取得较高的精度。 4.提出了一个偶应力正问题的h自适应EFG数值模型。EFG避免了有限元求解中因C1连续性要求造成的困难,自适应计算提高了计算精度,与均匀加密相比,显著降低了计算量。本文还建立了一个多宗量偶应力反问题的EFG和Gauss-Newton数值模型,可有效地对区域非均质的本构参数／载荷进行单一／组合反演,具有良好的抗噪性。 论文工作,有望为无网格方法的进一步研究／应用及相关问题的数值求解提供有价值的参考。