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

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

斜激波与内凹半圆柱面无粘相互作用

针对三维内转式进气道流动中激波入射内凹壁面主导的复杂波系干扰问题，采用斜激波入射内凹半圆柱面的简化构型，在来流马赫数为6的条件下，通过无粘数值模拟结合理论分析，研究了激波角βi=14°～29°的斜激波在内凹半圆柱面反射形成的三维流场。结果表明，流场对称面均会出现显著高于二维情况的逆压梯度。当βi≤25°时，从侧壁到对称面，斜激波经历了从马赫反射（MR）到规则反射（RR）的转变，形成了MR-RR型流场，转变点处产生的桥激波向对称面延伸，桥激波在对称面反射后产生的压力峰值高于二维斜激波反射；当βi≥25°时，斜激波在侧壁和对称面均发生马赫反射，形成了MR-MR型流场，两种马赫反射分界点处产生的桥激波向侧壁发展，侧壁气流在对称面相撞后产生的压力峰值高于正激波后的压力；当βi=25°时，流场存在MR-RR型和MR-MR型双解现象。通过降维分析理论，揭示了两类流场中转变点和分界点的形成机制，并厘清了桥激波的产生原因和初期演化特征。当βi≥18°时，无粘激波干扰所主导的侧壁气流加剧向对称面汇聚，并在对称面附近产生流向涡对。无粘分析获得的认识，有助于揭示内转式进气道中流动汇聚和流向涡对等现象的形成机...     

斜激波与弓形激波相互作用的彩虹纹影实验研究

斜激波与弓形激波相互作用形成的波系结构可以分为六类,其中IV型激波干扰产生了超声速"喷流",产生了极高的热载和压载。基于彩虹纹影显示技术,实验获得了六类激波干扰彩虹纹影照片,分析了各类激波干扰的波系结构。通过数值计算,获得了有无斜激波干扰的条件下,壁面附近温度与压强的分布,讨论了IV型激波干扰产生极高热载和压载的原因。研究表明,激波干扰发生以后,钝头体表明压载与热载急剧增加,尤其是IV型激波干扰,形成的超声速"喷流"射向钝头体壁面的超声速射流,会导致峰值压力剧烈升高,其数量相当于来流静压的近40倍,驻点附近的热流升高将近8倍。     

压力恢复系统扩压器激波串现象的数值模拟

通过求解由S-A一方程湍流模式封闭的三维雷诺平均Navier-Stokes方程,采用有限体积方法对激光器压力恢复系统变截面扩压器常温状态下的三维流场进行了数值模拟.数值计算和实验符合良好,较好模拟了扩压器中由激波/边界层干扰诱导的复杂流场的流场特性并再现了流场中的"激波串"和"伪激波"现象.本文对激波串的形成机理进行了分析,发现了几种变截面扩压器特有的激波串结构,并研究了扩压器扩张角对扩压效率的影响.     

马赫数对后掠激波和湍流边界层干扰特性的影响

本文介绍了尖前缘翼诱导激波和湍流边界层干扰流场壁面特性，着重强调马赫数影响。给出２．０≤Ｍ≤８．２、ａ≤３５°分离流场中，锥型干扰区内主分离线和再附线位置与无粘激波角β０和迎角ａ的相关式，证实无粘条件是控制锥型区尺度的主要因素，面高超声速与超声速干扰流中二次分离随激波强度的不同发展，表明干扰流场的细致结构与可压缩性有关。     

考虑粘性时激波反射及壁面热流计算

本文采用显式ＴＶＤ格式，求解薄层Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ方程，计算了激波与楔面相互干扰的流场情况，得出了马赫反射的各种等强度级及固壁上详细的物理分布，根据计算得出的流场，计算了壁面摩擦阻力及热流分布。     

吸气式飞行器进气道唇口三维激波/激波干扰

对吸气式飞行器进气道唇口处三维曲面激波/弓形激波干扰流场进行数值模拟,利用典型三维气动干扰试验对采用的数值计算方法进行验证。利用拼接网格技术及逆距离加权插值方法获得入口处流场的非守恒变量,作为激波干扰研究的入口边界条件。数值模拟表明,唇口处激波干扰流动的三维效应十分显著,曲面激波与弓形激波产生斜交,尽管唇口前缘半径很小,但Edney提出的6类激波干扰类型可能沿唇口展向方向同时存在；第Ⅲ和Ⅳ类激波/激波干扰的诱导使得唇口热流分布异常严酷；激波相交处形成斜向“伤疤”状局部高热流条带,峰值热流达到参考热流的4~6倍,可能引起唇口结构的局部烧蚀或破坏。      

