亚、超音速机翼非定常广义气动力系数之间的关系

 <正> 研究弹性飞机的动力学特性应计及非定常广义气动力。为了尽可能保留非定常运动历程,应首先计算由阶跃函数形式下洗产生的非定常气动力函数(指示函数或指数函数),再利用卷积形式的迭加原理处理任意函数形式下洗产生的非定常气动力。用振型法建立飞机弹性运动方程,可将飞机上任一质点的弹性位移近似地表示为     

机翼及振动襟翼的非定常气动力计算

 <正> 1.引盲 研究细长翼及其前缘襟翼大幅度振动所引起的气动力,有助于了解非定常流的一些特点,对未来高机动飞机的设计和飞行控制也有重要意义。计算这种情况的气动力,因涉及大尺度非定常分离流问题,需要能模拟自由涡面变化的位流方法,如常值偶极子格网法〔‘、涡格法山。〕和混合涡法〔‘〕。这些方法及其应用还在发展中〔,,。〕。为拓展位流方法的应用范国,本文就涡格法对单、双三角翼带振动前缘襟翼的气动力进行了探讨     

跨声速机翼非定常气动力的全位势粘位迭代计算

采用Ｃ－Ｈ型网格,守恒型非定常全位势方程的时间精确近似因式分解差分地计算二维,三维的跨声速非定常位势流,用准定常,准二维方法计算边界层位移厚度,通过粘位迭代得到的跨声速翼型,机翼的非定常气动力,所得结果与实验数据吻合很好。     

带操纵面机翼的非定常跨声速流数值解法

本文介绍一种带操纵面机翼的非定常跨声速流的有限差分计算方法。采用的方程是三元非定常跨声速修正小扰动位势方程,使用时间积分法,求解格式是近似因式分解交替方向隐式(ADI)格式。使用这种方法计算了F-5机翼在来流马赫数为0.9和0.925时的定常气动力和操纵面振荡的非定常气动力。计算结果与国外的NLR试验结果和XTRAN3S方法的计算结果进行了比较,表明计算是成功的。     

二维涡轮级的定常/非定常数值模拟

为了解定常与非定常计算差别对工程应用参考价值的影响,本文通过求争二维Navier-Stokes方程,针对一级涡轮叶栅分别进行了定常／非常数值模拟。结果发现,对于叶片表面压力分布,第一排导叶定常计算与非定常计算得到一致的结果,但是在动叶上差别较大,尤其在叶排小间距时。当叶排间距拉大后,定常计算结果趋近于非定常计算结果时间平均值。     

超声速和高超声速翼型非定常气动力的一种近似计算方法

本文叙述当地流活塞理论,根据这一理论导得超声速和高超声速尖前缘翼型的压强系数和振动导数公式。数值结果与实验数据和其它理论的结果作了比较。实例表明在活塞理论失效的情况(M_∞δ_(max)>1)下,本方法仍能给出合理的数值结果。     

氢气超音速燃烧非定常过程数值模拟

用有限差分法求解二维Ｎ－Ｓ方程和组分守恒方程,使用代数涡粘性湍流模型,两步化学反应模型和ＭａｃＣｏｒｍａｃｋ时间分裂法,模拟超音速气流中横喷氢气自动着火及火焰传播过程,清晰地描述了喷氢、火焰传播到主流、稳定燃烧过程中,温度、压力、速度和油气当量比的变化。计算结果表明,高温附面层和喷嘴前回流区能促使氢气自动着火,进口气流的静温是火焰向主流传播的重要影响因素。计算还证明了喷嘴后氢与空气两股流呈现平行流动,其间存在着强烈的湍流交换和化学反应边界层,超音速燃烧效率则主要取决于两股气流的混合速率。     

动态网格快速生成及组合体跨音速非定常气动力计算

采用代数方法、局部变形技术及快速动态弹性变形技术生成三维翼/体组合的贴体动态网格。用时间精确近似因式分解差分方法求解守恒型非定常全位势方程。考虑附面层影响,通过粘/位迭代得到跨音速翼/体的非定常气动力。所得结果与实验数据相比相当吻合。     

二维轴流压气机转—静非定常Euler方程数值计算

将Euler方程写在包括转子和静子的绝对坐标系下,通过在转—静之间设计一个可以随时间扭动的平行四边形网格计算了二维轴流压气机转—静非定常流场。由于整个方程统一写在动网格下,转—静之间没有使用插值。      

COMPACT FINITE VOLUME SCHEMES AND THEIR APPLICATIONS

