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

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

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

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

滚珠自转对高速球轴承环下润滑内部流动的影响

高速球轴承正常工作时,处于高速旋转状态,此时油相和气相因重力和离心力的作用在轴承环间剧烈运动。为了更准确地分析轴承环间的两相流动,采取VOF模型进行内部流场的模拟,采用多重旋转坐标系描述部件运动。建立球轴承环下润滑计算模型,分析了考虑滚珠自转因素下轴承内部的流动,并在此基础上探究了转速及供油量对轴承工作状态的影响。结果表明,对比整体模型与滚珠自转模型,发现滚珠自转使得轴承内部油体积分数增大,同时也使得滑油穿透间隙达到外环的能力增加;在考虑滚珠自转情况下,转速的不断增大,使得轴承内部的油相体积分数不断减小,在低转速情况下滚珠自转对流体运动影响较为明显,在高转速情况下公转速度对流体运动起到主导作用,滚珠自转对流体运动影响减弱;供油量不断的增大,使得滚珠自转模型内部的油相体积分数也在不断增大,而且滚珠自转运动会加强滑油在轴承内部的分布。     

类X?47B无人机菱形编队气动干扰数值模拟

设计了一种由4架类X?47B飞翼布局无人机（UAVs）组成的菱形编队。通过求解RANS方程的数值模拟方法,研究了菱形编队无人机气动干扰问题,详细分析了影响机理,定量给出了编队减阻效果。计算结果表明:头机气动性能基本保持不变。两侧僚机受上洗气流影响,其减阻效果明显。尾机主要受下洗气流影响,其阻力增大,对编队久航和远航不利。在重力配平条件下,两侧僚机飞行阻力的减小是由攻角减小和诱导阻力减小共同引起的。尾机在编队中飞行阻力的增大主要是攻角增大带来的阻力增加,诱导阻力增大仅带来了20%的阻力增量。从减小下洗气流对尾机的不利影响出发,对不同垂向间距的尾机升阻特性进行了研究,并参考雁群头鸟变换行为机制,给出了无人机菱形编队飞行建议。      

高速可压缩来流与凹腔驻涡相互作用及其物质输运特性研究

为了深入认知凹腔驻涡与高速可压缩来流之间的相互作用过程以及来流马赫数对上述动力学过程的内在影响规律,对长深比为1.0的单凹腔基础构型在常温、常压、Ma=0.3~0.7的高速可压缩来流下进行冷态非定常数值模拟,并基于拉格朗日视角研究来流与凹腔驻涡相互作用及其物质输运特性。结果表明：凹腔上方剪切层内旋涡结构仍然具有类似自由剪切层的周期性演化过程,主流与凹腔通过三次掺混进行物质输运;随着主流马赫数的增加,旋涡运动方式及凹腔-主流物质掺混机制都发生了本质上的改变,剪切层内旋涡和凹腔主涡涡量均呈现数量级式增长,剪切层旋涡结构脱落频率近似线性增长并逐渐卷吸至凹腔中部,进入凹腔的主流流体占比降低,但扩散程度和速度加快;此外,源于凹腔的流体流出唇口线后经剪切层的夹带作用回流占比增大,降低了凹腔-主流的整体质量交换水平;驻留参数整体衰减,高马赫数时仅为低速工况的1/2,显著降低了回流区对于流体驻留作用的贡献。     

亚格子模型对泵喷推进器宽带非定常力预测的影响研究

