使用GAO-YONG湍流方程组对扩压器流动的计算

采用基于同位网格的SIMPLE方法求解GAO-YONG不可压湍流方程组,对二维扩压器流动进行了数值模拟.通过界面速度动量插值法解决了压力锯齿波问题,并且采用了正交贴体网格和二阶QUICK格式离散对流项以提高计算精度.与实验结果以及BL(Baldw-Lomax)模型计算结果的比较表明,不需要任何经验系数及壁面函数的GAO-YONG不可压湍流模型方程组能够对有压力梯度的湍流流动做出很好的预测.计算发现,在机械能方程中引入平均流压力梯度的作用,对GAO-YONG湍流方程组正确模拟逆压力梯度流动起到了关键作用.      

塞式喷管三维流场的数值模拟

从曲线坐标下的三维平均雷诺的N-S方程出发,并用k-ε两方程湍流模型封闭方程组,采用二阶精度的精度无波动、无自由参数的耗散差分格式(NND格式),发展了模拟塞式喷管三维流场的数值程序,并对线性塞式喷管和瓦状塞式喷管在不同落压比下的三维流场进行了数值模拟研究.结果较好地预示了塞式喷管的三维流场特征,并且显示瓦状塞式喷管除了圆柱形结构比平板更能承载外,其塞锥压强分布可以进一步改善其受力情况.     

乘波构型飞行器磁流体进气道一体化概念设计

采用准一维的磁流体方程组对设计状态下内磁流体发生器进行参数设计,利用斜激波理论设计磁流体进气道一体化的外形.以该外形为基点,对安装五对电极的内磁流体发生器采用五方程模型进行三维数值模拟,并通过调节其磁作用数以满足准一维假设下设计的进气道外形.研究表明,磁流体进气道最小出口流动均匀参数可以作为判别内发生器最优磁作用数的准则,在该磁作用数下内磁流体发生器可以有效地提高进气道的性能.      

模拟可压缩流的高分辨率有限元法

采用经典的Galerkin方法对N-S方程进行空间离散,再以该半离散格式作为基本格式,参照Kuzmin的思想加入耗散与反耗散项,使所构造的格式具有局部极值不增(LED,Local Extremum Diminishing)性质.对上述半离散格式进行时间离散后所导出的稀疏线性代数方程组,采用了GMRES(Generalized Minimal Residual)迭代法进行求解.为验证所建立的格式及相应的程序,给出了激波管问题和绕圆柱与双椭球超声速流动问题的数值模拟结果.      

基于SA湍流模型的DES方法应用

使用基于SA(Spalart-Allmaras)湍流模型的分离涡模拟(DES,Detached-Eddy Simulation)方法,对圆柱绕流和超声速圆柱底部流动等高雷诺数大分离流动进行了数值模拟.无粘通量的计算采用二阶迎风格式,时间推进使用二阶隐式双时间步法.计算表明相对于雷诺平均N-S方程(RANS,Reynolds-Averaged Navier-Stokes)模拟,时空均为二阶精度的DES模拟能够分辨更为精细的分离区三维旋涡结构,对分离区主要湍流结构的捕捉和平均压力分布的预测取得了满意的结果.     

Euler方程数值模拟绕悬停旋翼桨叶的流动1)

本研究应用格心格式有限体积龙格-库塔时间推进法求解三维欧拉方程，模拟绕旋翼桨叶的流动.旋翼尾流的作用通过自由尾流分析来求解局部的诱导下洗速度以修正翼型攻角.选择O-H型式生成绕桨叶贴体的弦向和展向网格.通过算例计算为进一步寻求有效的Euler或NS方程数值解积累了经验.     

三维发动机进气道水滴撞击特性分析

用带数值扩散项的欧拉法对某型三维发动机进气道水滴撞击特性进行了研究.实现了可压流场中三维复杂形体的局部水滴收集系数的计算.数值扩散项的引入解决了欧拉法计算产生奇异结果发散的问题.通过对不同进气道流量的算例结果的比较,发现进气道流量对进气道唇口附近的水滴撞击区和撞击量都有影响.进气量大,最大局部水滴收集系数也大.通过对有无整流罩的计算结果进行比较,发现此型发动机整流罩对水滴撞击特性的影响很小.      

不同厚度翼型动态失速涡运动数值研究

在低马赫数下, 对三种不同相对厚度的NACA系列基本厚度翼型在俯仰振荡运动中的动态失速现象进行了数值研究.数值模拟时采用双时间法和LU-SGS隐式解法相结合的模式求解了非惯性坐标系下的非定常纳维-司托克斯(NS)方程组, 空间离散采用Roe格式,并结合刚性动网格生成技术.采用全湍流计算,通过引入BL(Baldwin-Lomax)模型计入湍流影响.计算出的气动力迟滞曲线与实验结果变化趋势符合较好,表明了数值方法的有效性,同时通过比较和分析不同厚度翼型在轻失速和深失速下的流场结构,发现绕翼型的失速旋涡产生和发展规律是明显不同的.      

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

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

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

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

乘波机外形设计

计算了由锥型流理论、吻切锥理论及吻切轴对称理论生成的不同乘波机外形.结果表明,圆锥绕流流场生成乘波机外形的方法虽然比较简便,且生成的外形具有较高的升阻比,但进气道进口流场具有锥型流的性质而不够均匀,给发动机的工作带来不利影响.应用吻切锥、吻切轴对称理论与特征线方法结合起来进行反设计,可以得到更一般的乘波机外形,并且能有效地改善进气道进口流场,提高容积效率.此外,这几种方法的计算速度都很快,对进一步的优化处理及机身-发动机一体化十分有利.      

