开题报告DEMLBM耦合算法的实现及其在岩土流固合问题上的应用.doc

攻读硕士学位论文选题报告 DEM-LBM耦合算法的实现在岩土流固耦合问题上的应用2014年5月22日 1 研究背景 耦合历来是一个存在广泛的。领域、石油开采领域都有着典型的流固耦合问题。在岩土工程中，的耦合有滑坡涌浪问题、渗透破坏、问题等而和固体相互作用的现象更是出现在各种各样的问题上：地下水开采中会遇到由于地下水的渗流和地层的相互作用而引起的地面沉降问题；水利工程中遇到库区水与大坝渗流发生的稳定渗流非稳定渗流、管涌破坏等关系大坝安全的基坑开挖的过程中，会遇到地下水基坑稳定安全问题；深埋管道也会遇到地下水的问题；在自然界中，也存在着滑坡、泥石流-固体耦合作用的灾害问题。所以对流固耦合问题进行试验研究，探究颗粒与流体的介观作用机理是具有意义的。2 研究现状 1 传统的流固耦合研究方法在，达西就提出了在渗流的水的流速，即的达西定律土体中的平均渗透速度与土体的渗透系数以及水力梯度算作岩土工程，流固问题研究的开端。 在，俄国的茹科夫斯基提出了渗流微分方程土体视为，考虑达西定律水中的渗流特性。 10年，理查森提出了差分方法求解渗流微分方程土体中水的渗流特性。 ，太沙基在岩土固结问题上提出了单向固结理论，考虑消散的，得到土体随的关系。 的流固耦合算法在工程具有一定的性，但是本质都是单相的。只考虑了两相中一相另一相的影响动态变化，而另一则产生动态。在流床中投入少量颗粒，颗粒流场剧烈的影响，而流场不会受到颗粒的影响。方法在流固耦合问题中的应用计算机能力的发展，两相耦合问题成为可能流固耦合算法已在有限差分法、有限单元中得到了。差分法和有限元法本质上是将固体和流体都是为连续介质进行计算。对于土体这种来说，的假设有一定的不适用性。 Cundall提出了离散元的方法以来，离散元在工程方面的应用得到了大力的发展在流固耦合方面，也有学者利用离散元与计算流体力学软件进行了研究。 （）FlUENT和IDEM耦合来进行渗透耦合理论研究并在、开采、油藏等方面进行了应用（）编写了CFD计算软件并与已有的DEM软件进行耦合进行了方腔流动等测试问题。 采用传统的CFD技术与DEM方法耦合会遇到很大的困难。原因在于强大的 总的，目前的与计算流体力学耦合是将流体视作连续介质，利用总体控制方程进行推演，计算速度上不具有优势。传统的采用数值积分方法的流体软件处理的具有复杂边界的流固变得异常困难。提出的LBM方法，是Qian等人提出了LBM方法的DnQm模型之后，了流体界的关注。LBM不同于其他的方法，直接由宏观上的Euler，由、动能和连续方程推演而出的NS控制方程。而是从分子动力学的角度出发，微观的动量、动能守恒定律通过一离散化方法，得到了方法。 方法空间、时间、都是离散的而且流体微团的运动遵循简单的碰撞定律，某处的仅取决于该处以及其临近的几个格点的状态算法使得解决了传统数值积分法模拟中的困难由于简单的原理，使得该方法在计算易于实现，且具有并行特性。但这样的简单的离散方法，却并没有丧失对于流体宏观的描述能力，能够有效宏观流场的。 总结该方法其他数值方法相比有： 2.流体的压力不需要解Possion方程直接通过流体密度和声速的状态方程便可得到 3.采用一个简单的离散速度集合来描述整个速度场的平均表象 4.能够很好的体现模型的相互作用关系,便于模拟多相流和多孔介质流 5.流动和碰撞过程是在流体的局部发生,因此便于并行处理 6.便于处理复杂边界 小结LBM方法，LBM方法是在微观离散基础上的，离散来研究的方法。该方法概念清晰，方程简单，易于编程实现，具有天然的特性，所以成为国际热点。普遍认为方法在多流方面不可比拟的优势1]。 。 流固耦合系统中，随着颗粒数量的增多，颗粒间的相互碰撞就越来越频繁，如何较好处理颗粒间的碰撞问题就越来越重要。然而，前面所提到的工作对于颗粒碰撞的描述都很粗糙并且缺乏实际的物理准确性，都受到了颗粒数量的限制，由此也限制了耦合方法的应用范围。Feng把颗粒间以及颗粒与墙之间的碰撞都采用排斥力进行处理。Niu采用琼斯势能法计算颗粒间的相互作用力。如此使得在前面的模拟中都加入了太多的人工参数，要模拟一个新问题必须首先进行反复实验，增加了过程的复杂度。此外，为了排除颗粒碰撞时的动荡现象，都设置了“安全区”，采用远程力代替颗粒间的接触力，以保证颗粒不接触。正如Yu和Xu在文章中所指出的，在此算法下，模拟流固耦合的关键主要在于如何解决固体颗粒问题。可见固体颗粒问题处理在流固耦合机理研究中的重要性。而DEM方法恰恰是处理固体颗粒问题的最好方法-格子玻尔兹曼耦合算法来外关于-格子玻尔兹曼DEM-LBM）耦合的论文发表 2007年，英国K.Han发表了DEM-LBM耦合湍流作用下不规则颗粒数值试验方法的。模拟中用到了Smagorinsky的湍流方法扩展LBM方法并且在流体