Heat equation on infinite domain

Click For Summary

Discussion Overview

The discussion revolves around the solution of the heat equation in an infinite domain, specifically examining the formulation of the solution using Fourier transforms. Participants explore the representation of the solution and the implications of the initial conditions on the properties of the Fourier transform.

Discussion Character

  • Technical explanation
  • Mathematical reasoning
  • Debate/contested

Main Points Raised

  • One participant presents a solution to the heat equation involving an integral representation and questions the correctness of their approach to express the solution in terms of complex exponentials.
  • Another participant acknowledges the technical correctness of the approach but suggests that evaluating the real integrals directly, as done in the referenced paper, may be more straightforward.
  • A participant expresses uncertainty about how to demonstrate that the function ##C(\omega)## is complex conjugate, questioning its relevance to solving the heat equation.
  • Further discussion leads to an exploration of the relationship between the Fourier transform of the initial condition and the properties of the function ##C(\omega)##, with a proposed demonstration based on the reality of the initial condition function.
  • Participants discuss the implications of the initial condition being real on the equality of the Fourier transforms at negative and positive frequencies.

Areas of Agreement / Disagreement

Participants generally agree on the technical aspects of the Fourier transform approach, but there is disagreement on the necessity of demonstrating the complex conjugate property of ##C(\omega)## in relation to solving the heat equation.

Contextual Notes

Participants note that the initial conditions and the function ##\Psi(x,t)## are real, which influences the properties of the Fourier transform. There is an acknowledgment of the complexity involved in transitioning from the definitions of ##c_1(\omega)## and ##c_2(\omega)## to the desired properties of ##C(\omega)##.

solanojedi
Messages
33
Reaction score
0
Hi everyone, I'm reading this paper about the solution of the heat equation inside an infinite domain: https://ocw.mit.edu/courses/mathema...quations-fall-2006/lecture-notes/fourtran.pdf
1) Please let me know if the following discussion is correct.
The solution ##\Psi(x,t)## reported is the following $$\Psi (x,t)= \int_{0}^{+\infty} [A(\omega) \cos(\omega x)+ B(\omega) sin( \omega x) ] e^{-\kappa \omega^{2}t} d\omega.$$ From second order ODE theory, the ##x## dependent part of the solution can be rewritten to get two complex exponentials $$\Psi (x,t)= \int_{0}^{+\infty} [c_{1}(\omega) e^{i\omega x}+ c_{2}(\omega) e^{-i\omega x} ] e^{-\kappa \omega^{2}t} d\omega.$$ The difference between the two expression should be that ##c_1(\omega)## and ##c_2(\omega)## are complex, while ##A(\omega)## and ##B(\omega)## are real (since the solution ##\Psi(x,t)## is real). Now, if I want to explicit the Fourier Transform inside the solution, I can decompose the solution in the sum of two integrals
$$\Psi (x,t)= \int_{0}^{+\infty} c_{1}(\omega) e^{i\omega x} e^{-\kappa \omega^{2}t} d\omega + \int_{0}^{+\infty} c_{2}(\omega) e^{-i\omega x} e^{-\kappa \omega^{2}t} d\omega$$ and changing the sign on the second integral we get $$\int_{0}^{+\infty} c_{1}(\omega) e^{i\omega x} e^{-\kappa \omega^{2}t} d\omega + \int_{-\infty}^{0} c_{2}(-\omega) e^{i\omega x} e^{-\kappa \omega^{2}t} d\omega = \int_{-\infty}^{+\infty} C(\omega) e^{i\omega x} e^{-\kappa \omega^{2}t} d\omega,$$ where ##C(\omega)=\begin{cases} c_{1}(\omega), ~~\omega \geq 0 \\ c_{2} (-\omega), ~~\omega<0 \end{cases}.## Hence we get the full Fourier transform inside the solution. ##C(\omega)## will be evalued from the initial condition of the problem. Is this approach correct?

2) if this approach is correct, how can I now show that ##C(\omega)## is complex conjugate? From Fourier theory, supposing the signal ##\Psi(x,t)## and the initial conditions are real functions, the transform ##C(\omega)## should be complex conjugate, but in this general situation I'm stucked with the two functions ##c_{1}## and ##c_2## while I should get something like ##C(-\omega)=C^{*}(\omega)##...

Thank you in advance!
 
Physics news on Phys.org
solanojedi said:
Is this approach correct?
Although what you have done is technically correct, I think it would have been more straight forward to just evaluate the real integrals as done in the paper you linked. You have shifted the work from evaluating those integrals to evaluating the form of ##C(\omega)## which will require a Fourier transform of the initial condition. Granted, this would be a clever approach if the initial condition was something simple like a delta function.
solanojedi said:
how can I now show that C(ω)C(ω)C(\omega) is complex conjugate?
Why do you want to show this? I don't see what it has to do with solving the heat equation.
 
Hi NFuller, thank you for your answer!
NFuller said:
Why do you want to show this? I don't see what it has to do with solving the heat equation.
I have to admit that it has nothing to do with it. :)
However, when I got to the point of defining the ##C(\omega)## as ##c_1(\omega)## and ##c_2(\omega)## I questioned myself about how these two different constants could be rearranged to get the situation where ##C(-\omega)=C^{*}(\omega)## (since the initial conditions and hence ##\Psi(x,t)## are real functions) and I wasn't able to solve it.

However, maybe a possible demonstration would be that since the definition of ##C(\omega)## is the Fourier transform of the initial condition ##\Psi(x,0)=f(x) ## we have $$C(\omega)=\int_{-\infty}^{+\infty} f(x) e^{-i\omega x} dx$$ and we can write $$C(-\omega)=\int_{-\infty}^{+\infty} f(x) e^{i\omega x} dx$$ and $$C^{*}(\omega)=\int_{-\infty}^{+\infty} f^{*}(x) e^{i\omega x} dx$$ that are equal having ##f(x)=f^{*}(x)##, i.e. ##f(x)## is real.
It seems to work for me.
 
solanojedi said:
that are equal having f(x)=f∗(x)f(x)=f∗(x)f(x)=f^{*}(x), i.e. f(x)f(x)f(x) is real.
Yes, it seems that if ##f## is real then this works.
 

Similar threads

  • · Replies 6 ·
Replies
6
Views
3K
  • · Replies 17 ·
Replies
17
Views
3K
  • · Replies 1 ·
Replies
1
Views
11K
  • · Replies 6 ·
Replies
6
Views
2K
  • · Replies 1 ·
Replies
1
Views
3K
  • · Replies 1 ·
Replies
1
Views
3K
  • · Replies 8 ·
Replies
8
Views
4K
Replies
1
Views
2K
  • · Replies 3 ·
Replies
3
Views
3K
  • · Replies 4 ·
Replies
4
Views
2K