12高精度旋翼非定常涡流场数值模拟方法.docVIP

第二十届（201）全国直升机年会论文 史勇杰 魏 鹏 徐国华 招启军 （南京航空航天大学直升机旋翼动力学重点实验室，南京，210016） 摘要：Caradonna-Tung模型旋翼悬停状态下的旋翼流场进行了模拟，通过分析旋翼尾迹特征及耦合边界处流场特征，表明所建方法充分发挥了CFD方法和粘性涡方法各自的优势，在旋翼流场计算方面的具有独特的优势。最后通过对比桨叶表面压强系数分布，进一步验证了该方法在旋翼流场分析中的有效性。 关键词： 0 引言 图1给出了本文耦合拉格朗日-欧拉法的示意图，该机算方法将 图 1 耦合拉格朗日-欧拉法示意图 1.1 旋翼尾迹求解模型 直升机飞行时，旋翼桨叶拖出的尾涡发生卷起现象并形成桨尖涡，这些结构紧凑且强度很高的桨尖涡堆积在桨盘下（后）方，成为旋翼尾迹的主要组成部分。拉格朗日粘性涡方法将该尾涡系离散为一系列无序涡元，并在拉格朗日描述体系进行数值求解来预测尾迹的卷起、畸变及耗散等过程。不考虑桨叶附近流场，直升机流场问题可认为不可压问题，离散涡元随时间的变化由运动学方程和涡量动力学方程进行求解 式中，下标为涡元索引号，表示位置，为环量值，为当地速度,为涡诱导速度，为时间，为运动粘性系数。 (1)中，，Biot-Savart定理计算，并引入快速多极子加速算法(FMM, Fast Multipole Method)[10]对该计算过程进行加速。而第二项表示的是涡量的粘性扩散效应，本方法采用粒子强度交换法(PSE, Particle Strength Exchange)[11]对其进行求解。 受桨尖涡诱导作用，旋翼流场存在着的强烈剪切流导致长时间数值计算后离散涡元会出现严重的疏密失调，导致涡元间无法满足保证计算精度的重叠性要求，从而致使计算精度下降。为此，本文采用文献[12]的涡元重布方法来调整涡元的疏密度，该方法根据当地速度应变张量判断局部流场的拉伸状况，当涡元所在位置的局部拉伸达到特定值时，将涡元分裂成两个相同强度的涡元填补由于拉伸导致的“空白”。然而，随着新涡元从桨叶脱落以及旧涡元的分裂，计算过程中涡元数目会迅速增长，导致计算效率降低；为合理控制因涡元数目增长导致的效率问题，一方面可以将远离旋翼的涡元删除，另一方面则需要对距离很近且方向较为一致（涡矢量夹角小于特定值）的涡元进行合并处理。 为捕捉桨叶近体空间的非定常流动 式中，守恒变量，为惯性系下三个方向上的速度，分别表示压强、密度和能量；为控制体单元体积，为面积，分别为对流通量喝粘性通量，为粘性项逻辑开关。 本文旨在研究所建方法的耦合有效性，着重于两计算域在交互过程中的信息传递，考虑到粘性项对此没有直接影响，在实际计算过程中一方面出于方便对数值耗散进行研究的缘故，另一方面出于效率方面的原因，并未打开粘性项逻辑开关。计算过程中所有的物理量均以悬停桨尖速度、密度和桨叶弦长作为参考值而无量纲化。控制方程的空间离散采用二阶中心方法，而时间推进上采用双时间方法。 对于本文所建方法，拉格朗日计算域与欧拉计算域之间的相互信息交换至关重要，合理的信息交换策略无论对计算精度还是对计算效率都有帮助，因此需要对其进行仔细处理。 2.1 欧拉计算域至拉格朗日计算域 CFD计算域至粘性涡计算域的信息传递，可采用两种方式：集中涡源法和分布涡方法。后者直接将CFD计算域计算的涡量分布传递给拉格朗日计算域，物理意义更为明确，但却容易将CFD计算域的数值耗散引入到拉格朗日计算域，同时也存在着计算工作量大的问题。为此，本文采用了集中涡源法，即先由CFD方法计算得到桨叶展向载荷分布，再根据库塔—如科夫斯基（Kutta-Joukowski）环量/升力定理来获得桨叶展向附着涡环量分布 其中是桨叶剖面升力系数，为附着涡环量值。 桨叶附着涡随空间及时间的变化是桨叶脱落尾涡的主要原因，其中附着涡随空间变化脱落的尾涡为尾随涡，随时间变化脱落的桨尖涡为脱体涡。因此在求解得到桨叶附着涡后，拉格朗日计算域新生涡元的强度即可由下面公式确定 式中，表示新生涡环量值，为桨叶剖面相对气流速度，包括来流速度、旋翼旋转速度以及桨叶挥舞运动速度。 2.2 拉格朗日计算域至欧拉计算域 粘性涡计算域至CFD计算域的信息传递，常用的有几种方式，第一种方法是采用涡诱导速度直接修正桨叶当地迎角；第二种是采用广义网格速度法将影响施加在全体网格格点中，即将涡线产生的诱导速度等效成网格畸变；第三种方式直接通过边界条件施加于CFD计算域外边界上。其中，第一种方式尽管效率很高，然而在迎角修正过程中会引入误差，精度并不高；第二种方法尾迹贡