三维复杂外形的非结构网格自动生成技术与应用

给出了一套三菜的非结构网络自动生成方法。任意外形由几种代表性的曲线和曲面来描述,将曲面变换为平面,在平面内由二维阵面推进法生成三角形网格,然后反变换为物形表面网格。空间的四面体网格由三维阵面推进法生成。四面体网格的优化采用了附加“关联质量”约束的节点松驰的局部的Ｄｅｌａｕｎａｙ变换技术。另外还发展了非结构风格的自适应技术。数值算例表明了方法的灵活性和通用性     

NND格式在非结构网格中的推广

在张涵信提出的无波动、无自由参数的差分格式（ＮＮＤ格式）的基础上，构造了适用于非结构网格的二阶精度ＮＮＤ有限体积格式，解决了现有非结构网格方法中为抑制激波附近的波动而必须引入含自由参数的人工粘性项的困难，并采用网格自适应技术以提高效率．通过对二维平板激波反射和前台阶在管道内的流动问题的计算，表明本方法可有效地用于Ｅｕｌｅｒ方程的求解．     

虚拟飞行中气动、运动和控制耦合的数值模拟技术

随着计算机科学和计算流体力学的发展, 以非定常数值模拟为核心, 开展气动/运动/控制等多学科耦合的“数值虚拟飞行” 模拟成为可能. 数值虚拟飞行有助于飞行器设计师在设计之初和整个设计过程中分析和评估飞行器的非线性飞行力学和稳定性与控制性能. 该文综述了国内外数值虚拟飞行中“气动/运动/控制” 耦合的一体化模拟技术的研究进展, 分析了其中的关键科学和技术问题, 重点介绍了气动/运动/控制耦合一体化计算方法, 并介绍了作者在一体化耦合计算方法方面取得的进展及初步应用情况. 最后探讨了数值虚拟飞行中的一些挑战性问题, 并对未来发展趋势进行了展望. 可以预期, 随着E 级计算的到来, 在不久的将来, 数值虚拟飞行将给新型飞行器设计带来革命性的变化.      

动网格生成技术及非定常计算方法进展综述

对应用于飞行器非定常运动的数值计算方法(包括动态网格技术和相应的数值离散格式)进行了综述.根据网格拓扑结构的不同,重点论述了基于结构网格的非定常计算方法和基于非结构/混合网格的非定常计算方法,比较了各种方法的优缺点.在基于结构网格的非定常计算方法中,重点介绍了刚性运动网格技术、超限插值动态网格技术、重叠动网格技术、滑移动网格技术等动态结构网格生成方法,同时介绍了惯性系和非惯性系下的控制方程,讨论了非定常时间离散方法、动网格计算的几何守恒律等问题.在基于非结构/混合网格的非定常计算方法中,重点介绍了重叠非结构动网格技术、重构非结构动网格技术、变形非结构动网格技术以及变形/重构耦合动态混合网格技术等方法,以及相应的计算格式,包括非定常时间离散、几何守恒律计算方法、可压缩和不可压缩非定常流动的计算方法、各种加速收敛技术等.在介绍国内外进展的同时,介绍了作者在动态混合网格生成技术和相应的非定常方法方面的研究与应用工作.     

计算流体力学中的验证与确认

计算流体力学(CFD)在航空航天等诸多领域的应用越来越广泛.特别是近年来, CFD在实际飞行器的设计中扮演着越来越重要的角色, 许多设计参数直接来源于CFD的计算结果.由此, 飞行器设计师对CFD提供结果的可信度提出了更高的要求.验证(verification)与确认(validation)是评价数值解精度和可信度的主要手段.本文综述了国内外开展CFD验证与确认研究的进展.在引言中论述了开展CFD验证与确认的重要性和必要性, 简述了国内外CFD验证与确认研究的历史和发展现状.第2节中讨论了CFD验证与确认的一些基本概念,以及这些概念定义的形成过程, 并指出了进行CFD验证与确认的基本步骤.第3节和第4节分别讨论了CFD验证与确认的方法, 如CFD验证中的精确解比较方法, 制造解比较方法, 网格收敛性研究;CFD确认中的层次结构, 流动分类法, 确认实验指南.在第5节中我们列举了几个CFD验证与确认的应用实例.最后, 对我国开展CFD验证与确认研究工作提出了若干建议, 包括:(1)开展流动分类法研究,(2)推行软件质量工程方法,(3)开展规范精细的实验, 建立国内的网络数据库.      

动态混合网格生成及隐式非定常计算方法

