Roe格式中不同类型熵修正性能分析

针对Roe格式中较为流行的3类熵修正:Muller型、Harten-Yee型和Harten-Hyman型熵修正,从理论分析入手, 辅以Euler方程为控制方程的激波管问题,前台阶流动和运动激波的双马赫反射3个数值实验,对不同熵修正的性能做了深入研究,得出如下结论: Muller型和Harten-Yee型熵修正方法作用于激波和非物理的膨胀激波2种情况:激波情况下引入的数值耗散,能根本改善"Carbuncle"现象,膨胀激波时其数值粘性也足以使膨胀激波得以耗散;而类Harten-Hyman型熵修正只对膨胀过程起修正作用,对激波情况无效,不能改善"Carbuncle"现象,该类熵修正的数值粘性较小,不足以使膨胀过程正确求解;直接采用 δ 值代替特征值的特征值修正方法效果好于传统的特征值修正方法.     

基于DSI插值的三角网格质量优化

通过对三角网格的单元顶点进行几何位置调整,提高了网格的质量,实现了网格的质量优化.几何位置调整是使用离散点光滑插值(DSI,Discrete Smooth Interpolation)实现的,针对在计算时影响质量优化的邻接边界的单元顶点,采用了在边界处补偿三角形的方法,消除了单元收缩,提高了网格的质量.与加权拉普拉斯算法进行了比较和分析,优于拉普拉斯算法;为了使三角网格在位置调整时保持原始网格的几何细节特征,在插值算法中施加了控制点约束.最后使用算例对算法进行了验证.      

翼型动态失速气动力二次峰值数值模拟研究

为揭示翼型动态失速状态下气动力二次峰值的发生机理,基于运动嵌套网格技术、有限体积方法、LU-SGS隐式格式和Roe-MUSCL格式建立了俯仰振荡翼型非定常流场的数值模拟方法。首先,基于所建立的数值方法对NACA0012翼型在深度动态失速状态下的气动特性进行模拟,计算结果与试验数据吻合良好,验证了数值模拟方法的准确性。然后,通过对NACA0012翼型动态失速状态流场的研究,揭示了气动力二次峰值的发生机理。最后,开展了翼型厚度、弯度和弯度位置等外形参数对气动力二次峰值的影响研究。结果表明,动态失速涡诱导形成的后缘涡是导致气动力二次峰值的关键因素;翼型外形参数的变化会引起动态失速过程中动态失速涡和后缘涡的变化,使得气动力二次峰值相对谷值的增量有规律地增加或减小,二次峰值位置有规律地前移或后移。     

激光功率与底面状态对选区激光熔化球化的影响

为研究激光功率与底面状态对选区激光熔化熔池流动的影响,基于离散单元法建立了选区激光熔化铺粉模型,采用粒径分布与实验相符的马氏体时效钢粉末分别铺展到平坦底面和增材底面上,将计算获得的粉末分布导入到基于有限体积法建立的选区激光熔化熔池计算流体力学模型中,研究激光功率和基板底面粗糙度对熔池流动和熔道表面形貌的影响。采用激光单道扫描实验验证铺粉模型和选区激光熔化模型。结果表明：随着激光功率的降低,单位长度的球化数量增加;由于增材底面使熔池润湿性变差,同时又对熔池流动行为产生扰动,使得增材粗糙底面上熔道的球化数量增加。选区激光熔化铺粉模拟及激光单道扫描模拟结果与实验结果吻合较好。本研究可为选区激光熔化工艺中工艺参数的选择提供理论指导。     

二元曲面可调进气道流量系数精确预测方法

