Tak Moonho, Park Duhee, Park Taehyo
Research Professor e-mail:
J Eng Mater Technol. 2013 Apr;135(2):210131-2101311. doi: 10.1115/1.4023776. Epub 2013 Mar 28.
On micro scale the constitutions of porous media are effected by other constitutions, so their behaviors are very complex and it is hard to derive theoretical formulations as well as to simulate on macro scale. For decades, in order to escape this complication, the phenomenological approaches in a field of multiscale methods have been extensively researched by many material scientists and engineers. Their theoretical approaches are based on the hierarchical multiscale methods using a priori knowledge on a smaller scale; however it has a drawback that an information loss can be occurred. Recently, according to a development of the core technologies of computer, the ways of multiscale are extended to a direct multiscale approach called the concurrent multiscale method. This approach is not necessary to deal with complex mathematical formulations, but it is noted as an important factor: development of computational coupling algorithms between constitutions in a porous medium. In this work, we attempt to develop coupling algorithms in different numerical methods finite element method (FEM), smoothed particle hydrodynamics (SPH) and discrete element method (DEM). Using this coupling algorithm, fluid flow, movement of solid particle, and contact forces between solid domains are computed via proposed discrete element which is based on SPH, FEM, and DEM. In addition, a mixed FEM on continuum level and discrete element model with SPH particles on discontinuum level is introduced, and proposed coupling algorithm is verified through numerical simulation.
在微观尺度上,多孔介质的结构会受到其他结构的影响,因此其行为非常复杂,很难推导理论公式以及在宏观尺度上进行模拟。几十年来,为了避免这种复杂性,许多材料科学家和工程师对多尺度方法领域中的唯象方法进行了广泛研究。他们的理论方法基于使用较小尺度上先验知识的分层多尺度方法;然而,它有一个缺点,即可能会出现信息损失。最近,随着计算机核心技术的发展,多尺度方法扩展到了一种称为并发多尺度方法的直接多尺度方法。这种方法无需处理复杂的数学公式,但需要注意一个重要因素:多孔介质中各结构之间计算耦合算法的发展。在这项工作中,我们试图在不同的数值方法——有限元法(FEM)、光滑粒子流体动力学(SPH)和离散元法(DEM)中开发耦合算法。使用这种耦合算法,通过基于SPH、FEM和DEM提出的离散元来计算流体流动、固体颗粒的运动以及固体域之间的接触力。此外,引入了连续体水平上的混合有限元法和不连续体水平上带有SPH粒子的离散元模型,并通过数值模拟验证了所提出的耦合算法。