基于格子Boltzmann方法的非线性渗流研究

论文摘要

非线性渗流或非达西渗流作为渗流力学中的一个重要分支，由于其在农业、能源、冶金、化工、材料、环境、生命科学以及医学等众多领域的广泛应用而引起国内外学者的普遍关注，并发展成为国际上研究的一个热点问题。通过对大量相关文献的分析可以看出，非线性渗流问题自身的复杂性大大限制了人们对其背后物理机理的认识，以至于对于某些问题至今也未能给出较为合理准确的数学模型或物理公式。目前对非线性渗流的研究主要基于三类方法：一是实验方法，作为一种传统的研究方法，其在渗流力学中仍然发挥着重要的作用，但是，在实际中该方法又受到很多条件的限制，比如实验周期长、花费大、受实验环境的影响等。二是解析或近似方法，这类方法通常会对问题先进行一定的简化或给出一些合理的假设，在此基础上，进一步通过数学的手段来求得非线性渗流的解析解或近似解。目前，这类方法中最具代表的就是近年来发展起来的基于分形理论的解析方法。另外一类就是数值方法，随着计算机技术的飞速发展，该方法越来越受到人们的广泛关注，并已成为解决实际的问题的一种重要手段。与实验方法相比，数值方法依赖于前期建立的数学模型，但它具有周期短、耗费小、不易受环境影响等特点；同时，该方法也可以弥补解析或近似方法所存在的不足。目前对非线性渗流的数值研究主要是借助于传统的数值方法（如有限差分法，有限体积法，有限元方法等）来数值求解孔隙尺度上的Navier-Stokes方程或表征体元（REV）尺度体积平均的Navier-Stokes方程。作为一种研究非线性渗流问题的有效手段，数值研究近年来备受关注，其应用范围也越来越广泛，但多孔介质的复杂结构使得这些方法在研究渗流问题时会遇到如下困难：（1）边界处理复杂，破坏数值方法自身的稳定性；（2）多孔介质中的微尺度效应难以体现；（3）计算量较大，并行效率低。鉴于此，发展高效的数值方法，并利用其研究非线性渗流规律，建立完善的数学描述是揭示非线性物理机理的理论基础。格子Boltzmann方法（Lattice Boltzmann Method，LBM）作为一种新兴的微观方法，它最初起源于格子气自动机（Lattice Gas Automata，LGA），但后来也被证实可以从连续的Boltzmann方程经过差分离散得到。格子Boltzmann方法与以连续微分方程为基础的宏观计算流体力学方法有着本质的不同，它是基于流体微观模型和细观动理论方程的方法，因此不受连续假设的限制，这为其能够模拟多孔介质中微细管道中的流动、刻画滑脱效应（Klinkenberg效应）提供了坚实的理论基础。此外，与传统的计算流体力学方法（基于Navier-Stokes方程的数值方法，比如，有限差分法、有限体积法、有限元方法等）相比，LBM还具有许多独特的优势，如计算效率高、边界条件容易实现、具有天然并行性等，更为重要的是其自身的微观特性使得它成为研究多孔介质中非连续流动提供了一种新的有效手段。本文首先发展了LBM的相关理论，在此基础上，利用LBM对单相非线性渗流力学中的一些前沿问题进行了探讨。论文的主要工作包括以下两个方面：（1）在格子Boltzmann方法的相关理论方面，论文首先提出了一种处理微尺度流动的反弹与Maxwell漫反射组合边界条件，并详细分析了该边界条件在单松弛模型与多松弛模型中的离散效应；其次，论文还给出了一种求解Poisson方程的新的格子Boltzmann模型，该模型可以有效克服已有模型的不足。这些理论的发展将为推动LBM在渗流领域的应用奠定必要的基础。（2）在非线性渗流的数值研究方面，论文首先系统研究了高速非线性渗流，详细分析了低雷诺数下产生非线性现象的物理机理，并进一步并给出了改进的数学模型；其次，基于（1）中的理论结果，论文详细研究了气体通过多孔介质的滑脱效应，利用文中发展的一些理论结果揭示了滑脱效应产生的物理根源；最后，论文还进一步考察了微尺度多孔材料中的电渗流，系统分析了各种因素对电渗流速度分布的影响，这些结果对实际中的工程问题具有重要的指导意义。总之，本文利用格子Boltzmann方法研究了孔隙尺度和REV尺度的非线性渗流，为推动该方法在渗流领域中的应用作出了许多有益的尝试。同时，针对针对高度非线性的Poisson-Boltzmann方程和复杂的多孔介质，我们分别构造了新的模型和边界条件；此外，我们还对格子Boltzmann方法进行了大量的数值实验，验证了该方法用于模拟渗流问题的可靠性。这些工作为后续研究的开展奠定了必要的基础。

论文目录

• 摘要
• Abstract
• 目录
• 1 绪论
• 1.1 研究背景
• 1.2 格子Boltzmann方法及其发展历史
• 1.3 格子Boltzmann方法的应用
• 1.4 研究内容
• 1.5 全文安排
• 2 格子Boltzmann方法的基本原理、模型和边界处理
• 2.1 Boltzmann方程、平衡态分布函数和BGK近似
• 2.2 从连续的Boltzmann方程到离散的格子Boltzmann方程
• 2.3 格子Boltzmann方程的一些基本模型
• 2.4 格子Boltzmann方程还原到Navier-Stokes方程的推导过程
• 2.5 边界条件的处理格式
• 2.6 本章小结
• 3 高速非线性渗流
• 3.1 研究背景
• 3.2 描述多孔介质流动的数学方法
• 3.3 数值方法与边界处理
• 3.4 数值结果与相关讨论
• 3.5 本章小结
• 4 微尺度单管内的气体流动及粗糙壁面的影响
• 4.1 研究背景
• 4.2 微尺度流动的格子Boltzmann方法
• 4.3 数值结果与相关讨论
• 4.4 本章小结
• 5 气体通过多孔介质的的滑脱效应
• 5.1 研究背景
• 5.2 相关理论及其分析
• 5.3 格子Boltzmann方法与边界处理
• 5.4 数值结果与相关讨论
• 5.5 本章小结
• 6 微尺度单管内的电渗流
• 6.1 研究背景
• 6.2 微尺度电渗流的研究现状
• 6.3 宏观控制方程
• 6.4 格子Boltzmann方法
• 6.5 数值结果与讨论
• 6.6 本章小结
• 7 表征体元尺度的非线性渗流
• 7.1 微尺度多孔介质中的电渗流
• 7.2 多孔介质方腔内混合对流
• 8 总结与展望
• 8.1 全文总结
• 8.2 研究展望
• 致谢
• 参考文献
• 附录1 攻读学位期间发表的学术论文
• 附录2 攻读学位期间参加的学术会议
• 附录3 攻读学位期间参与的科研项目