高超声速圆球模型飞行流场的数值模拟和实验验证

对高超声速圆球模型飞行流场进行数值模拟 ,分别采用空气完全气体模型、平衡气体模型以及热化学非平衡 1 1组元气体模型求解非定常轴对称N -S方程组。使用有限差分时间相关法捕捉激波 ,得到了定常流场的解。差分方程隐式部分采用了LU -SGS方法以避免矩阵运算 ,对化学反应和振动能量源项采用预处理矩阵以解决刚性问题。由计算结果处理得到的阴影图和干涉条纹图与再入物理弹道靶实验照片进行了对比分析 ,验证了实验中圆球飞行流场大部分区域接近于平衡状态。     

超声速化学反应流场迎风格式数值模拟及并行计算

基于有限体积迎风格式对超声速燃烧流场进行了的数值模拟.由于超声速燃烧流场绕流的复杂性,要求对多组分Euler/N-S方程求解的数值模拟方法应具有较高的计算精度及效率.本文引用辅助点方法建立了具有空间二阶精度的van Leer迎风矢通量分裂格式,并应用于超声速燃烧流场绕流的数值模拟.化学反应为氢气/空气十反应模型,采用考虑了化学反应特征时间的当地时间步长显式Runge-Kutta时间推进格式.对钝头体模型爆轰现象、后向台阶氢气喷射及二维内外流超声速燃烧流场模型进行了区域分裂技术的并行计算.计算结果与参考文献作了对比,得到了满意的结果.     

激波绕楔形物体的流场分析

使用数值计算方法和数值干涉技术并结合实验干涉图片,分析了激波过拐角绕楔形物体运动的流场.并采用数值彩色干涉技术,将数值计算结果转换为彩色条纹编码的数值干涉图,并与实验干涉图比较,可以辨识过激波发生跳跃后的干涉条纹的位移数,从而确定实验干涉照片的密度场.数值计算采用显式二阶迎风TVD格式.采用数值干涉技术将计算密度场转化成数值干涉图与实验干涉图进行了比较,检验了数值结果的可靠性.数值计算结果显示了包括激波、接触间断和旋涡在内的详细的流场结构,与实验相比,更为细致地刻画了流场演变过程的细节.通过数值计算结果观察到激波传播到旋涡中心附近并没有发生明显耗散,旋涡上半部的激波得到减速,下半部分的激波得到加速.     

AUFS格式在多组分反应流场模拟中的应用

对AUFS格式进行了讨论,并将其应用于高超声速多组分化学反应流场的数值模拟.超声速燃烧流场复杂性使得发展基于非结构网格的高精度逆风格式数值求解方法成为必然.本文中控制方程为轴对称多组分Eu-ler方程,空间分裂采用基于非结构网格的有限体积逆风格式AUFS方法,时间推进为考虑了化学反应特征时间效应的显式5步Runge-Kutta方法.化学反应模型采用经典的H2/Air的7组分8反应模型.对钝头体激波诱导亚爆轰及超爆轰流场进行了数值模拟,得到了与实验及文献相符合的计算结果.     

爆炸波与超声速飞行物体相互作用研究

本文数值模拟了超声速飞行物体与爆炸波相遇后产生的两波干扰流动，采用了激波捕捉、有限差分法和LU-TVD格式。同得到的实验结果对比，计算结果反映出了两波作用后的复杂流动，并且可清楚地看到飞行体头部产生的弱接触面，与实验符合较好，同时进一步模拟了飞行体以-5°攻角飞行时与爆炸波相遇后的流动。     

辐射及化学非平衡流耦合场计算方法

把一种新的三维热辐射计算方法应用到高超声速非平衡流场的数值模拟过程中:在非结构网格上对高超声速化学非平衡流场的数值模拟过程中,耦合了辐射输运方程来求解辐射热流。控制方程为含有化学反应源项和辐射源项的三维N-S方程,数值离散格式采用了Jam eson有限体积法,辐射输运方程采用有限体积法求解。化学反应模型为7组元7反应模型,并采用了化学反应特征时间步长。本文分别对RAM C-Ⅱ和MU SES-C模型进行了计算,计算结果与参考文献的结果吻合。     

超声速/高超声速风洞试验段结构形式对流场性能的影响研究

