Time-domain simulation for fractional relaxation of Havriliak-Negami type

The time-simulation of models described by the Havriliak-Negami response function is a challenging problem due to the absence of an explicit formulation for the corresponding differential operator in the time-domain. In this work, we discuss a convolution quadrature rule with convolution weights evaluated on the basis of the representation of the response function in the Laplace transform domain. We describe a general and straightforward technique for the computation of the weights and we present some numerical experiments to illustrate the effectiveness of the proposed approach.