矩形／非结构混合网格技术及在二维／三维复杂无粘流场数值模拟 …

建立了一套模拟复杂无粘流场的矩形／非结构混合网格技术，其中非结构网格仅限于物面附近，发挥非结构网的几何灵活性，以少量的网格模拟复杂外形，同时在以外的区域采用矩形结构网格，发挥矩形网格计算简单快速的优势，有效地克服全非结构网格计算方法需要较大内存量和较长ＣＰＵ时间的不足。     

非结构化多重网格流场数值计算

本文对二维非结构网格无粘流场计算的多重网格方法作了探讨。研究了非结构化多重网格的生成方法,插值算子的构造以及流场计算．以绕流ＮＡＣＡ００１２翼型无粘流场为算例,研究了多重网格法在非结构网格上的实现方法及其加速效果表明本文的方法有显著的加速收敛效果。     

基于各向异性非结构网格的超声速流动自适应计算

基于各向异性非结构网格,实现了自适应求解技术,将其应用于超声速楔形体绕流及超声速横向喷流问题的计算研究,对各向异性网格自适应计算过程的收敛性和求解精度进行了对比分析,展示了各向异性网格自适应算法在降低问题求解规模。提高各向异性流场分辨率等方面的优势及存在的问题。超声速楔形体无粘绕流场的自适应计算,清晰地捕捉到了激波附近的流场信息,激波前后的马赫数、压力、密度。温度等物理量与理论分析解吻合很好。超声速横向喷流流场存在激波、分离涡、边界层等流场结构的相互干扰。计算研究表明,单纯基于Hes-sian度量张量的各向异性网格生成及自适应算法不能有效模拟边界层内的流动情况,是将来需要进一步开展的研究挑战。     

超声速湍流流场的RANS/LES混合计算方法研究

采用对接/拼接网格技术,建立了基于分区混合和基于湍流尺度混合的双重RANS/LES混合计算模型,并对环翼低速绕流、翼型跨声速绕流和球锥带凹窗外形二维超声速绕流进行了初步的数值模拟.环翼和翼型绕流计算表明,该混合模型可给出较合理的湍流宏观平均量;球锥带凹窗外形二维超声速绕流计算表明,该混合模型可得到超声速瞬态湍流脉动流场,凹窗处存在复杂的旋涡结构和波系结构,呈现较大尺度的脉动.但该模型还需要进一步的考核验证.     

Jameson有限体积法对非结构网格推广的改进

本文将Ｄｅｌａｕｎａｙ法用于复杂物形的非线性结构网格生成,并利用黄明格提出的布点原理避免边界穿透。绕流的计算采用Ｅｕｌｅｒ方程和Ｊａｍｅｓｏｎ的有限体积法。本文对三维问题提出新的激波敏感器,对方法做了改进,并用于ＡＧＡＲＤ Ｍ６机冀、导弹模型和飞机模型的绕流计算。     

高超声速连续/稀薄流自适应界面推进重叠网格方法及应用研究

根据高超声速飞行器在过渡流域飞行时绕流流场具有连续流与稀薄流并存的特点,提出了一种适用于混合流场自适应界面推进的多尺度重叠网格方法。用两套独立生成不同尺度的结构网格和非结构网格覆盖整个流场。根据连续/稀薄流界面失效函数的截止值在结构网格内进行自适应"挖洞"。在挖洞后的稀薄流区利用非结构网格DSMC方法进行计算;在连续流区的结构网格内采用基于N-S方程的CFD算法。两种算法进行耦合计算和求解,在连续/稀薄流分界面进行流场信息的插值、交换和处理。根据流场实时耦合计算结果,利用界面失效函数自适应改变两套网格的重叠部分。最后,通过对超声速圆柱和高超声速带扩张角圆管绕流的数值模拟,验证了方法的可行性和有效性。     

开放式气动力数值模拟系统研究

气动力数值模拟系统是ＣＦＤ流场解算技术、网格生成技术、数据可视化技术和网络技术相结合的产物。描述了一种运行于计算机网格环境下多任务、多用户、交互式的数值模拟软件环境－ 气动力数值模拟系统,该环境能计算模拟亚、跨、超音速绕流的气动力,也能方便地集成模拟不同复杂气动布局飞机流场的解算程序。介绍了其网格生成、数据可视化、系统接口等核心部分的结构和功能。     

基于非结构重叠网格的二维N—S方程求解与应用研究

本文发展了非结构重叠网格技术.通过构造任意单元类型网格重叠的网格间边界确定算法、宿主单元搜索算法以及网格间插值算法,把非结构重叠网格推广应用于求解NS方程,并结合湍流模型模拟二维湍流流场.N-S方程的求解采用有限体积法,通量计算采用Osher格式,湍流模型采用SA模型.为验证本文方法的正确性和有效性,模拟了RAE2822翼型和某三段翼型绕流湍流流场,结果表明采用非结构重叠网格模拟复杂外形湍流流场不但网格生成效率高,而且容易推广应用于具有相对运动的非定常问题,具有较高的工程应用价值.     

非结构网格及混合网格复杂无粘流场并行计算方法研究

在分布式存储网络异构编程PVM环境下,本文将原串行计算程序改造为并行计算程序,提高了计算效率.利用非结构网格,对球体及航天飞机外形绕流流场进行了计算,取得了较好的性能结果.最后将该方法应用于混合网格的计算,求解了多种复杂外形的无粘绕流,如翼身组合体和DDF1C导弹,得到了满意的结果.     

基于结构网格的栅格翼绕流数值模拟

