求解叶轮机械三维粘性可压流的初步研究

给出了用“参数多项式方法”求解叶轮机械三维粘性可压缩流的基本思路、方程和方法,采用圆柱坐标系下的三次参数多项式。中还应用此方法以一个透平和一个压气机叶片为例计算了三维粘性流场。     

解不可压N—S方程的PCIMG法

本文提出了一种解不可压N-S方程的预测、修正、迭代、多重网格法。这种算法程序简单,结果精确及收敛快速。     

物理质团法关于可压流计算的程序验证

为了评估无网格物理质团法(PECM)对可压流计算的可信性,结合物理认识对该方法的设计思想和数值格式进行了先验分析,基于9个一维模型数值解与参照解的对比信息,通过后验分析方法对PECM进行了程序验证,参照解包括无网格再生核粒子法(RKPM)的数值解以及精确解。验证结果表明:各种物理参数以及状态方程的强烈间断会导致RKPM严重的失稳或失真,而PECM在包括极大密度比在内的各种强间断条件下仍保持稳定和收敛并具有较高的准确度。      

二维可压无粘流的自适应流量修正有限元解

从非定常Ｅｕｌｅｒ方程组出发,采用非结构三角形网格,研究二维无粘流动问题的自适应流量修正有限元解。将ＦＥＭ－ＦＣＴ与网格自适应相结合计算了ＮＡＣＡ００１２翼型的三个算例。计算结果与实验及其他计算结果比较,亚音速和超音速绕流符合很好;跨音速绕流情况存在一些差别。     

逆风因子法在非正交曲线坐标系中解不可压N-S方程中的应用

所谓的逆风因子法就是分别用二阶精度公式计算出对流项和扩散项的贡献,然后按此两项贡献的比值来确定一个逆风因子,来保证对流扩散方程数值解的精度和稳定性。这种逆风因子法是作者在用有限元素法解不可压N-S方程中提出来的。本文将逆风因子法应用于在非正交曲线坐标系中用有限差分法解不可压N-S方程。      

可压流的龙格-库塔时间步进求解法及其应用

研究了用经典四阶龙格－库塔法计算可压缩流场的可行性,并将其用到了电弧喷射推力器内部等离子体流场的数值求解中。应用情况表明,在结合了局域时间步长、隐式残值光滑加速收敛措施后,本格式能够成功地计算比较复杂的可压缩气体流场以及等离子体流场。有关计算结果揭示了气体流经电弧喷射推力器通道但无电流时形成的纯气动流场以及有电流通过时形成的等离子体流场的丰富的结构和一些重要的影响因素,为研究其过程机制提供了依据。     

一种不需要迭代的发动机辅助变量建模方法

现代航空发动机建模计算时常用牛顿迭代法求解非线性方程组,但由于发动机模型非线性太强且十分复杂而不易在全包线收敛。本文针对这一不足之处提出了一种不需要迭代就可以求解此非线性方程组的方法:引入航空发动机容积动力学中的变量及方程作为辅助变量及辅助方程,从而使整个非线性方程组闭合而可以直接求解。本文将此方法成功应用于某涡扇发动机,取得了良好的效果,验证了本方法的可行性。      

一种较简单的低循环疲劳寿命预测法

<正> 1。引言 局部应力应变法是目前预测结构裂纹起始寿命的一种行之有效的方法。根据局部应力应变法的原理,可以用计算的方法来预测构件裂纹起始寿命。在计算中,从载荷谱得到局部应变谱,可利用有限元法分析结构的局部应力应变响应。当然还可利用半经验公式(如Neuber法等)。在计算中若考虑材料的记忆特性和循环σ-ε曲线,这无疑对真实地反映受载构件的局部应力应变响应是有利的。但要做到这一点,用有限元分析所需的机时相当可观,因为它要对整个寿命期间的载荷历程进行分析计算,至少也要针对一个典型谱分析计算。 计算中材料的循环应力应变曲线一般可采用下式     

可压缩粘性流绕翼型的N─S方程自适应网格解

从Ｎ－Ｓ方程出发,采用ＬＵ－ＡＤＩ隐式分解方法求解绕翼型的可压缩粘性流动。湍流模式采用Ｂａｌｄｗｉｎ－Ｌｏｍａｘ代数湍流模式和ｑ－ω二方程微分模式以研究湍流边界层的非平衡效应;为更好地捕捉激波,对网格进行了自适应调整。通过对ＮＡＣＡ００１２和ＲＡＥ２８２２两种翼型在亚音速和跨音速不同马赫数和雷诺数下进行大量计算,表明该方法对改进湍流计算、提高激波分辨率有较好的效果     

