A model that predicts the physical changes that pulverized coal char particles undergo during combustion has been developed. In the model, a burning particle is divided into a number of concentric annular volume elements. The mass loss rate, specific surface area, and apparent density in each volume element depend upon the local particle conditions, which vary as a consequence of the adsorbed oxygen and gas-phase oxygen concentration gradients inside the particle. The model predicts the particle's burning rate, temperature, diameter, apparent density, and specific surface area as combustion proceeds, given ambient conditions and initial char properties. A six-step heterogeneous reaction mechanism is used to describe carbon reactivity to oxygen. A distributed activation energy approach is used to account for the variation in desorption energies of adsorbed O-atoms on the carbonaceous surface. Model calculations support the three burning zones established for the oxidation of pulverized coal chars. The model indicates two types of zone II behavior, however. Under weak zone II burning conditions, constant-diameter burning occurs up to 30% to 50% conversion before burning commences with reductions in both size and apparent density. Under strong zone 11 conditions, particles burn with reductions in both size and apparent density after an initial short period (< 2% conversion) of constant-diameter burning. Model predictions reveal that early in the oxidation process, there is mass loss at constant diameter under all zone II burning conditions. Such weak and strong burning behavior cannot be predicted with the commonly used power-law model for the mode of burning employing a single value for the burning mode parameter. Model calculations also reveal how specific surface area evolves when oxidation occurs in the zone II burning regime. Based on the calculated results, a surface area submodel that accounts for the effects of pore growth and coalescence during combustion under zone I conditions was modified to permit the characterization of the variations in specific surface area that occur during char conversion under zones II conditions. The modified surface area model is applicable to all burning regimes. Calculations also indicate that the particle's effectiveness factor varies during conversion under zone II burning conditions. With the adsorption/desorption mechanism employed, a near first-order Thiele modulus-effectiveness factor relationship is obeyed over the particle's lifetime. (C) 2007 The Combustion Institute. Published by Elsevier Inc. All rights reserved.