导弹垂直发射燃气排导系统两相流场数值分析

根据燃气排导系统内两相流场的分布特点,利用可压缩雷诺 N-S平均方程对流场进行计算。湍流模型采用标准 k-ε二方程模型,采用颗粒轨道模型对颗粒相进行模拟,利用 FLUENT就发射高度及燃烧室压强对三维流场参数分布的影响进行了重点分析,为优化发射结构设计提供了重要的理论依据。     

飞机/进气道一体化流场数值计算

对飞机进气道组合体进行了粘性流场的数值模拟,通过求解N-S方程并建立合理的边界条件来模拟不同飞行速度(亚、超音速)、不同攻角和发动机流量下的进气道性能,给出了进气道进口局部马赫数、攻角和侧滑角的分布以及进气道的总压恢复情况。     

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

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

基于N-S方程的直升机旋翼翼型反设计方法

为了在实际旋翼翼型设计中考虑粘性影响,采用N-S方程为主控方程,建立了旋翼翼型流场求解方法。使用以Poisson方程为控制方程的网格生成程序生成围绕旋翼翼型的N-S粘性贴体网格,该网格生成程序能在反设计迭代过程中自动更新过渡翼型的粘性网格。在翼型流场计算及网格生成基础上,采用了MGM方程作为翼型反设计方程,建立了一套能够考虑粘性影响的直升机旋翼翼型的反设计方法。应用该方法,分别对二维NACA系列翼型、OA系列翼型、超临界翼型进行了反设计分析,获得了满足要求的二维旋翼翼型,并与目标值(压力分布、升力系数、阻力系数、力矩系数)吻合良好。     

外吹式动力吹气襟翼N-S方程数值分析

提出了一种带动力发动机计算模型,以N-S方程为控制方程,首先数值模拟了涡轮动力模拟器(TPS),进排气流场,验证了发动机计算模型的可靠性,然后利用该数值方法模拟并分析了外吹式动力吹气襟翼(EBF)流场特性.     

单转子风扇的三维反问题气动设计

采用准三维流函数反问题设计方法与三维 N-S求解方法的相互迭代对单转子风扇进行气动设计。在叶片初步计算中得到叶栅进出口气流角沿径向分布 ,并将它作为本文的目标函数。采用准三维反问题求解方法 ,依次构造出各个 S1流面上的叶片几何形状和气流角分布。然后再采用 N-S方程的求解方法 ,对叶片进行全三维流场的数值计算。通过 N-S方程计算结果与目标函数的对比 ,重新修正叶片出口气流角分布 ,并作为下一次反问题设计的目标函数。经过反问题与 N-S方程求解的反复迭代 ,最终得到满足设计要求的叶型。     

全机复杂流场网格生成及N-S方程数值模拟

全机复杂流场网格生成对计算结果起着举足轻重的作用。本文首先通过合理的拓扑结构,分块对接的方法对复杂流场进行了网格划分,然后采用N-S方程有限体积法及分区对接法对翼—身—尾全机流场进行数值模拟。计算结果验证了网格生成的有效性和计算结果的合理性。     

等速上仰翼型绕流结构的数值模拟

通过数值求解涡量-流函数形式的N-S方程研究等速上仰翼型的分离流动结构。用ADI方法解涡量方程;用Poisson方程直接法解流函数方程。在计算得到的流场中除可看到实验中观察到的多涡结构外,还发现一些实验中未观察到的现象。将分离流动结构的计算结果与翼型气动力变化相联系,并结合涡动力学理论,定性地解释了上仰翼型产生高升力的原因。     

应用多块对接结构网格方法的直升机涵道尾桨气动特性分析

基于湍流Navier-Stokes方程,建立了一套直升机涵道尾桨流场及气动特性的CFD分析方法.在该方法中,针对直升机涵道尾桨流场的特点以及动量源方法对网格系统的要求,提出了一种多块对接网格生成方法,分块贴体网格采用求解Poisson方程的方法生成.通过结合涵道尾桨的运动方式、几何特征及气动特征,建立一个包含动量源项的N-S方程的涵道尾桨流场计算方法和迭代流程;为较好地捕捉涵道壁附近存在的旋涡及分离现象,采用了一方程S-A湍流模型.通过对ONERA M6机翼的绕流特性的计算分析,以及对涵道尾桨在悬停状态下的拉力、诱导速度、载荷分布等的计算分析,验证了该CFD方法的有效性.在此基础上,采用该方法进行直升机悬停、侧飞和前飞状态下涵道尾桨流场与气动特性的数值分析,得到了关于涵道尾桨流场和气动特性的一些有意义的结论.     

旋转机翼悬停气动特性研究

鸭式旋转机翼(CRW)是一种先进的高速直升机方案,旋翼同时也是机翼,是其最关键气动部件之一。采用结构分区拼接网格技术进行空间离散,分区建立参考系,通过求解多重参考系下的N-S方程来计算旋翼流场,首先以传统的Caradonna-Tung实验旋翼的亚、跨声速悬停流场分析为例验证该方法的可靠性,进而采用该方法对旋转机翼悬停流场进行了数值计算,旋翼拉力计算值和地面实验值吻合较好,结果分析表明旋转机翼的悬停流场有着不同于传统旋翼的流场特性。     