三维可压和不可压缩流动数值模拟的一种预处理多重网格方法研究(英文)

将矩阵预处理方法与多重网格方法结合,发展了适用于机翼及复杂外形气动数值优化设计的高效数值模拟方法和程序。在不改变定常流动解的前提下,对 Navier-Stokes 方程的时间导数项实施 Choi-Merkle 矩阵预处理,从而改善可压缩控制方程在低速情况下的系统刚性。Choi-Merkle 预处理最初用于二维低马赫数粘性流动数值模拟,本文将其推广应用到计算三维流动,并针对提高跨音速流动计算收敛性进行了相应的改进。预处理后控制方程采用中心格式有限体积进行空间离散,Runge-Kutta 混合多步方法进行显式时间推进求解,并应用 FAS 多重网格方法加速收敛。采用所发展的方法和程序数值模拟了绕 M6 机翼和 M100 机翼的低速和跨音速流动(马赫数从0.01 到 0.839)以及绕 F4 翼身组合体流动。计算结果与实验值吻合良好,验证了所发展的方法的正确性。算例研究表明所发展方法较大幅度地提高了三维流动数值模拟的效率,且同时适用于可压和不可压流动计算,在机翼及复杂外形气动优化设计方法具有应用前景     

非定常不可压粘性/无粘性耦合方程的一种分步分解方法

给出了数值求解初始变量不可压Navier-Stokes/Euler耦合方程的一种分步块LU分解方法.与传统的时间分裂法不同,该法无需压力中介边条件,从而避免了传统时间分裂法要求的复杂的压力中介边条件逼近.分步块LU分解方法可看做经典的Uzawa算法的改进,后者曾被成功应用于不可压Navier-Stokes/Euler耦合方程的求解.但本文显示分步块LU分解法比经典的Uzawa方法更经济.分析显示该法具有良好的稳定性和高精度,数值结果支持这一理论分析.     

一种简单的自适应延迟补偿方法

针对模型参考自适应控制对延迟鲁棒性较差的问题,提出了将动态特性接近被控对象的参考模型作为前馈补偿器,从而改进系统性能的补偿方法。仿真结果表明,该补偿方法效果明显,且对不确定的延迟扰动有一定的鲁棒性。     

一种基于三维粘性流场的乘波体生成方法

为了便于乘波体的优化设计工作,本文实现了一种基于任意三维粘性流场,在激波面上截取乘波体前缘线来生成乘波体外型的设计方法,并与NASA Langley研究中心的CFL3D软件的算例结果进行了比较检验,误差满足工程计算要求.由于在流场计算中直接计入了粘性,而不像大多数的研究方法那样通过无粘计算再进行粘性修正,而且对于高超声速乘波体设计,可以在程序中考虑真实气体效应、化学反应等因素,因此本文的设计方法更接近工程实际,更具有推广价值.     

绕流问题的区域收缩法及其数值解

本文提出了一种适于绕流问题的数值方法,它考虑了流场中不同区域内不同的物理现象,给出了外部边界条件的良好逼近。使用这种方法可以节省大量的计算机存贮单元和计算时间,也可以推广到三维绕流问题的数值计算。     

盘形凸轮廓线设计的一种简单算法

在分析直动从动件盘形凸轮机构廓线方程的基础上,提出了一种计算过程较为简单的算法,并应用MATLAB仿真软件通过实例对该算法进行了验证.     

一种基于TDDM裂谱信号的解模糊位搜索法

时分数据调制(TDDM)裂谱信号作为一种新调制信号,具有独有的电文翻转位模糊特性,导致了采用通常的搜索方法无法正确输出电文的现象。针对TDDM裂谱信号特性以及现有搜索方法的不足,提出了一种解模糊位搜索法。该方法主要分为合路快捕及解模糊位判决两个阶段。其合路快捕阶段,不但可以快速缩小搜索范围,而且可以避免由扩频方式而引起的漏捕现象;解模糊位判决阶段,可以去除电文翻转位的模糊度,正确输出电文的翻转位置。在理论推导的基础上,从性能及信噪比容限角度进行了仿真测试,测试结果表明,解模糊位搜索法对TDDM裂谱信号的电文翻转位判决精度可以达到100%。