### Author Topic: FE2  (Read 728 times)

#### Victor Ivrii

• Elder Member
• Posts: 1793
• Karma: 0
##### FE2
« on: December 13, 2016, 07:50:03 PM »
$\newcommand{\erf}{\operatorname{erf}}$
Solve  IVP for the heat equation
\begin{align}
&4u_t -  u_{xx}=0,\qquad &&0 <x<\infty,\; t>0,\label{2-1}\\[2pt]
&u|_{x=0}=0,\\
&u|_{t=0}= f(x)\label{2-2}
\end{align}
with $f(x)=xe^{-x^2}$.

Solution should be expressed  through $\displaystyle{\erf(z)=\frac{2}{\sqrt{\pi}} \int_0^z e^{-z^2}\,dz}$.

#### Victor Ivrii

• Elder Member
• Posts: 1793
• Karma: 0
##### Re: FE2
« Reply #1 on: December 15, 2016, 08:00:35 AM »
Wrong solution removed. Several students got correct solutions, including calculations. I hope that  one of them posts (or I will).

#### Luyu CEN

• Jr. Member
• Posts: 12
• Karma: 0
##### Re: FE2
« Reply #2 on: December 15, 2016, 06:36:09 PM »
This is 1D heat equation on half-line problem with Dirichlet boundary condition.

U(x,t)=\int
_{-\infty}^\infty G(x,y,t) F(y)\,dy.

where $F(y)$ is the function by taking the odd continuation of $f(y)$ and

G(x,y,t)=
\frac{1}{2\sqrt{k\pi t}}e^{-\frac{(x-y)^2}{4kt}}

Note that $f(y)$ is actually an odd function itself and substitute k with $\frac{1}{4}$, we have

U(x,t)=\int
_{-\infty}^\infty G(x,y,t) f(y)\,dy = \int
_{-\infty}^\infty \frac{1}{\sqrt{\pi t}}e^{-\frac{(x-y)^2}{t}} ye^{-y^2}\,dy

\begin{multline*}
\int
_{-\infty}^\infty \frac{1}{\sqrt{\pi t}}e^{-\frac{(x-y)^2}{t}} ye^{-y^2}\,dy
= \frac{1}{\sqrt{\pi t}}e^{-\frac{x^2}{1+t}}\int
_{-\infty}^\infty ye^{-\frac{(1+t)(y-\frac{x}{1+t})^2}{t}} \,dy
\\=\frac{1}{\sqrt{\pi t}}e^{-\frac{x^2}{1+t}}\int
_{-\infty}^\infty \bigl(\frac{x}{1+t}e^{-\frac{(1+t)(y-\frac{x}{1+t})^2}{t}} + (y-\frac{x}{1+t})e^{-\frac{(1+t)(y-\frac{x}{1+t})^2}{t}}\bigr)\,dy
\\= \frac{1}{\sqrt{\pi t}}\sqrt{\frac{t}{1+t}}\frac{\sqrt{\pi}}{2}\frac{x}{1+t}e^{-\frac{x^2}{1+t}}(erf(\infty)-erf(-\infty)) + \frac{1}{\sqrt{\pi t}}e^{-\frac{x^2}{1+t}}\sqrt{\frac{t}{1+t}}e^{-\frac{(1+t)(y-\frac{x}{1+t})^2}{t}}|_{-\infty}^{\infty} \\= \frac{xe^\frac{-x^2}{1+t}}{\sqrt{1+t}(1+t)}
\end{multline*}
Therefore, the solution is given by

u(x,t) = \frac{xe^\frac{-x^2}{1+t}}{\sqrt{1+t}(1+t)}

#### Victor Ivrii

Indeed, $u=\frac{x}{(t+1)^{\frac{3}{2}}}e^{-\frac{x^2}{(t+1)}}$ and there is a way to get it even almost without calculations. First, as mentioned, since $f(x)=xe^{-x^2}$ is an odd function, it is sufficient to solve Cauchy problem.
We know that $v(x,t)=t^{-\frac{1}{2}}e^{-\frac{x^2}{t}}$ satisfies our equation (\ref{2-1}).
Plugging $t:=t+1$ we conclude that $w(x,t)=(t+1)^{-\frac{1}{2}} e^{-\frac{x^2}{t+1}}$ also satisfies (\ref{2-1}) (since coefficients do not depend on $t$). Further, $w(x,0)=g(x):=e^{-x^2}$.
Next, $-\frac{1}{2}w_x$ also satisfies (\ref{2-1}) (since coefficients do not depend on $x$) and $-\frac{1}{2}w_x(x,0)=-\frac{1}{2}g'(x)$ which is $f(x)$. Bingo! $u=-\frac{1}{2}w_x=\frac{x}{(t+1)^{\frac{3}{2}}}e^{-\frac{x^2}{(t+1)}}$.