AUSM类格式在热流数值模拟中的比较研究

系统分析了AUSM类空间离散格式中AUSM+、AUSMDV及AUSMPW三种格式的构造方法及特点,并分别对钝锥、钝双锥的高超粘性绕流流场进行了数值模拟,重点分析了格式对热流数值模拟结果的影响.数值结果表明:AUSM类格式在高超音速流动模拟中具有高间断分辨率和数值稳定性,其中AUSM+和AUSMPW格式更适用于热流数值模拟,而AUSMDV格式可能会出现“红玉”现象.     

跨声速压气机数值模拟格式效应研究

针对跨声速叶轮机械复杂内部粘性流动,开发了全三维高精度RANS数值求解器。分别采用FVS格式与AUSM+格式结合多种限制器对Rotor 67跨声速压气机转子进行数值模拟,分析和比较了各格式的计算效果并与试验结果对比。结果表明:与FVS格式相比,AUSM+格式的数值粘性更小,边界层的模拟精度更高;Hemker限制器的综合表现最优,Van Albada限制器的粘性分辨率略低于Hemker限制器,Minmod限制器对流动分离现象的捕捉能力较差,Van Leer限制器容易引入色散误差。     

高分辨率格式与钝舵绕流数值模拟

 采用Beam-Warming/TVD格式，Roe二阶、三阶以及四阶通量差分分裂格式数值计算雷诺平均Navier-Stokes方程，湍流附加粘性采用Baldwin-Lomax模型计算，针对超声速绕钝圆舵流动进行了数值模拟比较研究，计算给出了流场的详细结构，物面压力分布与实验值符合良好，尤其是Roe三阶格式给出了较好的计算结果，两个超声速回流区以及湍流边界层分离也得到了很好的数值模拟。     

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

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

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

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

二元喷管三维粘性流场的数值模拟

本文用数值方法对二元喷管的三维粘性流场进行模拟,将 Mac Cormack有限差分格式拓广到全三维 N-S方程的求解。在数值迭代中使用加快收敛速度的当地时间步长,并根据格式压缩计算机存储量和计算量。进行了两个算例的数值模拟,计算结果与实验相符。      

高阶气体动理学格式在湍流数值模拟中的应用

本文回顾了高阶气体动理学格式在湍流数值模拟中的应用。与传统的Riemann求解器相比,气体动理学格式可以提供时空耦合的演化过程,这对发展高精度格式十分重要。因此,基于两步四阶时间离散和高精度WENO重构,发展了具有四阶时间精度的气体动理学格式。该格式有更高的数值精度和稳定性,并且具有更好的处理复杂流动问题的能力。目前,两步四阶格式已经成功地应用到低雷诺数湍流直接数值模拟和高雷诺数工程湍流RANS模拟中,包括低速槽道湍流、超声速均匀各向同性衰减湍流、二维亚声速翼型湍流和三维跨声速翼身湍流等。数值结果表明该格式对湍流直接数值模拟和湍流RANS模拟具有高数值精度和高数值稳定性。下一步将利用高阶气体动理学格式研究更具有挑战性的可压缩湍流问题,例如超声速湍流边界层和激波边界相互作用等。     

超声速横喷干扰湍流场数值模拟

本文探讨了超声速平板横向喷流干扰湍流流场的数值模拟方法。分别采用Roe、Harten-Yee及AUSM 格式对三维雷诺平均Navier-Stokes方程进行离散,湍流模拟应用Baldwin-Lomax代数模型,讨论了B-L模型中Fy最大值的不同选取方式。分析了数值格式、网格划分对计算结果的影响。将来流马赫数为5、喷流马赫数为3的实验结果与计算结果进行比较,选择了与实验符合较好的数值模拟方法,为进一步利用数值计算方法研究喷流干扰流场特性奠定基础。     

三角翼大攻角绕流场数值模拟研究

