跨音松弛法在叶栅绕流中的应用

文章介绍跨音松弛法应用于叶栅绕流的定常与非定常流场计算。在定常流中,采用弦线法向小扰动假定;在非定常流中,采用线化非定常分量小扰动假定,从而简化了数值求解过程。由于松弛法大大地节省机时和内存,便于工程上推广应用。     

火焰筒头部流场研究之二——二元轴对称等温廻流区计算

本文将火焰筒头部的等温流场,用二元轴对称不可压定常流来处理并用数值方法解出初步结果。基本方法是用涡量流函数法,双方程粘度模型,采用交替方向隐含格式差分求解,在边界及迭代初值提出自己的看法,初步计算结果与[10]之实验定性一致。     

低速流预处理Roe格式中的数值粘性

使用Roe格式计算低马赫数流动,需采用预处理技术.在预处理技术中,需要选择适合的预处理参数,但有时很难找到,因而得不到合理的计算结果.通过数值实验,发现使用预处理方法时,格式的数值粘性对低速流动的模拟有显著的影响.提出了在预处理Roe格式中加入可调参数,用以控制数值粘性对计算结果的影响.带有可调节数值粘性参数的预处理Roe格式的数值粘性更接近低速流动的真实耗散,反映了低速流动的真实物理现象.对方腔流动和圆柱绕流的低雷诺数定常流动所做的数值实验,得到了比较满意的结果,表明这一改进是有效的.对同一低雷诺数不同低马赫数的定常圆柱绕流所做的数值实验,均得到了满意的结果,表明了方法的适用性.     

浮空气囊定常及非定常气动特性计算

对于0~20 km范围内不同外形的浮空气囊,采用虚拟压缩方法求解不可压缩非定常N-S(Navier-Stokes)方程,数值模拟了5个高度下气囊的定常和非定常流场.定常状态时,随高度及气囊外形的变化,模拟了气囊气动力及背风区复杂流动分离.结果表明,当气囊厚度率逐渐增大时,其背风区存在复杂的空间流动分离.运用双时间步方法及动网格技术,成功实现了气囊平移或俯仰运动状态下非定常流场的数值模拟.考察了不同的平移方向及速度对气囊非定常气动特性的影响,得到了气囊平移过程中气动力及流场的变化情况.对于气囊的俯仰运动,获得了不同俯仰速度及不同最大俯仰角时气囊的气动力及流场.计算结果表明气囊的非定常运动对其流场和气动力影响显著.     

用GAO-YONG湍流模式对翼型跨音粘流的数值模拟

应用GAO-YONG可压缩湍流模式,数值模拟了NACA0012,RAE2822翼型的定常跨音速粘流算例.对流项采用三阶ROE格式,扩散项采用二阶中心格式,用多步Runge-Kutta显式时间推进法求解空间离散后的控制方程.计算结果很好地预测了翼型表面的压力系数的分布、激波的位置、马赫数等值线的分布等情况,并且对翼型表面激波与边界层相互干扰以及层流向湍流的转捩问题进行了分析计算.计算结果与实验值符合很好,表明GAO-YONG可压缩湍流模式应用合适的计算方法能够高精度模拟翼型跨音粘性流动问题,并且基于GAO-YONG可压缩湍流模式各向异性湍流粘性的机理,提供了一种预测转捩起始位置的判别方法.      

关于非结构网格生成中的背景函数

背景信息对于用阵面推进法生成非结构网格是一个关键的环节.文章基于众多文献中提到的背景信息提供方法,提出一种更为实用的背景信息提供技术.此方法不仅避免了不必要的机时,内存的消耗,又使网格具有光滑性,网格的生成具有灵活性,方便性的特点.基于生成的网格采用有限体积法研究无粘高超音速化学非平衡流的流场分布,生成的网格能够成功的运用于高超音速化学非平衡流的数值计算。     

固体火箭发动机燃烧室三维流动数值计算

