适于黏性计算的自适应笛卡儿网格生成及其应用

复杂外形的黏性网格生成已经成为计算流体力学(CFD)在工程应用中的主要瓶颈,高效而鲁棒的网格生成技术研究显得尤为迫切。使用基于物面外形自适应加密的笛卡儿网格方法生成计算网格,投影方法拟合壁面边界,将投影得到的柱形单元沿法向分层获得适于黏性计算的网格。应用特征恢复技术恢复凹面特征并改善凸面处的网格质量。流场解算中应用最小二乘法重构获得2阶精度,Roe格式计算无黏通量,Holmes-Connell方法计算黏性通量。加入Spalart-Allmaras一方程湍流模型以获得湍流解。上下对称高斯 赛德尔迭代(LU-SGS)格式隐式时间推进,对网格重排序保证格式的稳定性和平衡性。算例显示网格生成方法是快速的和鲁棒的,流场计算精度满足工程需求,因此该方法非常适合应用于实际工程计算。     

适于黏性计算的自适应笛卡儿网格生成及其应用

复杂外形的黏性网格生成已经成为计算流体力学(CFD)在工程应用中的主要瓶颈,高效而鲁棒的网格生成技术研究显得尤为迫切.使用基于物面外形自适应加密的笛卡儿网格方法生成计算网格,投影方法拟合壁面边界,将投影得到的柱形单元沿法向分层获得适于黏性计算的网格.应用特征恢复技术恢复凹面特征并改善凸面处的网格质量.流场解算中应用最小二乘法重构获得2阶精度,Roe格式计算无黏通量,Holmes-Connell方法计算黏性通量.加入Spalart-Allmaras一方程湍流模型以获得湍流解.上下对称高斯-赛德尔迭代(LU-SGS)格式隐式时间推进,对网格重排序保证格式的稳定性和平衡性.算例显示网格生成方法是快速的和鲁棒的,流场计算精度满足工程需求,因此该方法非常适合应用于实际工程计算.     

基于非结构背景网格的动网格方法研究

借鉴弹簧动网格方法和Delaunay背景图方法,提出一种基于非结构背景网格的动网格生成方法。先用少量的物面控制点生成相对较粗的非结构背景网格,再用弹簧法将动边界的变化量作用到该背景网格,然后利用计算网格与背景网格的映射关系插值生成对应的动网格,最后对物面网格点及其附近受影响的密网格点进行局部网格修正,获得最终计算所需的动网格。所提方法与直接弹簧法比较,由于弹簧法作用的背景网格远粗于实际的计算网格,提高了计算效率,此外,由于网格的运动受控于受弹簧法作用的背景网格,边界变形或运动能较为均匀地传递到内部网格上,使网格单元过渡光滑,提高了动网格的质量。给出若干算例,展示了方法的效率和可行性。     

基于局部拟合的压气机叶片加工坐标生成方法

建立了基于前缘局部拟合技术生成叶片加工坐标的方法。包括网格插值、造型逼近、局部拟合和叶型光顺等步骤。以某多级高效率高负荷轴流压气机第2级转子为例,通过几何分析与计算流体动力学分析检验了该方法。结果验证了所生成的加工坐标的几何准确性,且对叶片热转冷有限元分析,本文方法所需的网格数量比现有技术节省约75%,计算时间节省75%以上。本文方法相比现有技术可提高工作效率,更好满足工程需求。     

非结构网格生成技术的研究

本文描述了非结构网格生成的过程。在生成方法上综合Ｄｅｌａｕｎｅｙ和Ａｄｖａｎｃｉｎｇ Ｆｒｏｎｔ两种方法之长，构造了新的网格生成方法，使非结构网格生成更方便、快捷。网格生成过程分三个部分：１．表面网格的构造和最佳逼近；２．辅助点的生成与确定；３．空间四面体网格的生成。使用上述方法，本文构造限复杂构形物体的非结构网格。     

网格生成软件GRGES的设计与实现

