火星大气非平衡条件下驻点热流计算公式的探讨

针对目前绝大多数火星再入飞行器驻点热流公式未考虑热力学非平衡效应故而适用性及准确性存疑的问题,分析了典型的零攻角驻点热流公式Fay-Riddell公式在火星再入飞行器热力学非平衡驻点热流计算中的局限性;然后采用数值计算方法,对涵盖火星再入飞行器飞行状态的大范围计算状况开展了热力学非平衡模拟,并对模拟得到的驻点热流公式进行拟合,得到了适用于火星再入飞行器热力学非平衡条件下的零攻角驻点热流计算公式;利用Mars Pathfinder飞行数据对公式进行了验证,计算表明,对于验证的计算状态,该公式的计算误差小于10%,符合工程预估的可接受误差范围。     

玻璃纤维增强泰氟隆烧蚀理论

本文给出玻璃纤维增强泰氟隆材料的理模型和计算方法。得到了有液层机械流失的传质参数的关联公式，利用此公式可对此材料的驻点与非驻点烧蚀性能进行预测。     

超音速边界层转捩的数值预计

本文对Wilcox和Traci提出的湍流模型进行了合理的改进,使改进后的湍流模型既适用于转捩的后阶段也适用于转捩的初始阶段。对于零攻角圆锥超音速绕流的边界层转捩,本文作了系统的数值实验。数值计算结果与风洞实验数据比较,在各主要方面二者符合的程度是令人满意的。     

航天器再入大气层热力分析

以OREX（orbital reentry experiment vehicle）的飞行试验数据和相关的CFD数值模拟结果为基础,采用传热理论及相关公式,分析计算了OREX再入大气层过程中的轨迹,驻点处热流密度非平衡假设和平衡假设下的换热问题.所计算的轨迹、热流密度非平衡假设下计算的驻点温度和热流密度值同试验数据及相关的CFD数值结果取得了很好的一致,相应的计算方法可作为航天器驻点热力分析的通式.然而在平衡假设条件下,尽管计算所得到的驻点热流密度与之前的CFD数值结果差别不大,但这种情况下计算得到的温度与试验数据不符,这应归结于计算的热流密度结果对驻点处温度变化的不敏感.比较非平衡假设和平衡假设下的换热计算结果表明,对于航天器再入过程中的热力探讨不能仅仅满足于热流密度分析,对温度的考察或许更重要.      

中心进气旋转盘流动与换热的数值研究

本文采用混合长度模型及共轭计算的方法对中心进气转静系旋转盘、腔内的流动与换热进行了数值计算, 并与实验结果进行了比较, 证明混合长度模型对转静系旋转盘、腔的流动与换热的计算是可行的。计算结果显示: 冲击区努赛尔特数最大值不在驻点, 而是偏离驻点一段距离; 随着转速和冷气流量的增加, 盘面换热总体加强。      

基于轴对称比拟的高超声速复杂外形气动热预测方法研究

轴对称比拟法是高超声速飞行器气动热计算的一种有效方法。针对轴对称比拟法应用时公式推导繁琐、变量迭代复杂、计算量大的缺陷,直接在笛卡尔坐标系下采用三维线性方程拟合物面方程及流场变量,推导出相应的流线尺度因子计算公式,扩展了驻点区热流密度的计算方法,提出了一种驻点区下游流线推进格式。以改进的轴对称比拟法为基础,将边界层外无粘流场数值方法与边界层内气动热工程算法进行耦合,发展了一套适用于三维复杂外形飞行器的气动热计算方法。通过对球头钝锥和双椭球算例进行验证,结果表明:方法计算效率较高,适用范围较广,热流计算结果和实验数据吻合良好。     

基于流线跟踪法的气动热工程计算研究

