The Boundary Element Method (BEM) allows efficient solution of partial differential equations whose kernel functions are known. The heat equation is one of these candidates when the thermal parameters are assumed constant (linear model). When the model involves large physical domains and time simulation intervals the amount of information that must be stored increases significantly. This drawback can be circumvented by using advanced strategies, as for example the multi-poles technique. We propose radically different approach that leads to a separated solution of the space and time problems within a non-incremental integration strategy. The technique is based on the use of a space-time separated representation of the unknown field that, introduced in the residual weighting formulation, allows to define a separated solution of the resulting weak form. The spatial step can be then treated by invoking the standard BEM for solving the resulting steady state problem defined in the physical space. Then, the time problem that results in an ordinary first order differential equation is solved using any standard appropriate integration technique (e.g. backward finite differences). When considering the nonlinear heat equation, the BEM cannot be easily applied because its Green’s kernel is generally not known but the use of the PGD presents the advantage of rewriting the problem in such a way that the kernel is now clearly known. Indeed, the system obtained by the PGD is composed of a Poisson equation in space coupled with an ODE in time so that the use of the BEM for solving the spatial part of the problem is efficient. During the solving, we must however separate the nonlinear term into a space-time representation that can limit the method in terms of CPU time and storage, that is why we introduce in the second part of the paper a new approach combining the PGD and the Asymptotic Numerical Method (ANM) in order to efficiently treat the nonlinearity.< Réduire