Mathematics Stack Exchange is a question and answer site for people studying math at any level and professionals in related fields. Join them; it only takes a minute:

Sign up
Here's how it works:
  1. Anybody can ask a question
  2. Anybody can answer
  3. The best answers are voted up and rise to the top

I would like to know if the Rittner equation : $$\partial_{t}{\varPhi(x,t)=k\partial_{xx}{\varPhi(x,t)}}-\alpha{\partial_{x}{\varPhi(x,t)}-\beta{\varPhi(x,t)}+g(x,t)}$$

can be solved using the Lax pair method or the Fokas method.


share|cite|improve this question
up vote 3 down vote accepted

With $\alpha$ and $\beta$ being constants, this equation can be reduced to a heat equation and solved exactly. The first step is to remove the $\beta$ term. We do this by setting

$$\Phi(x,t)=e^{-\beta t}\Phi_1(x,t).$$

The equation becomes

$$\partial_t\Phi_1(x,t)=k\partial_{xx}\Phi_1(x,t)-\alpha\partial_x\Phi_1(x,t)+g(x,t)e^{\beta t}.$$

The next step is to remove the drift term. This can be done by setting


with $a$ a constant to be determined. We just note that $\partial_x(e^{ax}\Phi_2(x,t))=e^{ax}(\partial_x+a)\Phi_2(x,t)$ and so one gets

$$\partial_t\Phi_2(x,t)+c\Phi_2(x,t)=k(\partial_x+a)^2\Phi_2(x,t)-\beta(\partial_x+a)\Phi_2(x,t)+g(x,t)e^{\beta t}e^{-ax}$$

and so, choosing $a=\frac{\beta}{2k}$ and $c=-\frac{\beta^2}{4k}$ the equation becomes a heat equation

$$\partial_t\Phi_2(x,t)=k\partial_{xx}\Phi_2(x,t)+g(x,t)e^{\left(\beta-\frac{\beta^2}{4k}\right) t}e^{-\frac{\beta}{2k}x}$$

This equation can be solved exactly by the kernel

$$\Delta(x,t)=\frac{e^{-\frac{x^2}{4kt}}}{\sqrt{4\pi kt}}$$

Then, a full solution is given by

$$\Phi_2(x,t)=\tilde\Phi_2(t,x)+\int dx'dt'\Delta(x-x',t-t')g(x',t')e^{\left(\beta-\frac{\beta^2}{4k}\right) t'}e^{-\frac{\beta}{2k}x'}$$

being $\tilde\Phi_2(t,x)$ a solution of the homogeneous equation. Then, it easy to recognize that kernel for the initial equation is just

$$\Delta_0(x,t)=\Delta(x,t)e^{-\left(\beta-\frac{\beta^2}{4k}\right) t}e^{\frac{\beta}{2k}x}$$

share|cite|improve this answer

Setting $\psi(x,t) = e^{(\alpha^2/4k + \beta)t - \alpha x/2k} \Phi(x,t)$ and plugging in I get $$ \begin{align*} \psi_t(x,t) - k \psi_{xx}(x,t) &= e^{(\alpha^2/4k + \beta)t - \alpha x/2k}g(x,t) \\ &\equiv f(x,t). \end{align*} $$ If the forcing $f(x,t)$ is in $L^1(\mathbf{R}\times \mathbf{R})$ then we can use the heat kernel $$ K(x,t)=\begin{cases} (4\pi k t)^{-1/2} e^{-x^2/2kt},& t>0 \\ 0,& t\leq 0 \end{cases} $$ and a solution to the problem on the domain $\mathbf{R}\times\mathbf{R}$ is $\psi= K*f$. Unfortunately this puts very harsh restrictions on the original forcing term $g=g(x,t)$ because of the exponential growth of the factor $e^{(\alpha^2/4k + \beta)t - \alpha x/2k}$.

Similarly, if one is interested in the Cauchy problem with $\Phi(x,0)=\Phi_0(x)$, then the new problem has initial data $\psi_0(x) = e^{-\alpha x/2k} \Phi_0(x)$. Using heat kernel methods this can again lead to severe restrictions. Recall the solution to the homogeneous Cauchy problem would be $$ \psi(x,t)=K_t * \psi_0(x) $$ which is well defined in a classical sense if, for instance, $\psi_0 \in L^p(\mathbf{R})$ for $1\leq p\leq \infty$. This means the associated initial data $\Phi_0(x)$ would need exponential decay, which is a little artificial.

For the Cauchy problem for the Rittner equation on the full line (for simplicity consider the homogeneous problem $g=0$, the more general case is similar) one can use the Fourier transform. Setting $\psi = e^{\beta t} \Phi$ we have the Cauchy problem $$ \begin{align*}\psi_t - k \psi_{xx} + \alpha \psi_x &=0, \\ \psi(x,0) &= \Phi_0(x), \end{align*} $$ the solution to which is $$ \psi(x,t) = \frac{1}{2\pi} \int_{-\infty}^\infty e^{\mathrm{i} \lambda x -\lambda (k\lambda + \mathrm{i}\alpha) t} \hat{\Phi}_0(\lambda)\, \mathrm{d}\lambda. $$ The problem on the half-line is more subtle. The Fokas method provides a solution in this case (for given boundary data at $x=0$, $\Phi(0,t)=\Phi_1(t)$ say) in terms of an integral along a smooth contour in the complex $\lambda$-plane. Now if the initial data $\Phi_0(x)$ has sufficient exponential decay (which is a big restriction), then one can show that this contour can be deformed onto the real $\lambda$-axis, giving a Fourier type integral as in the case of the Cauchy problem on the full line. However, in general, this contour cannot be shifted.

This problem is treated in Fokas' book, A unified approach to boundary value problems. He emphasizes that this inability to shift the contour back to the real axis (in the half-line problem) reflects the fact there is no classical transform for the problem.

share|cite|improve this answer

Your Answer


By posting your answer, you agree to the privacy policy and terms of service.

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