二元矢量喷管内流特性的数值模拟

应用多重网格ＭａｃＣｏｒｍａｃｋ显示格式对不同矢量角下二元喷管的流动进行了模拟分析。计算时,在细网络上采用了雷诺平均Ｎ－Ｓ方程和Ｋ－ε湍流模型方程;在粗网格上则采用Ｅｕｌｅｒ方程。数值计算给出了较清晰的流谱。壁面压力分布的计算值也与实验结果吻合。通过计算结果分析了喷管几何参数对流场的影响。所得结果显示出多重网格算法高效可靠,可用于开展参数研究     

强迫振动法求解偏航阻尼导数

求解了非定常薄层近似Ｎ－Ｓ方程模拟强迫振动条件下的非定常扰动流场,并在Ｅｔｋｉｎ定义下给出偏航阻尼导数Ｃηγ－Ｃηβｃｏｓα的计算方法和数值积分公式,还给出Ｍａ＝１０,α＝９．８４７°,１９．６８°下的计算结果。     

超声速有翼飞行器的无粘分区数值模拟

本文采用ＭａｃＣｏｒｍａｃｋ两步显式格式，用空间推进法求解Ｅｕｌｅｒ方程来对有翼飞行器的超声速绕流流场进行数值模拟，应用特征边条件处理场面边界及装配激波的计算，在翼缘及表面斜率不连续处，采用局部解析处理，提高了算法的可靠性，复杂外形的计算采用分区处理技术，网格由代数形式生成，在大攻角的计算中，通过对横流速度型的修正改善了背风面压力分布特性，解决了推进中断问题，同时模拟出背风面涡。     

固体火箭发动机内流场分区耦合数值计算

采用分区耦合方法计算固体火箭轴对称燃烧室与喷管流场。对于低速的燃烧室流场，选用不可压流的Ｎ－Ｓ方程描述并用ＳＩＭＰＬＥＣ方法数值求解；对于高Ｒｅ数的喷管流场，则采用Ｅｕ－ｌｅｒ方程描述并用ＳＣＭ方法求解。计算时用燃烧室出流为喷管流场提供入口参数，同时用喷管流场压强分布反馈影响燃烧室流动状况。对耦合边界条件处理方法进行了探讨。对典型的侧壁加质燃烧室与喷管流场进行了计算，计算结果揭示了单独喷管流场计算难以反映的喷管收敛段近壁区的低速区域，与已有的燃烧室流场实验结果一致并反映了燃烧室与喷管流场之间的联系，较好地模拟了流动中的物理现象。     

火箭垂直发射燃气流场的数值模拟

数值模拟了火箭垂直发射燃气流场的非定常过程,在计算中采用了二阶精度的ＭＵＳＣＬ格式,在引入单调性限制条件和一定的耗散机制后,该方法适应于Ｅｕｌｅｒ坐标系下的计算。计算结果表明,该方法得到的流场激波结构十分清晰,计算结果和实验规律相一致。     

低密度小喷管流场的DSMC直接模拟

用ＤＳＭＣ方法数值模拟了低密度轴对小喷管的流场流动，获得了相应的结果，并将所得结果与用差分求解Ｎ－Ｓ方程方法获得的结果进行了对比，二者吻合较好。     

FATIGUE CRACK CLOSURE MEASUREMENT OF 2024-T3 SHEET SPECIMEN

ＦＡＴＩＧＵＥＣＲＡＣＫＣＬＯＳＵＲＥＭＥＡＳＵＲＥＭＥＮＴＯＦ２０２４－Ｔ３ＳＨＥＥＴＳＰＥＣＩＭＥＮＷｕＹｉｓｈｅｎｇ（ＩｎｓｔｉｌｕｔｅｏｆＭｅｃｈａｎｉｃｓ，ＣｈｉｎｅｓｅＡｃａｄｅｍｙｏｆＳｃｉｅｎｃｅｓ，Ｂｅｉｊｉｎｇ，Ｃｈｉｎａ，１００...     

求解Euler／Navier—Stokes方程的有限体积龙格库塔方法

