低速两段翼型的数值模拟

在多段翼型风洞实验中,很难观察翼型不同迎角下主翼、襟翼上的绕流及缝道流动。使用了O-H型混合的结构化网格,采用k-ε二方程湍流模型求解可压的N—S方程来模拟二元风洞中两段翼型的流场。计算结果与实验进行了比较,结果表明：计算与实验结果吻合良好,说明本方法可以较好地模拟两段翼型的绕流,并能很好地显示流场的变化情况。另外,从结果中可以看出：随着迎角的变化,主翼尾流和缝道间的流动对襟翼附面层分离有很大的影响作用。     

二维非定常欧拉方程的有限元解法

本文通过求解时间相关的欧拉方程，对振荡翼型的非定常二维无粘可压缩绕流作了数值模拟。求解方程时采用了有限体积积分的方法，并用ＶａｎＬｅｅｒ逆风通量分裂格式来计算控制体通量。隐式格式的采用，大大提高了计算效率；同时，通过插值使得空间离散达到了二阶精度。另外，本文采用基于非结构三角形网格的动态网格来处理刚性物体的移动，从而顺利地解决了复杂几何外形物体的非定常绕流的问题。文章最后给出了ＮＡＣＡ００１２翼型绕四分之一弦线作简谐俯仰振动的计算结果。一个周期内的瞬时压强分布以及升力和力矩系数都与实验结果吻合良好，     

高超声速验证飞行器助推分离段流场数值研究

为了模拟有相对运动的多体非定常绕流，发展了基于弹簧近似的非结构动网格方法及耦合动网格的Euler解算器，作为方法验证，模拟了俯爷振动的NACA0012翼型绕流，计算结果与实验及文献结果非常接近。采用二维近似外形，对高超声速验证习行器助推分离段流场进行了数值研究，得到了不同时刻由于多体相对运动形成的干扰流场结构以及分离过程的气动力参数。     

Gurney Flap绕流气动特性研究

采用C-H型多块结构网格和SST k-ω湍流模型求解二维定常雷诺平均NS方程,对NACA4412单段翼型后缘附近和NACA63(2)-215B多段翼型襟翼后缘附近的Gurney Flap绕流进行了数值模拟。讨论了局部网格对计算结果的影响,重点研究了Gurney Flap位置和高度对气动特性的影响。计算结果表明:对于NACA4412单段翼型,Gurney Flap位置对气动性能的影响较小,但随着高度的增加,零升攻角和失速攻角减小,最大升力系数增大,同时阻力系数和力矩系数增加;对于NACA63(2)-215B多段翼型,Gurney Flap对气动性能的影响类似于NACA4412单段翼型情形。     

网格拓扑和湍流模型参数对多段翼流动的影响

计算网格和湍流模型是影响CFD结果不确定的重要因素。基于NHLP-2D高升力翼型开展了数值试验对比研究,采用多块分区方法研究了适用于二维复杂构型的网格拓扑结构和网格分布方法,并在尽量排除网格影响因素的情况下,采用均匀设计方法初步分析了SST模型参数在低速多段翼型绕流流场中的影响大小。研究成果在工程应用中具有一定的参考价值。     

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

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

多重网格技术在复杂粘性流场计算中的应用及研究

流体力学控制方程的数值求解过程中,当网格加密、粘性效应强或来流马赫数低的时候,流场收敛非常缓慢。本文在多块网格、LUSGS隐式格式基础上采用多重网格技术,计算了NLR7301case#1两段翼型的低速绕流及DLRF4跨声速流场,在复杂二维、三维粘性流场求解中获得了满意的加速收敛效果。文中还比较了V循环、W循环的加速效率,为多重网格的工程应用奠下基础。     

非结构聚合多重网格法流场数值模拟研究

 在非结构网格上应用聚合多重网格技术来加速 Euler方程的收敛过程。设计了一种高效率的网格聚合方法。采用四重网格,在每一层网格上采用有限体积法进行计算。通过对多段翼型绕流的数值求解验证了多重网格加速收敛的高效性。     

多段翼型大迎角分离流动的DELAYED RANS/LES混合算法

多段翼型的大迎角绕流发生大范围附面层分离,具有明显的三维与非定常流动特性。RANS/LES混合算法继承了LES对流动分离区大尺度漩涡准确模拟的优点,避免了纯LES算法需求网格量巨大与亚格子模型壁面函数不成熟等问题,对分离流动的模拟效果优于RANS算法。以S-A湍流模型与Smagorinsky亚格子模型为基础,借鉴DDES的附面层延迟控制思想,构造了可用于对接网格、重叠网格的DELAYED RANS/LES混合算法。研究了GA（W）-1多段翼型的大迎角分离流动及其气动特性。     

