I want to know the analytical solution of a transient heat equation in a 2D square with inhomogeneous Neumann Boundary.

The transient heat conduction equation in a 2D square cavity :


and the boundary are:

$$\cases{T(0,y)=T_1\\T(L,y)=T_0\\\frac{\partial T(x,0)}{\partial y}=a\\\frac{\partial T(x,L)}{\partial y}=a}$$

The initial Temperature inside is $T_0$

If there is no analytical solution, then what is the solution of steady heat conduction (with $\nabla^2T=0$)?

Really appreciate anyone who can help!

  • 1
    $\begingroup$ $T$ depends on 3 variables: $T(t,x,y)$. That may be the cause of the confusion. $\endgroup$ Jun 21 '18 at 23:13
  • $\begingroup$ +1 for cute fox. No seriously the problem is interesting too. $\endgroup$ Jun 23 '18 at 9:09

This is very involved problem. To solve, you'll need to break the solution into successively smaller pieces. First, separate the steady-state and transient solutions, then split up the boundary conditions in order to use separation of variables.

To start off, I'm going to label

$$ T(x,y,t) = u(x,y,t) + v(x,y) $$

where $v(x,y)$ is the time-independent, steady-state solution, and $u(x,y,t)$ is the decaying, transient solution.

The steady-state solution should satisfy $\nabla^2v = 0$ and all the boundary conditions as listed. Since all boundaries are inhomogeneous, we need so split it up further

$$ v(x,y) = v_1(x,y) + v_2(x,y) $$

such that

\begin{matrix} \begin{cases} \nabla^2 v_1 = 0 \\ \\ v_1(0,y) = T_1,\ v_1(L,y) = T_0 \\ \\ \dfrac{\partial}{\partial y}v_1(x,0) = \dfrac{\partial}{\partial y}v_1(x,L) = 0 \end{cases} &&& \begin{cases} \nabla^2v_2 = 0 \\ \\ v_2(0,y) = v_2(L,y) = 0 \\ \\ \dfrac{\partial}{\partial y}v_2(x,0) = \dfrac{\partial}{\partial y}v_2(x,L) = a \end{cases} \end{matrix}

The homogeneous boundaries allow us to use separation of variables to solve each individual problem

  • The first one is easier since you can intuitively guess that it's constant in $y$ and linear in $x$. This turns out to be

$$ v_1(x,y) = T_0\frac{x}{L} + T_1\frac{L-x}{L} = \frac{(T_0-T_1)x}{L} + T_1 $$

  • For the second problem, the homogeneous boundary condition on $x$ returns a series solution of the form

$$ v_2(x,y) = \sum_{n=1}^\infty \sin\left(\frac{n\pi x}{L}\right)\left[A_n\cosh\left(\frac{n\pi y}{L}\right) + B_n \cosh\left(\frac{n\pi(L-y)}{L}\right)\right] $$

Then, applying the remaining boundary conditions on $y$ will give

$$ A_n = -B_n = \begin{cases} \dfrac{4aL}{n^2\pi^2\sinh(n\pi)}, & n \text{ odd} \\ 0, & n \text{ even} \end{cases} $$

For a more thorough explanation of why I'm using hyperbolic functions instead of exponentials, check out my answer for this similar problem. Short answer: It makes the math easier

The remaining transient solution is homogeneous on all boundaries and has initial conditions that cancel out the steady-state:

\begin{cases} \dfrac{\partial u}{\partial t} = \nabla^2 u \\ \\ u(0,y,t) = u(L,y,t) = \dfrac{\partial}{\partial y}u(x,0,t) = \dfrac{\partial}{\partial y}u(x,L,t) = 0 \\ \\ u(x,y,0) = T_0 - v(x,y) \end{cases}

Applying separation of variables once more and matching the homogeneous boundaries, we obtain

$$ u(x,y,t) = \sum_{n,m} c_{n,m} \sin\left(\frac{n\pi x}{L}\right)\cos\left(\frac{m\pi y}{L}\right)\exp \left[-(n^2+m^2)\frac{\pi^2}{L^2}t\right] $$

The initial condition for this last piece will have you solve a double Fourier series in $x$, $y$. You can use linearity to simplify calculations a bit, i.e.

\begin{align} u_1(x,y,0) &= T_0 - T_1 - \frac{T_0-T_1}{L}x = \sum_{n,m} c_{1(n,m)} \sin\left(\frac{n\pi x}{L}\right)\cos\left(\frac{m\pi y}{L}\right) \\ u_2(x,y,0) &= -\sum_{n=2k+1} \frac{4aL}{n^2\pi^2\cosh(n\pi)} \sin\left(\frac{n\pi x}{L}\right)\left[\cosh\left(\frac{n\pi y}{L}\right) - \cosh\left(\frac{n\pi(L-y)}{L}\right)\right] \\ &\qquad = \sum_{n,m} c_{2(n,m)} \sin\left(\frac{n\pi x}{L}\right)\cos\left(\frac{m\pi y}{L}\right) \end{align}

The first boundary function is constant in $y$ and the second is already a partial Fourier series, so they simplify to \begin{align} T_0 - T_1 - \frac{T_0-T_1}{L}x &= \sum_{n=1}^\infty c_{1(n,0)} \sin\left(\frac{n\pi x}{L}\right) \\ -\frac{4aL}{n^2\pi^2\sinh(n\pi)} \left[\cosh\left(\frac{n\pi y}{L}\right) - \cosh\left(\frac{n\pi(L-y)}{L}\right)\right] &= c_{2(n,0)} + \sum_{m=1}^\infty c_{2(n,m)}\cos\left(\frac{m\pi y}{L}\right) \end{align}

The last round of integration finishes the job

  • $\begingroup$ Thanks a lot for your help!! I can understand the steady solution process, but the transient results like the last two equations are hard for me. Do I have to solve the last two equations to obtain 'c_1(n,0) and c_2(n,m) '? and does the c(n,m) in u(x,y,t) equals to c1(n,m)+c2(n,m)? $\endgroup$ Jun 25 '18 at 22:28
  • $\begingroup$ Yes, $c_1+c_2=c$ for all indices $n,m$. You'll need to obtain the Fourier series of those LHS functions and then match them up. The integrals are a bit cumbersome, but they are doable. $\endgroup$
    – Dylan
    Jun 26 '18 at 4:45
  • $\begingroup$ You can also obtain the entire $u(x,y,0)$ Fourier series as a double integral, but I find splitting it up this way makes it a bit easier. Since the cosine series of a constant is that same constant, you can conclude that the coefficient of the linear $x$ piece, which I've marked $c_{1(n,m)}$, is zero for $m\ne 0$. $\endgroup$
    – Dylan
    Jun 26 '18 at 4:53

Consider the method of separation of variables. Small googling provides a solution to 2D heat equation. Perhaps you will have to adjust the solution a little so that it works for your exact boundary conditions.


In the end, you will likely arrive at an infinite sum of terms for your solution.

  • $\begingroup$ Really appreciate your answer! I know I may solve it with separation of variables, but the thing is I'm not familiar with it, and I have read some similar cases to figure out the process of separation of variables, but I'm still unable to solve this issues. That is the reason why I post the question. $\endgroup$ Jun 21 '18 at 20:41
  • $\begingroup$ This should be a comment, not an answer $\endgroup$
    – Dylan
    Jun 23 '18 at 8:26

Your Answer

By clicking “Post Your Answer”, you agree to our terms of service, privacy policy and cookie policy

Not the answer you're looking for? Browse other questions tagged or ask your own question.