A new numerical approach for solving two-phase flows is presented in the framework of the recently developed Consistent Particle Method (CPM). The spatial derivatives along fluid interfaces where abrupt discontinuity of fluid density exists are computed by including neighbor particles of another fluid and applying Taylor series expansion on specific pressure energy. In addition, a new h-adaptive particle selection scheme is proposed to overcome the potential problem of ill-conditioned coefficient matrix of PPE when particles are sparse and non-uniformly spaced. Using this method, water-air sloshing in a closed rectangular tank is studied with experimental validation.
