Simultaneous C-N-S removal from industrial wastewater was mathematically modeled.
Interactions between SRB, NRB, SOB, FB and MPA were depicted in the model.
Kinetic parameters for competitive coefficient between SRB and NRB were best estimated.
The developed mode was the first simulation tool for complicated C-N-S dynamics.
The model was versatile and feasible to predict relevant processes (e.g., SR-SO, DSR).