Benard对流问题五模类Lorenz方程组动力学行为及数值模拟.pdfVIP

第30卷第3期 大 连 交 通 大 学 学 报 V01．30No．3 DAI．IAN 2009年6月 JOURNAl。OF JIAOTONGUNIVERSlTY Jun．2009 文章编号：1673·9590(2009)03．0096-04 动力学行为及其数值模拟 王海燕1，王贺元2 (1．大连交通大学理学院，辽宁大连116028；2．辽宁工业大学数理系，辽-7锦州121001) 摘要：考虑底部温度高于顶部温度的两平行板间流层的Benard对流问题。该问题可以用Navier-Stokes方 程与热传导方程束描述，对其速度场和温度场等相关变鼍进行二维傅立叶展开，经过复杂运算得到几个新 五模类Lorenz方程组，对这些方程组的动力学行为进行了数值模拟． 关键词：Navier-Stokes方程；热传导方程；混沌 文献标识码：A 中图分类号：0175．14；0241．81 O 引昌 T(戈，=，t)一L． 这里的咒。在两板间呈线性递减，假定流体 Lorenz方程是将满足一定边界条件的Navier． 不可压缩，则方程(1)可表示为如下形式： Stokes方程和热传导方程进行傅立叶展开后截断 而产生的¨≈J．其实际背景来源于简化的大气对 去V 流模型，两平行板问流体从下板加热，热量从下至 上传递，当温度比较低时，热量通过热传导方式传 争=一辫+管筹+珊口㈣ 递，流体不发生运动。当温度增高时，液体就要发 函数沙(算，z，t)和0(石，z，t)在边界上满足 生对流，当温度继续升高时，液体会翻转而发生混 沌．在热压力条件下不可压缩的流体运动可以被 砂=面00=0d石 (4) 描述为如下偏微分方程组： 对沙(菇，z，t)和0(髫，石，t)傅立叶展开成如下形式： (击+u‘V)M=占g△r一吉VP+秽V2u 蜘而归毳A秘)sin竿sin竿 (言十Ⅱ·V)丁=后V2丁 (1) 吣而牡孙扣)cos半sin等 V·U=0 其中，M=／／,(戈，Y，z)为流体的速度场，T=r(戈，Y， (5) z)为流体的温度场，在这个方程式中常数g，8，秒， 在具体截断过程中，将m，，m：，n，，n：赋予具体的 和k分别表示重力加速度，热膨胀系数，运动粘性 值，使得沙(髫，z，t)和0(戈，z，t)成为具体的式子， 系数和热传导系数，△丁为两板间的温度差，P为 然后将沙(戈，：，￡)和9(菇，z，f)代入式(2)、(3)、 流体的压力场，在边界上U(戈，Y，：)=0． (4)中，经过复杂计算通过待定系数得到截断方程 仅考虑竖直截面上二维流动问题，引入函数 1新五模方程的形成 沙(x，z，￡)，它的梯度为速度场，同时引入p(x，z，t) 为流体的温度场，并且在静态条件下0(菇，：，t)= 依照上述截断的方法，在截取五模类方程组 收稿日期：2008—09—08 基金项目：辽宁省教育厅高等学校科学研究计划资助项目 作者简介：王海燕(1966一)，女，副教授，学士，主要从事分歧理论及其数值分析的研究 E-mail：wanghaiyansy@djtu．edu．cn． 万方数据