Analysis of fluid-structure interaction for two nearby underwater vehicles immersed in the sea is quite challenging because simulation of flow around them is very difficult due to the complexity of underwater vehicle shapes. The conventional approach using body-fitted or unstructured grids demands much time in dynamic grid generation, and yields slow convergence of solution. Since an analysis of fluid-structure interaction must be based on accurate simulation results, a more efficient way of simulating flow around underwater vehicles, without sacrificing accuracy, is desirable. An immersed boundary method facilitates implementation of complicated underwater-vehicle shapes on a Cartesian grid system. An LES modeling is also incorporated to resolve turbulent eddies. In this paper, we will demonstrate the effectiveness of the immersed boundary method we adopted, by presenting the simulation results on the flow around a modeled high-speed underwater vehicle interacting with a modeled low-speed one.