Fluid-Elastic Instability is the most important mechanism among the vibration excitations in heat exchanger tube bundles subjected to cross flow. Flows through the heat exchanger are mostly two phase flow like air-water, vapor-water or Freon. Tow phase numerical methods are so complicated because of some parameters like VOF and interaction between two phases. Experimental studies are so problem and costly. Therefore, numerical methods are so important for studying two phase flow. In this Article, threshold of vibration has been numerically predicted by simulation of incompressible, viscose, and unsteady cross flow through a tube bundle in normal triangular arrangement. Interactions between the fluid and the structure has been counted in a fully coupled manner. HEM was used for analyzing two phase flow, In HEM method treats no difference in velocity between gas and fluid. In this study, two phase flow with HEM was solved around a single flexible cylinder surrounded by rigid and flexible tubes of bundle. Eventually, the flow through tube bundle has been simulated and analyzed by monitoring critical reduced velocity. Result shows that with increasing VoF, amplitude decreased and the critical velocity increased.