Author Topic: FE2  (Read 4044 times)

Victor Ivrii

  • Administrator
  • Elder Member
  • *****
  • Posts: 2607
  • Karma: 0
    • View Profile
    • Personal website of Victor Ivrii
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

  • Administrator
  • Elder Member
  • *****
  • Posts: 2607
  • Karma: 0
    • View Profile
    • Personal website of Victor Ivrii
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
    • View Profile
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.
\begin{equation}
U(x,t)=\int
_{-\infty}^\infty G(x,y,t) F(y)\,dy.
\end{equation}
where $F(y)$ is the function by taking the odd continuation of $f(y)$ and
\begin{equation}
G(x,y,t)=
\frac{1}{2\sqrt{k\pi t}}e^{-\frac{(x-y)^2}{4kt}}
\end{equation}
Note that $f(y)$ is actually an odd function itself and substitute k with $\frac{1}{4}$, we have
\begin{equation}
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
\end{equation}
\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
\begin{equation}
u(x,t) = \frac{xe^\frac{-x^2}{1+t}}{\sqrt{1+t}(1+t)}
\end{equation}

Victor Ivrii

  • Administrator
  • Elder Member
  • *****
  • Posts: 2607
  • Karma: 0
    • View Profile
    • Personal website of Victor Ivrii
Re: FE2
« Reply #3 on: December 15, 2016, 07:42:07 PM »
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)}}$.
« Last Edit: December 15, 2016, 07:45:25 PM by Victor Ivrii »