Laminar, transient forced convection problem over a 2D backward facing step (BFS) at an inlet Reynolds number (Re) of 400 is investigated numerically using OpenFOAM. To increase the Nusselt number (Nu) along the bottom wall, active flow control is applied by zero-net-mass-flux (ZNMF) combinations of suction and injection through three thin slits which are placed on the top, step and the bottom walls in the vicinity of the BFS. The combinations of each jet velocity is determined by jet to inlet mean velocity ratios which are limited to integer numbers between-2 and 2 and satisfying ZNMF condition where negative and positive values indicate suction and injection, respectively. All 19 cases which satisfy these rules are investigated. Average Nusselt number, friction coefficient and recirculation zone lengths are calculated along the bottom wall from time averaged flow fields. Among 19 cases with each having different jet configuration, some cases converged to steady state solution while others indicated temporal effects and converged to periodic solutions. To understand these transient effects, velocity oscillation magnitude and Strouhal number which are monitored at a selected critical point are evaluated. It is shown that temporal interaction of chosen active flow control methodology has significant effect on enhancing mixing which results in an increase of Nusselt number. Among all cases, the best case concerning thermal improvement has an increase of 78.5% in Nu number while the best aerodynamic improvement is achieved for another case with a decrease of 81% in total recirculation zone length compared to the reference case where no control is applied.