脉冲式燃烧风洞起动特性数值研究

利用计算流体力学软件Fluent,采用隐式双时间步方法模拟非定常过程,控制方程的离散采用二阶迎风格式,湍流模型选择剪切应力输运(SST)k-ω模型,并求解组分扩散方程,对脉冲式燃烧风洞的起动过程进行了数值模拟,分析了喷管起动时流场结构的变化情况和自由射流区超声速流场的建立过程.计算结果为改进风洞起动特性提供了一定的参考依据.      

受DBD等离子体控制的低速流动数值模拟方法研究

基于正离子、电子和氧负离子的三组分空气DBD等离子体动力学模型和流体动力学方程,利用有限体积方法处理控制方程空间项,采用双时间步方法模拟非定常过程,对DBD空气等离子体动力学模型方程和流体动力学方程松耦合求解,数值模拟了高压交流电源激励的低速空气等离子体流动过程。仔细研究了一个周期中真实时间步长对模拟DBD等离子体形成过程的影响,给出了等离子体的形成历程以及等离子体对低速翼型边界层流动的影响,结果在定性上是合理的。     

翼型扑翼运动流场的数值模拟

基于非定常N-S方程,利用动网格方法和预处理方法对翼型在低速度下扑翼运动进行数值模拟。采用Jameson中心格式的有限体积法、双时间推进法求解非定常N-S方程;紊流模型采用基于SA紊流模型的DES模型。本文计算模拟翼型在低雷诺数下扑翼运动过程,计算所得推力系数结果与实验参考结果吻合较好,对扑翼飞行器设计有一定参考价值。     

用Newton子迭代方法计算前飞旋翼粘性绕流

通过求解Navier-Stokes方程数值模拟了直升机旋翼前飞非定常流场，为了同时保证计算的时间精确性和计算效率，时间推进格式采用了双时间推进方法，在该方法中，子迭代过程由十分高效的LU-SSOR方法完成，且使子迭代过程成为Newton子迭代，空间上应用中心平均的有限体积法进行离散，为了模拟前飞桨时间的相对运动，网格布局采用了运动嵌套网格方法，应用本文方法对-悬停流场进行了数值计算，计算结果与实验吻合较好，尽管缺乏实验数据的验证，对-有升力前飞状态的数值模拟结果是可信的。     

一种用于分离流动的网格自适应算法

对于可压缩粘性流动,提出利用流场速度的紊乱度作为指示变量进行网格自适应.Jameson中心格式有限体积法、五步Runge-Kutta时间推进法/双时间推进法求解定常/非定常N-S方程.基于雷诺平均N-S方程模拟紊流,选用SA一方程模型.在数值求解二维静态失速和动态失速问题过程中,加入网格自适应算法,提高数值模拟对流动分离特性的捕捉和分辨能力.算例结果表明在流场发生失速后,运用本文的自适应算法能够在增加少量网格单元的情况下明显提高计算精度.     

舰载垂直发射系统内燃气流运动

通过对舰载垂直发射系统的物理模型进行简化及结合计算流体力学知识,采用数值模拟的方法对垂直发射导弹的发射过程进行了仿真分析,得到了发射过程中燃气排导系统内燃气流场的动态特性,为研究导弹发射过程提供了参考和依据。     

直升机旋翼悬停流场的欧拉方程计算

描述了三维Ｅｕｌｅｒ方程数值模拟旋翼悬停流场的过程,欧拉方程的求解采用了风格中心有限体积法、五步Ｒｕｎｇｅ－Ｋｕｔｔａ时间推进格式。     

螺旋桨非定常滑流数值模拟

运用非结构动态嵌套网格技术,通过求解三维非定常Euler方程,对螺旋桨飞机的螺旋桨滑流进行了真实数值模拟.动态嵌套网格方法通过建立运动区子域网格来模拟螺旋桨的旋转运动,静止的机身包含在静止区网格中,有效解决了螺旋桨滑流模拟中的相对运动问题.采用双时间方法求解控制方程,在每一个物理时间步上运用经典四步Runge-Kutt...     

失速限制型风轮叶片流场的数值模拟

