数值模拟有外挂物的空腔流动

数值模拟了有外挂物存在的开式空腔的超音速流动。数值计算采用全隐式无分裂方法,湍流模型采用修正的B-L代数模型。研究了有外挂物存在的空腔流动的剪切层自持性振荡,并与无外挂物存在时的现象进行类比。压强与时间的变化曲线、空腔底部流向的静压系数分布以及空腔后壁面法向的静压系数分布表明,空腔内的压力脉动与外挂物所处位置是有关的。对于Rossiter振荡模态,数值计算结果能够与实验结果相吻合。     

基于RANS/LES方法的超声速底部流场数值模拟

分别采用基于两方程k-ω剪切应力输运(SST)湍流模型的延迟DES(DDES)、更改的DDES(MDDES)和改进的DDES(IDDES)方法,并引入可压缩修正,结合三阶MUSCL-Roe和五阶WENO-Roe两种空间离散格式,针对超声速底部的复杂流动现象,开展了数值模拟研究。计算结果表明本文方法能够捕捉到超声速底部流动中丰富的湍流结构,通过分析计算结果对超声速底部的流动机理有了进一步的认识,为下一步的超声速底部流动减阻改进和雷诺平均NavierStokes/大涡模拟(RANS/LES)方法在非定常高可压缩性流动中的应用提供了参考。通过对比分析不同空间离散格式的计算结果研究了数值耗散对计算的影响,五阶WENO-Roe格式的计算结果与实验结果吻合良好;对不同RANS/LES混合方法的计算结果进行了对比分析,结果表明IDDES方法在近壁区的表现优于DDES和MDDES方法。     

二维高亚声速空腔流激振荡的数值模拟研究

采用基于P .L .Roe的近似Riemann解的修正Osher&Chakravarthy (MOC)三阶TVD有限差分格式 ,数值求解二维雷诺平均全Navier Stokes方程 ,并用Cebeci Smith代数湍流模型 (对腔内区域作修正 )来模拟湍流效应 ,时间方向的积分采用四阶Runge Kutta方法 ,通过对GAMM超音速前台阶绕流的计算 ,验证了格式及程序的有效性 ,对高亚声速来流下的空腔流动作了数值模拟研究 ,取得了较好的效果     

非平衡流计算方法及其模拟激波诱导振荡燃烧

给出一种数值模拟高超声速化学非平衡流动的计算方法，采用流动项与化学反应生成源项解耦处理，即可以有效地解决非平衡流动计算中遇到的刚性问题，还具有公式推导简单和应用方便的优点。从薄层近似N－S方程出发，采用ENO差分格式，数值模拟了超音速H2／Air预混合气体钝体绕流场，较好分辨出流场中激波诱导振荡燃烧非定常流动现象，表明计算方法具有较高的精度。对计算方法的效率和网格数的影响也进行研究。     

后向台阶横向喷流干扰流场的数值研究

 用有限差分法求解层流全N-S方程,对超音速后向台阶绕流与台阶后横向喷流干扰流场进行了数值模拟。不带喷流的超音速后向台阶绕流数值模拟结果与实验及其他计算结果吻合较好。最后给出横向喷流压力比p_j/p_∞=3.033～7.033和喷口宽度w=0.1h～0.2h(h为台阶高度)时的数值模拟结果。得到了合理的流场结构;并可看到喷流压力比和喷口宽度对流场特性有显著影响,在计算参数范围内,底部压力与系数Q(=wp_j/hp_∞)几乎成线性关系。     

过渡区球锥气动力特性的桥函数估算与DSMC仿真研究

采用DSMC(DirectSimulationMonte-Carlo)方法模拟球锥体过渡区流动。论文给出了不同Knudsen数球锥绕流流场的温度与压力分布图,以及在不同Knudsen数、不同马赫数条件下球锥阻力系数和底部阻力系数的分布,并将阻力系数模拟结果与桥函数一体化估算结果进行了比较分析。     

用于粘性可压缩流动数值计算的SIMPLE方法

对SIMPLE方法加以推广,使之适用于可压缩粘性流动的数值计算。采用非交错网格技术,推出三维任意曲线坐标系下可压缩形式SIMPLE方法的计算方程。对亚音速、跨音速及超音速等5个流场进行了数值计算并与有关文献及实验数据进行比较。      

空腔流动特性数值模拟研究

