斜网格修正的低速计算方法及应用

采用无坐标变换的一种有限体积法对RANS方程进行离散求解,离散方法考虑了对网格非正交情况的二阶修正,通过构造一个动量插值实现了同位网格下的压力修正算法,并发展了相应的分区算法.利用以上方法对轿车绕流和三峡坝区主流摆动问题进行了数值计算,计算结果与实验吻合较好.     

三维非结构网格的生成及优化

叙述了一种新的生成三维非结构网格的Delaunay方法,该方法能够高效地生成优质四面体网格。由于采用了局部剖分的方法来保证边界的完整性,因而不需要添加新的节点到原始表面三角形中。该方法在四面体外接球的球心处插入新的节点,在插入全部内点后,采用删除边、边/面交换、网格光顺等方法对生成的网格进行了优化,基本能删除质量差的四面体。     

浸入式贴体网格边界方法

提出一种基于重叠网格的浸入式边界方法,用以模拟复杂的三维黏性跨声速流动。该方法采用一套固定的笛卡儿正交网格,用以主流的求解;采用一套可以移动的贴体等距面网格,用以拟合或者离散物面的作用力,通过空间插值,实现两套网格重叠部分的信息传递。分析了浸入式贴体网格边界方法的优势,介绍了贡献单元的寻找策略和物理通量的插值方法。流场求解采用Spalart-Allmaras带湍流模型的Navier-Stokes方程组,其中对流项采用流通矢量分裂和5阶WENO(weighted essentially non-oscillatory)-Z格式离散,黏性项采用6阶中心差分格式离散,时间项采用龙格库塔显式格式离散。数值验证算例表明:该方法具备高于4阶的空间求解精度,并适用于刚体动网格非定常流场模拟,且无需更新网格形状。等距面贴体网格生成过程简单,避免了繁琐的人工调节过程,与笛卡儿网格结合,可提供足够的壁面附近网格密度,同时有效减少了网格总量需求。     

基于Euler方程结合附面层修正方法的应用研究

采用欧拉加附面层修正的方法在结构网格上对民用飞机翼身组合体流场进行了数值模拟研究.采用一种适合翼身组合体模型的的网格生成技术来生成计算用的结构网格,该技术是先采用一种基于Higenstock源项修正的椭圆型偏微分方程生成二维网格,然后沿机身轴向合理的光滑性调整成三维网格.在求解Euler方程时采用中心离散格式和显式时间推进,并利用附面层修正的方法来考虑粘性的影响.为了验证此方法的正确性,本文对DLR-F4翼身组合体进行了数值模拟,结果表明此方法和试验结果能够很好的吻合.     

应用TVD格式分离求解不可压N-S流动

提出一种新的将连续方程与动量方程分离求解的TVD引用方法。对动量方程中的对流项采用三阶精度、迎风TVD格式离散, 粘性项和压力项采用中心差分格式离散。压力采用压力替代法求解。计算所用的网格是具有交错网格主要特点的修正非交错网格, 并采用一种压力修正法和G-S中心对称迭代法加速收敛。给出了二维管流和二维空穴流动在不同状态下的计算结果, 并与其它解法的计算结果或精确解进行比较。结果表明, TVD格式不仅能用于不可压流数值计算, 而且还能得到比非TVD差分法精度更高的结果。      

非交错网格求解传热及流体流动

在非交错网格系统上提出了一种新的以控制容积为基础的有限差分方法，来求解传热及流体流动问题。为了避免压力棋盘问题，压力梯度采用偏迎风压力离散格式进行处理，以保持压力－速度的耦合。利用本算法对数值传热的典型问题进行计算，并将计算结果与采用交错网格的PHOENICS通用程序所得结果相比较，二者吻合甚好。     

面向脏几何的曲面网格自动生成算法

脏几何是存在法向错误、穿插、狭缝、贴合等问题的几何模型,在此类模型的原始边界上直接生成网格无法避免上述问题。因此,在脏几何上无法直接生成拓扑相容的表面网格,这也将影响后续数值分析研究的开展。本文实现了一种基于四面体化方法的脏几何网格生成算法,其核心思想在于一个有效的四面体网格的表面一定是一个有效的曲面网格。本文对几何的离散表征进行法向修复,并在此基础上生成四面体网格,通过离散边界切割四面体解决了模型中存在的穿插问题,基于四面体质量提升解决狭缝问题,并应用曲面网格重构以提升表面网格质量。本文采用视图采样的方法实现法向修复,从多个方向观察该模型,使其在各个方向上的视图均满足质量要求;使用多级的能量阈值提升四面体优化效率;采用多级误差估计来避免曲面网格重构中可能出现自相交的问题。数值验证结果表明本文的脏几何曲面网格自动生成算法鲁棒、高效。     