为了满足二元曲面可调进气道模态转换马赫数范围(来流马赫数为2.2~3.2)的流量要求,针对唇口平移、转动和转动+平移三种调节方案,基于理论分析和基准进气道的流场,提出了一种流量系数精确预测方法,并通过数值计算进行验证且获得了进气道的总体性能。结果表明:调节后的进气道流量系数与预测值完全相等,而且无需多次试算,符合设计预期,可拓展应用于轴对称进气道。相对基准进气道,唇口前移时流量系数和压缩效率同时增加,来流马赫数为2.5时出口总压恢复系数相等而增压比增加了14.6%;在降低相同流量系数条件下,后移唇口使得增压比和压缩效率均降低,来流马赫数为2.5时出口总压恢复系数基本相等而增压比减小了12.9%,转动唇口使增压比进一步减小了9.1%,唇口后移方案性能更优。      

变几何径向透平损失模型及其分析

基于99 mm叶轮实验, 分析了变几何径向透平静叶的各种损失以及不同的损失模型.给出了变几何透平静叶冲角损失的定量分析.通过计算对比不同静叶喉部宽度下叶道边界层摩擦损失, 验证了Glass-man模型在计算叶片流道损失时, 与试验结果是相符的.基于NASA损失模型分析了径向透平静叶的尾缘损失, 损失随着出口边厚度呈线性增加, 减小出口边厚度是减小尾缘损失的有效方法.      

高超声速飞行器-进气道一体化热流数值计算

采用CFD(计算流体动力学)技术, 开展了飞行器前体/发动机一体化气动热环境分析.对层流区、转捩区和湍流区分别采用计算模型, 在湍流区利用压缩性修正的SSGZ-Jk-ε湍流模型, 在转捩区引入代数型转捩因子模型描述边界层由层流逐渐过渡为完全湍流的流动过程.计算了前体和内通道的表面热流, 并与实验结果进行了对比.结果表明所采用的计算方法可以较好地预测前体及发动机内通道热流率, 流动状态、几何结构及激波入射对热流值影响较大.      

SIMULATION OF COMPOSITE NON-LINEAR MECHANICAL BEHAVIOR OF CMCS BY FEM-BASED MULTI- SCALE APPROACH

The non-linear behavior of continuous fiber reinforced C/SiC ceramic matrix composites （CMCs） under tensile loading is modeled by three-dimensional representative volume element （RVE） models of the composite. The theoretical background of the multi-scale approach solved by the finite element method （FEM） is recalled first- ly. Then the geometric characters of three kinds of damage mechanisms, i.e. micro matrix cracks, fiber/matrix interface debonding and fiber fracture, are studied. Three kinds of RVE are proposed to model the microstructure of C/SiC with above damage mechanisms respectively. The matrix cracking is modeled by critical matrix strain en- ergy （CMSE） principle while a maximum shear stress criterion is used for modeling fiber/matrix interface debond- ing. The behavior of fiber fracture is modeled by the famous Weibull statistic theory. A numerical example of con- tinuous fiber reinforced C/SiC composite under tensile loading is performed. The results show that the stress/ strain curve predicted by the developed model agrees with experimental data.     

一种水平尾翼流动控制装置的实验研究

某型机采用上单翼低平尾布局,在着陆襟翼小迎角状态时平尾下翼面翼根部位出现局部气流分离,导致飞机振动,力矩特性出现异常变化。本文提出的解决方案是在平尾翼根前方0.12倍根弦长,下方0.30倍根弦长位置的机身上加装一对小展弦比负弯度小翼作为涡流发生器/导流片,一方面加速了后方分离区边界层与外流的能量交换,另一方面利用其上洗作用降低了平尾翼根区域的局部负迎角绝对值。通过数值计算和风洞实验研究表明,优化后的导流片使平尾分离区面积缩小50%,小迎角俯仰力矩拐点推迟4°以上,以较小的性能和结构重量代价解决了局部气流分离问题,拓展了飞机飞行边界。     

汽心泵特性分析与计算方法

本文给出汽心泵压头特性、汽心变化、泵内水力损失、集流器内压力分布和径向力等主要特性的计算方法和程序框图。该算法可用来估算新泵设计初始阶段的特性,亦可用来对现有泵的特性作近似验算。