气相爆轰波马赫反射非自相似性特征的实验

为了研究在CJ(Chapman-Jouguet)状态下气相爆轰波马赫反射的非自相似性特征,搭建了由驱动段、传播段以及观测段组成的矩形管道;利用纹影和烟膜实验方法分别对3种预混气发生马赫反射现象进行了实验研究.实验结果表明:对稳态或非稳态气体,测量了马赫反射三波点位置高度,验证了马赫反射三波点轨迹线为波动的曲线,即自相似性对爆轰波马赫反射失效;非自相似性的重要特征是三波点轨迹线在楔角初始位置时基本遵循无反应冲击波理论计算结果,随后逐渐偏离并向有反应冲击波转变,之后接近平行于反应冲击波理论计算并稳定地处于两个理论值之间.最后给出在初始压力为10kPa和楔角为30°的条件下,对于所研究的3种气体,其马赫反射三波点发生转变的位置高度分别约为0.8,1.05cm以及0.5cm,可见相同初始条件下非稳态气体的马赫反射发生转变比稳态气体提前.      

单边扩张燃烧室燃烧非对称及非稳态现象研究

为了探究单边扩张超声速燃烧室在不同当量比下出现的多种模态及其内在机理,通过实验和数值模拟,对该构型燃烧室在并联凹腔对称喷注的工况下进行了系统研究。实验中,隔离段来流入口马赫数为3.46,来流总温为1486K;采用乙烯作为燃料,当量比范围从0.26到0.51。结果表明:燃烧室工况随着当量比的上升而改变,且可分为三个模态。当量比较低时(0.26),形成只有剪切层内有火焰的对称燃烧。当量比较高时(0.51),形成上凹腔为凹腔内燃烧、下凹腔为剪切层燃烧的非对称燃烧。当量比适中时(0.4),火焰在上下凹腔间来回震荡。在非稳态过程中,预燃激波串沿燃烧室轴向做周期性往复运动,定量分析表明火焰震荡的频率约为170Hz。燃烧震荡来源于分离区喷注与燃烧的耦合效应,震荡机制与涡无关,可能是由声振机制引起的凹腔自激震荡。     

动网格技术在增升装置优化中的应用

将基于Denaulay背景图的动网格方法应用于增升装置二维缝道参数优化设计过程中。针对二维增升装置的构型特点和翼型运动的特点,采用了添加辅助点的双层Denaulay背景图动网格方法,并详细论述了辅助点添加、运动技术。为了进一步提高动态网格质量,提出了针对前后缘运动位置中比较极端的两种位置生成两个初始网格,并采用两种初始网格进行动态网格的生成。结果表明新动态网格方法可以为二维增升装置提供高质量的动态网格。     

一种新的涡扇发动机动力特性计算模型

以等熵流动原理为出发点,结合发动机短舱入口、出口及远场特征边界条件,建立了一种新的涡扇发动机动力特性计算模型.在基于非结构网格的并行计算软件中实现了模型,对轴对称超高涵道比涡扇模型(CRUF)、轴 、对称涡轮动力模型(TPS)以及某型民用飞机整机带发动机模型进排气流场进行了数值模拟,计算结果与传统的发动机动力特性模型计算结果以及实验数据进行比较,三者吻合较好,从而验证了构造的涡扇发动机动力特性计算模型的正确性、可靠性以及工程实用性.     

欧洲计算流体力学——个人观点

