1971年，CUNDALL提出的一种处理非连续介质问题的数值模拟方法，离散元
方法（Discrete Element Method,简称DEM），理论基础是结合不同本构关系（应力-应变关系）的牛顿第二定律。随后，这种方法被越来越广泛的应用于涉及颗粒系统地各个领域。通过求解系统中每个颗粒的受力（碰撞力及场力），不断地更新位置和速度信息，从而描述整个颗粒系统。
CFD-DEM耦合方法的基本思路是：通过CFD技术求解流场，使用DEM
方法计算颗粒系统的运动受力情况，二者以一定的模型进行质量、动量和能量等的传递，实现耦合。
该方法的优势在于，无论流体还是颗粒，都可以采用更适合自身特点的数值方法进行模拟，将颗粒的形状、材料属性、粒径分布等都考虑进来，更准确地描述颗粒的运动情况及其与流场的相互影响。
大多为气固（Gas-Solid Fluid）两相流(PFC)
1. 北京海基.流体－颗粒系统数值模拟的 FLUENT－EDEM解决方案, 2009 ２．Coupled particle and fluid flow modeling of fracture and slurry injection in
weakly consolidated granular media
We have investigated the use of coupled fluid flow and particle flow models to simulate fracture and dilation processes during waste injection.
３． GRANULAR ＦＬＯＷＳ ＡＮＤ ＮＵＭＥＲＩＣＡＬ ＭＯＤＥＬＬＩＮＧ ＯＦ ＬＡＮＤＳ
ＬＩＤＥＳ
４． Physically Representative Micro-Mechanical Models of Fluid-Filled
Granular Media
５． Ｍｏｃｒｏ－ｍｅｃｈａｎｉｃａｌ Ｓｉｍｕｌａｔｉｏｎ ｏｆ Ｗａｔｅｒ Ｆｌｏｗｔｈ
ｒｏｕｇｈ Ｇｒａｎｕｌａｒ Ｓｏｉｌｓ
６． DEM simulation of impact force exerted by granular flow on rigid structures
（Springer-Verlag 2011）
7. Yasuhiko Okada， Shear behaviour in numerical triaxial compression tests by 3D fluid-coupled DEM: a fundamental study on mechanisms of landslide initiation（2011）
2.郑建祥. 粘附性颗粒动理学及气固两相流体动力特性的研究. 工学博士学位论文,2008
模拟方法：大涡模拟（LES）
基于 Chapman 和 Cowling 的稠密气体动理学方法，建立了粘附性颗粒动理 学理论，提出了颗粒相本构方程
3.Tatemoto 运用DEM方法数值模拟粘附性颗粒在振动流化床内流动
2011 12.29（博士论文阅读采取点）
振动压实对道路材料空间组构及其力学性能演化的离散元模拟
1. 颗粒之间只是在接触时，才有力的作用，颗粒之间存在压应力和剪应力，
但其规律性比固态物质复杂的多，具有对边界面产生压力的性质，难以抵抗拉力，抗剪强度取决于围压的大小，无围压时，理想的颗粒系统无抗剪强度，具有压硬性、剪胀性。
2. 单个颗粒的尺寸、形状、物理性质以及颗粒分布决定了固体骨架的细观力
学性质。固体、液体、气体三个组成部分之间的比例关系和相互作用决定了道路材料的宏观力学性质。由试验分析确定各种材料的力学计算参数，只能根据试验试件整体表现出的力学性能参数为依据，难以考虑材料自身微观的差异，同时试验误差引起的数据离散程度大，导致材料参数的可靠度下降。
3. 一般是由于颗粒系统自身具有的特点(离散性，耗散性，形态的多样性等)
而造成，使得研究者们不能用传统的理论揭示颗粒系统的普遍规律。
4. 这是困为离散元理论的计算单元和我们研究的真实材料之存在很多而又
难于控制的因素，这些因素之间的非线性是非常显著的，而且相互之间的影响又是很明显的。比如与材料强度有关的细观参数就与颗粒之刚的接触模量、平行粘结模量、法向平行粘结强度和切向平行粘结强度有关，如果同时变化多个细观参数柬得lB材料宏观强度的变化规律是很难控制，而且是不太可能的。
5. 主要基于混台物多孔介质理论建立三相体数值模型。采用杨松岩教授提出
的“多孔介质的车构描述”的理论为离散元时步差分控制简化方程唧J，研究湿颗粒(考虑固、液、气三项体)振动压实下的瞬态响应
6. 由于数值模拟软件将基本颗粒假定为刚性体，在模拟中水被简化为流体域
上的体力
7. 颗粒与颗粒间的动摩擦转化为颗粒与水的动摩擦
8. 由固体颗粒和填隙流体形成的一种多孔介质材料
气固两相流技术发展状况
2012.2.12
复杂两相流动中颗粒碰撞的ＤＥＭ—LES/DNS耦合模拟研究
一、
a) 欧拉——欧拉方法（Eulerian- Eulerian approach）,即将连续相和颗粒相都看作可进入的连续流质和拟连续流质，用连续介质力学的连续性方程、动量方程和能量方程来分别描述两相的性质和运动。
b) 欧拉——拉格朗日方法（Eulerian- Lagrangian approach）,即对离散相采用跟踪计算颗粒的运动轨迹，用直接模拟或者随机模型的方法计算颗粒-颗粒相互作用等。
这两种途径都需要考虑流体-颗粒两相的相互作用，欧拉-欧拉方法最典型的即使双流体模型（Two-Fluid Model，TFM），而欧拉-拉格朗日方法典型的即是离散颗粒模型（Discrete Particle Model，DPM）。此外，对分散相的描述又分随机性方法（stachastic way）和确定性方法。
二、
ＤＥＭ方法在颗粒体系混合现象的理论研究和实际应用进行了综述和讨论，包括连续模型和离散模型、DEM的验证、非球形因素和粘结性因素的考虑等，同时他们指出了DEM模型在实际应用上受计算
时间和颗粒数目的限制，在解决此限制方面，可采取的途径将CFD计算和DEM耦合起来，如采用有限元方法计算群体速度，同时考虑每个组分的浓度导出积分型输运方程等。采用两维的硬圆碟无弹性碰撞模型模拟了稳态重力驱动颗粒流。
在计算机技术和离散颗粒模型的快速发展的基础上，颗粒尺度的颗粒-颗粒、颗粒-流体、颗粒-壁面之间的相互作用的研究工作在近期得到了快速发展。
研究发现当恢复系数为0.7和0.9、颗粒摩擦系数为0.3和0.6、
壁面摩擦系数为0和0.3时模拟结果和放射性颗粒跟踪实验测量结果非常接近。
对于颗粒流的数值研究，也可分为拟连续介质描述和拉格朗日
离散颗粒描述两大类，其中连续介质描述用于弹性颗粒快速颗粒流（rapid granular flow），因为此时颗粒流行为与流体流动行为十分相似。
由于ＤＥＭ方法是基于颗粒尺度行为的基本物理定律和材料物理特性的，碾磨过程中的功耗、转矩、装载、速度等重要参数以及碰击能、力传输链、壁面应力等关键力学量都可通过DEM模拟来预测和评估。
这两天你突击把看过的考虑粘性和流体作用的文献和我上面给你提的建议（考虑粘性或流体作用的离散单元法方面的文献(尤其是看是否有用PFC的)，主要关注粘性或流体作用如何加入到离散元模型中以及做的哪些具体的问题等）的文献总结整理一下，这周五我们讨论一下。
考虑粘性和流体作用的散体介质静、动力学行为分析
粘性 viscosity,
作用 1. 动词 affect 2. 名词 [影响] effect 3. 名词 [活动] action
散体物质：Discrete Material
静、动力学行为：the behavior of still and dynamics Force;
散体物质：属于软物质范畴，是自然界和工程界最为普遍的一种物质类型，具有多学科内涵和丰富的力学行为以及广泛的应用背景。
散体物质的多尺度特性更多的体现在不同层面上具有明显不同的物理本质和率相关性方面。（a）微观尺度（特质颗粒层面）上颗粒间的挤压与摩擦在细管尺度上形成分布不均匀且可以承受外部载荷的力链。（b）宏观尺度上又可能表现出变形或流动等特征。
散体物质力学行为不仅受到其颗粒的形状、尺寸、孔隙度、颗粒表面的粗糙度、颗粒的排列方式和接触形态等内部因素的影响，也受到加载条件和环境等外部因素的影响，包含了多个事件和空间尺度的相互耦合。
对于散体物质这种由大量颗粒组成的复杂系统，仅仅通过求解单一颗粒的运动方程来实现对宏观系统的定量描述是不现实的，因此我们必须探讨新的范式来包括跨尺度耦合的过程。统计处理和分析可能是不可或缺的，这是应为系统中的大量颗粒会表现出相应的统计规律，而且由此得到的统计规律可能是问题解决过程中的一个重要输入量。
2012.2.19 滑坡：landslide
泥石流：Debris flow
润滑的；喝醉了的：Lubricated 隔离，分离；种族隔离：Segregation