采用双时间步方法在旋转坐标系下求解非定常方程,构造了零偏航条件下风轮叶片流场的数值模拟方法,采用该方法对失速限制型风轮叶片NREL Phase-Ⅵ进行了模拟,其中网格采用重叠网格进行模拟,控制方程为可压缩N-S方程,同时加入了低速预处理技术,模拟的风速范围从13m/s到25m/s,这时候风轮叶片发生失速以限制扭矩输出.本文采用数值方法对风轮叶片在发生失速以后的主要流场特点进行了分析,从风轮扭矩、载荷分布等方面对数值计算结果进行了对比验证.     

自适应直角切割网格民机增升装置绕流数值模拟

提出并介绍了基于直角切割网格的分区面搭接技术,采用变长宽比网格方法,成功地进行了存在外形间断(剪刀叉)的民机增升装置的网格生成和绕流Euler方程数值模拟.根据外形的特点,运用"根"网格分区算法,降低了整个网格的生成难度;通过基于外形的自适应网格加密技术,详细地描述了外形上的大量缝道和凹槽,提高了网格质量;在分区交界面上,利用重叠面积切割和面搭接算法实现了两侧网格间的流场信息传递,并保证了通量守恒;采用中心有限体积方法,结合双时间推进算法,完成流场的Euler方程数值模拟,计算结果与实验数据吻合良好,说明所述方法的正确性.     

带气腔的离心压缩机旋转失速的三维数值模拟

使用商业CFD计算软件CFX求解三维雷诺平均的Navier-Stokes方程组, 结合出口气腔模型对某带无叶扩压器的离心压缩机的旋转失速现象进行数值模拟.首先使用定常计算得到了该离心压缩机的稳态性能曲线, 并和实验测量值进行了比较.然后引入出口气腔模型, 模拟离心压缩机内的旋转失速流动.在小流量下, 从矢量图和径向速度等值线图中观察得到了离心压缩机内部流场的非定常流动现象.还研究了气腔模型不同参数对失速流动的影响, 结果表明气腔体积越大, 计算得到的失速频率越低.      

液体火箭发动机有无隔板的热态解数值模拟

用数值方法模拟了液体火箭发动机燃烧室内有无隔板两种情况的热态解。气相控制方程用欧拉坐标系下的Navier-Stokes方程组描述,液相控制方程在Lagrangian坐标系下进行描述。      

涡轮集气腔出流流量分配和温度特性研究

采用模型实验和数值方法,针对特定的涡轮集气腔结构进行了出流特性研究。首先,在实验室条件下,进行全尺寸基准涡轮集气腔的模型实验研究并验证数值模拟方法;其次,在接近发动机真实气-热参数条件下进行系列的数值模拟,分析进出口压比、腔室高度、出流孔孔径和出流孔轴向位置等对出流孔流量分配和温度分布的影响。研究表明,除了集气腔出流流量分配的不均匀性之外,气流与热壁面之间的对流换热也导致明显的出流温度不均匀性;改变出流孔轴向位置距离、增加集气腔腔室高度是在不影响集气腔通流能力下改善出流均匀性的有效技术措施,当出流孔轴向偏置距离或集气腔腔室高度达到一定值后,其对改善出流流量分配均匀性的作用趋于微弱,但对于改善出流温度分布均匀性依然有较好的作用。     

径向开槽药柱发动机内流场数值模拟

采用时间相关的显示ＭａｃＣｏｒｍａｃｋ预估－校正两步格式和可压缩全Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ方程，计算了轴对称径向开槽药柱发动机内流场。计算中采用了贴体坐标系统和Ｂａｌｄｗｉｎ－Ｌｏｍａｘ代数湍流模型。边界参数计算采用参考平面上的特征线法和物理边界条件。通过大量算例的分析表明，方法有效，计算结果合理。当有集中质量加入时，燃烧室头部压强升高，喷射强度和开槽口的倾斜方向对头部压强也有较大的影响。     

槽壁试验段低超声速流场特性数值模拟

 在跨声速风洞中通过开槽和抽气可以建立低超声速流场,由于槽壁试验段设计参数多,流场结构复杂,为提高设计准确性,通过数值模拟对槽壁试验段低超声速流场特性进行了研究。首先根据槽壁试验段的一般设计准则进行了气动设计,给出了槽壁尺寸和不同马赫数所对应的抽气量。基于设计结果,通过数值模拟对流场特性进行了研究,计算表明：通过抽气可以建立均匀的低超声速流场,抽气量对试验段马赫数均匀区长度有较大影响;随后对不同气动外形进行了比对,结果表明:抽气口位置、壁板厚度、驻室容积、开闭比及槽壁外形等对试验段的气流质量有影响,开闭比和槽型的影响尤为显著;最后对槽壁的通流特性进行了分析,探讨了槽型对试验段流场影响的成因。     

