Thermal debinding is one of the most important steps in powder injection moulding process. Thermogravimetric analysis (TGA) was employed to analyze thermal debinding behaviour under argon atmosphere. Thermal debinding kinetics with different heating rates have been compared using 316L stainless steel feedstocks loaded typically at 60, 62, 64 and 66 vol.% for fine metallic powders D-50 = 3.4 mu m. The Kissinger and Ozawa methods have been used to estimate the kinetic parameters from thermogravimetric experiments. To set up the numerical simulations of thermal debinding stage using finite element methods, a coupled mathematical model for mass diffusion and heat transfer in deformable porous media have been developed. The basic steps of the proposed model consist of solving the following coupled problems: thermal degradation of binder coupled with heat transfer and deformation phenomena by finite element methods using Comsol Multiphysics (R) software. The obtained numerical simulation results are in proper agreement with experimental data. The proposed numerical simulations allow the determination of remaining binder distribution, temperature distribution and deformation fields in the component during the whole thermal debinding process at any time. (C) 2012 Elsevier B.V. All rights reserved.