涡核对自由涡计算的影响的研究.pdf

第二十届（2004）全国直升机年会论文 涡核对自由涡计算的影响研究 曾洪江 胡继忠 （北京航空航天大学直升机研究设计所，北京，100083） 摘要 在使用自由涡理论进行旋翼气动计算时，为了避免出现诱导速度无穷大的数值问题，引进 了涡核。在各种涡核模型中，都涉及到一个涡核大小选取的问题。在以前的研究中，涡核的选择 一般都是依靠经验或试验结果。本文比较系统地研究了涡核半径的选取对自由涡计算的影响。通 过把环量沿展向分布的计算结果与试验结果进行比较的方法对涡核半径进行优化，从而得出最优 涡核半径。通过比较使用优化和非优化涡核半径的计算结果，发现涡核选取对自由涡计算影响较 大。 1. 前言 在用自由涡理论进行旋翼气动计算时，由于涡丝可以在空间自由移动，涡丝与涡丝 之间就可能相交、缠绕。如果不对涡丝加入涡核修正，则基于 Biot-Savart 的诱导速度 计算可能出现数值问题。目前的涡核模型中既有选用固定涡核半径的，也有随时间变化 的涡核模型的[1], 到目前为止，不管是固定还是变化的涡核半径的选取，都是依靠经验 或实验观察结果。文献[2]中指出，桨尖涡、附着涡和尾随涡(不包括桨尖涡的内部涡， 下同) 的涡核半径应该取不同的值；在文献[1]中指出粘性涡核半径是时间的函数，在一 定的生命周期内近似与时间的平方根成正比；文献[3]中指出，在桨叶旋转过几度以后 桨尖涡就几乎完全卷起；在文献[4]中指出，相对于桨尖涡和尾随涡的影响而言，脱体 涡的影响是比较小的。根据上述研究结果，在本文中作如下假设：桨尖涡和尾随涡的环 量保持不变；桨尖涡和尾随涡的涡核半径与时间的平方根成正比；同时参考文献[2]， 为桨尖涡和尾随涡涡核半径指定一个初值，桨叶附着涡的涡核半径为常值。本文根据文 献[7]中提供的实验数据，对附着涡涡核半径、桨尖涡和尾随涡的初始涡核半径进行了 优化，从而得出针对该模型的最优涡核半径。由于自由涡的计算本身需要大量的计算机 时间资源，为获得优化的涡核进行的优化计算更会成倍增加，一种快速有效的自由涡计 算方法在这种优化计算中是非常必要的。本文首先发展了一种快速的自由涡计算方法， 在此基础上对涡核进行优化计算，进而把使用优化和非优化涡核的计算结果和实验数据 比较，得到涡核选取对自由涡的计算结果有较大影响这一结论。 2. 自由涡计算方法 本文计算方法的特点是，在文献[5]中提出的伪隐式方法的基础上，直接去除远尾 迹中引起尾迹发散的速度分量，并且在迭代过程中引入长周期概念，让新的环量和速度 场信息迅速传遍整个尾迹，从而快速得到收敛的尾迹结构。 553 (1) 迭代方法 在自由涡计算中，自由涡圈数取多少没有一个通用的准则，文献[6]中指出，尾迹 的轴向长度为桨盘直径的2 倍是合适的，由此可以大致确定自由涡的圈数。 迭代的思路是利用上一次求出的尾迹控制点诱导速度求出新的尾迹几何结构，在新 的尾迹下，根据桨叶上的控制点处满足边界条件求出新的环量，再求出新的诱导速度场， 最终得到尾迹和环量的收敛解。在通常的直接时间步进法中，每一时间步都进行上述迭 代过程，这样导致环量和整个速度场变化太快，容易引起数值发散。在文献[5]中引入 伪隐式方法，让环量和速度场在一个周期内（ψ从 0 变化到 2π，这样的周期在本文中称 为短周期）变化一次，从而增加了数值稳定性。在确定新的尾迹空间位置时，本文引入 了一个长周期的概念，一个长周期中ψ从 0 变化到 2π×n （n 为自由涡圈数）。每个长周 期内尾迹环量和诱导速度场保持不变，在每一个短周期内对尾迹位置进行迭代更新，连 续更新 n 次。引入长周期的目的是为了让每次更新以后的环量和速度场信息充分传递到 整个尾迹，因为环量是影响速度场的重要因素，速度场又决定尾迹控制点的位置，尾迹 上新生成的控制点位置又受到旧的控制点位置的影响。如果能够把这种影响充分传递到 整个尾迹，是有利于尾迹的稳定和收敛的。在每个长周期迭代结束以后，利用松弛因子 对新的尾迹进行修正，计算表明，这种修正是必要的，否则容易引起振荡，难以得到收 敛解。松弛因子修正公式为： n n−1 n r (ψ ,ζ ) (1−ω )r