The heat equation is a partial differential equation describing the distribution of heat over time. In one spatial dimension, we denote as the temperature which obeys the relation



where is called the diffusion coefficient. Problems related to partial differential equations are typically supplemented with initial conditions and certain boundary conditions. In this article, we go over how to solve the heat equation using Fourier transforms. It's recommended that you be familiar with their properties before proceeding.

Preliminaries

  • In this article, we use the following convention for the Fourier transform and its inverse. Note that the Fourier transforms are being applied to real space, not time.
  • Diffusion problems frequently encounter the error function, a special function defined as the antiderivative of the Gaussian. The normalization factor is such that the function has a range of

Steps

  1. 1
    Transform the equation into Fourier space. In this section, we outline the steps to finding the fundamental solution, a term whose name we will shortly come to understand.
    • Taking the Fourier transform of a derivative of order is the same as multiplication by Because the Fourier integral is independent of we can pull the derivative out of the integral and write
  2. 2
    Solve the resulting ordinary differential equation.
    • Solutions are decaying exponentials in The constant term is the initial conditions in Fourier space, denoted by
    Advertisement
  3. 3
    Transform back into real space.
    • The property of the Fourier transform we take advantage of here is convolution: multiplication in Fourier space corresponds to convolution in real space.
    • The term is the fundamental solution sought after, also known as the heat kernel. It is the solution to the heat equation given initial conditions of a point source, the Dirac delta function, for the delta function is the identity operator of convolution.
  4. 4
    Evaluate the inverse Fourier integral. The inverse Fourier transform here is simply the integral of a Gaussian. We evaluate it by completing the square. If one looks up the Fourier transform of a Gaussian in a table, then one may use the dilation property to evaluate instead.
    • This is the well-known fundamental solution to the heat equation. From here, we need only substitute initial conditions and evaluate the resulting convolution integral to obtain a solution.
  5. 5
    Find given initial conditions of the rectangular function.
    • The function written below is known by other names, including the gate function, or the unit pulse.
    • Now, we simply substitute this function into the convolution integral. Here, the form is particularly simple.
    • In the last step, we make use of the fact that
    • A plot of this function over time above shows that the "sharpness" of the function diminishes over time, eventually tending towards an equilibrium solution. This is what the heat equation is supposed to do - it says that the time rate of change of is proportional to the curvature of as denoted by the spatial second derivative, so quantities obeying the heat equation will tend to smooth themselves out over time. The steady-state solution where will therefore obey Laplace's equation.
    • Setting the initial conditions are plotted in blue, while is being plotted for values and for orange, green, and red plots, respectively.
  6. 6
    Find given initial conditions of the ramp function over a restricted domain. Specifically, where denotes the Heaviside step function. This is the ramp function over the domain Its solution is slightly more complicated. To find we need to split the integral into two pieces.




    • We see that the second integral is different from the first only by the lower boundary. Therefore, we will only detail the process for the first integral only. We make a substitution which splits this integral into two integrals which we may easily evaluate. Note that below refers to a substitution variable, not to the temperature density.



    • The second integral is found by a similar process.



    • Hence, our final answer is written as follows.



    • Setting the initial conditions are plotted in blue, while is being plotted for values and for orange, green, and red plots, respectively.
  7. Advertisement

Discussion

  • The heat equation we have been dealing with is homogeneous - that is, there is no source term on the right that generates heat.
    • We can show that the total heat is conserved for solutions obeying the homogeneous heat equation. That is, the relation below must be satisfied.
    • We simply substitute the convolution integral, exchange the order of integration, and then recognize that the integral in is simply 1.
    • Because is simply a dummy variable, we have shown that the total heat is conserved, as it should be.
  • A word should be said about the physicality of the solutions that we have obtained.
    • The initial conditions describe functions that have compact support. Intuitively, this means that the functions map to non-zero values within some limited domain, and map to zero elsewhere. This is a reasonable description for most materials.
    • However, the solutions are defined for and since the error function is a smooth function over the real line, does not have compact support, implying that the function takes on non-zero values everywhere. We know physically that heat transfer is limited by at least the speed of light, so the model cannot be applied when such conditions become a significant factor. Nevertheless, the solution does decay exponentially, so we may treat the "non-local" regions as an approximation to be neglected.


About This Article

wikiHow is a “wiki,” similar to Wikipedia, which means that many of our articles are co-written by multiple authors. To create this article, volunteer authors worked to edit and improve it over time. This article has been viewed 24,539 times.
47 votes - 88%
Co-authors: 2
Updated: December 2, 2021
Views: 24,539
Advertisement