本文采用数值方法模拟研究了空腔的长深比（L／D）对流动特性的影响。控制方程采用Navier-Stokes方程，计算方法采用有限体积法，湍流模型采用Baldwin-Lomax模型。通过分析三种典型空腔的流线图谱、压力云图、底部压力分布曲线，得到空腔流动特性与其几何参数和来流参数间的相应关系。从数值模拟的结果来看开式空腔流动的空腔内压力分布均匀，流场结构相对稳定，有利于内埋式武器安全平稳地分离。     

凹槽非定常超音速流动的研究

用SST的两方程紊流模型,采用时间和空间上二阶的方法,求解雷诺平均N-S方程,得到开放型凹槽的流动特性的变化。数值模拟了飞行马赫数为5时的测试结果,得到不同时刻流场参数的变化。凹槽内流体质量的增加和减少,有助于油气的混合和燃烧,对于维持燃料的超音速燃烧是非常关键的,对超燃冲压发动机燃烧室的设计具有重要的价值。     

基于RANS/LES混合方法的分离流动模拟

飞行器在大迎角、快速俯仰机动时,流场中含有大尺度、非定常的涡结构,传统雷诺平均Navier-Stokes （RANS）模型不能准确模拟流场结构,根据国际上相关研究的发展趋势,需要采用混合RANS/大涡模拟（LES）模型来对复杂分离流动进行准确模拟。本文对基于分区混合与湍流尺度混合的双重RANS/LES混合计算模型进行发展与应用。通过典型简化模型的静、动态湍流大分离流动,测试和验证所采用的脱体涡模拟（DES）类方法,重点研究改进的延迟DES （IDDES）模型在动态问题应用中的正确性和有效性,并对所采用的数值模拟方法和相应的计算软件的可靠性、鲁棒性以及精度进行了考核验证。典型算例包括超声速圆柱底部流动、跨声速方腔流动、NACA0015机翼深失速分离涡模拟等。计算表明：发展的IDDES类混合计算模型可有效解决对数层不匹配的问题;对于定态非定常分离流动,DES、DDES、IDDES等模型计算结果差别不大,随着流动的非定常特性增强,IDDES模型的优势逐渐显现;对于动态非定常分离流动,则需要采用IDDES类模型。     

液晶流动可视化方法研究拟似冲击波的内部超声速流

在一个压力一真空超声速风洞中,剪切应力敏感液晶流动可视化技术被应用来研究方管内马赫数2拟似冲击波（pseudo-shockwave）的超声速流动。它主要提供关于整个流动的定性信息,诸如湍流边界层分离、再附着位置以及流动的维数等。而且液晶也反映了表面流线,分离区内的涡流和管道流动的角效应。使用两种不同黏度的液晶分别进行实验,分析黏度对结果的影响。液晶实验的结果与纹影照片所得结果比较吻合,说明了液晶是一个非常有效的流动可视化工具。     

喷流流场旋涡激波干扰的数值模拟研究

本文采用隐式CSCM格式,通过数值求解NS方程,对平板上二维超声速主流中横向喷流流场进行了数值模拟,并与实验和其它计算结果进行了分析比较,结果令人满意。     

超声速流动中横向射流流场的影响参数

采用高精度的Weighted ENO格式，结合两方程湍流模型，准确模拟了二次射流形成的干扰流场，详细地描述了平板上单股射流干扰流场和喷管扩张段二次射流干扰流场中的激波、流动分离和旋涡运动，同实验结果进行了比较。探讨了射流／主流总压比、射流宽度。以及射流与来流夹角对射流穿透深度、分离距离等影响，揭示了二次射流推力短量控制干扰流场的控制机理。     

等离子体气动激励控制超声速边界层分离的实验研究

等离子体气动激励与超声速气流相互作用已成为高速流动控制领域的研究热点。激波与边界层相互作用现象广泛存在于超声速飞行器之中。本文进行了等离子体气动激励控制压缩角区和激波诱导边界层分离的实验,通过流场纹影显示和壁面静压测量,研究等离子体气动激励如何影响激波、激波如何影响边界层特性的科学问题。实验结果表明:施加毫秒量级表面电弧放电能够前移压缩角区的诱导斜激波,使分离区后移,分离区域增加,但激波强度减弱,流场总压增加;施加微秒量级表面电弧放电能够抑制激波诱导边界层分离,使分离区减小,流场总压减小。基于实验结果,认为毫秒量级表面电弧放电激励控制超声速气流的主要机理为放电过程的焦耳热效应;微秒量级表面电弧放电激励控制超声速气流的主要机理为焦耳热效应和冲击波效应共同作用。     