通过求解轴对称 N-S 方程,对Φ1 m 高超声速风洞马赫数3和6状态下的流场进行了模拟,计算结果与试验数据基本一致,验证了所用数值方法的可信性。在此基础上,对比研究了马赫数3和6状态下采用闭口等直圆截面和开口自由射流两种试验段结构形式的超声速/高超声速风洞在起动条件下的稳态流场性能。结果表明:采用闭口等直圆截面试验段和开口自由射流试验段的流场均匀区内速度场性能指标均满足相关标准要求;马赫数3喷管采用闭口试验段时,沿风洞轴向-300mm~900mm 截面范围内的流场均匀区直径均保持在Φ882mm 以上,均匀区面积较开口试验段增加了约31.57%;马赫数6喷管采用闭口试验段时,均匀区面积比开口试验段仅增加了约8.24%,流场品质略为提高。超声速条件下,闭口试验段的流场均匀区增加明显;但在高超声速条件下,闭口试验段的流场均匀区增加比较有限。     

运动激波与气泡串相互作用的多介质数值模拟

通过对激波和流体界面相互作用而诱导的大变形界面演化的数值模拟,验证了Level set 方法精确模拟多个流体界面的有效性.采用间断有限元Galerkin方法求解欧拉方程得到流场解,采用5阶WENO格式求解Level set方程追踪多流体界面,界面附近的边界条件由虚拟流体方法处理.对运动激波和两个气泡相互作用过程进行了数值模拟,得到了不同时刻的压力和密度等值线分布,并分析了计算域中两个气泡同是氦气泡,以及一个是氦气泡,一个是R22气泡情况下的计算结果.计算结果表明:利用多界面Level set方程可高质量地捕捉多个流体界面,处理3种多介质流场数值模拟问题.     

三维典型超声速燃烧室流场数值模拟

超声速燃烧是实现以吸气式冲压发动机为动力进行高超声速飞行的关键技术之一.本文对三维超声速燃烧室流场进行了基于有限体积方法的并行数值模拟,还对JPNal三维典型超声速燃烧室内的激波诱导燃烧流场进行了数值模拟.控制方程为三维多组元带燃烧反应模型的Euler方程,空间离散采用二阶精度VanLeer迎风通量分裂格式,时间推进为考虑了化学反应特征时间的5步Runge-Kutta方法,燃烧反应为氢气/空气十反应模型.为提高计算效率,采用了基于区域分裂技术的并行计算.最后对比分析了单、双气态氢气喷流对燃烧室内流场结构及燃烧特性的影响.     

高升力下二维构型阻力准确测量研究

翼型风洞试验阻力测量常使用尾迹流场测量积分求取阻力的方法,但各积分公式均建立在一定的假设基础上,有一定适用范围.在多段翼型流场N-S方程数值模拟和风洞试验的基础上,研究高升力情况下低速风洞阻力精确测量技术.通过N-S方程数值模拟求解多段翼型绕流场,分析尾迹流场的特点和常规风洞试验阻力计算公式推导时所作假设,提出新的更为准确的型阻计算公式;利用多段翼型绕流的数值模拟结果,积分表面压力和摩擦力求得翼型的气动特性,并利用计算得到的尾迹流场信息按照常规和新提出的风洞试验型阻计算公式计算阻力,将三者进行比较,检验提出的新型阻计算公式的准确性;通过风洞试验检验数值模拟得到的流场特点和新型阻计算公式.研究表明:在高升力条件下,传统型阻计算公式有很大的局限性,必须进行改进;提出的考虑尾迹区流动特点的新型阻计算公式能够得到更准确的阻力值.     

弹尾超声速轴对称喷流的正格式数值模拟(英文)

采用二阶正格式方法对超声速轴对称喷流流动进行数值模拟。将二维守恒方程的正格式方法发展到轴对称Euler方程组的求解,并对导弹尾部超声速伴随射流进行了数值计算。数值结果与实验照片所反映的流动特征吻合较好,与三阶精度间断有限元方法计算结果吻合,相比于二阶高分辨率TVD格式的计算结果,正格式方法的计算结果在间断点处具有较大的梯度变化。这表明该方法对激波具有较强的捕捉能力,在间断处不会出现数值振荡,对弹尾声音速伴随射流的数值模拟真实、有效。     

旋转冲压发动机冲压转子盘腔冷态流场数值模拟