ＣＯＭＰＡＣＴＦＩＮＩＴＥＶＯＬＵＭＥＳＣＨＥＭＥＳＡＮＤＴＨＥＩＲＡＰＰＬＩＣＡＴＩＯＮＳＪｕＨｏｎｇｂｉｎ（ＳｃｈｏｏｌｏｆＰｏｗｅｒａｎｄＥｎｅｒｇｙＥｎｇｉｎｅｅｒｉｎｇ，ＳｈａｎｇｈａｉＪｉａｏｔｏｎｇＵｎｉｖｅｒｓｉｔｙ，Ｓｈａｎｇｈａｉ...     

高温燃气进入舵机舱过程仿真与流动机理分析

采用计算流体力学(CFD)方法研究了火箭发动机工作拖尾段高温发动机燃气进入舵机舱的物理现象。结合导弹实际飞行弹道参数变化特点和超声速流场扰动不向前传递的空气动力学理论,提出了简化而不失真的非定常流场仿真方案,显著缩短了仿真周期;复现了某型导弹实际飞行时舵机舱先被\"抽气\"再进高温燃气的动态过程,并分析了高温发动机燃气进入舵机舱的流动机理,即在发动机工作段,导弹底端面压强低于舵机舱内压强,舵机舱被\"抽气\",在拖尾段随着燃烧室总压降低,喷口附近的马赫盘向导弹底端面移动,使导弹底端面压强增大且高于舵机舱内压强,高温燃气进入舵机舱烧毁电路致使导弹折断;明确了某型导弹折断故障产生的诱因,提出了改进措施和检测方法,并得到了大量飞行靶试的验证,解决了舵机舱热防护结构可靠性问题。     

基于动态非结构重叠网格法的直升机前飞非定常流场数值模拟研究

 发展了非结构重叠网格技术,提出适合于直升机旋翼运动的动态非结构重叠网格算法,该算法使得在采用非结构重叠网格模拟直升机流场时实现程序自动处理网格,不需任何人工干预。流动控制方程采用非定常可压缩Euler方程,空间离散采用Jameson有限体积法,时间推进采用双时间步长法,为加速收敛,在双时间步长内循环迭代中采用预处理技术。为验证本发展的预处理解算器以及动态非结构重叠网格的算法的可信性和有效性,开展一系列的数值实验,模拟了不同速度范围的NACA0012翼型绕流,Caradonna旋翼悬停和前飞流场以及Robin直升机机身-旋翼组合的前飞非定常流场,结果表明采用本文算法得到的结果是可信的,网格处理是自动高效的,可以应用于直升机全机流场模拟,为直升机设计提供可信指导数据。     

网格之间流场信息传递的高精度方法和验证

理论分析了插值方法引入误差的原因,提出一种利用时间推进方法和移动网格技术实现网格之间流场信息传递的思想,并构造计算方法。在介绍新的信息传递算法原理的基础上,根据三维情况下传递过程需要满足的时间推进步长稳定性条件,推导出最小推进步数的要求。最后,通过对激波诱导物体运动的非定常流动模拟,比较了信息传递算法和线性插值、加权最小二乘插值方法,计算结果表明,与插值方法相比,信息传递过程没有引入误差,能有效抑制非物理振荡,保持流场计算精度。     

一种新的动网格方法及其应用

 介绍了一种新型的动网格方法。该方法以Delauney图为基础进行整个网格的移动变形：首先是生成满足Delauney特性的Delauney图;然后把对网格移动的特性映射到Delauney图的移动上;最后把网格点在Delauney图以映射的方式进行重新定位而得到新的网格点坐标,从而得到移动后的网格。该方法不需要迭代,因而效率高,并且能适用于任意拓扑结构的网格。以翼身组合体的机翼颤振、翼身融合体的变体过程的非定常流动等为例来说明该方法的高效性及适用面广等特点。     

远场假设对喷流噪声预测中格林函数求解的影响

目前基于雷诺平均Navier-Stokes（RANS）的喷流噪声预测方法在格林函数求解时, 为简化求解过程,通常对喷流流动做平行流假设, 对观测点做远场假设。随着格林函数求解方法发展,近年来的研究表明平行流假设对下游观测点格林函数的计算会引起较大偏差,而目前远场假设对格林函数求解的影响仍不清楚。为研究远场假设对喷流格林函数求解的影响,以二维喷流为例, 采用计算气动声学方法（CAA）分别数值求解了观测点远场假设条件与实际条件下90°~150°方向喷流内伴随格林函数,进而分析远场假设对格林函数求解的影响。研究结果表明,对于不同方向的观测点,由观测远场假设导致的伴随格林函数求解偏差不尽相同,且对于越靠近喷流中心线方向的观测点,远场假设导致的偏差越大,其中150°方向观测点,采用远场假设后,格林函数计算结果最大偏差达到-15 dB以上。因此,对于靠近喷流中心线方向的噪声观测点而言,为避免预测偏差,应采用实际观测条件求解喷流格林函数。     