三维副翼铰链力矩计算

 数值模拟了带副翼偏转构型的二维、三维亚、跨、超声速来流情况的绕流流场，采用有限体积方法和AUSM+迎风格式数值求解N-S方程，四步龙格-库塔时间推进，湍流模型为Baldwin-Lomax和Spalart-Allmaras湍流模型。二维模拟压强分布与实验数据以及CFD软件模拟结果吻合较好。采用交错对接网格，数值模拟了不同副翼偏角和缝隙、来流马赫数、迎角和雷诺数的三维构型流场。通过舵面铰链力矩与实验数据的对比分析表明了该数值方法模拟此类复杂流场的可靠性和舵面铰链力矩计算的有效性。发现并研究了亚声速来流时铰链力矩随迎角的反向变化趋势，初步分析了副翼的缝隙和雷诺数效应。     

物面突变引起的二维分离流动的数值模拟

本文用一个单步二阶精度格式,对由物体截面突变引起的二维超声速分离流动进行了数值模拟。计算了近尾迹流场,并详细研究了来流M数和Re数对绕后台阶分离流场结构的影响,得到了与理论和前人计算一致的结果。另外,对超声速绕竖舨流动进行了数值计算,给出了细致的流场结构图画。     

低雷诺数下微型飞行器主动变形机翼非定常气动特性数值分析

以微型飞行器的气动外形设计为研究背景,通过数值求解N-S方程分析了主动变形柔性机翼的气动特性。空间离散格式采用中心有限体积法,时间推进采用双时间推进方法,其中子迭代过程由隐式LU-SSOR方法完成。在动网格技术的基础上,模拟分析了薄翼面上作行波运动的流场。计算结果发现,该行波翼型的气动特性相比于刚性翼型有明显改善,有效地抑制了大攻角下大尺度流动分离,同时升阻比提高了约35%,起到了显著的增升减阻的作用。本文的工作对于进一步的柔性机翼优化具有良好的理论研究与数值模拟的参考价值。     

跨声速HTA涡轮导向器全三维流场计算分析

利用CFD软件CFX-TASCflow对具有试验数据的跨声速HTA涡轮导向器流场进行了计算，结果表明，计算的叶面等熵马赫数分布和出口气流角与试验值吻合良好；而计算的总压恢复系数与试验数据存在一定的差别。     

基于 N-S 方程的飞机风挡雨刷气动力计算

风挡雨刷在飞行过程中易出现被气流吹离复位位置、飘起现象,阻碍飞行员视野,严重影响飞行安全。选取某型机风挡雨刷作为优化对象,采用了结构、非结构混合网格的办法,在机身壁面附近生成结构化网格,在雨刷机构表面生成非结构化网格,用 Interface 面将两部分结合生成雨刷最终计算模型。利用 Fluent 软件,基于 N-S 方程和 S-A 湍流模型,采用时间平均法对不可压缩流的参数进行时均化,计算了不同飞行状态下,风挡雨刷在复位位置以不同角度停放及增加扰流板时所受气动力,对风挡雨刷的外形及停放角度进行优化,比较分析了优化前后雨刷所受气动力。仿真计算结果表明,雨刷合理的停放角度、在刷臂上加装扰流板能有效引导风挡雨刷迎风面气流的方向,从而减小雨刷迎风面所受气动力。通过试飞验证了数值计算的可靠性及优化方案的可行性,为风挡雨刷预紧力的设置及复位位置的确定提供依据。     

湍流模型对喷流计算精度的影响研究(英文)

采用三维RANS代码MBNS3D对包含了内流、外流及后体干扰的三种简化尾喷流,进行了数值模拟研究。通过壁面压力分布及喷流截面速度型与实验数据的对比,研究了SA和SST两种湍流模型及其修正(SACC,SADES,SSTCC)对喷流计算精度的影响。本文计算表明,SSTCC模拟效果最为精确;利用SA模型进行脱体涡模拟,在同样的网格下优势不明显。湍流模型可压缩修正改善了模拟喷流的膨胀趋势,但对附面层分离位置的预测,在本文的工作中,还不如原始模型。两类模型相比,SST模型的结果好于SA模型。     

某型民机着陆构型粘性气动力计算

采用ICEM对某型民机着陆构型进行多块结构化点对点网格剖分,运用基于N-S方程的WISE MAN PLUS软件对该构型低速粘性复杂流场进行数值模拟.在进场条件下,计算的升力和试验值误差在3.8%以内,俯仰力矩和试验误差在3.6%以内;计算与试验失速迎角差别在1°左右,最大升力系数较试验小6.3%.同时还给出了一些流场显示结果,对复杂流场进行了机理分析.计算与风洞试验结果符合较好,满足工程需要.     

旋成体大迎角分离流的Euler方程数值模拟

用三维Euler方程计算旋成体的大迎角绕流,我们发现:亚声速自由流时,计算得到的是无分离流动,与实际的分离流不相符合,因而计算的法向力大大小于有分离的实验值;而当自由流为超声速,且横向自由流M数较高时,由于背风表面出现激波,无粘Euler方程的数值计算能自动捕捉到背风面上的分离涡,计算法向力与实验值接近。为了对亚声速自由流用Euler方程数值模拟大迎角旋成体背风面上分离涡,根据问题的物理图像,在假设的分离线上,提出强加Kutta条件的方法,成功地计算出有分离流动,使计算的法向力与实验值接近。     

后向台阶横向喷流干扰流场的数值研究

 用有限差分法求解层流全N-S方程,对超音速后向台阶绕流与台阶后横向喷流干扰流场进行了数值模拟。不带喷流的超音速后向台阶绕流数值模拟结果与实验及其他计算结果吻合较好。最后给出横向喷流压力比p_j/p_∞=3.033～7.033和喷口宽度w=0.1h～0.2h(h为台阶高度)时的数值模拟结果。得到了合理的流场结构;并可看到喷流压力比和喷口宽度对流场特性有显著影响,在计算参数范围内,底部压力与系数Q(=wp_j/hp_∞)几乎成线性关系。