基于LBM-LES方法的多段翼型流场数值模拟

运用Lattice Boltzmann(LB)方法与Large Eddy Simulation(LES)结合的方法数值模拟研究了二维翼型绕流流场,其中,粒子速度模型采用D2Q9模型,松弛时间中对应的湍流粘性系数由Smagorinsky模型计算,网格生成方法采用无网格技术。对NACA0012翼型进行数值模拟,验证其有效性;为进一步研究其模拟复杂物体绕流的能力,以三段翼型30P30N为研究对象进行数值研究,计算出了给定攻角条件下的气动特性,并与实验值进行比较。计算结果表明,LBM-LES方法的计算结果能够与实验值一致,是一种可行、有效且较新的数值模拟方法。     

动网格技术在增升装置优化中的应用

将基于Denaulay背景图的动网格方法应用于增升装置二维缝道参数优化设计过程中。针对二维增升装置的构型特点和翼型运动的特点,采用了添加辅助点的双层Denaulay背景图动网格方法,并详细论述了辅助点添加、运动技术。为了进一步提高动态网格质量,提出了针对前后缘运动位置中比较极端的两种位置生成两个初始网格,并采用两种初始网格进行动态网格的生成。结果表明新动态网格方法可以为二维增升装置提供高质量的动态网格。     

增升状态下的全机位流数值模拟

研究了势流理论下的全机增升装置复杂外形气动力数值计算方法,该方法以低阶面元法为基础,具有计算量小、速度快的优点.复杂流动的粘性效应主要通过附面层粘性修正技术加以体现.采用该方法对某增升装置翼身组合体外形压力分布特性进行计算,并和风洞实验数据相比较,验证了方法的正确性和实用性.在此基础上对某型飞机真实复杂外形的压力分布特性进行了计算,为复杂外形飞行器增升装置的初步设计、方案选型与评估提供了非常理想的工程技术手段.     

低速增升装置失速特性预测研究

民机低速增升装置失速特性是具有挑战性的空气动力学难题之一,决定着飞机起、降阶段机翼处于流动分离或尾迹干扰状态时的操纵稳定性,是气动设计必须关注的重要方面。利用基于N-S方程的自研多块结构网格计算软件和混合网格计算软件,对NASA高升力梯形机翼标准模型开展了失速特性计算研究,系统评估了计算软件、计算网格、湍流模型等主要因素对失速特性的影响,获得了有效预测高升力构型失速特性的计算方案,并将其成功应用于工程实际,为增升装置失速特性计算评估提供了有意义的参考。     

An automatic isotropic/anisotropic hybrid grid generation technique for viscous flow simulations based on an artificial neural network

Based on the author’s previous research, a novel hybrid grid generation technique is developed by introducing an Artificial Neural Network(ANN) approach for realistic viscous flow simulations. An initial hybrid grid over a typical geometry with anisotropic quadrilaterals in the boundary layer and isotropic triangles in the off-body region is generated by the classical mesh generation method to train two ANNs on how to predict the advancing direction of the new point and to control the grid size....     

A 3D hybrid grid generation technique and a multigrid/parallel algorithm based on anisotropic agglomeration approach

A hybrid grid generation technique and a multigrid/parallel algorithm are presented in this paper for turbulence flow simulations over three-dimensional (3D) complex geometries. The hybrid grid generation technique is based on an agglomeration method of anisotropic tetrahedrons. Firstly, the complex computational domain is covered by pure tetrahedral grids, in which anisotropic tetrahedrons are adopted to discrete the boundary layer and isotropic tetrahedrons in the outer field. Then, the anisotropic tetrahedrons in the boundary layer are agglomerated to generate prismatic grids. The agglomeration method can improve the grid quality in boundary layer and reduce the grid quantity to enhance the numerical accuracy and efficiency. In order to accelerate the convergence history, a multigrid/parallel algorithm is developed also based on anisotropic agglomeration approach. The numerical results demonstrate the excellent accelerating capability of this multigrid method.     

各向异性网格自适应计算在超燃模拟中的应用

