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

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

二维任意域内基于节点的局部网格生成算法

凸域内基于节点的局部网格生成算法,克服了基于节点的有限元方法的网格生成可能产生的不一致性。将该基于节点的局部网格生成算法的适用范围拓展到二维任意域。另外,提出了通过使用约束Delaunay路径来划分任意域的区域划分算法,该算法使得在并行实现网格生成的过程中各处理器之间无需通信,从而大大提高了节点给定情形下有限元方法网格生成的并行效率。     

基于Euler方程结合附面层修正方法的应用研究

采用欧拉加附面层修正的方法在结构网格上对民用飞机翼身组合体流场进行了数值模拟研究.采用一种适合翼身组合体模型的的网格生成技术来生成计算用的结构网格,该技术是先采用一种基于Higenstock源项修正的椭圆型偏微分方程生成二维网格,然后沿机身轴向合理的光滑性调整成三维网格.在求解Euler方程时采用中心离散格式和显式时间推进,并利用附面层修正的方法来考虑粘性的影响.为了验证此方法的正确性,本文对DLR-F4翼身组合体进行了数值模拟,结果表明此方法和试验结果能够很好的吻合.     

矩形/非结构混合网格技术及在二维/三维复杂无粘流场数值模拟中的应用

建立了一套模拟复杂无粘流场的矩形／非结构混合网格技术,其中非结构网格仅限于物面附近,发挥非结构网格的几何灵活性,以少量的网格模拟复杂外形;同时在以外的区域采用矩形结构网格,发挥矩形网格计算简单快速的优势,有效地克服全非结构网格计算方法需要较大内存量和较长ＣＰＵ时间的不足。混合网格系统由修正的四分树（２Ｄ）／八分树（３Ｄ）方法生成。本文将ＮＮＤ有限差分与ＮＮＤ有限体积格式有机地融合于混合网格计算,消除了全矩形网格模拟曲面边界的台阶效应,同时保证了网格间的通量守恒。该算法被成功地运用于ＮＡＣＡ００１２翼型的跨、低超声速绕流,二维弯管内运动激波的非定常流,以及三维捆绑火箭的超声速绕流的数值模拟。计算结果表明本方法在模拟复杂无粘流场方面的灵活性和高效性。     

保角变换在复杂外形网格生成中的应用

保角变换在二维翼型网格生成中已经得到了大量的应用。本文描述保角变换方法生成三角翼翼身组合体，全机标模的计算网格，最后，对实用 残击机外形的网格生成也做了有效的尝试。     

固体火箭发动机绝热层温度场的有限元计算方法

利用有限元法计算了固体火箭发动机绝热层在移动边界条件下的二维温度场.采用碳化层-热解面-原始材料的二维碳化烧蚀模型;推导了将热解气体对流项作为源项的有限元计算方法;采用当量对流换热系数和当量热流的方法处理复杂边界条件.采用无限插值法获得移动边界条件下的三角形网格,提高了网格生成速度和网格质量.计算结果表明,利用有限元法计算固体火箭发动机绝热层的温度场收敛性和稳定性都较好.      

三维多块结构椭圆型贴体网格的生成

通过求解椭圆型方程生成三维贴体网格，采用分块技术生成多块结构网格体系。通过实验验证了方法对任意复杂网络区域生成三维多块结构网格能满足流场计算的要求，其质量是相当满意的。     

数值保角转绘新技术及应用

本文发展了一种系统的一次性数值保角转绘方法,可以将一含多个角点的任意曲边多边形的外部区域,一次保角转绘为单位圆的外部区域。由于方法中应用了快速富氏变换技术(FFT),计算速度很快。这种方法可应用于计算绕任意物体的二维不可压位流,计算二维物体的表观质量系数,生成有限差分计算中所需的正交保形贴体坐标网格。对几种外形的计算结果表明,在网格生成方面,所需计算机时与Thompson方法相比,少一个数量级,网格正交性处处保证。此外,对许多问题,变换平面内的主管方程特别简单。     

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

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

二维DRAGON网格自动生成技术的研究

对重叠网格技术中新出现的Cut-Paste方法进行了研究和改进,结合网格洞边界标识的缝合与重组织技术,首次将其应用在DRAGON网格生成中,提出了二维DRAGON网格的自动生成方法,最后编制出通用程序,给出了多体网格的复杂示例.实际计算表明,本文方法具有自动化程度高、可靠性好、DRAGON网格质量高等特点.     