运用结构化网格求解三维Euler方程,计算得到边界层外缘无粘流场气流参数;利用无粘流场气流参数和表面流函数的方法计算了飞行器无粘表面流线分布;在理论和半经验公式的基础上,计算了定比热比和变比热比情况下驻点热流密度,非驻点区域采用参考焓、局部相似性等方法来确定飞行器表面的气动加热,实现了数值算法与工程算法的耦合.上述方法用于求解高超声速钝双锥的表面热流分布,计算结果与经典的热流公式和实验结果进行对比,平均精度为10%左右,满足高超声速飞行器概念研究和初步设计的需要.      

超声速小攻角下的圆锥转捩攻角效应

采用高阶差分格式求解三维可压缩滤波Navier-Stokes方程,对超声速零攻角和小攻角圆锥边界层转捩进行了空间大涡模拟研究.在三维Tollmien-Schlichting(T-S)波扰动作用下,计算得到了流场非定常转捩过程.计算结果清晰地展现了非定常大尺度结构演化特征和攻角效应引起的非对称转捩现象,流场发卡涡、横流涡等转捩主控大尺度结构与实验一致.     

反舰导弹三维流场数值模拟

三维流场的数值模拟是国外军事技术发达国家在飞行器设计过程中所采用的重要方法.本文建立了导弹的真实三维几何外形,然后结合实验数据,采用有限体积法,对国产某典型反舰导弹在高亚音速、不同攻角时的三维流场进行了数值模拟.计算出了导弹在不同工况下的阻力系数、升力系数和力矩系数.把计算所得结果与实验值进行了比较,用理论公式对计算结果进行了对比分析,并简单地分析了偏差的原因和解决的可行措施.文中所用的残差光顺技术和网格自适应技术能有效地提高计算的收敛速度和精度.     

热流CFD计算中格式和网格效应若干问题研究

通过数值实验、理论分析、试验数据和经验公式的综合研究,探讨了CFD热流计算的格式效应及网格效应:研究了Roe的FDS格式、Van Leer的FVS格式、AUSM 格式及中心差分格式的热流计算分辨率,探讨了不同格式对热流计算精度的影响,认为AUSM 格式在热流计算精确性方面具有优势;在热流绝对值的模拟上,网格的影响很大,能否计算出正确的驻点热流值将是获得精确物面热流解的关键;通过大量的数值计算、理论分析及其同试验数据的对比研究,对可能影响热流计算的网格因素进行了细致的分析,揭示了各因素对热流计算结果的影响规律。     

机翼大攻角气动力数值计算

本文描述了一个计算任意平面形状的机翼在无粘、亚音速、大攻角情况下的气动力数值方法,计算大攻角下气流分离时的定常非线性气动载荷。考虑翼面具有脱体涡分离,对旋涡强度和自由涡的位置进行松驰透代,取得收敛解。计算结果与实验数据很接近。     

超声速溢流条件下二元超声速进气道附加阻力计算

通过几何关系导出了单楔、双楔和三楔超声速进气道在零攻角及有攻角时的附加阻力系数计算公式,并进一步得出了N-1楔进气道的附加阻力系数计算通用公式.通过算例研究分析了飞行马赫数、飞行攻角及进气道总转角对附加阻力的影响.分析表明,进气道在超声速溢流条件下,附加阻力只同飞行马赫数、攻角及各楔面转角有关;附加阻力随飞行马赫数增加而减少,随攻角及进气道外压缩面总转角增大而增大;对于加速爬升用冲压发动机而言,设计时应注意减少附加阻力,并结合弹道、气动外形合理选择攻角.     

翼型和叶栅前后缘局部流场的新精化解法

本文提出一种利用前后驻点自动定位来求前后缘局部流场精化解的方法 ,即在全流场解的基础上 ,把局部流场解分为两部分的叠加 ,其中之一是驻点邻区的近似解析解 ,它吸收了原解中变化急剧的部分 ,使剩余部分成为变化比较平缓、具有更好的解析性态的分布 ,因而易于求准确。算例表明 ,本文方法为提高亚声绕统计算精度 ,特别是带攻角绕流工况 ,提供了有效的手段 ,本文思想也可同变分差分法结合应用     