本文基于N-S方程采用Jameson的显式中心有限体积格式和Roc的上风格式对三角翼大攻角绕流场进行初步数值模拟，空间网格为H-0型网格。计算结果表明两种计算格式均可较好的模拟大攻角绕流场。本文计算结果为今后开展全机大攻角绕流场数值模拟提供了计算经验。     

结构动力学方程的显式与隐式数值计算

从理论上详细推导了动力学方程的显式隐式时间积分数值计算方法,针对经典中心差分格式在求解阻尼结构计算效率降低的问题,研究提出了一种新的显式积分格式,理论证明该数值积分格式能够处理各种复杂阻尼状态而不降低计算效率。对于不同的时间步长通过设定调节加速度参数,能得到比经典显式计算格式和NEWMARK逐步时间积分更高的精度,最后通过数值算例对三种计算格式的理论分析进行了验证。     

基于动力学分析的发动机测振截面选取

研究了发动机测振截面选择的建模方法和优化算法。确定测振截面选择思路,建立转子-支承-机匣整机模型,计算正常和故障条件下转子动力学特性,并依据总体方案规定的传感器数目,确定测振物理量较大的转子截面或支承位置,完成振动传感器布置方案的初步确定。进一步提出测振方案优化方法,通过对实测数据进行分析,考虑实测信号的信噪比和信息冗余性,确定各测点振动数据有效信息量,判定测振方案是否合理,并选择机载测点。以某型发动机为例,设计监测点布置方案。结果表明:机载测点与该型发动机实际机载传感器测点位置相同,验证了测振方案的建立方法、优化方法和机载传感器选择方法正确可靠。单元敏感度越高,反映转子振动越明显。振动相似度越低,反映发动机整机振动特性越全面。研究结果为航空发动机测振截面布置提供重要依据。      

矩阵的奇异值分解 在红外光谱预处理中的应用

研究了矩阵的奇异值分解在红外光谱预处理中的应用。对奇异值序列,采用间距最大最高峰原则,选择有效秩的阶次;对原始光谱进行奇异值降噪处理、一阶微分处理,用相关系数法选择建模谱段为3 100~2 650 cm~(-1)和1 600~1 250 cm~(-1);采用PLS法建立运动粘度的校正模型,并与卷积平滑等预处理方法进行比较。模型校正集的相关系数(R_c)为0.977 9,标准偏差SEC为0.100 5,预测集的相关系数(R_P)为0.941 2,标准偏差SEP为0.154 7。研究结果表明:奇异值分解和一阶微分相结合可有效去除光谱噪声和基线漂移的干扰,提高PLS分析模型的准确度。     

三角形耦合本振阵列天线波束形成技术

以对称三角形耦合振荡器阵列为本振阵列,分析了耦合本振振荡器的相位控制原理,推导了耦合本振稳定同步形成均匀的平面相位分布,且反对称地调节边界单元的自由振荡频率,可改变本振信号的线性相位分布,确定了加权矢量,得出波束形成方向图,并通过仿真验证了三角形耦合本振阵列天线阵接收波束形成的可行性。     

麦克风阵列安装偏差引起的旋转声源定位误差

对旋转声源定位中3种常见阵列安装偏差即角度偏差、x轴向偏移偏差、z轴向偏移偏差导致的定位误差进行了详细研究。结果表明:当角度偏差在-10°～10°范围内,声源定位位置误差随角度偏差呈近似线性变化,最大定位位置误差为0.089 m,定位强度误差呈随机变化,其中最小定位强度误差为0.97 dB,最大为4.69 dB;当x轴向偏移偏差在-0.1～0.1 m范围内,声源定位位置误差随x轴向偏移偏差呈近似线性变化,最大定位位置误差为0.098 m,定位强度误差呈随机变化,其中最小定位强度误差为0.91 dB,最大为4.94 dB;当z轴向偏移偏差在-0.1～0.1 m范围内,最大定位位置误差为0.01 m,正偏移偏差引起的定位强度误差总体小于负偏移偏差,其中最小定位强度误差为0.81 dB,最大为4.51 dB。研究结果可为麦克风阵列在实际应用中控制测量误差提供指导。      