建立了一种基于动态混合网格的非定常数值计算方法. 混合网格由贴体的四边形网格、外场 的多层次矩形网格和中间的三角形网格构成. 当物体运动时，贴体四边形网格随物体运动而 运动，而外场的矩形网格保持静止，中间的三角形网格随之变形；当物体运动位移较大，导 致三角形网格的质量降低，甚至导致网格相交时，在局部重新生成网格. 新网格上的物理量 由旧网格上的物理量插值而得. 为了提高计算效率，采用了双时间步和子迭代相结合的隐式 有限体积格式计算非定常Navier-Stokes方程. 子迭代采用高效的块LU-SGS方法. 利用该 方法数值模拟了NACA0012振荡翼型的无黏和黏性绕流，得到了与实验和他人计算相当一致 的结果.     

基于非结构/混合网格模拟黏性流的高阶精度DDG/FV混合方法

间断Galerkin有限元方法(discontinuous Galerkin method, DGM) 因具有计算精度高、模板紧致、易于并行等优点, 近年来已成为非结构/混合网格上广泛研究的高阶精度数值方法. 但其计算量和内存需求量巨大, 特别是对于网格规模达到百万甚至数千万的大型三维实际复杂外形问题, 其计算量和存储量对计算资源的消耗是难以承受的. 基于“混合重构”的DG/FV 格式可以有效降低DGM 的计算量和存储量. 本文将DDG 黏性项离散方法推广应用于DG/FV 混合算法, 得到新的DDG/FV混合格式, 以进一步提高DG/FV混合算法对于黏性流动模拟的计算效率. 通过Couette流动、层流平板边界层、定常圆柱绕流, 非定常圆柱绕流和NACA0012 翼型绕流等二维黏性流算例, 优化了DDG 通量公式中的参数选择, 验证了DDG/FV 混合格式对定常和非定常黏性流模拟的精度和计算效率, 并与广泛使用的BR2-DG 格式的计算结果和效率进行对比研究. 一系列数值实验结果表明, 本文构造的DDG/FV混合格式在二维非结构/混合网格的Navier-Stokes 方程求解中, 在达到相同的数值精度阶的前提下, 相比BR2-DG格式, 对于隐式时间离散的定常问题计算效率提高了2 倍以上, 对于显式时间离散的非定常问题计算效率提高1.6 倍, 并且在一些算例中, 混合格式具有更优良的计算稳定性. DDG/FV 混合格式提升了计算效率和稳定性, 具有良好的应用前景.      

庆祝中国力学学会成立50周年大会暨中国力学学会学术大会'2007会议简介

1 会议概况 值此中国力学学会50华诞之际,为回顾中国力学的光荣历程、展现我国力学的优秀成果、弘扬中国力学的优秀传统、展望新世纪力学学科的发展趋向,中国力学学会成立50周年暨中国力学学会学术大会'2007(Chinese Conference of Theoretical and Applied Mechanics-2007,CCTAM'2007)于2007年8月20～22日在北京召开.大会由中国力学学会主办,39个单位协办.     

HyperFLOW软件非结构网格亚跨声速湍流模拟的验证与确认

计算流体力学(computational fluid dynamics, CFD)数值模拟在航空航天等领域发挥越来越重要的作用,然而CFD数值模拟结果的可信度仍然需要通过不断地验证与确认来提高.本文给出了从制造解精度测试、简单到复杂外形湍流模拟网格收敛性研究等三个方面开展CFD软件验证与确认的方法,并对自主研发的CFD软件平台HyperFLOW在非结构网格上模拟亚跨声速湍流问题的能力进行了验证与确认.首先通过基于Euler方程和标量扩散方程的制造解精度测试,分别验证了HyperFLOW在非结构网格上对Euler方程和黏性项的求解精度,结果表明其能够在任意非结构网格上达到设计的二阶精度.其次,通过NASA Turbulence Modeling Resource中的湍流平板、二维翼型近尾迹流动、二维Bump等几个典型的亚声速湍流算例的网格收敛性研究,量化考察了数值结果的观测精度阶和网格收敛性指数,并与国外知名CFD解算器CFL3D,FUN3D的计算结果进行了对比,验证了HyperFLOW对简单湍流问题的模拟能力,且具有良好的网格收敛性和计算精度(阶).最后,通过NASA Common Research Model标模定升力系数的网格收敛性研究和升阻极曲线预测,验证了软件在复杂外形亚跨声速湍流流动数值模拟中也具有良好的可信度.     

基于动态混合网格的不可压非定常流计算方法

鱼类、昆虫等运动速度较低,对它们的数值模拟需要解决不可压问题.虚拟压缩方法通过在连续性方程中加入压强对虚拟时间的偏导数,从而把压力场和速度场耦合起来,解决了不可压缩流的计算问题.基于动态混合网格技术,利用双时间步方法耦合虚拟压缩方法来解决非定常不可压缩流的计算问题.为了加快每一虚拟时间步内的收敛速度,子迭代采用了高效的块LU-SGS方法,并且耦合了基于混合网格的多重网格方法.利用该方法数值模拟了不同雷诺数下的静止圆柱、振荡圆柱的绕流,得到了与实验和他人计算一致的结果.