 Solution

Use the Fourier coefficientsThey are continuous for $0\le y\le 1$ and bounded by $2\Vert u{\Vert }_{\infty }$. On compact subintervals of $0<y<1$, the equation is uniformly elliptic with smooth coefficients. Local elliptic regularity at the flat vertical sides with zero boundary values permits differentiation and integration by parts there. ThusFor completeness, this identity also follows in the distributional sense without assuming boundary derivatives initially. Integrate on $(\epsilon ,\pi -\epsilon )$ against $\sin (nx)$ and a test function of $y$ with support away from zero and one. The interior gradient estimate near a zero Dirichlet side gives $|u_x|\le C{\epsilon }^{-1}\sup |u|$ on a nearby ball of radius comparable to $\epsilon$; continuity and the zero side value make that supremum $o(1)$ uniformly over the support in $y$. Therefore the boundary products $u_x\sin (nx)$ and $u\cos (nx)$ tend to zero. Move the $y$ derivatives to the test function before taking $\epsilon \downarrow 0$. This proves the ordinary differential equation, which makes each $b_n$ smooth in $0<y<1$ and yields the same classical identity.

The Cauchy-Euler differential equation has indicial rootsFor $n\ge 1$ the negative root is strictly negative. Boundedness as $y\downarrow 0$ forces $B_n=0$. Continuity at $y=1$ identifiesBecause $p_n^{+}>n$, the series ${\sum }_{n\ge 1}{A}_n{y}^{p_n^{+}}\sin (nx)$ converges absolutely and uniformly for $0\le y\le \rho <1$. On compact subsets with $y>0$, differentiated series converge as well. At each fixed $0<y<1$, its Fourier coefficients equal those of $u(\cdot ,y)$; completeness of the Fourier sine basis in $L^2(0,\pi )$ and continuity in $x$ make the two functions equal. ThereforeThis is also the requested sum over integers: take all coefficients with $n\le 0$ to be zero. Negative indices give redundant Fourier modes, and the $n=0$ sine vanishes. No pointwise convergence of the ordinary Fourier series on the top boundary is needed.

An important consequence is a forced zero trace at a quadratically degenerate boundary. Indeed the absolute sum is bounded by $2\Vert u{\Vert }_{\infty }{\sum }_{n\ge 1}{y}^n=2\Vert u{\Vert }_{\infty }y/(1-y)$ for $y<1$, which tends to zero as $y\downarrow 0$, uniformly in $x$. Thus $u(x,0)=0$ on the entire lower side. It follows also from $b_n(0)=0$ for every sine coefficient and completeness.