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

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

共轴式双旋翼悬停流场和气动力的CFD计算

考虑到共轴式双旋翼流场特征高度复杂、桨叶承受非定常气动载荷的特点,为更好地预测共轴式双旋翼的气动特性,把非结构嵌套网格方法和网格的自适应技术相结合,发展了一套适合于共轴式双旋翼流场数值计算的求解器.在该求解器中,采用非结构嵌套网格方案来模拟桨叶之间存在的相对运动,自适应网格技术用来捕捉尾迹对流场和气动特性的影响,求解惯性坐标系下的非定常N-S主控方程来模拟流场的非定常特性.应用该求解器,首先计算了有试验结果可供对比的一副试验旋翼的诱导速度场分布,在此基础上,计算了共轴式双旋翼的桨尖涡轨迹和拉力分布特性,并与单旋翼的计算结果进行了对比和分析,得到了一些有意义的结论.     

单级跨音速压气机三维非定常粘性流动计算

本文以数值求解雷诺平均非定常N-S方程方法完成了单级跨音速压气机动静叶相干三维非定常粘性流动计算,给出了高分辨率NND格式的时间推进有限差分解。对某跨音速单级压气机进行了定常和非定常流的详细数值计算,设计转速下定常粘性流场计算所得到的性能参数与实验结果吻合较好。在此基础上,对最高效率状态进行了动/静叶相互干扰的非定常流动计算,并对结果进行了分析。      

大尖削机翼非定常跨音速绕流的数值解法(英文)

本文给出一种计算飞机机翼上定常和非定常跨音速气动力的数值解法,使用了一种特殊设计的坐标变换,是用时间精确交替方向隐式(ADI)有限差分算法来求解非定常跨音速修正三元小扰动位势方程。给出了F5战斗机机翼的数值结果并与XTRAN3S,ATRAN3S及试验结果进行了比较,表明本方法是有效的和经济的。     

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

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

基于Navier-Stokes方程/自由尾迹/全位势方程的旋翼流场模拟混合方法

为充分考虑旋翼尾迹对流场的影响和减少尾迹的数值耗散，建立了一个基于Navier-Stokes方程／自由尾迹分析／全位势方程的旋翼流场求解的新的混合方法。该方法的求解域由三部分组成：一是围绕旋翼桨叶周围的粘性区域．采用可压Navier-Stokes方程来捕捉近场信息，包括激波及尾迹；二是离桨叶较远、粘性可以忽略的等熵流区域，以全位势方程来描述其流动；三是在无粘区域中嵌入自由尾迹模型，模拟桨尖涡从粘性区域进入势流范围的发展变化。为便于流场分区求解和信息传递，采用了嵌套网格方法，并给出了不同区域之间的信息传递方法。以两叶的Caradorma＆Tung模型旋翼和四叶的UH-60A直升机旋翼为算例，计算给出了旋翼桨叶表面的压强分布以及桨尖涡的位置，并与可得到的试验数据及无尾迹模型方法的计算结果进行了对比，表明本文的混合方法能够很大程度地减少旋翼尾迹的数值耗散。     

计算三维振动机翼上跨-亚音速压力分布的局部线化面元法

本文以小扰动位势流理论为基础,采用局部线化假设,建立了三维振动机翼的跨音速下洗积分方程;给出了跨-亚音速流核函数的计算方法;用推广的偶极子栅格法计算了矩形谐振动机翼上的非定常压力分布。结果表明,本方法可以较准确地计算非线性的跨-亚音速非定常流动,关于跨-超音速流动的数值解法,尚需作进一步的研究。     

跨音速压气机动静叶排相干的三维非定常流动数值分析

采用时空二阶精度ＮＮＤ显式差分格式求解三维Ｅｕｌｅｒ方程，对跨音速压气机级中动静叶排相干形成的非定常流动进行了数值分析。除动静叶排流动的控制方程体系、数值方法外，着重讨论了动静叶排计算域、网格以及边界条件的处理。对某压气机第１级动静叶排相干的非定常流场的计算，给出了叶片表面的马赫数分布、压力波动幅值，以及叶片槽道中瞬态密度等值线等，并分析了结果的合理性。     

基于欧拉方程的机翼跨音速颤振的频域计算

机翼跨音速颤振的频域计算方法是以给定机翼模态分布下机翼上各点的模态值作为运动幅值，以三维非定常Euler方程为控制方程，采用有限体积法和双时间推进，求解三维机翼简谐运动下的非定常气动力。所求得的气动力作为已知值运用于颤振方程，利用v-g法进行求解。对得到的一系列的阻尼、速度和频率进行了线性插值，从而得到颤振速度和颤振频率。     

