Numerical computation to predict the pulverized-coal diffusion flame using Lagrangian treatment of the gas-solid flow and combustion is implemented. The authors analyse and discuss the gas-particle momemtum, mass and energy transfer processes. A modified scheme for existing computation model is proposed by the authors, which is effective in promoting two phase coupling process and can overcome divergency in computation for high pulverized-coal concentration combustion.