基于非结构动态网格的翼型积冰过程数值模拟

讨论了积冰气象条件下翼型积冰过程的非结构动态网格数值模拟方法.外部无黏绕流流场采用基于非结构网格的分布式并行计算技术,空间离散采用VanLeer迎风有限体积格式,时间推进为5步显式龙格-库塔方法.采用四阶龙格-库塔方法求解过冷水滴运动轨迹方程,使用了一种合理、简单的物面水积冰量计算方法,提高了计算效率.使用非结构动态网格技术对积冰后的翼型网格进行快速调整.在对NACA0012翼型积冰形状进行校核验证后,对某型飞机在积冰气象条件下飞行时翼型积冰过程进行了数值模拟,得到了与参考文献结果相符的冰形特性.     

基于离散曲率计算的三角网格模型优化调整

采用面积夹角加权的三角网格模型顶点法矢及三角片质心权值对Taubin的三角网格模型离散曲率计算方法进行了改进,在此基础上提出了一种新的三角网格模型优化调整方法。用该方法调整三角网格模型,在模型上曲率变化较平缓的平坦区域及曲率变化较剧烈的特征区域都能取得较好的调整效果。     

三维Euler方程的自适应非结构网络计算

介绍一种三维自适应非结构网格的Euler方程求解技术。采用有限体积中心差分的格点格式对Euler方程进行空间离散,多步龙格-库塔时间推进,并采用当地时间步长、焓阻尼修正等加速收敛。采用Delaunay三角剖分技术生成三维流场的四面体网格。网格自适应技术中的网格局部加密也采用Delaunay三角剖分技术。最后用ONERAM6机翼的跨音速绕流计算说明本文方法的有效性。     

基于径向基函数与混合背景网格的动态网格变形方法

发展了一种基于径向基函数(RBF)插值和混合背景网格映射的动态网格变形方法。混合背景网格由运动物面附近的各向异性网格和偏离运动物面的各向同性网格组成。物面变形或运动后,包含在各向异性背景网格中的所有物面网格点的位移由指定的物面运动规律得到,而各向同性背景网格点的位移通过径向基函数插值获得。然后变形后的计算网格利用计算网格与混合背景网格之间固定的映射关系代数插值得到。各向同性网格距离物面一定的距离,消除了物面保形对径向基函数插值的要求,从而可以减少径向基函数基点的数量,以提高插值的效率和可靠性。最后,通过二维和三维算例对该方法的变形网格质量、参数影响、网格拓扑依赖性、复杂外形变形能力及网格变形效率进行了测试。测试结果表明,该方法能够高效生成大变形下高质量的变形网格。     

三维非结构自适应多重网格技术

发展了一种基于几何外形的非结构自适应多重网格技术。根据流场参数的变化梯度确定加密边,由加密准则对粗网格进行自适应剖分得到分布合理的较密网格。通过预先生成的初始极密表面网格来将边界网格的加密点投影到边界上,使得边界保持初始外形。为了提高流场求解速度,将自适应加密后的一系列网格作为多重网格,采用了多重网格算法进行计算。M6机翼的流场计算表明该方法可以准确捕捉激波位置,并且表面的压力系数分布与实验值相吻合。     

基于混合网格的三维Navier-Stokes方程并行算法

提出了一种基于混合网格的三维Navier-Stokes方程的并行计算方法。Navier-Stokes的求解采用了基于面的有限体积方法,该方法适用于任何网格类型。采用一方程Spalart-Allmaras模型来计算紊流黏性。并行计算采用区域分裂的方法,利用METIS网格分区系统实现了各节点的加载平衡。节点间的数据交换通过调用MPI库函数来实现,采用非阻断通讯的方式来减少数据交换时间。充分利用FORTRAN90的动态存储特性来减少对内存的需求。最后,通过对绕DLR-F6外形(翼身组合体+挂架+发动机短舱)黏性流动的数值模拟,验证了该并行程序的准确性,高性能并行计算以及处理复杂几何外形的能力。     

