强制偏流喷管的三维网络生成

运用泊松方程实强了制偏流喷管内的三维网络生成，方程中的源项采用了与边界形状和边界网格点分布无关的数学形式，源项内的自由参数由横向网格线与边界的交角及横向网格线在边界处的局部曲率两个约束条件来确定，实现了由边界上的网格分布直接控制内部的网格点，使生成的网格具有良好的正交性和合理性。     

圆柱管多分支联接三维分块贴体网格的耦合生成

利用分块耦合技术生成了两圆柱管成任意角度,管径不相等的三分支联接的三维贴体网格,所介绍的网格生成方法在分块生成网格时能保证穿过分界面上的网格线是光滑连续的。此外,用方法生成网格时间通过边界上网格的疏密或者稳定边界上网格线之间的夹角来控制内部区域的网格疏密。     

贴体与边界正交网格及大迎角涡流计算

本文应用空间理想非均匀线性弹性力学平衡关系生成三维贴体与边界正交网格。通过弹线性的刚度来控制网格的密度，并且自动调整边界节点上外力分布使网格线与边界正交。文中还用欧拉方程有限体积法三步显式时间推进及焓修正方法求解了绕三角翼的大迎角流动。计算结果表明这种方法能够较好模拟跨声速胶体涡流动。     

复杂空气动力外形网格生成技术

通过对用于网格生成的Ｌａｐｌａｃｅ方程进行抛物化差分离散,发展了一种快速抛物型方程贴体及与边界正交网格生成方法,另外,采用源项修正技术求解用于网格生成的Ｐｏｉｓｓｏｎ方程,发展了一种可实现对网格线与边界间距及正交性进行双重控制的椭圆型方程方法,运用这种两种基本方法以及这两种方法的混合,分别生成了真实歼击机全机外形,带十字叉翼导弹以及双发,下单翼运输机外形的计算网格。并已成功地秀于跨声速得杂流场的娄     

机翼机身跨声速绕流的计算

本文采用变换与数值网格生成方法，对由椭圆抛物而描述的每一个展向曲面进行网格生成，引入一种用于控制网格疏密分布的函数变换方法，其特点是易于对不同族网格线分别控制，根据物面梯度自适应地调整网格的疏密。用有限体积法计算守恒型全速势方程跨声速绕流，对M-100翼身组合体进行计算。理论计算结果与实验结果表明，本文方法有效。     

两种椭圆型方程求源项方法在喷管内流场网格生成中的应用

用椭圆型方程生成二维网格，发展两种Hilgenstock求源项方法，实现对网格与边界间距及与边界正交的控制。     

用欧拉方程计算机翼与翼身组合体绕流

本文用保角转绘加剪切变换生成三维机翼与翼身组合体Ｃ－Ｈ型网格，并在这类网格拓扑上，用有限体积法研制出可供分析机翼与翼身组合体绕流的三维欧拉方程计算程序。本方法的特色是改进了机翼表面网格点的分布，使机翼后缘与网格线一致；采用了当量机身，使机岙表面到对称面光滑过渡；嵌入了边界层粘性修正。改善了计算结果。数值计算表明，对粘性影响较大的超临界绕流，边界层修正效果显著。     

双边界代数法在喷管内流场网格生成中的应用

用一种双边界代数方法生成二维网格，通过在超限插值函数中引入边界的一阶导数实现网格与边界是的正交，采用几何一致性条件保证物理区域不变形，发展一系列拉伸方程控制网格在边界或给定位置的间距。应用此方法生成了满意的喷管内流场网格，网格生成效率明显优于微分方程网格生成方法。     

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

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

基于超网格的重叠网格守恒插值方法

保证重叠网格边界数据插值的守恒性是计算流体力学面临的一大难题。基于超网格方法与格心格式有限体积法,发展了一种新型混合重叠网格守恒插值方法。详细研究了用网格切割技术在重叠网格边界构造局部超网格,用网格求交算法合理地扩大贡献单元模板,建立了一种适用于任意网格类型的隐式并行重叠网格守恒插值方法,以超网格为媒介可实现重叠网格二阶精度的守恒插值。数值结果表明,本文方法对二阶分布的流场变量具有严格的守恒性,相比三线性插值方法和逆向距离权插值方法,本文方法减小了数值误差,提高了重叠网格边界的插值精度,加快了计算收敛速度,改善了重叠区域网格尺度相差较大时流场的光滑性和连续性。     

翼-身组合体绕流的Euler方程数值模拟