起落架缓冲支柱柱塞组件故障分析及改进研究

某起落架缓冲支柱修理时发现柱塞组件的螺母螺纹发生咬合变形,为了探究该故障发生的原因,根据柱塞组件的结构特点与工作原理,构建螺纹咬合变形故障的故障树,基于历史试验数据、螺纹剩余强度分析以及产品状态检查等方法对故障树底事件逐一分析并进行改进研究。结果表明：螺母设计缺陷和螺母未装配到位是导致螺纹发生咬合变形的原因;将螺距改为...     

等离子体对翼型流动分离控制历程的PIV试验研究

采用粒子图像测速（Particle Image Velocimetry,PIV）技术,研究了介质阻挡放电等离子体激励对NA—CA0015翼型表面流动分离的控制特性及控制效果随时间历程的变化规律。结果表明,激励电压存在一个阈值,当电压小于阈值时,控制无效或效果不明显;当电压接近阈值时,控制表现出不稳定性并最终趋于稳定;当电压大于阈值时,控制效果稳定且显著,气流能够很好地重附在翼型表面。     

基于确定接触迹线的点接触齿轮副构建及实验

提出一种点接触齿轮副,在已知齿面上根据设计需要确定接触迹线,推导其共轭曲线,将齿廓曲线沿共轭曲线扫掠后形成轮齿齿面,构建与已知齿轮正确啮合的配对齿轮,并推导该齿轮副的滑动系数计算方法;以渐开线内齿轮为已知齿轮,设计以抛物线齿面作为轮齿齿面的配对齿轮,进行齿轮样机制造和效率实验研究。结果表明,该齿轮副在啮合过程中时刻保持点接触状态,啮合点沿着理论接触迹线移动,与理论分析一致;滑动系数的大小取决于齿轮两端选取啮入点和啮出点所对应的渐开线参数值范围;实验后齿轮副实际接触迹线与理论接触迹线一致,且效率稳定在97.2%~98.5%。      

聚四氟乙烯接枝GMA的反应研究

使用叔丁基锂(t-BuLi)与乙二胺(EDA)的混合溶液活化聚四氟乙烯(PTFE),并在活化的PTFE上接枝甲基丙烯酸缩水甘油酯(GMA)。使用傅立叶红外吸收光谱(FTIR)对活化与接枝后的PTFE进行了分析,结果表明:本工艺成功将GMA接枝到PTFE上。通过扫描电子显微镜(SEM)观察了接枝表面,使用TGA法测定接枝聚合物的接枝率达到5.83%。最后研究了接枝反应时间和温度对PTFE-g-GMA接枝率的影响。结果表明,接枝率随着反应时间的延长逐渐增加,反应时间超过8h后基本恒定;在70℃反应时接枝率出现最大值。     

航空发动机涡轮机匣传热分析技术研究

通过数值仿真和基础试验系统研究了涡轮机匣的传热特点,建立了高压涡轮机匣传热分析方法。总结机匣传热分析的主要技术环节,开展了机匣关键部位的CFD流动换热数值分析,给出了机匣表面换热分布规律,通过换热试验修正了机匣表面换热分析方法;通过机匣的过渡态温度分析,建立了涡轮机匣传热分析工程方法。     

直升机机身气动外形的低阻优化设计

直升机机身阻力是飞行阻力的主要来源之一,通过对机身外形的优化设计,能够实现直升机的高效低阻飞行.首先,把机身划分为头部、中段和尾梁三段,对其外形轮廓线进行CST参数化表示;其次,采用拉丁超立方法选取试验设计样本点,计算各样本点的阻力系数,构造Kriging代理模型,估计模型预测的精度;最后,选用序列二次规划算法对其进行优化,并对优化后的机身模型进行了风洞试验.通过计算分析可知:所建立的Kriing代理模型能够精确预测阻力系数值,优化后得到了机身的设计参数;机身阻力系数减小了15.3％,理论值与试验值吻合良好.