6.2. Separation of variables: Misc equations

$\renewcommand{\Re}{\operatorname{Re}}$ $\renewcommand{\Im}{\operatorname{Im}}$ $\newcommand{\erf}{\operatorname{erf}}$ $\newcommand{\dag}{\dagger}$ $\newcommand{\const}{\mathrm{const}}$ $\newcommand{\arcsinh}{\operatorname{arcsinh}}$

6.2. Separation of variables: Misc equations

  1. Schrödinger equation
  2. 1D wave equation
  3. Laplace equation in half-strip
  4. Laplace equation in rectangle

In the previous Section 6.1 we considered heat equation.

Schrödinger equation

Consider problem \begin{align} & u_t= iku_{xx},&& t>0,\ 0< x< l,\label{eq-6.2.1}\\[3pt] & (\alpha_0u_x-\alpha u)|_{x=0}=(\beta_0u_x+\beta u)|_{x=l}=0. \label{eq-6.2.2}\\[3pt] & u|_{t=0}=g(x)\label{eq-6.2.3} \end{align} where either $\alpha_0=0$, $\alpha=1$ and we have a Dirichlet boundary condition or $\alpha_0=1$ and we have either Neumann or Robin boundary condition and the same at $x=l$.

Let us consider a simple solution $u(x,t)=X(x)T(t)$; then separating variables we arrive to $\frac{T'}{T}=i k\frac{X''}{X}$ which implies $X''+\lambda X=0$, \begin{equation} T'=-i k\lambda T \label{eq-6.2.4} \end{equation} (explain, how). We also get boundary conditions $(\alpha_0 X'-\alpha X)=(\beta_0 X'+\beta X)(l)=0$ (explain, how).

So, we have eigenvalues $\lambda_n$ and eigenfunctions $X_n$ ($n=1,2,\ldots$) and equation (\ref{eq-6.2.4}) for $T$, which results in \begin{equation} T_n=A_ne^{-ik\lambda_n t} \label{eq-6.2.5} \end{equation} and therefore a simple solution is \begin{equation} u_n=A_ne^{-ik\lambda_n t}X_n(x) \label{eq-6.2.6} \end{equation} and we look for a general solution in the form \begin{equation} u=\sum_{n=1}^\infty A_ne^{-ik\lambda_n t}X_n(x). \label{eq-6.2.7} \end{equation} Again, taking in account initial condition (\ref{eq-6.2.3}) we see that \begin{equation} u=\sum_{n=1}^\infty A_nX_n(x). \label{eq-6.2.8} \end{equation} and therefore \begin{equation} A_n=\frac{1}{\|X_n\|^2}\int_0^l g(x)X_n(x)\,dx. \label{eq-6.2.9} \end{equation}

Remark 1.

  1. Formula (\ref{eq-6.2.6}) shows that the problem is well-posed for $t<0$ and $t>0$.
  2. However there is no stabilization.
  3. What we got is a finite interval version of analysis of Subsection 5.3.2.

1D wave equation

Consider problem \begin{align} & u_{tt}= c^2u_{xx},&& 0< x< l,\label{eq-6.2.10}\\[3pt] & (\alpha_0u_x-\alpha u)|_{x=0}=(\beta_0u_x+\beta u)|_{x=l}=0. \label{eq-6.2.11}\\[3pt] & u|_{t=0}=g(x),\qquad u_t|_{t=0}=h(x).\label{eq-6.2.12} \end{align} Actually we started from this equation in Section 4.1 but now we consider more general boundary conditions. Now we have \begin{equation} T''= -k\lambda T \label{eq-6.2.13} \end{equation} and we have \begin{equation} T_n=\left\{\begin{aligned} & A_n \cos (\omega_n t)+ B_n \sin (\omega_n t) && \omega_n = c\lambda_n^{\frac{1}{2}}\text{as } \lambda_n>0,\qquad\\ & A_n + B_n t && \text{as } \lambda_n=0,\\ & A_n \cosh (\eta_n t)+ B_n \sinh (\eta_n t) && \eta_n = c(-\lambda_n)^{\frac{1}{2}} \text{as } \lambda_n>0, \end{aligned}\right. \label{eq-6.2.14} \end{equation} and respectively we get \begin{equation} u =\sum _n T_n(t)X_n(x) \label{eq-6.2.15} \end{equation} and we find from initial conditions \begin{align} & A_n=\frac{1}{\|X_n\|^2}\int_0^l g(x)X_n(x)\,dx, \label{eq-6.2.16}\\ & B_n = \frac{1}{\|X_n\|^2}\int_0^l h(x)X_n(x)\,dx \times\left\{\begin{aligned} & \frac{1}{\omega_n} && \text{as } \lambda_n>0,\\ & 1 && \text{as } \lambda_n=0,\\ & \frac{1}{\eta_n} && \text{as } \lambda_n<0. \label{eq-6.2.17} \end{aligned}\right. \end{align}

Laplace equation in half-strip

Consider problem \begin{align} & \Delta u:=u_{xx}+u_{yy}=0,&& y>0,\ 0< x< l,\label{eq-6.2.18}\\[3pt] & (\alpha_0u_x-\alpha u)|_{x=0}=(\beta_0u_x+\beta u)|_{x=l}=0, \label{eq-6.2.19}\\[3pt] & u|_{y=0}=g(x).\label{eq-6.2.20} \end{align} To make it uniquely solvable we need to add condition $|u|\le M$.

Again separating variables $u(x,y)=X(x)Y(y)$ we get \begin{equation} Y''= \lambda Y \label{eq-6.2.21} \end{equation} and therefore assuming that $\lambda>0$ we get \begin{equation} Y=Ae^{-\sqrt{\lambda}y}+Be^{\sqrt{\lambda}y} \label{eq-6.2.22} \end{equation} We discard the second term in the right-hand expression of (\ref{eq-6.2.22}) because it is unbounded. However if we had Cauchy problem (i.e. $u|_{y=0}=g(x)$, $u_y|_{y=0}=h(x)$) we would not be able to do this and this problem will be ill-posed.

So, $u=Ae^{-\sqrt{\lambda}y}$ and (\ref{eq-6.2.20}) and \begin{equation} u_n=A_ne^{-\sqrt{\lambda_n}y}X_n(x) \label{eq-6.2.23} \end{equation} and assuming that all $\lambda_n>0$ we get \begin{equation} u=\sum_n A_n e^{-\sqrt{\lambda_n}y}X_n(x) \label{eq-6.2.24} \end{equation} and (\ref{eq-6.2.19}) yields \begin{equation} A_n=\frac{1}{\|X_n\|^2}\int_0^l g(x)X_n(x)\,dx. \label{eq-6.2.25} \end{equation}

Remark 2.

  1. If there is eigenvalue $\lambda=0$ we have $Y= A+By$ and as we are looking for a bounded solution we discard the second term again; so our analysis remain valid as all $\lambda_n\ge 0$.
  2. If we have a problem \begin{align} & \Delta u:=u_{xx}+u_{yy}=0,&& y>0,\ 0< x< l,\label{eq-6.2.26}\\[3pt] & (\alpha_0u_x-\alpha u)|_{x=0}=\phi(y), \label{eq-6.2.27}\\[3pt] & (\beta_0u_x+\beta u)|_{x=l}=\psi(y), \label{eq-6.2.28}\\[3pt] & u|_{y=0}=g(x)\label{eq-6.2.29} \end{align} with $g(x)=0$ we could reduce it by the method of continuation to the problem in the whole strip and solve it by Fourier transform (see Subsecion 5.3.5).

  3. In the general case we can find $u= u_{(1)} +u_{(2)}$ where $u_{(1)}$ solves problem with $g=0$ and $u_{(2)}$ solves problem with $\phi=\psi=0$ (explain how it follows from the linearity).

  4. One can replace Dirichlet boundary condition $u|_{y=0}$ by Robin boundary condition $(u_y-\gamma u)|_{y=0}=g(x)$ ($\gamma\ge 0$) but there is an exceptional case: there is an eigenvalue $\lambda_0=0$ and as $y=0$ we have Neumann boundary condition.

  5. In this exceptional case (usually as we have Neumann b.c. everywhere--as $x=0$, $x=l$, $y=0$) a required solution simply does not exists unless $\int_0^l g(x)X_0(x)\,dx=0$.

Laplace equation in rectangle

Consider problem \begin{align} &\Delta u:=u_{xx}+u_{yy}=0,&& 0< y< b,\ 0< x< a,\label{eq-6.2.30}\\[3pt] &(\alpha_0u_x-\alpha u)|_{x=0}=(\beta_0u_x+\beta u)|_{x=a}=0, \label{eq-6.2.31}\\[3pt] & u|_{y=0}=g(x), \label{eq-6.2.32}\\[3pt] & u|_{y=b}=h(x).\label{eq-6.2.33} \end{align}

Then we get (\ref{eq-6.2.21}) and (\ref{eq-6.2.22}) again but with two b.c. we cannot diacard anything; we get instead \begin{align} & A_n && + B_n&& =g_n,\label{eq-6.2.34}\\[3pt] & A_n e^{-\sqrt{\lambda_n}b} && + B_n e^{\sqrt{\lambda_n}b} && =h_n\label{eq-6.2.35} \end{align} where $g_n$ and $h_n$ are Fourier coefficients of $g$ and $h$ respectively, which implies \begin{align*} & A_n= \frac{e^{\sqrt{\lambda_n}b}}{2\sinh (\sqrt{\lambda_n}b)}g_n -\frac{1}{2\sinh (\sqrt{\lambda_n}b)}h_n) , \\[3pt] & B_n= -\frac{e^{-\sqrt{\lambda_n}b}}{2\sinh (\sqrt{\lambda_n}b)}g_n +\frac{1}{2\sinh (\sqrt{\lambda_n}b)}h_n \end{align*} and therefore \begin{equation} Y_n(y)=\frac{\sinh (\sqrt{\lambda_n}(b-y))}{\sinh (\sqrt{\lambda_n}b)} g_n+ \frac{\sinh (\sqrt{\lambda_n}y)}{\sinh (\sqrt{\lambda_n}b)} h_n. \label{eq-6.2.36} \end{equation} One can see easily that $\frac{\sinh (\sqrt{\lambda_n}(b-y))}{\sinh (\sqrt{\lambda_n}b)}$ and $\frac{\sinh (\sqrt{\lambda_n}y)}{\sinh (\sqrt{\lambda_n}b)}$ are bounded as $0\le y\le b$.

Problem 1. Investigate other boundary conditions (Robin, Neumann, mixed).

Remark 3.

  1. There is an exeptional case: there is an eigenvalue $\lambda_0=0$ and as $y=0$ and $y=b$ we have Neumann boundary conditions. Then solution does not exist unless $\int_0^a h(x)X_0(x)\,dx -\int_0^a g(x)X_0(x)\,dx=0$.

  2. We can consider general b.c. with $(\alpha_0u_x-\alpha u)|_{x=0}=\phi(y)$, $(\beta_0u_x+\beta u)|_{x=a}=\psi(y)$. Then we can find $u= u_{(1)} +u_{(2)}$ where $u_{(1)}$ solves problem with $g=h=0$ and $u_{(2)}$ solves problem with $\phi=\psi=0$ (explain how it follows from the linearity). The second problem is also "our" problem with $x$ and $y$ permutted.

  3. Assume that we have Neumann b.c. everywhere--as $x=0$, $x=a$, $y=0$, $y=b$. Then solution does not exist unless \begin{equation} \int_0^a h(x)\,dx -\int_0^a g(x)\,dx + \int_0^b \psi(y)\,dy - \int_0^b \phi(y)\,dy=0 \label{eq-6.2.37} \end{equation} which means that the total heat flow is $0$. How from two assumptions we can get one? Well, we just need to consider $\phi=g=0$, $\psi=\frac{y}{b}$, $h=-\frac{x}{a}$ (explain why) but there is a solution $u=\frac{y}{b}-\frac{x}{a}$ for that.

$\Leftarrow$  $\Uparrow$  $\Rightarrow$