旋翼/机身气动干扰的数值模拟

通过求解粘性非定常可压缩N-S方程,发展了旋翼/机身气动干扰的数值模拟方法,主要技术措施包括:(1)通过ADT数据结构,建立了高效的多块对接网格动态重叠技术;(2)在可压缩方程求解中引入低速预处理方法,克服低速流场求解时的刚性问题;(3)加入了多重网格方法,提高收敛效率;(4)湍流模型方程求解需要反复计算壁面最近距离,采用Wigton的优化算法很大程度上提高了计算速度。利用建立的计算方法对佐治亚大学(GIT)旋翼/机身干扰模型进行了数值模拟,分别比较了机身表面的平均压力分布、旋翼的下洗流场以及非定常压力脉动,对干扰现象及其流动机理进行了分析。     

时间推进积分方程法求解非定常跨声速流动

本文提出了一种求解非定常跨声速流动的新方法——时间推进积分方程法,此法克服了时间线化积分方程法的限制,能较好地模拟激波的运动。本文首先用一维波(?)问题——模型问题阐明此法的基本思想,然后将它应用于二维非定常跨声速流动中。本文还首次引入拟速度位的概念,使时间推进积分方程式得到简化,尾涡条件和Kutta条件更易处理。数值计算表明时间推进积分方程法是合理可靠的。     

二维非定常跨音速时间线化积分方程的推导

1.引言 求解非定常跨音速流动的主要困难是非线性问题。对于微幅翼型振动问题(同时引起激波的微幅振动)可做时间线化简化处理,得到时间线化微分方程定解问题。时间线化积分方程是由时间线化微分方程推导出来的。     

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

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

低超声速来流中机翼跨声速非定常气动力数值解法

本文给出一种低超声速来流中机翼跨声速非定常流的计算方法。该方法使用跨声速小扰动理论,正确地计入了沿流向平面流动特征线的压力波的运动边界条件。流动方程是使用时间精确交替方向隐式(ADI)有限差分解法求解。试算了M=1.1的矩形翼和M=1.1、1.34的F5机翼的定常和非定常气动力,并与国外的风洞试验和计算结果作了比较,符合很好。表明本方法是有效的。     

计算绕薄翼型跨音速非定常流的积分方程法

 本文从小扰动方程出发推导出绕薄翼型的跨音速非定常流动的积分方程,为使此方程适用于具有激波的流场引入人工粘性,并对其数值求解。通过算例讨论了人工粘性和计算域大小对计算结果的影响,与实验及其它数值计算结果比较,表明本方法的准确度令人满意。计算量小,收敛性好是此法的特点。     

三角翼匀速上仰跨声速非定常流场数值模拟

用有限差分法求解Ｎ－Ｓ方程，数值模拟了细长三角翼作匀速上仰运动时的距声速绕流流场。给出了非定常流场结构，涡破裂特性和非定常升力，阻力特性，研究了上仰运动对动态失速的影响。     

任意回转面亚、跨音流场计算方法

本文所设计的任意回转面亚、跨音流场的计算软件 ,能根据叶型和进出口条件分析流动特性 ,并自动选择精确和节省机时的合适的数值解法。用户在给出叶型坐标和定解条件后 ,可以得到叶面马赫数分布或压力分布 ,全流场的等马赫数或等压图     

欧拉方程计算非定常流动的一种数值方法

提出了利用欧拉方程计算非定常流动的一种数值方法。在利用中心差分的有限体积法和双时间推进法求解非定常欧拉方程过程中，通过在固定物面边界上满足完全的非定常边界条件，计算出非定常流场。算例结果表明，本文方法具有高效、简便的特点，而且能描述激波的非定常运动和变化情况。     

跨声速压气机转子叶尖非定常流动的数值模拟

为了进一步加深对转子叶尖区域非定常流动现象的认识，针对某一跨声速单转子轴流压气机，采用三维非定常数值方法开展了详细研究。对单转子在不同进气条件以及不同工作流量下分别进行非定常模拟，来研究该转子叶尖区域复杂流动结构以及叶尖非定常流动的发展过程。结果表明:不同进气条件下，转子叶尖区域流场结构形式表现不同.当进口非轴向进气时，叶尖非定常流动呈单通道周期形式；而轴向进气时，叶尖区域出现了类似于“旋转不稳定”的沿周向传播的非定常流动现象，且随着工作流量的减小，叶尖区域流场结构也相应发生改变.