It is well known experimentally that the circuit matching network plays an important role in the mode transition behavior of inductively coupled plasmas. To date, however, there have been no reports of numerical models being used to study the role of the matching circuit in the transition process. In this paper, a new two-dimensional self-consistent fluid model that couples the components of an equivalent circuit module is developed to investigate the effects of the equivalent circuit on the mode transition characteristics of an inductively coupled, hydrogen plasma. The equivalent circuit consists of a current source, impedance matching network, reactor impedance, and plasma transferred impedance. The nonlinear coupling of the external circuit with the internal plasma is investigated by adjusting the matching capacitance at a fixed input current. The electron density and temperature as well as the electromagnetic fields all change suddenly, and the E to H mode transition occurs abruptly at a certain matching capacitance as the impedance matching of the external circuit is varied. We also analyze the fields and the plasma characteristics during the transition process, especially for the case of the capacitive E mode. V C 2015 AIP Publishing LLC.