After taking a topics course in applied mathematics (partial differential equations), I found that there were equations that I should solve since I would later see those equations embedded into other larger-scale equations. This equation was Laplace’s equation (future post). Once I solved this equation, I realized that it becomes a differential operator when acted upon a function of at least two variables. Thus, I could solve equations such as the Schrödinger equation using a three-dimensional laplacian in spherical-polar coordinates (another future post) and the three-dimensional heat equation. I will be solving the latter.

The heat equation initial-boundary-value-problem is therefore

subjected to the boundary conditions

and the initial condition

Now, this solution is not specific to a single thermodynamic system, but rather it is a more general solution in a mathematical context. However, I will be appropriating certain concepts from physics for reasons that are well understood (i.e. that time exists on the interval ). It is nonphysical or nonsensical to speak of negative time.

To start, consider any rectangular prism in which heat flows through the volume, from the origin to the point .At a time , the overall heat of the volume can be regarded to be a function . After a time period has passed, the heat will have traversed to the point from the origin (think of the heat traveling along the diagonal of a cube). The goal is to find the heat as a function of the three spatial components and a single time component. Furthermore, the boundary conditions maintain that at the origin and the final point, the heat vanishes. In other words, these points act as heat sinks. (A sink is a point where energy can leave the system.)

To simplify the notation, we define Thus the initial-boundary-value-problem becomes

where . Also, this definition also reduces the three-dimensional laplacian to a second-order partial derivative of u. The boundary and initial conditions are then

Now, we assume that the solution is a product of eigenfunctions of the form

Taking the respective derivatives and dividing by the assumed form of the solution, we get

Now, Eq.(4) can be equal to three different values, the first of which is zero, but this solution does not help in any way, nor is it physically significant since it produces a trivial solution. The second is , so we can apply to the time dependence equation which then becomes

whose solution is

The third case where allows us to write the spatial equation upon rearrangement as

whose solution is

where . Next we apply the boundary conditions. Let in :

and let in (7.1) to get

Solving for above gives

Substituting into Eq.(7.1) and simplifying gives

To find as many solutions as possible we construct a superposition of solutions of the form

Rewriting as and using the solution for the time dependence, and also let the coefficients form a product equivalent to the indexed coefficient , we arrive at the solution for the heat equation:

Now suppose that in Eq.(11). In this case we get the initial heat distribution given by :