Cole–Hopf transformation

Template:Short description The Cole–Hopf transformation is a change of variables that allows to transform a special kind of parabolic partial differential equations (PDEs) with a quadratic nonlinearity into a linear heat equation. In particular, it provides an explicit formula for fairly general solutions of the PDE in terms of the initial datum and the heat kernel.

Consider the following PDE:

$$u_t-a\mathrm{\Delta }u+b\Vert \mathrm{\nabla }u{\Vert }^2=0,u(0,x)=g(x)$$

where $x\in {ℝ}^n$, $a,b$ are constants, $\mathrm{\Delta }$ is the Laplace operator, $\mathrm{\nabla }$ is the gradient, and $\Vert \cdot \Vert$ is the Euclidean norm in ${ℝ}^n$. By assuming that $w=\varphi (u)$, where $\varphi (\cdot )$ is an unknown smooth function, we may calculate:

$$w_t={\varphi }^{\prime }(u)u_t,\mathrm{\Delta }w={\varphi }^{\prime }(u)\mathrm{\Delta }u+{\varphi }^{″}(u)\Vert \mathrm{\nabla }u{\Vert }^2$$

Which implies that:

$$\begin{array}{cc}{w}_t={\varphi }^{\prime }(u)u_t& ={\varphi }^{\prime }(u)(a\mathrm{\Delta }u-b\Vert \mathrm{\nabla }u{\Vert }^2)\\ & =a\mathrm{\Delta }w-(a{\varphi }^{″}+b{\varphi }^{\prime })\Vert \mathrm{\nabla }u{\Vert }^2\\ & =a\mathrm{\Delta }w\end{array}$$

if we constrain $\varphi$ to satisfy $a{\varphi }^{″}+b{\varphi }^{\prime }=0$. Then we may transform the original nonlinear PDE into the canonical heat equation by using the transformation:

Template:Equation box 1This is the Cole-Hopf transformation. With the transformation, the following initial-value problem can now be solved:

$$w_t-a\mathrm{\Delta }w=0,w(0,x)=e^{-bg(x)/a}$$

The unique, bounded solution of this system is:

$$w(t,x)=\frac{1}{(4\pi at{)}^{n/2}}{\int }_{{ℝ}^n}{e}^{-\Vert x-y{\Vert }^2/4at-bg(y)/a}dy$$

Since the Cole–Hopf transformation implies that $u=-(a/b)\log w$, the solution of the original nonlinear PDE is:

$$u(t,x)=-\frac{a}{b}\log \left[\frac{1}{(4\pi at{)}^{n/2}}{\int }_{{ℝ}^n}{e}^{-\Vert x-y{\Vert }^2/4at-bg(y)/a}dy\right]$$

The complex form of the Cole–Hopf transformation can be used to transform the Schrödinger equation to the Madelung equation.

Applications

• Aerodynamics
• Stochastic optimal control
• Solving the viscous Burgers' equation
• Madelung equation