Assuming $X$ and $Y$ follow a bivariate normal distribution with correlation $\rho\in[0,1)$, zero-means, and variances $\sigma_x^2$ and $\sigma_y^2$, respectively. Then, the pdf of the random variable $Z = \sqrt{X^2+Y^2}$ follows the Hoyt distribution, that is
\begin{align*}
f_Z(z)
&=\frac{z}{\sqrt{1-\rho ^2} \sigma_x \sigma_y}
\exp\!\left(-\frac{z^2 \left(\sigma_x^2+\sigma_y^2\right)}{4 \left(1-\rho ^2\right) \sigma_x^2 \sigma_y^2}\right)
I_0\!\left(\frac{z^2 \sqrt{4 \rho ^2 \sigma_x^2 \sigma_y^2+\left(\sigma_x^2-\sigma_y^2\right)^2}}{4 \left(1-\rho ^2\right) \sigma_x^2 \sigma_y^2}\right),~ z>0,
\end{align*}
where $I_0$ is the modified Bessel function of the first kind.
Indeed, if $\sigma_x = \sigma_y = \sigma$ and $\rho = 0$, we have the Rayleigh distribution: $f_Z(z) = \frac{z}{\sigma^2} e^{-z^2/(2\sigma^2)}$.
Additional comment: Note that we can always find $X'$ and $Y'$ with variances $\sigma_x'^2$ and $\sigma_y'^2$, respectively, that are uncorrelated and satisfies $\sqrt{X'^2+Y'^2} = \sqrt{X^2+Y^2} = Z$. Then, the pdf of $Z$ can be written in the simpler form
\begin{align*}
f_Z(z)
&=\frac{z}{\sigma_x' \sigma_y'} \exp\!\left( -\frac{z^2}{4}\left(\frac{1}{\sigma_x'^2} + \frac{1}{\sigma_y'^2}\right) \right) I_0\!\left( \frac{z^2}{4} \left(\frac{1}{\sigma_x'^2} - \frac{1}{\sigma_y'^2}\right)\right),~ z>0.
\end{align*}