An investigation has been done to identify the heat source model suitable for numerical modeling of the Gas Metal arc Welding (GMAW) welding process. In the present investigation, a numerical model was constructed to simulate the GMAW welding process through transient thermal analysis for different thermal heat source models for weld deposition of Mild steel considering semi-cylindrical weld bead geometry. The Gaussian heat source with radial coordinate system and double ellipsoidal thermal heat source have been considered to develop a Finite Element Analysis (FEA) model. Validation of numerical analysis results has been confirmed through experimental investigation. Numerical modeling confirms the utilization of both heat sources, i.e. Gaussian and Double ellipsoidal for the prediction of the weld bead deposition profile. However, the numerical modeling developed using a double ellipsoidal heat source has more resembles experimental investigation due to the utilization of bead profile data. However, the Gaussian heat source is found to be suitable in the case of the unavailability of experimental data. The experimental studies show a variation in the weld bead geometry, along its welding direction due to multiple environmental factors. The microstructural investigation also shows the development of a fine columnar grain structure at the welding bead & weld penetration zone, and the grain structure became coarser at the Heat Affected Zone (HAZ).