涡格法中Multhopp离散的超收敛性的进一步研究

在中弧线为四阶以下多项式的二维薄翼绕流问题的研究中，Ｍｕｌｔｈｏｐｐ离散超收敛性的证明展示了涡格法中优化的离散化格式的潜力。本文把这一研究扩展到中弧线为任意Ｎ阶多项式的二维薄翼绕流的一般情况。通过误差分析，证明了Ｍｕｌｔｈｏｐｐ离散对该问题可以得到离散化误差为零的翼面涡密度、升力和俯仰力矩，唯一的要求是使用的离散元数目应大于［（Ｎ＋１）／２］。     

模拟有限长涡片旋转的一种新数值模型

本文把涡格法中的超收敛性推广到自由涡片的演化问题，涡片的演化由ＢｉｒＫ－ｈｏｆｆ－Ｒｏｔｔ方程所控制，该方程具有强烈的非线性，难以用解析法求解。在数值模拟中通常离散涡方法。本研究表明对自由涡片的旋转问题在精确解和离散涡方法结果之间的误差主要来自点涡近似。使用基于超收格式的新数值模型可以显著地提高数值解的精度。     

升力面尾涡卷起的地面效应

本文就二维平板从理论上证明了改进涡格法的超收敛性。用该方法研究了近地升力面尾涡卷起的地面效应。尾涡面上的自由线涡方向应与当地流线方向一致,通过迭代满足这一非线性边界条件的确定尾涡的形状。地面效应以镜象涡系模拟。一系列计算结果表明迭代迅速收敛。本文讨论了升力面距地高度对尾涡位置和形状的影响,计算表明当该高度小于一定距离后,尾涡呈现不稳定。     

超声速空气介质自由扩散混合流涡的演化

在不同谐波扰动促发方式下,观察二维超声速空气/空气自由扩散混合流中涡的形成及演化.考虑含空气组分扩散效应的可压缩Navier-Stokes方程,对流项采用3阶迎风紧致格式离散,输运项采用6阶对称紧致格式离散,非定常时间推进采用3阶紧致存储显式Runge-Kutta方法.在无相差情形下,获得了大尺度基频涡结构的饱和、一次对并、二次对并、三涡对并等现象.在Ma_c=0.3低对流Mach数下,基频涡较饱满,但流向尺度较小.受空气真实气体特性影响,在Ma_c=0.8高对流Mach数下未发现涡/小激波结构.      

突扩后突缩的二维流场计算

本文模拟固体火箭发动机药柱后端流场,建立了简化的计算模型,采用压力耦合方程半隐式解法,对突扩后突缩的二维粘性流场进行了计算分析.在计算过程中,对网格划分及控制体的选择、复杂边界条件的处理、差分格式的选取、方程的离散化及方程的解法等做了分析研究.计算结果表明:在突扩台阶后得到的二次涡与Finash,F.和Ghaddar,N.K.等人的结论吻合,进一步证实了Pan,F.和Acriuos,A.的嵌套涡族理论.     

带运动顶盖的封闭方腔内湍流流场的计算

本文运用κ-ε模型求解了带运动顶盖的封闭方腔内的二维湍流流场,其流动雷诺数为10000,对乘方定律格式(PLDS)和二次上风差分格式(QUICK)进行了比较,采用了PISO求解方法、截面连续性修正并对k、ε方程的源项作特殊处理,使计算的收敛性和稳定性明显提高。计算出的主涡、二次涡位置及流场速度分布值与实验吻合良好。计算所用网格为12×12,在Honeywell DPS 8/49计算机上运行的CPU时间对PLDS与QUICK分别为1.9分和1.8分。     

有翼飞行器大攻角分离涡模拟

基于势流方程,采用涡格法和三维离散涡相结合,模拟具有升力面分离涡的大攻角弹冀、弹体、尾翼组合体绕流。该方法适用于高亚音速、小攻角或大攻角有翼飞行器的压力分布、法向力和压力中心等气动特性的数值计算,并具有显示三维涡流场的能力。     

ENO-AUSMPW格式对激波涡干扰的数值模拟

高阶ENO格式在复杂流场计算中具有重要的地位,它理论上可以达到任何阶的精度。本文采用二维高阶ENO插值方法,构造了高阶AUSMPW格式,并对包含复杂激波结构的激波涡干扰流场进行了数值模拟。计算结果表明了高阶ENO插值下的AUSMPW格式具有较高的激波分辨率和较低的数值耗散。     

端板法翼型实验的数据处理与洞壁干扰修正方法

 用端板法测量气动力数据时,端板附近有自由涡拖出。通过实验观察假设了端板翼型实验的涡系模型。由涡系模型可计算出机翼加端板后的有效展弦比,再按举力线理论可将实验结果转换成二维(翼型)数据;按镜象法也可作升力洞壁干扰修正。对于端板法中间剖面的测压结果亦应作二维转换和洞壁干扰修正,方法类似。经过这样的数据处理后,翼型数据准确可靠。     

叶栅内冷气射流场结构的数值模拟

采用具有 TVD性质的三阶精度 Godunov格式对涡轮叶栅背弧表面进行了冷气喷射的数值模拟。射流由冷气孔喷出后在空间的发展过程中逐渐形成典型的卵形涡对, 在压力场的作用下射流被压弯并逐渐向背弧表面靠拢。射流在出口附近具有刚体的特性, 并在射流前方形成了鞍点及马蹄涡。同时本文还验证了大冷气喷射速率条件下, 卵形涡消失于主流而非边界层内的结论。      

