I was hoping that by solving the pde the normal way which is part of the question that it might give a clue as to how the change of variables method should should turn out.

Solving the original PDE I get

. Taking the positive solution of y

which gives

The IC's are u(x,0)=0, therefore

. Let therefore

I am not sure of this is right either...