多孔介质振荡流格子Boltzmann模拟.pdfVIP

‘ 中国工程热物理学会 传热传质学 2007年学术会议 编号；073378 I t 多孔介质振荡流格子一Bozman．n模拟 孟繁孔李志信 热科学与动力工程教育部重点实验室，清华大学航天航空学院，北京100084 联系电话：(oloEmail：ii圣丛@鱼塾g鱼丛垒：璺鱼坠：鲤 摘要多孔介质内振荡流及其换热是脉管制冷机和热声热机回热器内的基本物理过程，研究 振荡流特点对提高整机性能的优化设计具有重要意义。应用格子-Bol乞，埘a助方法对多孔介质内 的振荡流进行了数值研究，表明多孔介质振荡流的最大压差随最大雷诺数的增加而迅速增加， 在所计算的范围内，最大阻力系数与最大雷诺数呈反比例关系。在低孔隙率条件下，压差与流 速之间的相位差变化较小，孔隙率较大时，相位差随孔隙率的增加而线性降低。 关键词： 格子一Boltzm锄方法多孔介质振荡流相位差 ． I 引言 近年来，脉管制冷机和热声热机由于具有众多优点而取得迅速发展。热声热机具有 3， 无机械运动部件、以惰性气体为工质等优点，在近二十年中得到广泛的研究和关注n 进一步提高热声热机的效率是将热声热机推广到工程实际应用的关键。脉管制冷机及热 声热机的回热器一般为有填料的多孔介质结构，多孔介质内的振荡流及其换热是回热器 内的基本物理过程。了解回热器内的振荡流的压降及相移特点，对脉管制冷机和热声热 机的优化设计具有十分重要的作用晗’。 目前已有一系列针对多孔介质内的振荡流及其换热的研究工作。Zhao等b1实验测量 了丝网状的多孔介质振荡流的阻力特性。Ju等H1实验研究了脉管制冷机的回热器内的振 荡流过程，分析了不同因素对振荡流阻力以及压力与速度问相移的特点。Choi等旧。也采 用实验手段对制冷机回热器内振荡流进行了类似的研究。Le0119等№1实验研究了金属泡 沫结构内的振荡流，分析了铝金属泡沫内振荡流的阻力特性。目前，对复杂多孔介质内 振荡流及其换热进行数值模拟的较少，一方面，传统的CFD方法难于处理多孔介质内的 复杂几何边界，另一方面，构造出接近于实际的多孔介质也面临很多挑战。随机生成是 最简单的多孔介质构造方法，但是构造出的多孔介质微观结构特点不能控制：重构方法 是常用方法之一，但是确定修正函数关系式l卜分复杂；wang等川结合多孔介质构造的孔 隙生长模型和丛生理论，提出了可控制多孔介质微观结构的四参数生成方法(QsGs)。 该方法的构造过程更接近于很多孔介质的生成过程。 具有程序简单、易于并行、容易处理复杂边界等优点，已经被广泛的应用于流动和换热 问题的研究旧·。本文利用LBM易于处理复杂边界的优势，对多孔介质内的振荡流进行数 值模拟，分析了多孔介质振荡流的压降与相移的特点。 2格子一BoItzmann模型 流动与换热等输运过程，宏观参量可以通过对分布函数的统计而得到。常用的格子结构 方向的格子速度为： 652 O． f=O (1) f=1，2，3，4 cI={c(cos(f—J)石／2，sin(f—1)万／2) l√易(cos【(f一5)万／2+石／4】，sin【o一5)万／2+万／4】)f=5，6，7，8 其中，c=厕，月为气体常数，L为流体平均温度。 ‘图l D2Q9型格子结构 对碰撞项采用BGK近似，离散后的分布函数演化方程为：