# Derivation of the Heat Equation

### ODEs vs PDEs

I began studying ODEs by solving equations straight off the bat. I started with the simplest one I could think of,$$\frac{dy}{dx}=0,$$and moved onto progressively more difficult ones, like$$\frac{dy}{dx}=x^2.$$While learning, I developed a 'toolkit' of techniques that could solve equations that met very specific criteria. For example, one 'tool' claimed that an equation of the form$$\frac{dy}{dx}+P(x)y=f(x)$$had a solution$$y=\frac{\int e^{\int P(x)\,dx}f(x)\,dx}{e^{\int P(x)\,dx}}.$$I picked up several other tools that could tackle every kind of ODE that you could imagine. However, there were never enough; the set of

*linear*ODEs alone,$$\sum_{k=0}^{n}a_k(x)\frac{d^ky}{dx^k}=f(x),$$has an infinitesimal subset of equations that can be solved by analytic methods, and another one that requires a numerical approach instead. Moreover, when it comes to nonlinear ODEs, what lies before our eyes, more often than not, is a vast ocean full of esoteric sea creatures totally disjoint from all of our equation-solving findings, so to speak.

I thought that I could try an

*a posteriori*approach toward PDEs, and although it could work, it would definitely never provide the necessary theoretical groundwork required to reach higher levels of understanding. To circumvent this issue, I decided to opt for a physically-based introduction to them, as mathematical legends such as Laplace found lots of inspiration from the observational works of Newton and similar applied scientists.

Now, without further due, let me present the

**heat equation**to you!

### The Heat Equation

The first PDE that we are going to study is called the heat equation. Let us derive it by considering a one-dimensional model rod of length $L$ with a perfectly insulated lateral surface. Visually:

**Fourier's law of heat conduction**:$$\varphi(x,t)=-K_0\frac{\partial u}{\partial x}$$Substituting it into the heat equation yields$$c(x)\rho(x)\frac{\partial u}{\partial t}=K_0\frac{\partial^2u}{\partial x^2}+Q(x,t)$$For the sake of simplicity, let us treat $c$, $\rho$, and $K_0$ as constants. Hence,$$\frac{\partial u}{\partial t}=k\frac{\partial^2u}{\partial x^2}+\frac{Q(x,t)}{c\rho}$$where $k=K_0/c\rho$. Letting $Q=0$ yields the canonical heat equation; this is what we are mostly going to work with from here on.

This is how the heat equation is derived. The equalities we presented are physical in nature, but since this is a mathematical blog, we are just going to accept them as true. In the next entry, we are going to go over the

**method of separation of variables**to solve the heat equation with initial and boundary conditions$$\begin{align}\text{IC}&:u(x,0)=f(x)\\\text{BCs}&:\left\{\begin{matrix}u(0,t)&=0\\u(L,t)&=0\end{matrix}\right.\end{align}$$I will talk more about initial and boundary conditions in the next entry as well. Until then, have fun!

## No comments: