Heat/Diffusion Equation on the whole real line

  • Thread starter Thread starter eventhorizonof
  • Start date Start date
  • Tags Tags
    Line
Click For Summary
SUMMARY

The discussion focuses on solving the heat equation, specifically the partial differential equation defined by u_t = u_xx, with boundary conditions u(x,0) = exp[-((x-14)^2)/4]. The general solution for the whole real line is presented as u(x,t) = 1 / sqrt(4kt) ∫ exp[-((x-y)^2)/4kt] * PHI(y) dy. It is established that while the integral cannot be avoided, a substitution can simplify the computation by transforming the argument of the exponential, ultimately leading to a solvable integral of the form ∫ e^{-p^2} dp, which equals √π.

PREREQUISITES
  • Understanding of partial differential equations, specifically the heat equation.
  • Familiarity with integral calculus and Gaussian integrals.
  • Knowledge of boundary conditions in mathematical physics.
  • Experience with substitution techniques in integrals.
NEXT STEPS
  • Study the derivation of the heat equation and its general solutions.
  • Learn about Gaussian integrals and their applications in solving PDEs.
  • Explore substitution methods for simplifying integrals in mathematical analysis.
  • Investigate the implications of boundary conditions on the solutions of PDEs.
USEFUL FOR

Mathematicians, physics students, and engineers dealing with heat transfer problems, as well as anyone interested in advanced calculus and differential equations.

eventhorizonof

Homework Statement



boundary conditions:

u(x,0) = exp[-((x-14)^2)/4

Homework Equations



u_t = u_xx

x [tex]\in[/tex] R


The Attempt at a Solution



the problem says there is a way for this partial differential equation to be solved without computing any integrals. i know the general solution for the whole line is:

u(x,t) = 1 / sqrt(4kt) [tex]\int[/tex] exp[-((x-y)^2)/4kt] * PHI(y) dy


where PHI(y) is the boundary conditions.

is there some kind of substitution or trick to be used here to not have to compute any integral?

Thanks.
 
Physics news on Phys.org
I'm under the impression that you are mixing a few things on the notation. You use k in the resulting integral, but that is nowhere defined. I know that sometimes the partial differential equation is written as:

[tex]a^2 \cdot \frac{\partial^2 u}{\partial x^2} = \frac{\partial u}{\partial t}[/tex]

or:

[tex]k \cdot \frac{\partial^2 u}{\partial x^2} = \frac{\partial u}{\partial t}[/tex]

In your case a² or k is equal to 1. However you still have it in the integral. There is also something wrong with a pi. I assume you mean that the solution has to be:

[tex]u(x,t)=\frac{1}{2a\sqrt{\pi t}} \cdot \int_{-\infty}^{\infty} f(y) e^{-\frac{(x-y)^2}{4a^2t}} dy[/tex]

or in your case:

[tex]u(x,t)=\frac{1}{\sqrt{4 \pi t}} \cdot \int_{-\infty}^{\infty} f(y) e^{-\frac{(x-y)^2}{4t}} dy[/tex]

Now as for the question, You can't do without the integral. However you can avoid having to calculate it explicitly by changing the argument of the exponential. This means that you need to write it as:

[tex]u(x,t)=\frac{1}{\sqrt{4 \pi t}} \cdot \int_{-\infty}^{\infty} e^{-\frac{(y-14)^2}{4}} e^{-\frac{(x-y)^2}{4t}} dy[/tex]

And now transform the argument of the exponential to something like:

[tex]A+B\cdot (y-C)^2[/tex]

The part with A can be put outside the integral and you need to use the following substitution:

[tex]\sqrt{B}(y-C)=p[/tex]

Giving you a final integral as:

[tex]\int_{-\infty}^{\infty}e^{-p^2}dp[/tex]

Which is equal to [tex]\sqrt{\pi}[/tex]. So you do not need to calculate it, just rewriting it in another form.
 

Similar threads

  • · Replies 3 ·
Replies
3
Views
2K
  • · Replies 6 ·
Replies
6
Views
7K
  • · Replies 4 ·
Replies
4
Views
2K
Replies
2
Views
2K
Replies
0
Views
1K
  • · Replies 1 ·
Replies
1
Views
2K
Replies
5
Views
2K
  • · Replies 12 ·
Replies
12
Views
2K
  • · Replies 5 ·
Replies
5
Views
1K
Replies
12
Views
5K