首先,本文说明了研制网格生成软件对于促进计算流体力学研究和应用的重要意义。其次,本文论述了数值网格生成的基本概念和基本步骤,提出了网格生成软件GRGES的设计思想,简单介绍了GRGES的总体结构和五种网格生成方法。最后给出了三个网格生成算例,它们均可用于计算实际流场。     

用代数法生成飞机翼型进气道等复杂形状物体的网格

本文论述了研制网格生成软件对于促进计算流体力学研究和应用的重要意义。着重介绍了网格生成的多面方法,编制了多面法网格生成程序,生成了飞机翼型、进气道等复杂形状物体的网格.该方法可快速生成二维或三维边界正交网格,同时还可有效地控制区域内网线的分布和边界上网格点的分布.     

旋翼桨尖涡生成及演化机理的高精度数值研究

为了细致捕捉直升机旋翼桨尖涡的生成和演化过程,建立了一个基于高精度网格和高阶通量计算格式的旋翼桨尖涡计算流体力学(CFD)求解方法。在该方法中,流场求解选取旋转坐标系下的Navier-Stokes方程为控制方程;空间离散采用迎风Roe格式,并采用低耗散的5阶WENO(Weighted Essentially Non-Osciltatory)格式进行对流通量的计算;时间推进则采用双时间法,在伪时间步上使用隐式LU-SGS(Lower-Upper Symmetric Gauss-Seidel)推进格式;应用嵌套网格方法实现桨叶网格和背景网格的数据交换。应用所建立的方法对悬停状态的旋翼桨尖涡流场进行了高精度模拟,在桨叶网格上精细地捕捉到了桨尖涡的具体生成过程,在背景网格上捕捉到了更多圈数的桨尖涡尾迹,并对桨尖涡的演化机理进行了深入研究。结果表明:建立的高精度数值方法能够有效地对旋翼桨尖涡的生成和演化过程进行细致模拟;悬停状态下旋翼桨尖气流在上下表面压力梯度的作用下经历了边界层增厚、逐渐卷起形成涡核以及最终脱离桨叶形成桨尖涡的过程。     

ALE方法中一种新的二阶保界守恒重映算法

在用拉格朗日格式求解流体力学问题时,随着时间的推进,计算网格会扭曲变形,影响格式精度,甚至导致计算中断。因此,要在网格变形较大时进行网格重分和物理量重映,以保证网格质量和格式精度。针对间断有限元方法求解流体力学问题的二阶拉格朗日格式,给出了一种守恒重映算法。该重映算法包括两步:第一步是用已有重映方法计算新网格上的单元平均值,并用相应修补算法对单元平均值进行调整,保证单元平均值的保界性;第二步是由已得到的新单元平均值重构出新网格上分片一次多项式,再使用Van Leer限制器对新网格上的梯度进行限制,使之不出现新的极值。最后用数值算例表明了该重映算法的保界性和二阶收敛性。     

飞机外挂物离机过渡过程的有限元算法.

研究飞机外挂物在按给定轨迹下沉旋转的离机过渡过程中，流场非结构网格的调整和再生成方法以及过渡过程欧拉方程有限元数值解法。中指出外挂物运动时流场网格调整或再生成仅需在围绕在外挂物的一小区内进行以减少计算量，而且固定网格的欧拉方程的有限元数值解法也可适用于运动网格情况，只要将固定网格的守恒通量项Fk用运动网格的Ek来代替，并建立了Fk和Ek之间的关系。末给出了带外挂物飞机的实例计算结果。     

一种适用于飞行器复杂外形的表面与空间网格生成方法

本文通过建立反映机身横截面外形特征的控制点概念，及在控制点间分配网格点数方式，给出了一种非常实用、交互式的表面网格生成方法，所生成的网格失真度小，物体保形好，生成适用于CFD计算需要的表面网络花费的时间少，适用性强。基于无限插值理论，引入本文提出的网格正交控制、物面法向量控制、加权平均光顺措施，给出了一种改进的空间网格生成方法，有效克服了传统无限插值网格生成方法在复杂外形网格生成方面的缺陷，有效改善了网格的法向疏密性、贴体性及周向网格的均匀性，对真实飞行器部件网格生成非常迅速。以某飞行器复杂机身表面与空间网格生成，及无人机翼身组合体空间网格生成为例，检验了本文网格生成方法的实用性和有效性。     

