全速位积分方程跨音速绕流的数值解

从全速位方程出发,利用Green公式将其化为激波捕获积分方程和激波装配积分方程,然后离散进行数值解。流场出现激波时,对激波捕获积分方程应用上风技术捕捉到激波,然后应用激波装配技术计算,得到了满意的结果。经算例考核,该方法具有计算区域小,收敛快和CPU时间较少等优点。     

积分方程法计算翼型的跨音速绕流

从跨音速小扰动方程出发推导积分方程的过程中,本文用任意形状的封闭曲线C_Q（其极限趋于零）挖去奇点Q,最终得到无奇性（指无穷奇性,不包括Cauchy奇性）的积分方程。 对于跨音速流中的圆头翼型的前缘问题,提出了一种解决办法。 证明了Nixon给出的反演公式对于超临界有激波的小扰动流动也成立。 关于积分方程法中的人工粘性方法,对Sachdev和Lobo提出的方法做了改进。 最后给出了NACA0012翼型在有无升力和有无激波各种状态下的计算结果。比较表明,本方法的计算结果与其它方法的计算结果符合得较好,且计算量很小。     

密切弯曲激波乘波体技术的应用及有效性分析

将密切技术设计乘波体的应用推广至弯曲激波外锥流场中。针对不同激波形状(ICC)约束条件,在凸、凹激波曲外锥流场中,生成了四种构型的密切弯曲激波乘波体,采用数值模拟及理论分析的手段开展了密切弯曲激波乘波体技术应用的可行性验证及有效性分析。研究结果表明:(1)基于密切凸、凹激波外锥流场的乘波体乘波压缩面、出口截面压力分布、激波形状(ICC曲线)均与设计吻合,最大偏差小于5%,说明密切凸、凹激波外锥流场的方法设计乘波体是可适用的;(2)利用宽高比为0.5的超椭圆方程作为ICC控制曲线生成的乘波体,流场压力的理论解与数值解偏差可控制在1.6%以内,而宽高比为2的超椭圆方程作为ICC控制曲线生成的乘波体,流场压力偏差可控制在4%以内。     

计算机翼跨声速绕流的全速位积分方程数值解

本文利用格林定理将全速位方程转换为全速位积分方程，应用高斯定理对沿激波面的显式积分进行简化。将全速位积分方程的解表示成为面元项和场元项之和，然后离散进行数值求解。本文将激波捕捉方法和激波装配技术结合起来，对机翼跨声速绕流进行数值实验，计算结果令人满意。     

跨声速机翼抖振初始迎角N-S方程定常计算分析

 应用代数法和椭圆型方程优化相结合的方法生成翼身组合体块结构网格。采用有限体积空间离散法和五步Runge-Kutta显式时间推进法求解N-S方程。基于雷诺平均N-S（RANS）方程,选用Spalart-Allmaras（SA）一方程模型模拟紊流。根据N-S方程定常计算结果,采用升力曲线、俯仰力矩曲线、后缘压力发散、跨声速激波位置以及机翼表面极限流线等几种判据,对跨声速机翼的抖振初始迎角进行了合理的预测分析。     

粘性高超声速激波—激波碰撞的计算和分析

本文用特片变量型的隐式ＮＮＤ格式计算了五种类型的高超声速激波－激波干扰流场，得到了清晰的流场结构图画。激波和剪切层等分辨得很清楚。计算的物面压力与实验甚为吻合。借助于激波极曲线，对第４类激波干扰结构进行了分析，解释了干扰引起物面压力增高的原因。     

跨声速非匀熵全速位方程的有限元数值解(英文)

本文改进了传统全速位方程的有限元数值解,既保留了全速位方程简单的优点又考虑了出现激波后流场不是匀熵的特点,对密度和后缘库塔条件进行修正,改善了计算结果。文中给出翼型和机翼的算例,其结果与传统全速位方程的结果相比,激波的位置偏往上游,强度较弱,这个趋势使结果较接近于欧拉方程的解。     

后缘连续变弯度对跨声速翼型气动特性的影响

针对后缘连续变弯度对跨声速翼型气动特性的影响进行了研究。首先不考虑翼型后缘连续变弯度,基于搭建的优化设计系统对跨声速翼型进行气动减阻优化设计,通过添加不同的约束优化得到两种跨声速翼型:无激波翼型和超临界翼型。然后在这两种翼型的基础上,以后缘偏转角度为设计变量、以阻力系数最小为目标,针对不同的升力系数分别进行优化设计,并根据优化结果深入分析后缘连续变弯度对这两种翼型极曲线特性的影响机理。优化结果表明:无激波翼型与超临界翼型相比,其设计点处的气动特性较好,但鲁棒性较差;升力系数小于设计升力系数时,应用后缘连续变弯度后,无激波翼型的极曲线特性明显提高,减阻最高达到3.9%,而超临界翼型的极曲线特性提高不明显;升力系数大于设计升力系数时,应用后缘连续变弯度后,无激波翼型和超临界翼型的极曲线特性都明显提高,减阻分别达到2.4%~18.1%和1.7%~13.2%。     

高超声速化学非平衡粘性激波层方程的数值解

