A meshless local Petrov-Galerkin (MLPG) method is proposed to obtain the numerical solution of nonlinear heat transfer problems. The moving least squares scheme is generalized, to construct the field variable and its derivative continuously over the entire domain. The essential boundary conditions are enforced by the direct scheme. The radiation heat transfer coefficient is defined, and the nonlinear boundary value problem is solved as a sequence of linear problems each time updating the radiation heat transfer coefficient. The matrix formulation is used to drive the equations for a 3 dimensional nonlinear coupled radiation heat transfer problem. By using the MPLG method, along with the linearization of the nonlinear radiation problem, a new numerical approach is proposed to find the solution of the coupled heat transfer problem. A numerical study of the dimensionless size parameters for the quadrature and support domains is conducted to find the most appropriate values to ensure convergence of the nodal temperatures to the correct values quickly. Numerical examples are presented to illustrate the applicability and effectiveness of the proposed methodology for the solution of heat transfer problems involving radiation with different types of boundary conditions. In each case, the results obtained using the MLPG method are compared with those given by the FEM method for validation of the results.