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

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

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

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

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

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

风力机远尾流的计算研究

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

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

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

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

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

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

为了研究凹坑形非光滑表面的减阻效果与流动控制机理,将半球形凹坑单元布置在直升机尾舱门外表面。采用基于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.     

基于混合笛卡儿网格方法的可压流动数值模拟

以可压缩黏性流动的数值模拟为研究背景,发展了一套自适应混合笛卡儿网格(AHCG)方法以及基于有限体积方法的雷诺数平均Navier-Stokes(RANS)的数值求解方法.为更好地模拟边界层的黏性流动在近壁面处采用贴体结构网格,剩余计算区域自动生成与之相重叠的笛卡儿网格,并同时发展了基于流场特征的笛卡儿网格自适应技术.结合ADT(alternating digital tree)算法显著减少了网格生成中“挖洞”和“贡献单元”搜索的消耗机时,50万左右的网格数目下,搜索耗时为0.062s,仅为普通遍历方法的1/1847.通过二维圆柱与两段翼型绕流的数值算例显示,定常AHCG方法能够准确地预测物面压力分布与升阻力系数并且具备处理复杂外形的能力;同时通过二维非定常圆柱绕流问题与NACA0015矩形机翼翼尖尾涡的捕捉算例显示,结合了动态自适应网格加密的非定常AHCG方法尤其适用于旋涡主导流动.      

基于混合笛卡尔网格的翼型绕流数值模拟

在对重叠网格技术研究的基础上,一套以笛卡尔网格为主体、贴体的结构网格为特殊固壁边界的混合笛卡尔网格(Hybrid Cartesian Grid,HCG)生成方法得到了发展,同时基于流场特征的笛卡尔网格自适应技术和有限体积方法的雷诺平均Navier-Stokes(RANS)的数值求解方法也得到了开发。通过对经典跨音速RAE2822翼型、双NACA0012翼型和高升力两段NLR-7301翼型进行绕流数值模拟,获得了较好的结果,验证了方法的可行性与精确性,为研究二维复杂外形流动问题提供了一种解决方法。