一种基于双重空间网格结合的高效DSMC实现方法

基于双重空间网格结合策略，发展了一种高效率的贴体DSMC（direct simulation Monte Carlo）方法。通过将仿真分子在物理空间结构网格的位置坐标映射于计算空间的直角网格中，并在计算空间中完成分子所属网格单元的定位以及分子与边界是否发生作用的判断，从而结合结构网格的贴体性和直角网格的高效率计算的优点，提高DSMC方法的贴体性和计算效率。基于双重空间网格结合策略，通过直接映射和间接映射法分别建立DSMC程序，对微尺度收缩扩张喷管气体流动和超声速圆柱绕流进行模拟。数值结果表明：两种方法均很好地模拟出微喷管因尺度缩小导致的黏性效应和速度滑移现象以及超声速圆柱绕流发生的激波现象，具备有效性。与传统结构网格方法对比，两种方法的计算效率平均分别提高了3.85倍和2.85倍，具备高效性。     

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

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

一种适合于旋翼涡流场计算的非结构自适应嵌套网格方法

针对旋翼涡流场的特点,提出了一种新的非结构自适应运动嵌套网格方法。在该运动嵌套网格方法中,贴体网格采用非结构网格,背景网格使用可自适应的笛卡尔网格,以便发挥非结构网格的形状描述能力和笛卡尔网格的空间自适应优势,方便网格生成过程并提高计算精度,同时克服在物体表面附近自适应笛卡尔网格生成工作量大和过程烦琐的不足。在此基础上,使用求解N-S主控方程的方法,分别对旋翼的悬停和前飞流场进行了数值模拟,计算结果表明本文的非结构自适应运动嵌套网格方法对提高旋翼涡流场的细节捕捉能力和计算精度是很有效的。     

基于共模型技术的大展弦比飞机气动/结构耦合性设计

在考虑弹性变形情况下,以大展弦比飞翼飞机为例,对气动/结构耦合设计方法进行研究提出了几何/结构/气动的共模型技术.首先是基于共模型技术生成表面气动网格和空间气动网格以及结构有限元网格,并实现气动载荷在结构模型和气动模型间的数据传递,最后进行气动/结构耦合计算分析,在进行气动/结构耦合计算分析时要进行多次迭代,直到结构变...     

TRIP2.0_SOLVER的开发与应用

本文介绍了“亚跨超CFD软件平台”2.0版本流场解算器部分(TRIP2.0_SOLVER)的开发工作和应用情况,重点是结构网格部分的开发工作,包括对接网格、拼接网格、重叠网格等多种网格拓扑结构。TRIP2.0_SOLVER采用有限体积法离散雷诺平均的N-S方程,数值模拟绕流飞行器的空间流场和飞行器的气动特性,本文从软件已经具备的基本功能,用户界面的开发、软件测试和和多种网格拓扑结构的应用等四个方面总结了亚跨超CFD软件平台流场解算器(TRIP2.0_SOLVER)的工作进展情况。     

基于多重网格方法的跨声速颤振数值模拟研究

使用多重网格的方法对二维翼型和三维机翼的跨声速颤振进行了数值求解.流场控制方程为N-S方程,湍流模型采用SA模型网格.计算网格采用结构网格,空间离散采用中心格式,使用双时间法进行时间推进.流场与结构之间的数据通过径向基函数(RBF)插值方法进行交换.通过耦合求解流场控制方程和结构运动方程得到Isogai(case A)机翼模型和Agard445.6机翼的跨声速颤振边界.通过与使用单重网格时得到的结果相比较得出,多重网格方法能够节省计算时间,提高颤振分析计算效率.     

基于网格变形技术的涡轮叶片变形传递

为实现学科耦合,需要在涡轮叶片多学科设计优化的气动模型和结构模型之间传递温度、气压和变形等载荷.研究结构变形向气动网格传递的实现方法.参数空间插值方法解决耦合面网格不一致时的变形插值问题;网格变形技术用来调整气动网格内部节点的位置,保证变形传递后的网格质量.某涡轮叶片变形传递实例的表明:该方法能够用于多学科设计优化问题中的变形传递.      

空间曲面的非结构网格生成

介绍一种空间曲面的 Delaunay三角形网格生成技术。采用参数表示法将空间曲面变换为平面参数域,引入相关的度量矩阵后将平面参数域进行三角剖分,最后得到空间曲面的三角剖分。空间曲面的边界点作为已知点被给定,内点采用一类自动生成的方法计算,例如按边插点法,Voronoi边插点法等。该方法可以很好地控制曲面网格的质量,使用方便,可以直接用于物理问题数值模拟中的曲面网格生成,也可以用于三维非结构网格生成中提供边界面三角剖分。文中给出的算例说明了本方法的应用情况。     

