基于全欧拉方程的二维平行剪切层声波产生和辐射的数值模拟

本文基于全欧拉方程对二维平行剪切层声波的产生和辐射进行了计算气动声学数值模拟。空间离散采用了频散相关保持有限差分格式,时间积分采用了低频散低耗散的龙格库塔法。数值模拟结果表明,当声源强度较小时,基于全欧拉方程的计算结果与解析解及基于线化欧拉方程的计算结果均符合得很好,当声源强度较大时(约140dB),线化欧拉与全欧拉的计算结果有较大的偏差,这说明此时非线性的影响不可忽略,线化方程不再适用。      

轴向磁通磁场调制型电机的解析分析

提出了一种新型二维解析分析方法用于计算轴向磁通磁场调制型电机的气隙磁场分布和电机特性。对电机每个子域列出拉普拉斯方程或泊松方程并结合边界条件,求解出每个子域的矢量磁位,进而得到气隙磁场分布、空载反电势、电磁转矩和轴向磁拉力的解析表达式。为了评估解析方法的计算精度,对二维解析分析结果和三维有限元仿真结果进行了比较,结果表明二维解析计算结果的计算误差在10.3%以内。由于解析分析在保证计算精度的前提下大大缩短了计算时间,因此该方法适用于轴向磁通磁场调制型电机初始设计阶段电机参数的确定。     

计算带有层流分离的粘性/无粘流相互作用的一个方法

基于粘性/无粘流相互作用模型和边界层积分方程,推导出计算二维和三维不可压层流粘性/无粘流相互作用的方程组。该方程组可以用于带有中等程度的分离和再附情况下。文中给出了计算方法和算例。所得的二维计算结果与其他计算和实验结果相符。本法具有简便迅速的优点。     

预估蒸汽轮机叶片颤振的三维方法

发展了一种可考虑叶片间相角变化的三维失速颤振预测方法。二维计算的最小气动阻尼比三维值低, 说明取特征截面的二维方法偏于保守, 且特征截面的取法具有任意性, 故二维方法仅能作为叶片颤振初步预估之用。为提高颤振预估的精度, 需用三维计算对二维计算结果进行校核。三维和二维计算均表明叶片间相角对蒸汽轮机叶片气弹稳定性影响很大, 叶片颤振预测必须考虑这一关键因素的影响。      

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

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

航空发动机整机二维气动热力数值模拟

研究了用于航空发动机整机数值仿真的含粘性力项的二维欧拉流动模型,进行了整机流场数值仿真。用代数和椭圆形偏微分方程混合方法分别生成叶栅通道和燃烧室内贴体坐标网格;分别针对压气机转子、静子和涡轮叶栅,定义了三种不同损失系数;采用高阶精度Godunov格式和时间推进法,交替使用显/隐格式,收敛速度快,计算准确度高。计算了某双轴涡扇发动机设计点的整机性能与部件性能,计算结果与设计数据相吻合,分析整机仿真计算得出的流场参数分布,可以更好的了解流场细节。数值试验表明,建立的程序可以预测发动机稳态条件下的整机性能以及用于发动机的设计与改进,奠定了二维数值试验台建设的基础。     

基于HLLC格式的理想MHD方程数值模拟

针对理想MHD(magnetohydrodynamics)方程数值求解困难的问题,基于原始的HLLC (Harten-Lax-Van Leer Contact wave)近似黎曼解方法,发展出一种新的适用于理想MHD问题的通量计算格式。控制方程采用有限体积法离散,时间推进采用隐式的LU-SGS（lower-upper symmetric Gauss-Seidel）格式,并且采用双曲型散度清除技术来抑制磁场散度的累积。通过一维激波管问题的数值模拟表明,HLLC-MHD格式能准确的分辨并捕捉复杂的磁流体力学波系结构,保证数值计算的精度,最大的数值计算误差不超过10%;通过二维的Rotor问题的数值模拟研究表明,HLLC-MHD格式能够应用于多维理想MHD问题的数值模拟,并且能够准确捕捉多维情况下磁流场中的阿尔文波;通过对比有、无散度清除的Rotor问题计算结果表明,双曲型散度清除技术可以将磁场散度峰值从50降低到2,有效抑制高磁场散度区域的散度累积,同时也会将误差传播到低散度区域,并且引起边界处散度的累积,影响计算的稳定性。      

上海标二厂的计算法也适用于MJ螺纹

上海标准件二厂（简称上海标二）的螺纹毛坯直径计算方法（详见贵刊１９９７年第１期《辗压螺纹毛坯直径计算》一文）多年来一直用于计算所有对称圆柱普通螺纹，新近笔者对其是否能用于计算ＭＪ螺纹毛坯直径作了验证，得出的结论是：上海标二的计算式也适用于计算ＭＪ螺纹...     

翼上螺旋桨构型耦合气动特性及翼型优化设计

