Hi everyone! I want to apologise for bumping this old thread, but I was recently going through my old posts and found this again. After spending a couple of days having another go at it, I managed to find a solution on my own :) I've decided to share how I went about it here in case anyone else is curious. I suspect it's a rather inefficient way of approaching the problem, but it's just the way I found that got me the correct answer first.
So, before we get started: when I started redoing this problem I made some changes to how I defined the radii for convince sake. So going forward I'll be using the following:
$$q = r_{q} = \text{ radius of the quarter-circle}, \text{ Point } Q = \text{Centre of quarter-circle with radius }q$$
$$b = r_{b} = \text{ radius of the big circle}, \text{ Point } B = \text{Centre of circle with radius }b$$
$$s = r_{s} = \text{ radii of the two smaller circles}$$ $$\text{ Point } S = \text{Centre of the small circle, with radius s, touching the horizonal}$$
As I mentioned in the original post, getting ##q## in terms of ##b## is pretty straight forward. By drawing a line between ##Q## and ##B##, we can create a right angled triangle with side lengths ##b## and hypotenuse ##\overline{QB}##. From Pythagoras, we can calculate the hypotenuse to be ##\sqrt2 b##. From inspection, we can then see that ##q=( \sqrt2 b + 1 )b##.
Now, onto the new stuff. The first step I took this time around was to construct another right angled triangle with hypotenuse ##\overline{SB}##, like so:
This triangle has a hypotenuse of ##b+s##, a height of ##b-s## and a base length of ##q-b-s-a## (where ##a## is the relatively small length needed to be consider to take the curvature of the quarter circle into account).
From here, we'll now have to start thinking of the image on the Cartesian plane. This means that point ##Q## will be at ##(0,0)##, point ##B## will be at ##(b,b)##, and point ##S## will be at ##(q-s-a, s)##.
The next step is to find a way to express ##a## in terms of ##b## and ##s##. The method I ended up going with was using the points we previously determined to calculate the distance between ##S## and ##B##:
$$\text{Distance} = \sqrt{(x_{2}-x_{1})^2+(y_{2}-y_{1})^2}$$
$$\text{Distance} = \sqrt{(q-s-a-b)^2+(s-b)^2}$$
We know that ##q=( \sqrt2 b + 1 )b## and we know that the hypotenuse has a length of ##b+s##, therefore we can substitute these in:
$$b+s = \sqrt{(( \sqrt2 b + 1 )b-s-a-b)^2+(s-b)^2}$$
After some quick simplifying and rearranging, we can solve for ##a##:
$$a = \sqrt{2}b-2\sqrt{b\cdot s}-s$$
It was around here that I got stuck for a good long while. I couldn't seem to make any more progress trying towards getting a solution algebraically, so instead I decided to graph everything on Desmos and just pick a random value for ##b##. Using this I then picked a value for ##s## which generated a graph of closest fit to what was originally pictured. At this point, I had the following equations for the graphs:
$$x^{2}+y^{2}=(( \sqrt2 b + 1 )b)^{2} |\left\{x\ge0\right\}\left\{y\ge0\right\}$$
$$\left(x-b\right)^{2}+\left(y-b\right)^{2}=b^{2}$$
$$\left(x-b\left(\sqrt{2}+1\right)+s+a\right)^{2}+\left(y-s\right)^{2}=s^{2}$$
I decided to start playing around with changing the horizonal translation of the smaller circle. I first tried to graph ##x^{2}+\left(y-s\right)^{2}=s^{2}##, but that wasn't too interesting. I then shifted it to the right a little with ##\left(x-s\right)^{2}+\left(y-s\right)^{2}=s^{2}## and that's when the eureka moment hit.
Just from inspection, it looked like the distance between ##\left(x-s\right)^{2}+\left(y-s\right)^{2}=s^{2}## and ##\left(x-b\left(\sqrt{2}+1\right)+s+a\right)^{2}+\left(y-s\right)^{2}=s^{2}## would fit an exact number of circles. So I tried adding another circle with the equation ##\left(x-3s\right)^{2}+\left(y-3s\right)^{2}=s^{2}## and it fit perfectly. After adding the remaining circle in, this is what I had discovered:
All the copies of circle ##S## fit in perfectly. Thus, I was able to quickly determine that radius of the quarter circle must be ##8s+a## (I also realized that the equation I'd been using up to this point for circle ##S## was rather messy and can be neatly simplified to ##\left(x-s\right)^{2}+\left(y-7s\right)^{2}=s^{2}##).
This left me with: ##8s+a = q##. Substituting in ##a = \sqrt{2}b-2\sqrt{b\cdot s}-s## and ##q=( \sqrt2 b + 1 )b##, I had:
$$8s + \sqrt{2}b-2\sqrt{b\cdot s}-s = \left(\sqrt{2}+1\right)b$$
And finally, after one last rearrangement...
$$s=\frac{b}{49}\left(4\sqrt{2}+9\right)$$
The last step
Of course, one final step remains: actually calculating the percentage area remaining:
$$\frac{\frac{\pi}{4}q^{2}-\left(2\pi s^{2}+\pi b^{2}\right)}{\frac{\pi}{4}q^{2}}\cdot100$$
$$=\frac{\frac{\pi}{4}\left(\left(\sqrt{2}+1\right)b\right)^{2}-\left(2\pi\left(\frac{b}{49}\left(4\sqrt{2}+9\right)\right)^{2}+\pi b^{2}\right)}{\frac{\pi}{4}\left(\left(\sqrt{2}+1\right)b\right)^{2}}\cdot100$$
$$=19.09\%$$
Conclusion
To anyone still reading, I just wanted to say thank you. This may be very simple and basic to most people on here, but for me it means so much more than you can imagine. I spent a lot of hours attempting to solve this question, and while I very much appreciated all the help I received here, it was all a fair bit beyond my understanding for it to be of any help to me. I'm just so happy I managed to solve this on my own using a method that I was able to figure out on my own and that made sense to me. I hope this is helpful to someone else out there trying and struggling with some geometry :)