介绍了70年代初到20世纪末欧洲科学家对CFD发展的主要贡献,并重点关注了对国际流体力学数值方法会议(ICNMFD)和国际计算力学学术会议(ISCFD)的成功举办具有重大贡献的科学家.ICNMFD始于1969年(新西伯利亚),由参与了美国和苏联太空竞赛的两国科学家联合创办,而ISCFD则作为ICNMFD在亚太地区的同类会议初创于1995年(东京).由Oshima教授和本人提议,得益于上述两个会议全体科学委员会成员的共同努力,ICNMFD和ISCFD自2000年开始,合并成为每两年举办一次的国际计算流体力学会议(ICCFD).绝大部分本文讨论所涉及的学者通常都在受邀进行大会报告后,成为了前述三大国际会议中的一个或多个科学委员会的成员.作者对文中讨论所涉及的研究人员的圈定负有全权责任,上述人员的选择更多地基于应用而非理论计算流体力学,因此尚存在不完善之处.此外由于文中所讨论的学者大都是作者在40多年计算流体力学研究生涯中的私交好友,因此并没有将年轻一代的科学家考虑进来,但这并非由于他们的贡献不重要.另外值得注意的是欧洲的计算流体力学同美国的计算流体力学有非常紧密的联系,本文讨论所涉及的研究学者有些最早从欧洲开始涉足计算流体力学而现在则生活和工作于美国(例如Glowinski,Roe,Van Leer等人),而有些则是在美国取得了博士学位(例如Launder,Napolitano),他们都同美国的科学家进行过合作,因此很难清晰地对欧洲和美国的计算流体力学进行界定.对那些现已移居美国的欧洲科学家们,本文优先考虑了他们移居之前的贡献.最后但并非是最不重要的,作者通过电子邮件一一联系了文中所涉及的科学家,希望他们本人能够给出其自身贡献的评价草稿,从而帮助作者以最权威的视角在文中相应章节对其工作进行评述.本文向提供自述评价的诸位致以谢意,同时也希望其他同仁能够对本人不恰当的擅自评述给予谅解.     

4116航天润滑油摩擦特性的影响因素

采用球盘试验机测定了4116航天润滑油在不同工况下的摩擦因数,尤其成功实现了高低温环境下的测试.应用正交试验设计方法,研究了滚动速度、滑动速度、最大工作压力和入口油温4个因素对润滑油摩擦特性的影响,并绘出各因素对摩擦因数的影响曲线.研究结果表明:滑动速度和滚动速度对摩擦因数有显著影响,入口油温和最大工作压力对摩擦因数影响不显著.摩擦因数随着滑动速度的增加而增加,但是增加的幅度随滑动速度的增加而逐渐减小;随着滚动速度的增加,摩擦因数迅速减小,当滚速高于20m/s时趋于稳定;入口油温低于0℃时,其对摩擦因数有较大影响;随最大工作压力的增加摩擦因数缓慢增加.      

W型地下井发射环境数值模拟与分析

为研究导弹地下井发射环境,采用计算流体力学方法对W型井导弹发射过程的流场进行了数值模拟。通过将仿真数据与理论计算及试验结果进行对比,验证了计算模型和数值方法的准确性。采用动网格技术处理流场计算区域随导弹运动的变化过程。仿真结果表明,在整个发射过程中燃气流具有反射、引射和充分发展3个阶段,造成发射环境参数有较大变化。分析了压力和温度参数的变化规律。该数值算法和仿真结果可推广应用于发射井的设计、优化及改进。     

用非结构直角网格和欧拉方程计算运载火箭绕流

采用八叉树结构，物面附近用投影法，生成运载火箭绕流计算的非结构直角网格。并将Jameson的有限体积法推广用于这种网格的欧拉方程计算。算例表明，网格生成的时机化费很少，质量好；欧拉方程解算的收敛质量也好。     

固体推进剂非稳态燃速测量及燃速对压强变化响应滞后

本文扼要介绍了一个能直接、快速地测量固体推进剂瞬时燃速的激光系统，并用于测量压强变化条件下的瞬时燃速响应滞后。实验结果显示，在压强振荡、降压和升压过程中，瞬时燃速对压强变化响应都有滞后。滞后时间（ｔ＊）随压强变化速率增加而减少。但ｔ＊－ｄｐ／ｄｔ曲线不是一条直线，而是一条近似的双曲线     

三维热化学非平衡流动非结构网格DSMC方法及其应用

研究了过渡流域三维热化学非平衡流动DSMC方法实现的过程。以四面体非结构网格为基本单元，提出一种新型的高效搜索算法，该算法不仅可以快速跟踪模拟分子在网格之间的迁移，而且可以准确判别分子与物面是否相互作用，避免了原有算法中分子表面反射的非确定论判据。设计了适合三维DSMC方法的动态局部时间步长技术，节约了计算时间。将碰撞距离的思想引人到非结构网格上来，有效地消除了网格尺度小于三分之一平均自由程的限制。利用Fortran90的动态分配内存技术编制了适用于任意外形的通用计算程序。最后对全尺寸航天飞机高超声速绕流进行了数值试验，计算结果验证了该算法的可行性及高效性。