Method of Characteristics with 3 Partials Use the method of characteristics to solve:
$$yu_x-xu_y+u_z=1, \ u(x,y,0) = x+y.$$
I feel comfortable solving the using the method of characteristics in two dimensions, but am having trouble extending it to three. 
I know I should start by solving the system:
$$\frac{\partial x}{\partial s} = y$$
$$\frac{\partial y }{\partial s } = -x $$
$$\frac{\partial z }{\partial s} = 1$$
$$\frac{\partial u}{\partial s} = 1.$$ 
I am having issues solving the first two equations, as well as knowing how to proceed from there...
 A: $$yu_x-xu_y+u_z=1$$
The change to cylindrical coordinates would simplify the calculus. Nevertheless, we will solve the PDE in Cartesian coo1dinates, just to show that it is not very difficult.
The differential characteristic equations are : 
$\quad \frac{dx}{y}= \frac{dy}{-x}=\frac{dz}{1}=\frac{du}{1}$ 
A first family of characteristic curves comes from $dz=du \quad\to\quad u-z=c_1$
A second family of characteristic curves comes from :$\quad\frac{dx}{y}= \frac{dy}{-x} \quad\to\quad x^2+y^2=c_2$
A third family is a bit more difficult to find : $ \frac{dz}{1}=\frac{dx}{\sqrt{c_2-x^2}} \quad\to\quad \tan^{-1}\left(\frac{x}{\sqrt{c_2-x^2}} \right)-z=c_3 $
On the characteristic curves $c_1,c_2,c_3$ are independent. Elsewhere they are related by an implicit equation : 
$$\Phi\left(u-z\:,\:x^2+y^2\:,\:\tan^{-1}\left(\frac{x}{y} \right)-z\right)=0$$
where $\Phi$ is any differentiable function of three variables.
An equivalant manner to express the relationship is :
$$u-z=F\left(x^2+y^2\:,\:\tan^{-1}\left(\frac{x}{y} \right)-z\right)$$
where $F(X,Y)$ is any differentiable function of two variables.
This is the general solution of the PDE in Cartesian coordinates :
$$u(x,y,z)=z+F\left(x^2+y^2\:,\:\tan^{-1}\left(\frac{x}{y} \right)-z\right)$$
Or in cylindrical coordinates :
$$u=z+F\left(\rho^2\:,\:\tan^{-1}\left(\frac{\cos(\theta)}{\sin(\theta)} \right)-z\right) =z+F\left(\rho^2\:,\:\frac{\pi}{2}-\theta-z\right) $$
Now, we consider the boundary condition :
$$u(x,y,0)=x+y=\rho\left(\cos(\theta)+\sin(\theta)\right)=F\left(\rho^2\:,\:\frac{\pi}{2}-\theta\right)$$
Thus, the function $F(X,Y)$ is determined :
$$F(X,Y)=\sqrt{X}\left(\cos(\frac{\pi}{2}-Y)+\sin(\frac{\pi}{2}-Y) \right)=\sqrt{X}\left(\sin(Y)+\cos(Y) \right)$$
Bringing it back into the above general solution, with $X=\sqrt{x^2+y^2}$ and $Y=\tan^{-1}\left(\frac{x}{y} \right)-z$ , the particular solution of the PDE according to the boundary condition is :
$$u(x,y,z)=z+\sqrt{x^2+y^2}\left[\sin\left(\tan^{-1}\left(\frac{x}{y} \right)-z \right) +\cos \left(\tan^{-1}\left(\frac{x}{y} \right)-z \right) \right]$$
A: Here is a solution that proceeds from the OP's formulation. We take the initial conditions to be $$x(0)=x_0,\;y(0)=y_0,\;z(0)=0,\;u(0)=x_0+y_0.$$ The ODEs for $z(s),u(s)$ may be immediately integrated to obtain $$z(s)=s,\; u(s) = s+x_0+y_0.$$ Note that $s=z$ for all $z$. For $x,y$ we observe that $x''(s) = y'(s)=-x(s)$, and thus \begin{align} x(s) &= x(0) \cos s+x'(0) \sin s\\ &=x_0 \cos s+y_0\sin s,\\ y(s) &= x'(s) \\ &= -x_0\sin s+y_0\cos s.\end{align}
Expressing $s,x_0,y_0$ in terms of $x,y,z$, we obtain
\begin{align}
s&=z\\
x_0 &= x\cos z-y\sin z\\
y_0 &= x\sin z+y\cos z
\end{align}
and therefore
\begin{align}
u(x,y,z)
&=s+x_0+y_0\\
&=(z)+ (x\cos z-y\sin z)+(x\sin z+y\cos z)\\
&=z+(x+y)\cos z+(x-y)\sin z.
\end{align}
