NND格式在航天飞机头部段N—S方程求解中的应用

本文采用实质上为二阶精度的无波动、无自由参数的耗散差分格式(NND格式)求解了非定常的二维和三维完全N-S方程,对类似于航天飞机头部外形的超声速粘性气体绕流流场进行了模拟。在计算中,使用了代数网格生成技术,对脱体激波采用了装配法。作为验证算例的二维圆柱绕流计算结果同文献[4]的数据进行了比较。二维凹陷外形绕流的计算结果同无粘流结果进行了对比,由于粘性流动出现了分离,二者的壁面压力是有差异的。最后,对三维外形的绕流也进行了模拟,计算结果提供了流场的细节。     

高升力系统外形的数值模拟计算

分析和讨论了高升力系统外形无黏和黏性流动的物理特征,表明了流动现象的复杂性和求解的艰难性。深入研究了二维和三维高升力外形空气动力数值模拟计算的雷诺平均N-S(RANS)方程方法、计算结果、以及它们和实验数据的比较。二维结果表明压强分布、速度型等计算值和实验值吻合很好,但最大升力系数两者仍有差异。三维计算最大升力系数的结果与二维类似。还分析了影响计算结果的各种因素,包括计算网格、湍流模型、外形的几何模拟的逼真度等。特别指出了转捩模型研究对提高计算准确度的重要性。文中解释了研究高升力系统外型的数值模拟方法,以及同时发展能估算C_Lmax的工程半经验方法的重要性和迫切性。     

一种新的分区重叠洞点搜索方法—感染免疫法

本文分析了分区重叠计算中洞点搜索的传统模式及春应用上的缺陷，针对凹形重叠区难以精确寻点的问题，提出了一种新的洞点搜索方法：感染免疫法。文中给出了二维和三维情况下的搜索算例作为验证，实际应用表明，方法无论是寻点能力还是在通用性方面都优越于传统模式。     

微尺度拉伐尔喷管冷态流场及其推进性能

为了研究微尺度超声速喷管的推进性能,采用无滑移的连续介质模型,对二维和三维微尺度拉伐尔喷管的冷态流场进行数值模拟。对不同算例的模拟结果进行比较,给出二维微喷管对几何结构和工作条件的要求,结果表明,在绝大部分工作条件下,20°扩张角是较好的外形设计选择;对三维和二维模拟结果进行比较,给出二维微喷管对蚀刻深度的要求,结果表明,蚀刻深度和喷管喉部宽度的比值超过10时,微喷管具备很好的二维特性。     

三维非结构贴体直角网格的N-S方程数值模拟

基于八叉树数据结构,采用直接分层光顺的方法,解决了非结构直角网格在物面不正交贴体的难点,发展了任意复杂外形的三维非结构贴体直角网格,并且可以根据流场特征进行自适应处理。在此基础上采用最小二乘法对N-S方程进行了MUSCL有限体积法的非结构网格的粘性求解,提出了带约束边界条件的最小二乘法的处理办法。计算表明,三维非结构贴体直角网格的粘性求解,可以适应各种高超声速的复杂外形,计算结果与试验吻合较好。     

多芯片组件的热管理

较详细地介绍了二维和三维多芯片组件的几种冷印方法,给出了导热计算公式;并提出了热设计的若干问题。     

三维组件布局建模与优化设计

 从工程应用背景出发,在研究现有组件布局问题的基础上,将Sphere-trees概念引入三维装填布局优化设计问题。采用有限球族进行非规则任意形状三维组件外形轮廓的统一描述,而非直接采用组件的实际外形轮廓进行装填布局设计,建立了任意形状组件装填优化设计干涉检查的系统性计算方法。通过几个典型算例,验证了该方法与序列二次规划法相结合解决有限空间内组件布局优化问题的适用性。     

基于阻力气动特性计算的飞艇艇身外形研究

采用SIMPLE算法、S-A湍流模型求解了三维不可压雷诺平均N-S方程,数值模拟了绕飞艇的低速不可压粘性流动,计算了不同艇身外形的阻力特性.本文方法计算结果与实验结果或其他文献的计算结果符合良好,可以用于飞艇的气动特性计算.通过对6个不同外形的相同体积飞艇阻力特性计算,得到了最佳艇身外形,表明本文方法可以用于艇身外形的选型设计.     

几种典型迎风格式的分析与比较

重点分析了AUSM 格式的构造方法,同时将其与Van Leer, Roe格式等构造方法进行比较和分析,通过求解一维激波管问题,以及二维和三维跨音速流场的计算比较了AUSM 格式,VanLeer格式,Roe格式和中心格式的数值模拟结果和网格适应性,以及在经典四步Runge-Kutta显式时间推进和LU-SGS隐式时间推进下上述四种格式的计算效率.获取了对于基础研究和工程应用很有意义的结论.     

复杂外形气动力计算的阻力修正方法

基于气动力计算的位流线化理论,发展了一种亚音速附面层阻力修正方法。采用联合流场求解方法获得复杂外形的表面压力分布,在此基础上,将绕复杂外形物体沿流向面元划分为片条,在片条上,把三维附面层转换为二维附面层,对二维附面层进行数值求解,得到附面层粘性阻力修正。算例表明,该方法计算量小,与实验和N-S数值解结果吻合较好。为复杂气动外形飞行器初期设计、方案选型及评估提供了较为理想的工程设计手段。     

向量化多重网格的Mac Cormack显格式解轴对称底部干扰流场

本文将多重网格法与MacCormack显格式的算法相结合,计算了轴对称底部流场中斜向喷流与超声速外流粘性干扰流场,提高了计算效率。流场的分析表明分辨率较好,计算合理。在YH-1上将程序向量化,大大节省了计算时间。     

