Let $R>0$ be a large constant, such that for any $x \in \Omega$, $\Omega \subset B_R(x)$. Consider the following problem in $\mathbb{R}^2$: $$ -\varepsilon^2 \Delta u=1_{\{u>a\}} \text { in }\, B_R(0) \quad u=0 \text { on }\, \partial B_R(0), $$ where $a>0$ is a constant. Then, the above equation has a solution $U_{\varepsilon , a}$, which can be written as $$ U_{ \varepsilon, a}(y)= \left\{\begin{aligned}& a+\frac{1}{4 \varepsilon^2}(s_{ \varepsilon}^2-|y|^2), && |y| \leq s_{ \varepsilon}, \\ & a \ln \frac{|y|}{R} / \ln \frac{s_{ \varepsilon}}{R}, && s_{ \varepsilon} \leq|y| \leq R,\end{aligned}\right. $$ where $s_{\varepsilon}$ is the constant, such that $U_{\varepsilon, a} \in C^1(B_R(0))$.
So, $s_{\varepsilon}$ satisfies $$ -\frac{s_{\varepsilon }}{2 \varepsilon ^2}=\frac{a}{s_{\varepsilon} \ln \frac{s_{\varepsilon}}{R}} . $$
From $$ \frac{s_{\varepsilon } \sqrt{\ln \frac{R}{s_{\varepsilon }}}}{\varepsilon }=\sqrt{2 a}, $$ we see that if $\varepsilon>0$ is small, the above equation is uniquely solvable for $s_{\varepsilon}>0$ small. Moreover, we have the following expansion for $s_{\varepsilon}$ : $$ s_{\varepsilon }=\frac{\sqrt{2 a} \varepsilon}{\sqrt{|\ln \varepsilon |}}\Big(1+O\Big(\frac{\ln |\ln \varepsilon |}{|\ln \varepsilon |}\Big)\Big) . $$ Above are some calculations in the paper [Planar vortex patch problem in incompressible steady flow]. But I can not find out the last estimate. Any comments are welcome.