基于迎风型有限体积法的非结构化网格流场格心数据高质量体绘制方法

在计算流体力学领域,对复杂流场进行数值模拟时常用非结构化网格划分空间,并广泛采用迎风型有限体积法求解流场.该方法将数值解存储在非结构化网格单元(控制体)中心,产生格心数据,并假设物理量在相邻单元邻接面上形成"间断".然而,在可视化此类数据时,目前的重构方法破坏了间断约束,导致"流间断"信息丢失,严重制约了绘制质量.为此...     

高超音速流场的隐式TVD格式并行计算

本文采用区域分割技术和拼接网格的并行策略，发展了一个适合于分布式存贮多机系统的ＴＶＤ隐式有限体积并行算法；并在ＰＶＭ并行环境下，对三维高超音速绕流流场实现了多机并行计算，通过负载平衡等方法得到了较高的加速比（在二处理机系统上加速比为１∶８４，在四处理机系统上为３∶４４）。     

Φ50mm口径气体炮炮口制退器的设计与效率分析

针对气体炮发射过程中受后坐力作用影响瞄准精度的问题,基于兵器标准推荐数值确定了炮口制退器的主要尺寸参数及其变动范围,设计了适用于Φ50 mm口径的一级气体炮反作用式炮口制退器,理论分析了制退器的制退效率,获得了高压气室压力、弹体质量、侧孔倾角以及侧孔直径等试验参数对制退效率的影响,并通过有限元流场仿真验证了制退效率,该方法和制退器的结构形式为其他应用场景的气体炮炮口制退器的设计提供了参考。     

超声速飞行器流场的并行数值计算及效率分析

基于有限体积方法、TVD差分格式和显式Runge-Kutta迭代方法的框架,针对超声速/高超声速飞行器绕流流场,在超级并行计算机上完成了2~64个CPU并行数值计算工作。通过测试程序在超级计算机上的并行效率,并将并行程序应用于航天飞机绕流流场计算,检验了计算程序进行大规模并行计算的性能。结果表明,在负载平衡的条件下,程序在该超级并行计算机上达到了不同程度的超线性加速比,并行效率最高达到了126%,远远高于微机Cluster并行平台上的结果,适合复杂流场的大规模并行计算。     

三维实验固体火箭冲压发动机燃烧室湍流反应流数值模拟

在对固体火箭冲压发动机的研究中 ,建立合适的数学模型对其内部工作过程进行模拟是冲压发动机设计和试验的重要环节。本文在三维湍流反应流情况下 ,对某实验固体火箭冲压发动机建立了数学模型并进行数值求解 ,得出了一些有用的结论 ,为冲压发动机的设计和实验提供了有用的理论研究基础     

某新型减阻试验台理论仿真

针对现今水洞、风洞减阻试验准备过程繁琐、周期长、费用昂贵等问题,根据同轴旋转粘度计原理.提出了针对脊状结构减阻、涂层减阻等减阻技术的新型试验台.通过选取V型和U型两种脊状结构,对试验台模型进行模拟仿真计算,将具有脊状结构的试验台壁面剪切力值和湍流强度与光滑试验台数据进行比较,得到了较好的减阻效果.从而说明运用本试验台进行减阻试验的方案是切实可行的.     

二维水翼空化流动的数值模拟

基于RANS方程,结合不同湍流模型与空化模型对粘性流中二维NACA66012模型水翼附近流场进行数值计算与模拟.计算得到固定攻角、固定来流速度时升力系数和阻力系数随空化数的变化规律以及固定空化数时升力系数随攻角的变化规律,并对相应流场进行分析.在二维水翼空化流动的准稳态模拟中,重点对比分析了不同湍流模型、空化模型对空化流动模拟结果的影响,并将计算结果与实验数据进行比较.研究结果表明,在二维水翼附近空化流动的数值模拟中,采用k-ω湍流模型和均相流空化模型计算得到的结果与实验数据吻合较好.最后,采用分离涡模型对不稳定的空化流动进行了瞬态模拟,并得到空化流动的变化规律.     

二维激波/湍流边界层干扰的混合LES/RANS模拟

为了降低高雷诺数条件下大涡模拟方法的计算量,将两方程k-ωSST湍流模型与Yoshizawa一方程亚格子模型通过一个衔接函数相结合,构造一种混合大涡/雷诺平均NS方程模拟方法(混合LES/RANS)。使用这种方法及AUSM+-up格式对20°压缩斜坡的马赫2.85流动进行模拟,并考察了固定入口和在入口添加白噪声两种湍流入口边界条件对于结果的影响。模拟结果再现了边界层的分离、再附以及分离激波等现象,计算得到的分离区要显著大于试验结果,对于这种混合模拟方法的缺点进行了分析,并提出了可能的改进方法。     

一种基于低扩散通量分裂方法的WENO格式研究

将低扩散通量分裂格式(LDFSS)和加权基本无振荡格式(WENO)相结合,构造出一种混合格式,其中WENO格式用于物理量重构,而LDFSS用于通量分裂。采用这种格式对Riemann问题、钝头体高超声速无粘绕流流场进行了计算,并对超声速平板湍流边界层进行了混合LES/RANS模拟,计算结果表明:相对于采用Lax-Fridrichs分裂的WENO格式来说,这种混合格式对于激波和接触间断的分辨率更高,并且在标量保正性方面更优,收敛性更好;而相对于采用带有Minmod限制器的MUSCL方法进行物理量高阶重构的LDFSS格式来说,这种混合格式在混合模拟的计算中能够更好地反映湍流流场的脉动特性,计算得到的湍流速度脉动量的统计值更加准确。     

脊状表面湍流边界层高阶统计量特性

通过对不同风速下V型脊状表面及光洁平板表面湍流边界层内偏斜系数和平坦系数的测试,对比分析了脊状表面边界层内偏斜系数和平坦系数的分布规律。试验在一小型专用风洞中开展,流场测试中使用恒温式IFA300智能型流动分析仪,测试模型则采用有机玻璃材质的矩形平板结构,且试验中模型表面脊状结构的方向与流向一致。最终研究结果表明,脊状表面湍流边界层内偏斜系数和平坦系数分布规律与平板表面基本一致,但脊状结构的存在降低了距壁面无量纲高度y′10区域(包括整个粘性底层和过渡区的一部分)的偏斜系数和平坦系数,而对边界层中过渡区以外区域则影响不明显。由此可以推断,脊状结构主要影响边界层流场的近壁区。