This paper seeks to implement fluid simulation in tetrahedral mesh using hybrid method. Since Jos Stam introduced Stable Fluid in computer graphics community, FDM (finite difference method) in regular grid (hexahedral mesh) became widely used. The drawback of FDM in regular mesh has to do with the fact that treating curved boundary presents a challenge. To handle irregular boundary, tetrahedral mesh was used as a simulation domain. For advection, modified PIC/FLIP method was used to take advantage of low numerical dissipation. Unlike the traditional linear interpolation in PIC/FLIP, SPH kernel function based velocity interpolation method was experimented.