颗粒流体两相流具有非线性非平衡的特征，其中的多尺度结构对于系统的流动、传质、反应等都具有显著的影响。最近随着计算机技术的发展，离散颗粒方法因其能准确的捕捉颗粒流体两相流中的两相结构而受到越来越广泛的关注。然而现有离散化方法在流固耦合方面往往采用连续介质方法的处理方式，忽略了离散化方法给出的对应颗粒位置的非均匀信息，使得颗粒与流体在不同的尺度上进行计算，影响了计算的准确度;而对于流体计算方面，传统TFM方法受到自身并行性能的限制，计算速度较难满足大规模的计算需求，需要寻求具有更佳并行性能的方法。同时，虽然介尺度结构及其相关能量关系对于系统有显著的影响，但是目前局部区域上基于介尺度结构的能量关系仍然不明，系统局部所遵循的稳定性条件仍然欠缺。　　鉴于此种情况，为了实现高效、准确的颗粒流体两相流模拟这一长远目标，本文建立了应用于离散化方法的基于结构的曳力模型，结合并改进了MP-PIC方法，对于流体相计算建立了基于LBM算法的并行效率优异的TF-LBM方法，并对于系统能量关系进行了初步的尝试与探索。　　论文第二章首先在开源软件MFIX的基础之上实现了MP-PIC的程序计算，为了验证程序，对于一些鼓泡床和循环流化床提升管算例进行了模拟，结果显示与实验结果吻合，但采用传统均匀曳力会产生较大的偏差。而直接采用TFM层次的EMMS曳力如前所述会影响计算的准确度。　　为解决这一准确性问题，论文第三章提出了面向DPM的EMMS曳力模型，模型引入了颗粒位置所代表的非均匀信息，并因此将控制单元内气固相划分为一个稀相及多个密相，采用压降平衡条件对各相中流体进行分配，以此回归拟合，得到改进的曳力系数。结果显示，相比于均匀的曳力系数，面向DPM的EMMS曳力模型的计算结果准确度有了显著的改善。　　为实现快速模拟的目标，论文第四章提出了TF-LBM算法，它能够在保留LBM并行性能的前提下，在TFM的尺度上对于流体相进行计算，从而可望显著加快流体相的计算速度。通过一些简单的算例，对于本方法也进行了初步的验证。　　最后，为探索介尺度结构对系统的影响，论文第五章对于气固流态化系统中的能量关系进行了分析。本工作加深了对于系统能耗的组成和EMMS模型中各变量的物理意义及关系的理解，通过具体的分析讨论了EMMS稳定性条件不能在局部成立的可能原因，对于未来建立局部稳定性条件与介尺度结构的关联提供了一定的基础。　　论文最后总结了所得到的主要结论，并对于模型的前景以及进一步研究的方向进行了展望。