超声速前机身及翼－身组合体的全速势方程数值计算

通过数值求解全速势方程，计算了超声速来流的前机身及翼－身组合体。当流场存在亚声速区时，在此区域内采用中心差分格式，迭代求解，并引入多重网格技术，加快收敛；当流场中某一区域沿某一方向是超声速时，在此区域内采用沿该方向的推进解法。计算结果表明，本文的方法可靠，结果准确，可以向工程应用方面推广。     

给定下游边界的超声速流场逆向求解方法

为了实现超声速流场出口参数的可控设计,借鉴特征线方法在古尔赛特与柯西问题两类边值条件中的应用,提出了一种根据下游边界条件求解其依赖域的逆柯西问题求解方法,并在此基础上发展了一种均匀来流条件下的超声速流场逆向求解方法。对轴对称内收缩流场、二维平面流场、轴对称外锥流场三个典型算例校核表明,逆向求解方法得到的计算结果与正向特征线方法一致,由此验证了该逆向求解方法的可行性,并有效拓展了超声速流场的设计思路。     

尖拱弹身横向喷流数值模拟

本文通过数值求解层流Navier—Stakes方程来模拟超声速流场中横向超声速喷流的干扰流场，采用分块对接网格和“O”型网格技术，精确模拟喷口截面形状并生成高质量的贴体计算网格，对尖拱弹身的横向超声速喷流的干扰流场进行了数值模拟，研究了喷口附近的流场结构、涡系结构、波系结构及流动分离等流场特性，得到喷流流场的干扰效应，计算结果与实验值吻合较好。     

强激光与高超声速球锥流场干扰数值模拟研究

本文采用有限差分技术，数值求解带能量源项与空气11组分化学动力学模型的轴对称粘性流体力学方程组，采用类氢原子理论模型计算空气对激光的吸收系数，建立了模拟强激光与超声速钝体流场干扰的计算软件，得到了高能浙江与超声速球锥流场干扰数值模拟的结果。     

超声速和高超声速横向喷流数值模拟

为了解横向喷流的干扰影响和喷口附近的流场结构,采用分块对接网格和“0”型网格技术,数值求解N—S方程来模拟超声速和高超声速流场中横向喷流的干扰流场。对两种尖拱弹身外形的超声速和高超声速喷流干扰流场进行了数值计算,计算结果与风洞实验数据吻合一致。在此基础上,开展了某型导弹多喷构型干扰流场的数值模拟,得出一些横向喷流数值模拟的结论。     

超声速气流中液滴破碎的大涡模拟

为了研究超声速气流中的液滴破碎过程,采用基于CLSVOF界面追踪方法的大涡模拟方法。该方法中用可压流和不可压流求解器分别计算超声速气流流场和液体流场,并且在分别计算气相或液相时,另一相作为该相计算时的边界条件参与计算,从而实现可压流求解器和不可压求解的融合。通过去散度化液体速度外推解决界面处大液气密度比所带来的动量误差大的问题。利用该方法对来流Ma=1.358的超声速气流中的水滴破碎进行了数值模拟。模拟结果展示了超声速气流中液滴破碎的几何结构变化细节和基本破碎形态,且与实验吻合较好。数值模拟结果揭示了Rayleigh-Taylor不稳定和气动剪切力对液滴破碎的作用过程。破碎过程中,液滴在垂直来流方向上的最大无量纲宽度(Dmax/D0)约为4.5,无量纲破碎时间(tbU∞/D0)约为29.85。     

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

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

超声速化学反应流动的LU-SGS伪时间迭代空间推进求解

介绍了求解超声速多组分有限速率化学反应流动的伪时间迭代(lower-upper symmetric-Gauss-Seidel,LU-SGS)方法.空间推进求解多组分抛物化(parabolized Navier-Stokes,PNS)方程时,在一个推进面上采用修正的LU-SGS方法迭代至收敛,把得到的结果作为初值赋给下一个推进面.沿推进面依次迭代直至求解完整个流场.采用伪时间迭代LU-SGS方法求解化学反应PNS方程,计算结果的准确性和时间迭代求解完全Navier-Stokes(N-S)方程相当,求解效率提高一个数量级.      

基于各向异性非结构网格的超声速流动自适应计算

基于各向异性非结构网格,实现了自适应求解技术,将其应用于超声速楔形体绕流及超声速横向喷流问题的计算研究,对各向异性网格自适应计算过程的收敛性和求解精度进行了对比分析,展示了各向异性网格自适应算法在降低问题求解规模。提高各向异性流场分辨率等方面的优势及存在的问题。超声速楔形体无粘绕流场的自适应计算,清晰地捕捉到了激波附近的流场信息,激波前后的马赫数、压力、密度。温度等物理量与理论分析解吻合很好。超声速横向喷流流场存在激波、分离涡、边界层等流场结构的相互干扰。计算研究表明,单纯基于Hes-sian度量张量的各向异性网格生成及自适应算法不能有效模拟边界层内的流动情况,是将来需要进一步开展的研究挑战。     

氢/碳氢燃料超声速燃烧的数值模拟

采用隐式上、下三角分解（Lower-Upper Decomposition）的算法，结合对组分连续方程中化学反应源项的点隐式处理，求解了多组分的N－S方程组，得到了二维和三维条件下的超声速掺混及燃烧的数值模拟结果，通过与水平喷射氢气及预燃煤油等算例的验证，计算结果与实验数据基本符合。在此基础上进行了垂直喷射氢气的三维超声速燃烧的数值模拟，计算结果显示了氢气在超声速空气中掺混、燃烧的过程，该程序可进一步用来模拟超燃冲压发动机燃烧室内的复杂流场。     

带座舱飞船高超声速粘性绕流的分离形态及流场结构分析