高负荷跨声速涡轮激波损失机理及控制技术研究

为了指导高性能、高负荷跨声速涡轮的设计,对其叶栅内的流场结构、尾缘波系结构、减小激波损失的机理及其控制技术进行了分析研究。结果表明:跨声速涡轮尾缘流场结构复杂,存在分离膨胀波、分离激波、基底区、再附激波、尾迹、吸力面反射波甚至激波边界层相互干扰等流动现象。通过采用收缩-扩张通道,喉道后采用直线型吸力面,减小吸力面尾缘弯折角、尾缘厚度和尾缘附近局部修型等措施,从而减弱激波强度,减小激波损失。     

激波绕流与激波碰撞的数值模拟

本文应用基于Poe平均构造的最小耗散的扩散与反扩散守恒差分格式，数值模拟单激波绕流三我形障碍物，及等强度激波碰撞和不等强度激波碰撞的二维非定常复杂流场。计算结果表明，该格式具有很高的单波精度及较好的激波与接触间断分辨度，可直接用于复杂流场的数值模拟。     

高马赫数下激波湍流边界层干扰数值研究

应用GAO—YONG可压缩湍流方程组数值模拟了入射斜激波／平板湍流边界层相互干扰现象,计算了来流马赫数为5．0,激波入射角度分别为15．876°、23．287°两种不同激波干扰强度下的流场。计算程序中的对流项、扩散项分别采用二阶ROE格式和二阶中心差分格式离散,并用多步Runge—Kutta显式时间推进法求解空间离散后的控制方程。计算较好地模拟了高马赫数下的激波／湍流边界层干扰的流场结构,位移边界层厚度,动量损失厚度等,也比较准确地预测了平板壁面压力、摩阻系数等气动力参数的分布。     

燃烧固体药柱内腔表面的典型裂纹及其力学行为

根据燃烧环境中高装填固体推进剂药柱内应力的分布特征,提出了药柱内腔表面的两种典型裂纹, 即横向表面裂纹和纵向表面裂纹, 并分别采用方坯药柱和圆盘药柱表面预制的裂纹进行了模拟。比较燃气作用下两种药柱的有限元分析结果, 明确了挤裂模式对于燃烧环境中固体药柱表面裂纹问题的适用性。这些结论与实验观测现象非常吻合, 对实际发动机工作环境中固体药柱表面裂纹的力学行为预示和故障分析有着重要的指导意义。     

顶部间隙对超声速膨胀器流动特性的影响

隔板与机匣之间留有间隙,间隙的存在势必会对超声速膨胀器的内部流场和总体性能产生影响,为了获得超声速膨胀器内部间隙流动的流动细节,采用三维雷诺平均Navier-Stokes方程和标准k-ε湍流模型,就顶部间隙对超声速膨胀器流动特性的影响进行了数值研究。结果表明:膨胀流道出口斜激波导致吸力面压力高于压力面,隔板尾缘附近部分泄漏流体经间隙流回压力面侧;间隙的存在导致吸力面进口及中、后部近下端壁压力上升,而压力面前缘附近压力下降,对比同一隔板位置,间隙高度每增加1%喉部高度,超声速膨胀器隔板载荷系数最高下降2.6%;端壁损失和斜激波损失降低,但产生了泄漏损失,三维流道内总的流动损失增加,膨胀器效率降低,本文研究范围内效率最多下降8.8%;马蹄涡、泄漏涡及二者之间的相互作用是顶部区域的主要涡系结构;前缘附近气流经间隙流到吸力面侧和尾缘附近泄漏流体越过间隙重新流回压力面侧是间隙内气流的主要运动形式。     

AP/HTPB复合推进剂低压燃烧特性数值研究

为研究AP/HTPB复合推进剂低压燃烧特性,使用二维三明治模型进行了数值仿真研究。该模型使用简化的四步总包反应机理,考虑燃烧过程中的燃面退移,计算时大粒径AP取平均粒径,小粒径AP采用“均质黏合剂假设”,认为其与HTPB组成“均质黏合剂”。在0.02MPa-0.08MPa的压强下进行了仿真,并与高压（5MPa）下的燃烧仿真进行对比。结果表明,低压下固相AP对黏合剂的传热导致了交界面处燃面的下陷,下陷程度随压强减小而增大;随着燃烧的进行燃面各点燃速趋于一致,燃面形状不再改变。低压下气相反应表现为预混燃烧,初焰为主导反应。低压下燃速仿真与实验值基本一致,压强指数小于高压下的压强指数。     

叶片全三维反问题优化设计方法

