用三维N—S方程计算钝体绕流

本文采用有限体积方法求解三维低速Ｎ－Ｓ方程，研究了三维钝体和类汽车物体的绕流。用中心差分和迎风差分的混合差分格式构造代数方程组，采用了隐式和显式相结合的方法使主对角元素占优以提高迭代求解过程的稳定性。利用压力校正算法，采用面松驰迭代来实现方程的求解，湍流的计算采用了标准的双方程模型，本文给出了钝体和类小汽车外形的计算结果，钝体的计算结果和实验结果进行了对比，与其它方法的计算结果也进行了对比，它们之     

用三维N－S方程计算钝体绕流

本文采用有限体积方法求解三维低速Ｎ－Ｓ方程，研究了三维钝体和类汽车物体的绕流。用中心差分和迎风差分的混合差分格式构造代数方程组，采用了隐式和显式相结合的方法使主对角元素占优以提高迭代求解过程的稳定性。利用压力校正算法，采用面松弛迭代来实现方程的求解。湍流的计算采用了标准的双方程模型。本文给出了钝体和类小汽车外形的计算结果。钝体的计算结果和实验结果进行了对比，与其它方法的计算结果也进行了对比，它们之间具有合理的一致性。     

高超音速化学非平衡钝劈绕流数值计算

 自1970年Davis提出粘性激波层方法以来,用数值方法求解高超音速轴对称钝体绕流问题国内外已做了大量工作,但至今仍未见到关于平面问题的计算结果。对于平面问题,虽然方程在形式上比轴对称简单,但由于二维效应,激波层较厚,用文献[2]的方法向下游区推进有困难。另外,在驻点线上采用极限关系式虽能克服方程的奇性,但驻点解对流向步长Δξ有依赖。     

二维非均匀射流遇障碍时波系的变化

本文研究了非均匀射流在二维管道内的钝头体绕流问题。所采用的Godunov格式适合于不规则区域上双曲守恒系统的流场计算,数值计算的结果给出该流场的主要特征和激波结构,对于认识这一物理现象和了解波的运动规律是有帮助的。     

气液两相绕流钝体稳燃器的冷态数值研究

提出并建立了基于可变时间间隔平均方法的多尺度两相湍流模型,通过算例验证,证实了所推导的气液两相多尺度数学模型、所标定的模化参数及描述气液两相作用机理的合理性,及其计算的时间平均阻力系数、脉动升力系数误差、斯特劳哈尔数和回流区大小误差分别为1.45%,0.323%,2.17%,2.33%,计算结果明显优于标准k-ε和重整化群(RNG)k-ε模型.用多尺度湍流模型对气液两相绕流6个不同的钝体稳燃器进行数值计算表明:船形与锥形是综合性能较优的两个钝体稳燃器结构.其中,船形钝体比锥形钝体稳燃器回流区大10.53%,其时间平均阻力系数比锥形稳燃器大4.776%,方均根脉动升力系数比锥形稳燃器小44.73%,通过全方面综合比较,船形是综合性能最优的钝体稳燃器结构.     

强旋-V型钝体组合燃烧装置绕流实验研究

为了探明强旋流对钝体绕流旋涡脱落规律的影响,利用PIV速度场仪和高频动态压力传感器等测量仪器研究了强旋-V型钝体组合稳燃装置中旋流器出口安装位置对旋涡脱落频率和尾流流场的影响规律。实验结果表明,当旋流器出口位置安装在钝体内部时,可获得较小的旋涡脱落频率。随着旋流器安装位置向钝体出口方向逐渐移动时,旋涡脱落频率增大。当旋流器出口置于钝体外部时,将破坏钝体剪切旋涡的几何特征,从而导致钝体旋涡脱落机制的改变。另外,当旋流体积流量控制在一定范围(0~1m~3/h)内,增大旋流流量有助于减小钝体绕流剪切旋涡的脱落频率,对比无旋流时的V型钝体绕流,当旋流流量为1m~3/h时,旋流器出口位置位于钝体尾缘内部1倍开口直径处的V型钝体绕流St数下降幅度最大,约为34.96%。一定范围内的旋流同样有利于扩大回流区,在旋流流量为0.5m~3/h时回流区长度达到最大,约为钝体开口直径的2.3倍。但当旋流体积流量过大时,将导致下游流动恶化,钝体的火焰稳定性能可能会急剧下降。     

方柱-强旋组合旋涡脱落机制研究

