土壤下渗问题的格子玻尔兹曼模拟_张东辉.pdfVIP

第38 卷第6期 河 海 大 学 学 报 (自然 科 学 版 ) Vol.38 No.6 2010年11月 Journal of Hohai University(Natural Sciences) Nov.2010 DOI:10.3876/j.issn.1000 1980.2010.06.014 土壤下渗问题的格子玻尔兹曼模拟 1 2 1 1 张东辉 , 芮孝芳, 马哲树 ,孔祥雷 (1.江苏科技大学船舶与海洋工程学院, 江苏 镇江　212003;2.河海大学水文水资源学院, 江苏 南京　210098) 摘要:将一种新的格子玻尔兹曼模型(简称LB 模型)应用于土壤水流下渗过程的探讨.在恰当的时 间和空间多尺度化方案基础上, 给出了Richards 下渗方程的LB 模型所对应的宏观量和平衡态分布 函数形式.通过对扩散方程和线性Richards方程的分析, LB 模型的模拟结果与分析解相吻合, 并详 细探讨了弛豫系数、网格步长和时间步长等参数对计算误差的影响.与Philips 解的比较表明, 该LB 模型可成功应用于非线性Richards 下渗方程的求解, 并在计算稳定性和处理非线性等方面展现出 很好的优点. 关键词:格子玻尔兹曼模型;土壤下渗;理查兹方程;非线性 + 中图分类号:S152.7 2　　　文献标志码:　　　文章编号:1000 1980(2010)06 0671 06 应用连续介质模型来分析流体流动问题, 实际上是通过分析微元体的质量守恒、动量定理、能量守恒推 得N-S 方程组, 然后应用有限差分法等方法对微分方程进行离散, 得到各结点的代数方程组进行求解.这类 计算方法比较直观, 目前应用非常广泛.但如果微分方程是强非线性的, 数值稳定性便成为很大的问题, 往往 [1] 需要通过变换和设计特殊的数值格式来分析具体问题 . 与上述思路不同, 研究流体的流动也可从微观离 散模型出发, 即从更微观的尺度———大量分子运动的 [2] 主控方程出发, 向上“集成”来探讨流体流动特性 , 如 [2] 图1 所示.统计热力学中的玻尔兹曼方程 是描述系 统中大量粒子分布函数f (x, v , t )的演化方程, 它是联 系微观与宏观的桥梁, 其简化形式为 f f F f ＊ t +e x +m e = Ψ(f ,f ) (1) 图1　两种推导宏观微分方程的途径 式中:e ———微观粒子的速度, e=(e , e , e );x ———粒 x y z Fig.1　Two kinds of processes deriving 子的位置矢量, x =(x, y , z );F ———作用在分子上的 macroscopic differential equation 外力(也称源项);Ψ(f ＊,f )———粒子与粒子之间的碰撞散射项, 表示粒子分布函数f (x , e, t )由于粒子碰撞 而发生的改变率.玻尔兹曼方程的解是不同时间的粒子分布函数f (x , e, t ).一旦求得f (x , e, t )的具体表达 式, 可通过它与质量、速度、含水率等宏观量之间的关系, 来推得后者的变化特点. [3] 格子玻尔兹曼方法采用网格离散的方法求解玻尔兹曼方程 .如果不考虑外力项的作用, 式(1)在x 处、α方向上可离散为[2]