三维机翼结冰模拟

针对机翼结冰问题,利用欧拉法求解空气水滴两相流流场,通过Fluent软件用户自定义函数功能编程求解三维机翼表面水滴局部收集系数以及结冰情况.在计算中.将水溢流方向沿机翼弦向和展向分解,以此对Messinger结冰热力学模型进行改进,使之更符合三维情况,并在此基础上进行三维传质传热分析;修正传统的二维表面对流换热系数求解...     

小型平板回路热管蒸发器内传热与流动的三维数值模拟

针对小型平板回路热管蒸发器内的流动与传热,建立了多区域耦合的数学物理模型,并应用FLUENT软件进行了三维数值模拟。结果表明:蒸发器传热特性在不同热负荷下呈现出较大的差异,其温度分布不仅取决于热负荷,更依赖于毛细芯表面发生的两种传热机制,即毛细蒸发和热传导。相比高热负荷(Q=120W)和低热负荷(Q=40W),中等热负荷(Q=80W)下蒸发器各个部位的温度均较低。三种不同热负荷下,毛细芯反向导热均大于侧壁漏热,补偿腔内与毛细芯相邻处易出现高温区。冷凝回流液在补偿腔内的流动形成两个涡,这种流动特点有利于降低毛细芯的温度。当热负荷与系统冷凝能力匹配时,整个系统流动与传热特性最优。      

利用OpenGL建立运动目标体素模型的方法

三维体素模型的建立是运动目标序列图像产生的基础和关键步骤，三维体素模型就是利用计算机系统将运动目标的几何特性在计算机显示器上直观地表示出来。介绍了应用SGI公司的OpenGL三维图形库和微软的MFC基类库在VC 5.0的编程环境中建立运动目标三维体素模型的方法和过程，并以较复杂的静态飞机模型为例进行了模拟实现，得到较好的效果。     

燃烧过程辐射换热离散坐标模型的应用

对离散坐标(DO)辐射传热模型分别在长方体空间和旋流燃烧室内进行了验证与应用。长方体空间的计算结果表明:DO模型的纯辐射传热预报结果比六热流(SHF)模型准确,与计算精确的区域法的计算结果吻合较好,而SHF模型的预报结果却有较大的误差。对于旋流燃烧室内温度场的模拟预报也表明,DO模型的计算结果和实验值相吻合,明显高于用SHF模型计算所得的温度值。另外,DO模型在计算时间和边界处理等方面比区域法优越;在收敛速度和计算准确度等方面也比工程常用的六通量模型优越。      

三维坑湍流分离流动的粘流／无粘流干扰计算

本文对坑的三维分离流动做了低速粘流与无粘流的相互作用计算。对三维边界层反方法进行了分析和讨论。用数值试验的方法验证了在H和α作为已知量的情况下,三维边界层反方法的积分方程是双曲型的,并提出了一种近似数值特征线法进行求解。无粘流采用低速位流面元法。计算表明所用方法可计算出三维效应很强(即横向变化很大)的三维边界层分离流动。     

渤海冬夏季环流的数值研究

用三维陆架海模式对渤海的冬夏季环流进行了研究。指出：就垂向平均而言、渤海的冬夏季环流结构基本一致,且与垂向平均的潮余流相似。渤海环流有显著的三维结构,冬季,上层主要为风漂流,下层主要为风漂流,下层为补偿流;夏季;潮,风和密度对环流的贡献较大,表现为复杂的三维结构。     

可压混合层涡并放热效应的涡动力分析

 基于双元单步不可逆燃烧动力学简化模型 ,采用紧致差分算法直接求解二维可压 Navier-Stokes方程 ,研究了时间发展非预混扩散火焰平面反应混合层流动。拟序结构由数值线性稳定性分析方法设定的粘性可压最不稳定小扰动初场促生。随放热强度增大 ,得到了基频涡卷撕裂现象。涡动力分析表明 ,放热主要通过膨胀和斜压这两种强非线性过程抑制大尺度剪切涡的演化     

边界层分离流边值问题的算法研究

本文给出了边界层方法边值问题的一种快速求解方法,该方法可用于二维与三维问题,它是在正方法边界层算法的基础之上增加了一个矩阵算法,该矩阵很容易逆,总的ＣＰＵ时间仅增加约２０％,使用该算法对反方法问题求解和进行数值稳定性分析都将十分方便。     