一种求解涡轮平面叶栅的跨音松驰法

一、引　　言跨音松驰法是七十年代初发展的一种新方法。文献 [1 ]是一篇关于跨音松驰法的奠基性文章 ,文章提出了求解小扰动势方程的混合差分格式。文献 [2 ]发展了求解全位势方程的旋转混合差分格式。近年来 ,求解跨音全位势方程的方法有了更全面的发展 [4、5]。作者于 80年发展了一种求解大弯度二维叶栅的跨音松驰法[3] 。本文将该法推广到求解涡轮平面叶栅跨音带激波流场 ,直接求解全位势方程 ,进一步提高边界区差分格式的精度 ,并注意初场的给法。本法具有节约机时和内存的优点。数值算例表明 ,对于头部不大的涡轮平面叶栅 ,应用本方法…     

Integral-projection method and singularities of its application in mating shells with inconsistent parameters of approximation

We describe a numerical integral-projection method used by the authors for the approximate solution of systems of interrelated two-dimensional linear boundary-value problems in mechanics of composite shell systems. The method is based on discretization in each shell substructure of a two-dimensional problem along one of coordinates using a projection-grid variant of the Galerkin-Petrov method and its subsequent transformation to a system of ordinary differential equations; by integration and introduction of sought functions as unknown derivatives, the system is reduced to a system of integral equations being solved by the method of mechanical quadratures. The method is characterized by the fact that its application requires no additional conditions of conformity with discretization parameters of substructures being mated.     

Real-time solution of nonlinear potential flow equations for lifting rotors

Analysis of rotorcraft dynamics requires solution of the rotor induced flow field.Often,the appropriate model to be used for induced flow is nonlinear potential flow theory(which is the basis of vortex-lattice methods).These nonlinear potential flow equations sometimes must be solved in real time––such as for real-time flight simulation,when observers are needed for controllers,or in preliminary design computations.In this paper,the major effects of nonlinearities on induced flow are studied for lifting rotors in low-speed flight and hover.The approach is to use a nonlinear statespace model of the induced flow based on a Galerkin treatment of the potential flow equations.     

锯齿形转捩片触发高超声速进气道边界层转捩的大涡模拟

为促进锯齿形转捩片在高超声速进气道中的应用，以地面风洞条件下的二元进气道为研究对象，采用高精度大涡模拟方法对锯齿形转捩片在三级压缩楔面上触发的边界层转捩现象开展了研究。数值方法基于隐式亚格子模型，空间离散采用高精度通量限制型紧致格式，时间推进采用显式Runge-Kutta方法。数值模拟清晰捕捉到了边界层转捩的空间发展演化过程，并获得了统计平均流场以及流场脉动特征。数值模拟结果表明转捩片能够有效触发进气道压缩面边界层转捩；通过与等熵压缩面及单楔面数值模拟结果的对比分析，获得了转捩片触发边界层转捩的内在机理，为后续研究工作奠定了基础。     

高亚声速二维风洞收缩段的设计

本文提供了一个仿射变换式,它将亚声速位势方程精确变换成Laplace方程。应用这个理论,本文给出了一个设计二维高亚声速风洞收缩段的精确方法。在马赫数为0.9时,求得的收缩段形状曲线与由计算机提供的位势方程数值解相比较,二者极为符合。此外本方法计算简单,精度高,并具有广阔的应用范围。     

跨音速大迎角Euler方程数值分析

应用有限体积法离散三维Euler方程,选用三步显式格式进行时间推进求解;并通过当地时间步长、残值光顺、焓修正等方法加速收敛。对一维波动方程的von Neumann稳定性分析说明这种三步格式的最大Courant数是2。对三角翼跨音速大迎角流动的数值分析表明文中的三步格式及焓修正方法优于人们经常采用的四步、五步格式以及焓阻尼技术。     

轴向超声速通流串列构型弦长比对气动性能影响的数值研究

为避免高马赫数、大攻角来流引发的叶片颤振,将串列叶片技术引入到超声速通流风扇叶栅中,对其进行串列改型及气动性能研究。利用准二维数值模拟,对串列叶片前、后排叶片的弦长比参数进行了详细的对比研究。结果表明：影响气动性能的关键因素是后排叶片进口压力侧激波的落点,在本文研究条件下,随着弦长比的减小总压损失呈减小的趋势,当弦长比由0.99减小到0.43时,设计攻角下,15°折转角叶型总压损失可减小27%,30°折转角叶型总压损失可减小38%。进一步的研究表明,通过减小弦长比可有效控制后排叶片前缘斜激波在相邻叶片吸力侧的落点以实现损失降低,并且这种降低效应在小弯角叶型上比大弯角叶型更容易实现。     

用于气动声学计算的非均匀网格紧致差分格式

为克服传统紧致差分格式在数值求解非均匀网格问题时产生的寄生波,构造了一种新的高精度紧致差分格式.通过泰勒展开分析方法,详细给出了格式系数的通用形式;利用傅里叶分析方法,分析了数值耗散、色散误差.以3对角6阶精度紧致差分格式求解均匀扰动网格问题为例,计算表明:色散值和耗散值随扰动因子的增加而更加趋近于精确值;当扰动因子大于0.213时,格式不稳定,当扰动因子小于等于0.213时,格式渐近稳定;对一维对流波和二维波传播的模拟计算所得数值解与精确解吻合,体现了该格式在求解非均匀网格问题时的优越性.