亚音速机翼和机身大迎角下非线性气动载荷的两种计算方法

本文计算无粘亚音速流中任意平面形状薄翼和旋转体机身在大迎角下计入气流分离时的定常非线性气动载荷。采用两种方法:(1)定常流的松弛法;(2)陡然起动流的时间相依法。方法(1)对旋涡强度和自由涡的位置进行松弛迭代,取得了收敛的结果。方法(2)计算机翼陡然起动后一系列时刻的旋涡强度和自由涡的位置,直到结果趋近于常数解。亚音速流的压缩性效应是假设小扰动计算的。两种方法计算的结果很接近,和可以找到的实验数据也很接近。     

三维Euler方程的自适应非结构网络计算

介绍一种三维自适应非结构网格的Euler方程求解技术。采用有限体积中心差分的格点格式对Euler方程进行空间离散,多步龙格-库塔时间推进,并采用当地时间步长、焓阻尼修正等加速收敛。采用Delaunay三角剖分技术生成三维流场的四面体网格。网格自适应技术中的网格局部加密也采用Delaunay三角剖分技术。最后用ONERAM6机翼的跨音速绕流计算说明本文方法的有效性。     

三维复杂域的Delaunay三角剖分

介绍一类自动生成内点的Delaunay三角剖分算法 ,其中包括重心插点法 ,外心插点法和Voronoi边插点法。该类方法适用于二维和三维复杂计算域的非结构网格生成 ,所需人工干预少 ,使用方便 ,计算效率高 ,生成的网格质量也很好。文中并给出几个算例说明本方法的应用。     

Delaunay非结构网格生成之布点技术

Ｄｅｌａｕｎａｙ 非结构网格生成技术是 Ｃ Ｆ Ｄ 领域广为应用的有效技术之一。如何布置流场中的计算点是 Ｄｅｌａｕｎａｙ 三角剖分的一个关键问题。本文以二维非结构网格生成为例， 介绍三种布点策略， 即重心布点、外接圆心布点和 Ｖｏｒｏｎｏｉ 边布点， 算例表明三种方法均可生成满意的非结构网格， 其中 Ｖｏｒｏｎｏｉ 边布点方法生成的网格质量更好一些。三种方法均可方便地推广到三维问题。     

三维Euler方程非结构自适应网格投影和光顺技术研究

研究了三维Euler方程非结构自适应网格的投影和光顺技术.以ONERA M6机翼为计算模型,非结构网格为计算网格,守恒变量的范数为自适应加密准则,对采用一般自适应计算策略得到的网格和数值结果进行了分析,并发展了一种考虑网格投影和光顺的自适应计算策略.数值试验结果表明:通过对自适应所生成网格的投影与光顺,实现了自适应计算过程中的几何保真性,获得了相对光滑的压力分布曲线.     

非定常薄翼脱体涡数值计算

本文考虑无粘、不可压缩流中平板薄翼在大迎角下气流沿边缘分离出现脱体涡时的非定常气动载荷数值计算。采用离散涡格法,从某一初始备件开始,计算一系列时刻机翼的非线性气动载荷以及脱体涡的几何形状,从而了解脱体涡的形成过程。当保持边界条件不随时间改变,则随着计算时刻的增加,结果趋近于稳态解。算例计有矩形翼、箭形翼。     

三维非结构网格的欧拉方程聚合多重网格法

介绍一种基于三维非结构网格的聚合多重网格技术,用于求解Euler方程定常解,加速收敛效果十分明显。求解Euler方程的基本算法采用有限体积中心差分格点格式。聚合多重网格技术按一定的规则将细网格的若干控制体积聚合在一起,组成一个较粗的网格。数值实验表明本文多重网格法十分有效。     

机翼大攻角气动力数值计算

本文描述了一个计算任意平面形状的机翼在无粘、亚音速、大攻角情况下的气动力数值方法,计算大攻角下气流分离时的定常非线性气动载荷。考虑翼面具有脱体涡分离,对旋涡强度和自由涡的位置进行松驰透代,取得收敛解。计算结果与实验数据很接近。     

复杂组合体气动力的高阶面元法数值计算

本文论述了采用高阶面元法计算包括机身、机翼、尾翼、挂架、外挂物在内的复杂组合体的亚音速气动力。复杂组合体分解为若干部件,每一部件以若干网片描述,并剖分为有限数目的面元。本方法中以一个或多个平面多边形近似曲面面元。在每一个面元上布置二次偶极子奇性分布,利用边界条件建立确定偶极子强度的线性方程组,从而接得偶极子强度并计算所需要的气动力参数。