一种抑制激波-边界层相互作用的新型无源被动控制

为了探索一种抑制激波-边界层相互作用的新型无源、被动控制方法,采用有限体积法计算流体力学程序,分别对压缩拐角和激波入射平板两种典型流动的附面层分离进行了无源被动控制的数值研究.数值结果表明,新型无源被动控制方法消除了20°压缩拐角流动所产生的分离泡;对激波入射平板所诱导的附面层分离,采用的自适应无源控制方法可将分离区长度减小为无控制时的58％、总压恢复可比无控制时提高7.2％、x=0.85m截面最大回流速度较无控制时减小了69.3％.从而验证了无源被动抑制控制激波-边界层相互作用的可行性.     

减小翼型激波阻力的鼓包流动控制技术

针对2020年使用的N+2代民用飞机的翼身融合(BWB)布局发展需要,以减小激波阻力为目标,采用计算流体力学(CFD)方法,开展弱化激波、减小激波阻力的鼓包流动控制技术研究.提出了λ形激波结构“强干扰”和等熵压缩“弱干扰”两种鼓包激波减阻流动控制原理,给出了两种鼓包基本形状设计方法和工程应用的可行性分析,指出λ形激波结...     

双锥Bump压缩面设计及气动特性

采用反设计方法研究双锥Bump压缩面设计技术,该方法实质是融合运动间断边界在运动网格条件下对轴对称欧拉方程组求解.为避免激波捕获法对激波型面位置求解误差,应用计算域分块的方法预估第2道激波,在此基础上采用流线追踪法生成Bump压缩面;结合平面机身设计一双锥Bump实例,运用计算流体力学仿真手段对其进行黏性数值模拟.研究结果表明:①该模型流场结构仍保持较强的附面层扫掠能力;②在来流马赫数为2.0条件下,相比传统的正圆锥乘波体Bump压缩面设计,新型设计方法可使Bump外压缩系统总压恢复系数提高0.04左右,为Bump进气道性能提高奠定基础.      

RANS/LES在超声速突起物绕流中的应用研究

安装在超声速/高超声速飞行器表面的突起物如机翼、控制舵等通常会导致复杂的激波/边界层干扰,对突起物的局部气动特性甚至飞行器整体的气动特性产生较大的扰动.在采用计算流体力学(CFD)数值模拟此类问题时,传统的求解雷诺平均Navier-Stokes(RANS)方程方法由于不能准确预测湍流脉动流场并且精度有限,在应用上受到一定的限制.本文在研究B-L (Baldwin-Lomax)内层模型和Smagorinsky亚格子模型优缺点的基础上,提出了一种新型的RANS/LES(Large Eddy Simulation)混合模型,并进行了算例验证,证实了该方法的可行性.在此基础上,对火箭表面突起物的干扰流场进行了数值模拟研究,细致地刻画了突起物附近的激波/边界层干扰、剪切层失稳和底部分离涡形成的非定常过程,获得了突起物及火箭表面上的压力脉动历程并进行了频谱分析.研究发现,相对于突起物底部的非定常分离流动,突起物前缘的激波和边界层相互干扰的非定常过程是突起物周围压力脉动的主导因素,这种高频的压力脉动可能对火箭内设备的正常工作产生不利的影响.     

激波结构内流动问题的气体运动论描述

从介观Boltzmann速度分布函数理论出发,发展计及分子粘性碰撞截面与扩散碰撞截面,可描述各流域一维气体流动问题的Boltzmann简化速度分布函数方程及其气体运动论数值计算方法。通过对不同Knudsen数下非定常激波管流动及不同马赫数定常正激波结构问题数值模拟,研究分析不同流区的激波突跃变化过程以及近连续流、稀薄过渡流特有的分子输运现象,揭示不同马赫数、不同分子模型的激波内流动与传热变化规律,证实基于Bo-ltzmann模型方程的气体运动论数值计算方法用于激波结构内流动研究的准确可靠性。     

基于本征正交分解和代理模型的流场预测方法

 为了实现流场的快速求解,基于本征正交分解(POD)与代理模型提出了一种全新的流场预测方法。其基本原理为:首先采用本征正交分解将一定数量的样本流场分解为同等数量的基模态流场;然后用少数包含了绝大部分样本流场特征的基模态流场拟合所有的样本流场;最后用代理模型建立起决定样本流场的输入参数与拟合系数之间的近似函数关系。针对几何外形不同的二维翼型定常流场预测结果表明:在亚声速情况下,预测误差收敛的模态数量不超过20个,继续增加模态的使用数量不能明显提高预测精度;在跨声速情况下,预测误差收敛的模态数量为26个,当使用的模态数量达到前10个时,继续增加模态的使用数量能提高绝大部分流场区域的预测精度,但同时会在激波附近引入"噪声"激波特征而降低该区域的局部预测精度。在这两种情况下,预测流场所需时间均不到高精度计算流体力学(CFD)方法的1/200。     

