基于N-S方程和自由转捩预测耦合求解的钝后缘翼型气动性能计算

大厚度钝后缘翼型由于结构和气动性能上的优点,被用于作为大型风力机叶片设计的内侧翼型。而由于其大厚度的特点和风洞实验阻塞度的限制,大雷诺数的风洞实验数据很少,给此类翼型的设计和使用带来困难。本文基于雷诺平均Navier-Stokes(RANS)方程和自由转捩预测耦合求解方法,进行了大厚度钝后缘翼型的气动性能计算研究;针对由于钝后缘后的涡脱落而造成的翼型表面压力分布的脉动,发展了有时均效应的转捩模型以考虑这种周期脉动的时均效应,并对传统的耦合求解方法进行了修改。通过对风力机翼型DU97-W-300的钝后缘改形DU97-Flat翼型的气动性能计算,分析了网格数对计算结果和计算效率的影响。并将计算的DU97-W-300翼型的气动性能和实验结果进行了比较,表明通过本文发展的耦合求解方法可以在更少的计算网格时得到比参考文献更吻合风洞实验结果的气动性能,为此类翼型的设计和使用提供数值计算基础。     

基于黏性伴随方法的跨声速机翼气动优化设计

 进行了基于黏性伴随方法和Navier Stokes方程的跨声速机翼气动优化设计研究。分别推导了适用于三维跨声速机翼气动反设计和减阻设计的黏性伴随方程、边界条件和梯度求解表达式,并研究了伴随方程的数值求解方法。通过将网格生成、流场计算、黏性伴随方程数值求解、梯度求解和拟牛顿优化算法等几方面的有效结合,发展了一种跨声速机翼气动优化设计方法。为了提高计算效率,将多重网格方法应用到方程的数值求解中来加速收敛。跨声速机翼反设计和减阻设计算例验证了本文所发展的方法的正确性。采用本文的方法进行优化设计,一般通过20~30次迭代就能得到满意的结果。     

RANS方程和附面层方程耦合求解转捩位置的方法

研究了基于线性稳定性分析的RANS（Reynolds-Averaged Navier-Stokes）方程和附面层方程耦合求解转捩位置的方法。RANS方程计算得附面层外边界的速度分布作为附面层方程计算的输入,求解附面层方程得到稳定性分析需要的附面层参数。由于附面层的解在分离点处的奇异性,用交互式附面层方程求解方法。用基于线性稳定性分析的e^N方法计算转捩位置,并考虑了Tollmien-Schliching波和层流附面层分离造成的转捩。迭代求解RANS方程和附面层方程,直到转捩位置收敛。RANS求解中使用转捩过渡区模型,避免了点转捩模型引入的数值扰动。通过对NA-CA0012翼型的计算并和实验结果及Xfoil计算的结果进行比较,吻合较好,有较好的工程实用价值。     

多目标自然层流翼型反设计方法

 进行了基于扰动放大N-因子的目标压力分布设计方法的多目标自然层流（NLF）翼型反设计方法研究。流场分析和转捩位置计算用XFOIL程序,大大减少了计算花费。用N-因子设计方法进行有NLF范围要求和满足气动约束的目标压力计算,压力恢复段的压力分布用Stratford分离准则来进行设计。用基于响应面方法的优化方法来进行反设计计算,使用不含二阶交叉项的二次多项式模型的响应面模型,大大减少了构造模型所需的试验次数;设计空间内试验点的选取满足D-优化准则。根据设计目标的设计状态,进行了多目标翼型反设计。计算结果表明,设计结果的层流范围和设计目标基本吻合,该方法可以用于NLF翼型的多目标反设计中。