Abstract: In order to study the filling process of the plastic encapsulation, a stabilized three dimensional simulation method was presented. In this method, NS equations combined with Kamal conversion equation were introduced as governing equations, and FAN(Flow Analysis Network) technique was used to track melt fronts for each time step. To avoid numerical surge during calculation and achieve stabilized scheme, GLS(Galerkin Least Square) scheme was used in the momentum equation, GLS\\GGLS(Gradient Galerkin Least Square) scheme was used in the energy equation, and SUPG(Streamline Upwind Petrov Galerkin) scheme was used in the conversion equation. The experiment proves this method has a good stabilization and accuracy.