基于CAESES的超临界二氧化碳（sCO2）轴流透平叶片优化设计研究

这篇文章中采用CAESES进行sCO2轴流透平设计的方法，发表于2019年7月17日至21日的ASME透平会议。

 简 介 

传统发电厂采用蒸汽作为工质，通过透平产生动力，超临界二氧化碳（以下简称sCO2）循环使用的是温度和压力均高于临界点（超临界状态）的CO2，在这种状态下，CO2表现出介于气体和液体之间的特性，并且具有较高的密度和体积热容，这种状态下的特性为高循环效率提供了巨大潜力。由于工质的能量密度更高，因此可以减小组件尺寸，从而减小占地面积和成本。sCO2也被认为是一种安全的介质，其资源十分充足且使用收益高，因此，从效率和成本角度来看，sCO2发电有潜力取代蒸汽发电。

本文对用于废热回收应用的新型sCO2轴流透平设计进行了探索，文中基于10兆瓦的案例进行介绍。

文中采用Kulfan Class Shape Transformation（CST）变换方法进行二维轴流叶型轮廓变形优化设计，并在设计优化过程中同时考虑叶片的气动效率及应力情况。

 轴流透平设计原理 

首先基于尺寸、性能、运行工况等设计需求，项目中使用了Triveni Turbines开发的一维均线计算内部工具进行设计计算，均线计算的结果构成了二维叶片轮廓设计的基础。

项目中采用CAESES进行轴流透平的几何建模，调用二维/准三维流动求解器MISES用于方针分析，并采用印度科技学院（IISc）内部开发的Matlab脚本进行前后处理。通过CAESES软件的自动优化平台封装了整个过程，用以优化透平叶片的气动性能。

有多种方法能够用于轴流透平叶片二维截面形状的参数化建模，下图展示了一种通用方法，该方法可以直接控制有意义的参数，例如气流角、楔角、前尾缘半径、厚度等，但是这种方法在几何变化的灵活性方面存在局限性，因此很难得到突破性的进展。其他方法比如采用B样条曲线建模和采用类函数变换（CST）方法等，可以带来了更大的灵活性和局部形状控制，从而使得优化过程中出现新的、出乎意料的设计。

透平叶片截面形状参数化定义

本研究使用基于B样条曲线和CST方法进行几何参数化、仿真和优化，具体步骤如下：

1. 导入现有二维基础轮廓线；

2. 将其分割为压力面和吸力面；

3. 使用B样条曲线和CST方法自动匹出新的压力和吸力面曲线；

4. 创建设计变量控制B样条曲线控制点和CST参数变化，从而控制叶型轮廓线变化；

5. 创建几何约束，例如某些区域的内切轮廓面积和的曲率单调性，保证曲线变化时的光顺性；

6. 建立软件连接，导出文件，运行Matlab脚本及MISES软件；

7. 在CAESES中使用DAKOTA全局优化算法，进行透平叶片气动性能优化。

透平叶片的表面曲率对亚音速机翼的压力分布有重大影响，因此CAESES中创建了独特的曲率识别功能，并将其作为约束，从而在保证参数变化灵活性的同时，保持曲率分布合理过渡。下图显示了拟合后叶片型线轮廓的曲率分布，其中位置0表示压力面的尾缘，位置1表示吸力面的尾缘。

叶型表面单调曲率分布

通过曲率分布的导数能够评估曲率的单调性，将曲率分析划分为不同的区域，并设置约束条件来判断这些区域的单调性，下图显示了违反曲率约束的示例（以红色显示）。CAESES中能够采用光顺程序进行曲线曲率分布的自动光顺（仅适用于B样条曲线）。

 轴流透平叶片的空气动力学优化