In this paper we present a novel mathematical model to describe the permeation of a fluid through a polymeric matrix, loaded with drug molecules, followed by its subsequent desorption. Both phenomena are enhanced by temperature. We deduce energy estimates and stability estimates for the weak solution of the model, showing that this solution of the problem is stable in bounded time intervals. Numerical simulations illustrate how the coupling effects, of viscoelastic properties and thermal external assistance, can have a central role in the design of drug delivery devices.