论文部分内容阅读
在实际工程应用当中,经常会遇到带有移动边界的复杂的非稳态流动问题,例如空气弹性变形,流固耦合问题,多体分离问题等。计算流体动力学中的一个关键问题就是如何处理移动的边界。计算区域随时间改变,计算区域的网格也应作出相应的调整。网格的调整方法包括网格的移动与重剖分。网格移动方法只改变节点的位置,节点之间的拓扑连接关系不变。网格的重剖分是对求解区域重新划分网格。网格的重剖分需要重新生成一次网格,并且需要将原网格上的信息映射到新的网格上,这两步的计算量较大,并且容易产生插值误差,对后续的计算不利。因此,首选采用网格的自适应移动方法。本文分析研究并改进了几种主要的网格移动方法。首先,在绪论中介绍了课题的研究背景及国内外的研究现状,并分析了几种当今流行的网格移动算法。在第二章中分析并实现了弹簧近似方法,并对ball-vertex方法进行了简单的介绍。弹簧近似法仅能处理网格的中小变形问题。针对网格大变形问题,第三章和第四章分别分析并实现了基于Laplacian方程的网格移动方法和弹性体方法。基于Laplacian方程的网格移动方法是采用有限元法求解Laplacian方程,将需要移动的网格直接作为有限元离散网格,分别求出每个单元节点三个坐标方向的位移量;将位移叠加到上一步的节点坐标上,即可得到更新后的网格。本文提出将一个体积系数引入到基于Laplacian方程的网格移动方法中,使移动物体可转动更大角度而不会出现网格的交叉重叠问题。基于Laplacian方程方法只需求解三次n阶的线性代数方程组,n为网格节点数,而弹性体法需求解的方程组为3*n阶。所以,在求解线性代数方程组的过程中,此方法有更小的计算量。与基于Laplacian方程方法和弹簧近似法相比,弹性体法的计算量大,但处理网格大变形的能力更强。弹性体法中包括弹性模量和泊松比两个参数,对单元形状的控制能力强。本文对弹性体法进行了改进,进一步提高了该方法解决网格大变形的能力。