失速限制型风轮叶片流场的数值模拟

采用双时间步方法在旋转坐标系下求解非定常方程,构造了零偏航条件下风轮叶片流场的数值模拟方法,采用该方法对失速限制型风轮叶片NREL Phase-Ⅵ进行了模拟,其中网格采用重叠网格进行模拟,控制方程为可压缩N-S方程,同时加入了低速预处理技术,模拟的风速范围从13m/s到25m/s,这时候风轮叶片发生失速以限制扭矩输出.本文采用数值方法对风轮叶片在发生失速以后的主要流场特点进行了分析,从风轮扭矩、载荷分布等方面对数值计算结果进行了对比验证.     

不同马赫数的无粘和粘性流动高阶精度隐式计算方法

 应用隐式时间推进法对不同马赫数的无粘和粘性流动进行数值分析,给出了基于预处理方法的高阶精度隐式求解方法。方程离散采用改进的高阶精度对流迎风分裂格式。该方法通过求解曲线坐标系可压缩Euler和 Navier-Stokes方程,对低速到超音速范围内的无粘和粘性流动的典型问题进行了数值分析。计算结果与文献的数值结果或实验数据对比分析表明,该方法对低速到超音速范围内的无粘和粘性流动问题进行数值分析是可行而有效的。     

预处理AUSM+格式在全速域反应流模拟中的应用

对预处理修正的AUSM+(advection upwind splitting method+)格式进行了详细讨论, 并将其应用于全速域范围内的化学反应流场数值模拟.基于多组分Navier-Stokes方程, 结合控制方程时间导数项的预处理, 采用预处理修正的AUSM+格式发展了一种模拟全速域范围内化学反应流场的有效算法.验证算例包括超声速、亚声速以及极低速情形的化学反应流场, 计算结果表明基于预处理AUSM+格式发展的算法能够有效地用于全速域范围内化学反应流动问题的模拟.      

零质量射流推迟翼型失速的数值模拟

 通过在具有分离流动的翼型表面局部地引入零质量射流干扰,可以较小的能量代价实现对翼型绕流的全局性主动控制,从而达到抑制分离、推迟翼型失速的目的。本文发展了零质量射流与翼型绕流干扰的数值分析方法。在格心有限体积法的基础上,空间离散采用AUSM+up格式,时间推进运用含牛顿型LU-SGS子迭代的全隐式双时间法,且引入了预处理方法和多重网格法等技术。以TAU0015翼型失速主动控制为算例,通过在翼型上表面前缘附近引入零质量射流控制,使翼型失速特性得到改善。将计算结果与实验结果做了比较,验证了本文发展的数值模拟方法及程序的正确性。还以NACA0015翼型为例研究了射流动量系数、射流频率、射流偏角等关键参数对失速控制效果的影响,得出一些关于如何选择控制参数的结论。     

导弹垂直发射系统三维超声速射流的数值模拟

对导弹垂直发射系统(VLS,Vertical launching system)内的三维湍流流场进行了数值模拟.采用了单元中心有限体积法和三阶迎风TVD差分格式离散Navier-Stokes方程,采用双时间步和预处理过程提高计算收敛速度,引入对角化的双时间ADI算法求解控制方程.仿真结果准确地模拟了发射系统内的激波形成过程,真实地反映出VLS中的复杂流动.研究表明,压力室底部射流中心的波后压强以大约1 200Hz的频率振荡.      

预处理法求解定常/非定常混合网格的全速流场

运用时间导数预处理法在混合网格上求解三维定常/非定常N-S方程。预处理后的N-S方程,用二阶迎风格式的有限体积法离散,隐式、LU-SGS迭代求解,双时间步推进求解非定常流。本文以一维预处理欧拉方程为基础进行分析,提出了一种特征变量的边界条件以适应预处理后的特征系统。从极低速到跨声速范围的众多算例表明,本文的边界条件处理稳定有效,预处理法在各种速度下都具备同样快速的收敛性,在较宽的使用范围内均能得到理想的计算结果。     

预处理AUSMDV格式在微型扑翼绕流数值模拟中的应用

由于微型扑翼的低速低雷诺特性,对其粘性绕流的数值模拟存在较大困难。传统的中心空间离散格式不具有格式耗散而需要加入带有自由度的人工耗散,给计算带来不确定因素;此外,在低速低雷诺条件下,RANS方程组的条件数太大,导致计算收敛缓慢甚至不收敛,需要引入预处理方法加速收敛。本文使用具有低耗散特性的AUSMDV空间离散格式结合预处理方法对扑翼的非定常粘性绕流进行了数值模拟,非定常推进采用了含双时间迭代的LU-SGS隐式时间推进,湍流模型采用了BL模型。计算结果表明,加入预处理能够克服AUSMDV格式难以收敛的缺点,并且与中心格式相比,基于预处理的AUSMDV格式能够更准确地模拟扑翼的非定常粘性绕流,通过将本文计算结果与文献结果对比验证了本文方法的有效性。在此基础上研究了扑翼相关参数对气动特性的影响,对进一步了解扑翼气动力产生机理,指导扑翼的气动设计具有参考价值。     

低速叶型气动反问题设计方法

 低马赫数不可压流动中声速与流速大小差别巨大,采用基于可压缩流动控制方程的计算格式求解流场时,由于数值黏性的污染,解的精度低且收敛性差,通常可使用时间预处理技术来解决这一问题。在基于控制理论的优化方法中,共轭方程的Jacobian矩阵和流动方程的系数矩阵相似,因而在低流动马赫数下,求解共轭方程存在着与求解流动方程相同的数值污染和数值刚性问题。首先推导了带有预处理的Roe格式,然后发展了适合全速度流动的共轭方程求解方法,最后选取翼型和叶栅两个典型算例进行了验证。计算结果表明所发展的方法可很好地用于低马赫数时的气动反问题设计。     

基于预处理的非稳态流场计算方法

采用预处理方法可提高Navier-Stokes方程或Eu ler方程计算低马赫数流场的准确度和收敛性能。采用双时间步和预处理相结合的方法计算非稳态流场,保证了采用预处理方法后计算结果的时间精确性。在每一物理时间步,控制方程虚拟时间导数项采用LU隐式时间格式,对流通量采用AUSM 和MUSCL格式进行求解。通过对前台阶流和圆柱绕流的计算,验证了该方法的广泛适用性和计算结果的可靠性。     

预处理方法在叶轮机械三维数值模拟中的应用

应用预处理方法求解三维可压缩流动Navier-Stokes方程, 发展了一种适用于叶轮机械的既可以求解跨声速流动也可以求解不可压流动的数值模拟方法, 对NASA Rotor67和北航低速大尺寸压气机算例进行了计算, 计算结果与实验数据相吻合, 表明本文采用的预处理方法具有较大的工程实用性.