Does a solution for this PDE exist ? If so is it possible to find it using Laplace tranform?

131 Views Asked by At

Find the solution for the PDE where $u = u(x,t)$

$$ u_x + xu_t = t $$

$u_x=\frac{\partial{u}}{\partial{x}}$ and $u_t=\frac{\partial{u}}{\partial{t}}$

Given $u(0,t) = 0$ and $u(x,0)=0$

I was asked this in one of my exams to solve using Laplace Transform but it ended up resulting in a solution with $erfi(z)$ whose inverse LT happens to be extremely complex. I'm a fresher taking my first DE class.

2

There are 2 best solutions below

1
On BEST ANSWER

Let me sketch here the solution of the PDE $$ u_x+xu_t=t. \tag{1} $$ using the Laplace transform. Transforming both sides of $(1)$ and using the initial condition $u(x,0)=0$, we obtain $$ U_x(x,s)+sx\,U(x,s)=\frac{1}{s^2}. \tag{2} $$ The ODE $(2)$ can be solved using standard methods; its general solution is $$ U(x,s)=U_0(s)e^{-\frac{1}{2}sx^2}+\int_0^{x}\frac{1}{s^2}e^{\frac{1}{2}s(\xi^2-x^2)}d\xi. \tag{3} $$ Taking the inverse Laplace transform of $(3)$, we obtain $$ u(x,t)=f(t-\frac{1}{2}x^2)\,\Theta(t-\frac{1}{2}x^2)+\int_0^x(t+\frac{1}{2}(\xi^2-x^2))\Theta(t+\frac{1}{2}(\xi^2-x^2))\,d\xi, \tag{4} $$ where $f(t)={\cal L}^{-1}(U_0(s))$ and $\Theta$ is the Heaviside step function. To compute the integral in $(4)$, we notice that the argument of $\Theta$ in the integrand is positive if $\xi<-\sqrt{x^2-2t}$ or $\xi>\sqrt{x^2-2t}$. The first possibility can be ignored, since it is outside the range of integration; the second yields \begin{align} \int_0^x(t+\frac{1}{2}(\xi^2-x^2))\Theta(t+\frac{1}{2}(\xi^2-x^2))\,d\xi &=\int_{\sqrt{x^2-2t}}^x(t+\frac{1}{2}(\xi^2-x^2))\,d\xi \\ &=tx-\frac{1}{3}x^3+\frac{1}{3}(x^2-2t)^{3/2}. \tag{5} \end{align} We notice, however, that the last term in $(5)$ only makes sense if $x^2-2t\geq 0$, so we shall enforce this condition with the Heaviside function $\Theta(x^2-2t)$. Then, combining $(4)$ and $(5)$ we obtain $$ u(x,t)=tx-\frac{1}{3}x^3+\frac{1}{3}(x^2-2t)^{3/2}\,\Theta(x^2-2t)+f(t-\frac{1}{2}x^2)\,\Theta(2t-x^2). \tag{6} $$ We notice that $(6)$ satisfies the initial condition $u(x,0)=0$ for any $f$ (as it should, as this condition was used in the computation of the Laplace transform of the PDE). On the other hand, the function $f$ is determined by the boundary condition $u(0,t)=0$: $$ 0=u(0,t)=f(t), \tag{7} $$ so $(6)$ simplifies to $$ u(x,t)=tx-\frac{1}{3}x^3+\frac{1}{3}(x^2-2t)^{3/2}\,\Theta(x^2-2t). \tag{8} $$ This is precisely the solution found by @JJacquelin in his answer.

0
On

I didn't solve it with the Laplace transform but with the method of characteristics. In order to compare the results my calculus is shown below.

$$u_x+xu_t=t$$

Charpit-Lagrange characteristic ODEs : $$\frac{dx}{1}=\frac{dt}{x}=\frac{du}{t}$$ A first characteristic equation comes from solving $\frac{dx}{1}=\frac{dt}{x}$ : $$t-\frac12 x^2=c_1$$ A second equation comes from solving $\frac{dx}{1}=\frac{du}{t}=\frac{du}{c_1+\frac12 x^2}$ : $$u=c_1 x+\frac16 x^3+c_2$$ The general solution of the PDE on implicite form $c_2=F(c_1)$ is : $$ u=\left(t-\frac12 x^2\right)x+\frac16 x^3+F\left(t-\frac12 x^2 \right)$$ $$\boxed{u(x,t)=tx-\frac13 x^3+F\left(t-\frac12 x^2 \right)}$$ $F$ is an arbitrary function to be determined according to the specified conditions.

First condition $u(0,t)=0=t*0-\frac13 0^3+F\left(t-\frac12 0^2 \right)\quad\implies\quad F(t)=0$. $$u(x,t)=tx-\frac13 x^3 $$

Second condition $u(x,0)=0*x-\frac13 x^3+F\left(0-\frac12 x^2 \right)\quad\implies\quad F\left(-\frac12 x^2 \right)=\frac13 x^3$ .

With $X=-\frac12 x^2\quad\implies\quad F(X)= \frac13 \left(-2X \right)^{3/2}$

The function $F(X)$ is determined. We put it into the general solution :

$$u(x,t)=tx-\frac13 x^3+\frac13 \left(-2(t-\frac12 x^2) \right)^{3/2} $$

So they are two conditions and accordingly two different functions $u(x,t)$ which are solutions of the PDE and satisfy the conditions respectively. This means that the solution is a piecewise function made of the two functions each one valid in a different domain.

The bondary between the two domains starts from $u(0,0)=0$. The equation of the bondary is $$t-\frac12 x^2=0$$

enter image description here

Of course the piecewise function can be expressed on closed form thanks to the Heaviside's step function. This is the result found by Gonçalo : $$u(x,t)=tx-\frac{1}{3}x^3+\frac{1}{3}(x^2-2t)^{3/2}\,\Theta(x^2-2t).$$ The equation of the bondary between the two domains corresponds to the argument of the Heaviside function.