In a diesel engine, an inside of a cylinder is divided into intra-cavity and extra-cavity regions. Ideal heat release rate waveform models, each formed of an isosceles triangle in which each oblique line gradient is a reaction rate, an area is a reaction amount and a base length is a reaction period with a reaction start temperature as a base point, are generated respectively for a vaporization reaction, low-temperature oxidation reaction, thermal decomposition reaction and high-temperature oxidation reaction of injected fuel for each region. An ideal heat release rate waveform of the reaction modes is generated by smoothing the ideal heat release rate waveform models through filtering and combining the ideal heat release rate waveforms, and is compared with an actual heat release rate waveform obtained from a detected in-cylinder pressure. A reaction mode having a deviation larger than or equal to a predetermined amount is diagnosed as being abnormal.