采用k-ε双方程湍流模型,以SIMPLE计算程式求解翼柱型装药固体火箭发动机燃烧室内三维非定常不可压流N-S方程.在建立翼柱型装药简化模型的基础上采用边界标志法来表达燃面推移,采用从二维到三维的初场给定方法结合多重网格法求出了多个时间步的非稳态流场结构.计算结果表明燃烧室内旋涡运动呈现一定的空间与时间分布,周向加质导致燃气通道横断面上轴向速度分布的极不均衡.      

间隙非线性气动弹性颤振控制

研究了亚音速不可压来流中二元机翼气动弹性颤振主动控制问题.采用Theodorsen准定常气动力建立俯仰方向含有间隙非线性的气动弹性动态方程,然后基于状态依赖Riccati方程推导了非线性颤振控制律.假设只有俯仰角位移和控制面偏转位移可以直接测量,状态空间中的其它变量通过所设计的状态观测器进行了估计.仿真结果显示,观测器可以快速准确地对非直接测量变量进行估计,系统状态变量与控制变量都能够迅速地收敛于零点,表明所设计的控制律可以有效地实现对间隙非线性二元机翼颤振的抑制.      

用于叶轮机械流场计算的交错网格有限体积法

在对各种交错网格离散格式的奇偶不耦合数值误差进行了详细的傅里叶分析和比较后,采用一种紧致的交错网格有限体积法将纳维-斯托克斯方程和k-ε湍流模型方程紧密的耦合在一起求解,发展了一套三维定常流程序.通过对NASA Rotor 67设计和非设计工况内部流场的数值模拟结果与实验的比较、分析,初步验证了程序的准确性和可靠性.      

对旋轴流风机非定常流场数值模拟

利用全隐格式的非稳态压力修正求解方法,分析和研究了对旋叶轮级间流场非定常流动情况,上游尾迹与下游叶片位势作用对流场的影响,并与动态实验结果进行了有益的对比.结果表明:非定常流场的计算更能反映叶轮机实际运行性能,如果能很好地组织对旋叶轮的非定常流动,还可以进一步提高其压比和效率.     

绕大钝头短体外形跨音速大迎角湍流流动计算

采用数值模拟求解N-S方程方法计算与分析了绕大钝头短体外形的跨音速大迎角时的湍流流动.数值方法中选用Van Leer 分裂格式离散无粘通量项,并构造了一种Limiter函数以保证TVD(Total Variation Diminishing)性质,采用中心差分格式来离散粘性通量项.计算中采用了B-L湍流模型.算例定性地表明绕流的三维分离具有准周期性的非定常特性.     

几何和流动参数对上翘后体阻力系数的影响

利用N-S方程对9个上翘后体模型进行了气动力计算,主要研究后体几何参数和流动参数对上翘后体阻力系数的影响.研究结果表明, 后体的压差阻力系数分别随上翘角、收缩比的增加及迎角的减小而明显增加;后体摩擦阻力系数分别随后体的长细比的增加和雷诺数的减小而增加;后体越扁平,其压差阻力系数越大;在跨音速时,波阻系数也与上翘角有关,上翘角增加会导致波阻系数进一步加大.     

轴流压气机内导叶/转子干涉对转子流场的影响

叶轮机械内的转静干涉现象多年来一直是一个研究的热点,然而现有的对于转 静干涉的研究多局限于对流动现象的描述,缺乏理论指导,因此无法在叶轮机械的设计体系 中包含这一非常重要的非定常流动现象,限制了叶轮机械性能的进一步提升.利用流动稳定 性及感受性理论来指导对这一问题的研究,能够找到利用这一现象提升叶轮机性能的规律. 对于一台单级低速轴流压气机内的导叶/转子干涉现象进行了实验研究.利用二维粒子图像 速度仪测量得到了导叶/转子干涉的定量流动细节,及其对于转子流场的影响.结果表明, 转子流场在定常来流条件下也表现出显著的非定常性.在引入导叶/转子干涉之后,转子流 场结构发生显著变化.动叶尾迹强度明显减弱,叶背流动分离得到抑制,叶尖泄漏涡也得到 了周期性的抑制.除了导叶尾迹对动叶流场的干涉作用以及气流角的周期性变化以外,动叶 非定常流动对于不同的导叶/动叶干涉激励频率的不同响应是引起压气机性能变化的主要原 因.      

