有限元三维网格自动生成

讨论了为提高有限元法求解流动问题的精度和效率而寻求自动生成三维网格的方法。给出了利用代数变换、输入极少量区域边界数据即可插值生成三维网格的方法。最后举出了自动生成三维网格的例子。     

Mesh generation: Art or science?

Mesh generation has evolved to the point where highly complicated domains can be covered by a variety of mesh types including hexahedral, tetrahedral and overset meshes. The application of these methods to computational aerodynamics has become a routine exercise and numerical predictions over complete aircraft now complement experimental results obtained from wind tunnels. This paper surveys the main developments that have taken place and traces the evolution of mesh generation over the last 35 years. This is followed by an assessement of the accuracy of Navier Stokes codes that are currently in use for predicting the drag of an aircraft at transonic cruise. The relationship between solution accuracy, mesh size and mesh type is examined in some detail and the implications for further research are discussed.     

NLR-7301高升力翼型绕流评估计算

探讨了在二维高升力多段翼型绕流计算中附面层网格采用各向异性三角形网格的技术。通过对经典的NLR-7301两段翼型绕流的计算,从数值计算角度评估了自研混合网格雷诺平均N-S方程求解软件对该类型附面层网格的解算能力。通过与确认试验数据对比,在附面层采用各向异性三角形网格进行二维高升力多段翼型绕流计算,其结果具有较高的可信度。在此基础上,采用同样策略,评估了多段翼型绕流计算中翼型后缘的几何处理、湍流模型选取等关键CFD要素对计算结果的影响。     

复合材料帽形加筋结构稳定性研究

建立了复合材料帽形加筋结构的参数化有限元模型,给出了两种进行层板结合处的模拟计算方法,通过计算得到了模型的屈曲载荷并总结了加筋结构几何参数对结构稳定性的影响,给出了模型的收敛性分析,对比了模型网格密度对于计算结果的影响,找到最佳的网格密度,提高了模型计算的效率和准确性,为实际的结构设计提供了重要参考,为结构分析提供了依据。     

气动弹性计算中网格变形方法研究进展

网格变形是气动弹性计算中实现计算网格随运动边界变形的主要方法。在总结目前网格变形方法发展现状的基础上,对近几年常用的网格变形方法,即弹簧法、弹性体法、超限插值法、Delaunay背景网格法、径向基函数插值法和温度体法等做了简要的总结。根据各方法构建模型的不同,将它们分成物理模型法、数学插值法和混合方法3类,简要介绍了各方法的基本思想和研究进展,重点比较了各方法的网格变形特性（变形能力、变形质量和变形效率）和优缺点。总结了当前流场边界与结构边界的数据传递方法,对当前气动弹性计算中遇到的网格变形的难点问题作了简要评述并对未来网格变形方法的发展方向进行了探讨。     

高精度差分格式及多尺度流场特性的数值模拟

简要介绍了作者近年发展的高精度差分格式,给出了两类格式精度的分析方法,分析了粘性面差分逼近所引入的耗散量低亏效应。文中按群速度将格式分为快型,慢型和混合型三类。对每类格式分析了高频波的传播特性。     

基于Optistruct的结构静动力拓扑优化设计

基于Optistruct平台,研究了二维平面薄板及三维结构的静力、动力以及静动力联合的拓扑优化,给出了典型结构的算例,对比了优化参数取值不同时所对应的优化结果,探讨了优化过程中的网格依赖性以及拓扑优化的效率等问题.     

用隐式多重网格法计算实用外形无粘绕流

网格展弦比的好坏严重地影响多重网格的效率。文中细致地处理了多重网格算法的细节,了网格展弦比的影响,发展了一种计算实用外形绕流的隐式多重网格解法,计算效率提高了５倍。     

逆流冷却罗茨鼓风机涡流与排气脉动的数值分析

运用计算流体动力学方法对逆流冷却在削弱罗茨鼓风机排气过程中存在的排气脉动方面进行了数值模拟,以动网格方法实现了计算域的变化并保证与转子旋转角度的匹配,在预进气与基元容积的动、静交接区域采用了非一致网格界面.对涡流与脉动流动瞬态特性的分析证明,逆流冷却可以避免排气口处的回流现象;在一个运行周期内,逆流冷却鼓风机的小脉动幅值仅为空冷鼓风机的20%～25%,这种小幅值脉动在降低排气气动噪声方面起重要作用.      

非结构混合网格鲁棒自适应技术

数值格式、湍流模型和计算网格是影响CFD数值模拟精度的3个主要因素。结合流场信息的网格自适应技术具备动态优化计算网格的能力,被NASA列为未来CFD发展的一项关键技术。本文针对非结构混合网格,发展了网格单元分布优化、表面网格几何投影和空间网格协调匹配3项关键技术,建立了高鲁棒性几何保真的网格自适应系统。首先,为了提高自适应方法的鲁棒性和通用性,发展了基于标准面网格的多面体网格单元分布优化方法。其次,发展了仅依赖表面网格信息的局部曲面重构技术,采用参数点映射方法实现了新增表面网格点的几何投影,消除了自适应系统对几何CAD系统的依赖。再次,采用改进的距离函数方法实现了空间网格与投影后表面网格的快速匹配。最后,结合基于流场特征的自适应探测器,采用二阶格式的有限体积方法,开展了30P30N三段翼绕流和三角翼大迎角绕流的网格自适应数值模拟。结果表明,通过网格自适应对网格单元的分布进行优化后,流场求解的收敛性和模拟精度都得到了显著提高。     

基于雷诺平均Navier-Stokes方程的表面传热系数计算

