旋成体表面沟槽减阻试验研究

在旋成体模型表面粘贴铝基沟槽蒙皮,研究沟槽表面对模型阻力的影响。给出基于来流单位雷诺数的沟槽尺寸无量纲公式,用于设计沟槽齿高和齿宽,依据试验风洞的单位雷诺数范围,选定3种尺寸的沟槽进行减阻试验。结果表明,在旋成体表面湍流流动区域顺流向布置沟槽,当沟槽齿高和齿宽的无量纲值h＋、s＋均小于25时,在小迎角下具有减阻作用,当s＋约为15时减阻效果最明显,可减小约3%～4%的阻力。     

细长旋成体背风面非对称流场的DES模拟

利用脱体涡模拟DES方法对细长旋成体非对称绕流进行数值模拟,比较了基于S-A湍流模型的DES与RANS方法对背风面非对称分离涡的数值模拟能力,发现在大迎角非对称涡未破裂情况下,DES与RANS均能模拟出与实验相符的非对称绕流;当迎角增大,背风面非对称涡发生破裂时,RANS无法准确捕捉到背风面流场的非定常性,而DES能准确预测非对称涡的飘起与破裂,并与实验值接近。计算结果表明:与RANS方法相比,DES方法具有更好的模拟大攻角,大分离流动的能力,尤其在非对称流场的分离涡破裂模拟方面具有明显优势,能够更真实模拟出细长旋成体背风面分离涡破裂之后的非定常流动特征。     

Re数对细长旋成体非对称涡及气动力特性影响的实验研究

通过对细长拱柱旋成体大迎角绕流不同截面测压结果分析，探讨了绕流舶数对非对称涡结构和气动特性的影响，得出舶数不仅影响分离线位置和绕流流态结构，而且影响边界层的绕流特征及其分离涡的强度，非对称性的出现与细长体两侧边界层的绕流特征和分离涡的强度不等存在密切的关系。特别是在同种流态下，两侧边界层的绕流特征和分离涡强度不等是造成侧向力的主要原因；在两侧不同的流态下，转捩不对称是产生大侧向力的主要原因。     

粗糙带起始位置对旋成体绕流过临界雷诺数模拟效果的影响

使用低速闭口风洞,对头部带有确定人工扰动的尖拱型旋成体,在模型两侧θ=±52.5°粘贴条状粗糙带的条件下,粗糙带的轴向起始位置在x/d=0.75之前变化,固定迎角为40°,在雷诺数ReD=0.67～9.4×105变化范围内,通过表面测压研究粗糙带起始位置对模拟过临界Re数下旋成体绕流的影响.结果表明大迎角下,人工转捩带的轴向起始位置会影响到头部部分截面的压力分布,但对后体的影响很小.粗糙带的起始位置越靠后,所受影响的轴向位置越远.通过分析得出粗糙带起始位置与受影响的轴向位置x/D之间的关系曲线,由此确定模拟过临界雷诺数下流动时合适的粗糙带起始位置.     

细长旋成体大迎角非对称流动特性的试验研究

简要介绍了细长体大迎角流动非对称性的试验结果，分析了Re数、湍流度和安装条件等因素对大迎角流动非对称性的影响，探讨了大迎角流动非对称性的产生机制。文章最后着重阐述了该研究得出的一些支持细长体大迎角流动非对称性产生机制的空间动力不稳定性观点的理由。     

船尾形状对旋成体马格努斯效应的影响

尖头旋成体和船尾形状是子弹、炮弹及火箭弹等抛射体上常用的布局形式。研究表明船尾布局具有减小底部阻力、增大射程的作用,但此时旋成体的马格努斯效应增大,对运动稳定性产生不利影响。为了解释这种流动现象,对三维旋转弹流场进行了数值模拟,对从亚声速到超声速下的旋成体马格努斯力和力矩进行了分析,重点对标准形状和船尾形状两种底部进行了比较。结果表明,相对于标准形状,在所有来流下船尾形状都起到了增大马格努斯效应的作用,并且马格努斯力和力矩与船尾角成正比。为了揭示其流动机理,选择代表性计算状态对两种布局马格努斯力矩系数分布、边界层厚度分布和边界层位移厚度分布进行了对比分析,结果表明,在亚跨声速下船尾马格努斯效应由绕拐角的加速流动引起,使当地压力系数幅值增大;在超声速下船尾马格努斯效应由船尾段的气流膨胀引起,使旋成体左右两侧的边界层位移厚度畸变增大。上述两种效应都使马格努斯力矩增加,对于亚声速流动来说,该效应发生在柱段与船尾段连接位置;对于超声速流动来说,该效应发生在连接点以后的船尾段上。当来流速度在声速点附近时,上述两种效应都可能发挥作用,使船尾形状的旋成体马格努斯效应大幅增加。     

旋成体大迎角分离流的Euler方程数值模拟

用三维Euler方程计算旋成体的大迎角绕流,我们发现:亚声速自由流时,计算得到的是无分离流动,与实际的分离流不相符合,因而计算的法向力大大小于有分离的实验值;而当自由流为超声速,且横向自由流M数较高时,由于背风表面出现激波,无粘Euler方程的数值计算能自动捕捉到背风面上的分离涡,计算法向力与实验值接近。为了对亚声速自由流用Euler方程数值模拟大迎角旋成体背风面上分离涡,根据问题的物理图像,在假设的分离线上,提出强加Kutta条件的方法,成功地计算出有分离流动,使计算的法向力与实验值接近。     

细长旋成体大迎角绕流的背涡结构与压力脉动特性

