轴对称内锥流动中马赫盘的形成与演化

针对轴对称内收缩直锥流场中的马赫反射问题，采用无黏数值模拟结合理论分析，在来流马赫数6及直锥入口半径R恒定的条件下，以模型壁面前缘角度θ_w和壁面长度w/R为几何参数，研究了内锥形激波汇聚增强过程以及马赫盘下游流动对马赫盘位置的影响。结果表明，壁面前缘角度θw通过改变轴对称内锥形入射激波汇聚增强过程中达到von Neumann强度的位置，限定了马赫盘可能存在的位置范围；马赫盘位置还受其下游流动的影响，在不同的θ_w和w/R条件下，马赫盘下游流管内声速喉道的形成机制可以划分为依赖或不依赖壁面尾缘膨胀波两种类型；对于不依赖于壁面尾缘膨胀波的情况，声速喉道形成之前流管内的压力仅与反射激波的波后非均匀压力相匹配，马赫盘的位置不随w/R而变化；对于依赖于壁面尾缘膨胀波的情况，声速喉道形成之前流管还需要匹配壁面尾缘膨胀波引起的压力变化；根据反射激波下游与马赫盘下游流管的压力匹配关系，可以求解马赫盘位置。     

二维激波管有障碍物的激波反射流场计算

本文应用最小耗散的混合格式及MacCormack二阶格式加θq╱2耗散修正格式计算了一维激波管Riemann问题、二维激波管中激波与楔的反射流场。这两个格式的特点是对任何初值其耗散项永远不等于零,因而,计算得到的解为物理解,而且除间断外,均具有一致二阶精度,计算出的单波很精确,并且激波分辨度高。本文计算了楔角θ~*=30°、45°,激波马赫数M_s=5.29、10的情况,给出了流场的非定常过程及流场的自模解,一维结果与Riemann问题理论解、二维结果与实验结果的比较,说明计算是成功的。     

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

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

DMD方法对流场特征捕捉能力的研究

为论证动力学模态分解(DMD)方法对流场主要特征的快速识别能力,分别计算分析了有、无激波产生的2组翼型中相同工况下的翼型振荡运动。在S-A方程湍流模型下,得出了翼型表面的流场压力数据,并完成DMD分解。在DMD模态云图及特征频率中,可以识别出翼型周围流场的流动结构与主要的动力学信息。结果表明,有、无激波算例中分别捕获了翼型的非定常流动特征,前4阶DMD模态反映了一系列的流场特征信息:一阶模态可以基本匹配均匀流场的形态;二阶DMD模态的频率则与翼型预设的振荡频率恰好匹配,剩余的三、四阶模态则正好反映了这一振荡频率的倍频。在有激波流场算例中,后3阶模态均准确地捕捉到了激波间断处;无激波的流场算例中,后3阶模态则可以有效地捕捉到大压差区域。     

应用γ-Reθ湍流模型和PSE方法研究高超声速流动

采用考虑转捩的γ-Reθ湍流模型研究了高超声速复杂流动,结合PSE(parabolic stability equa-tion)稳定性分析方法和e-N方法,通过给出边界层最不稳定点的参数,用以控制γ-Reθ的转捩经验关联式,减少对经验参数的依赖.计算了锥体高超声速流场,得到的壁面换热量和边界层转捩位置与实验吻合较好.计算了平板激波/边界层干扰和高超声速拐角压缩流场,得到了准确的波系结构、压力分布和分离区的大小.结果表明,该方法与传统湍流模型相比具有较大的优越性.     

来流参数对混合管道流场结构的影响