为了研究亚格子模型对泵喷推进器非定常流动与宽带非定常力预报结果的影响,采用分块结构化网格建立了模型尺度下艇后泵喷推进器的计算模型,并进行了大涡模拟数值计算。从艇尾非定常流场特征量和推进器转子脉动载荷两个方面对比了三种不同亚格子模型计算结果的差异,并分析了泵喷内部流动与转子非定常力间的内在联系。研究结果表明：三种亚格子模型均能得到含有叶频宽带峰的轴向推力谱,且整体趋势相近。泵喷转子上游的湍流强度和尺度的分布对亚格子模型较敏感,其中Smargrinsky-Lilly模型得到的湍流强度较强,尺度较大,该模型下的湍流谱在非平衡区量级较大,但由低频向高频的衰减较快,并且预测到的分离区范围大,导叶尾缘脱落涡分散,导致转子上游来流空间分布不均程度较强。对于转子叶片上的载荷脉动,Smargrinsky-Lilly模型预测的推力谱中线谱成分明显,并且叶频处宽带谱峰“陡峭”,而WALE模型和KET模型的结果宽带特性较强,对比非定常推力测试结果可知,WALE模型和KET模型更适于该问题宽带非定常力预测。     

风扇/增压级带间隙三维粘性流场数值模拟

使用三维数值模拟软件NUMECA(Numerical Mechanics Application)对双涵道风扇/增压级带转子叶尖间隙的粘性流场进行了三维定常流动数值模拟,获得了其设计转速下的流场特性.结果显示,随着外涵背压的升高,风扇顶部激波后阻塞区周向和轴向尺寸增大,激波变形加剧,继续提高外涵压比,则外涵效率和进出口流量下降.内涵背压的升高,使内涵各转子叶片顶部的泄漏流动影响范围向上游扩散.内涵背压过高,效率和进出口流量同样会有下降的趋势.同时通过和实验结果的比较表明,数值模拟软件模拟的风扇/增压级的三维定常粘性流场能够较真实的反映出转子叶尖间隙泄漏流场的结构特征,以及间隙对流场参数的影响.      

适于旋流杯下游流场的改进型半分析关系式

利用典型(CFM-56燃烧室)旋流杯结构下游的流场实验数据,参照经验关系式,提出了物理意义更明确、拟合更准确、且更便于对大量实验数据进行相关性研究的半经验关系式,并与旋流杯下游紧靠出口处附近的轴向速度及旋流速度测量数据进行了对比,发现各特征数据(如峰值位置、峰值衰减速度、回流区宽度及其沿下游的发展、回流速度及其衰减等)与实验数据非常相符,从而初步验证了提出的半经验关系式.同时,给出了公式中不同的经验常数取值带来的变化,分析了与几何及工况参数之间可能的联系,为进一步验证该关系式打下基础.分析表明,各经验常数的变化之间相互独立,且物理意义明确,可以方便地用于进一步验证.      

层流平板摩擦阻力的数值计算

为提高层流摩擦阻力数值计算的精度,采用数值求解N-S方程的方法,对层流平板摩擦阻力进行了计算.计算中采用4种先进上风格式,同时配合以14种不同的壁面法向网格分布和5种常用限制器,并将计算结果与布拉修斯解进行对比.结果表明,van Leer的FVS(Flux Vector Splitting)格式数值耗散大,对网格分布和限制器敏感,不适合用于摩擦阻力计算.限制器和计算格式的粘性分辨率越高,摩擦阻力的计算精度越高.在计算量允许的情况下,可以减小壁面法向第1层网格高度,以提高计算精度,但应同时兼顾网格长宽比.对于摩擦阻力计算,边界层内应至少布置4个网格点.      

脉冲射流强化圆柱喷流混合的三维数值模拟

对利用脉冲射流强化圆柱喷流混合的流场进行了数值模拟,这项工作将高频高强度信号射流射入亚音剪切层以强化喷流混合.利用有限体积法和重整化群(RNG)k-ε湍流模型求解N-S方程.分别进行了低亚音马赫数及高亚音马赫数喷流直接排入空气的数值模拟.主要计算3个状态:未受激励的状态,斯德鲁哈尔数St分别为0.2和0.4时的状态.数值模拟结果显示,脉冲激励引起喷流在与激励源平行的平面迅速扩展,在与之垂直的平面收缩,当脉冲射流激励的St在0.2附近时,相关的混合较好,而且少量的射流流量可以极大地增加速度的衰减率.