15.2.A. Heat Equation

Form: ∂𝑢/∂𝑡 = 𝛼² ∂²𝑢/∂𝑥², 0 < 𝑥 < 𝐿, 𝑡 > 0

The heat equation is a partial differential equation (PDE) that describes how heat diffuses through a given region over time. In this case, the equation describes the behavior of temperature (u) as a function of both space (x) and time (t). Here, ∂u/∂t denotes the rate of change of temperature with respect to time, while ∂²u/∂x² represents the curvature of the temperature distribution with respect to position. The constant α² refers to the thermal diffusivity of the material, quantifying how quickly heat moves through the medium.

Imagine placing a metal rod in a fire on one end. Initially, only the end in the fire is hot. Over time, the heat travels along the rod, resulting in an increase in temperature at points further along the rod. This process of heat traveling is described mathematically by the heat equation.

Boundary Conditions: 𝑢(0,𝑡) = 0, 𝑢(𝐿,𝑡) = 0

Boundary conditions are crucial in solving PDEs because they define how the solution behaves at the boundaries of the domain. In this heat equation, the boundary conditions indicate that at position x=0 and x=L, the temperature u is fixed at zero for all times t. This means the ends of the rod are kept at a constant temperature (which is zero in this case), often representing a perfect heat sink.

Think of a metal rod that is held at both ends with ice. The ice keeps the ends of the rod at a temperature of 0°C. As you heat the middle of the rod, the heat will travel towards the ends, but they will always remain at 0°C because of the ice, illustrating the boundary conditions.

Initial Condition: 𝑢(𝑥,0) = 𝑓(𝑥)

The initial condition for the heat equation specifies the temperature distribution along the entire rod at the initial time (t=0). The term u(x,0) = f(x) means that the temperature function u at time t=0 is determined by a specific function f, which describes the initial temperature profile along the length of the rod. This sets the starting state from which the evolution of temperature is calculated.

Continuing with the metal rod analogy, imagine that you initially heat the middle portion of the rod to a high temperature while the ends are at 0°C due to the ice. The initial condition describes this setup, often varying based on how and where you apply heat to the rod.

Solution Approach:

1. Separation of Variables: Assume 𝑢(𝑥,𝑡) = 𝑋(𝑥)𝑇(𝑡)

The method of separation of variables involves assuming that the solution can be expressed as a product of two functions, one depending only on x (spatial function X(x)) and the other only on t (temporal function T(t)). This helps transform the PDE into two ordinary differential equations (ODEs), making them easier to solve individually.

Think of baking. If you want to understand how to bake a cake (the final outcome), you can separate the process into two tasks: preparing the batter (the spatial part) and baking it (the time part). By focusing on one task at a time, you simplify the process and clarify your next step.

2. Substituting in the PDE and separating variables gives: 1/𝛼²𝑇(𝑡) (𝑑𝑇/𝑑𝑡) = (1/𝑋(𝑥))(𝑑²𝑋/𝑑𝑥²) = −𝜆

After separating variables, we end up with the relationship where both sides of the equation equal a constant, denoted by -λ. This results in two separate equations: one for T(t) and one for X(x). Each of these equations can be analyzed and solved using techniques appropriate to ordinary differential equations, leading to easier handling of the complex problem originally described.

Consider trying to find the height and weight of several apples on different trees. By treating each tree independently, you might measure each tree's height and weight separately, allowing for more straightforward comparisons.

3. Solve two ODEs: 𝑑²𝑋/𝑑𝑥² + 𝜆𝑋 = 0 𝑑𝑇/𝑑𝑡 + 𝛼²𝜆𝑇 = 0

Each ODE can be solved using techniques appropriate to their forms. The first equation for X(x) is a second-order linear differential equation typical in problems involving oscillations, while the second equation for T(t) is also linear and shows exponential decay related to the constant λ. The solutions yield functions that describe how temperature is distributed in space and how it evolves in time.

If we liken solving these equations to tuning musical instruments, the first equation would be like figuring out how to adjust the pitch of a string by altering its length and tension (X(x)), while the second would be similar to how the sound fades away over time (T(t)). Both need their own methods to find the perfect results.

4. General solution: 𝑢(𝑥,𝑡) = ∑(𝐵 sin(𝑛𝜋𝑥/𝐿)e^(-𝛼²(𝑛𝜋/𝐿)²𝑡)) Where 𝐵 is obtained from the Fourier sine series of 𝑓(𝑥).

The general solution to the heat equation is expressed as an infinite series involving the sine function and an exponential decay term. The term B refers to coefficients derived from the initial condition function f(x) using Fourier sine series, which reflect how the initial temperature distribution contributes to the evolving solution over time.

Imagine you are creating music by layering different sounds (like sine waves) to create a rich harmony. Each sound’s volume and presence depend on how it fits with the initial melody you composed, which is represented here by the initial temperature distribution.

--