为了探究来流参数对混合管道流场结构的影响,通过试验和数值计算方法研究了在不同马赫数、静压来流参数条件下的背景流场和激波串流场结构。试验采用高速纹影和压力测量设备分别进行流场的观察和测量,采用喷管支板结构将来流分为上侧二次流、一次流和下侧二次流。试验与数值结果均表明,背景流场中的背景波系、混合层均由支板后缘处产生,背景波系在壁面间反射并与混合层发生干扰作用,背景流场结构较为复杂。一、二次流的静压关系影响混合层的初始发展方向,导致θ从-9°～7°变化,而马赫数决定起始激波与混合层的夹角。依据前缘激波的反射类型、亚声速区面积等特征,混合管道的激波串流场可以分为规则反射型、临界反射型、正向马赫反射型和反向马赫反射型四种典型结构。临界反射型的波后流场在垂直流动方向上具有"亚-超-亚-超"间隔流动的特点。     

高分辨率激波/边界层干扰时间演化过程分析

采用基于纳米示踪的平面激光散射（NPLS）技术,探索了流场超高帧频成像测试研究的试验系统,主要解决了多个单脉冲激光器并联后的稳定性和合束、阵列CCD相机的整体设计和布局以及测试系统的同步精确控制等问题,通过对系统的时序和分系统调试,实现了测试系统的精确控制。基于此系统,在单位雷诺数为6.30×10⁶/m的条件下,在马赫数为3.4的超声速低噪声风洞中开展了θ=20°激波发生器入射激波与来流壁面湍流边界层干扰相关试验研究。在试验条件下获得了序列连续时间相关的激波与湍流边界层干扰的瞬态流场精细结构图像,并分析了其流场结构的时空演化特性。     

动态间断装配法模拟斜激波壁面反射

基于非结构动网格技术和边界装配思想提出了动态间断装配法,该方法能够应用于求解含有间断的流动问题。无论入射激波还是反射激波都是作为边界进行处理,激波运动速度由兰金-许贡纽(Rankine-Hugoniot)关系确定。激波作为动网格的一部分,其运动由动网格技术实现。采用该方法模拟了超声速二维流场中激波与壁面相交问题,并且与捕捉法进行比较,二者的流场结构符合良好,但是在细节上还是存在明显差异。通过对流动结构的分析,得出采用装配方法得到的流场要优于捕捉方法的结论。激波壁面反射的问题模拟,也说明了边界激波装配方法对于复杂的激波相交问题是具有处理能力的。     

超燃进气道激波/湍流边界层干扰

针对超燃进气道湍流边界层/激波干扰引起的分离问题,采用基于5阶WENO数值格式的大涡模拟(LES)方法开展流场湍流非定常预测,旨在分析进气道湍流化技术实现进气道起动的可行性。研究表明,平板激波/湍流边界层干扰(STBLI)问题,LES方法能够清晰、可靠预测反射、分离激波形成过程及激波与充分发展湍流边界层的相互干扰,定量结果与试验一致;进气道研究方面,层流状态下,激波干扰产生强分离,导致进气道堵塞,而采用湍流化控制后试验和计算均表明流场分离明显减小,流场稳定且无明显堵塞现象,进气道可以起动,总压恢复系数达到要求,该结果表明,利用强湍流化减弱分离,实现进气道起动思想是可行的。      

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

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

电传飞机系统动态特性识别

提出了一套识别飞机系统动态特性的方法。对于短时激励信号,可同时采取时域法和频域法;对长时激励信号仅采取频域法。方法验证采用了变稳飞机电传系统试飞数据和专门的仿真试验数据。结果表明:该方法可解决电传飞机系统识别和评定中的许多问题,并为更进一步的应用研究打下基础。     

P-S_a-S_m-N曲面的求法

 本文论述了在疲劳可靠性分析和设计中求P-S_a-S_m-N曲面的必要性,探讨了P-S_a-S_m-N曲面一般形式的选择原则,提出了矩法和最大似然法两种疲劳试验方案,分别给出了根据试验所得的疲劳性能数据求P-S_a-S_m-N曲面的方法,并用实例加以说明。     

论有限元素理论中SUPG方法与PG方法在解高速粘流时的等效性