用代数方法生成翼-身组合体H-O型网格,并用有限体积法研制出翼-身组合体绕流三维Euler方程计算程序。该方法的特点是改进了机身与机翼表面网格点分布,机翼后缘有后掠时也能保证后缘与网格线一致。程序除能按常规提供横流截面展向压强分布外,还能提供弦向压强分布。对NASA TND-712的翼-身组合体模型的计算结果与实验符合很好。     

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

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

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

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

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

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

基于背景网格簇的动网格生成方法

提出了基于背景网格簇实现网格变形的新方法。背景网格簇由任一内边界节点和远场边界角点生成的背景网格构成,并随边界运动而变化。通过保持计算网格节点在每一个背景网格中所在单元的面积坐标不变(三维时保持体积坐标不变)而求出一组期望位置,加权平均该组期望位置,确定出该计算网格节点的新位置。其中权值与内边界节点和计算网格节点距离的倒数相关。新方法引入的背景网格簇可以有效地改进单个背景网格的不足之处。算例结果表明:基于背景网格簇的动网格生成方法实现简单,与弹簧近似法相比,新方法因不需要迭代求解方程组而非常高效,且拥有更强的网格变形能力;与Delaunay图映射法相比,该方法背景网格的单元个数极少,因此易于定位,且不会出现背景网格单元交叉的现象,网格变形能力更强,变形后的网格质量更好。     

实际机身跨音速流场分析计算

本文应用修正的跨音速小扰动(MSD)方程和疏密网格交替迭代法,将Jameson格式在非等步长网格中的推广形式与Engquist-Oslier格式相结合,在机身侧边局部横向网格线上使用横向线松弛法,经数值实验得到了一个收敛速度比较快的计算方案。许多算例均证明本方案迭代次数少,计算结果与实验结果符合得也很好。     

任意旋成面叶栅的跨声速绕流计算

本文中的方法是把时间推进有限体积法和流线迭代法结合起来,用以解决任意旋成面叶栅跨声速绕流的计算问题。在旋成流面上,网格是由一族流线及一族固定的节距线所组成的。流线的位置预先并不知道,在计算过程中需要逐步加以修正。采用流线作为一族网格线可使网格线的分布更为合理,并使差分格式简化。同时,为了加速计算过程的收敛,在本计算中采用逐步加密网格技术。     

粘性流体对流换热问题用直角坐标网格的解法

研究了非结构化直角坐标网格全贴体求解粘性不可压流体自然对流换热问题。该方法在于利用直角坐标网格但通过在边界附近保留不规则控制体,使得算法是完全贴体的。利用此数值方法对二维倾斜空腔内的自然对流及上固体边界移动时的对流换热问题进行计算,并与STAR-CD的计算结果进行比较。对比结果表明,该方法计算的结果与STAR-CD的结果很接近,可以满足一定精度要求。      

基于四叉树网格下的浸入边界-格子Boltzmann方法

提出了一种基于非均匀四叉树网格下的浸入边界-格子Boltzmann方法.在不同层网格界面上,为了保证物理量的连续性,需要在时空方向上进行插值.由于四叉树同层网格在空间上步长相等、且相邻时间层上时间步长相等,所以在时空方向上采用平均值插值,这样做的优点是时空方向插值既可以达到二阶精度,又可以简化计算过程,节省资源,对任意边界加密下的网格,格子Boltzmann方法的实现比较容易.为了充分利用均匀笛卡尔网格的优势,物面边界的处理采用了速度修正法,与传统浸入边界-格子Boltzmann方法中的直接力法、动量交换法相比,无滑移边界条件得到了较好的保证.同时,把大涡模型加入到浸入边界-格子Boltzmann方法中,实现了在四叉树网格数据结构、边界处理技术、大涡模拟几种模型相结合下,绕障碍物的较大雷诺数流动的模拟.通过不可压缩粘性流中圆柱绕流算例验证,结果与其它方法结果吻合良好.     

基于动态边界控制的优化方法在高升阻比气动外形优化中的初步应用

本文是利用动网格技术和有限体积法来研究基于动态边界控制的气动优化的方法.该优化方法是通过在求解非定常流动的过程中计算优化的目标函数随设计变量变化的梯度.而利用参数化的结构动网格技术能有效、快速地生成与优化过程出现的外形对应的网格并进行计算,具有较高的效率.本文以典型的高升阻比外形优化为例,用非定常NS方程为流动控制方程进行优化,取得了初步的结果.