三维复杂域的Delaunay三角剖分

介绍一类自动生成内点的Delaunay三角剖分算法 ,其中包括重心插点法 ,外心插点法和Voronoi边插点法。该类方法适用于二维和三维复杂计算域的非结构网格生成 ,所需人工干预少 ,使用方便 ,计算效率高 ,生成的网格质量也很好。文中并给出几个算例说明本方法的应用。     

基于梯度下降的水滴收集率计算方法

拉格朗日方法在计算三维复杂外形的水滴撞击特性时，一般需要通过插值获得物面网格上的水滴收集率，而这一过程容易造成计算结果失真。建立了一种基于梯度下降的水滴收集率计算方法，通过确定撞击物面网格角点的水滴轨迹实现物面网格中心处的水滴收集率计算，避免插值运算。该方法的计算思路如下：根据撞击点到目标点的距离，采用梯度下降法自适应地调控水滴的初始释放位置，使撞击点落在壁面网格角点附近的给定阈值范围内；采用邻点检索算法，快速确定迎风面内所有撞击网格角点的水滴轨迹；基于流管定义，计算出网格中心处的水滴收集率。采用该方法对典型三维外形的水滴收集率进行了计算，并与相关文献结果进行对比。结果显示计算结果与文献数据吻合良好，有效克服了现有技术中因插值导致的结果失真问题，能够为飞机结冰及防除冰技术研究提供参考。     

旋翼流场计算嵌套网格并行装配方法改进研究

针对嵌套网格系统下进行旋翼流场计算时的网格装配问题,给出了一种改进的并行化网格装配模型。在该模型中,使用压缩存储的辅助结构网格进行并行"洞切割",生成插值边界;应用基于插值点坐标范围的局部可自适应辅助结构网格定位,缩小了处理器之间进行并行贡献单元搜索时的交互数据范围及生成辅助结构网格的内存消耗,结合局部范围的Neighbor to Neighbor(N2N)搜索重启和穷举方法,消除了网格分区、物面、外场等边界面对贡献单元搜索过程的影响,且方法对网格的拓扑没有限制,适合于结构、非结构、直角网格等网格块的装配过程。在流场计算方面,采用并行化的隐式双时间格式数值求解Navier-Stokes方程,用于旋翼数值模拟问题。理论分析和计算测试表明,在使用粗化的辅助网格时,改进的计算方法能够达到传统方法贡献单元搜索速度近似10倍,在低内存开销情况下降低网格装配过程对流场计算时间占比。在内迭代大于7次时,网格装配时间小于流场计算的1%。该方法是对旋翼多周期计算时反复网格装配过程的很好的效率改进。     

面向中小企业的动态联盟及其运行模式研究

将网格化方法应用到飞机S型进气道模型的铺放轨迹规划中,即利用UG网格生成器网格化进气道CAD模型,并生成网格化数据,以该数据为基础建立网格单元和节点之间的拓扑关系,在每个网格单元中计算铺放轨迹和轨迹点的法向量.通过仿真试验证明该方法具有程序执行速度快、生成的铺放轨迹曲线光滑、无曲线突变情况等特点,具有一定的工程应用价值.提出利用UG的二次开发功能开发适用于该轨迹生成方法的网格化工具.     

无人飞行器外形布局设计及其气动特性计算分析

设计了一种新的无人飞行器气动外形布局,应用基于二次曲线的模线设计方法完成了飞行器气动外形建模,并对所建模型进行气动特性计算分析。将飞行器几何外形节点数据导入Gambit软件,生成气动特性计算所需的网格文件。应用Fluent软件,根据不同的飞行条件,选择合适的计算流体力学分析方法,对无人飞行器外形布局进行气动特性计算。主...     