强旋流和钝体绕流均可构造低速回流区以利于高速气体稳定燃烧,本文通过在钝体中心开孔通入旋流气体模拟研究了旋流流动对钝体绕流旋涡脱落的影响规律。计算结果表明,在旋流作用下,钝体绕流涡脱规律更加复杂,涡脱由上下边界层失稳主导变为由旋流管中的旋流涡主导,并且形成的涡尺度增大,在脱落过程中会由于旋涡的破碎产生"点潭"以增强周围涡强度。同时,柱体涡和旋流涡的相互作用可以减小流动阻力,增加回流区长度和回流质量流量,从而强化传热传质,起到强化燃烧增加燃烧稳定性的效果。     

化学非平衡边界层绕流计算

本文介绍用化学非平衡边界层方程组,数值计算高超声速钝头体绕流流场、组元浓度和电子密度分布情况。计算中考虑了电离空气的七种主要成分以及相互间的化学反应,比较了变熵、普朗特数及Lewis数的影响,讨论了完全催化壁和非催化壁两种情况,计算结果与文献及试验进行了比较,结果令人满意。     

机翼跨声速绕流的Galerkin有限元解法

本文给出了翼型和机翼跨声速绕流的Galerkin带权余数有限元解法。利用一系列适当形状的元素将计算区域离散,在元素内采用线性插值函数。对于超临界流动,文中采用了适合于有限元法的超声速上风技术,具有捕捉激波的能力。将这种技术应用于线松弛迭代解法中,在计算钝头翼型和机翼绕流时,都获得了成功。     

高超声速球锥粘性激波层压力分裂法数值求解

本文将压力分裂法引入细长钝体高超声速粘性激波层数值求解中,较好地解决了这种流场远后身部分的推进求解问题。针对非解析球锥体绕流数值求解的困难,作者提出了实用型初始激波形状计算方法,达到了简化数值求解步骤,提高数值计算效率的目的,数值求解结果与实验结果相符合。     

二维超音速湍流混合特性研究

着重讨论顺流喷射情况下的超音速混合的流场特征及影响因素。简要介绍了顺流喷射混合的研究现状。同时，探讨了不同的来流参数对流场特征的影响及混合机理。计算采用了一种二阶ＴＶＤ格式有限体积法及Ｂ－Ｌ湍流模型求解二维湍流Ｎ－Ｓ方程。结果表明：所采用的计算方法能成功地抑制非物理振荡，适用于二维超音速内流的湍流混合特性研究。     

渐变反推力平面叶栅端壁流动仿真与试验

针对某型反推力装置中的渐变反推力平面叶栅,开展了性能预测和流动试验的研究与验证。结合数值仿真和性能试验,运用表面油流等测量手段,探索了反推力叶栅流线分布和流动的基本结构,揭示由不同弯角组成的反推力叶栅旋涡结构和流动特征。试验与仿真结果表明:反推力叶栅端壁呈现出3维非定常有旋的流动特征,其发展和演变的过程与涡轮叶栅中的流动过程相似,依据仿真结果和流道收敛的特征,确立了反推力平面叶栅的基本测量方法,并通过了试验结果的验证,对今后的反推力叶栅数值仿真和试验工作有很大的借鉴意义。     

突发事件下大规模空中交通流量管理的组合优化模型

为了解决危险天气和军航活动对管制运行的影响，使用动态网络流方法对突发事件下短期空中交通流量调度问题展开研究。首先，结合中国民航管制的特点，分析了突发事件对流量管理的影响。同时，根据航路航线网络及其高度层的特点，给出了网络和高度层的数学描述，并根据机型将航空器分成大、中和小3种交通流，介绍了使用多品种流描述3种机型的必要性；其次，根据网络拥挤程度随时间和流量变化的特点、危险天气随机变化的特点、空中等待和地面等待费用差异的特点构建了3个优化目标，考虑机场容量、扇区容量、航班连续性和扇区连续性约束建立了多目标优化模型；再次，针对军航活动时需要协调空域的特点，多品种流模型求解时间复杂度高，不能适应短期流量管理的缺陷，改进了逐步宽容约束法，设计了一种阶段式求解的近似算法；最后，以西南空管局管制的空域为例，利用实际流量数据，设计了3个场景进行仿真。结果表明，所提模型和算法能有效求解突发事件下的短期流量调度问题，算法效率比起传统算法在大流量下更具优势。     

火箭级间分离喷流干扰数值模拟与风洞试验研究