对流场中含激波问题,若用有限元法来数值模拟,目前较有效的方法是SUPG方法和PG方法。本文讨论了这两种方法在数值求解N-S方程组时具有等效性,并对它们的应用前景作了简单的评估。     

预浸料制备方法影响复合材料湿热稳定性的原因分析

分别用热熔法和溶液法制备 T3 0 0 /540 5预浸料并在相同的工艺条件下制成相应的复合材料层压板。比较 2种复合材料的树脂含量、内部缺陷 (用 C-扫描技术 )、空隙率 (用光学显微镜法 )和纤维 /基体界面粘结(以 [0 2 /90 1 2 /0 2 ]层压板的横向裂纹密度和断口形貌评价 ),发现热熔法复合材料的孔隙率较低和纤维 /基体界面粘结较好是其在水煮条件下稳定性优于溶液法复合材料的主要原因     

能量法和加权残值法的联合应用——构造有限元的新途径

 本文将能量法和加权残值法结合起来(记为EWR方法),构造非协调位移元,特点是(1)做法简便——用加权残值法来放松单元间位移连续条件。可使能量泛函中不再包含不协调位移引起的能量贡献项,于是单元只有位移一类变量。(2)使用可靠——由于将加权残值法与能量泛函统一考虑.权残方程可赋予明确的物理意义。一方面。适当选择权函数,权残方程的满足就意味着非协调元收敛准则的满足。另一方面,变换权函数,可将现有各种位移型单元纳入EWR方法的统一形式中。     

基于随机场局部平均的随机有限元法在实特征值问题中的应用

 本文将基于随机场局部平均的随机有限元法应用于结构分析中的实特征值问题,给出了2次摄动求随机实特征值和特征向量的1、2阶矩的有限元公式。算例表明,本文方法随单元数的收敛速度比现有的方法快得多,在较少单元数下即可获得满意的结果。     

粘弹性结构振动特性分析的摄动求解技术

 本文根据粘弹性结构的特点,对其振动特性进行简化分析,构造了这类结构振动特性分析的摄动求解方法,给出了低阶摄动求解公式。还讨论了零阶摄动有重特征值时的求解。     

应用能量－Casimir方法研究充液卫星系统的运动稳定性

利用能量－Ｃａｓｉｍｉｒ方法研究了带４个挠性粱的充液卫星系统在无外力矩状态下,关于绕铅垂轴稳态转动的非线性稳定性条件,该条件考虑了液体的涡旋、弹性梁的振动、卫星主刚体的旋转以及流－弹－刚之间的耦合,此外还考虑了离心力与Ｃｏｒｉｏｌｉｓ惯性力的影响。推广了Ｒｕｍｊａｎｔｓｅｖ和Ｍａｒｓｄｅｎ的部分结果,为带挠性粱的充液航天器的运动稳定性总体设计提供了可靠的理论依据。     

克服欧拉方程奇异性的双欧法

论述一种解算飞机姿态角的新方法，它利用正、反欧拉方和是精华区倒挂关系进行分区接替运算，充分发挥二者的长处，把精华区扩展到全域，消除了奇异区，比其他方法具有解算准确，使用简便的显著优点，飞行仿真实践表明它根除了奇异性，是的全姿态方程，其原理对其他工程和物理领域都有广泛推广应用价值。     

极大似然法对比试验研究及其试验数据处理

 基于中值Ｓ－Ｎ曲线的三参数幂函数式，根据极大似然原理，分别给出了测定Ｐ－Ｓ－Ｎ曲线的单点试验、成组试验和单点－成组试验结果的Ｐ－Ｓ－Ｎ曲线公式。并进行了对比试验，使用文中方法统计分析试验数据，发现单点法、单点－成组法及成组法测定出的Ｐ－Ｓ－Ｎ曲线较接近。由于单点法能够大量节省试样，且能保证一定的精度；因此它适合构件和模拟件Ｐ－Ｓ－Ｎ曲线的测定