In the IBM, an Eulerian description of the equations governing the fluid flow is used for modeling the fluid dynamics, and a Lagrangian description of the curvilinear boundary structural mechanics is used for objects immersed in the fluid. In this method,the boundary effect on the fluid is taken into account by spreading the surface force into the bulk fluid and treating it as a body force.In the present work, the immersed body is assumed to consist of mass less fibers, so that all of the force generated by distortions of the boundary can be calculated easily and transmitted directly to the fluid. Fig. 1 depicts a 2D fluid domain including a single closed immersed boundary fiber. In order to calculate the Lagrangian forces arising from the elastic energy of the immersed body, we need to specify the material points of the boundary. Thus, we use a Lagrangian description of the boundary. Suppose that X(s, t) is an arbitrary point on the immersed boundary where the parameter s (0 ≤ s ≤ l0) indicates the Lagrangian coordinates along the boundary and t denotes time. Here, l0is the length of the immersed boundary at equilibrium. Based on the principle of virtual work, the Lagrangian force F is derived using the variational derivative of the elastic energy functional E[X(s, t)], i.e.