一种基于聚类分析的二维激波模式识别算法

在激波捕捉求解器计算的可压缩无黏流场基础上,提出了一种探测并识别二维激波干扰模式的新算法,从3个层面详细介绍了该算法的实施流程。首先,采用基于当地流场参数设计的传统激波探测方法,辨识出激波附近的一系列网格单元;其次,通过经典的K-means聚类算法将这些激波单元划分成许多簇,并根据簇的相邻信息定义每个簇的类别;最后,设定相关准则对某些紧邻的簇进行合并,进而确定各个激波干扰点的位置,记录各条激波分支所对应的簇,采用Bézier曲线拟合算法分别对其聚类中心进行拟合以获取更加光滑的激波线。数值试验表明,该算法不受网格类型的限制,不仅可以保证最终拟合的激波线具有较高的位置精度,还可以清晰地识别出流场中多激波干扰的模式,同时对分析非定常流场中激波的运动与演化过程也提供了一种有效的可视化手段。     

Numerical simulation of the three-dimensional flow-path in a ram-rotor

A kind of three-dimensional flow-path of the ram-rotor was designed in this paper according to the design method of the typical supersonic intake of aircraft and missiles,especially the design method of the compression ramp,the throat and the subsonic diffuser.The design Mach number of the flow-path was based on the averaged relative Mach number,regardless of the change of incoming Mach number along radial direction of the flow-path.Then,the flowfield of the designed flow-path was simulated numerically with computational fluid dynamics software at design point.The purpose of simulation was to obtain the details of the flowfield,including the structure of the shock wave system,position of flow separation zone,and so on.The performance of the ram-rotor was also evaluated.The numerical results show the structure of the shock wave system in the flow-path is similar to that in the supersonic intake.To improve the overall performance of the ram-rotor,the reasonable compromise between the total pressure ratio and the isentropic adiabatic efficiency should be reached.The structure optimization of the flow-path and the configuration of the shock wave system have great influence on the performance of the ram-rotor.      

High-enthalpy hypersonic project at ONERA

This article deals with the research performed in the framework of the internal ONERA Hyperenthalpic Hypersonic Project. This project involved fifteen researchers specialised in different domains, which are the following: fundamental studies on shock/wave and shock/boundary layers interactions; study of the laminar/turbulent transition; real gas effects and associated physical modelling; the development of non-intrusive optical diagnostic methods usable in cold and hot wind tunnels; flow rebuilding in a hot shot wind tunnel such as F4; conception, realisation and validation of computational fluid dynamics (CFD) solvers.     

联合使用激波捕获—拟合法求解激波间断的跨声速流场

吸取激波捕获和激波拟合两种方法各自的优点,发展了一种求解有激波的跨声速流场的方法,并对回转叶栅中定常跨声速流场做了计算。计算中,通过捕获法确定初始激波位置,然后经反复修正,拟合出确定的流场通道激波。计算表明,这种方法可自动获得清晰的激波,流场中气流各参数分布合理。这一方法可适用于复杂边界和不同进口M数,计算时间仅比势函数方法多一倍左右。     

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

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

超燃冲压发动机唇口气动热计算研究与分析

采用数值模拟方法,对超燃冲压发动机进气道唇口的气动热进行了计算和分析.依据唇口与前体激波所处相对位置的不同,分三种情形进行了计算研究.结果表明,三种情形下,进气道前缘驻点热流均大于头部.通过算例验证,证明该方法和结论的正确性.研究表明,激波干扰以及前体激波压缩后超声速来流使唇口热流增大,热防护和结构设计时应当予以充分考虑.      

变截面管道中激波诱导的两相流动

本文研究扩张型与收缩型变截面管道中激波诱导的非定常两相流动。对于稀相气固悬浮体采用了双连续介质模型,对于两相流动按照准一维近似处理。控制方程则是利用算子分裂技术和二阶GRP方法数值求解。文中给出了波后气固两相的流动结构,并讨论了粒子对两相流场的影响。     

高超声速计算中的气体动理学格式

回顾了高超声速连续流部分的计算流体力学(CFD)方法,总结了近些年兴起的气体动理学格式。阐述了该格式的构造机制,强调了将物理规律直接用于构造数值方法的思路。结合一些应用实例,例如激波相互作用、激波边界层相互作用以及边界层分离等高超声速问题,说明了这种构造思路给数值模拟带来的优点。从高超声速的发展历程来看,气体动理学格式的构造过程包含了更基础的物理规律,而且具有多尺度的特性。这些特性有助于研究复杂的高超声速问题。介观或者微观角度直接构造数值方法的发展趋势为高超声速计算工具指出了可能的发展方向。     

