Abstract We develop a self-consistent theoretical formalism to model the dynamics of heat transfer in dissipative, dispersive, anisotropic nanoscale media, such as metamaterials. We employ our envelope dyadic Green’s function method to solve Maxwell’s macroscopic equations for the propagation of fluctuating electromagnetic fields in these media. We assume that the photonic radiative heat transfer mechanism in these media is complemented by dynamic phononic mechanisms of heat storage and conduction, accounting for effects of local heat generation. By employing the Poynting theorem and the fluctuation-dissipation theorem, we derive novel closed-form expressions for the radiative heat flux and the coupling term of photonic and phononic subsystems, which contains the heating rate and the radiative heat power contributions. We apply our formalism to the paraxial heat transfer in uniaxial media and present relevant closed-form expressions. By considering a Gaussian transverse temperature profile, we also obtain and solve a system of integro-differential heat diffusion equations to model the paraxial heat transfer in uniaxial reciprocal media. By applying the developed analytical model to radiative-conductive heat tranfer in nanolayered media constructed by layers of silica and germanium, we compute the temperature profiles for the three first orders of expansion and the total temperature profile as well. The results of this research can be of interest in areas of science and technology related to thermophotovoltaics, energy harvesting, radiative cooling, and thermal management at micro- and nanoscale.