迷宫螺旋泵和密封的数值分析

 对迷宫螺旋体内的流动提出了若干假定;建立了流动模型;简化了Navior-Stokes方程和连续方程。选用轴及套筒上两组斜交螺纹所围成的蜂窝状空间作为计算区。在计算区内截取两组沿螺纹方向的斜交截面,采用涡量流函数法求出截面上的平面流动解。根据两组截面交线上速度唯一的条件,用一组平面解的已知速度作为给定数据输给另一组平面求出其平面解。如此将两组平面解交替迭代,最后得到的收敛值即该区的空间流动解。将所有蜂窝体内的压力和流量综合起来,就得到迷宫螺旋体的泵特性和密封特性。     

两相喷管流动研究进展

本文综述了两相喷管流动研究的进展,着重介绍了两相喷管流动控制方程的数值解以及两相喷管流动研究的实际应用。     

混合气体流动的一种数值模拟方法

本文建立一个数值模拟完全气体混合流动的理论模型.该模型首先应用混合气体的Euler方程和每种气体组分的质量分数方程来控制流动.为了消除混合网格内气体组分界面附近出现的非物理振荡,我们假定混合气体的每种组分达到了动力学平衡状态然而尚未达到热力学平衡状态.这种思想导致需要另外给定每种气体组分的总能量方程.为使用高分辨格式来求解这组双曲型偏微分方程并且简化对所需要的Jacobi矩阵的推导,混合气体的压力方程也被耦合起来.Godunov型的波传播方法被采用来离散求解所获得的控制方程.从典型算例结果来看,一维问题的数值解与精确解一致,二维问题的数值解与理论分析一致.这说明本文的理论模型是合理的.     

对称驱动方腔定常流动的多解现象研究

对方腔内顶盖正弦速度驱动的定常流动进行了数值研究.采用非线性动力系统数值分岔分析方法,讨论了Re<2000时Navier-Stokes方程定常解的多解共存现象,计算了对称和非对称定常解的分岔曲线和相应的流动模式,给出了非对称解出现的临界参数Re=320,稳定的对称和非对称解同已有文献的结果一致,并获得了一组新的不稳定的非对称解.随着Re数的增加,侧壁出现流动分离并有小涡出现,随之发生流场拓扑结构分岔,但并不伴随N-S方程解的分岔.     

非零压差边界下间隙流动雷诺方程近似解析解

为了获得间隙密封、滑动轴承等的液膜压力分布解析表达式,将间隙流动雷诺方程的解析解表示为特解和通解之和,其中特解表示为周向坐标函数和轴向坐标函数之和,通解表示为周向坐标函数和轴向坐标函数之积,结合分离变量法和非零压差边界条件,推导获得间隙流动雷诺方程的近似解析解。基于近似解析解计算间隙密封的液膜承载力和偏位角等参数,并与采用有限差分法的计算结果对比,其中液膜承载力偏差小于5%,偏位角偏差小于3°,而计算效率提高了数百倍。      

等速上仰翼型绕流结构的数值模拟

通过数值求解涡量-流函数形式的N-S方程研究等速上仰翼型的分离流动结构。用ADI方法解涡量方程;用Poisson方程直接法解流函数方程。在计算得到的流场中除可看到实验中观察到的多涡结构外,还发现一些实验中未观察到的现象。将分离流动结构的计算结果与翼型气动力变化相联系,并结合涡动力学理论,定性地解释了上仰翼型产生高升力的原因。     

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

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

跨、超音速流动的区域分解方法与并行算法

研究二维跨、超音速无粘性流动的 Euler方程区域分解方法、并行算法及其应用。通过内边界耦合条件实现相邻子区域解的光滑过渡,以得到总体流场的数值解。发展了一种多块区域之间守恒型的有效内边界耦合方法,对二维翼型跨音速流动和钝头体超音速流动等进行了分区数值求解,分区计算结果与其他单区计算结果作了比较,并讨论了多种区域分解数目的分区计算效率。并行计算采用纯结点并行编程方式和“先进先出”的同步控制等待机制,利用 PVM并行环境对二维绕翼型跨音速流动做了二区和四区分区并行计算。     

平面叶栅跨声速非等熵流的时间推进解

本文根据Viviand等人提出的求解平面流动问题的拟不定常时间推进法的思路,推导了完整的求解平面叶栅无粘非等熵拟不定常流动的方程系统,应用经过改进了的有限控制面积时间推进格式进行了数值计算。计算结果与精确解或试验结果吻合较好,且能有效地节省计算机时。     

粘性跨音速喷管流场计算

本文继用隐式近似因子分解法成功地计算无粘跨音速喷管流场之后,用同样的方法,结合任意曲线坐标系,通过非定常方程在相当长时间后的定常解,针对矩形截面喷管和轴对称喷管粘性跨音速流场,求解了可压缩层流薄层N-S方程,获得了核心流,特别是边界层中的流动参数.对于矩形截面和轴对称两种喷管进行了计算,其结果和实验数据相当一致.将本文的方法应用到粘性两相流动计算,可望较多的节省机时.     

高超声速球锥粘性激波层压力分裂法数值求解

本文将压力分裂法引入细长钝体高超声速粘性激波层数值求解中,较好地解决了这种流场远后身部分的推进求解问题。针对非解析球锥体绕流数值求解的困难,作者提出了实用型初始激波形状计算方法,达到了简化数值求解步骤,提高数值计算效率的目的,数值求解结果与实验结果相符合。     

