In our study, the Ni/YSZ ReaxFF reactive force field was developed by combining the YSZ and Ni/C/H descriptions. ReaxFF reactive molecular dynamics (RMD) were applied to model chemical reactions, diffusion, and other physicochemical processes at the fuel/Ni/YSZ interface. The ReaxFF RMD simulations were performed on the H2/Ni/YSZ and C4H10/Ni/YSZ triple-phase boundary (TPB) systems at 1250 and 2000 K, respectively. The simulations indicate amorphization of the Ni surface, partial decohesion (delamination) at the interface, and coking, which have indeed all been observed experimentally. They also allowed us to derive the mechanism of the butane conversion at the Ni/YSZ interface. Many steps of this mechanism are similar to the pyrolysis of butane. The products obtained in our simulations are the same as those in experiment, which indicates that the developed ReaxFF potential properly describes complex physicochemical processes, such as the oxide-ion diffusion, fuel conversion, water formation reaction, coking, and delamination, occurring at the TPB and can be recommended for further computational studies of the fuel/electrode/electrolyte interfaces in a SOFC.
Keywords: ReaxFF; SOFC; reactive molecular dynamics; triple-phase boundary.