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

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

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

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

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

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

用交替方向迭代法计算高超音速钝体的粘性激波层

本报告用交替方向Spline隐式迭代求解简化的Navier-Stokes方程组,以给出粘性激波层解。方程建立在参考线坐标系上,引入坐标变换以增加靠近物面处流场的分辨度。在x方向上用隐式迭代,在y方向上用Spline逼近修正激波。主要由物面曲率变化的情况控制激波松弛的次数。计算了球-锥,球-锥-柱等组合体的粘性激波层,并与国内外的一些结果进行了比较。     

高超声速钝体绕流的奇异摄动理论

对于二维和轴对称的物体的无粘高超声速绕流,当比热比Υ→1和来流Mach数M→∞时,激波将贴近物面,Chernyi等人应用薄激波层理论给出了尖头物体绕流问题的渐近展开解,它属于正则摄动理论。但对于钝头体,Chernyi等人的理论遇到了困难,因为该理论中对各物理量的量阶估计在物面附近是不适用的。本文在对物面附近各物理量的量阶重新进行估计的基础上,以Chernyi等人的渐近展开式作为外解,并在物面附近给出了各物理量的新的内解的渐近展开式,从而用奇异摄动理论中匹配的渐近展开法,得到了钝头体绕流问题的解。文中还给出了渐近展开式中首项的解析表达式结果。     

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

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

管内正激波/湍流边界层干扰的数值模拟

本文采用时均Ｎ－Ｓ方程和Ｂａｌｄｗｉｎ／Ｌｏｍａｘ代数湍流模型计算了典型拉伐尔喷管内正激波与湍流边界层干扰流场。计算与实验结果的比较表明，方法可准确地预测激波结构，激波与边界层干扰区流动特征，激波位置、激波前马赫数和壁面压力分布等。     

可压缩球涡和平面激波相互作用的数值研究

行进了可压缩球涡和平面激波相互作用的数值研究。一方面，用非定常Ｅｕｌｅｒ方程数值模拟了一种可压缩球涡的自身演化，同时，利用Ｒａｎｋｉｎｅ－Ｈｕｇｏｎｉｏｔ关系，在流场中嵌入运动激波，求解了激波－球涡相干的流动过程。     

差分计算中激波上、下游解出现波动的探讨

本文探讨了差分计算中激波上、下游解出现波动的原因。研究指出:在NS激波方程中,附加二阶耗散项,可平滑激波;但附加三阶弥散项或四阶耗散项,在一定条件下,可引起解在激波上、下游波动。因此在求解NS激波方程时,若采用二阶或三阶差分格式,由于格式弥散项和耗散项的存在,将引起解的波动。     

高超声速运动中的熵层问题

本文研究了钝锥及钝柱在高超声速飞行中的熵层问题,、分析的基础是把激波层分成熵层和外层两个区域。在熵层和外层内,我们分别建立了气体运动的微分方程组和它们的解。利用两层交界面上的压力相等的条件, 两层内的解被联结起来, 同时,给出了确定流场各物理量的公式。通过分析, 阐明了激波层内压力过度膨胀、二次激波形成及激波出现拐点是熵层对外层流动作用的结果。和电子计算机的数值结果及实验结果比较,理论给出的表面压力、表面M数以及熵层内的密度和熵是符合的,这是以往的理论作不到的。     

运动激波绕射过程中的粘性效应

应用 TVD格式求解二维非定常全 N- S方程 ,以研究运动激波绕射现象。考虑粘性后 ,绕凹角流动时 ,压缩拐角附近表面压力分布与实验结果更为符合 ;绕凸起物时 ,出现反射激波因边界层干扰根部加宽 ,前、后角点附近有局部分离与再附 ,并引起局部压力上升 (与无粘流场模拟相比 )等现象。结果表明 ,全 N- S方程可较好地模拟实际粘性流中复杂激波障碍物干扰流场特性     

压气机叶栅中的激波附面层相互作用

压气机叶栅中激波附面层相互作用是十分复杂的问题,由于相互作用引起分离是决定跨音速压气机性能的重要因素之一。然而迄今为止,尚未深入进行过压气机叶栅激波附面层相互作用的研究,发表的文献极少。作者经过计算和分析,说明压气机叶栅流中主要的相互作用形式是叶栅槽道中激波和湍流附面层的相互作用。作者通过分析指出,研究压气机叶栅激波附面层相互作用,不能直接应用Pearcey分离准则。作者并提出了适用于压气机叶栅的分离准则的函数关系为f(M_1,p_1/p_(L.E.),P_(r.E.)/p_2,Re_0)=O。     

马赫数6柱-裙激波/边界层干扰直接模拟