低速实壁风洞测压实验洞壁干扰修正研究

 本文用壁压信息矩阵法对不同绕流特性的。测压数据进行修正,均获满意结果。 文中证明了洞壁对模型驻点的压力不产生干扰;分析了测压实验洞壁干扰特性与规律,对不同修正公式作了对比;进而说明不能用对速压修正的方法消除洞壁干扰     

FLUENT软件的二次开发及其在火箭气动计算中的应用

针对火箭设计过程中的气动特性计算,在计算流体力学软件FLUENT的基础上进行二次开发,从后台启动FLUENT进程,生成参数文件,采用批处理功能自动控制计算过程,减少用户工作量,极大地提高了设计效率。以两种不同的火箭构型为例,计算了在攻角α=4,°Ma=0.4,0.6,0.8,1.0,1.2,2.0,4.0,6.0状态下的气动特性。给出了火箭升力系数、压心系数和阻力系数随马赫数的变化规律。计算结果与实验数据吻合良好,精度满足设计要求,计算的气动数据可以为火箭的初步设计提供参考和依据。     

飞机亚音速非对称运动流场数值计算

本文以飞机亚音速对称运动（仅小攻角）流场数值计算为基础。采用涡格法理论开发飞机亚音速非对称运动（小攻角、小侧滑角、小角速度）流场特性数值计算。本文特点：从毕－萨公式出发推导出发推导的涡格法流场诱速公式比面涡法简便；采用多坐标系，详细推导了各坐标系转换公式；将复杂的干扰流场简化为一阶干扰和二阶干扰效应，并仔细推导了相应的流场诱速计算公式。已在ＩＢＭ－４３４１计算机上建立了流场计算实用性程序，既可单独     

翼型大攻角绕流的数值模拟

以求解二维N-S方程数值模拟NACA0012翼型大攻角状态的可压绕流特性;N-S方程是在贴体坐标系中给出的,以代数方法生成C型网格系统。采用LU-ADI格式计算,为提高格式的稳定性在隐式和显式部分分别添加了2阶和4阶人工粘性项。应用BaldwinLomax湍流二层代数模型模拟了大攻角时前缘分离涡的形成,旋涡对流及其非定常现象。在某些Mach数和攻角下NACA0012翼型的湍流解具有周期性。通过比较,本文数值计算结果同实验及国外相应的数值计算结果基本吻合。     

在正交流线坐标系中数值计算二维升力翼型跨音速绕流

本文将文献[1]提出的流线迭代法进一步应用于计算正交流线坐标系中有攻角时的二维跨音速绕流。用该算法计算了α=1°情况下的6％双圆弧翼型和NACA0012翼型,计算结果与实验数据符合较好。     

人椅系统绕流流场的数值模拟

以某种成熟型号的弹射座椅为研究对象,对其进行合理简化、改进,获得了接近实体的人椅系统几何模型.借助CFD软件FLUENT,采用控制容积法,对N-S方程、壁面函数、标准湍流方程进行了求解.在验证计算数据与实验数据相吻合的基础上,考察了马赫数为0.2,0.6,1.2,攻角为0°60°,偏转角为0°30°条件下可压缩人椅系统绕流流场情况,得到了人椅系统绕流流场的流速、压力、马赫数、表面压力及气动力参数分布情况和变化规律.结果表明:高马赫数、零攻角、大偏转角对人体造成的损害程度最大.      

涡轮叶片型面换热系数计算

本文用数值方法计算了无气膜冷却涡轮叶片上包括前驻点的整个型面上的换热系数。驻点区采用相似方程计算, 并用圆柱绕流近似表示驻点附近流场。较成功地模拟了边界层转捩和过渡区的换热情况。计算中考虑了来流湍流度的影响。