I imagine this is already found but I cannot find the proof.
The formula also works for $f_1=\sqrt{2}^{1/x}$ to get many multiples of $\pi$
Can you also prove that when $f_1=\sqrt{\frac{2}{\pi}}$ then $f_n$ tends to $\sqrt{n}-.5$ ?
Example for $f_1=2$
Example for $f_1=\sqrt{\frac{2}{\pi}}$


Leave $f_1$ arbitrary; $f_{n}f_{n-1}=f_{n-1}f_{n-2}+1$ implies $f_{n}f_{n+1}=n$ (by induction). Let $$ f_n=\frac{\Gamma\left(\frac{n+1}2\right)}{\Gamma\left(\frac n2\right)}a_n; $$ then $a_{n}a_{n+1}=2$, thus $a_{2n+1}=a_1$ and $a_{2n+2}=2/a_1$ for all $n$. Now the known limit $$ \lim_{x\to+\infty}\frac{\Gamma(x+a)}{x^a\ \Gamma(x)}=1\qquad(a\in\mathbb{R}) $$ easily gives $\displaystyle\color{blue}{\lim_{n\to\infty}\frac{f_{2n+1}}{\sqrt n}}=a_1\color{blue}{=f_1\sqrt\pi}$ and $\displaystyle\color{blue}{\lim_{n\to\infty}\frac{f_{2n}}{\sqrt n}}=\frac2{a_1}\color{blue}{=\frac2{f_1\sqrt\pi}}$.
(An elementary alternative here is to use the Wallis product.)
For $f_1=\sqrt{2/\pi}$ we get $a_n=\sqrt2$ identically, and the asymptotics of $$ f_n=\frac{n}{\sqrt{2\pi}}\mathrm{B}\left(\frac{n+1}2,\frac12\right)=n\sqrt\frac2\pi\int_0^1\frac{x^n\,dx}{\sqrt{1-x^2}}=n\sqrt\frac2\pi\int_0^\infty\frac{e^{-nt}\,dt}{\sqrt{e^{2t}-1}} $$ can be obtained using Watson's lemma. The (correct) result is $$ f_n\asymp n^{1/2}\color{red}{-\frac14 n^{-1/2}}\color{gray}{+\frac1{32}n^{-3/2}+\frac5{128}n^{-5/2}+\dots} $$