AbstractParticle methods for high-viscous free-surface flows are of great use to capture flow behaviors which are intermediate between solid and liquid. In general, it is important for numerical methods to satisfy the fundamental laws of physics such as the conservation laws of mass and momentum and the thermodynamic laws. Especially, the angular momentum conservation is necessary to calculate rotational motion of high-viscous objects. However, most of the particle methods do not satisfy the physical laws in their spatially discretized system. The angular momentum conservation law is broken mostly because of the viscosity models, which may result in physically strange behavior when high-viscous free-surface flow is calculated. In this study, a physically consistent particle method for high-viscous free-surface flows is developed. The present method was verified, and its performance was shown with calculating flow in a rotating circular pipe, high-viscous Taylor–Couette flow, and offset collision of a high-viscous object.