Blast furnace (BF) remains the dominant ironmaking process worldwide. Central coke charging (CCC) operation is a promising technology for stabilizing BF operations, but it needs reliable and quantified process design and control. In this work, a multi-fluid BF model is further developed for quantitatively investigating flow-thermal-chemical phenomena of a BF under CCC operation. This model features the respective chemical reactions in the respective coke and ore layers, and a specific sub-model of layer profile for the burden structure for the CCC operation. The simulation results confirm that the gas flow patterns and cohesive zone's shape and location under the CCC operation are quite different from the non-CCC operation. Under the CCC operation, the heat is overloaded at the furnace center while the reduction load is much heavier at the periphery regions; the profiles of top gas temperature and gas utilization show bell-shape and inverse-bell-shape patterns, respectively. More importantly, these differences are characterized quantitatively. In this given case, when the CCC opening radius at the throat is 0.35 m, the cohesive zone top opening radius is around 0.50 m, and the isotherms of CCC operation become much steeper (~80 deg) than those of non-CCC operation (~60 deg) near BF central regions. In addition, it is confirmed that carbon solution-loss reaction rate can be decreased significantly at BF central regions under CCC operation. The model helps to understand CCC operation and provides a cost-effective method for optimizing BF practice.