采用各向异性网格自适应求解技术,将其应用于DLR超燃冲压发动机燃烧室中的超燃模拟.开展了3个算例,包括采用滑移壁面条件的燃烧室冷流场模拟、采用无滑移壁面条件的冷流场模拟及采用无滑移壁面条件的反应流模拟.模拟中,各向异性网格自适应计算捕捉到了如激波、射流、边界层、火焰面等具有各向异性特征的大梯度区域,并利用各向异性网格进行了很好的加密.对比利用各向同性网格的初始流场计算,各向异性网格自适应计算使基于滑移条件、无滑移条件的冷流计算及反应流计算的网格单元数量分别下降了36.2%,36.4%和36.8%,有效降低了计算规模,而且流场大梯度区域的计算结果更准确,辨析度更好.结果表明:对于像超燃这类具有各向异性特征的问题,各向异性网格系统比各向同性网格系统有更好的计算效率及准确性,同时也表明基于Mach数场构造的各向异性网格系统可以有效应用于超燃计算.      

流动参数对合成射流控制叶栅流动分离的影响

采用大涡模拟方法、结构化网格建立了低压高负荷透平Pak B叶栅的非稳态数值分析模型,研究了不同流动参数对合成射流控制叶栅流动分离的影响.控制前随着雷诺数的减小和气流攻角的增大,叶栅流动分离区域变大,在气流攻角为5°下发生分离未在尾缘前再附的情况.合成射流控制后,不同流动参数下的流动分离都得到了有效的控制,并且在射流偏角为30°时,合成射流控制效果最好.合成射流使叶栅吸力面的流动分离位置推迟,再附位置前移,分离泡尺寸减小,叶栅吸力面的逆压梯度段缩短,吸力面边界层表面的剪切层在向下游迁移的过程中,没有发生充分的抬升,避免了大尺度涡旋的形成,并且很快地黏附于壁面,进而有效地控制了流动分离.      

国家数值风洞(NNW)工程中的黏性自适应笛卡尔网格方法研究进展

相比传统的结构和非结构网格,笛卡尔网格具有自动化和高质量兼顾的优势,是未来网格技术发展的重要方向。面向笛卡尔网格生成软件开发需求,依托国家数值风洞（NNW）工程,介绍了黏性自适应笛卡尔网格技术相关课题的研究工作进展,重点从网格生成技术、自适应方法、黏性物面边界处理方法等方面展开。网格生成技术方面,基于全线程树数据结构的思想,开展了叉树数据结构的优化;针对物面离散单元信息的快速检索,构建了稳定、平衡、高效的叉树式物面网格数据结构;对于网格单元类型的判断,发展了染色方法和包围盒方法以提高效率;针对大规模网格量,发展了高性能并行计算技术。自适应技术方面,发展了基于几何特征和流场解特征的自适应方法,建立了流场结构捕捉判据,并开展了三维复杂构型的应用;为了降低全局网格规模,发展了各向异性自适应方法,构建了各向异性特征识别的综合判据;发展了边界层法向射线加密技术,以减少边界层内的网格数目。物面边界处理方面,发展了重叠笛卡尔网格方法,构建了重叠网格尺度匹配技术和数据交互方法;开展了壁面函数方法在笛卡尔网格框架下的适应性应用,以放宽壁面网格尺度限制、降低边界层网格规模,并通过典型算例验证了方法的可靠性。     

低压涡轮叶型边界层相互作用的数值模拟

基于Lantry-Menter转捩模型,应用商用流体计算软件求解非定常雷诺平均N-S方程组,对进口雷诺数为2.5×104,来流湍流度为2.5%,尾迹折合频率为0.92状态下的尾迹/PAK-B叶型边界层的相互作用进行了数值模拟。数值计算揭示的尾迹在叶栅通道中的输运特性、尾迹诱导卷起涡形成、尾迹诱导边界层转捩等物理机制大致符合相关实验定性描述。     

Mesh deformation on 3D complex configurations using multistep radial basis functions interpolation

The Radial Basis Function (RBF) method with data reduction is an effective way to perform mesh deformation. However, for large deformations on meshes of complex aerodynamic configurations, the efficiency of the RBF mesh deformation method still needs to be further improved to fulfill the demand of practical application. To achieve this goal, a multistep RBF method based on a multilevel subspace RBF algorithm is presented to further improve the efficiency of the mesh deformation method in this research. A whole deformation is divided into a series of steps, and the supporting radius is adjusted in accordance with the maximal displacement error. Furthermore, parallel computing is applied to the interpolation to enhance the efficiency. Typical deformation problems of the NASA Common Research Model (CRM) configuration, the DLR-F6 wing-body-nacelle-pylon configuration, and the DLR-F11 high-lift configuration are tested to verify the feasibility of this method. Test results show that the presented multistep RBF mesh deformation method is efficient and robust in dealing with large deformation problems over complex geometries.