Consider unilateral series in a single variable where the exponent is a polynomial, and the coefficients are periodic. Quadratic polynomials correspond to partial theta series. Such series converge inside the unit disk. We compute limits and asymptotic expansions of those series as the variable tends radially to a root of unity. The method is suitable for automation and part of the computations is based on the idea of the $q$-integral.