A numerical prediction method has been proposed to predict non-linear free surface oscillation in a three-dimensional container. The fluid motions are numerically predicted with Navier-Stokes equations discretized in a Lagrangian scheme with sufficient numerical accuracy. The profile of a free surface is precisely represented with three-dimensional body-fitted coordinates (BFC), which are regenerated in each computational step on the basis of the arbitrary Lagrangian-Eulerian (ALE) formulation. The computational method was applied to non-linear sloshings and transitions from sloshing to swirling motions. The predicted free surface motions were visualized as sequential image files and animations to understand their dynamic futures