动量源方法2资料.pptVIP

动量源方法Momentum Source Method V22悬停状态时流场 将旋转的桨叶简化成无限薄的作用盘，桨叶对气流的作用被等效为时间平均的动量源项加到Euler方程的右侧，即 S。将桨叶沿展向离散成许多微段，可以认为在此微段上翼型弦长c、翼型厚度、负扭转、安装角等参数均为常数。取其中长度为dr的微段为例进行说明，如下图所示。 * * * * * * 目前，直升机旋翼流场的数值模拟主要采用两种方法：一是涡流理论方法，二是计算流体力学方法。其中，利用涡流理论方法求解旋翼流场已经取得了很大的进展，并在实际中得到广泛应用。但涡流方法的主要不足在于它对旋翼流场细节特征的模拟不够准确。而CFD方法若采用Euler方程或Navier-Stokes方程作为控制方程，涡量可以作为解的一部分存在，因而可以很好的捕捉旋翼尾迹，从而更准确的捕捉流场信息。 一、概述 在利用CFD数值模拟旋翼流场时，主要有两种思路： （1）围绕各片桨叶生成贴体网格，整个旋翼网格系统为运动嵌套网格。在该网格中求解Euler方程或Navier-Stokes方程来模拟旋翼流场； （2）利用作用盘（Actuator Disk）理论，将旋转的桨叶等 效为一个作用盘。动量源方法属于作用盘方法的一种，它的基本思想是，桨叶对气流的作用被等效为时间平均的动量源项添加到控制方程（Euler或Navier-Stokes）中。这样，桨叶对气流的作用通过气流动量的变化来表征。 嵌套网格方法 为降低计算代价和提高模拟旋翼下洗流场的效率，一些研究者已进行了很多尝试。Rajagopalan和Chaffin采用以桨盘代替旋翼的方法对旋翼流场进行了计算，在计算效率上取得了明显的进展。采用动量源方法最大的特点是桨叶对气流的作用以动量源的形式来代表，忽略桨叶附近的细节流动，把周期性的流动通过时间平均方法转化为“准定常”流动。采用上述方法在保证旋翼下洗流场的本质属性的前提下，舍去求解旋翼桨叶周围流场的流动细节给计算精度提出的要求。同时，不用采用围绕桨叶生成贴体网格，简化了网格生成。 二、带动量源项的三维Euler方程及其离散 分别为流体的密度和压强 表示单位体积的总能 表示笛卡尔坐标系下的三个速度分量 为动量源项 为动量源项在笛卡尔坐标系下的三个分量 将计算域分成有限个互不重叠的网格单元，这样流场空间被离散化，在每一个网格单元中求解Euler方程。任取其中的一个网格单元K，在网格单元K中Euler方程可表示为： 因为网格单元很小，可以认为守恒变量W在网格单元K中 处处相等，所以有： i 表示四面体K的一个网格面，n 为面 i 的外法向矢量， 引入变量Z 可以得到： 最终可以得到空间离散后的Euler方程： 为网格单元K的动量源项，将在下一节中介绍它的求法。 为减小解的振荡，在方程的右边加入人工耗散项，可得： 时间方向上采用五步Runge-Kutta法对控制方程进行求解， 并引入了当地时间步长，以加速解的收敛。 三、动量源项的求解 桨叶微段 桨盘平面 在Euler方程的不断迭代中，可以求得任意时刻桨叶微段处气流的绝对速度（惯性坐标系中） 桨叶微段局部坐标系的旋转速度为： 可得到桨叶微段处气流的相对速度为： 由公式 可以得到桨叶微段处的当地音速 进而可以得到当地马赫数 由当地马赫数与桨叶微段处的气流迎角，可以确定翼型的 桨叶微段上的力为： 将计算的力分解到桨叶微段局部坐标系中，并忽略径向力 表示桨叶旋转一周所需时间 表示单位时间内单片桨叶旋转了多少圈 单位时间内单片桨叶作用于半径 r处，角度为 的矩形微面的力为： 单位时间内N片桨叶作用于该矩形微面的力为： 如果采用四面体网格，所以在桨盘处的微面为三角形微面，为此，将作用于矩形微面的力矢量按面积平均，再分配到三角形微面中，可得： 图2 桨盘平面 由力的牛顿第三定律可知，桨叶微段作用于气流的作用力为： 桨叶微段受力为： 即“动量源项” 将 分解到计算域坐标系中，可得动量源项 在不断迭代求解Euler方程的过程中，动量源项也在不断的变 化，求得的动量源项添加到桨盘下方紧邻桨盘的网格单元中 上述动量源项的求解过程可用下面的图形表示： 四、一些采用动量源方法计算的结果 涵道间隙分别为0.01R，0.03R时的流场轴向速度分布