超声速底部喷流干扰流场数值模拟

数值模拟了不同马赫数,不同喷流压比下的轴对称超声速底部喷流干扰流场,采用LU隐式算法进行数值求解并引入了Baldwin-Lomax代数湍流模型.采用分区网格将弹身与底部区域合为一个整体进行计算,得到了清晰的流场结构和弹体表面及底部的压力分布,并与试验结果进行了比较.数值模拟结果表明超声速底部喷流干扰流场结构复杂,有、无喷流时底部流场有很大不同, 喷流对底压分布有明显影响,进而对轴向力系数影响显著.     

风洞分流锥及孔板整流的数值模拟研究

采用数值模拟技术研究了2m×2m超声速风洞设计引导试验大开角扩散段配置一个中心倒锥和两层球冠状孔板的内流场,中心倒锥和两层球状孔板的不同组合共有五套。孔板有几百至上千个开孔,是模拟的难点,发展了一种孔板流动CFD边界条件模型。另一个难点在于风洞管道内流入口和出口边界条件的准确处理。计算表明,分流锥可将流动有效抑制在大开角洞壁附近,防止洞壁附近的扩张分离,但在分流锥底部将产生大尺度的分离涡环结构,经过孔板的整流后,总压有较大损失,但流线趋于平直,可达到预期效果。     

液晶显示方法研究管内马赫数4拟似冲击波流动

在一个压力-真空型超声速风洞中, 剪切应力敏感的液晶表面流动显示技术被用来研究方管内马赫数4拟似冲击波的超声速流动.实验是在日本室兰工业大学的超声速风洞中进行.它提供了关于测试表面流动结构的定性信息,诸如湍流边界层分离和再附, 以及边界层流动的维数等.液晶显示结果与纹影实验照片和油膜方法进行比较,证实了它是一种非常有效的表面流动显示工具.     

等离子体合成射流对钝头激波的控制与减阻

为了验证等离子体合成射流对超声速流动的流动控制和减阻效果,在考虑热完全气体效应的情况下,工程拟合等离子体热力学属性和输运系数,利用能量源项模型对超声速平板和球头等典型流场结构进行了数值模拟,并取得了与实验较为一致的结果。研究结果表明：对于马赫数为2的超声速平板流动,等离子体合成射流能有效干扰边界层的发展,并诱导产生一系列大尺度结构;对于马赫数为3的球头流场,等离子体合成射流能显著改变激波脱体距离与球头的阻力特性;在放电后第1个周期内,合成射流能使球头平均阻力降低6.3%,而在射流峰值情况下使阻力降低32.0%,激波脱体距离增加2倍,实现了激波控制和流动减阻的预期效果。     

超声速湍流流动中平板-楔模型分离、传热和烧蚀试验研究

本研究之目的是为了模拟可控有升力飞行器周围的加热和压力场,研究受热表面的烧蚀现象,以便发展地面的实验能力。在电弧加热器上,利用超声速湍流平板技术,进行控制翼模型分离、传热和烧蚀实验。结果表明,由于翼角前有一条横缝,分离效应减弱,导致产生分离的最小翼角增大。由于分离流动影响,在翼上及其周围压力和热流升高,烧蚀速度大大增加。本文给出了有关的相关公式。     

超声速条件下多体干扰与分离试验研究

采用根据国外公开文献设计的类CAV模型,在0．6m×0．6m 跨/超声速风洞中开展了多体干扰与分离网格测力试验研究,初步获得了典型多体飞行器分离过程中的气动特性变化规律。试验结果表明,载荷模型气动特性受分离位置变化影响非常明显。载荷模型沿轴向分离时,气动力（矩）逐步接近自由流中气动力（矩）值,载荷模型法向位置改变会引发其气动力（矩）值发生更为剧烈的变化。引发这种现象的原因有两个：一是尾迹和头激波的发展改变了不同轴向位置处载荷模型的表面流态,从而影响了其气动特性;二是母机模型底部流动具有明显的非对称膨胀特征,不同法向位置处流速大小和方向差异明显,导致载荷模型气动特性随法向位置变化更为剧烈。