改进的翼型积冰数值模拟方法

基于微团概念,以雨滴运动的基本假设出发,对翼型积冰的拉格朗日方法进行了改进.新方法克服了原方法计算效率低、冰形同雨滴分布密度及网格密度相关的缺点,继承了原方法的鲁棒性,不存在欧拉方法的收敛性问题.为验证方法的正确性,计算了NACA0012翼型在0°和4°攻角下的霜冰积冰情况、表面水收集量和局部撞击系数.结果表明新方法减少了计算量,预测的冰形和实验结果一致,并分析了积冰对气动性能的影响.      

混合解析/数值方法及其在湍流数值模拟上的应用

针对带小时间尺度的源项的方程描述的流动问题,提出了混合解析/数值方法。混合解析/数值方法的基本思想是:分裂原始方程组为对流-扩散部分的偏微分方程和源项的常微分方程。偏微分方程采用合适的数值方法求解,而常微分方程采用解析方式积分。模型方程的理论误差研究表明,混合方法提高了源项处理的精度,降低了混合方法的整体数值误差。分析同时表明,基于时间分裂的算法在求解含源项双曲系统的定常类型问题,会存在数值振荡。为此发展了非分裂方式的混合解析/数值方法,在湍流模型数值计算中提高了数值稳定性,而且加快了计算的收敛速度。     

带有双层壁扩压器的波瓣喷管混合流场的数值计算

采用三维贴体曲线坐标网格,边界网格加密且正交,在整个计算区域进行全场计算。为避免因波瓣造成的网格强烈的非正交引起解的发散,采用了Chen和Kim修正的k-ε湍流模型及同位网格SIMPLC计算方法,对带有双层壁扩压器的波瓣喷管混合流动进行了数值计算和分析。对波瓣及双层壁,采用大粘性的方法解决流固耦合。计算结果表明:在双层壁间有外界冷气流被引射进入,形成了壁面的冷却气流;自波瓣出口截面沿流向产生的环流速度场,强化了主次流的掺混,速度分布渐趋均匀。计算结果与实验数据符合良好,二者在离开波瓣52 mm的混合管内相对主流速度的最大误差为21.65 %。     

叶轮机非定常气动设计的缘线匹配技术

本文提出一种新的叶轮机气动设计自由度——缘线匹配技术。缘线匹配是指相邻两叶排中前排叶片后缘线与后排叶片前缘线的空间关系,更准确地说是其相位角关系。叶片尾缘线表征了叶片出口气流的分布型,它囊括了尾迹、激波、二次流等所有影响,另一方面,叶片前缘线代表了其对上游流场的势干扰。定常气动设计体系只能部分考虑叶片缘线的相互影响,大部分则被忽略了,然而在非定常框架下,相邻叶排缘线所代表的相互影响是显著的。本文首次提出了叶轮机非定常设计的缘线匹配技术,并通过理论分析及直列涡轮叶栅的非定常数值模拟结果展示了缘线匹配技术在未来提高叶轮机气动性能、气弹性能、气动噪声和热传导性能上的潜在能力。作者认为,缘线匹配技术将使叶轮机的非定常气动设计具有真实而可操作的内容,是叶轮机非定常气动设计的核心技术之一。      

二维非结构网格的生成及其Euler方程解

 采用推进阵面法生成三角形非结构网格,用多重网格技术求解泊桑方程以取得结构背景网格.提出了一种高效的边界剖分法,并用它来生成初始阵面.另外,还提出了一种新颖的数据结构——数组链接表,用它来存贮网格数据.用有限体积法在非结构网格中求解了二维可压缩Euler方程,给出了单段翼型和多段翼型的算例,并与实验结果进行了比较和分析.     

爆炸波与物体干扰流场的数值模拟

本文构造一类新型的计算效率较高的ＥＮＯ差分格式。利用该格式非定常Ｅｕｌｅｒ方程出发,对爆炸波扫过物体时形成的二维和三维干扰流场进行了数值模拟,计算结果与实验数据符合较好。     

地球流体力学数值方法的若干研究

本文简要综述了三十多年来中科院大气物理所发展的一些地球流体力学数值方法及有关问题的研究。其中包括用物理观点和数学分析相结合的方法阐述造成计算紊乱和计算不稳定的机理;论证计算稳定性与算子非负性和能量守恒性之间的密切关系;对国际上流行的瞬时能量守恒的格式找到非线性计算不稳定的特解;设计出一批内部协调并且具有完全平方守恒性的差分格式（包括隐式的、显式的和高精度的）;发展了多种经济省时算法等等。同时也简要     

椭圆形截面螺旋管道内二次流动的摄动解

从曲线直角坐标系中的Ｎ－Ｓ方程出发,采用摄动方法求解了椭圆形截面螺旋管道内粘性流动。结果表明：一阶挠率项仅影响椭圆形截面上的二次流动,并不影响截面上轴向速度;当轴向压力梯度较小时,二次流动主要受挠率的影响,随着轴向压力梯度的增加,曲率的作用越来越强。     

导弹配平飞行数值模拟及控制优化

传统基于建模仿真的控制系统检验方法受建模质量的影响,在非定常效应较强时很难得出准确结果。为此,采用基于CFD的虚拟飞行数值模拟方法检验和优化控制系统。首先介绍了数值模拟方法,进行了纵向开环控制飞行模拟,结果与气动力建模仿真结果、风洞试验结果吻合较好,共同验证了计算方法和气动力模型。针对迎角拉起设计了PID控制器,阶跃响应仿真结果表明控制器具有良好的动态性能,然后进行了数值模拟验证,结果表明迎角超调量超过50%、调节时间1.7s。通过试探法对控制参数进行了优化改进,将超调量减小到14%,调节时间减小到1.18s,动态性能显著提高。分析表明,受控时尾舵角速度达到200°/s,非定常效应强,线性气动力模型不能精确描述这些过程,造成了仿真结果的差异;虚拟飞行数值模拟技术在控制律检验和优化上可以起到关键作用。