高阶精度间断Galerkin有限元方法研究

在二维非结构网格上,对高阶精度间断Galerkin有限元方法求解跨音速欧拉方程进行研究。运用间断有限元理论,采用施密特正交化多项式基函数对流场解进行近似描述,使用HLLC近似黎曼解方法计算网格单元边界处的数值通量,积分项通过高斯积分求解,时间推进采用经典四步Runge-Kutta方法,并引入斜率限制器,抑制数值振荡。对NACA0012翼型跨音速无粘流动进行数值模拟,数值结果表明：间断Galerkin有限元方法具有较高的精度,较小的数值耗散和较强的激波捕捉能力。     

基于TFI的局部网格变形方法研究

本文发展了一套基于TFI（Transfinite Interpolation,超限插值法）技术的局部网格变形方法,该方法不同于以往的整体网格变形思路,将网格变形控制在有限的区域内,在保证网格质量的同时,不影响其他区域的网格,因此能够较好地配合其他动网格方法解决复杂流动问题。以带后缘襟翼的NACA0012翼型为例,对翼型俯仰振荡耦合襟翼偏转运动的动态流场进行了数值模拟,计算结果与实验值吻合良好,证明了本文发展的局部网格变形方法的可行性。     

跨音涡轮中三维流动的数值试验

 本文从积分型的基本方程组出发,以我们的关于平面叶栅和任意旋成面叶栅的工作为基础,将文献（3、4）中的计算方案推广到三维流动情况。编制了相应的计算机程序,进行了一些数值试验,得到了初步的计算结果。 数值试验表明,拟三维计算结果可以导致相当大的误差,从而说明了完全三维流动计算的必要性。另一方面,如果取2000个左右的计算点,则在国产的013机上算一个例子,约需100分钟。从实际应用的观点来看也是可行的。     

解跨音速升力翼型的有限元法

 在应用解全速位方程的最小压强积分有限元法求解绕升力翼型的跨音速流动时,将不可压流中求解绕升力翼型的耦合单位环量流动和无环量流动的解法推广到可压流中。为了确定环量,本文所用Kutta条件是:在后缘处,气流流向平行于后缘角二等分线。因有限元法对网格无正交性要求,因而可在椭圆变换前后进行剪切和延伸变换。这种网格生成法易于构成适用于复杂形状的有限元网格。通过计算并将其结果与文献中的数据比较,表明这种方法应用方便且有较快的计算速度和较高的计算精度。     

近水面飞行气动特性数值计算方法研究

根据对水面波的数学描述，给出了其对流场影响的精确解析解。针对地(水)面效应飞行器的特点，分析了近水面波飞行时的流场特征，利用线性叠加原理，提出了飞行器与水面波耦合影响的非定常速度势修正项概念。分析了非定常气动特性的数值计算方法。对某实用型地效飞行器进行了计算，并对计算结果与试验结果进行了分析总结。     

等截面矩形隔离段内流场的三维数值模拟

针对三维矩形等截面隔离段中由激波／附面层干扰诱导的复杂流场，采用有限体积方法，OC－TVD格式结合B－L代数湍流模型隐式求解了三维雷诺平均N－S方程，数值模拟了三维流场内复杂的激波串和旋涡结构，计算结果表明所采用的数值方法能够很好地捕捉隔离段内流场的复杂流动现象。     

用弱非线性理论研究二次稳定性问题

本文把研究三波共振的弱非线性理论应用于平板边界层和平面Poiseuille流三维扰动的二次稳定性问题。其数值结果与文献[1]、[2]用Floquet形式的二次稳定性理论得出的结论是一致的。     

单级跨声速压气机进气畸变数值模拟

基于彻体力模型的基本思想,发展了一套能够反映进气畸变对风扇/压气机气动性能影响的三维数值计算程序CSAC.首先详细研究了彻体力模型的具体建模方法,然后发展了相应的三维数值计算程序.利用该程序针对某单级跨声速压气机在均匀进气下的流场进行了计算验证,计算结果与雷诺平均Navier-Stokes(RANS)计算吻合得很好,流场主要参数分布与相关实验数据也有很好的一致性.利用该程序分别计算分析了进口存在稳态周向总压畸变及周向总温畸变的全环三维畸变流场,结果表明:压气机进口畸变流场存在大尺度、强三维特征,将导致压气机总压比及稳定裕度下降.所有计算结果均显示该模型能够在降低对计算资源和工程经验依赖的同时,有效分析进气畸变对压气机特性的影响.      

二次空气系统盘腔一维非定常计算方法

通过对盘腔结构的一维算法展开研究,提出适用于具有非标准结构、包含多个进出口的复杂盘腔的一维建模规则,并采用特征线法,实现盘腔结构的一维非定常流动与换热的计算,通过将广义的变截面、有摩擦、考虑对外热交换的盘腔一维可压缩非定常流动与换热的偏微分方程组转换为常微分方程组,进而获得盘腔内部沿流动方向压力、流量、温度等系统参数的特征线数值解,采用商业计算流体软件CFX对其进行验证.结果表明:所提出的盘腔建模与计算规则能够直观且准确地描述盘腔一维非定常流动及换热特征,仿真结果与CFX计算结果最大偏差约为9.4%,可应用于具有非标准结构的多接口复杂盘腔的模拟,可以作为二次空气系统盘腔结构一维非定常计算研究的有效手段.