基于自适应混合网格的脱体涡模拟

基于混合网格和CGNS（CFD General Notation System）数据结构,建立了一种各向同性加密/稀疏的网格自适应方法。在悬空点的后处理中,让含有悬空点的单元转化为任意多面体,从而简化了自适应单元剖分模版,同时自适应网格单元之间可完全相容,自适应生成的网格能够直接用于可处理任意多面体的流场求解器。将该自适应方法与脱体涡模拟（DES）算法相结合,开展了65°后掠三角翼大迎角流动的数值模拟应用,并与初始网格的模拟结果进行了详细比较。对比表明：采用网格自适应方法适当增加局部网格量,能够以较小的成本迅速提高三角翼背风区的空间分辨率,增强数值模拟对小尺度涡系结构的解析能力,从而弥补了基于混合网格的脱体涡模拟中常用二阶格式计算的空间分辨率相对偏低、不利于湍流多尺度结构精细模拟的不足。     

复杂外形飞行器多块结构贴体网格生成

通过求解椭圆型方程生成三维贴体坐标网格,采用分块技术生成多块结构网格体系。每个子区域转换成计算空间内的立体域,具有自身的曲线坐标系,从而所形成的整个网格系统不会过分扭曲和不规则,这将有利于对于复杂边界的真实飞行器绕流流场的求解。文中讨论了控制函数对网格生成的影响。     

鲁棒的结构网格自动化重叠方法

网格技术是目前数值模拟中的关键技术之一。重叠网格是一种放宽拓扑要求、减小网格生成难度的网格技术。本文以结构重叠网格为基础,分别针对挖洞、寻点以及洞面优化方法进行了研究和改进,同时完成物面网格重叠,形成了一套鲁棒的、自动化的网格重叠系统。在挖洞方面,结合“最小洞映射”方法,提出“复合式挖洞”方法,节省内存开销;在寻点方面,通过构建格心虚网格,保证搜索空间的连续性,同时结合“有效搜索”思想,排除部分对寻点无贡献的网格点,进而减少ADT叉树节点;在洞面优化上,改变填补判别法则并引入两类受保护洞内点,确保两层插值边界建立,提高鲁棒性;在物面网格重叠上,利用物面投影法完成坐标修正,实现物面附近网格流动变量的准确传递。为验证本文方法,分别对定常翼身组合体DLR-F6绕流和非定常机翼挂载分离过程进行了数值模拟,计算结果与实验结果吻合良好,表明该结构重叠网格系统对多物体间定常、非定常扰流具有较好的数值模拟能力和较高的模拟精度,具有较高的工程应用价值。     

某栅格翼导弹减阻外形设计及其阻力特性数值研究

介绍和分析了国内外栅格翼的研究现状和技术难点。为了详细了解栅格翼气动布局栅格内部的流动特性,以及其与阻力设计满足飞行器性能要求的栅格翼,本文开展了几种不同单栅格形状外形的气动特性数值计算。针对传统数值计算中遇到的栅格翼计算网格生成困难、工作量大的问题,利用非结构空间直角网格对于复杂外形易于生成的特点对单个栅格的外形进行了初步选型研究。对于常规导弹和栅格翼导弹的气动性能进行了对比。     

先进复合材料格栅结构(AGS)应用与研究进展

 具有空间点阵的先进复合材料格栅结构（AGS）是一种新型高效的结构形式，在航空、航天结构中具有良好的应用前景。介绍了AGS优良的各项性能以及在美、俄等国航空、航天工程中的应用现状，着重描述了混合工艺法、模具膨胀工艺和拉挤 互锁格栅结构几种现行工艺方法，并对当前国内外针对格栅结构的力学性能和结构健康监测的研究进行了概述。最后提出了AGS进一步研究的工作展望。     

Numerical simulation of aerodynamic interaction for a tilt rotor aircraft in helicopter mode

