应用于非线性热传导方程的Lattice Boltzmann新方法.pdfVIP

应用于非线性热传导方程的 Lattice Boltzmann新方法 周根明 张东辉 （江苏科技大学，船舶与海洋学院，镇江， 210023 Tel:01511Email: dhzhang20@126.com） 摘要：本文通过对非线性热传导方程的探讨，给出了 Lattice Boltzmann 新方法所对应的宏观量以及平 衡态分布函数，由此还原出对应的非线性热传导方程。新的 Lattice Boltzmann 方法展现出极大的灵活 性和普适性，具有很好的应用前景。 关键词：热传导，Lattice Boltzmann 方法，多重尺度，非线性 0.引言 工程中遇到的大量传热问题, 从根本上说都是非线性的。当导热系数和温度之间存 在关联时，对应的热传导方程属于非线性方程，其解不再满足叠加性。目前求解非线性 方程采用的方法包括小扰动法、变分法、积分法和基尔霍夫变换法等等[1, 2, 3]。Lattice Boltzmann 方法是一种处理微分方程的新方法，对于热传导类型的方程，线性情况的 Lattice Boltzmann分析可见文献[4,5]，能否使此方法应用于真正的非线性热传导方程？ 目前国内外未见文献发表。从更广的角度来看，热传导方程是属于抛物型偏微分方程。 在大自然的物理过程中，扩散、渗流等问题与热传导现象存在共性之处，即“流(flux)” 是源于“势梯度”的推动，两者之间存在如下广义关联： dx dT kT dx dT TDq m== )( (1) 这里的 T 可以代表温度、浓度等变量，q 则是对应的热流、粒子流等，而 D 则是 导热系数、扩散系数等。如果 D 与物理量无关，即 m=0, 那么我们称此问题为线性问 题；反之 m≠0 则是非线性问题。对于扩散现象，非线性扩散意味着扩散系数是浓度的 函数；而对于热传导现象，则是导热系数与温度有关。 严格而言，如果边界条件也是温度的非线性函数关系，这类问题也可概括入非线性 热传导问题。但这里并不对其进行讨论。 这里所讨论的非线性热传导问题是： ))(( 2 x T TT xt T ∂ ∂ ++ ∂ ∂ = ∂ ∂ γβα (2) 0≤x≤1, 其边界条件是： 0)0,( 0),1( 1),0( = = = xT tT tT 中国工程热物理学会 传热传质学 学术会议论文 编号: 083573 *资助项目：国家自然科学基金之资助项目 这是一维非稳态导热问题, 已化为无因次形式，两端温差保持不变。如果假设 ρ=c =1, 那么导热系数与温度之间的关系是： 2TT γβαλ ++= (3) 为了应用 Lattice Boltzmann方法, 可使上述方程变换为： ) ) 32 ( ( 32 x TTT xt T ∂ ++∂ ∂ ∂ = ∂ ∂ γβ α (4) 文章中为便于比较和讨论，假设 α=1.0. 1. Lattice Boltzmann方法 Boltzmann离散方程则为下列形式： )),(),(( 1 ),(),( txftxftxftexf eqααααα τ εε −−=−++


(5) 为了应用新的 Lattice Boltzmann 方法解决上述问题，这里采用三速模型，只需考 虑前三重时间尺度。这里引入三个时间尺度： )20( −== ktt kk ε 那么时间和空间坐标的多重尺度表达为： )( )( 0 3 2 21 0 ε ε εε Ο+= Ο+++= xx tt tt 空间坐标的微分形式不变，时间坐标的微分其多重尺度的形式是：       Ο+ ∂ ∂ = ∂ ∂ + ∂ ∂ + ∂ ∂ + ∂ ∂ = ∂ ∂ )( )( 0 3 2 2 10 ε εεε xx O tttt 粒子分布函数同样可展开为： ( )5)2(2)1()0( εεε αααα Ο