针对栅格翼这种复杂的构造形式,采取有效的多块结构网格生成策略,对栅格翼导弹生成高质量贴体的结构网格。在此基础上,通过求解N-S方程,对栅格翼导弹亚、跨、超声速绕流流场进行了数值模拟。计算结果与风洞试验结果吻合较好,能够准确地捕捉到栅格框内复杂流场的流动特性,从而结构网格生成策略与计算方法得到验证,为栅格翼导弹绕流流场的精确数值模拟提供了技术途径。     

非结构网格生成方法的改进及气膜冷却三维静子流场的求解

讨论五面体网格下,涡轮静子叶片含多排冷却孔时三维粘性流场雷诺平均Navier-Stokes方程组的数值求解问题。考虑到叶片含多排冷却孔的工程背景,本文发展了生成五面体非结构网格的快速有效方法即仅在一个拟S₁面上生成非结构网格,而沿叶高方向仍采用结构网格的生成办法。文中还采用了新的网格优化措施以改善生成网格的质量。流场的计算采用了非结构网格下的高分辨率迎风格式。几个典型算例表明:本文发展的网格生成方法快速有效,而且网格生成质量较高;发展的计算格式能够获得较为满意的三维流场数值解,并具有较高的激波分辨率。      

任意平面域的三角形网格和混合网格生成

本文以阵面推进法为基础,给出了一种可对任意平面域进行三角形网格和混合网格自动生成的方法。该方法包含了易于生成与调控的矩形背景网格技术,简便的外形输入和初始阵面划分技术,加速查寻的数据结构,以及节点松弛和基于Ｄｅｌａｕｎａｙ准则的对角线交换技术。网格生成的实例表明,方法能快速可靠地生成任意平面域的高质量三角形网格和适于粘性计算的三角形／四边形混合网格,具有较好的通用性与实用性。     

复杂外形的动态混合网格生成方法

本文将层推进方法、四分树方法和阵面推进法相结合,建立了二维复杂外形的动态混合网格生成方法.首先用层推进方法在运动物体附近生成贴体的四边形网格,然后用四分树方法在外场生成自适应的矩形网格,而中间由三角形网格过渡.三角形网格由阵面推进法生成,并由节点松弛和"对角线交换"方法进行优化.当物体运动时,贴体的四边形网格随物体运动而运动,而外场的矩形网格保持静止.中间过渡层的三角形网格随物体运动而变形.当物体运动位移较大导致出现质量较差的单元甚至导致网格相交时,则在局部重新生成网格.利用上述动态混合网格生成方法,生成了多种含相对运动的复杂外形的动态混合网格,网格质量令人满意.     

用非结构直角网格和欧拉方程计算飞机绕流

 采用八叉树结构 ,生成复杂外形绕流计算的非结构直角网格。物面附近用投影方法 ,使网格贴体。并将Jameson的有限体积法推广用于这种网格的欧拉方程计算。对歼击机模型的绕流计算表明 ,网格生成的机时花费很少 ,总体质量好 ,因而欧拉方程解算的收敛质量也好。     

复杂积冰翼型气动性能分析

为了克服复杂外形难以划分结构化网格的问题, 采用结构/非结构混合多块网格对复杂积冰翼型的气动性能-进行了分析计算.在靠近积冰边界的内层采用了非结构网格, 在非结构网格外采用结构化网格以节省存储空间和计算时间.不同类型网格之间采用了点点对接的方式以简化通量重构的处理方式.整个流动区域采用SST湍流模型, 求解了雷诺平均N-S方程.数值计算结果表明结构/非结构混合多块网格在复杂积冰翼型气动性能计算上的可行性, 积冰对翼型的气动性能造成了严重的影响.      

混合网格方法在栅格翼数值模拟中的应用研究

栅格翼是采用一种蜂房结构设计的升力面和控制面,和传统舵面相比具有铰链力矩小、能在较大攻角下保持升力等特点,并得到了越来越广泛的应用。由于栅格尾翼的特殊结构,网格生成具有很大的难度。本文采用了结构、非结构混合网格的办法,在弹体和栅格翼的物面附近区域,生成适合粘性计算的大长宽比的结构网格,弹体网格和栅格翼网格之间采用非结构网格进行填充,满足非结构网格和结构网格交接面的完全对接。本文基于结构/非结构网格体系采用有限体积方法求解NS方程,对不同舵偏角下的栅格翼构形进行了数值模拟,并通过实验结果对数值方法进行了验证。     

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.     

非结构网格上新型的NND有限元格式

张涵信发展的ＮＮＤ差分格式是由中心差分格式,二阶迎风格式和一阶迎风格式混合组成的杂交型格式。众所周知,和中心差分格式相对应的是Ｇａｌｅｒｋｉｎ有限元格式。通过对中心型有限元格式加修正项的方法本文成功地构造出二阶迎风型有限元格式和一阶迎风型有限元格式。     

自适应笛卡尔网格超声速黏性流动数值模拟

复杂外形/流场的高质量网格的生成往往需要占用大量人力资源,而自适应笛卡尔网格方法能够自动化生成高质量网格,具有很好的工程实用价值和应用前景。基于笛卡尔网格方法,采用叉树数据结构进行数据的存储和访问,分别从几何特征和流场解特征出发进行网格的自适应加密和粗化,发展了一种二维情况下自动、高效的自适应笛卡尔网格生成方法。从浸入边界方法出发,结合虚拟镜像对称方法和曲率修正技术进行黏性物面边界条件的处理,同时建立了多值点问题的处理技术,发展了一种在笛卡尔网格下可有效模拟黏性物面边界条件的方法。针对自适应笛卡尔网格非均匀的特点,发展了悬挂网格的处理方法,并构建了适用于自适应笛卡尔网格的黏性数值求解器。通过典型算例的考核,验证了所发展的自适应笛卡尔网格生成技术和构建的数值求解器具有较高的精度和可靠性。     

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

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