The dominant role of electricity generation and environment consideration have placed strong requirements on coal fired power plants, requiring them to improve boiler combustion efficiency and decrease carbon emission. Although neural network based optimization strategies are often applied to improve the coal fired power plant boiler efficiency, they are limited by some combustion related problems such as slagging. Slagging can seriously influence heat transfer rate and decrease the boiler efficiency. In addition, it is difficult to measure slag build-up. The lack of measurement for slagging can restrict conventional neural network based coal fired boiler optimization, because no data can be used to train the neural network. This paper proposes a novel method of integrating non-dominated sorting genetic algorithm (NSGA II) based multi-objective optimization with computational fluid dynamics (CFD) to decrease or even avoid slagging inside a coal fired boiler furnace and improve boiler combustion efficiency. Compared with conventional neural network based boiler optimization methods, the method developed in the work can control and optimize the fields of flue gas properties such as temperature field inside a boiler by adjusting the temperature and velocity of primary and secondary air in coal fired power plant boiler control systems. The temperature in the vicinity of water wall tubes of a boiler can be maintained within the ash melting temperature limit. The incoming ash particles cannot melt and bond to surface of heat transfer equipment of a boiler. So the trend of slagging inside furnace is controlled. Furthermore, the optimized boiler combustion can keep higher heat transfer efficiency than that of the non-optimized boiler combustion. The software is developed to realize the proposed method and obtain the encouraging results through combining ANSYS 14.5, ANSYS Fluent 14.5 and CORBA C++. (C) 2014 Elsevier Ltd. All rights reserved.