A methodology has been developed to achieve rapid two‐dimensional self‐consistent simulation of plasma transport and reaction in an inductively coupled source of arbitrary geometry and with arbitrary plasma and surface chemistries. In this modular finite element fluid simulation the reactor was divided into bulk plasma and sheath. The bulk plasma was assumed quasineutral and the electrons were assumed to be in Boltzmann equilibrium. Separate modules computed the power deposition into the plasma, electron temperature, charged species densities, and neutral species densities. Simulation results agreed favorably with available experimental data, taken in a chlorine plasma in a Gaseous Electronics Conference reference cell, without using any adjustable parameters. Rapid convergence makes the simulation tool especially attractive for technology computer‐aided design (TCAD) applications. ©1996 American Institute of Physics.