谐波平衡法在低速非定常流模拟中的应用

谐波平衡法是一种有效的周期性非定常流的计算方法.采用基于可压缩流的谐波平衡方程在计算低速不可压流动时,会由于对流通量计算格式中的数值粘性污染,降低解的精度和收敛性.采用预处理技术,使得基于可压缩流的谐波平衡方程可以直接用于低速周期性非定常流的计算中.选取典型的不可压方腔驱动流和低雷诺数圆柱绕流为例进行了时间推进法和谐波平衡法的计算对比.计算结果表明预处理后的谐波平衡方程适合于低速流的计算,在谐波平衡法中采用较少阶数的谐波计算就可以还原出几乎准确的非定常流场.     

高阶熵条件格式下Euler方程与 N-S方程的混合算法

针对在工程问题中使用N-S方程计算量大的问题,提出只在粘性作用显著的局部区域求解N-S方程,在流场的其余部分采用Euler方程计算的办法来模拟具有复杂外形的流动问题,并提出了一种基于高阶熵条件格式的算子分裂算法.应用所构造的算法对绕迫击炮弹的亚、跨超声速流动进行了模拟计算.计算结果与试验数据的对比表明,所提出的方法是切实可行的.      

高超声速三维化学非平衡流动的数值模拟

高温非平衡化学反应对高超声速再入飞行器绕流流场的气动热特性有较大影响.通过采用有限体积NND(Non-oscillatory, containing No free parameters and Dissipative scheme)格式求解三维N-S方程,分别计算了小攻角条件下球锥的完全气体粘性流场和化学反应非平衡气体流场.比较并分析了化学非平衡作用对三维高超声速钝头体绕流流场的影响.所得结果与国外文献结果符合的较好,证明了计算方法与程序的有效性和实用性.      

直升机发动机入口绕流的数值模拟

对直升机发动机入口处流场采用有限体积法,通过生成贴体结构化网格求解Navier-Stokes(N-S)方程, 进而对MI-171直升机低空悬停和前飞状态进行数值模拟,计算了该型直升机发动机进气道处的流场特性.然后以广义尾流入手,并引入尖涡涡核作用的半经验修正的旋翼自由尾流模型对该机发动机进气道处流场进行计算分析.对沿发动机轴线剖面以及垂至于该轴线剖面的分析结果表明,CFD(Computational Fluid Dynamics)数值模拟方法与旋翼自由尾流分析方法均可以对涡诱导流场进行预测.     

基于涡方法生成大涡模拟进口条件的数值计算

探索改进涡方法来生成大涡模拟的非定常进口条件.改进方法中,为了避免局部漩涡数量过多或者过少,采用密度分布方式放置漩涡场.并用Langevin方程控制漩涡运动,模拟实验方法中的蜂窝器,使改进后的涡方法生成的脉动速度场更加符合湍流的特征.在已知雷诺平均的流场结果下,利用涡方法产生漩涡场,进而生成能满足大涡模拟所需要的非定常进口流场.为了检验改进后的涡方法在生成脉动速度过程中的效果,在槽道中进行了对比数值试验,并借助于直接数值模拟数据做对比,对比分析槽道进出口的平均速度、涡量以及雷诺应力统计,证明改进后的涡方法在生成大涡模拟进口条件下是非常有效的.      

二维微槽道电渗流的增强混合

微尺度流动中的混合一直是微流系统研制中的一个大问题,对壁面zeta势非均匀分布情况下的二维微槽道电渗流进行了数值分析,重点分析了这种情况下的微流动混合增强.基于电渗流的Helmholtz-Smoluchowski滑移速度模型,应用SIMPLE方法计算了4种滑移速度在整个壁面上的分布不均匀一致时的微槽道电渗流动的流场及被动标量在其中的混合.计算结果表明,滑移速度的非均匀一致导致微槽道电渗流中的流线发生剧烈的转折,从而使得微槽道流动的混合效率有很大的提高.流场中流线的转折越剧烈、转折的次数越多,微槽道流动混合效率的提高就越大.     

发动机进气道水滴撞击特性分析

针对某发动机进气道的防冰系统设计,对进气道的水滴撞击特性进行了研究.在对进气道流场进行计算的基础上,采用差分法对水滴运动方程进行了数值求解,得到了水滴运动轨迹,从而确定了水滴的撞击极限、总收集系数和局部水收集系数等水滴撞击特性参数;此外,还研究了飞行高度、飞行速度及水滴半径对水滴撞击特性的影响;研究发现水滴撞击极限、总收集系数和局部水收集系数随飞行速度的减小、飞行高度的增加或者水滴半径的增加而增大.这些为进气道的防冰系统设计奠定了基础.     

高温化学非平衡计算中源项处理方法

高超声速三维化学非平衡绕流流场的数值模拟中,在计算初期容易发生组元密度出现负值的非物理现象,另外源项的刚性是影响计算稳定性和收敛速度的主要原因。根据源项线性化理论,采用2种源项线性化方法处理化学反应源项.一种为根据流场内化学反应物理规律构造的线性化方法,该方法能抑制在计算过程中组元密度出现负值现象,提高了计算稳定性,加速了收敛速度.同时采用了时间预处理矩阵的线性化方法,较好地解决了非平衡化学反应与流场耦合的刚性问题。证明了所构造源项线性化方法相容性.数值实验表明该方法有效地避免了源项计算中密度出负的问题,加快了计算的收敛速度,从而提高了计算效率.