理想流体力学大作业(hesssmith方法求附加质量)

《理想流体力学大作业(hesssmith方法求附加质量)》由会员上传分享，免费在线阅读，更多相关内容在学术论文-天天文库。

1、理想流体大作业应用hess—smith方法求附加质量一物理背景当船舶在海洋中运动或潜艇在深海中运动吋，可以看成是物体在无界流中的运动，其所受到的力和力矩可以川物体的速度，加速度和附加质量决定，因此，附加质量是物体最重要的水动力特性。对于附加质量的求解就显得非常重要。巾于附加惯性力的作用，物体在理想流体屮的变速运动相当于物体自身质量上增加了一个附加质量而在真空屮运动，换句话说，理想流体增大了物体的惯性，是物体很难加速也难减速。附加质量有物体上的单位运动绝对速度势来决定，与流体密度，物体形状和运动方向有关，而与时间无关。二理论

2、依据物体的附加质砧可巾下时表示dnni=P氏，柴dsQ，j=1>2—6)<2.1)(2.2)其中，叭，(py•分别是总速度势的分量，ny表示物法线方向在各个方向分量，^7~表示物体在i方14运动引起其在j方叫的附加质:U:，，那么，只要求出在物体表而速度势在法线方向的各个分量及法向量的分量，物体的附加质量便可求出。在无界流中，物体的总速度势可以分解为定常速度0是和扰动速度势cp。其在物面上满足定解条件盖=—(2.3)应用分布源来模拟物面上的扰动势，则有

3、分量表达式勞=2崎)+<7⑷是(七)dsq(2.5)由物面条件得到2tt(j(p)+ffs_£(士)dSg=-v^n(2.6)所以只要求fli分布源密度，就可以求fli物体表而的速度势。为了计算附加质量，要求中.位速度势叭在物面上的值，其满足式(2.6)。从而得到分布源的单位速度势密度满足的积分方程。只要求解出单位速度势密度，便可以求得附加质量。三数值模型主要通过在物而上的分布源来模拟物体在无界流中所产生的扰动势，应用数值离散的方法列代数方程组求解积分方程，即用离散量来表示连续变量的方法。首先，将物体表面离散，划分成由n个

4、四边形单元近似组成。S=Z△勺7=1物而上各个单元的法向可由节点求出。其次，可以列出求解0；•的线性代数方程组其屮a'Pi从而，速度势密度a,•可解，单位速度势可求，附加质量可以计算得到叫几何模型对球体，椭球和圆柱进行网格的划分，应用不同数量的㈧边形网格来近似物体的表面，通过编写程序对其进行网格的划分，并按顺时针输出个节点的坐标值，相应的程序如图五计算参数主要对半径r=10的球壳，短轴长度为10,长短轴比为5的椭球，直径为20,柱体长和截面直径比为5的圆柱进行计算，对其进行网格的划分并求其附加质量系数。六结果讨论(1)球壳

5、对球壳进行网格划分，R出不同网格数对应的计算得到的附加质量系数曲线如图7OO..5oO附加质量系数0.45IIIIIIIIII050100150200250300350400450网格数0.7m335.66550..50.O附加质量系数450400OO00O0015054•O网格数有图付以看出，随着网格数的增加，附加质量系数趋于常数，且mil和m33的值是相同的，这主要是由于球壳是高度对称的原因(2)旋转椭球对旋转椭球进行网格划分，W出不同网格数对应的汁算得到的附加质量系数曲线如图附加质量系数0.0850.080.0750

6、.070.0650.06m335151c.11..0n11c附加质量系数0.90100200300400网格数500(3)柱体0.134对柱体进行网格划分，画出不同网格数对应的计算得到的附加质以系数曲线如陶mil附加质量系数0.1330.13250.1320.131100200300400500600网格数数nni为0,所以它们的附加质:W:系数nni逐渐减小。