本文讨论了高超声速粘性激波层方程数值计算时差分格式引起的物理失真问题。具体分析了全隐格式格式粘性的影响,并作了数值试验。为了验证隐式结果的可靠性。在超声速激波风洞中测量了钝锥的表面压力分布,并与计算结果作了比较,两者基本一致。 本文采用隐式有限差分法数值计算了高超声速化学非平衡粘性激波层绕细长球锥的流动。计算时采用连续方程和法向动量方程耦合求解的方法以解决细长体远后身区计算中的问题。应用网格技术和加强系数矩阵主对角元素优势的方法提高了化学非平衡流计算的雷诺数范围。文中给出了高超声速化学非平衡流的计算结果,并与其它文献的结果作了比较。     

跨声速非等熵势函数方程的数值求解

基于动量方程的特点,本文将只与温度有关的“等熵密度”与熵的变化分离,而由散度型动量方程直接计算出气体经过激波的熵增。与经典势函数方程迭代,可以方便地求解非等熵跨声速流动。从对任意回转面叶栅的跨声速流动计算来看,这种非等熵势函数得到的激波位置可比等熵势函数前移1～2个网格,强度有所减弱,和试验结果较接近。由于各流线的熵不相等,对Kutta条件和出口边界条件作了修正。本文建议的方法也可用于计算流函数的熵增。     

Reflection of a moving shock wave over an oblique shock wave

The reflection of a moving shock wave over a wedge immersed in a still gas and the reflection of a wedge induced steady shock wave over symmetrical and asymmetrical reflecting surfaces have received intensive considerations since more than 70 years ago. Here we consider a different shock reflection problem—reflection of a moving shock wave over an initially steady oblique shock wave induced by a wedge immersed in supersonic flow. For the flow condition we considered, five moving triple points, with each connecting an incident shock wave, a reflected shock wave and a Mach stem, are identified. By using the reference frame co-moving with each triple point, the type of each shock wave of this triple point is clarified. The present study is significant in that it treats a new shock reflection problem leading to a new shock reflection configuration and showing potential applications in supersonic flow with unsteady shock interaction.     

一种形状可控的激波增强管道型线设计新方法

激波管所产生的非定常运动激波,若强度和形状能够按照一定的设计要求进行可控条件下的调节,将可望为燃料点火燃烧试验等提供具有独到优势的研究手段。基于激波动力学理论,针对激波管中所产生的平面运动激波,通过设计特定的上下壁面收缩型线,使初始平面运动激波,经收缩段（包括光滑凹形曲线段、斜直线段和光滑凸形曲线段）的变形和强度增加,再以平面波面形状进入较小截面直管段的连续转变过渡,得到了强度增加的平面激波。进一步对所设计的典型型线分别采用数值计算和试验的方法,考核分析激波运动过程中的形状变化,验证了理论方法的可靠性。在此基础上,分析了型线设计的关键参数对激波增强幅度的影响,结果表明,相对于传统激波管方法,本文中所提出的收缩截面方法能更显著地增加平面激波强度;另外,还考察了初始入射激波马赫数对壁面型线和运动激波波面形状的影响,结果表明,对于较强的初始入射激波来说,壁面型线对入射激波强度依赖较小,也就是说,当实际入射激波马赫数即使稍偏离设计状态时,仍然能得到近乎完美的平面形状增强激波。     

非等强度多道冲击波作用下空泡溃灭机制分析

航空航天、生物医学等领域均存在冲击波与空泡作用问题,冲击波尤其是多道冲击波作用下空泡溃灭包含着复杂的多相瞬变行为与物理现象。基于自主搭建的可压缩多相流并行数值平台,对液体中单一/多道冲击波与空泡的作用过程进行了数值模拟。通过对冲击波作用下空泡内波系性质,以及多道冲击波与空泡作用后的系列反射波系在液体中的相互作用过程进行详细解析,分析了不同冲击波作用下空泡的形变演化过程,探究了空泡的溃灭机制与特性。研究发现,相较于单一冲击波的作用,多道冲击波作用下空泡内部及周围液体流场中的波系结构更为复杂,然而无论是单一冲击波或是多道冲击波工况,其最终诱发的空泡溃灭形态均十分相似。而且,当空泡发生溃灭前所受到多道冲击波的总强度与单一冲击波强度相等时,空泡的溃灭波强度也十分相近,为多道冲击波与空泡作用效果的定量化评估提供了理论基础和依据。     

高马赫数来流条件下斜激波与流向涡对相互作用

针对斜激波入射流向涡对形成的复杂流动干扰现象,在来流马赫数6条件下,利用后掠斜坡式涡流发生器生成流向涡对,采用数值模拟结合激波风洞实验,研究了斜激波与流向涡对相互作用引起的激波变形和旋涡演变。结果表明,斜激波与流向涡对相互作用的剧烈程度大致分为强、中、弱三种。在强相互作用下,激波面向上游凸起的程度最剧烈,呈现类“T”字型特征,其对称面处接近于一道正激波,波后出现回流区和局部流动滞止,使得流向涡对发生破碎;在中等相互作用下,对称面处存在局部正激波,但没有出现回流区;在弱相互作用下,不再出现局部正激波。进一步地分析发现,斜激波与流向涡对相互作用产生的流向涡对破碎现象,不符合经典的激波作用下单个流向涡破碎理论的预测。采用流向涡对对称面处的气流参数,对经典的流向涡破碎理论进行修正后,能够预测本文流向涡对在斜激波作用下是否发生破碎。本文为掌握激波作用下流向涡对破碎现象的规律,提供了重要参考。     

