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
: