CFD数据交易标准及其软件实现

论述推行统一的CFD数据交换标准的必要性 ;报告在研制通用航空气动力数值模拟系统ANSS2 0 0 0过程中 ,建立内部数据接口文件规范的尝试 ;推荐国际上主流的空气动力学和CFD数据库标准———CGNS ;讨论CGNS在ANSS2 0 0 0中的实现 ;大力呼吁同步使用与国际接轨的CFD数据交换标准。     

叶轮机械三维不可压Euler方程的直接解法

本文用近似因式分解交替方向隐式格式直接求解叶轮机械中原参数三维不可压Euler方程组。在连续方程中引入“拟压缩性”之后,将原来类型不确定的控制方程变成双曲型方程,给定初边值之后用时间推进法求得稳定解,对一个单级压气机转子内部流场进行了计算,结果与实验值比较符合。     

求解Euler方程的无网格与有限体积混合算法研究

研究了无网格方法与有限体积法相结合的、求解Euler方程的混合算法.与单一的无网格方法相比较,由于在大部分计算区域采用了有限体积法,使得发展的算法在运算效率上能与有限体积法相当;同时在物体附近嵌入了无网格区,使得在处理几何外形上更具灵活性.有关算法涉及的无网格与网格区域搭接处理,通过在交界面处引入辅助卫星点和网格单元构成边界信息,实现了区域间的流动信息传递.Euler方程的空间导数分别在两区域用有限体积法和无网格法离散近似,时间方向都采用四步显式Runge-Kutta格式推进求解,数值模拟了喷管内流和绕翼型外流,并分别与整体有限体积法和整体无网格方法进行了比较.算例展示出,提出的混合算法能有效捕捉激波间断,且两区域等值线过渡光滑,算法效率如预期与有限体积法相当,表明该方法是可行的.     

含内热源球壳的复合传热特性及其表面温度均匀性研究

应用计算流体动力学方法,对带内热源的空腔球壳所涉及的含导热、自然对流和辐射换热的耦合传热问题进行了数值建模,得到了球壳空腔内、外侧流体的温度场、速度场。在数值模拟基础上,提出并试验验证了能有效提高球壳外表面温度均匀性的电加热绳方案。     

当量比对连续旋转爆轰燃烧室特性的影响

为了探究当量比对甲烷-空气连续旋转爆轰燃烧室（CRDC）特性的影响,利用二维可压缩欧拉方程对CRDC进行了数值研究,分析了爆轰波的发展过程和贫燃熄火过程,对比了不同工况下CRDC特性参数的变化情况。结果表明:CRDC起爆后燃烧场在由不稳定状态到相对稳定状态的过程中发生了2次碰撞,当进气当量比较低时,CRDC未能完全发生2次碰撞过程就已经熄爆。随着进气当量比的降低,爆轰波传播速度、轴向平均速度、出口平均温度、出口平均总压均呈下降趋势;增压比随当量比降低而减小的根本原因在于旋转爆轰燃烧过程和等压燃烧过程的熵增差减小,使吉布斯自由能增量差减小。CRDC的燃料驻留时间处于亚毫秒量级,燃烧热效率保持在99%以上。      

旋转整流罩积冰过程中水膜流动和脱离的数值研究

为了更准确地模拟旋转整流罩积冰过程,在现有三维积冰与冰层表面薄水膜流动耦合模型基础上,基于功平衡分析的方法引入了旋转部件表面水膜脱离模型,并发展了相应的计算方法,给出了水膜脱离的判定依据:当气流曳力做的功和由于离心力使水膜增加的潜能之和大于黏附功时整流罩表面的水膜会发生脱离。对旋转整流罩积冰进行数值模拟,计算结果与实验结果吻合得较好,验证了该模型的合理性和计算方法的可行性。之后分析了转速和来流速度对整流罩表面水膜脱离和积冰的影响,结果表明:转速和来流速度越大,水膜发生脱离的比例越大。在研究范围内,转速为3000r/min和6000r/min时,因水膜脱离导致积冰总量分别减少13.4%和15.8%;来流速度为40、50m/s和60m/s时,因水膜脱离导致积冰总量分别减少为12.2%、13.4%、14.2%。      