A rotor CFD solver is developed for simulating the aerodynamic interaction phenomenon among rotor, wing and fuselage of a tilt rotor aircraft in its helicopter mode. The unsteady Navier–Stokes equations are discretized in inertial frame and embedded grid system is adopted for describing the relative motion among blades and nacelle/wing/fuselage. A combination of multi-layer embedded grid and ‘‘extended hole fringe" technique is complemented in original grid system to tackle grid assembly difficulties arising from the narrow space among different aerodynamic components, and to improve the interpolation precision by decreasing the cell volume discrepancy among different grid blocks. An overall donor cell searching and automatic hole cutting technique is used for grid assembly, and the solution processes are speeded up by introduction of Open MP parallel method. Based on this solver, flow fields and aerodynamics of a tilt rotor aircraft in hover are simulated with several rotor collective angles, and the corresponding states of an isolated rotor and rotor/wing/fuselage model are also computed to obtain reference solution.Aerodynamic interference influences among the rotor and wing/fuselage/nacelle are analyzed,and some meaningful conclusions are drawn.     

A compressing and saving method for structured grid data based on block construction

The increasing grid data in CFD simulation has brought some new difficulties and challenges, such as high storage cost, low transmission efficiency. In order to overcome these problems,a novel method for compressing and saving the structured grid are proposed. In the present method,the geometric coordinates of the six logical domains of one grid block is saved instead of all grid vertex coordinates to reduce the size of the structured grid file when the grid is compressed. And all grid vertex co...     

计算网格对气动力气动热数值模拟影响的研究

以二维圆柱钝头超声速流场为研究对象,应用计算流体力学方法研究了计算网格在气动力及气动热数值模拟中的影响.计算程序中对流项、扩散项分别采用AUSM格式和中心差分格式离散,并用多步Runge-Kutta显式时间推进法求解空间离散后的控制方程.计算结果表明,壁面法向第一层网格间距、网格长宽比以及扩张度等因素都会对壁面摩阻、热流的计算产生影响,尤其是壁面法向第一层网格间距影响最为突出.边界层的计算结果与实验数据进行了仔细对比,包括壁面压力分布,摩阻系数Cf分布,斯坦顿数St分布等.高马赫数滞止气流导致壁面附近产生巨大速度梯度以及温度梯度,这些梯度的准确计算需要合适的壁面法向第一层网格间距以及合理的网格长宽比、扩张度.     

基于非结构网格有限差分法的扎染算法

应用有限差分法时遵循贴体坐标系下离散等价方程的离散准则,格式中仅用到当地网格点的坐标变换系数,据此可以构建新的非结构网格有限差分法,数值算例表明在空间二维情况下用连接离散点的3条网格线构造出来的一阶迎风格式可以稳定收敛。把这种非结构网格有限差分法推广到三维,在4条网格线基础上进行离散计算,然后提出一种局部区域多重网格、重复计算的新算法。首先在包含物体的整个区域采用直角坐标系下的均匀网格进行计算,然后消除物体内部点对物体外部流场区域的影响,最后在物面和直角网格之间很小的局部区域填补非结构网格进行计算,计算过程类似于中国传统的扎染工艺。处理流场内存在曲面物体边界问题时,相较于常规结构网格差分算法,扎染算法尽管在计算中包含需要特殊处理的无用网格点,也存在重复计算的局部区域,但是以这些计算能力的浪费为代价,换取了编程简便、内存量小、网格生成快捷、易于扩展网格规模等优势。定性分析表明,这种扎染算法非常适合研制大规模并行计算,建立的计算流体力学应用软件可以充分发挥数十万至百万处理器核心的超级计算机的效能。     

二维及三维贴体网格的一种生成方法

本文讨论了二维及三维情况下应用Ｌａｐｌａｃｅ方程生成贴体网格的方法，并应用ＡＦ１格式求解方程以大大加速迭代收敛的速度，对于二维及三维情况，分别讨论了相应用Ｄｉｒｉｃｈｌｅｔ边界问题解法。对于二维情况，还讨论了翼型的Ｏ－、Ｃ－和Ｈ－型三种典型的网格生成方法。着重讨论了有关割缝上的网格点的处理方法，以保证网格的光滑性。三维情况下，以简化航天飞机外形为例，讨论了网格的生成方法，尤其着重讨论了边界点的处理     

Z源NPC三电平变换器恒功率并网控制策略研究

Z源中点钳位式(NPC)三电平变换器把逆变和升压两种环节结合在一起,结构简单、工作效率高,但是针对这种结构输出的无功功率研究较少。提出了一种Z源NPC三电平变换器的恒功率并网控制策略。首先对这种变换器的结构及原理进行分析；再把恒功率控制引入到这种变换器的并网控制中,并采用空间矢量脉宽调制(SVPWM)实现了变换器的并网控制。MATLAB/Simulink软件仿真结果证明了这种变换器能够提高电压升压比。该控制策略可使变换器输出的有功功率、无功功率稳定于设定值,且并网电流谐波较低。