Abstract: To solve the problem of unreasonable boundary conditions of the existing SPH non-Newtonian fluid simulation,a boundary handling method is proposed.First,a unified modeling of Newtonian and non-Newtonian fluids was implemented by Cross model.Then,by calculating the weighted sum of solid boundary particles,a normal vector of the boundary force was set for the fluid particles at the boundary.Finally,the velocity of fluid particles near the boundary was iteratively updated using Coulomb friction boundary conditions,which can be integrated into the predictive-corrective algorithm framework.Experiments show that the friction coefficient of the fluid-solid boundary can be adjusted according to the user’s experimental needs.Compared to simulation under free-slip boundary conditions,the non-Newtonian fluid simulated by our method is able to increase in viscosity with the decrease of velocity and exhibit physical properties of solids.