吸气／喷气叶栅的三维贴体网格生成技术

对某吸气/喷气叶栅的复杂的多联通域边界,采用孔斯双参数样条曲面模型,通过一维问题边界条件控制疏密度,选择合理的参数来强调正交性;分别在双参数域和绝对坐标系下求解微分方程,得到二维、三维网格。应用上述方法,为某叶栅叶片前缘和端壁开多排孔方案生成了三维贴体网格,获得了理想的效果。     

基于非结构动网格的非定常气动力计算

对用于非结构动网格生成的弹簧近似方法进行了研究。通过采用顶点弹簧方法,分析研究了弹簧倔强系数的取值,同时通过引入挤压倔强系数和边界修正,对标准弹簧近似方法进行了改进。改进后的方法可以大大提高网格变形能力和网格质量。应用本文发展的非结构动网格生成方法耦合求解基于Arbitrary Lagrangian-Eulerian(ALE)描述的欧拉方程,计算了谐和振动NACA0012翼型及矩形机翼的跨声速非定常气动力,计算结果与参考文献提供的结果及实验结果吻合良好。     

薄壁结构的加筋布局优化设计

采用拓扑优化方法研究了薄壁结构的加筋布局优化问题.考虑到航空航天领域薄壁加筋结构的大量使用与结构形式的复杂性,提出了一种适用于有限元自由网格剖分的加筋设计新方法--几何背景网格法.该方法一方面通过定义几何背景网格的大小,实现了加筋设计域内任意离散网格沿加筋高度方向的布局参数化定义;另一方面,通过背景网格曲线坐标系下的定义,实现了三维曲面薄壁壳结构的加筋布局设计.采用该方法分别对平面结构与曲面薄壁结构进行了以结构刚度最大为目标的加筋优化设计.数值结果表明,所提出的方法能有效获得合理的加筋布局.     

薄壁结构的加筋布局优化设计

采用拓扑优化方法研究了薄壁结构的加筋布局优化问题。考虑到航空航天领域薄壁加筋结构的大量使用与结构形式的复杂性,提出了一种适用于有限元自由网格剖分的加筋设计新方法——几何背景网格法。该方法一方面通过定义几何背景网格的大小,实现了加筋设计域内任意离散网格沿加筋高度方向的布局参数化定义;另一方面,通过背景网格曲线坐标系下的定义,实现了三维曲面薄壁壳结构的加筋布局设计。采用该方法分别对平面结构与曲面薄壁结构进行了以结构刚度最大为目标的加筋优化设计。数值结果表明,所提出的方法能有效获得合理的加筋布局。     

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.     

网格生成和数值模拟的讨论

讨论了多种网格生成方法及其数值模拟。给出一种生成多部件复杂外形飞行器分区网格和Ｅｕｌｅｒ方程分区解。讨论了前缘分离涡为主体的旋涡流场。作了跨声速流分区的并行计算的初步尝试。给出一种二维非结构网格的生成方法。     

二维结构/非结构混合网格的生成及其流场模拟

采用分区求解方法在二维结构/非结构混合网格中对粘性绕流作了数值模拟。利用层面推进法生成绕物体的O形结构网格,利用推进阵面法在结构网格外生成非结构网格。在结构网格区内,采用B-L湍流模型,求解雷诺平均N-S方程。在非结构网格区内求解Euler方程。两区的解通过在交界面上满足通量守恒条件进行耦合。数值算例表明了方法的有效性和实用性。     

多面体网格在CFD中的应用

介绍了多面体网格技术在计算流体力学(CFD)中的应用.考虑到粘性计算的需要,将混合网格技术引入多面体网格中,发展了一套多面体网格生成及求解技术,并用于M6机翼、DLR-F6复杂外形的气动特性模拟中.计算结果表明,多面体网格技术具有准确度较高、效率较高及适用性广等优点.     

用非结构动网格技术计算确定再入飞行器配平攻角

提出了用数值计算直接确定飞行器配平攻角的自由振动法.该方法采用非结构动网格技术,使飞行器在气动力作用下绕通过质心的轴线自由振动,最终衰减到配平攻角.构造了一种MUSCL类型有限体积格式,求解考虑了动网格效应的三维非定常Euler方程,采用弹簧近似技术实现非结构网格的变形运动.作为例证,计算确定了载人飞船返回舱的配平攻角,计算结果与实验结果接近.该方法外形适应能力强,效率较高,具有较好的应用前景.