PurposeThe purpose of this paper is to recommend a validated numerical model for simulation the flue gases heat recovery recuperators. Due to fulfill of this demand, the influences of ash fouling characteristics during the transient/steady-state simulation and optimization of a 3D complex heat exchanger equipped with inner plain fins and side plate fins are studied.Design/methodology/approachFor the particle dispersion modeling, the discrete phase model is applied and the flow field has been solved using SIMPLE algorithm.FindingsAccording to obtained results, for the recuperator equipped with combine inner plain and side plate fins, determination of ash fouling characteristics is really important, effective and determinative. It is clear that by underestimating the ash fouling characteristics, the achieved results are wrong and different with reality.Originality/valueFinally, the configuration with inner plain fins with characteristics of: di =5 mm, do = 6 mm, dg = 2 mm, dk = 3 mm and NIPFT = 9 and side plate fins with characteristics of: TF = 3 mm, PF = 19 mm, NSPF = 17·2 = 34, WF = 10 mm, HF = 25 mm, LF = 24 mm and ß = 0° is introduced as the optimum model with the best performance among all studied configurations.