可压缩粘性流动笛卡尔网格虚拟单元方法研究

以二维高雷诺数可压缩粘性流动问题为背景,提出了一种全新的笛卡尔网格虚拟单元方法。基于壁面函数基本假设,构造了壁面函数-虚拟单元方法(WF-GCM),用于定义湍流壁面边界条件。引入参考点的概念计算虚拟单元上的基本变量与湍流变量值,定义了"非贴体"笛卡尔网格下的湍流壁面边界条件,并通过壁面函数模型修正近壁面单元与界面单元。基于自适应笛卡尔网格体系,采用发展的具有二阶精度的格心格式有限体积求解器,数值模拟了跨音速RAE2822翼型绕流问题与超音速圆柱绕流问题,计算结果与实验值吻合良好,显示了WF-GCM对高雷诺数可压缩粘性问题是有效的。     

外挂武器对无人直升机纵向气动特性影响研究

针对外挂武器对无人直升机纵向气动特性影响的问题,采用雷诺时均Navier-Stokes(Reynolds average Navier-Stokes,RANS)方法对加装武器系统前后,不同前飞速度、武器安装角和挂载状态下的气动特性进行了数值计算。然后将大速度前飞状态时的气动特性与加装外挂武器之前的风洞试验结果进行了对比分析。结果表明,采用悬臂梁外挂方式加装武器对阻力有显著影响,前飞速度和武器发射安装角变化对纵向气动特性几乎没有影响。外挂武器安装位置和数量变化对无人直升机的纵向气动特性影响很小。研究结果可为武装无人直升机选择合适的武器外挂和发射方式提供参考。     

壁面函数在激波诱导分离流动中的应用

以数值模拟激波-附面层干扰引起的流动分离问题为研究背景,发展了基于有限体积方法的雷诺平均Navier—Stokes（RANS）方程的流场数值模拟方法。利用壁面函数模型得到壁面剪切应力,通过修正壁面粘性通量,构造了一种新的湍流边界处理方法,并将其耦合到RANS方程和SSTk-ω湍流模型的数值求解中;同时,针对激波诱导引起的附面层流动分离问题,提出一种附面层网格加密技术,能够自适应加密分离区内附面层网格,使得在流动分离区域也能够使用壁面函数模型。数值算例表明,壁面函数模型能够降低数值模拟结果对网格的依赖性;同时也验证了壁面函数耦合附面层网格自适应方法,在处理激波诱导引起的附面层流动分离问题时的有效性和准确性。     

二维多段翼型网格生成与数值模拟研究

为提高多段翼型的网格生成效率和数值模拟精度,发展了一套自适应混合笛卡尔网格(AHCG)生成方法和基于有限体积方法的雷诺数平均Navier-Stokes (RANS)的数值求解技术.混合笛卡尔网格由围绕物体几何外形的贴体结构网格和填充流场其他区域的笛卡尔网格构成,两套网格之间的信息传递由“贡献单元”提供,且“贡献单元”由基于ADT(Alternating digital tree)技术的搜寻方法获得.为更准确地捕捉流场信息,采用了基于流场特征的网格自适应技术.数值模拟结果显示,AHCG方法能够准确且高效地模拟高升力多段翼型绕流问题.     

基于有限体积格式的自适应笛卡尔网格虚拟单元方法及其应用

基于自适应笛卡尔网格,结合虚拟单元方法处理浸入边界,发展了一种有限体积格式实现模式,数值模拟了二维无粘可压缩流动。本文发展的方法相对于传统的有限差分方法具有执行效率快,易于获得守恒性等优点。最后,利用一些标准算例和多物体流场的计算问题验证了算法的有效性。     

风力机远尾流的计算研究

基于工程上广泛使用的Park模型构造了一种新的远场尾流计算模型,同时又提出了一种CFD与制动盘理论相结合的方法对单个风力机远尾迹区的流动状况进行了数值模拟.新的尾流计算模型是由一维Park模型进行三角函数修正得到的,CFD方法是结合制动盘理论使用FLUENT提供的Fan边界数值求解整个尾流流场.给出了单台风力机尾流中风轮中心的轴向速度分布以及风力机下游不同位置处的径向速度分布,并与风洞实验结果以及Park模型的计算结果进行分析对比,表明本文的两种计算方法都能够较好地捕捉尾流的流场特性.由于新的远场尾流计算模型较之原Park模型能够更准确地反映径向尾流分布,且计算量非常小,因此可作为工程开发的有效工具.制动盘理论结合CFD的方法也能够在计算量相对较小的情况下,准确的反映远尾流的流场信息,这两种方法都可为风场微观选址和风力机排布提供参考依据.     

基于非光滑表面的直升机尾舱门减阻技术

为了研究凹坑形非光滑表面的减阻效果与流动控制机理,将半球形凹坑单元布置在直升机尾舱门外表面。采用基于k-ω剪切应力输运模型(k-ω shear-stress-transport,k-ω SST)的计算流体力学(Computational fluid dynamics,CFD)方法,对比分析了非光滑表面机身和光滑表面机身的阻力特性、后体局部流场和表面压力分布。根据凹坑单元内局部流场和压力分布的分析结果,得到了凹坑形非光滑表面的减阻机理:凹坑单元内形成的低速漩涡,能够延缓机身后体流动分离;半球形的凹坑造型和基本恒定的压力分布,可减小对阻力直接做出贡献的低压区面积,两者的共同作用降低了直升机阻力。     

ADAPTIVE HYBRID CARTESIAN GRID METHOD FOR VORTEX-DOMINATED FLOWS

An efficient compressible Euler equation solver for vortex-dominated flows is presented based on the a- daptive hybrid Cartesian mesh and vortex identifying method. For most traditional grid-based Euler solvers, the excessive numerical dissipation is the great obstruction for vortex capturing or tracking problems. A vortex identif- ying method based on the curl of velocity is used to identify the vortex in flow field. Moreover, a dynamic adaptive mesh refinement （DAMR） process for hybrid Cartesian gird system is employed to track and preserve vortex. To validate the proposed method, a single compressible vortex convection flow is involved to test the accuracy and effi- ciency of DAMR process. Additionally, the vortex-dominated flow is investigated by the method. The obtained re- sults are shown as a good agreement with the previous published data.     

多机起降旋翼间气动干扰数值模拟研究

以多架直升机近距离起降为研究背景,采用基于运动嵌套网格技术的计算流体力学(Computational fluid dynamics,CFD)方法开展旋翼间的气动干扰研究,重点分析了旋翼间距和环境风速对旋翼升力系数和扭矩系数的影响。数值分析结果表明,旋翼之间的气动干扰对上风位置旋翼的影响较小,但会导致下风位置旋翼的升力系数降低,扭矩系数增加,且其影响随着旋翼间距的增加而降低,随着环境风速的增加而增加。     

SST-DDES模型在大分离流动问题中的应用

采用基于SST k-ω湍流模型耦合延迟脱体涡模拟技术(Delayed detached-eddy simulation,DDES)的SSTDDES模型,数值求解Navier-Stokes方程。并针对低速流动问题,发展了基于隐式双时间步方法的低速预处理技术,用于数值模拟低速、高雷诺数湍流大分离流动问题。数值模拟了低速大迎角三角翼绕流及大迎角6∶1椭球绕流,观察到与物理现象相一致的旋涡特征,且得到的定量数值结果与实验数据相吻合。