小振幅振动翼型升力迟滞环变向的成因分析

在低雷诺数条件下研究了SD7037翼型小振幅强迫振动时，振动轴位置和减缩频率对翼型气动特性的影响。通过升力系数随时间变化曲线与迎角随时间变化曲线的对比，发现引起升力系数随迎角变化的迟滞环变向的原因主要是升力系数随时间变化曲线的相位发生变化，并从数学上证明了这个相位是关于减缩频率的函数。而且，当升力系数随时间变化曲线的相位变大时，会使得对应的升力系数随迎角变化的迟滞环曲线从逆时针向顺时针方向变化，这中间必然会经历一个"直线"过程，就好像迟滞现象消失了一样。进一步通过对比流场，发现振动轴位置和减缩频率对此相位的影响机制不同，振动轴位置主要改变翼型的有效迎角，而减缩频率的增大则影响了周围流场结构，使得附加质量带来的反作用力变大，进而提高升力，振动轴向前移动和减缩频率的增大都会增大升力系数曲线的相位。     

一种适合迭代求解的反馈力浸入边界法

对Goldstein提出的反馈力浸入边界法进行了新的思考，改进了其对力源项的计算，拓展了该浸入边界法的使用范围。传统的反馈力浸入边界法在进行力源项的计算时，含有对速度误差的时间积分项，只能用于含时间项的Navier-Stokes （N-S）方程的求解，且在显式时间推进时有严格的时间步长限制。本文改进的方法则直接通过迭代过程中的速度误差求和来计算力源项，避免了时间相关的参数，使其不仅能适合非定常隐式时间推进，还能与不含时间项的定常N-S方程求解方法结合。为了验证改进方法的可靠性，对二维静止圆柱绕流、静止流体中的振荡圆柱、运动椭圆翼以及三维静止圆球的流场进行了计算，计算结果均与文献结果符合较好，表明本文改进的方法是有效的。得出的结论为：可以直接基于迭代次数进行反馈力源项的计算，改进的反馈力浸入边界法不仅可与非定常N-S方程结合，进行隐式求解，还可以与定常N-S方程结合用于定常流动的模拟，可将改进的方法运用到更多的流动问题当中。     

螺旋桨/机翼耦合下的目标螺旋桨滑流设计

针对分布式电推进(DEP)构型等具有多螺旋桨特征的飞行器,发展了通过优化螺旋桨滑流来达到提高机翼升阻比的方法。提出了一种可以获得目标诱导速度分布的螺旋桨设计方法,基于面元法发展了一套可以快速计算螺旋桨机翼干扰的气动程序Prop-wing,基于Kriging代理模型建立了一套高效的优化方法获得最优的螺旋桨诱导速度分布提高机翼升阻比。优化结果显示当拉力保持相同时,螺旋桨桨毂附近的轴向诱导速度越大,下游机翼的升阻比越大。在不对螺旋桨功率进行限制时,优化后的螺旋桨使得下游的翼段阻力相比较安装最小能量损失设计的螺旋桨的翼段减少了1875%,而翼段升阻比提升达到了2563%,当优化螺旋桨功率被限制后,翼段升阻比提升为962%。虽然升阻比的提升需要付出螺旋桨效率下降的代价,但是研究还是给分布式动力滑流的利用提供了一种思路。