针对半环形式翼上螺旋桨构型，研究了螺旋桨-机翼耦合流场特性，并以短距起降（STOL）状态最优升阻特性为目标对机翼翼型进行全局优化。首先，针对螺旋桨-气动面耦合构型，通过动量源法与真实桨叶模型CFD的计算对比，分析动量源法用于该构型设计分析的可行性。其次，为得到有利于桨-翼耦合特征的新翼型，建立了翼上螺旋桨构型自由型面变形（FFD）参数化模型，采用遗传算法对翼上螺旋桨构型机翼翼型进行全局寻优设计，分析了优化翼型参数及流场变化规律。最后，将优化翼型用于三维半环形机翼，分析其流场特性与二维计算结果的异同，验证二维翼型优化的有效性。结果表明：真实桨叶多重参考系（MRF）方法不能准确计算翼上螺旋桨构型下的流场结构，而动量源法计算结果与真实桨叶滑移网格非定常方法较为吻合；采用二维动量源CFD方法进行翼型的遗传算法优化是有效的，受半涵道的保护，二维优化翼型的优势在三维构型中得到了有效继承；翼上螺旋桨构型的翼型优化应当着重关注翼面曲率变化，在本文计算状态下，通过增加桨盘附近翼面曲率、保持附着流动来加强Coanda效应，有效实现了气动增升，优化后机翼升力提高了22.51%，显著减弱桨盘后高压区并产生二次吸力峰值，同时保持了机翼负阻力特性。     

用二方程湍流模型对湍流射流场的数值计算

在许多工程问题中,利用Prandtl混合长度模型计算复杂类型的湍流射流问题已经不能满足要求。本文利用二方程湍流模型,以Raw Sinha等提出的计算射流及尾迹流平均流场的隐式有限差分方法为基础,建立了可用于二维及轴对称可压缩湍流射流场的计算方法。最后,以两个不同类型的同轴湍流射流场作为算例进行了流场计算,其结果与实验符合良好。     

亚、超音喷管中声传播的数值模拟

采用计算气动声学方法研究了管道近音速区的声传播和声波与激波的相互作用问题。空间离散采用高阶DRP格式,时间离散采用低耗散低频散的龙格库塔方法,边界条件采用Giles的特征变量一维非定常无反射边界条件。计算结果和解析解的精确吻合表明,计算中应用的这种高阶的DRP格式有很高的分辨率和很好的稳定性,适合声传播的计算。      

Interface flux reconstruction method based on optimized weight essentially non-oscillatory scheme

Aimed at the computational aeroacoustics multi-scale problem of complex configurations discretized with multi-size mesh, the flux reconstruction method based on modified Weight Essentially Non-Oscillatory (WENO) scheme is proposed at the interfaces of multi-block grids. With the idea of Dispersion-Relation-Preserving (DRP) scheme, different weight coefficients are obtained by optimization, so that it is in WENO schemes with various characteristics of dispersion and dissipation. On the basis, hybrid flux vector splitting method is utilized to intelligently judge the amplitude of the gap between grid interfaces. After the simulation and analysis of 1D convection equation with different initial conditions, modified WENO scheme is proved to be able to independently distinguish the gap amplitude and generate corresponding dissipation according to the grid resolution. Using the idea of flux reconstruction at grid interfaces, modified WENO scheme with increasing dissipation is applied at grid points, while DRP scheme with low dispersion and dissipation is applied at the inner part of grids. Moreover, Gauss impulse spread and periodic point sound source flow among three cylinders with multi-scale grids are carried out. The results show that the flux reconstruction method at grid interfaces is capable of dealing with Computational AeroAcoustics (CAA) multi-scale problems.     

高精度有限差分格式的色散优化及耗散控制

本文综述了我们在高精度有限差分格式的色散优化和耗散控制方面的研究进展。首先,我们提出了半离散有限差分格式色散和耗散相互独立的充分条件,实现了差分格式色散和耗散特性的独立调节。在此基础上,提出了色散最小、耗散可控的高精度差分格式,称为MDCD格式。MDCD格式已经得到了广泛应用,取得了很好的计算效果,但其主要缺点是耗散的调节依赖于经验。为了解决这一问题,我们进一步提出耗散的自适应调节方法。具有自适应耗散特性的高精度有限差分格式的基本特征是,差分格式的耗散能够随解的局部尺度自适应调节。为了构造这类格式,我们提出了一种新型的尺度识别器,它能够以等效无量纲波数的形式来定量衡量数值解的局部长度尺度。在此基础上,设计差分格式耗散参数与尺度识别器得到的等效无量纲波数之间的关系,从而构造了一类色散最小、耗散自适应的差分格式,称为MDAD格式。为了计算含有间断的问题,同时保持在光滑区的良好耗散特性,我们利用尺度识别器对一种经典的激波探测器进行改进,提出了一种新的激波探测器,并将自适应耗散格式与对应的WENO格式相混合,得到自适应耗散混合格式。近似色散关系显示该混合格式兼具高分辨率和鲁棒性。多个含间断流场的标准算例测试结果显示,自适应耗散混合格式具有良好的分辨率和激波捕捉能力。     