三维跨声速流粘性—无粘相互作用问题的分析和计算

本文对跨声速粘性-无粘相互作用问题进行了分析,论证了如何选择湍流边界层反模式的输入量。在FL027多重网格基础上作了局部加密。计算例子表明,本文建议的半逆方法适用于弱激波情形,可用于实际设计之中。     

液体火箭发动机推力室冷却通道温度分层数值研究

为了研究冷却剂温度分层的形成机理及其对流动和换热的影响, 应用雷诺应力模型(RSM)对液体火箭发动机推力室再生冷却通道的流动与传热进行了三维数值模拟, 冷却剂为气氢, 考虑其物性随温度和压力的变化.所得结果表明:冷却剂在非流动方向会出现温度分层现象, 随着冷却剂的不断受热, 温度分层现象越明显, 由于喉部二次流加强了冷却剂间的混合, 在喉部区域温度分层被减弱, 温度分层对冷却剂温升及压降影响较小, 严重影响气壁温度的估算.      

基于传热反问题方法的N_2O/C_2H_4预混推进剂燃烧室热流测量研究

测量液体火箭发动机的热载荷是获取燃烧室内部信息的重要方法。为了获取N_2O/C_2H_4预混推进剂燃烧室内壁的热载荷,建立了液体火箭发动机的热流计算的反问题方法,该方法基于对燃烧室壁面温度场的直接求解,通过对轴向多个位置测量温度的反演计算得到燃烧室内壁热流和温度。研究表明:应用文中建立的传热反问题方法能够较为准确地获得热流随时间及空间的分布;热电偶的位置对计算准确性有明显的影响,与理论深度偏差在0.2mm以内的随机深度偏差可导致超过4%热流反演误差;N_2O/C_2H_4预混推进剂燃烧室热流及温度沿轴向逐渐降低,表明燃烧室内的反应释热过程主要在燃烧室头部附近发生。     

基于再入轨迹和气动热环境的返回舱烧蚀研究

针对再入全过程合理预测热防护罩表面材料烧蚀深度和温度的动态变化问题,提出融合再入轨迹、气动热以及Newton-Raphson和三对角矩阵算法（TDMA）构建动态烧蚀的方法。该方法建立直入式和跳跃式三自由度再入轨迹,应用修正的牛顿流体理论估算气动参数,以及修正的Fay-Riddell和Sutton-Grave理论计算驻点区域的热流密度,利用一维非线性热传导方程模拟了热防护材料的烧蚀过程。仿真结果表明：此方法实现了再入全过程热防护材料烧蚀深度和温度连续动态变化的预测,同样适用于更为复杂结构飞行器的动态烧蚀预测,与热平衡积分法（HBI）相比其结果可靠合理,为进一步优化热防护系统（TPS）提供了一定的参考依据。     

基于再入轨迹和气动热环境的返回舱烧蚀研究简

 针对再入全过程合理预测热防护罩表面材料烧蚀深度和温度的动态变化问题,提出融合再入轨迹、气动热以及Newton-Raphson和三对角矩阵算法（TDMA）构建动态烧蚀的方法。该方法建立直入式和跳跃式三自由度再入轨迹,应用修正的牛顿流体理论估算气动参数,以及修正的Fay-Riddell和Sutton-Grave理论计算驻点区域的热流密度,利用一维非线性热传导方程模拟了热防护材料的烧蚀过程。仿真结果表明：此方法实现了再入全过程热防护材料烧蚀深度和温度连续动态变化的预测,同样适用于更为复杂结构飞行器的动态烧蚀预测,与热平衡积分法（HBI）相比其结果可靠合理,为进一步优化热防护系统（TPS）提供了一定的参考依据。     

求半透明材料导热系数的复合换热反问题

在激光脉冲法测不透明材料导热系数的基础上,本文提出了一种已知半透明材料辐射特性,求其导热系数的复合换热反问题方法。      

推进动力系统燃烧不稳定性产生的机理、预测及控制方法