超声速进气道可压及不可压流动数值模拟

采用Chio-Merkle预处理矩阵对可压NS(Navier-Stokes)方程组进行时间预处理,分析了预处理方法的物理和数学背景.用有限体积方法,结合LU-SGS(Lower-Upper-Symmetric-Gauss-Seidel)隐式时间积分和AUSM+(P)(Precondition Advection Upstream Splitting Method)格式、中心差分2种空间离散格式,求解该预处理NS方程组.通过圆弧凸包流动和二维方腔顶盖驱动流动的数值试验表明,该方法克服了传统时间迭代方法模拟低速流动时的刚性问题,加速了收敛过程,可以同时有效地模拟可压及不可压流动.用此方法与块结构化网格技术相结合,进行了在不同飞行马赫数、攻角和出口反压条件下,混压式轴对称超声速进气道三维流场数值模拟.计算结果表明:该方法能准确地捕获复杂激波系,清楚地揭示进气道三维流动现象;进气道性能随工作条件的变化规律,与理论分析一致.      

利用CFD技术计算飞行器动导数

利用非定常流场数值计算方法模拟飞行器强迫俯仰振荡仅能得到俯仰力矩系数对迎角变化率和俯仰角速度的动导数之和,而动稳定性分析需要单独的动导数数值.为解决这个问题,利用滑移网格模拟强迫振荡运动,得到俯仰力矩系数对迎角变化率和俯仰角速度的动导数之和.利用旋转参考坐标系模拟定常拉升,得到俯仰力矩系数对俯仰角速度的动导数.利用旋转参考坐标系模拟匀速滚转,得到滚转力矩系数对滚转角速度的动导数.对有翼导弹和水上飞机进行了纵向和横向动导数的计算.计算结果与试验数据、文献数据以及其他方法得到的结果具有较好的一致性,表明提出的方法可用于复杂外形飞行器动导数计算.     

轴流压气机尾流撞击效应机理的数值模拟

通过精确给定转子进口边界条件,以模拟上游静子尾迹对转子流场的影响,这样利用单排的计算结果模拟转/静干涉引起的非定常流动对流场气动参数的影响,因为相邻叶片排的叶片数不同,从而避免传统方法中常用的叶片数约化存在的难以正确模拟非定常激励频率的问题,对计算时间和计算能力的要求也可减少.针对一个低速转子的数值模拟结果显示,如果能合理地组织来流尾迹的扫过频率和尾迹亏损的幅度,对转子气动性能的改善是很显著的,从而为风扇/压气机中控制和利用非定常流动提供了一个可行的方向.     

基于润滑理论的二维积冰数值模拟

将基于润滑理论的水膜-冰层模型推广到在任意二维截面上建立的正交曲线坐标系,得到描述二维型线表面水膜流动和积冰的偏微分控制方程.采用一种隐式-显式有限差分格式建立控制方程组的代数形式,给出求解非稳态控制方程组的数值方法.为验证模型和数值解法的有效性,在典型的航空积冰和传输线积冰环境下对翼型和传输线表面的积冰算例进行数值模拟.将水膜-冰层模型的冰形计算结果与传统Messinger模型的模拟结果以及冰风洞试验结果进行对比,低温明冰条件下采用当前方法计算的翼型表面冰形接近于Messinger模型的模拟冰形曲线;相对高温条件下的计算结果比传统Messinger模型更为精确.低速积冰环境下Messinger模型难以模拟的传输线表面积冰同样可以采用水膜-冰层模型进行有效 预测.