Abstract: An implicit algorithm of Bi-penalty approximation with orthogonality projection for the numerical simulation of Bingham fluid flow problems is proposed in this paper.A Newton fluid flow with two kinds of artificial viscosity subjected to the inequality constraint is introduced to approximate the Bingham fluid flow This approach can effectively simulate the Bingham fluid flow with floating rigid cores or fixing rigid cores.