基于OpenFOAM的三维H2/Air连续旋转爆轰流场数值模拟

为进一步研究旋转爆轰流场特征,基于开源计算流体动力学软件OpenFOAM,采用9组分19步的基元化学反应模型,对H2/Air连续旋转爆轰流场进行了三维数值模拟,得到了旋转爆轰波稳定传播时燃烧室内部流场的详细结构,研究了燃烧室头部激波的传播特性,分析了旋转爆轰燃烧室的压力增益性能。结果表明:旋转爆轰波后的第一道反射激波在由燃烧室外壁面向内壁面传播过程中反射激波的高度增加并在靠近内壁面附近与滑移线交汇形成局部高温高压区域;旋转爆轰波在外壁面位置处相位约落后于内壁面0.003rad～0.15rad,其相位差随燃烧室曲率差的增加而增大;燃烧室头部反射激波数目受到曲率差和进气总压的影响,燃烧室曲率差增大,反射激波数目减少,进气总压增大,反射激波数目增多;燃烧室压力增益保持在0.3以上,在进气总压一定的条件下,压力增益随着燃烧室曲率差的增大有增加的趋势。研究结果揭示了三维旋转爆轰流场的精细结构和燃烧室头部激波的传播规律。     

激波-边界层-分离流相互干扰三维湍流的数值模拟

本文采用数值方法求解时间相关三维可压缩雷诺平均Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ方程组,模拟激波—边界层—分离流相互干扰三维湍流流动。湍流模型为Ｂａｄｗｉｎ－Ｌｏｍａｘ两层代数模型,改进后用于三维内流问题。采用单元中心有限体积法离散流场控制方程,ＶａｎＬｅｅｒ矢通量格式计算无粘通量,中心差分法计算粘性通量,ＬＵＳＧＳ时间推进格式计算定常流场。本文以二元跨音速扩压器内三流动为算例,数值模拟较强激波—边界层—分离流相互干扰维湍流流动,并与实验结果进行了比较。数值模拟结果,在激波强度、分离点位置和再附点位置等方面,与实验结果吻合较好。     

翼/短舱/吊架气动特性及其动力影响数值研究

应用计算流体动力学(CFD)数值模拟方法对DLR-F6模型分别在通流及动力条件下翼/短舱/吊架(W-N-P)组合体部分的气动干扰特性进行研究.结果表明:组合体相互干扰引起的干扰通道附增激波及短舱后部逆压导致的涡流增加了短舱外罩阻力,降低了机翼升力;动力条件较通流条件,F6短舱外罩阻力降低,机翼升阻力特性在零度攻角时有所提高.      

定后掠角密切锥乘波体的生成和设计方法

对定后掠角密切锥乘波体（OCWRCAS）的生成方法和考虑黏性的设计方法进行了研究。定后掠角乘波体的前缘具有特定的后掠角,能够在上表面产生稳定分离涡从而改善乘波体的气动性能。本文首先在传统密切锥乘波体生成方法的基础上给出了定后掠角密切锥乘波体的生成方法;从前缘后掠的几何特征中提取了后掠角、激波角和前缘曲线程度等设计变量,并研究了设计变量的取值范围;以遍历设计空间的思路对两类定后掠角密切锥乘波体进行了设计分析,研究了升阻比、体积效率随设计变量的变化规律,然后在设计空间内进行了多目标寻优;最后使用计算流体力学方法对定后掠角乘波体的乘波特性和涡升力特性进行了验证。结果表明,由本文生成方法得到的定后掠角密切锥乘波体具有明显的乘波特性并且能够在较高的升阻比时保证一定的体积效率;定后掠角前缘能够在一定的迎角下在上表面产生稳定的分离涡,产生涡升力。     

SFD方法模拟高超声速进气道中激波与边界层干扰

采用ＳＦＤ方法,在ＥＸＣＥＬ软件上利用ＮＮＤ格式数值模拟了激波与平板边界层之间的相互作用以及激波与管道壁面边界层的相互作用,即高超声速进气道隔离段内激波串流动现象。该方法的特点是简单、直观而且不用常规编程,计算结果基本反映了流场的基本结构,表明ＳＦＤ方法是一种具有潜力的计算方法。     

Evaluation of wave drag reduction by flow control

An analytical expression is proposed to estimate the wave drag of an aerofoil equipped with shock control. The analysis extends the conventional approach for a single normal shock wave in the absence of control, based on the knowledge that all types of successful shock control on transonic aerofoils cause bifurcated λ-shock structures. The influence of surface curvature on the λ-shock structure has been taken into account. The extended method has been found to produce fairly good agreement with the results obtained by CFD methods while requiring negligible computational effort. This new formulation is expected to be beneficial in the industrial design process of transonic aerofoils and wings where a large number of computational simulations have to be performed.