钝体燃烧器后能量逆转和紊流参数的数学模拟

本文用非定常数值法研究了钝体燃烧器后的紊流参数分布和能量逆转,并对Boussinesq的雷诺应力模式和双方程k-ε紊流模型提出了改进方案。这样不仅计算出了能量逆转区,而且紊流参数分布也与实验相符。     

喷打用红色陶瓷表面装饰墨水的制备与性能

制备了一种红色陶瓷表面装饰墨水——与传统陶瓷表面装饰的新技术(连续式喷墨打印法)配套的新型墨水。它以陶瓷颜料为基本着色料，用溶胶一凝胶法制备，以连续式喷墨打印墨水的质量标准为参照标准。经测试，所制墨水的电导率、表面张力、粘度等理化性能完全符合喷墨打印墨水的要求，适合用喷墨打印法陶瓷装饰工艺的需要。研究了分散剂、粘接剂、pH值、固含量等工艺因素对该墨水性能的影响。该墨水制成时为墨绿色，经煅烧，各组分反应生成颜料，显现粉红色。     

超声速流中激波/湍流附面层干扰数值模拟

采用修正的B／L湍流模型以及多块结构化网格求解了二维N-S方程。分别对超声速流和高超声速流中的激波／湍流附面层干扰进行了数值研究。本文首先研究了进口马赫数为2．96的超声速流。计算结果准确预测了入射斜激波在平直壁面引起湍流附面层分离的流动特征：分离点的反射激波、分离包引起的膨胀扇以及再附点的反射激波。计算的壁面压力分布与实验值吻合较好，计算的分离区长度与实验值比较有一定误差。本文还对进口马赫数为9．22的高超声速流中压缩角引起的激波／湍流附面层干扰进行了数值研究。计算结果与实验结果吻合较好。     

灰色系统理论中二维数据序列均值生成的讨论

灰色系统理论中非等时序序列或含有空穴的序列,通常可以用均值生成方法来填补,生成等时序序列。本文研究发现,用非紧邻均值生成方法生成的内点与实际存在有明显的误差,为此讨论了均值生成和有偏生成算法。利用数据序列点的局部凹凸特性来表示数据序列点的特征,提出了采用数据序列点连线的状态选取生成系数的方法。通过对一维数据序列传统均值生成方法的比较,给出了二维数据序列的定义,讨论了在二维数据序列下的均值生成方法,提出了有偏生成算法中生成系数分解的选取方法。     

高强驻波声场中的湍流

反射端声压级高于140dB以后驻波管中的声流结构会发生明显变化。利用流动显示和LDV（Laser Doppler Velocimetry）技术对该问题进行了实验研究。流动显示结果表明，在160dB左右的高强驻波声场中随着旋涡结构的合并、破裂等现象的出现，声致流动有明显的湍流特性。LDV测量结果显示，随着驻波管中声压级的提高，一些测点的速度频谱明显加宽，出现高频速度脉动，而另一些点的速度频谱则没有显著变化。径向湍流强度测量结果显示，此时管中流动已是复杂的三维结构,而非轴对称的二维结构。     

倒置喷管湍流流场计算

本文对倒置喷管的湍流跨音速流场进行了数字模拟，其算法基于ＹａｎＬｅｅｒ的矢通量分裂格式，求解时采用了线性Ｇａｕｓｓ—Ｓｅｉｄｅｌ迭代法，并结合局部时间步长技术对收敛过程进行加速，对雷诺数为１．６６×１０６的倒置喷管的流场进行了数值计算，结果与实验相符合。     

湍流燃烧模拟中化学反应的加速算法研究进展

为研究湍流燃烧数值模拟中化学反应机理计算的加速方法，讨论了动态自适应化学（Dynamic Adaptive Chemistry，DAC）方法和Krylov子空间近似的指数格式的应用情况。在湍流火焰大涡模拟中，使用DAC简化可以加速化学反应计算。然而，在并行燃烧数值模拟中，处理器核心的负载极度不平衡，加速效果有限。而Krylov子空间近似的指数格式的加速效果可以作用于每个处理器核心，更有利于整体计算效率的提高。在同等精度下，相比于隐式格式耦合DAC和MTS加速方法，Krylov子空间近似的指数积分格式对化学反应计算的加速效果更为显著。     

数值模拟方法对NASA CRM模型阻力预测的影响

在复杂工程外形的数值模拟中，网格类型、规模和分布、湍流模型、数值格式等都会不同程度影响模拟结果。如何评估这些模拟方法对模拟结果的影响，并识别对模拟结果影响较显著的因素，对于CFD的发展方向有积极的借鉴指导意义。为了综合研究不同因素对商用运输机外形阻力预测的影响，以AIAA第六届阻力预测会议外形NASA CRM为研究对象，同时考虑了网格、湍流模型、无黏通量格式和体心梯度求解方法等因素对阻力预测的影响。分析中采用枚举法和正交试验设计两种方法，并使用了基于聚类分析的定性敏感性分析方法和基于Mckay主影响分析的定量方法，识别出对阻力预测影响较大的因素，这为数值模拟方法的发展指明了方向。     

三类气动导纳数值识别方法的适应性研究

针对气动导纳函数的数值识别方法,借助于CFD,在简谐脉动来流、湍流和竖向阶跃来流下对平板断面和箱梁断面的导纳函数函数进行研究。首先,在无断面存在的空流域内详细研究了简谐脉动来流、湍流和竖向阶跃来流的传播特性及其数值计算方法。其次,对有断面存在的情况进行了数值计算。最后,识别得到了平板断面和箱梁断面在三种不同来流下的气动导纳函数。结果表明:对平板断面,三种方法识别得到的气动导纳函数与Sears函数吻合良好,验证了三种数值计算方法的可行性;对箱梁断面,简谐来流和湍流下识别的气动导纳差别不大。相比之下,完全基于线性叠加原理的阶跃来流方法产生了实质性的偏差,表明该法不宜用于钝体断面。计算效率方面,湍流的计算效率适中且对任意断面适用;简谐脉动来流的计算效率最低,适用于气动导纳与风场无关和弱相关的断面;竖向阶跃流方法具有计算时间短的优势,但它仅能用于气动导纳与风场完全无关的断面。     

高雷诺数湍流非预混火焰及NO生成的大涡模拟

采用两种亚网格湍流燃烧模型,即化学建表方法结合假定概率密度模型和稀疏拉格朗日过滤密度函数方法,对高雷诺数湍流非预混火焰Flame D进行数值研究,定量比较不同亚网格模型的差异,并对火焰特征、污染生物生成特性进行分析。结果表明,两类亚网格燃烧模型预测的温度及大组分分布相近,稀疏拉格朗日过滤密度函数方法可以更好地模拟CO质量分数分布。不同的假定概率分布均可合理描述湍流与火焰的相互作用,之间的差别主要体现在NO分布,Dirac函数远高估了NO生成,而Top-hat函数则略低估了NO生成,Beta函数表现最优。Flame D的高温区及NO质量分数均主要分布在当量混合线及富燃侧附近。受高温伴流的影响,NO质量分数与温度一直保持高度正相关,峰值主要集中在标量耗散率很小的区域。不同截面上,反应物中的O_2和生成物中的H_2O均与NO高度相关。