实施方案
[0032] 下面结合附图对本发明做进一步说明。
[0033] 如图1所示，基于能量梯度理论的离心压缩机弯道改进方法，具体步骤如下：
[0034] 步骤1、模拟离心压缩机模型级内的流动物理参数；
[0035] 步骤2、计算离心压缩机模型级的弯道内部流场的能量梯度函数值K，具体计算过程如下：
[0036] 根据能量梯度理论(参见论文《Mechanism of Flow Instability and Transition to Turbulence》，窦华书，International Journal of Nonlinear Mechanics期刊，2006年，5卷41期)，离心压缩机模型级内的能量梯度函数值K的计算公式为：
[0037]
[0038] 式(1)中， 为流体总压，p为流体静压；k为空气的比热容，通常取1.4；Ma为气流的马赫数；n为流体流动的法线方向，s为流体流动的流线方向；H为流体的能量损失；
[0039] 流体在离心压缩机模型级弯道内的流动属于压力驱动流动，其 和 的计算如下：
[0040]
[0041]
[0042] 其中：α表示弯道内流体在x方向的速度v与流体速度矢量之间的夹角。
[0043]
[0044]
[0045]
[0046]
[0047]
[0048]
[0049]
[0050] 表示P在x方向的一阶偏导， 表示P在y方向的一阶偏导， 表示U在x方向的一阶偏导， 表示U在y方向的一阶偏导， 表示μ在x方向的一阶偏导， 表示μ在y方向的一阶偏导， 表示α在x方向的一阶偏导， 表示α在y方向的一阶偏导， 表示U在x方向的二阶偏导， 表示U在y方向的二阶偏导， 表示U关于x、y的混合偏导。
[0051] 步骤3、根据弯道各位置能量梯度函数值K的大小，找到弯道内流动最不稳定的位置，将弯道对称设置的后半段圆弧形子午型线改进成直角导角形子午型线，直角导角形子午型线包括圆弧段及分别与圆弧段两端连接的两条直线段；重新计算弯道流动最不稳定位置的能量梯度函数值K，优化出弯道的直角导角子午线的圆弧段最佳半径。
[0052] 步骤1所述的离心压缩机模型级内的流动物理参数的具体模拟过程如下：
[0053] 采用CFD(计算流体动力学)对离心压缩机模型级内的流动进行数值模拟，模拟过程中控制方程采用定常三维可压缩的雷诺平均纳维-斯托克斯方程和连续性方程，同时使用Spalart-Allmaras湍流模型(单方程湍流模型)计算离心压缩机模型级内的流动；网格采用结构化的网格，并利用有限体积法对结构化网格下的控制方程在空间上进行离散；然后，在计算域上施加边界条件，进行模拟计算，得到离心压缩机模型级内整个流场的物理参数，包括气流速度、总温、压强和湍流粘度。所述计算域上施加的边界条件为：入口边界条件设定的是总温、总压入口，出口边界条件设定的是质量流量出口。
[0054] 步骤3所述的优化出弯道直角导角子午型线的圆弧段最佳半径，确定最佳半径的标准如下：在该半径下的弯道直角导角子午型线，使得弯道内流动最不稳定位置的能量梯度函数值K最小，稳定性能最优，则该半径为最佳半径。
[0055] 如图2所示，当弯道外壁半径为54.37mm时，计算弯道内部各位置的能量梯度函数值K，发现弯道外壁附近气流不稳定。将弯道对称设置的后半段圆弧形子午型线改进成直角导角形子午型线，分别取弯道直角导角子午型线的圆弧段半径为5mm、10mm、15mm、20mm、25mm、30mm、35mm和40mm进行数值模拟，结果表明，当直角导角子午型线的圆弧段半径取
5mm、10mm、15mm、20mm、25mm、30mm和35mm逐渐变大时，弯道内部特别是弯道外壁面附近的能量梯度函数值K逐渐减小，说明流动稳定性增加；但直角导角子午型线的圆弧段半径增加到
35mm后，再继续增加，改进后弯道的子午型线与原模型相比，内部能量梯度函数值K改进不明显，所以直角导角子午型线的圆弧段半径取35mm为最佳。如图3所示，直角导角子午型线的圆弧段半径为35mm时离心压缩机整级的多变效率较原模型有了显著提升。
