采空区自然发火是井下重大灾害之一，也是煤矿安全生产领域重点研究的问题。本书在分析煤自燃特性参数及采空区内产热与散热规律的基础上，利用数值模拟技术提出了采空区自燃发火的能量迁移理论，阐述并揭示了采空区自然发火的多场耦合致灾过程和机理，建立了三维多场耦合数学模型；确定了四面体单元控制体的合理圈划方案，提出了一种具有较高精度的有限体积法的新算法；利用新算法离散自然发火模型，编制了“采空区自然发火三维仿真系统”，结合现场参数，模拟得到了注氮前后及停采状态下的采空区各场分布情况，分析了推进速度、遗煤厚度及工作面风量等对自然发火的影响；提出了一种以工作面上隅角涌出量来定量评估采空区自然发火风险的新方法，提高了采空区自然发火早期预报的准确性，还自主研发了一种新型干冰快速升华装置，即干冰相变发生器，成功处置了井下的早期自然发火。

前言
第1章 绪论
1.1 煤矿煤自燃灾害现状
1.2 国内外研究现状
1.2.1 煤自燃的相关理论研究
1.2.2 煤自燃过程的实验研究
1.2.3 采空区自然发火数值模拟研究
1.2.4 研究发展趋势
第2章 煤低温氧化实验方法
2.1 煤低温氧化实验系统
2.2 煤自燃的特性参数
2.2.1 标准耗氧速率
2.2.2 标准CO生成速率
2.2.3 标准放热强度
2.3 遗煤粒度分布特征
2.4 氧气浓度的影响
2.4.1 煤样制备及实验过程
2.4.2 实验结果及分析
2.4.3 讨论及验证
2.4.4 采空区耗氧速率与c0生成速率
2.5 粒度的影响
2.5.1 煤样及实验过程
2.5.2 实验结果及分析
2.5.3 混合煤样耗氧速率
2.5.4 讨论
2.6 挥发分的影响
2.6.1 煤样及实验过程
2.6.2 实验结果及分析
2.6.3 讨论
2.7 本章小结
第3章 采空区自然发火多场耦合致灾机理
3.1 采空区自然发火的能量迁移理论
3.1.1 采空区移动坐标系
3.1.2 采空区自然发火的能量迁移理论
3.1.3 采空区最高温度预判方程
3.2 采空区自然发火的多场耦合机理
3.2.1 采空区中“场”的概念
3.2.2 采空区内的产热与散热规律
3.2.3 采空区多场耦合作用的机理及过程
3.3 采空区自然发火的多场耦合数学模型
3.3.1 方向导数与梯度
3.3.2 采空区流场数学模型
3.3.3 采空区渗流参数
3.3.4 采空区氧气浓度场模型
3.3.5 采空区冒落煤岩固体温度场模型
3.3.6 采空区气体温度场数学模型
3.3.7 采空区自然发火的多场耦合三维模型
3.4 本章小结
第4章 三维模型离散方法
4.1 导热微分方程
4.1.1 控制方程
4.1.2 边界条件
4.2 有限单元法离散
4.2.1 导热的有限单元方程
4.2.2 插值函数
4.2.3 第一类边界及内部单元的离散化
4.2.4 第二类边界单元的离散化
4.2.5 第三类边界单元的离散化
4.2.6 有向曲面的热通量
4.2.7 导热积分方程的推导
4.2.8 边界条件
4.3 有限体积法离散
4.3.1 四面体网格节点的控制体
4.3.2 导热的有限体积方程
4.3.3 插值函数
4.3.4 第一类边界及内部单元控制体的离散化
4.3.5 第二类边界单元控制体的离散化
4.3.6 第三类边界单元控制体的离散化
4.4 有限单元法与有限体积法对比
4.5 本章小结
第5章 自然发火模型的离散及求解
5.1 模型中的参数离散处理
5.1.1 采空区气体密度
5.1.2 采空区孔隙率、渗透率
5.2 采空区流场模型的离散
5.2.1 控制体与插值函数选取
5.2.2 采空区流场方程离散
5.2.3 流场边界处理
5.2.4 采空区速度方程的离散
5.3 采空区氧气浓度场模型的离散
5.3.1 控制体与插值函数选取
5.3.2 采空区氧气浓度场方程离散
5.3.3 氧气浓度场边界处理
5.4 采空区温度场模型离散
5.4.1 控制体与插值函数选取
5.4.2 冒落煤岩固体温度场离散
5.4.3 气体温度场离散
5.4.4 温度场边界处理
5.5 采空区解算区域及网格划分
5.5.1 解算区域
5.5.2 网格划分
5.5.3 节点编号及坐标
5.6 线性方程组的求解
5.6.1 节点线性方程的总体合成
5.6.2 系数矩阵的压缩与存储
5.6.3 系数矩阵的求解方法
5.7 程序设计
5.8 本章小结
第6章