前飞速度和升力偏置量对共轴刚性旋翼气动特性影响分析

共轴刚性旋翼前飞状态的气动特性主要由工况环境中的来流速度、密度和桨叶的翼型配置、弦长分布和扭转分布等气动布局参数决定。气动布局参数的综合影响决定了共轴刚性旋翼的的升力偏置量。了解前飞速度和升力偏置量对前飞性能的影响规律有利于设计更适合于高速飞行的共轴刚性旋翼。因此,本文通过求解可压雷诺平均N-S(Reynolds-averaged Navier-Stokes,RANS)方程对4 m直径的由两副2片矩形桨叶旋翼构成的共轴刚性旋翼模型的前飞流场进行了数值模拟,获得了不同前进比下的气动力并对不同升力偏置量下的旋翼性能进行了对比。数值模拟结果表明,随前进比增大,桨叶展向拉力分布更加趋于合理,拉力中心向桨叶中段移动,可以充分给桨尖卸载;旋翼升力主要由前行侧桨叶提供,升力偏置量过大容易产生激波诱导失速,不利于高速前飞。     

计算气动光学研究进展

计算气动光学是采用数值计算的方法研究高速空气动力流场对光波传输和光学成像影响的一门交叉学科,在红外成像制导设计等工程应用中发挥着重要的作用。总结了计算气动光学研究的三种基本计算方法,包括基于CFD简化方法的光学统计估算、基于RANS流场的计算和基于LES/DNS瞬态流场的计算,对这三种计算方法的基本理论、计算过程进行了介绍,综述了新近的一些计算气动光学研究进展,最后指出了未来的研究重点和方向。     

一种基于求解椭圆型方程的结构动网格生成方法

基于椭圆型网格生成法,实现了一种简单高效的贴体结构动网格生成方法,可用于具有移动边界问题的非定常流动数值模拟。该方法提出,在网格变形过程中,Poisson方程需要的控制网格间距和正交性的源项可以通过提取已知的静态网格源项直接得到,并在整个动网格生成过程中保持不变。因此,在椭圆型网格生成中需要通过外迭代确定源项的过程可以得到省略,而且该方法不需要人工指定参数。这使得方法具有高效和易于嵌入到已有程序中的特点。数值模拟结果证明,采用这种方法获得的网格能够较好地保持静态网格原有的正交性和光滑性,在相同迭代步数约束下,网格求解效率低于传统弹簧模拟法,但鲁棒性优于弹簧模拟法。     

反推状态下大涵道比涡扇发动机气动稳定性预测与评估

为了预测与评估反推状态下,反推气流再吸入对大涵道比涡扇发动机气动稳定性的影响,采用反推气流扰流流场三维CFD数值模拟、发动机整机稳定性计算分析以及反推状态下发动机进气畸变台架试验相结合的方法,开展了反推气流对大涵道比涡扇发动机气动稳定性影响的研究。通过三维CFD数值模拟手段,捕获了反推状态下发动机进口流场的畸变程度。在此基础上,通过采用稳定性计算程序预测了发动机的气动稳定性,并进一步通过发动机台架试验,验证了预测结果。CFD计算结果表明,随着相对来流马赫数的减小,反推气流被发动机重新吸入的可能性不断增大,当相对来流马赫数减小到0.05时,外侧发动机进口的流场畸变情况变得最为严重。进气畸变情况下的整机稳定性计算分析以及发动机台架试验结果表明,在所考核的目标状态,若只存在因反推气流再吸入引起的进口流场畸变,是不会导致发动机失稳的。     

径向入射激波聚焦实验和数值模拟研究

为得到不同入射激波马赫数对凹腔内径向入射激波聚焦过程的影响,对马赫数为1.25和2.18的径向入射激波在实验段内反射聚焦的过程进行了实验和数值模拟研究。用高速CCD拍摄了凹腔中气流流场的纹影图片,并采用保持强稳定性的龙格—库塔格式、加权基本无振荡算法和块结构动态自适应网格加密对激波反射聚焦的过程进行了数值模拟,实验和数值模拟结果吻合较好。复杂的激波形状及其之间的相互作用(如激波反射,激波聚焦和不同激波的相互作用)在实验和数值模拟中都能很好地观测到。通过比较不同马赫数入射激波的反射聚焦过程,发现较高的马赫数下聚焦后能产生射流,而较低的马赫数下却没有。同时发现,入射激波马赫数能影响二次激波出现的位置和强度,进而影响聚焦过程流场的紊乱程度。     

氢燃料燃烧加热风洞中水蒸气相变效应的数值研究

