Two Green's function-based numerical formulations are used to solve the time-dependent nonlinear heat conduction (diffusion) equation. These formulations, which are an extension of the first paper, utilize two fundamental solutions and the Green's second identity to achieve integral replications of the governing partial differential equation. The integral equations thus derived are discretized in space and time and aggregated in a finite element sense to give a system of nonlinear discrete equations that are solved by the Newton–Raphson algorithm. The mathematical simplicity of the Green's function of the first formulation facilitates its numerical implementation. The performance of the formulations is assessed by comparing their results with available numerical and analytical solutions. In all cases satisfactory and physically realistic results are obtained.