旋转冲压发动机是一种没有叶片和活塞的基于冲压压缩技术的新概念发动机,冲压转子是其惟一的核心转动构件。通过应用商用CFD(NUM ECA)软件,基于三维定常N-S方程和B-L湍流模型,对旋转冲压发动机冲压转子盘腔冷态流场进行了数值模拟,着重分析了其旋转冲压压缩效果和流场结构。分析结果表明,冲压转子能有效实现超声冲压压缩,进气道激波结构类似于超声进气道,但略有差异;盘腔内流场结构是带有周边小旋涡回流和低压中心区的大旋涡结构,这种流场结构有利于燃油供应和燃烧组织。     

嵌套网格粘性流动数值模拟用于风洞洞壁干扰研究

运用嵌套网格技术和Navier Stokes数值模拟对机翼半模和翼身组合体试验时风洞的四壁干扰进行综合模拟、评估和修正。计算格式在空间上采用中心有限体积离散 ,在时间上采用多步Runge Kutta时间步长格式进行积分。计算结果证明了该方法的可行性和优越性。     

激波运动与反射流场的ENO算法

在ＴＶＤ格式的基础上，以修正限制函数的方法，将通常的ＴＶＤ格式改造为性能更为优越的ＥＮＯ格式，并利用这种有限差分格式求解欧拉方程和薄层Ｎ－Ｓ方程。分别数值模拟了激波管中运动激波与三种楔形物体相互作用流场；给出了激波在不同物体上的反射及波涡干扰过程中非定常流场的变化情况；算出的波系结构与在ＣＡＲＤＣ小型激波管中做的相关实验结果一致。     

不同进口马赫数下超声速燃烧流场特性研究

对有无楔板超燃冲压发动机模型内横向氢气喷流超声速燃烧流场进行了数值模拟,分析了进口马赫数对超声速燃烧流场特性及特征参数分布的影响特性。采用有限体积法求解多组元Navier-Stokes(N-S)方程,对不同进口马赫数下的燃烧流场进行了数值模拟,细致对比了流场激波结构、喷流穿透深度、燃烧阵面,燃烧效率及总压恢复系数等参数随进口条件的变化特征。结果表明:无论是否存在楔板结构,喷口后流场压强均随着进口马赫数的增加而减小,并且随进口马赫数的增加,氢气喷流穿透深度减小,楔板对喷流穿透深度基本无影响。较无楔板结构而言,设置楔板结构可以缓解燃烧室内流场对马赫数变化的敏感度,使燃烧更为稳定。在同一进口马赫数条件下,楔板布局有明显的促燃作用及激波点火效果,在一定程度上可增加此类发动机工作的马赫数范围,但以总压恢复系数略微降低为代价。     

NUMERICAL SIMULATION OF THREE DIMENSIONAL INTERACTING TURBULENCE FLOW FIELD OVER PROJECTILE WITH LATERAL JETS

INTRODUCTIONIn recent years,the study of lateral jets hasaroused more interest of aerodynamicist all overthe world.On the one hand,the interacting flowfield caused by lateral jet includes complicatedshock wave patterns and has important values tostuding flow phenomena.It has been not veryclear to this day and it is very difficult for nume-rical simulation.On the other hand,lateral jethas important applications for the field of aero-nautics,astronautics and weapons:for example,lateral jet …     

含化学反应流动的TVD格式及其在燃气射流中的应用

以Harten标准TVD格式为基础，结合固体火箭燃气射流的特点，以数学方法系统地推导了适合于高温、高压和高速流体流动的数值格式，给出了压力偏导数的合理计算公式；利用化学动力学知识，对火箭燃气射流流场中存在的多组分、含有限速率的化学反应问题进行了论述，阐明了化学反应质量源项的求解方法。以12组分9反应方程模型为例，利用编制的计算程序，对某火箭燃气自由射流流场进行了模拟。通过对结果的分析，肯定了数值格式的正确性。     

电弧加热器超声速湍流平板烧蚀流场变化研究

电弧加热器超声速湍流平板烧蚀试验技术是研究防热材料烧蚀特性的重要手段。为研究超声速湍流平板烧蚀过程中流场变化情况，采用数值求解二维N-S方程的方法进行试验流场模拟。从模拟结果看，未烧蚀模型外形流场模拟得到的模型表面参数结果与试验结果吻合很好。然后对烧蚀过程中的模型外形进行了流场模拟，并与试验流场进行对比，根据模拟结果分析了试验过程中模型表面压力和热流密度分布的变化。根据分析可知，如果平板模型烧蚀量最大的位置在初始高热流区内，可以采用该烧蚀量计算烧蚀速率。