高超声速激波/边界层干扰比超声速工况下具有更强的可压缩效应,再附之后会形成极高的局部力/热载荷,严重影响飞行器飞行安全。而激波/湍流边界层干扰区附近流动的三维特性使得流动更加复杂而难以预测。采用直接数值模拟对高超声速条件下的柱-裙激波/湍流边界层干扰进行了详细研究,特别是对Görtler涡结构对分离泡、物面压力和热流造成的展向差异开展了定性和定量分析。研究发现,干扰区附近的分离泡结构呈现出明显的三维效应,且Görtler涡展向分离位置截面的分离泡要明显小于再附位置,而这两个截面上分离泡的运动基本同步,没有明显的延迟或超前现象。物面压力和热流在展向出现显著的不均匀性,展向再附位置的平均压力和热流要比展向分离位置分别高13%和16.2%,脉动压力和热流比展向分离位置分别高28%和20%。截面流向速度特征正交分解结果显示两个位置上的能量都集中在剪切层附近,并且展向再附位置上低频模态占有更高的能量。在采用模态重构流场分析分离区面积发现,展向分离位置重构误差更小,而展向再附位置上的重构误差收敛更快。     

三维激波／附面层相互作用的起始分离预测

给出了一种预测三维激盘／附面层相互作用诱导的流动起始分离的方法。研究结果表明：激波与附面层相互作用所诱导的二次流动是影响起始分离的重要因素。当来流相对马赫数大于１．５时，强的激波／附面层相互作用可能导致跨音风扇转子叶尖区域的流动分离。     

一类关联转子分离尺度与损失的工程模型

将分离与损失相关联，提出了利用转子进出口参数来确定跨音风扇转子内部分离流动的方法。应用该方法对几个转子流场计算的结果表明，在转子的叶尖区域，由于强的激波与附面层相互作用，使得流动在激波后不久便发生分离。叶尖漏流涡可能对激波与附面层的相互作用诱导流动分离有较大影响。     

拉瓦尔喷管外发生激波反射工况详细分析

平面拉瓦尔喷管的环境压力与入口总压的比值从设计工况到第三极限工况(管口处为正激波)逐渐增大时,喷管外的激波现象不同.根据已知的激波反射理论,将设计工况到第三极限工况间的工作状况进一步细分为强斜激波工况、马赫反射工况、双解区和正规反射工况、并将由气流偏转角所表示的发生正规反射和马赫反射的条件用环境压力与总压的比值来表示,找出了强、弱斜激波的环境压力与总压比值的临界值,得到这些临界值随着喷管面积比的变化趋势.最后针对面积比为5的平面拉瓦尔喷管,由理论分析得到其管外发生正规反射、马赫反射、强斜激波反射时对应的压比分别为:0.0841,0.0959和0.2005,并用数值模拟来验证理论分析的正确性.      

Shock control on a swept wing

Shock and boundary layer control by contour bumps and local boundary layer suction have been investigated experimentally and numerically on a transonic swept wing. Additional 2-D numerical investigations were performed for the airfoil, corresponding to the wing. The investigations were primarily stimulated by the question concerned with the influence of sweep on the bump effectiveness. This influence has been found to be rather small; the drag reduction by the bump is slightly lower for the swept wing than for the airfoil. A location of the bump in the shock region has shown its effectiveness for reducing shock strength and hence wave drag. A position of the bump downstream of the shock wave has been shown to reduce viscous drag and to postpone buffet-onset to higher lift coefficients. Furthermore, the results indicate that boundary layer suction is a powerful device for drag reduction, but the effectiveness decreases with increasing Reynolds number. Higher effectiveness of suction can be attained, when it is coupled with a contour bump. The parameters height and position (relative to the shock) of the bump, optimized in terms of drag, depend on the shock strength; an influence of the boundary layer thickness upstream of the shock on the optimal bump parameters has not been found. A possibility to control an adaptive bump, mounted on an aircraft wing, is to employ the trailing edge pressure.     

超声速弱欠膨胀冲击射流流场结构

为了对超声速弱欠膨胀冲击射流的流场结构细节进行研究,使用大涡模拟方法对其进行了数值模拟。利用三阶迎风和四阶对称紧致格式对无量纲化轴对称可压缩滤波N-S方程进行空间离散,时间上推进采用的是三阶精度的TVD型Rugge-kutta法。亚格子尺度模型采用的是修正Sm agorinsky涡粘性模型。通过与经典的冲击射流实验比较,证明了程序的可靠性。数值模拟得到了剪切层以及壁面射流中的涡结构和主射流中的激波结构,并且在此基础上对涡合并和板前激波和涡干扰现象进行了深入研究。发现涡合并现象主要出现在流场的上游,越往下游出现的几率越小;涡和板前激波的相互作用会引起激波位置和强度以及冲击平板上冲击区的压强的显著变化,同时也会导致涡的变形。     

考虑附面层影响的二元混压式进气道设计

采用等激波强度的方法,考虑附面层修正,设计了一种飞行马赫数Ma=3.0的二元混压式进气道.通过数值仿真,模拟了激波-边界层的相互影响,研究了附面层抽吸对内流场的影响,获得了进气道内部复杂的流场分布,以及不同背压下进气道的起动特性.计算表明所设计的进气道性能较好,附面层抽吸对稳定正激波有明显的作用,提高了进气道抗反压能力.给出的方法可用于二元混压式进气道的初步设计和验证.      

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

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