针对典型风洞试验状态,采用有限体积法求解三维雷诺平均Navier-Stokes(RANS)方程,数值模拟了给定级间距下火箭级间分离喷流干扰流动情况,获得了清晰的流场结构和有/无喷流干扰时两级气动特性变化规律,并对仅主机喷流和主机与四游机同时喷流的情况进行了计算,数值模拟研究结果为风洞试验方案的设计提供了参考依据。在典型状态数值模拟研究的基础上采用了冷喷流模拟技术对多个状态下火箭级间分离喷流干扰特性进行了风洞试验研究,获得了两级在有/无喷流、同轴变攻角情况下的气动力系数,研究结果表明,级间喷流干扰会对两级气动特性带来较大影响。相同状态下计算得到的两级气动特性变化与试验符合较好,表明了针对风洞试验状态的预先数值模拟对试验方案的设计验证可提供有益的参考。     

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

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

有翼飞行器大攻角分离涡模拟

基于势流方程,采用涡格法和三维离散涡相结合,模拟具有升力面分离涡的大攻角弹冀、弹体、尾翼组合体绕流。该方法适用于高亚音速、小攻角或大攻角有翼飞行器的压力分布、法向力和压力中心等气动特性的数值计算,并具有显示三维涡流场的能力。     

内埋式弹舱流场特性及武器分离特性改进措施

为改善内埋弹舱的流场特性以及内埋武器分离特性,采用在弹舱前缘悬细金属条的方法对弹舱流场进行流动控制,并在高速风洞中进行了试验研究。通过分析舱底静态压力试验结果以及脉动压力试验结果,研究了武器模型处于不同分离位置时流动控制对弹舱流场特性的影响;通过测量武器模型力和力矩,研究了流动控制对武器分离特性的影响。研究结果表明：武器模型处于不同分离位置时,该流动控制方法对弹舱底部静态压力分布以及总声压级分布的影响是相似的;当弹舱的流场类型为过渡／闭式穴流动时,采用该流动控制措施能有效降低舱内的静态压力梯度,并能有效改善武器的分离特性;当弹舱的流场类型为开式穴流动时,采用该流动控制措施能有效抑制舱内产生的气动噪声。     

一种多级轴流压气机特性预测工程方法

基于多级轴流压气机部件设计点气动参数与相应几何尺寸,发展了一种计算发动机高、低压压气机部件特性近似方法, 并将其应用于某型涡扇/涡喷发动机风扇和高、低压压气机部件特性预测,计算结果与实验数据比较表明,本方法作为一种近似方法具有可接受的工程精度,尤其适合于预测发动机压气机部件小转速特性,为涡扇/涡喷发动机起动过程分析与计算提供模型基础。     

复杂旋翼流场的耦合欧拉-拉格朗日数值方法

针对影响旋翼流场求解精度的关键因素“桨叶复杂近体流动”和“尾迹涡畸变”,结合计算流体力学(CFD)方法和黏性涡方法,发展了一套适合于复杂旋翼涡流场分析的耦合欧拉-拉格朗日数值方法:为捕捉桨尖三维效应、激波等细节流场特征,在桨叶近体区域采用CFD方法对其进行求解;针对高雷诺数旋翼流场中桨尖涡的紧凑结构特点,引入黏性涡方法建立了高分辨率的尾迹求解模型;两计算域间的信息交换采用了集中涡源法和边界修正法.应用所建立的计算方法,以旋翼CFD标准验证试验(Caradonna-Tung旋翼)为算例,对尾迹影响明显的悬停状态进行了数值模拟,通过对比耦合边界处流场特征及桨叶表面压力系数分布,验证了方法的有效性.此外,还从旋翼尾迹捕捉精度、涡量耗散特征及计算时间等方面对不同计算方法进行了对比分析,结果表明耦合方法可充分发挥CFD和黏性涡方法各自的优点,在旋翼流场数值模拟方面具有独特的优势.     

航空发动机1 次表面换热器流动换热性能分析

为了研究航空发动机1次表面换热器流动换热性能,基于传热单元数法并结合其结构特性,建立了换热器热力学设计方法并对经典热力学公式进行了对比分析。同时针对适用于航空发动机的4种不同结构形式1次表面换热器(直通道逆流型和1 5°,30°及45°叉流型),在真实工况下的流动换热特性开展了数值模拟研究。通过对比不同结构换热器在不同工况下的流动换热特点,可以为一次表面换热器芯体核心部件的优化设计提供设计依据和方法。基于数值计算结果,对比分析了不同交错角度θ对换热器的换热性能与流动特性的影响。结果表明:对于直通道逆流换热器,整个换热器内部温度有规律均匀分布;对于叉流换热器,由于波纹板片呈一定角度交替放置,内部流动复杂,局部存在明显的涡流强化换热,气体流动通道内的速度、温度分布极不均匀。随着交错角度的不断增大,叉流换热器的换热性能不断增强,但其冷热两侧压降也大幅增大。