驻室式处理机匣——内部流动及扩稳机理的数值研究

通过定常模拟的方法,研究了驻室式处理机匣对高负荷、高马赫数孤立转子性能的影响.结果表明:相同的处理机匣结构在不同的转速下取得的效果存在较大差异,在设计转速时,处理机匣并没有使风扇转子的失速裕度得到改善,并且效率下降约2.8%;在非设计转速时,处理机匣扩稳作用得到充分体现,转子稳定性得到较大程度的加强,并且没有导致转子性能的下降.通过对处理机匣内部流场的详细分析,揭示了驻室式处理机匣的扩稳机理.      

处理机匣对风扇稳定性影响的数值研究

在对驻室式处理机匣内部流动及扩稳机理进行数值研究叫基础上,采用数值模拟方法进一步分析处理机匣对风扇性能的影响。通过取消处理槽上方的驻室,提高了处理机匣对设计转速的扩稳能力。根据叶尖间隙内的相对速度分布和处理机匣内的回流量,分析了处理机匣的扩稳特性。研究结果表明：处理机匣产生扩稳作用的原因在于处理机匣能使转子叶尖进口气流加速,增加流体能量,从而抑制由叶尖泄漏流发展成泄漏涡．讲而改善流动环境。     

直升机旋翼桨叶翼型积冰的数值模拟

本文对直升机旋翼桨叶翼型积冰进行了数值模拟.建立了空气、水滴两相流控制方程,通过求解气体欧拉方程得到空气流场,求解水滴控制方程得到翼型表面的水滴收集率.根据结冰表面的对流换热特性,使用积冰数值模拟的热力学模型确定翼型表面的积冰量,利用冰层时间推进法模拟积冰过程,采用积冰法向生长模型预测积冰形状;对NACA0012翼型在不同环境温度下的积冰进行了预测,最后对预测结果进行了分析.     

具有气液同轴喷嘴的液体火箭发动机燃烧室压力振荡及其声学特性研究

为探究液氧/煤油液体火箭发动机气液同轴喷嘴模型燃烧室具有良好稳定性的原因,采用非稳态雷诺平均（URANS）方法数值研究了其燃烧不稳定性和声学特征。两相燃烧条件下,燃烧室压力振荡幅值约为室压的10%左右、最大不超过25%,且以纵向和横向振型为主。一周六径隔板对横向振型具有很强的抑制作用,但对纵向振型影响较小。与液-液撞击式液氧/煤油发动机模型燃烧室相比,本文研究的燃烧室中煤油液滴没有发生超临界蒸发现象,第三邓克尔数较小,诱发燃烧不稳定性的激励源较弱。进一步通过数值定容弹激发了燃烧室多模态声学特征压力振荡,并得到了其振荡特征频率、幅值和衰减率。结果表明,气喷嘴具有四分之一波长喷嘴特征,能显著减小目标振型的幅值,而集气腔对纵向振型具有很强的抑制作用,同时对其他振型也有程度不同的抑制效果。因此,较弱的燃烧不稳定性激发机制以及隔板、气喷嘴和集气腔对纵向和横向振型很强的抑制作用,使得该液氧/煤油发动机气液同轴燃烧室具有很好的稳定性。     

固体火箭发动机内流场分区耦合数值计算

采用分区耦合方法计算固体火箭轴对称燃烧室与喷管流场。对于低速的燃烧室流场，选用不可压流的Ｎ－Ｓ方程描述并用ＳＩＭＰＬＥＣ方法数值求解；对于高Ｒｅ数的喷管流场，则采用Ｅｕ－ｌｅｒ方程描述并用ＳＣＭ方法求解。计算时用燃烧室出流为喷管流场提供入口参数，同时用喷管流场压强分布反馈影响燃烧室流动状况。对耦合边界条件处理方法进行了探讨。对典型的侧壁加质燃烧室与喷管流场进行了计算，计算结果揭示了单独喷管流场计算难以反映的喷管收敛段近壁区的低速区域，与已有的燃烧室流场实验结果一致并反映了燃烧室与喷管流场之间的联系，较好地模拟了流动中的物理现象。