采用有限体积法数值求解控制二维绕流的雷诺平均Navier-Stokes(RANS)方程组,计算了光滑和粗糙NACA0012翼型以及圆柱表面的局部表面传热系数.分析了近壁面网格间距、湍流模式和表面粗糙度模型对数值计算结果的影响.结果表明:切应力输运(SST)湍流模型能够区分层流和湍流边界层的对流传热特性,并能预测转捩的发生;采用Spalart-Allmaras(S-A) 扩展模型能够计算粗糙壁面的对流传热系数,但采用忽略转捩函数的S-A模型不能有效计算层流边界层的传热系数.当近壁面网格间距接近10-5量级的黏性子层时,在光滑和粗糙壁面都能得到准确的传热系数分布.结合合适的近壁面网格间距,湍流模式和表面粗糙度模型可以得到与实验数据十分接近的表面传热系数曲线.通过与求解不可压缩RANS方程得到的结果比较后发现,不可压缩RANS方程主要忽略了压缩和黏性耗散效应,这种效应可以通过绝热升温项的形式并入总体热分析.      

非结构变形网格和离散几何守恒律

数值模拟流固耦合问题或多体分离问题的非定常流动时,常采用基于任意拉格朗日-欧拉（ALE）方程的有限体积法,涉及到变形网格和离散几何守恒律。在对变形网格算法进行综述时,按照构造思想分为物理比拟、椭圆光顺、插值、运动子网格（MSA）及其混合法共5类,分别介绍了基本原理、研究现状和适用范围,通过算例比较表明,径向基函数（RBFs）和运动子网格相结合的混合方法既有很好的变形能力,也有较高的计算效率,值得进一步发展和推广。在介绍了离散几何守恒律（DGCL）概念之后,采用二维几何模型进行分析,指出其机理是离散过程中体积增量与网格面元扫过的体积不相等造成的,把目前国内外应用的算法分为面积修正法、给定速度的面积修正法、速度修正法和体积修正法共4类,对其应用范围和存在的问题进行讨论,认为提出的体积修正算法既可以保证流固界面条件,也可以用于时间多层格式。     

非稳态挤压刚-粘塑性有限元分析中的误差估计和双重网格技术

 本文在刚-粘塑性有限元法中,应用误差估计和双重网格技术,并通过所设计的前后处理程序,有效和直观地模拟了管材非稳态挤压过程,用计算机绘出了分析网格和材料网格的变形图以及应变等值线图。     

关于 CFD 高精度保真的数值模拟研究

本文由四部分组成。第一部分研究了利用高阶格式、网格、物理模型和高性能计算机求解 NS 方程时,计算数据的高精度保真性的要求以及高精度保真的指标。第二部分研究了利用大规模计算开展大涡模拟和直接数值模拟存在的问题。第三部分分析了在高雷诺数下 NS 方程计算方法和网格的关系。最后提出了建立计算数值验证、确认的新方法。     

Delaunay graph-based moving mesh method with damping functions

The fluid–structure interaction and aerodynamic shape optimization usually involve the moving or deforming boundaries, thus the dynamic mesh techniques are the key techniques to cope with such deformation. A novel dynamic mesh method was developed based on the Delaunay graph in this paper. According to the Delaunay graph, the mesh points were divided into groups. In each group, a factor ranging from 0 to 1 was calculated based on the area/volume ratio. By introducing a proper function for this factor, this method can control the mesh quality with high efficiency. Several test cases were compared with other dynamic mesh methods regarding mesh quality and CPU time, such as radial basis function method and Delaunay graph mapping method.     

锥齿轮副误差测量成组夹具设计

本文重点讨论了锥齿轮副齿圈跳动误差、背锥及面锥相对零件中心线的位置误差、齿根余量误差和齿轮啮合中心距误差测量专用夹具解决的技术难点,提高锥齿轮副的误差检测的精度,为测量夹具提供思路。     

复合固体推进剂颗粒的三维网格自动生成

为了分析圆柱颗粒对复合固体推进剂性能的影响,必须对复合固体推进剂中三维圆柱颗粒进行网格划分。根据"力平衡"网格划分原理设计了一种适合于三维圆柱颗粒的网格自动剖分方法。编程运行结果表明,方法能自动完成划分,网格质量较高,满足要求。     

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

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

航空发动机主燃烧室火焰筒壁冷却的研究

概述了近年来国内外航空发动机燃烧室冷却方式的发展。重点论述了火焰筒壁温预测方法、三维壁温数值计算的控制方程、贴体正交网格和壁面函数的应用、计算辐射热流的蒙特卡罗法、三维计算域处理技术以及有高温涂层的壁温计算。调研表明,壁温预测技术对于未来高温升燃烧室设计而言,是有效、经济和实用的方法,对提高燃烧室寿命、飞行安全性具有重要意义。     

非结构混合网格自适应并行技术

计算流体力学（CFD）模拟实际工程问题所采用的网格规模可达千万量级，并行技术是减少计算时间的有效方法。耦合流场信息的网格自适应技术能有效动态优化计算网格，被NASA视为一项亟待发展的CFD关键技术。混合网格自适应系统包含网格分布优化、表面网格投影和空间网格匹配等关键技术。针对以上3项关键技术分别建立了高效的并行算法。首先，提出了"先唯一后同一"的两步法策略实现了网格单元分布优化过程的并行相容性；其次，基于局部曲面拟合思想，实现了曲面重构和新增物理网格点投影的完全并行；再次，提出了空间网格匹配技术的半并行算法，快速解决了网格单元交错问题。为了提高后续流场计算的并行效率，发展了基于并行重分区-网格数据迁移方法的动态负载平衡技术，并采用圆柱激波流场自适应模拟对动态负载平衡技术进行初步验证。最后，采用三角翼自适应加密测试了自适应系统的并行效率。结果表明，建立的混合网格自适应系统并行效率较高，且相比流场求解耗费总时间的比例低于1%。