The equations governing heat and fluid flow in weld pools for the TIG fusion welding process are presented and this coupled system is solved numerically using finite differences. Electromagnetic forcing terms, buoyancy forces, shear forces on the pool surface due to the variation in surface tension with temperature and an additional uniform magnetic field applied normal to the workpiece are all included in our model and results are displayed indicating the relative importance of these four mechanisms.