多孔介质中热对流二次分岔数值的分析.pdfVIP

北京理工大学学报990102 北京理工大学学报 JOURNAL OF BEIJING INSTITUTE OF TECHNOLOGY 1999年 第19卷 第1期 No.1 Vol.19 1999 多孔介质中热对流二次分岔的数值分析 徐世英　冯长根　刘赵淼 摘　要　目的　进行了多孔介质中热对流二次分岔的数值分析。方法　采用了扩展方程方法，此方法无 需求解原动力学方程和辅助参数。结果　得到了长宽比为1～ 矩形区域内多孔介质中热对流二次分 岔点的Rayleigh数及相应的流场和温度场分布。结论　所给方法对确定二次分岔点是有效的。 关键词　多孔介质； 二次分岔； 扩充方程组； 对称性 分类号　O357.3 Numerical Analysis of Secondary Bifurcations of Thermal Convertional Flow in Porous Media Xu Shiying　Feng Changgen　Liu Zhaomiao Abstract　Aim media. Methods auxiliary parameter. Results Key words 　　多孔介质中热对流广泛存在于地下贮油、天然气开采、化学反应工程及新型材料的制作中。 近年 来，关于方腔中的多孔介质对流分岔机理的研究取得一些成果。文献［1］用摄动法给出了方腔中的多 介质中热对流的主分岔解的流场和温度场；文献［2］用3次样条及时间相关法数值模拟了热对流流态的 分岔解分支的渐近展开式。 用数值模拟方法可以观察到解的变化，但计算量大，特别是由于给出的初值 不同而很难确定所得结果是哪一个分岔解，也难以确定分岔解是一次的还是二次的，即很难分析解的结 构。 作者构造了一种数值分析方法直接确定多孔介质中热对流二次分岔点和相对应的模态。 1　物理问题的数学描述 　　将等厚度的多孔介质水平安放，在下面缓慢加热，若上下表面温差很小，则多孔介质处于静止传热 过程，若上下表面温差很大，或者说Ra达到一定数值时，则多孔介质内流体出现回流。随着Ra的增加， 回流的数目和流态发生改变。设流场为边长1×1的矩形， X轴和Y 轴分别和矩形的下边界和左边界重合， 沿Z 的条件下，其量纲为1的动力学方程是 file:///E|/qk/bjlgdxxb/bjlg99/bjlg9901/990102.htm（第 1／5 页）2010-3-22 15:29:26 北京理工大学学报990102 　　　　　　　　　 　　　　　　(1) 式中　Ψ为流函数； θ为温度。边界条件为 　　　　　　　　　　　　　 　　　　　　　　　　　(2) 　　引进坐标变换　　　　　　　　　　　　X= lx,Y = y， l为变换后的长宽比；u为温度差函数。使得 　　　　　　　　　　　　　　　　　　　　 θ= θ+ u，　　　　　　　　　　　　　　　　(3) 0 式中　 θ=1-y相应于方程组(1)静态传热解。考虑方程(1)(2)驻态解得 0 　　　　　　　　　　 　　　　　　　　　(4) 边界条件：整个边界满足 Ψ=0； 并且y=0,1;u=0;x=0,1;u/x=0。 2　确定二次分岔点及其模态的扩充方程 　　在矩形区域上，取x轴和y轴步长为h=1/(2n+1)，n为正整数。作两族与坐标轴平行的直线 　　　　　　　　x= ih,　y= jh,　i=0,1,2,…,2n+1,　j=0,1,2,…,2n+1。 采用5点中心差分格式并考虑边界条件对方程(4)