变截面旋转裂纹梁横向振动特性的研究

对损伤结构进行动力特性分析是进行无损检测的重要基础,而对于旋转梁结构的裂纹损伤动力特性研究,却鲜有文献涉及。以变截面旋转裂纹梁为研究对象,对其横向振动特性进行研究,提出一种求解变截面旋转裂纹梁横向振动特性的新方法。首先利用扭转弹簧模拟裂纹效应,建立含裂纹梁局部柔度模型,然后采用Frobenius方法求解振动方程,得到方程的级数解析解,并研究裂纹位置和深度对振动频率的影响,分析不同损伤程度、不同转速工况下梁的前两阶固有频率变化情况。结果表明:本文方法是有效的,转速和损伤程度的变化并非独立影响梁的固有频率,两者间具有耦合作用机理,对于变截面梁同样成立。     

感性耦合夹层等离子体隐身天线罩电磁散射分析

设计了一种石英夹层感性耦合等离子体（ICP）隐身天线罩模型,采用有限元与Z变换时域有限差分（ZT-FDTD）联合仿真的方法,建立了ICP放电的流体模型,得到不同气压及功率下与电磁散射相关的电子密度空间分布,在此基础上建立Z变换时域有限差分法模型,对石英夹层等离子体隐身天线罩的宽频段后向散射进行计算,同时利用微波干涉法及XFDTD软件对算法及程序进行验证。结果表明：感性耦合等离子体可产生较高密度等离子体,能有效实现雷达散射截面（RCS）的减缩,在气压为2 Pa时,碰撞衰减较弱,等离子体密度分布较均匀,衰减带宽集中在等离子体振荡频率附近,功率增加会使衰减峰值向高频方向移动,气压为20 Pa时,碰撞衰减增强,且等离子体密度分布有较大梯度,衰减带宽有效增加,同时RCS曲线的波动特性加强。     

超临界压力下碳氢燃料裂解与流动传热模拟的快速算法

为了快速预测发动机冷却通道内碳氢燃料在考虑热裂解时的流动传热特性,基于正癸烷热裂解反应机理,建立了一套模拟正癸烷裂解吸热和超临界压力传热现象的快速算法。采用三维物性库查表算法计算裂解反应混合物的热物性,同时对组分输运方程进行简化,简化后只需求解1个组分输运方程。通过与现有的实验和数值结果进行比较,检验了快速算法的计算效率和可靠性。结果表明,快速算法与求解全组分输运方程的算法精度相当,但计算效率提升了约20倍。最后采用该算法对三维矩形截面冷却通道内的超临界压力正癸烷裂解与传热过程进行数值模拟,进一步考察了本文快速算法的计算精度及其工程应用价值。     

剪切稀化非牛顿射流撞击液膜破碎直接数值模拟

液态射流撞击是液体火箭推进系统中广泛采用的一种燃料雾化方法,其破碎特征直接影响燃料最终的掺混及燃烧效率。采用直接数值模拟（DNS）工具,研究了低雷诺数（Re_l=41）和中等韦伯数（We_l=163）条件下剪切稀化非牛顿射流撞击液膜破碎的问题,着重分析了对角液膜的三维结构、破碎特征和非牛顿特性等。研究结果表明：在所研究的射流参数下,该非牛顿撞击液膜破碎属于Open Rim类别,破碎过程具有三维特性并伴随液丝与边缘的融合、液丝向液滴的转变等时域流动特征。液体的总表面积随时间不断增长,但单位表面积随液膜破碎的发生而下降,液膜扩张半角随时间逐渐增加并趋于恒定值43°,而后部液膜的长度几乎不随时间发生变化。此外,撞击液膜表现出明显的剪切稀化特性,液体内部最低黏性系数仅为零剪切黏性系数的1/5。