Large amplitude sloshing in tanks is simulated by the least square particle finite element method （LSPFEM） in this paper. The least square finite element method （LSFEM） is employed to spatially discrete the Navier-Stokes equations, and to avoid the stabilization issues due to the incompressibility condition for equal-order interpolation of the velocity and the pressure, as usually used in Galerkin method to satisfy the well-known LBB condition. The LSPFEM also uses the Lagrangian description to model the motion of nodes （particles）. A mesh which connects these nodes is constructed by a triangulation algorithm to avoid the mesh distortion. A quasi a-shapes algorithm is used to identify the free surface boundary. The nodes are viewed as particles which can freely move and even separate from the main fluid domain. Finally this method is used to study the large amplitude sloshing evolution in two dimensional tanks. The results are compared with those obtained by Flow-3d with good agreement.