为明晰氢燃料燃烧加热风洞中的水蒸气相变效应,采用同时考虑凝结和蒸发相变模型的多组分输运可压缩流动数值模拟方法,研究了超声速喷管内的水蒸气凝结过程以及凝结气流再蒸发对典型压缩构型(二维斜劈、三维斜劈和圆锥)流场参数及气动力的影响。结果表明:氢燃料燃烧加热风洞喷管气流中,水蒸气在较低总温条件下确实存在较大程度的凝结相变,但同时在斜劈激波压缩后也存在较强的蒸发相变;凝结会使喷管流场参数发生较大偏离,但再蒸发效应又会使气流经过斜劈压缩后趋近理想状态,趋近的程度受到压缩程度的影响。本文计算条件下,喷管内水蒸气相变导致出口参数偏差大于10%,但经过12°斜劈压缩后,流场参数以及表面气动力可恢复到理想状态的95%以上,此时相变的影响基本可以忽略。     

发散磁场中等离子体加速和推进性能数值研究

稳态螺旋波等离子体推力器中,源室放电获得的一定能量的等离子体经过磁喷管加速产生预定的推力和比冲。为了分析在发散磁场约束下,等离子体的运动受约束磁场和内能变化的影响规律及其推进性能,引入了考虑电子和离子不同响应的二维轴对称数值模型。计算了入口中心磁感应强度[B0]为100～500G、电子温度[Te0]为3.0～10.0eV时等离子体的运动。结果表明,入口等离子体的内能增加,[B0]保持100G不变时,其最终膨胀的绝对速度增加,比冲从400s提高到约为580s;内能变化对比冲[Isp]的影响大于磁感应强度。不考虑等离子体与磁场相互作用情况下,文中计算的磁场范围可以最大限度地将内能转化为等离子体的轴向定向动能;为了提高[Isp],应适当增加电离段等离子体获得的能量,且可以适当降低对产生约束磁场的电流线圈输入能量要求。     

单/双涵道模式转换动态过程的数值研究

为研究变循环发动机单/双涵道模式转换动态过程中流场参数的变化规律,基于本项目组设计的模拟模式转换过程的试验台,利用动网格技术对该试验台进行了二维数值模拟。结果表明:在模式转换动态过程中,内外涵流量和总压的变化滞后于阀门缓慢关闭的状态;并且阀门变化速率越快,其内涵总压和流量的动态响应越滞后;而静压的变化规律则相反。此外,还研究了在阀门关闭过程中,不同初始涵道比对这种流场参数变化滞后性的影响。发现初始涵道比越小,滞后性越强;并且这种滞后性在阀门关闭至22°位置时达到最大。     

低温流体汽蚀的数值计算及可视化实验研究

为研究热效应对低温流体汽蚀的影响,选取文氏管为研究对象,进行数值计算及可视化实验研究。低温汽蚀计算模型考虑能量守恒的影响,物性参数随温度变化。使用高速摄影拍摄汽蚀发生位置及汽蚀区大小。热效应对常温水汽蚀的抑制很小为1.14%。液氮考虑热效应计算与实验结果误差为2.64%,吻合较好。由于低温流体的热物性,汽蚀发生伴有显著的温降及压降,汽蚀过程减弱,汽蚀区减小,最大体积分数下降20%,热效应对低温汽蚀的抑制作用明显。低温汽蚀受动量和能量交换的共同作用,液相、气相相互拖拽能力增强,汽蚀核心区由壁面扩展到整个流道;汽蚀区和气液界面的气泡浓度变化梯度更加平缓,气液界面不清晰。     

缺陷对高温合金蜂窝板弯曲力学性能的影响

采用预埋缺陷的方法,制备含有面芯脱焊缺陷的高温合金蜂窝板,进行了三点弯曲试验,研究了缺陷形状、大小、位置和蜂窝芯取向对蜂窝板三点弯曲损伤模式和承载能力的影响.研究发现,含面芯脱焊缺陷的蜂窝板在缺陷受压部位产生反向鼓包,并沿着宽度方向扩展为完全破坏;随着缺陷尺寸的增大,蜂窝板承载能力逐渐降低;当缺陷位于三点弯曲受拉面时,蜂窝板具有更高的三点弯曲极限载荷,且矩形缺陷试样的弯曲极限载荷要高于相同缺陷面积的圆形试样.最后,利用LS-DYNA对含圆形缺陷高温合金蜂窝板的三点弯曲性能进行了数值模拟分析,得到其三点弯曲过程中的应力分布状态.