战术导弹横向喷流数值模拟

为了研究喷流对导弹不同升力面的干扰影响,通过数值求解N-S方程来模拟超声速外流场中横向喷流的干扰流场,采用分块对接网格和"O"型网格技术,精确模拟喷口截面及弹翼形状,生成高质量的贴体计算网格.通过几种升力面和弹身组合的超声速横向喷流干扰流场的数值模拟,计算的马赫数为2和4.5,计算结果显示具有升力面的喷流干扰比光弹身的喷流干扰大.研究和分析了喷口附近流场的涡系结构和波系结构,并将导弹喷流干扰力放大因子及喷流引起的压心变化的计算结果与目前国外公布的试验值进行比较分析,二者显示了良好的一致性.     

有限体积法中面积分离散格式的精度分析

讨论了应用顶点中心格式对欧拉方程进行有限体积离散时的精度,详细推导了顶点中心格式在平滑网格和扭曲网格上的截断误差,同时与中心格式和顶点格式进行了比较.格式精度分析采用泰勒级数展开的方法.分析结果表明:顶点中心格式和顶点格式类似,总能达到至少是一阶的精度,并且可以通过加密网格的办法减少误差,而且顶点中心格式更适合于分区计算.通过圆弧亚音速绕流和NACA0012翼型绕流流场的数值计算,验证了分析的结论.     

基于低温辐射计的探测器绝对光谱响应率测量

探测器的绝对光谱响应率是衡量光电探测器性能优劣的重要参数,直接反映其光电转化的能力。通过对低温辐射计,功率稳定系统光路调整的介绍,基于低温辐射计,对探测器可见光范围的7个波长的绝对光谱响应率进行了高精度测量,分析了测量结果。     

镁粉的高温水反应特性研究

采用管式炉模拟发动机的高温环境(600～900℃),研究了镁粉与水的反应特性,得到了反应率与温度、镁粉粒度的关系和添加剂对反应的影响.结果表明,Mg/H2O反应的反应率随温度增加而增大,随粒径减小而增大,添加剂对反应有促进作用,Mg/H2O反应的表现活化能为24.67 kJ/mol.     

基于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,有效抑制高磁场散度区域的散度累积,同时也会将误差传播到低散度区域,并且引起边界处散度的累积,影响计算的稳定性。      

氧浓度对同轴乙烯层流扩散火焰中碳烟形貌及粒径演变的影响研究

为探究氧体积分数对乙烯同轴射流扩散火焰中碳烟颗粒生成及演变过程的影响,采用SiC纤维沉积法和热泳探针采样法,对层流扩散火焰中不同径向和轴向位置处的碳烟生成特性进行了研究。研究表明,氧体积分数增加,使火焰同一位置温度升高。在氧体积分数低于31%、火焰高度低于30mm时,SiC纤维上碳烟沉积物形态由火焰中心位置处的表面光滑的类液态演变为粗糙的凹凸块状,随温度升高及氧化作用加强,渐渐被氧化为疏松的团簇状、絮状,最终过渡为致密的团簇状和纤维网状。相同氧体积分数下,碳烟颗粒平均粒径随火焰高度增加呈现先增大后减小的趋势。在氧体积分数21%、火焰高度30mm处初生粒子直径达到最大值,为41.8nm。在火焰根部,碳烟颗粒平均粒径随氧浓度升高而增大,而火焰较高位置处则呈现了相反趋势。     

RTM成型聚酰亚胺复合材料的热氧老化特性

采用树脂传递模塑(RTM)工艺制备了U-3160碳纤维增强HT-350RTM聚酰亚胺树脂基复合材料(U-3160/HT-350RTM),研究了不同老化温度、老化时间下U-3160/HT-350RTM复合材料的失重率的变化规律,建立聚酰亚胺复合材料老化失效特征与老化时间/老化温度的关系,并通过微观形貌分析阐述了其在热氧老化过程中的失效机理。结果表明:在一定温度下复合材料的失重率变化符合三次多项式的变化规律,复合材料的老化在材料近表面尤为明显,由于氧分子作用,聚酰亚胺树脂发生降解导致孔隙率增加,因此温度越高、老化时间越长老化加速现象越明显。     

基于自适应混合网格的脱体涡模拟

基于混合网格和CGNS（CFD General Notation System）数据结构,建立了一种各向同性加密/稀疏的网格自适应方法。在悬空点的后处理中,让含有悬空点的单元转化为任意多面体,从而简化了自适应单元剖分模版,同时自适应网格单元之间可完全相容,自适应生成的网格能够直接用于可处理任意多面体的流场求解器。将该自适应方法与脱体涡模拟（DES）算法相结合,开展了65°后掠三角翼大迎角流动的数值模拟应用,并与初始网格的模拟结果进行了详细比较。对比表明：采用网格自适应方法适当增加局部网格量,能够以较小的成本迅速提高三角翼背风区的空间分辨率,增强数值模拟对小尺度涡系结构的解析能力,从而弥补了基于混合网格的脱体涡模拟中常用二阶格式计算的空间分辨率相对偏低、不利于湍流多尺度结构精细模拟的不足。     

后掠机翼欧拉方程数值模拟

本文对后掠机翼用保角转绘加剪切变换生成C-H型网格,并在这类网格拓扑上,用Jameson的三维欧拉方程有限体积法,四步Runge-Kutta时间推进格式,研制出可供分析三维后掠机翼亚、跨、超声速绕流的计算程序。本方法的特色是改进了机翼表面网格点的分布,使机翼后缘与网格线一致。在数值计算中,通过附面层位移厚度计算,引入了附面层粘性修正,改善了计算结果。     

基于体积分数形式的多介质流动数值模拟

把非守恒形式的体积分数方程与Euler方程组耦合，用交错中心型格式求解流体力学方程组。对一维情况下两种气体的相互作用和二维气泡的变形及塌陷问题进行了数值模拟，并把模拟结果与γ-model方法和NND（Non—oscillatory and non—free—parameter dissipation difference scheme，无振荡无自由参数耗散）格式的模拟结果进行比较，达到了抑制界面两侧非物理振荡和保持体积分数正性的目的。