基于数值求解Navier-Stokes方程,对叶轮机械叶片全三维反问题优化设计方法进行研究.反问题给定叶片表面的静压分布反求叶型,假设叶表网格点有虚拟移动速度,迭代过程中根据叶表的实际静压和目标静压之差来修改叶型.对单排叶栅进行了反问题设计,并给出了初步的优化算例,结果能很好地满足目标静压分布,验证了该方法的有效性,计算的时间和精度也适合进行工程应用.      

转静干涉对叶片非定常表面压力的影响

为了研究转静干涉对叶片表面压力的影响,通过数值模拟的方法,计算了一级低速轴流压气机流场,得到了充分收敛的周期性非定常流场解。在此基础上分析了叶片表面压力及气动力的变化规律,并通过改变背压,研究了流量对叶片表面压力和气动力的影响。最后在低速轴流压气机实验器上进行了实验验证,结果表明计算和实验的结果吻合得较好。     

基于0.6 m量级三声速风洞的压敏漆试验技术

基于0.6 m暂冲式三声速风洞，建立了压敏漆测压系统，解决了各分系统的同步控制问题。研究了涂料喷涂影响、图像滤波和系统测量稳定性及精准度等技术细节，并将该系统首次应用于未来大型客机减阻与激波控制的机翼表面压力测量中，获得了基本外形和鼓包外形机翼表面的压力分布、激波位置及形态。检验了设计鼓包在设计状态和稍偏离设计状态下的激波控制效果及其对上翼面压力分布和升力特性的影响。研究结果表明：涂料喷涂质量不佳造成的表面粗糙度和厚度变化会显著影响压力分布，喷涂质量需严格控制。窗口直径8像素迭代3次的高斯滤波对压力波动的平滑效果较好且不会失真。建立的压敏漆系统与压力传感器的压力系数测量均方根偏差在0.022以内，压力均方根偏差小于620 Pa，测量精准度较高。设计鼓包在设计状态及稍偏离设计状态下，均能够有效减弱激波强度，保证机翼升力变化很小，从而提高机翼的升阻比。     

超声速有翼飞行器的无粘分区数值模拟

本文采用ＭａｃＣｏｒｍａｃｋ两步显式格式，用空间推进法求解Ｅｕｌｅｒ方程来对有翼飞行器的超声速绕流流场进行数值模拟，应用特征边条件处理场面边界及装配激波的计算，在翼缘及表面斜率不连续处，采用局部解析处理，提高了算法的可靠性，复杂外形的计算采用分区处理技术，网格由代数形式生成，在大攻角的计算中，通过对横流速度型的修正改善了背风面压力分布特性，解决了推进中断问题，同时模拟出背风面涡。     

叶根开槽叶栅对角区分离的控制

吸力面/端壁角区的分离是轴流压气机流场中固有的现象。本文采取在叶根处开槽的方法,利用压力梯度从叶片压力面向叶片吸力面引入一股射流增加分离区的能量,从而减缓分离。通过数值模拟的方法分析了在不同攻角下槽的位置、大小和形状对扩压叶栅性能的影响,计算结果初步表明,在非设计工况下适当位置、大小和形状的槽可以有效地减小角区分离。     

后置静叶周向弯曲对表面非定常压力脉动的影响

在目前的叶轮机械中,普遍存在着旋转动叶与后置静叶干涉产生的非定常流动现象。这种流动现象会导致静叶表面产生非定常性压力脉动。这种脉动可分为旋转动叶的速度亏损尾迹引起的叶片表面的周期性脉动和气流本身的湍流脉动及其他涡团导致的随机脉动两个方面。采用一工业用单级通风机,对其动静叶的干涉现象进行研究。通过在后置静叶表面沿弦长及叶高方向布置动态压力传感器的方法,研究静叶的周向弯曲（倾斜）角度对静叶表面的非定常力的影响。研究结果表明：采用周向弯曲静叶会改变动叶尾迹所引起的后置静叶周期性脉动的幅值在静叶表面的分布特点,静叶顶部压力脉动的幅值与原始径向静叶相比能够有效地降低,同时与中部的压力幅值相比有一定提高,但整体呈现下降趋势;另一方面,静叶表面随机性的压力脉动强度受周向弯曲的影响较小。     

NASTRAN在压力系数分布形态分析中的应用

针对机翼压力系数分布形态分析工作存在辨识困难的问题，通过对压力系数数据格式的变换处理，利用MSC．NASTRAN for Windows2．0的图形显示功能，将压力系数分布的分析辨识工作由分析数据文件变为查看在电脑屏幕上可视化的分布图形，极大地提高了压力系数分布分析工作的效率。