The fact that $D$ is a disc tells you that you should use polar coordinates. Therefore, let $x = r \cos t$ and $y = r \sin t$ with $r \in [0,1]$ and $t \in [0, \frac \pi 2]$, in order to transform your integral into
$$\iint \limits _{[0,1] \times [0, \frac \pi 2]} \frac {r \cos t r \sin t } {\sqrt {1 - r^2 \sin^2 t}} r \Bbb d r \Bbb d t .$$
Notice that the derivative with respect to $t$ of the expression under the square root is $-r^2 2 \sin t \cos t$, which is almost what you already have in the numerator (save for a factor of $-2$). Threrefore, the computation can be continued as
$$\int \limits _0 ^1 -r \left( \int \limits _0 ^{\frac \pi 2} \frac {\frac {\partial} {\partial t} (1 - r^2 \sin^2 t)} {\sqrt {1 - r^2 \sin^2 t}} \Bbb d t \right) \Bbb d r = \int \limits _0 ^1 -r \left( \sqrt {1 - r^2 \sin^2 t} \ \Bigg| _0 ^{\frac \pi 2} \right) \Bbb d r = \int \limits _0 ^1 -r \left( \sqrt {1 - r^2} - 1\right) \Bbb d r = \\
\int \limits _0 ^1 r \left( 1 - \sqrt {1 - r^2} \right) \Bbb d r = \frac {r^2} 2 \Bigg| _0 ^1 + \frac 1 2 \int \limits _0 ^1 (1-r^2)' \sqrt {1-r^2} \ \Bbb d r = \frac 1 2 + \frac 1 2 \frac 2 3 \sqrt{1 - r^2} \Big| _0 ^1 = \frac 1 2 - \frac 1 3 = \frac 1 6.$$
The fact that the result should be $\frac 1 6$ can be also seen by using Fubini's thorem: your integral is
$$\int \limits _0 ^1 x \left( \int \limits _0 ^{\sqrt {1 - x^2}} \frac y {\sqrt {1-y^2}} \ \Bbb d y \right) \ \Bbb d x = \int \limits _0 ^1 x \left( - \sqrt{1 - y^2} \Bigg| _0 ^{\sqrt{1-x^2}} \right) = \int \limits _0 ^1 x ( -x + 1) \ \Bbb d x = - \frac 1 3 + \frac 1 2 = \frac 1 6.$$