基于多尺度无网格galerkin方法的流动问题模拟-应用数学专业论文.docx

摘要无网格 Galerkin 方法(EFG)是近年来迅速兴起的一种数值方法。该方法构造 形函数时只需要具体的节点信息，而不要网格，因此可显著减少因网格畸变带来 的困难。在涉及大变形、自由面、局部区域高梯度等问题的研究中，EFG 方法 显示出具有发展前景的明显优势。目前，EFG 方法在计算流体中的应用研究还比较少。一般在流动问题的求解 中，EFG 方法会遇到对流占优、速度-压力失耦、自由面难以追踪等问题。同时，在微尺度/纳米尺度下，EFG 方法也不能模拟出流体运动的微观效应。针对上述 问题，本论文为不可压流动问题的求解发展了相应的 EFG 方法。其主要工作和 结论如下：(1)针对 EFG 方法研究中形函数只有矩形、圆形影响域的研究现状，将其形 函数拓展到了任意形状的影响域。数值实验表明：当影响因子趋于 1 时，EFG 形函数具有类似于有限元形函数的单元分片性质，且在节点处近似满足 Delta 函数属性。(2)针对 EFG 方法求解对流占优对流-扩散方程时出现的数值伪振荡现象，提 出了变分多尺度无网格 Galerkin(VMEFG)方法，即全局 VMEFG(GVMEFG)方法和局部 VMEFG(LVMEFG)方法，后者又包含 LVMEFG 数值(LVMEFG_N)方法和 LVMEFG 分析(LVMEFG_A)方法。数值实验表明：这些 VMEFG 方法均可很好 地消除由对流项占优引起的数值伪振荡，但 GVMEFG 方法的计算量极大。(3)针对 EFG 方法求解不可压流动问题时产生的速度-压力失耦现象，提出了 LVMEFG 方法，即 LVMEFG_N 方法和 LVMEFG_A 方法。数值实验表明：这些 LVMEFG 方法在速度-压力采用等低阶基函数近似时，能有效地解决速度-压力的失耦问题。然而，对于不可压流动问题，LVMEFG_N 方法的计算量很大。(4)针对温度场在 EFG 方法求解时因对流项占优引起的伪振荡现象，将求解 等温不可压流动问题的 LVMEFG_A 方法拓展到了非等温不可压流动问题的求解。数值实验表明：LVMEFG_A 方法确实能消除温度场在 EFG 方法求解时因对 流项占优引起的伪振荡现象。(5)为了准确地模拟自由面问题，将 LVMEFG_A 方法与自由面捕捉技术相结合。数值实验表明：基于 LVMEFG_A 方法和自由面捕捉技术可方便准确地追踪 水波晃动现象中的自由面位置。(6)为了反映微尺度/纳米尺度下的物理现象，针对等温/非等温不可压流动问题，将分子动力学模拟与 LVMEFG_A 方法基于“区域分解”技术相结合，建立了一种耦合分子动力学的“真正的”多尺度 EFG 方法。数值实验表明：在流体力学领域中，基于“区域分解”技术耦合分子动力学模拟与 LVMEFG_A 方法， 能保证两者在重合区域耦合很好，数值解几乎重合。关键词：无网格 Galerkin 方法，变分多尺度方法，分子动力学，不可压流动， 自由面，多尺度方法ABSTRACTThe element free Galerkin (EFG) method is a new kind of numerical method, and it has been developed rapidly in recent years. When the EFG shape function is constructed, only node is required instead of mesh, so the influence of the mesh distortion is reduced evidently. Thus, this kind of method has some advantages with wide development prospects over the simulation of the problems which refer to large deformation, free surface and so on.Currently, EFG method is less involved in the field of fluid mechanics. When the flow problem is simulated, EFG method will suffer from many difficulties such as the advection-dominated, the coupling of velocity and pressure and the tracking of free surface and so on. Additionally, for the microscale or nanoscale f