A two-dimensional axisymmetric numerical model, including the influence of the cathode and the free surface of the weld pool, is developed to describe the heat transfer and fluid flow in gas tungsten arc (GTA) welding. In the model, a boundary-fitted coordinate system is adopted to precisely describe the cathode shape and deformed weld-pool surface. The current continuity equation has been solved with the combined arc plasma-cathode system, independent of the assumption of current density distribution on the cathode surface, which was essential in the previous studies of arc plasma. It has been shown that the temperature profile, the current, and the heat flux to the anode show good agreement with the experimental data. Moreover, the current and the heat-flux distributions may be affected by the shape of the cathode and the free surface of the weld pool.