多孔介质内溶解与沉淀过程地格子Boltzmann方法模拟.PDF

物理学报 Acta Phys. Sin. Vol. 64, No. 15 (2015) 154701 多孔介质内溶解与沉淀过程的格子 Boltzmann方法模拟 张婷 施保昌 柴振华 1)(理学院, 武汉科技大学, 武汉 430081) 2)(数学与统计学院, 华中科技大学, 武汉 430074) ( 2014 年10 月26 日收到; 2015 年2 月28 日收到修改稿) 本文采用格子Boltzmann 方法模拟了多孔介质内的溶解和沉淀现象, 并分析了雷诺数、施密特数、达姆科 勒数对多孔介质孔隙结构及浓度分布的影响. 结果表明: 对于多孔介质内的溶解(沉淀) 过程, 当雷诺数越大 时, 孔隙率越大(小), 平均浓度值越小(大); 当达姆科勒数或施密特数较小时, 溶解和沉淀过程均受反应控制, 此时反应在多孔介质的固体表面较为均匀的发生; 当达姆科勒数或施密特数较大时, 溶解和沉淀过程均受扩 散控制, 此时反应主要发生在上游及大孔隙区域. 关键词: 多孔介质, 溶解/沉淀, 格子Boltzmann 方法 PACS: 47.11.Qr, 47.56.+r DOI: 10.7498/aps.64.154701 来模拟存在对流的反应扩散系统. 在这些相对较 1 引 言 少的采用LGA 或者LBM 来模拟化学反应的文献 中, 大部分都关注于对流体与流体之间发生的化 多孔介质内的溶解和沉淀现象广泛存在于能 学反应进行研究78 , 很少有考虑流体- 固体之间的 源、环境和化工等领域, 如石油的开采过程中、二 化学反应问题. Wells 和Janecky 等人采用LGA 模 氧化碳的捕集和地下埋存以及地下水资源污染等 拟表面反应(包括溶解和沉淀), 其基本思想是, 当 等均涉及到这类问题. 随着计算机硬件和软件技术 溶解发生时, 固体结点变为流体结点; 当沉淀发生 的迅速发展, 以数值方法为基础, 运用计算机进行 9 数值计算的方法成为了研究这类对流反应扩散问 时, 流体结点变为固体结点. Keleman 等 采用类 题的重要手段. 近年来发展起来的格子Boltzmann 似的观点用格子Boltzmann 方法(LBM) 也模拟了 方法 (lattice Boltzmann method, LBM), 由于其微 溶解现象. 2000 年, He 等10 将流体-壁面间的表 观本质和介观特点, 为研究这类复杂问题提供了可 面反应与扩散机理相耦合, 明确地将化学反应与格 行的框架, 且已被广泛应用于各个领域13 . Wells 子Boltzmann 方法结合起来, 其中表面化学反应的 等人早在1991 年就用lattice gas automata (LGA) 发生主要是通过边界条件来体现. 2002 年, Kang 模拟了矿物表面和孔隙网络内溶质的传输和化学 等11 基于He 的模型建立了一类格子Boltzmann 反应4 . Janecky 5 也采用类似的方法模拟了地球 模型, 用来模拟多孔介质内两种水合物(分别为反