Proporcionaré una demostración directa no totalmente rigurosa. Me encantaría ver los detalles completamente rigurosos.
Recuerda la función hipergeométrica de Gauss $${}_2F_1 (a,b;c;z) = \displaystyle\sum_{k=0}^{\infty} \dfrac{a^{\overline{k}}b^{\overline{k}}}{c^{\overline{k}}} \dfrac{z^k}{k!},$$ donde $\xi^{\overline{k}}$ denota el factorial creciente $$\xi^{\overline{k}}=\dfrac{\Gamma(\xi+k)}{\Gamma(\xi)}=\xi (\xi+1) \ldots (\xi+k-1).$$ En particular, nota que tomamos $\xi^{\overline{0}}$ como $1$. Esto muestra que $${}_2F_1(a,0;c;z) = 1,$$ una función constante. Para nuestros propósitos, $$\displaystyle\lim_{|z| \rightarrow \infty} {}_2F_1(a,0;c;z)=1.$$
Wolfram alpha informa que la antiderivada necesaria es%20dx) $$\displaystyle\int \sqrt{1 + ax^b} \mathrm{d}x = \dfrac{xb {}_2F_1\left( \frac{1}{2}, \frac{1}{b};1+\frac{1}{b};-ax^b \right)}{b+2}+\dfrac{2x\sqrt{ax^b+1}}{b+2}+C.$$ Entonces, con $a=t^2$ y $b=2t-2$, el teorema fundamental del cálculo dice $$\begin{array}{ll} \displaystyle\lim_{t\to \infty}\int_0^1 \sqrt{1+t^2x^{2t-2}} \mathrm{d}x &= \displaystyle\lim_{t \to \infty} \left[\dfrac{(2t-2){}_2F_1\left(\frac{1}{2},\frac{1}{2t-2};1+\frac{1}{2t-2};-t^2\right)}{2t} + \dfrac{2\sqrt{t^2+1}}{2t}{}\right] \\ &= \displaystyle\lim_{t \to \infty} \left[ \left(\dfrac{t-1}{t} \right) {}_2F_1\left( \dfrac{1}{2}, \dfrac{1}{2t-2};1+\dfrac{1}{2t-2};-t^2 \right) + \dfrac{\sqrt{t^2+1}}{t}\right]. \end{array}$$ Ya que $$\displaystyle\lim_{t \to \infty} \dfrac{t-1}{t} = \displaystyle\lim_{t \to \infty}\dfrac{\sqrt{t^2+1}}{t}=1,$$ para finalmente responder la pregunta, debemos discutir el comportamiento límite de la función ${}_2F_1$.
Un enfoque riguroso consideraría fórmulas asintóticas de la serie hipergeométrica y deduciría el límite a partir de ellas. Nos relajaremos en rigor y simplemente tomaremos $t \rightarrow \infty$ en todos los parámetros a ${}_2F_1$ sugiriendo que el valor límite debería ser $${}_2F_1 \left( \dfrac{1}{2}, 0; 1; -\infty \right)=1.$$ Esto significa que tendríamos $$\displaystyle\lim_{t \to \infty} \displaystyle\int_0^1 \sqrt{1+t^2 x^2} \mathrm{d}x = 1+1=2,$$ como se deseaba.