Algebraic Reynolds stress model (ARSM) is often employed in practical turbulent flow simulations. Most of previous works on ARSM have been carried out for incompressible flows. In the present paper, a new ARSM model is suggested for compressible flows. The model adopts a compressibility factor function involving the turbulent Mach number and the gradient Mach number. Compared to incompressible flow, explicit solution for ARSM for compressible flow can hardly be obtained due to dilatation terms. We propose approximate representations for these dilatation-related terms to obtain an explicit procedure for compressible flow turbulence. The model is applied to compressible mixing layer, supersonic flat-plate boundary and planar supersonic wake flow. It is found that the model works very well yielding results that are in good agreement with the DNS and the experimental data.