It is known that the reduction process influences the initial performances and durability of nickel-yttria-stabilized zirconia composite anode of the solid oxide fuel cell. In the present study, the reduction process of nickel-yttria stabilized zirconia composite anode is simulated based on the phase field method. An three-dimensional reconstructed microstructure of nickel oxide-yttria stabilized zirconia composite obtained by focused ion beam-scanning electron microscopy is used as the initial microstructure for the simulation. Both reduction of nickel oxide and nickel sintering mechanisms are considered in the model. The reduction rates of nickel oxide at different interfaces are defined based on the literature data. Simulation results are qualitatively compared to the experimental anode microstructures with different reduction temperatures. (C) 2015 Elsevier B.V. All rights reserved.