基于两步优化策略的结构网格负载平衡

国家数值风洞（NNW）工程旨在发展完全自主知识产权的计算流体力学（CFD）软件,结构网格负载平衡问题研究是该工程中的一个重要组成部分。本文发展了两步优化策略以求解结构化网格的负载平衡问题。第1步优化采用传统的贪婪算法,完成对大块网格的剖分和以进程计算时间为指标的网格块分配;第2步采用遗传算法（GA）,目标函数兼顾进程计算时间和通信时间,在第1步优化结果的基础上,对网格块在进程上的分配开展二次优化。为准确计算GA的目标函数,构建了一套计算时间和通信时间的建模方法,包括样本生成、模型建立和模型验证,整体方法具有一定的通用性。根据负载平衡问题以及两步优化策略的特点,对GA的编码、交叉、变异和种群初始化进行了研究,详细分析了交叉操作的递归问题及解决方法。算例验证说明建立的进程计算时间和通信时间模型具有较高的计算精度,能够用于GA的目标函数计算;两步优化策略能够在第1步优化的基础上进一步改善优化结果,从而减少CFD问题的整体计算时间,对于计算量巨大的工程问题具有较大的实用价值。     

绕HERMES外形高超声速无粘流动的数值模拟

本文采用时间相关法和隐式推进迭代方法求解了绕ＨＥＲＭＥＳ外形高超声速无粘流动。将抛物型方程和代数方程相结合生成计算网格，既节省了网格生成的时间，又能得到较合理的网格分布。在时间相关法计算中，采用了显式ＮＮＤ格式；在隐式推进迭代计算中，采用了推广的隐式ＮＮＤ格式。     

网格之星:国家数值风洞的通用型结构网格生成软件

结构网格具有边界层模拟准确、计算效率高等优势,在当前的计算流体力学(CFD)研究和工程应用中被大量采用。然而,复杂外形的结构网格生成非常耗时,通常需要在图形化软件上借助频繁的交互操作来实现,对CFD研究和应用产生了重要影响。针对此问题,提出了附面层推进、网格块自动补齐、网格分布快速调整、破损数模网格优化等网格生成加速技术。依托国家数值风洞工程,研发出一款通用型结构网格生成软件——网格之星(NNWGridStar)。与主流商业软件Pointwise相比,网格之星的网格生成速度更快。经过大量复杂算例验证,网格之星满足工程应用要求。自软件发布以来,已有80余家科研院所下载使用,为实际工程任务提供了重要的支撑作用。     

经典跨声速翼型RAE2822数据分析

经典跨声速翼型RAE2822风洞试验数据长久以来被广泛用于CFD计算方法和软件的验证与确认，但是数据的正确使用或者说合理使用仍存在一些需要研究和注意的问题，包括计算网格、风洞试验数据修正、中弧线修正、翼型几何定义和建模，以及摩擦阻力系数和边界层速度剖面的原始定义等。在开展CFD研究之前，必须首先对计算方法进行验证，尤其是要先尽可能消除计算结果对计算网格的依赖性；经过对目前可开放使用的计算网格的不足之处进行分析，研制了一套高品质的计算网格并获得了预期的一阶网格收敛性；通过计算软件交叉验证，进一步确保所用计算软件的可信度。在将CFD计算结果与翼型风洞试验数据进行比对时，通常需要对马赫数和攻角进行修正，且如何修正是一个需要持续研究的问题；翼型中弧线修正是一种有效的方法，但需要考虑流动参数的影响。原始翼型几何构型采用有限离散点定义，计算网格生成过程中需要采用插值方法布置型面网格点，不同插值方法对于翼型前缘附近流动的数值模拟有细微影响。多数相关研究工作只比对压力分布；少数研究工作会比对摩擦阻力系数、边界层及尾迹速度剖面，但在比对时，需要注意原始风洞试验相关参数定义与CFD计算常用定义的区别。