基于热-流-固耦合方法火焰筒壁温三维数值模拟

准确地确定火焰筒壁温分布是对其进行热结构分析的前提和基础.基于ANSYS/CFX程序平台,采用基于燃烧模型的热-流-固耦合有限元方法,对某环形回流燃烧室的整个流场、温度场进行数值模拟,得出火焰筒壁温三维温度场分布.计算中考虑了燃烧热源和热辐射,燃烧和流场变化对换热的影响,以及燃气、煤油和固体物性参数随温度的变化.数值模拟结果表明,将流场和固壁进行耦合求解,能更精确反映流场和温度场的整个形态,可以模拟出较为合理的流场和温度场分布以及相应的流动换热特性.所做工作对火焰筒结构设计具有重要的工程应用价值.      

3点逐步r型自适应网格算法

为了根据某个流场参数梯度较合理地分布网格同时又不增加网格点的数量,提出一种基以总变差为基础的r型3点逐步自适应网格算法,通过3点逐步调整网格点的位置以减小网格线的扭曲,生成质量较好的计算网格.以绕二维圆柱和三维双椭球的超声速流动为算例,用Euler方程模拟,以流场的压力总变差来进行自适应,通过分别用所提出的自适应方法生成的自适应网格与用原始网格计算所得的流场结果对比,计算结果表明,用自适应网格计算所得的激波比用初始网格计算所得的激波薄,采用自适应方法所计算出的流场数值解具有更高的分辨率.      

面齿轮齿面的自适应采样方法

面齿轮齿面的数字化是齿面检测的关键技术之一.针对在三坐标测量机上进行的齿面采样,提出一种根据给定精度确定采样网格点数量的方法:对采样网格边界线进行初步自适应,利用截平面法得到初始采样网格;同时根据给定的采样网格点数量,借助基于形状的采样算法,对初始采样网格进行迭代,生成自适应采样网格,最终实现面齿轮齿面数字化检测采样的自适应规划.     

关于非结构网格生成中的背景函数

背景信息对于用阵面推进法生成非结构网格是一个关键的环节.文章基于众多文献中提到的背景信息提供方法,提出一种更为实用的背景信息提供技术.此方法不仅避免了不必要的机时,内存的消耗,又使网格具有光滑性,网格的生成具有灵活性,方便性的特点.基于生成的网格采用有限体积法研究无粘高超音速化学非平衡流的流场分布,生成的网格能够成功的运用于高超音速化学非平衡流的数值计算。     

多翼导弹复杂流场及气动力的数值计算

将改进的Jameson格式应用到多翼导弹的三维流场计算中,并使用分区技术和三维椭圆微分方程数值网格生成技术以精确模拟此复杂构形的复杂流态,计算结果与美国导弹标模的实验吻合很好.一系列攻角下导弹三维流场的数值模拟表明,修正Jameson中心差分格式能够适用于复杂外形、复杂流态、亚/跨/超/高超音速的流场计算,而且精度和分辨率都比较高,结果可信.     

直升机发动机入口绕流的数值模拟

对直升机发动机入口处流场采用有限体积法,通过生成贴体结构化网格求解Navier-Stokes(N-S)方程, 进而对MI-171直升机低空悬停和前飞状态进行数值模拟,计算了该型直升机发动机进气道处的流场特性.然后以广义尾流入手,并引入尖涡涡核作用的半经验修正的旋翼自由尾流模型对该机发动机进气道处流场进行计算分析.对沿发动机轴线剖面以及垂至于该轴线剖面的分析结果表明,CFD(Computational Fluid Dynamics)数值模拟方法与旋翼自由尾流分析方法均可以对涡诱导流场进行预测.     

基于CE/SE法的翼型流场激波捕捉计算

在二维四边形结构网格下,针对时空守恒元和解元(CE/SE, space-time Conservation Element and Solution Element)方法捕捉激波需要双倍空间网格点的问题,重新对守恒元和解元进行了定义,计算点同时包含单元中心和网格节点,并推导得到了新的CE/SE方法的计算公式.以此为基础,结合非敏感克朗数(CNIS, Courant Number Insensitive Scheme)计算格式和当地时间步长方法,对激波翼型流场进行数值模拟,并与原方法及AGARD报告进行对比.结果表明,对原CE/SE方法的改进是有效的,可明显提高激波分辨能力,并且激波前后无明显的数值震荡发生,新方法适合应用于翼型流场中捕捉激波.     

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

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

混合重叠网格插值方法的改进及应用

