A penalty immersed boundary method for viscoelastic particulate flows

•A numerical method is introduced to investigate particulate viscoelastic fluid.•The method is tested by a convergence study.•The effect of the elasticity number and channel wall is presented.•A chain is formed with multiple particles falling in an Oldroyd-B fluid.