Application of high-order accurate methods to solve the computational gas dynamics and aeroacoustics problems

The application of the low dissipative high-order accurate scheme for numerical solution of the Euler and Navier-Stokes equations in the two-dimensional statement is considered. The scheme of calculating space derivatives is the seven-point central-difference approximation of the fourth order with the coefficients optimized to reduce the dispersion errors. The optimized six-step Runge-Kutta method is used for evaluating the time derivatives. By means of filtration, spurious pulsations are suppressed and shocks are processed. The results of verifying the calculation scheme realized on the stationary problem of flow around the model turbine blades are presented.     

气动声学直接数值模拟的空间差分格式分析

流场内声波的直接数值模拟对数值格式的耗散、色散特性提出了严格要求,首先介绍了两种适于声学模拟空间差分格式——紧致差分格式与色散关系保持(DRP)格式的基本思想,然后分析了格式的色散及其与四阶龙格-库塔法结合时的耗散特性。根据Fourier误差分析和数值模拟发现:无论从格式可有效求解的波数范围还是从与龙格-库塔时间推进结合后的耗散特性来看六阶紧致格式的可有效求解范围要大于DRP格式;不过,数值计算表明紧致格式的计算量远远超过DRP格式,对计算机资源提出了较高的要求。      

矩形/非结构混合网格技术及在二维/三维复杂无粘流场数值模拟中的应用

建立了一套模拟复杂无粘流场的矩形／非结构混合网格技术,其中非结构网格仅限于物面附近,发挥非结构网格的几何灵活性,以少量的网格模拟复杂外形;同时在以外的区域采用矩形结构网格,发挥矩形网格计算简单快速的优势,有效地克服全非结构网格计算方法需要较大内存量和较长ＣＰＵ时间的不足。混合网格系统由修正的四分树（２Ｄ）／八分树（３Ｄ）方法生成。本文将ＮＮＤ有限差分与ＮＮＤ有限体积格式有机地融合于混合网格计算,消除了全矩形网格模拟曲面边界的台阶效应,同时保证了网格间的通量守恒。该算法被成功地运用于ＮＡＣＡ００１２翼型的跨、低超声速绕流,二维弯管内运动激波的非定常流,以及三维捆绑火箭的超声速绕流的数值模拟。计算结果表明本方法在模拟复杂无粘流场方面的灵活性和高效性。     

正反顺序加载下安装边螺栓预紧力分散度的数值模拟与试验验证

在航空发动机机匣安装边装配过程中,不同位置螺栓之间会产生弹性交互作用,不同加载方案会使装配后的螺栓预紧 力存在一定的分散度。针对这一问题,建立了3层带止口机匣安装边3维接触有限元模型,研究不同加载方案下螺栓预紧力的施 加方法。基于有限元模型数值模拟了3层机匣安装边分别在顺序加载方案、交叉加载方案和正反顺序2轮加载方案下螺栓预紧力 的分散度。搭建了相应的试验装置,分别对3种加载方案的预紧力分散度进行试验验证。结果表明：采用正反顺序2轮加载后的 螺栓预紧力分散度最小,在3层止口机匣安装边上采用正反顺序2轮加载方案后的螺栓预紧力分散度较小,在有限元仿真中螺栓 预紧力分散度为2.43%,在试验中螺栓预紧力分散度为4.98%。采用此方案加载后,止口机匣安装边能达到更好的密封效果和较 好的刚度对称性。     

非结构混合网格高雷诺数流动计算中一种改进的LU-SGS隐式格式(英文)

The lower-upper symmetric Gauss-Seidel (LU-SGS) implicit relaxation has been widely used because it has the merits of less dependency on grid topology,low numerical complexity and modest memory require...     

非结构网格高超声速化学非平衡MHD流场数值模拟

 在非结构网格上对考虑化学非平衡效应的二维高超声速磁流体绕钝头体流动进行了数值模拟。控制方程由二维理想磁流体动力学（MHD）方程和组元连续方程两部分组成,化学动力学模型为5组元17反应模型。MHD方程空间离散采用AUSM格式,时间推进采用显式5步龙格-库塔格式,并通过弱耦合的方式与化学反应控制方程结合在一起。计算模型为二维钝头体,外加磁场为偶极子场,磁场源位于钝头体内部。在高超声速来流条件下,对有、无磁场干扰,是否考虑化学反应下的4种工况进行了数值计算,得到了满意的结果,并与有限的参考文献进行了对比。结果表明本文发展的方法能准确地模拟考虑化学非平衡效应的高超声速MHD流场。