A random variable $X$ has a probability density function (PDF) given by $f(x) = 3x^2$ for $0 \le x \le 1$. If a new random variable $Y$ is defined as $Y = 2X + 1$, what is the variance of $Y$?