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.