One-Dimensional Heat Equation: Fourier Series and Boundary Conditions

Overview

This video examines solving the one dimensional heat equation for a rod, showing how the temperature evolves in time and space under boundary and initial conditions. It starts with a simple sine temperature profile and explains why sums of sine waves remain solutions, paving the way to Fourier series as a powerful tool.

Key insights

• Heat equation links time change to spatial curvature via the second derivative
• Simple sine profiles yield exponential decay in time
• Boundary conditions enforce flat ends, leading to eigenfunctions with quantized frequencies
• Any initial temperature distribution can be expressed as a sum of these eigenfunctions

Introduction

The discussion centers on the heat equation in one dimension, modeling heat flow along a rod. The PDE relates how the temperature at a point changes in time to how curved the temperature distribution is in space, via the second spatial derivative. Beyond the equation itself, boundary and initial conditions shape the actual evolution, and Fourier's method provides a way to pick the specific solution that matches the initial state.

Why sine waves are a natural starting point

To keep constants manageable, the video first tests a temperature profile at time zero of the form sin x. The second derivative of sin x is -sin x, so the PDE reduces to a simple scaling in time. This leads to a family of solutions where the temperature is a sine wave in space and decays exponentially in time, preserving the same spatial shape while the amplitude decreases. This insight points to exponentials in time as the time evolution mechanism for these modes.

Boundary conditions and the real evolution

In a physical rod with no heat flow across its ends, the endpoints must be flat with respect to the rod, corresponding to a zero spatial slope at the boundaries. This Neumann boundary condition means the spatial derivative ∂u/∂x vanishes at x = 0 and x = L for all t > 0. A naive sine profile may not satisfy this boundary restriction, so the analysis proceeds to modify the spatial form to meet the ends while still solving the PDE.

Introducing harmonics through frequency adjustment

The right hand side of the PDE introduces a dependence on the square of the spatial frequency. Increasing the spatial frequency sharpens curvature and accelerates decay in time, producing a quadratic relationship between frequency and decay rate. The lowest nontrivial frequency that keeps the right boundary flat occurs when the spatial input spans up to PI divided by the rod length, giving a base frequency. Higher harmonics are obtained by whole number multiples of this base frequency, with the constant function (frequency zero) also forming a valid solution under the boundary conditions.

Building the infinite family of solutions

With the boundary condition in place, an infinite family of functions is available that satisfy both the PDE and the end conditions. Each member looks like a sine or cosine wave in space and an exponential in time, and their sums are also solutions due to linearity. The power of this approach is that any initial temperature distribution can be approximated by a (potentially infinite) sum of these eigenfunctions, each decaying at a rate determined by its frequency. In a future discussion, this leads to a general solution constructed by Fourier series, combining all the relevant harmonics to fit the initial state precisely.

Takeaways and the road ahead

The key lesson is that solving the heat equation often means reducing a complex initial state to a sum of simple, well-behaved pieces whose evolution is easy to track. Boundary conditions and the interior PDE work together to select the appropriate set of modes, and Fourier series provides a practical recipe for handling arbitrary initial conditions. The upcoming material will show how to assemble these pieces into a complete solution for any initial temperature profile on the rod.