The mixed oil and gas including water and sand are extracted from well to offshore structure. This mixed fluid must be separated for subsequent processes by using wash tanks or separators. To design such a system, a proper numerical-prediction tool for multiphase fluids is required. In this regard, a new moving particle simulation (MPS) method is developed to simulate multiliquid-layer sloshing problems. The new MPS method for multifluid system includes extra search methods for interface particles, boundary conditions for interfaces, buoyancy-correction model, and surface-tension model for interface particles. The new particle interaction models are verified through comparisons with published numerical and experimental data. In particular, the multiliquid MPS method is verified against Molin et als (2012) experiment with three liquid layers. In case of excitation frequency close to one of the internal-layer resonances, the internal interface motions can be much greater than top free-surface motions. The verified multiliquid MPS program is subsequently used for more nonlinear cases including multichromatic multimodal motions with larger amplitudes, from which various nonlinear features, such as internal breaking and more particle detachment, can be observed. For the nonlinear case, the differences between with and without buoyancy-correction and surface-tension models are also demonstrated.