In this research, interaction between eel and vortex shedding behind a bluff body is studied numerically, aiming to optimize the oscillation of eels by coupling with the altering vertices. Effects of different factors such as length of eel, width of bluff body and flow speed on vibration of eels are also aimed to be investigated. A 2-dimensional Fluid-Structure Interaction Model that simulates the motion of eel under vortex shedding behind a bluff body is constructed by using the computational fluid dynamics (CFD) software – COMSOL Multiphysics. Simulations based on laminar flow regime are performed and the results show that undulating motion of the eel is successfully modeled with the Fluid-Structure Interaction Model along with Arbitrary Lagrangian-Eulerian (ALE) method.