A computational model is developed to represent the coupled behavior of elementary chemistry, electrochemistry, and transport in the vicinity of solid-oxide fuel cell three-phase boundaries. The model is applied to assist the development and evaluation of H 2 charge-transfer reaction mechanisms for Ni-yttria-stabilized zirconia ͑YSZ͒ anodes. Elementary chemistry and surface transport for the Ni and YSZ surfaces are derived from prior literature. Previously published patterned-anode experiments ͓J. Mizusaki et al., Solid State Ionics, 70/71, 52 ͑1994͔͒ are used to evaluate alternative electrochemical charge-transfer mechanisms. The results show that a hydrogen-spillover mechanism can explain the Mizusaki polarization measurements over wide ranges of gas-phase composition with both anodic and cathodic biases. Establishing elementary charge-transfer reaction pathways is one of the most difficult aspects of understanding the fundamental behavior of solid oxide fuel cells ͑SOFCs͒. Most modeling efforts represent charge transfer in terms of activation overpotentials and the Butler-Volmer equation. The objective of the present paper is to develop an elementary representation of the transport, chemistry, and electrochemistry in the vicinity of an anode three-phase boundary ͑TPB͒ between gas-phase H 2 , a Ni anode, and a yttria-stabilized zirconia ͑YSZ͒ electrolyte. The approach is based upon quantitative models and previously reported patterned-anode experiments by Mizusaki et al. Mizusaki et al.2 developed model anode structures using thin nickel films on single-crystal YSZ. The nickel films were patterned lithographically into an array of lines. By varying the linewidth and spacing, it was possible to control the TPB length per unit area. Others have followed Mizusaki's lead, and several experimental studies have been reported. [3][4][5][6][7] Nevertheless, Mizusaki's 1994 paper stands as the most complete data set available. The experiments were carefully designed to impose both cathodic and anodic potentials and to systematically vary the H 2 partial pressure while holding H 2 O partial pressure constant and vice versa.In addition to electrochemical charge transfer, the model simultaneously represents the competing processes of surface-diffusive transport and nonfaradaic chemistry on the Ni and YSZ surfaces. The reaction pathways and rates for thermal surface reactions are drawn from independent experimental and theoretical literature. The reaction mechanisms are constructed hierarchically, beginning with the nonfaradaic processes, and then adding charge-transfer reactions at the TPB, without modifying the rates or equilibrium constants for the nonfaradaic reactions. In this way, it is possible to develop robust reaction mechanisms that are applicable over wide ranges of operating conditions.The model is used to evaluate the validity of two alternative charge-transfer pathways. One is an oxygen-spillover mechanism in which the rate-limiting reaction involves atomic oxygen from the YSZ surface crossing over to the Ni...