FEA-based Gaussian density heat source models were developed to study the effect of convective and radiative heat sinks on the transient temperature field predicted by the available approximate analytical solution of the purely conduction-based Goldak’s heat source. A new complex 3D Gaussian heat source model, incorporating all three modes of heat transfer, i.e., conduction, convection and radiation, has been developed as an extension of the Goldak heat source. Its approximate transient analytical solutions for this 3-D moving heat source were derived and numerically benchmarked with the available measured temperature & weld pool geometry data by Matlab programming with ~5 to 6 times faster than FEA-based simulation. The new complex 3D Gaussian heat source model and its approximate solution could significantly reduce the computing time in generating the transient temperature field and become an efficient alternative to extensive FEA-based simulations of heating sequences, where virtual optimisation of a melting heat source (i.e. used in welding, heating, cutting or other advanced manufacturing processes) is desirable for characterisation of material behaviour in microstructure evolution, melted pool, microhardness, residual stress and distortions.