Find the Steady State Solution of the Heat Equation
Suppose we have an insulated wire, a plate, or a 3-dimensional object. We apply certain fixed temperatures on the ends of the wire, the edges of the plate, or on all sides of the 3-dimensional object. We wish to find out what is the steady state temperature distribution. That is, we wish to know what will be the temperature after long enough period of time.
We are really looking for a solution to the heat equation that is not dependent on time. Let us first solve the problem in one space variable. We are looking for a function \(u\) that satisfies
\[u_t=ku_{xx},\]
but such that \(u_t=0\) for all \(x\) and \(t\). Hence, we are looking for a function of \(x\) alone that satisfies \(u_{xx}=0\). It is easy to solve this equation by integration and we see that \(u=Ax+B\) for some constants \(A\) and \(B\).
Suppose we have an insulated wire, and we apply constant temperature \(T_1\) at one end (say where \(x=0\)) and \(T_2\) on the other end (at \(x=L\) where \(L\) is the length of the wire). Then our steady state solution is
\[u(x)= \dfrac{T_2-T_1}{L}x+T_1.\]
This solution agrees with our common sense intuition with how the heat should be distributed in the wire. So in one dimension, the steady state solutions are basically just straight lines.
Things are more complicated in two or more space dimensions. Let us restrict to two space dimensions for simplicity. The heat equation in two space variables is
\[\label{eq:3} u_t=k(u_{xx}+u_{yy}), \]
or more commonly written as \( u_t=k \Delta u \) or \( u_t=k \nabla^2 u \). Here the \( \Delta \) and \( \nabla^2 \) symbols mean \( \dfrac{\partial^2}{\partial x^2}+ \dfrac{\partial^2}{\partial y^2}\). We will use \( \Delta\) from now on. The reason for using such a notation is that you can define \( \Delta\) to be the right thing for any number of space dimensions and then the heat equation is always \(u_t=k \Delta u\). The operator \( \Delta\) is called the Laplacian.
OK, now that we have notation out of the way, let us see what does an equation for the steady state solution look like. We are looking for a solution to Equation \(\eqref{eq:3}\) that does not depend on \(t\), or in other words \(u_t=0\). Hence we are looking for a function \(u(x,y)\) such that
\[ \Delta u = u_{xx}+u_{yy}=0.\]
This equation is called the Laplace equation\(^{1}\). Solutions to the Laplace equation are called harmonic functions and have many nice properties and applications far beyond the steady state heat problem.
Harmonic functions in two variables are no longer just linear (plane graphs). For example, you can check that the functions \(x^2-y^2\) and \(xy\) are harmonic. However, if you remember your multi-variable calculus we note that if \(u_{xx}\) is positive, \(u\) is concave up in the \(x\) direction, then \(u_{yy}\) must be negative and \(u\) must be concave down in the \(y\) direction. Therefore, a harmonic function can never have any "hilltop" or "valley" on the graph. This observation is consistent with our intuitive idea of steady state heat distribution; the hottest or coldest spot will not be inside.
Commonly the Laplace equation is part of a so-called Dirichlet problem\(^{2}\). That is, we have a region in the \(xy\)-plane and we specify certain values along the boundaries of the region. We then try to find a solution \(u\) defined on this region such that \(u\) agrees with the values we specified on the boundary.
For simplicity, we consider a rectangular region. Also for simplicity we specify boundary values to be zero at 3 of the four edges and only specify an arbitrary function at one edge. As we still have the principle of superposition, we can use this simpler solution to derive the general solution for arbitrary boundary values by solving 4 different problems, one for each edge, and adding those solutions together. This setup is left as an exercise.
We wish to solve the following problem. Let \(h\) and \(w\) be the height and width of our rectangle, with one corner at the origin and lying in the first quadrant.
\[\begin{align} & \Delta u = 0 , & & \label{dirich:eq1} \\ & u(0,y) = 0 & & \text{for } 0 < y < h,\label{dirich:eq2} \\ & u(x,h) = 0 & & \text{for } 0 < x < w,\label{dirich:eq3} \\ & u(w,y) = 0 & & \text{for } 0 < y < h,\label{dirich:eq4} \\ & u(x,0) = f(x) & & \text{for } 0 < x < w.\label{dirich:eq5}\end{align}\]
Figure \(\PageIndex{1}\)
The method we apply is separation of variables. Again, we will come up with enough building-block solutions satisfying all the homogeneous boundary conditions (all conditions except \(\eqref{dirich:eq5}\)). We notice that superposition still works for the equation and all the homogeneous conditions. Therefore, we can use the Fourier series for \(f(x)\) to solve the problem as before.
We try \(u(x,y)=X(x)Y(y)\). We plug \(u\) into the equation to get
\[X''Y+XY''=0.\]
We put the \(X\)s on one side and the \(Y\)s on the other to get
\[ - \dfrac{X''}{X}= \dfrac{Y''}{Y}.\]
The left hand side only depends on \(x\) and the right hand side only depends on \(y\). Therefore, there is some constant \( \lambda \) such that \( \lambda = \dfrac{-X''}{X}= \dfrac{Y''}{Y}\). And we get two equations
\[ X'' + \lambda X = 0, \\ Y'' - \lambda Y=0.\]
Furthermore, the homogeneous boundary conditions imply that \(X(0)=X(w)=0\) and \(Y(h)=0\). Taking the equation for \(X\) we have already seen that we have a nontrivial solution if and only if \( \lambda= \lambda_n = \dfrac{n^2 \pi^2}{w^2}\) and the solution is a multiple of
\[ X_n(x) = \sin \left( \dfrac{n \pi}{w}x \right). \]
For these given \( \lambda_n \), the general solution for \(Y\) (one for each \(n\)) is
\[\label{eq:14} Y_n(y)=A_n \cosh \left( \dfrac{n \pi}{w}y \right) + B_n \sinh \left( \dfrac{n \pi}{w}y \right). \]
We only have one condition on \(Y_n\) and hence we can pick one of \(A_n\) or \(B_n\) to be something convenient. It will be useful to have \(Y_n(0)=1\), so we let \(A_n=1\). Setting \(Y_n(h)=0\) and solving for \(B_n\) we get that
\[ B_n = \dfrac{ - \cosh \left( \dfrac{n \pi h}{w} \right)}{ \sinh \left( \dfrac{n \pi h}{w} \right)}.\]
After we plug the \(A_n\) and \(B_n\) we into \(\eqref{eq:14}\) and simplify, we find
\[ Y_n(y) = \dfrac{ \sinh \left( \dfrac{n \pi (h-y)}{w} \right)}{ \sinh \left( \dfrac{n \pi h}{w} \right)}.\]
We define \(u_n(x,y)=X_n(x)Y_n(y)\). And note that \(u_n\) satisfies \(\eqref{dirich:eq1}\) - \(\eqref{dirich:eq4}\).
Observe that
\[ u_n(x,0)= X_n(x)Y_n(0) = \sin \left( \dfrac{n \pi}{w}x \right). \]
Suppose
\[ f(x)= \sum_{n=1}^{ \infty}b_n \sin \left( \dfrac{n \pi x}{w} \right).\]
Then we get a solution of \(\eqref{dirich:eq1}\) - \(\eqref{dirich:eq5}\) of the following form.
\[ u(x,y)= \sum_{n=1}^{ \infty}b_n u_n(x,y) = \sum_{n=1}^{ \infty}b_n \sin \left( \dfrac{n \pi }{w}x \right) \left( \dfrac{ \sinh \left( \dfrac{n \pi (h-y)}{w} \right)}{ \sinh \left( \dfrac{n \pi h}{w} \right)} \right).\]
As \(u_n\) satisfies Equation \(\eqref{dirich:eq1}\) - \(\eqref{dirich:eq4}\) and any linear combination (finite or infinite) of \(u_n\) must also satisfy \(\eqref{dirich:eq1}\) - \(\eqref{dirich:eq4}\), we see that \(u\) must satisfy Equations \(\eqref{dirich:eq1}\) - \(\eqref{dirich:eq4}\). By plugging in \(y=0\) it is easy to see that \(u\) satisfies \(\eqref{dirich:eq5}\) as well.
Example \(\PageIndex{1}\)
Suppose that we take \(w=h= \pi\) and we let \(f(x)= \pi\). We compute the sine series for the function \( \pi\)(we will get the square wave). We find that for \(0<x< \pi\) we have
\[ f(x)= \sum_ {\underset{n~ {\rm{odd}} }{n=1}}^{ \infty} \dfrac{4}{n} \sin(n x).\]
Therefore the solution \(u(x,y)\), see Figure \(\PageIndex{2}\), to the corresponding Dirichlet problem is given as
\[ u(x,y)= \sum_ {\underset{n~ {\rm{odd}} }{n=1}}^{ \infty} \dfrac{4}{n} \sin(n x) \left( \dfrac{ \sinh(n( \pi-y)) }{ \sinh(n \pi) } \right).\]
Figure \(\PageIndex{2}\): Steady state temperature of a square plate, three sides held at zero and one side held at \(\pi\).
This scenario corresponds to the steady state temperature on a square plate of width \( \pi\) with 3 sides held at 0 degrees and one side held at \( \pi\) degrees. If we have arbitrary initial data on all sides, then we solve four problems, each using one piece of nonhomogeneous data. Then we use the principle of superposition to add up all four solutions to have a solution to the original problem.
A different way to visualize solutions of the Laplace equation is to take a wire and bend it so that it corresponds to the graph of the temperature above the boundary of your region. Cut a rubber sheet in the shape of your region—a square in our case—and stretch it fixing the edges of the sheet to the wire. The rubber sheet is a good approximation of the graph of the solution to the Laplace equation with the given boundary data.
Footnotes
[1] Named after the French mathematician Pierre-Simon, marquis de Laplace (1749–1827).
[2] Named after the German mathematician Johann Peter Gustav Lejeune Dirichlet (1805–1859).
Find the Steady State Solution of the Heat Equation
Source: https://math.libretexts.org/Bookshelves/Differential_Equations/Book:_Differential_Equations_for_Engineers_%28Lebl%29/4:_Fourier_series_and_PDEs/4.09:_Steady_state_temperature_and_the_Laplacian