高阶精度非线性格式WCNS-E-5在二维流动中的应用研究

本文采用具有5阶精度的加权紧致非线性显式格式(WCNS-E-5)对定常与非定常二维流动进行数值模拟,研究表明该格式对各类间断有很好的分辨捕捉能力,而且对强间断如激波的计算,即使在高马赫数与高雷诺数条件下它仍具有很好的收敛性与可靠的计算结果.此外,WCNS-E-5在粗网格条件下也体现出优越性.类如WCNS-E-5的高精度激波捕捉方法将为以后开展湍流数值模拟工作提供坚实的技术保证.     

叶轮机械中三维周期性非定常流场测量

本文改进了三维非定常速度场的单斜丝测量技术,发展了三维非定常压力场的试验-数值确定方法。采用上述测量方法,借助于位移结构、高速数据采集系统和锁相采样平均技术,成功地测出了叶轮机转子出口三维周期性非定常流场。      

任意截面细长弯曲管道内的可压缩势流解

利用摄动法给出了任意截面细长弯曲管道内的可压缩定常势流的解。引入了一个表示管道细长程度的小参数ε,并将流场的解对小参数进行渐近展开。文中给出了对密度和轴向速度具有两项展开式、对横向速度具有一项展开式的结果。理论分析表明,在细长管道的情况下,问题的解可降低一维。同时,在上述近似条件下,求横向速度分量相当于求解一个平面Laplace方程的Neumann边值问题。     

实用外形亚跨超流场数值模拟研究

本文给出了实用外形亚跨超流场数值模拟若干新研究成果，包括实用的表面及空间网络生成技术，高效快速、稳定的流场求解方法，即ＬＵ分解方法、利用这套方法研制了实用程序，并对前机身、全机、前机身／进气道等实际问题，进行了较大规模的试算。试算结果与实验数据进行了比较，符合良好，为进一步开发工程上实用的软件打下了良好的基础。     

电弧风洞内球锥全目标绕流场计算分析

针对电弧风洞高焓气体来流条件,通过数值求解三维化学非平衡Navier-Stokes方程,分析无烧蚀球锥全目标绕流场真实气体效应.考虑空气电离化学反应系统,包含7个化学组元和6个反应方程.差分算法是隐式NND格式.依据实验来流条件,计算了不同壁温及不同壁面催化特性下的流场结构和电子数密度,观察了来流非均匀性的影响.结果表明.尾迹区对模拟条件较敏感,其电子数密度随壁温降低而升高,完全催化壁的电子数密度约比非催化壁低近一个量级,实验状态的计算应采用风洞实际来流条件,     

基于小波奇异分析的流场计算方法及应用

首次将三维小波以及二维小波奇异性分析的思想引进三维以及二维复杂流场的数值计算,发展了一种高效率、高精度、高分辨率的方法.该算法的核心是获取流场中物理量在不同空间点的Hlder指数α,而该指数α的获取又依赖于小波变换以及高维(即二维或三维)小波分析技术.在三维与二维欧氏空间中,为进行小波多分辨分析,需要在尺度空间与小波空间分别引进尺度基与小波基.对二维问题,尺度基与小波基的基底要由1个尺度函数与3个小波函数组成,而三维时要由1个尺度函数和7个小波函数组成.借助于小波奇异分析所找到的流场中奇异点区域与光滑区域,便可分别选用高分辨率、高精度的WENO(weighted essential non-oscillatory)格式与高精度中心差分格式进行流场的离散求解.一系列二维(即,①二维前台阶问题的Euler流;②二维双马赫反射的Euler流; ③著名的二维Riemann问题;④跨声速RAE2822翼型二维绕流; ⑤跨声速VKI-LS 59二维涡轮叶栅绕流)与三维(即,⑥跨声速轴流压气机转子NASA(National Aeronautics and Space Administration) Rotor 37三维黏性绕流;⑦跨声速风扇转子NASA Rotor 67三维黏性绕流)算例表明:该方法的计算效率比传统的WENO格式有较大的提高,大部分典型算例能够提高3~5倍,而且可以获取复杂流场中高分辨率的激波结构以及涡系分布,可以得到与有关实验数据较为吻合的数值结果.      

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

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

Detonation waves induced by a confined wedge

An auto-ignition detonation phenomenon can be initiated by a wedge confined in a channel. A new, self-sustaining, normal detonation wave engine concept is proposed. The detonation processes were numerically modeled with a simplified two-dimensional wedged channel flow that was deemed to emulate a real three-dimensional configuration. The results showed that within certain ranges of incoming flow Mach number or wedge angle, detonation could be self ignited in the channel. The study furthermore investigated the detonation waves based on three different detonation initiation positions. Different configurations of the detonation waves were observed and analyzed. The performance of the different detonation wave configurations was estimated and compared.     

二维超／高超声速进气道流场数值模拟

对超／高超声速三级压缩进气道流场进行了数值模拟,来流马赫数为4、6,进气道内流动为层流状态,根据二维Navier-Stokes方程,采用二阶精度Roe格式进行离散。按照流场特点,合理地设计网格分布及调整不同黏性范围的熵修正,防止了壁面附近过大的数值耗散,使计算结果更加合理。在进气道模型的各级压缩折转角处,获得了清晰的激波结构,在进气道内部的各种波系的相交、反射和激波诱导的边界层分离等现象都得到合理的描述。计算得到的压力分布,在各级压缩斜板上同简单波理论结果十分接近。用本文方法计算了另一个二级压缩进气道,沿上、下壁面的压力分布与试验比较符合得较好。