Abstract A mathematical model system was derived to describe the simultaneous removal of phenol biodegradation with chromium(VI) reduction in an anaerobic fixed-biofilm reactor. The model system incorporates diffusive mass transport and double Monod kinetics. The model was solved using a combination of the orthogonal collocation method and Gear's method. A laboratory-scale column reactor was employed to validate the kinetic model system. Batch kinetic tests were conducted independently to evaluate the biokinetic parameters used in the model simulation. The removal efficiencies of phenol and chromium(VI) in an anaerobic fixed-biofilm process were approximately 980 mg/g and 910 mg/g, respectively, under a steady-state condition. In the steady state, model-predicted biofilm thickness reached up to 350 μm and suspended cells in the effluent were 85 mg cell/l. The experimental results agree closely with the results of the model simulations.