It is one of challenging topics to carry out large-scale simulations for multiphase flows such as gas-liquid or and solid-gas two-phase flows. We have developed a weakly compressible flow solver and applied AMR (Adaptive Mesh Refinement) fine meshes to the gas-liquid interfaces. Dynamic domain partitioning based on space-filling curves has been employed and a few simulations on liquid film dynamics and tsunami flows with a lot of debris are carried out on 10-100 GPUs of the GPU supercomputer TSUBAME 3.0.