在亚临界流动范围内,通过尖拱细长旋成体在无侧滑状态下的脉动压力测量,对大迎角下细长旋成体非对称背涡结构沿轴向的压力脉动特征进行了研究分析,揭示了压力的脉动幅度与频率沿旋成体轴向的变化规律及其与背涡结构之间的对应关系,并对压力脉动特征沿轴变化的原因进行了初步探讨.     

轴对称旋成体流动结构的数值研究

利用N-S方程模拟了旋成体的绕流过程。旋成体流场由不同层次的涡结构组成,包括主涡、二次涡t、ertiary涡、螺旋涡、羊角涡、U型马蹄涡。文中清晰地展示了这些结构,简要分析了各类结构的特性和功能,强调了亚结构和三维结构在流态发展与流场演变中的作用。     

二维欧拉方程在非结构网格上的自适应有限元算法

发展了一种易于和自适应网格加密技术相结合的非结构三角形网格生成技术。这种方法的特点是域内节点可任意规定;根据欧拉方程解的误差估计,在流场某些区域增加节点后,能重新快速生成新的具有高质量三角形单元的非结构网格。欧拉方程采用Ｇａｌｅｒｋｉｎ空间离散和Ｒｉｃｈｔｍｙｅｒ两步显式时间推进相结合的有限元法求解,并以沿流线熵的变化量作为物理判据进行网格自适应加密,提高解的精度。     

用有限解析法解N—S方程

本文介绍有限解析法的一种新构想,并将本方法用于求解不可压缩粘性流体流动的N-S方程,采用贴体坐标,雷诺数达1000,计算结果表明本方法是有效的。正确的。     

用保角转绘坐标解二维和轴对称可压流Euler方程

本文推导出用保角转绘正交曲线坐标表示的二维与轴对称Euler方程,方程中的未知量以曲线坐标方向的速度分量表示,并有不含未知量导数的源项。方程组的形式与直角坐标系中的接近,因而简化了计算。本文将该方程组用于计算二维任意翼型和旋成体轴对称绕流,采用分辨率较高的Eberle特征平均显式有限体积差分算法求解,在不同区域中用不同的显式推进步数加速收敛。本方法已分别对翼型和旋成体轴对称绕流编制了计算程序,可用于自由流为亚声速、跨声速和超声速的绕流计算。大量计算结果表明,本方法精度好,激波分辨率高,编程序简单。     

论有限元素理论中SUPG方法与PG方法在解高速粘流时的等效性

对流场中含激波问题,若用有限元法来数值模拟,目前较有效的方法是SUPG方法和PG方法。本文讨论了这两种方法在数值求解N-S方程组时具有等效性,并对它们的应用前景作了简单的评估。     

不可压Navier—Stokes方程组的SUPG有限元数值解

本文从定常不可压NaVier-Stokes方程组出发,构造了SUPG加权剩余公式。为保证数值解的精度,本文对速度取八节点插值,保留了摄动项中的二阶导数项。从用本文方法所做的算例来看,计算结果是令人满意的。     

用三维欧拉方程计算机身与机翼绕流

本文用Jameson的三维欧拉方程有限体积法、四步Runge-Kutta时间推进格式,计算机身和大后掠细长机翼的三维可压缩绕流。对钝头机身,用C-O型网格;尖头机身与大后掠细长机翼用H-O型网格。本文介绍钝头旋成体、带座舱前机身和三角翼绕流的计算结果,显示流场等值M线分布,计算压强系数分布与实验比较,以及三角翼大迎角分离涡等。由于用了隐式残值光顺等加速收敛措施,有效地减少推进步数,节省机时。     

梁壳组合结构的一种非线性有限元计算方法

提出了梁壳垂直组合结构在有限元计算中存在的问题,并以几何大变形问题为依托提供了解决该问题的计算原理和有限元方法。数值算例证明了该理论和方法的正确性。该法不仅适用于梁壳组合结构的大变形非线性计算,也适用于静动力小变形问题的线弹性分析计算。     

无粘跨声速绕流的有限元法和界积分方程的某些进展

本文对无粘跨声速绕流的有限元法和边界积分方程法的目前进展,特别是对提高跨声速流算法效率起着主要作用的成果,例如普适元素和块结构有限元网格生成技术有限元解法的多重网格技术、SUPG有限元法、GRMES算法和可压流的完全边界积分方程法等方面作了回顾和评介。     

轴对称结构的有限元离散法

介绍一种利用高精度等参环元对轴对称结构进行有限元离散的方法（也可用于平面结构的有限元离散）。这种方法可方便地离散由多种材料构成的复杂轴对你结构，可单向和双向对局部单元加密，可对边界移动和结构间隙问题作有效处理，且可通过离散程序自动优化单元排序以减小总刚度矩阵的规模。离散后的结构可直接用于有限元分析。实践应用表明，该方法离散快速、准确。     

非结构网格下Euler方程的高分辨率高精度解

提出了一种非结构网格下求解Ｅｕｌｅｒ方程的高分辨率高精度迎风格式。以Ｒｏｅ的矢通量差分分裂为基础,吸收了ＮＮＤ格式的优点,使其具有捕捉激波和滑移线的良好性能;在时间方向上采用Ｊａｍｅｓｏｎ的Ｒｕｎｇｅ－Ｋｕｔａ积分,并结合局部最大时间步长和残差光滑技术加速收敛。最后成功地完成了二维平板激波反射、跨音速Ｌａｖａｌ喷管内的流动和ＧＡＭＭ超音速前台阶绕流等算例,显示了该方法的有效性