This article provides a computational analysis of the reaction of fully-dense layered aluminum and copper oxide systems. After the detailed presentation of the 2D nonstationary model implementing both oxygen and aluminum diffusion, the propagation of the reaction front in an Al/CuO thin film was studied. The model qualitatively reproduces the dependency of the reaction front progression rate spatially as a function of the fuel concentration. Calculations also evidence the inverse evolution of flame front width with respect to the reaction front velocity. A procedure to estimate the heat loss generated by the fact that the reactants and products may vaporize prior to reaction completion was proposed by imposing a flame temperature limit close to Cu vaporization point. This work also shows that microscopic fluctuations in the instantaneous reaction front velocity can be observed for reactant diffusion activation energy (Ea) of 125 kJ/mol, before quenching for greater Ea. Finally, we demonstrate the potential of this new 2D nonstationary model to investigate the thermal effect of additives such as metallic impurities in the Al/CuO thin film that can lead to the flame front corrugation at the microscale. The simulations show that a metallic particle acts first to boost the reaction velocity as its high thermal conductivity helps the upfront heating. Then, the particle being also a heat sink, a local slowing down of the front velocity is observed. Nomenclature Symbol Meaning Unit C Total concentration of the material (equivalent to a density) kg.m-3 Ci Concentration of the species i kg.m-3 Cp Average specific heat capacity J.kg-1 .K-1 Cp,i Specific heat capacity of the species i J.kg-1 .K-1 Di Diffusion coefficient of the species i m².s-1 Activation energy of the Arrhenius law defining the reaction rate of the copper oxide decomposition J.mol-1 , ℎ Activation energy of the Arrhenius law defining the reaction rate of the copper oxide decomposition J.mol-1 Fs Enthalpy flux related to species flux kg.m-2 .s-1 ℎ , Enthalpy of formation of the species i J.kg-1 0, ℎ Pre-exponential factor of the Arrhenius law defining the reaction rate of the copper oxide decomposition s-1 0, Pre-exponential factor of the Arrhenius law defining the diffusion of the species i m².s-1 L Simulated length µm Lt Total thickness µm n Number of bilayer PHrx Power of reaction W.m-3 PPC Required power for phase change W.m-3 R Gas constant (R=8.314) J.K-1 .mol-1 ri Reaction rate of the species i kg.m-3. s-1 Heat conductivity of the species i W.m-2 .K-1 av Average heat flux going through the heat front W.m-2 av.loss Average heat flux lost by limiting the temperature with Tmax W.m-2