在不稳定热释放和声波相互耦合作用下而产生的燃烧不稳定现象,表现为振荡燃烧的形式,常发生于火箭发动机、航空发动机加力燃烧室以及地面燃气轮机等装置当中,会对结构引起很严重的破坏。本文总结了当前关于燃烧不稳定问题的机理以及控制方法的研究进展。目前来看,关于燃烧不稳定的机理研究方面,包含实验研究以及数值方法研究。实验方面,有rijketube热声不稳定实验,还有模拟真实燃烧室环境的燃烧实验。数值方面的工作,包括线化的热声不稳定模型,以及对火焰进行描述的解析模型和大涡模拟等方面的工作。同时,人们尝试了各种抑制不稳定的方法。控制方法包含两大类,即:主动控制方法以及被动控制方法。主动控制方法在理论研究方面取得了重大成果,然而,由于其实现上需要复杂响应系统、执行机构,因此,距离实际工程应用还很遥远。而被动控制方法,例如,亥姆赫兹共振腔,以及穿孔板等装置,在工程上得到了很好的应用。     

膨胀循环推力室再生冷却换热的数值模拟

为了解液体火箭发动机膨胀循环推力室再生冷却换热特性,采用数值模拟方法,研究了冷却剂流动方式、推力室圆柱段长度、圆柱段室壁加肋和气壁面粗糙度等因素对冷却通道压降、冷却剂温升、壁面热流密度和温度分布等换热特性的影响.计算过程中采用k-ε双方程湍流模型.计算结果表明:采取顺流冷却要比逆流冷却的冷却通道压降低,但同时冷却剂温升也低;对于室壁加肋结构,在肋个数相同而只改变肋高度的情况下,总换热量正比于总换热面积.      

一种高效的多重网格三维N-S方程计算方法及其应用

介绍一种基于 Jameson中心差分和 Runge-Kutta时间推进的三维 N-S方程计算的有限体积方法 ,计算点选取体积单元中心 ,为了提高收敛效率和程序稳定性 ,采用了当地时间步长和多重网格方法 ,特别是针对粘性计算中计算网格形状比高的特点 ,隐式残差光顺采用当地变系数的方法 ,并专门设计了考虑计算网格形状比的人工粘性模型。应用本文方法对收扩喷管的三维流场和NACA0 0 1 2翼型二维流场进行了计算 ,结果表明程序具有良好的收敛能力和较高的稳定性。     

燃气非平衡流再生冷却流动传热数值模拟

为准确预测液体火箭发动机推力室身部再生冷却换热状况,采用数值模拟方法,对燃气、推力室壁和超临界气氢进行三维流动和换热耦合计算。采用6组分9步反应动力学模型计算燃气的非平衡化学反应,采用DO模型计算燃气辐射换热,考虑超临界气氢物性随温度和压力的变化。获得了室壁温度场、燃气及冷却剂流场。结果表明,Redlich-Kwong方程、Peng-Robinson方程、Lucas法、TRAPP法能分别准确计算超临界氢的密度、定压比热容、粘度、导热系数,采用燃气非平衡流计算所得流场值更符合实际情况。     

基于特征正交分解的跨声速流场重构和翼型反设计方法研究

在二维流场的重构问题中应用特征正交分解(Proper Orthogonal Decomposition,POD)数据处理方法。利用CFD技术计算得到的流场快照对气动力模型进行降阶,然后利用基于POD的降阶模型(Reduced Order Model,ROM)对所需的流场参数进行重构,在快照范围内可以得到高精度的结果,且具有一定的外插能力。在翼型反设计问题中该方法仍然是成功的,通过修正快照向量,利用基于POD降阶模型的数据重构方法,由已知的翼型表面压力分布通过反设计就能够高效精确地得到对应的最优翼型形状,这极大地简化了翼型反设计问题。本文分别在跨声速范围对RAE 2822翼型的流场重构和Korn翼型及NACA 63212翼型的反设计进行了验证,证明了基于POD的流场重构和翼型反设计方法的有效性和高效性。     

基于欧拉方程的机翼设计方法及应用

以Takanashi提出的三维机翼设计理论为基础，研究与发展了一个基于欧方程和“正－反迭代、余量修正原理”的机翼设计方法。用改进的无限插值方法生成绕机翼的O－O贴体网格，采用三维欧拉方程作为充动分析计算的基本方程，该设计方法已用于某无人机机翼和一个超临界机翼设计，设计结果达到了预期的目标。