土体饱和渗流的lattice boltzmann method数值模拟研究

《土体饱和渗流的lattice boltzmann method数值模拟研究》由会员上传分享，免费在线阅读，更多相关内容在学术论文-天天文库。

1、长沙理工大学硕士学位论文土体饱和渗流的LatticeBoltzmannMethod数值模拟研究姓名：樊火申请学位级别：硕士专业：岩土工程指导教师：王桂尧;蒋中明20090501摘要二十世纪八十年代发展起来的LatticeBoltzmannMethod(LBM)从最初的计算流体力学fCFD)应用领域逐渐向其他诸多领域渗透，己成为一个国际研究热点。LBM基于微观离散模型和某些细观动力学方程，以离散的速度集、自底向上(Down．top)建立模拟系统，兼具Euler法和Langrange法的特征。与传统

2、的数值模拟方法相比，它具有算法简单、编程容易、易于处理复杂边界条件、天然适应大规模并行计算等优点，这也使得它非常适合研究渗流问题。本文以LBM为切入点，在理解应用LBM解决渗流问题的关键在于建立适合渗流运动的演化方程后，对常规LBM演化方程和局部平衡态分布函数进行修正，即在常规LBM演化方程的右端增加多孔介质对流体的牵制作用力项(该力项可以写成离散速度集的幂级数形式)，并将多孔介质的孔隙率引入局部平衡态分布函数；接着确定了该力项的具体表达式和其合理的离散形式。奠定了用LBM解决渗流问题的理论基础

3、。进而，在消化和吸收近年来LBM理论和应用等方面相关成果的基础上，初步开发了利用LBM理论模拟土体饱和渗流的软件平台(LBMSeepage)。该软件平台的前处理模块既具有较为完备的绘图功能，又提供了与其他图形文件交互的接口，能方便地完成渗流场几何剖面的绘制；该软件平台的核心模块的核心部分采用目前LBM中应用最普遍的LBGK模型且边界处理采用非平衡外推格式，既加快了收敛速度，又能使演化方程在时问和空问上逼近二阶精度；该软件平台的后处理模块采用应用广泛的数据可视化软件包，极大地扩展了LBMSeepa

4、ge的功能。为验证LBMSeepage的软件平台的正确性，首先，用该软件验证了Darcy定律，进而用该软件在表征单元体(REV)尺度上模拟多种情况下的地下水运动和成层土样的竖向渗流运动，并与室内实验结果进行对比验证；此外，还用该软件在孔隙(Pore)尺度上对水在人造粗糙裂隙中的流动规律和变孔隙率的多孔介质渗流特性进行了模拟，并与室内实验结果进行对比验证。研究表明该软件平台具有一定的灵活性、有效性和可扩展性，有望成为研究渗流运动的有力工具。关键词：渗流；多孔介质；LBM；数值模拟；软件平台ABST

5、RACTTheLatticeBoltzmannMethod(LBM)developedin1980shasgraduallyinfiltratedfromtheinitialCFDapplicationstootherfieldsandhasbecameoneoftheinternationalresearchfocuses．LBMbasedonthemicroscopicdiscretemodelandsomemicro-dynamicequations，andthesimulationsys

6、temisestablishedusingdiscretespeedsetsinthewayofDown-top，whichhasthecharacteristicsofEulermethodandLangrangemethod．WiththetraditionalCFDmethods，LBMhassomeadvantages，suchassimplealgorithm，easyprogramming，easytohandlecomplexboundaryconditions，adaptnatu

7、rallytolarge—scaleparallelcomputingandSOon，whichmakesitveryappropriatetostudytheseepageproblems．Inthispaper，LBMasthebreakthroughpoint，aRerunderstandingthekeyofusingLBMtosolveseepageproblemsliesinestablishingsuitablefortheevolutionequationofseepagemov

8、ement，theevolutionequationsandlocalequilibriumdistributionfunctionofconventionalLBMfortheamendment，thatiS，theconstraintforce，porousmediaimposedtofluid，isaddedintherightsideoftheevolutionequationofconventionalLBM(theconstraintforcecanbewritteninpowers