冷气喷射法控制激波强度的数值研究

对冷气喷射时激波受到的影响进行了多方案数值研究.在高压级静叶吸力面反射激波生成点前后5个不同位置上采用相同总压及相同温度的冷气喷射,分析了不同位置的冷气喷射对激波强度和方向影响.结果表明:在吸力面激波折射点附近喷射冷气是一种行之有效的控制激波强度的方法.冷气喷射位置位于激波折射点附近时能够对激波的强度产生影响,在接近激波折射点前部位置注入冷气对流动有积极作用,减弱了激波.      

基于LES方法的平板非定常激波/湍流边界层干扰研究

以高超声速发动机进气道湍流分离控制为应用背景,采用大涡模拟(LES)方法进行马赫数为3.0(唇口附近马赫数约为3.0)的激波/湍流边界层干扰(SWTBLI)流场机理研究.利用扰动循环引入的方法,先得到充分发展湍流场,然后根据斜激波关系式引入激波的方法进行激波/湍流干扰模拟.研究结果显示:充分发展湍流场在激波作用下产生逆...     

波转子非定常泄漏流动机理

针对制约波转子性能的泄漏问题,通过提取波转子中与泄漏相关的主要流动现象并建立波转子单通道泄漏模型,对非定常泄漏流动机制进行了详细的数值研究,对本文给出的非定常泄漏损失预测模型进行了数值验证。结果表明:转子通道中存在不同程度的压力波动,波动幅值与间隙宽度有关;连续反射膨胀波、周期性出现的弓形激波及其反射激波是压力波动的根本原因;间隙内部泄漏过程存在3个主要的流动阶段;泄漏过程中通道激波传播速度不变、波后时均压力不变;在一定间隙宽度范围内,激波马赫数、激波静增压比与无量纲间隙宽度均呈线性关系,当间隙宽度从0增大到0.08时,激波马赫数衰减7.3%,激波静增压比衰减10.1%;泄漏流动通过泄漏产生的主膨胀波对激波传播过程施加影响,通道激波衰减本质上是理想激波与主膨胀波叠加效应的结果;泄漏损失预测模型与数值结果吻合良好。     

爆轰波-激波迎面作用后的稳定波系

本文为对当量比乙炔氧气混合气体中爆轰波与激波正面对撞产生稳定波系的实验和理论研究.实验主要以高速扫描摄影获取两波对撞的x-t纹影图,一维理论分析则基于三种热完全的组分求解两波对撞的稳定解并探寻它们的规律.实验发现透射波系包括一道激波和爆轰波,以及紧随爆轰波后的稀疏波区,这种波系情况与一维理论分析中CJ解一致.透射CJ爆轰与人射爆轰相比马赫数降低,而相对波前来流的传播速度有轻微提高,但在实验室坐标下其速度显著降低.透射波系受初始压强影响不大;初始温度提高使得爆轰波速度降低,而透射激波速度增加;对波系起实质影响作用的是入射激波强度,激波越强,则整个透射流场呈现偏向激波的趋势;理论分析还指出,稀疏波区的出现不可避免,当激波强度趋于声波时稀疏波区趋于消失,激波越强则疏波区趋于扩大.     

掠型设计调控跨声速离心压气机内部流动机理

跨声速离心压气机面临进口激波控制的难题,探索前缘掠型对跨声离心压气机影响机理对于进一步提升离心压气机性能具有重要意义。以两台跨声速离心压气机为研究对象,利用经过校核的数值模拟手段探究了不同前缘掠型下压气机的性能及流场变化。结果表明:设计工况下叶轮入口的激波位置直接决定了前缘掠型设计的作用效果。掠型设计通过改变导风轮叶顶载荷分布,进而影响其通流能力,同时改变了其入口激波结构和流动损失。当叶轮入口叶顶激波呈吞入状态时,采用前缘后掠会减弱槽道激波强度,从而减小导风轮叶顶附近损失;当叶轮入口叶顶激波呈脱体状态时,采用前缘前掠会减弱前缘激波强度,使得激波更加附体,从而减小导风轮叶顶附近损失。      

燃气射流起始冲击波形成机理的实验研究

利用实验研究的方法对燃气射流起始冲击波的形成机理进行研究。目的在于研究其发生、发展的规律,为寻求降低其危害性的研究提供参考。实验中使用了激光Ｍｏｉｒｅ偏拍仪对起始冲击波场的形成初期进行了光学显示,同时还使用了压阻式压力传感器对起始冲击波超压进行了测量。通过对测量结果的分析,认识了燃气射流起始冲击波的物理本质且总结出其形成和发展的一般规律,它在近场的高峰值超压对火箭武器的承载设备具有破坏作用。