提出了一个求解定常Ｅｕｌｅｒ／Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ方程的计算方法,求解的主管方程是积分形式的Ｅｕｌｅｒ方程或雷诺平均的Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ方程,计算方法是有限体积多步龙格库塔方法,通过在格式中加入自适应的二阶加四阶人工粘性项以保证数值计算的稳定性,用空间变时间步长、焓阻尼来加速定常问题的收敛速度。     

M∞→1时圆球绕前体流场的数值模拟

本文根据来流马赫数Ｍ∞选取坐标变换函数，将Ｍ∞→１时的低超声速圆球绕流前体流场变换到矩形的计算区域，忽略粘性影响，采用时间相关法，用ＴＶＤ有限差分格式求Ｅｕｌｅｒ方程的定常解，得到了Ｍ∞＝１．０５、１．０１和１．００５的流场分布。结果与弹道靶的实验吻合较好。     

低速大迎角三角翼旋涡流动的有粘怀无粘数值模拟

通过人工压缩性方法求解了三维定常不可压Ｅｕｌｅｒ方程和雷诺平均Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ方程，应用Ｂｅａｍ－Ｗａｒｍｉｎｇ近似因子分解格式及其对角化形式，采用Ｂａｌｄｗｉｎ－Ｌｏｍａｘ代数湍流模型。计算模型为７０°三角翼，分析了有粘与无粘流动旋涡特性的影响，尤其是对Ｅｕｌｅｒ方程模拟大迎角分离涡、涡破裂的能力及其局限性提供了数值验证，计算结果与实验结果相一致。     

平头增阻再入体俯仰动态特性计算与流动机理分析

采用自由振荡法数值模拟了一种平头增阻再入体的高超声速俯仰动态特性,沿用稳定性导数模型作为动态特性的指标,通过把稳定性导数视为攻角的函数来考虑非线性影响,并应用奇异分解线性最小二乘法辨识稳定性导数。对于本文研究对象,计算结果表明,平头外形在较高马赫数和小攻角范围内存在动不稳定现象,随着再入时马赫数的降低,动不稳定的攻角范围不断缩小,直至演化为动稳定的。文中对平头外形出现这种动态特性现象的空气动力学原因进行了分析,有必要作进一步的确证与研究。     

平／钝头外形俯仰动态特性对落点姿态的影响分析

采用自由振荡法数值模拟了平头、钝头外形的超声速俯仰振荡时间历程,并应用奇异分解线性最小二乘法辨识稳定性导数,得到动导数随马赫数和攻角非线性变化的规律.平头外形在较高马赫数和中小攻角范围内存在动不稳定现象,随着再入时马赫数的降低,动不稳定的攻角范围不断缩小,直至演化为动稳定的;小钝头外形则没有动不稳定现象.蒙特卡洛弹道仿真表明,平头外形初始时刻俯仰角的扰动会经历一个先增长然后衰减的过程,最终落地时俯仰角小于1°,满足设计要求,但其落地速度有所降低,可能影响该外形对靶标的侵彻效果.     

超音速、高超音速飞行器动导数的高效计算方法

将超音速、高超音速气动力工程计算方法综合起来,推广到非定常气动力的计算,并用于求解任意外形飞行器作强迫运动时的非定常气动力.在此基础上进一步发展了一种适合于超音速、高超音速大迎角状态的飞行器动导数计算方法和程序.计算结果与国际标准实验数据作了比较,证明了该程序的正确性和精度.研究了减缩频率和基本迎角对动导数的影响,结果表明基本迎角对动稳定性导数影响较大,而减缩频率在小迎角范围内对动导数影响很小.     

超声速和高超声速机翼俯仰导数—当地流活塞理论解法

本文给出计算超声速和高超声速尖前缘三维机翼俯仰导数的一个近似解析方法。该方法是根据当地流活塞理论和片条理论导得的。在激波附着于机翼前缘情况下,本方法可用于计算任意迎角,颇为一般的平面形状、尖前缘剖面机翼的俯仰导数。中小平均迎角的数值结果与有限的试验数据和其它理论作了比较,示例表明本方法的结果令人满意。     

超,高超声速机翼俯仰导数的欧拉方程摄动解法

