理想流体力学课程设计(Hess Smith方法求附加质量).docVIP

一、物理背景 二、理论依据 用表示无界流中的物体表面，来流为均匀流，其未扰动速度或无穷远处的速度为 (2.1.1) 用表示定常速度势，它在物体外部空间域中适合拉普拉斯方程，在物面上适合不可进入条件，在无穷远处，应该与均匀来流的速势吻合即 （物体外） (2.1.2) （物面上） (2.1.3) （） 其中,单位法线向量指向物体内部。 在速度势中分出已知的均匀来流项，记 (2.1.4) 这里的是扰动速度势，应适合以下定解条件： （物体外） (2.1.5) （物面上） (2.1.6) （无穷远处） (2.1.7) 易知过物面的通量为零，即 所以远方条件(2.1.7)可进一步具体化为 （） (2.1.8) 用表示点和之间的距离，对函数和在物面外部和远方控制面的内部之空间域内用格林公式，当点在上述空间域内时 (2.1.9) 从的远方条件(2.1.8)可知，上积分趋于零，式(2.1.9)成为 (2.1.10) 其中, 是物面外的任意一点。 在物体的内部域中构造一个合适的内部解，它在内部适合拉普拉斯方程，在物面上适合某种物面条件，其具体形式将在下面给出。对于上述物体外部的点函数在物体内部域中没有奇点，在内部域中对函数和用格林公式，得到 (2.1.11) 式(2.1.10)和(2.1.11)中的是物体外部同一个点，把两式相减，得到 (2.1.12) 在物面上取适合下述两种物面条件，得到两种的定解条件，一种是： (2.1.13) 定解问题(2.1.13)是拉普拉斯方程的第一类边值问题，它的解是存在且唯一的。取式(2.1.12)中的内部解为式(2.1.13)所决定的函数，则式(2.1.12)成为 　　　　　　　　　 (2.1.14) 其中 　 　(2.1.15) 式(2.1.14)表示扰动速度势可以用物面上的分布源表示，其中分布源密度是未知函数，将由扰动势的物面条件（2.1.6）来决定。 物体的附加质量，表示物体沿方向运动引起的方向的附加质量，公 式如下: (2.2.1) 利用式(2.1.14)，再结合物面条件，得到 (2.3.1) 这就是分布源密度所适合的线性积分方程。 把积分方程(2.3.1)转换成线性代数方程组，即用离散量代替连续变量。把物面分成小块，记 (2.3.2) 用平面四边形或三角形来近似代替小曲面。具体做法如下，取第小块的四个顶点坐标之算术平均值，得到中心点的坐标。计算对角线向量的向量积（指向与曲面法线指向相符合），用表示该方向上的单位向量，形成以为法线且通过中心点的平面，再把四个顶点向该平面作投影，以四个投影点为顶点组成平面四边形，用代替原来的小曲面，称为单元。通常把小范围内的分布源密度作为常数，因此只要分割不太粗，可以认为在单元上为常数，记作，从而 (2.3.3) 因此物面上的积分可以用个平面四边形（三角形）上积分之和来近似，即 (2.3.4) 上式左端的未知量是连续型变量，而上式右端的未知量是个离散量。为了求解这个未知数，须要个方程。取积分方程(2.3.1)中的动点为个单元的中心点，称之为控制点，即控制物面条件使之成立的点。用近似式(2.3.4)代替积分方程(2.3.1)的左端，便可以写出的阶线性代数方程组： (2.3.5) 其中 称为影响系数，即第个单元上的分布源在第个控制点上的影响。 求解线性代数方程组(2.3.5)得到的值以后，便可以得到速度势在控制点处的值，即 (2.3.6) (2.3.7) 另外，物面上的诱导速度为 (2.3.8) 其中表示求和是不计这一项。，这里的曲面法线指向物体内部。 三、模型 ，每一个四边形面面元为。为了简化计算，将面网格投影到各自对应的平面上，使曲面网格变为平面网格。投影的方法为： 取四个顶点坐标之平均值，作为中心点的坐标。计算对角线连线向量的向量积并使得积的方向与流域法向相同。用表示