一种计算非定常二维流动的无网格算法

主要目的是发展一套求解非定常流动的无网格算法。计算区域的离散方面,提出了一种按区域进行填充布点的点云自动生成方法;发展了一种点云的运动技术来实现离散点云对物面边界的随体运动;在点云离散的基础上,采用最小二乘法求解矛盾方程的方法来求取空间导数,进而获得数值通量;采用双时间方法进行时间离散推进,其中物理时间迭代采用二阶隐式格式,伪时间迭代采用四步龙格一库塔显式格式,为了加速收敛,在伪时间迭代中采用了当地时间步长和隐式残值光顺等加速收敛措施。最后,利用本文算法模拟了NACA0012翼型和NACA64A010翼型的跨音速非定常流动,并将计算结果与实验测量结果进行了对比分析,验证了上述方法的正确性和实用性。     

微型飞行器的仿生流体力学——昆虫前飞时的气动力和能耗

 用数值模拟方法研究了昆虫前飞时的气动力和需用功率。由N S方程的数值解提供速度场和压力场,从而得到涡量、气动力和力矩 (惯性力矩用解析方法计算 )。基于流场结构,解释了非定常气动力产生的原因;基于气动力和力矩,得到需用功率。悬停飞行中揭示出的 3个非定常高升力机制 (不失速机制,拍动初期的快速加速运动,拍动后期的快速上仰运动 )在前飞时仍然适用 (即使在快速前飞时,V∞ =2～ 2.5m/s,失速涡也不脱落 )。在低速飞行时 (V∞ ≈ 0.5m/s)平衡重量的升力既来自于翅膀的下拍运动也来自于上挥运动,并主要由翅膀的升力贡献;克服身体阻力的推力主要来自于翅膀的上挥运动,由翅膀的阻力贡献。在中等速度下 (V∞ ≈ 1.0m/s),升力主要来自于下拍运动,其中一半由翅膀升力贡献,一半由翅膀阻力贡献;推力主要来自上挥运动,也是一半由翅膀升力贡献,一半由翅膀阻力贡献。在快速飞行时 (V∞ ≈ 2.0m/s),升力主要来自于下拍运动,主要由翅膀阻力贡献;推力来自上挥运动,主要由翅膀升力贡献。悬停时,下拍和上挥做功同样大;前飞时,下拍做功较上挥大得多 :V∞ =0.5,1.0和 2.0m/s时,下拍做的功分别是上挥的 1.6,2.6和 3.5倍。     

叶轮机非定常气动设计的缘线匹配技术

本文提出一种新的叶轮机气动设计自由度——缘线匹配技术。缘线匹配是指相邻两叶排中前排叶片后缘线与后排叶片前缘线的空间关系,更准确地说是其相位角关系。叶片尾缘线表征了叶片出口气流的分布型,它囊括了尾迹、激波、二次流等所有影响,另一方面,叶片前缘线代表了其对上游流场的势干扰。定常气动设计体系只能部分考虑叶片缘线的相互影响,大部分则被忽略了,然而在非定常框架下,相邻叶排缘线所代表的相互影响是显著的。本文首次提出了叶轮机非定常设计的缘线匹配技术,并通过理论分析及直列涡轮叶栅的非定常数值模拟结果展示了缘线匹配技术在未来提高叶轮机气动性能、气弹性能、气动噪声和热传导性能上的潜在能力。作者认为,缘线匹配技术将使叶轮机的非定常气动设计具有真实而可操作的内容,是叶轮机非定常气动设计的核心技术之一。      

基于非结构动网格的非定常流数值模拟方法

给出了基于非结构动网格的非定常流动问题的数值模拟方法,重点说明两种动态非结构网格方法的原理以及在非定常流模拟中的算法.一种是基于非结构网格的重叠网格方法,该方法结合了非结构网格方法和重叠网格方法的优点,以物面距为准则参数对网格节点进行分类,进而建立重叠网格的网格问插值关系,通过采用动态非结构重叠网格算法,可以处理多体间...