通过求解薄层近似的Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ方程，数值模拟了带座舱飞船高超声速粘性绕流，首次系统地给出了攻角为５°、１０°、２０°、３０°和４０°情况下的物面分离形态、横戴面流线形态、纵剖面流线形态以及旋涡沿涡轴的演变形态。分析和比较表明，本文采用数值可视化方法给出的流场拓扑结构与相应的拓扑规律完全吻合。在飞船后体流场横截面流线中，本文发现了稳定的极限环流态，该极限环涡心处沿涡轴方向的速度是超声速     

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

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

非线性壁板颤振分析

利用一种改进的计算流体力学与计算结构动力学(CFD/CSD)耦合方法研究了由气动和结构几何非线性引起的壁板颤振问题。在非定常气动力计算中,考虑了通量分裂格式、隐式时间推进方法和几何守恒律;二维和三维壁板的结构几何非线性建模则采用了有限元的协同旋转理论,并利用一种近似能量守恒算法求解结构的非线性响应。流场和结构求解器采用二阶松耦合方法联立求解,并将其应用于壁板在超声速、跨声速和亚声速的颤振计算中。当考虑结构几何非线性和气动非线性时,出现了典型的极限环振荡现象,并对颤振边界和极限环振荡幅度进行了对比分析。     

三维超声速流动的压力反问题

为了进一步探索三维超声速流道的设计方法,采用一种预设壁面压力分布计算壁面型线的思想,并结合双特征线方法提出一种全三维超声速流动压力反问题的求解方法。在三维超声速流场设计中,可直接根据来流条件和壁面压力分布求解壁面的三维坐标,通过空间步进的方式,使得解在一系列解平面上推进,从而使得所设计的型面与预设的壁面压力分布相容。通过Prandtl-Meyer膨胀波的理论解验证了该格式的设计精度。根据预设的压力分布,设计了圆形和椭圆形入口的三维超声速喷管,并将设计方法与数值模拟进行对比验证。验证结果表明:所设计的流场与CFD计算得到的等值线符合得较好,因此基于双特征线的压力反问题求解方法具备三维超声速气动设计的能力,并具有纯三维、高精度、壁面压力分布可控的优势,对未来高超声速气动设计应用将起到重要的支撑作用。     

激波对激光束远场强度分布影响的数值研究

研究了半径0．1m、焦距10^6m、波长分别为10．6μm、1．315μm的激光束通过具有激波结构的超声速绕楔流场后的远场强度分布。得到了沿光束传播方向积分的流场的光程差分布以及激光束在穿过超声速绕楔流场并在无扰来流中传播1km后光束的强度分布。结果表明：强的激波结构使激光束远场强度分布发生了向光程增大方向的明显偏折，对于波长不同的激光束得到光束通过流场后的远场强度分布不同。     

裙后向台阶超声速绕流特性数值研究

本文研究了在中部带有裙体后向台阶的轴对称体的超声速绕流特性。采用空间推进法数值模拟超声速无粘流场，对于轴对称后向台阶的分离区采用Ｃｈａｐｍａｎ－Ｋｏｒｓｔ理论模型进行处理，计算结果给出了与实验值吻合的物面压力分布以及超声速流场的流动图画。     

外压式超声速进气道流场的数值模拟

采用CFD方法对某型外压式超声速进气道流场进行了二维数值计算,叙述了网格生成、边界条件确定方法,控制方程求解收敛技巧以及计算结果的分析.数值计算结果与理论相符较好.     

三维超声速压力反问题的特征线求解技术

为了进一步探索三维超声速气动反问题的求解方法,基于双特征线理论及Butler解法,研究了三维超声速流线压力反问题的适定性。为了确保解的唯一性,提出在限定壁面膨胀/压缩方向时存在壁面压力与三维坐标的一一映射关系。基于该映射关系,提出了三维压力反问题的双特征线求解技术(iMOC-3D求解器)。采用Prandtl-Meyer膨胀波、Busemann进气道的理论解,对iMOC-3D求解方法的膨胀、压缩过程进行了精度评估,误差均为1×10~(-4)量级。为了进一步验证设计方法的可靠性和易控性,设计了进口为矩形和三角形的超声速喷管;通过设计壁面压力分布,完成了均匀膨胀的轴对称喷管设计,并将设计结果与数值模拟进行对比验证。研究表明:预设三维流线下游未知点的压力值,存在多个流动方向满足该压力条件,即该问题的解不唯一,因此三维超声速流线压力反问题是非适定的。对比验证表明:所设计的流场与CFD计算得到的等值线符合得较好,流场参数的最大误差为1%。因此,所提出的双特征线解法具有一定的可靠性,有望为三维超声速流道设计提供新思路。     

横向喷流对鸭式导弹滚转特性影响研究

通过数值方法求解三维可压缩雷诺平均N-S方程,对导弹横向喷流的干扰流场进行了数值模拟,计算结果与实验数据吻合较好,基本验证了计算方法在横向喷流复杂流场数值模拟方面的有效性。在此基础上,对固定尾翼鸭式布局导弹亚、跨、超声速流场进行了数值模拟,并且计算分析了横向喷流对鸭式布局导弹滚转控制特性的影响,计算结果表明,横向喷流可以有效提高鸭式布局导弹的滚转控制能力。     

稳态周向总压畸变对小型风扇影响的全周数值模拟

采用全周非定常数值模拟的方法，求解两种畸变角时，进口周向总压畸变对某型风扇流场的影响，并详细分析了总压畸变沿风扇的传递情况。计算和分析结果表明，该方法可以清晰描述流场参数沿周向分布情况；总压畸变会引起总温畸变，并且进口畸变强度相同时，畸变角度及动静叶间的相互作用对畸变传递有较大影响。