Wave-body interaction is simulated using a developed code based on the flux-difference splitting scheme for immiscible and incompressible fluids and the hybrid Cartesian/immersed boundary method. A free surface is captured as a moving contact discontinuity within a fluid domain and an approximated Riemann solver is used to estimate the inviscid flux across the discontinuity. Immersed boundary nodes are identified inside an instantaneous fluid domain near a moving body, then dependent variables are reconstructed at those immersed boundary nodes based on interpolation along local normal lines to the boundary. Free surface flows around an oscillating cylinder are simulated and the computed wave elevations are compared with other reported results. The generation of a solitary wave by a moving wave-maker is simulated and the time histories of wave elevations at two different points are compared with other results. The developed code is applied to simulate body motion of an elastically mounted circular cylinder as a solitary wave passes the body. The force acting on an elastically mounted cylinder is compared with the force acting on a fixed cylinder. Grid independency of the computed body motion is established based on a comparison of results using three different-size grids.