应用片条理论，将文献（１）的尖前缘任意翼型的俯仰稳定性导数的计算方法作、推广、用于计算任意平面形状的三维机翼。在各种特殊的情况下，其计算结果与现有的理论都符合得很好。在激波附体的条件下，本文的计算方法可用于计算任意迎角。任意平面形状机翼的超、高速声速俯仰稳定性导数。     

翼型大攻角绕流的数值模拟

以求解二维N-S方程数值模拟NACA0012翼型大攻角状态的可压绕流特性;N-S方程是在贴体坐标系中给出的,以代数方法生成C型网格系统。采用LU-ADI格式计算,为提高格式的稳定性在隐式和显式部分分别添加了2阶和4阶人工粘性项。应用BaldwinLomax湍流二层代数模型模拟了大攻角时前缘分离涡的形成,旋涡对流及其非定常现象。在某些Mach数和攻角下NACA0012翼型的湍流解具有周期性。通过比较,本文数值计算结果同实验及国外相应的数值计算结果基本吻合。     

反流区对复合高速直升机旋翼气动特性的影响

针对复合式直升机高前进比旋翼反流区严重的特点,建立了适合于高前进比旋翼气动特性的分析方法,以H-34旋翼为例计算了该旋翼在高前进比状态下的气动性能,并与已有的风洞试验数据进行对比验证.在此基础上,进一步分析了反流区对高前进比旋翼气动性能以及对桨叶剖面迎角、升力系数和阻力系数的影响.结果表明:反流区越大,对旋翼的气动性能...     

扇翼飞行器翼型附面层控制数值模拟

基于扇翼飞行器翼型特殊的几何形状及流场特性,在原有翼型的弧形槽下方和后缘加装控制阀门,通过调节阀门开启及开启尺寸的大小,利用弧形槽低压涡所产生的吸力对翼型后缘的附面层进行一定的控制,达到增升减阻的效果。通过采用计算流体力学的方法对其机理及阀门开启尺寸的影响进行了详细计算和分析,研究表明当阀门开启的尺寸为10 mm时,修改翼型的最大升力系数、失速迎角及相同迎角下的升力系数和推力系数均大于基本翼型;随着阀门开启尺寸的增大,修改翼型的最大升力系数和失速迎角均减小,但是在失速前,修改翼型在相同迎角下的升力系数大于基本翼型。此方法可以改变先前通过增大横流风扇的转速来提高其气动性能的做法,减小了能量的消耗,增大了整个飞行器的航程,为扇翼飞行器能够早日投入实际运用奠定了一定的理论基础。     

Review of theoretical achievements for starting flow problem for all Mach numbers

The typical behavior of unsteady flow and force evolution in a number of applications, such as aero-elastics, gust-wing interaction, flapping flight and flight maneuvering, can be understood using the starting flow model. Starting flow model is obtained either by setting rapidly an angle of attack for a wing moving at constant speed, or by accelerating a wing to a constant speed while gaining an angle of attack. In the limiting case of impulsively starting flow, the wing is assumed to gain suddenly an angle of attack in an initially uniform flow. Theories have been developed for impulsively starting flow at small angle of attack long before and at large angle of attack only recently, especially for incompressible and supersonic flow. This paper intends to provide a state-of-art overview of the typical flow phenomena, force evolution characteristics and developed theories for impulsively starting flow at any angle of attack and for both lower speed flow (vortex dominated) and high speed flow (compressible wave dominated). This review also provides some new topics that deserve further studies.     

使用RANS/LES混合方法对钝前缘三角翼进行数值模拟(英文)

阐述了建立在一方程湍流模型基础上的分离涡模型(DES)和延迟脱体涡模型(DDES),它们很好地模拟了壁面附近的小尺度流动。该方法在大攻角分离流区域,使用了亚格子应力的Smagorinsky大涡数值模拟,模拟亚音速旋涡流动。使用大攻角下的钝前缘65°大后掠三角翼模型,研究了CFD对涡的发展、破裂以及复杂涡的发展和演变的模拟能力。研究结果说明了方法对模拟分离流动是可行的。