混合重叠网格间合理插值是保证流场计算正确的基础之一。本文针对重叠区网格尺寸匹配度较差时插值误差较大问题,发展了一种新型混合重叠网格插值方法,通过使用二阶精度插值、按单元尺度区分和扩充模板,改善了插值区网格尺寸匹配度较差时的插值精度。本文方法适用于任意单元类型的混合网格重叠,对各类单元处理透明,实现简单。计算结果表明,采用本文方法在网格重叠区流场变量传递正确有效,插值区网格尺寸匹配度较差时,相比原始方法,等值线过渡更为光滑,变量经过插值区耗散更小,计算与试验值符合更好。     

燃烧室头部进气规律对煤粉燃烧效率的影响

用实验方法研究了火焰筒式煤粉燃烧室头部进气规律对燃烧效率的影响。通过对6种结构的头部配气方案的分析比较,得出了煤粉燃烧室合理的头部配气方案,为今后煤粉燃烧室的设计和工作状态的选择提供了重要的实验依据。     

基于虚网格的格心ADT搜索法

针对重叠网格寻点问题,提出了一种直接在格心网格下操作,基于虚网格的ADT(Alternating Digital Tree)搜索方法.通过对计算网格进行扩展,建立虚边界网格,解决了格心网格因覆盖区域不完整而无法直接建立ADT数据结构的问题.搜索结果即为包含解的格心单元集合,可直接对结果列表遍历以得到合理贡献单元,故完全摒弃了可靠性差的Stencil Walk方法.由虚网格的定义,使寻点在边界附近的处理更为灵活,可以准确给出边界附近贡献单元的有效信息,同时简化了虚网格体系的构建.扩展了重叠边界类型,构建的搜索空间能完整覆盖网格范围,解决了对称面重叠问题.算例研究表明:该方法可靠性好,边界处理能力强,有效提高了重叠网格方法对大型复杂网格的解算能力.     

多段翼型贴体网格的生成方法

本文给出了一种多段翼型绕流的分区贴体网格生成方法。采用泊松方程的解来给出贴体网格的坐标,泊松方程中的源项通过靠近边界的网格间距和正交性来决定,本文对分区网格的区域边界的选择作了讨论,采用生成的分区贴体网格计算了绕多段翼型位流的解,其结果与实验所得是很一致的。     

B样条应用于超限插值的网格生成方法

采用B样条曲线曲面生成技术与传统网格生成方法中超限插值技术,并将两者结合起,形成一种计算网格的生成方法.它发挥了两者的优点和长处:能充分反映复杂外形的完全信息,提高对外形的适应程度;可控制切矢的边界条件,提高生成网格的质量.算例表明了上述特点.      

双支板超燃冲压发动机燃烧特性研究

为研究分级喷注超燃冲压发动机火焰稳定、燃烧状态及火焰传播特性,以双支板超燃燃烧室为基本构型,开展了当量比连续调节试验研究。模拟低飞行马赫数5.5工况,燃烧室入口马赫数为2,总温1436 K,试验表明:燃烧室单独上游喷注熄火当量比为0.19,该值不受下游燃烧的影响;单独下游喷注熄火当量比为0.46,上游火焰会削弱下游当量比变化对壁面压力的影响,并且会使下游熄火当量比值降低。通过调节上游当量比可实现燃烧状态的转换,转换过程存在迟滞。模拟高飞行马赫数6.5工况,燃烧室入口马赫数为3,总温1 899 K,试验表明:随着总温的增加,单独上游喷注可实现点火和稳焰,上游火焰发生抬举,燃烧室抗反压能力增强,可喷注更多燃料。     

层流平板摩擦阻力的数值计算

为提高层流摩擦阻力数值计算的精度,采用数值求解N-S方程的方法,对层流平板摩擦阻力进行了计算.计算中采用4种先进上风格式,同时配合以14种不同的壁面法向网格分布和5种常用限制器,并将计算结果与布拉修斯解进行对比.结果表明,van Leer的FVS(Flux Vector Splitting)格式数值耗散大,对网格分布和限制器敏感,不适合用于摩擦阻力计算.限制器和计算格式的粘性分辨率越高,摩擦阻力的计算精度越高.在计算量允许的情况下,可以减小壁面法向第1层网格高度,以提高计算精度,但应同时兼顾网格长宽比.对于摩擦阻力计算,边界层内应至少布置4个网格点.      

三维椭圆型方程网格生成中的新源项修正法

构造一种新的源项修正方法用于三维椭圆型方程网格生成，可以实现对网格与边界正交及网格与边界间距的双重直接控制．     

复杂外形网格生成技术及应用

以某型双立尾战斗机为例,给出了实际复杂型号飞行器结构网格生成的一种分块方案及实现过程,用一种矢性三次多项式插值的方法生成相邻块公共交界面的网格,用Hilgenstock方法修正源项,求解椭圆型方程生成块内空间网格,使网格线与边界面夹角、与边界面间距控制为期望的值.在所生成的双立尾战斗机分块结构网格中求解Euler方程,获得了合理的全机气动力系数、翼面压力分布和绕流图谱结果.