分层嵌套重叠网格自适应树结构动态组装方法

采用重叠网格可以有效地进行复杂流动的大规模数值模拟,特别是包含运动部件(如旋翼、投弹)的动态模拟。本文将树结构的自适应直角网格用于重叠网格组装过程中的切割和贡献单元的搜索,大大加快重叠网格的组装速度。通过二叉树自适应直角网格对物体外形进行离散,实现切割过程的快速定位;采用八叉树自适应直角网格对流场区域进行离散,高效地搜索贡献单元。使用基于壁面距离准则的重叠区域最小化方法和分层嵌套重叠策略,能提高重叠网格组装的效率和质量。对于具有运动部件的动态重叠网格问题,采用多个二/八叉树减少组装过程中信息更新的冗余计算,从而大幅度减少重叠网格组装的时间消耗。实际算例的重叠网格组装结果说明本文发展的重叠网格组装方法具有很高的计算效率,可以满足运动边界复杂流动问题的动态计算要求。     

基于STL文件的笛卡尔网格生成方法研究

采用从STL文件中直接获取模型几何数据的方法,利用ADT数据结构快速剔除相同顶点,提出一种新的“推进查找”方法替代以往网格单元是否在物面内部的判定方法,避免了对物面的自封闭要求,无需进行几何重构,直接完成网格生成,更进一步提高了自适应笛卡尔网格生成速度和自动化程度,为实现CAD与CFD的之间的数据对接提供了一种新的方式。     

曲线物面边界条件在DGM中的应用研究

基于二维非结构网格,采用高精度间断Galerkin有限元方法数值求解定常Euler方程。为了有效克服传统的反射物面边界条件在物面处易于生成伪熵降低计算精度的缺陷,提出一种曲线物面边界条件,从而实现对物理模型而不是数值模型进行数值模拟。对经典圆柱亚音速无粘绕流进行数值模拟,结果表明:采用曲线物面边界条件之后,流场解的精度得到很好的提高;此外,在非常稀疏的网格上,通过提高基函数的阶次仍可以得到高精度的数值解。     

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....     

三维自适应非结构网格的Euler方程解

 将 Ausm+ 迎风格式应用于三维非结构网格中求解 Euler方程。对单元变量进行重构以获得空间高阶精度,对时间域采用多步龙格库塔法推进,并采用了当地时间步长和隐式残差光顺技术来加速收敛。采用多点择优推进阵面法生成复杂曲面的三角形网格,利用推进阵面法生成四面体网格。采用网格自适应技术对网格进行局部加密,以减少总体网格数目,从而提高计算效率。最后给出了绕 ONERA M6机翼的跨音速流动及绕麻雀 导弹的超音速流动算例,结果表明了本方法的有效性。     

非结构化网格生成中背景网格对源的响应

背景网格对源的响应方法对非结构网格的生成具有重要意义,如果背景网格参数能恰当地反映出源的变化,则所生成的非结构网格就能满足网格设计的要求。将背景网格的三角形单元分为7个子域,根据源在单元中的位置对网格进行局部加密或变形,使每一个源点都与一个背景网格结点相对应,再通过求解热传导方程,使背景网格参数分布光滑,基本反映网格设计的要求,从而提高了非结构化网格的生成质量。     

A parameterized geometry design method for inward turning inlet compatible waverider

Intensive studies have been carried out on generations of waverider geometry and hypersonic inlet geometry. However, integration efforts of waverider and related air-intake system are restricted majorly around the X43A-like or conical flow field induced configuration, which adopts mainly the two-dimensional air-breathing technology and limits the judicious visions of developing new aerodynamic profiles for hypersonic designers. A novel design approach for integrating the inward turning inlet with the traditional parameterized waverider is proposed. The proposed method is an alternative means to produce a compatible configuration by linking the off-the-shelf results on both traditional waverider techniques and inward turning inlet techniques. A series of geometry generations and optimization solutions is proposed to enhance the lift-to-drag ratio. A quantitative but efficient aerodynamic performance evaluation approach （the hypersonic flow panel method） with lower computational cost is employed to play the role of objective function for opti- mization purpose. The produced geometry compatibility with a computational fluid dynamics （CFD） solver is also verified for detailed flow field investigation. Optimization results and other numerical validations are obtained for the feasibility demonstration of the proposed method.