# Blow-up for the quasilinear heat equation $u_t= u \ u_{x x}$ or the related $w_t= \left(w_x e^w\right)_x$

What kind of approaches can be used to study the following quasilinear parabolic pde for a scalar function $u=u(x,t)$ ? $$u_t= u \ u_{x x}$$

The physical problem where this pde comes from dictates that the Cauchy problem of interest corresponds to an initial condition which is a third-order polynomial with no-constant term $$u(x,0) = x \ (x^2-s x + p), \ \ \ x \in \mathbb{R}^+, \ \ \ s p \neq 0$$ Note that $u(x,0)$ is initially $\propto x$, and that $u(x,t)$ will remain so during the evolution as long as the pde makes sense. Indeed, and this is the crux of the matter, numerical experiments and heuristic arguments strongly indicate that there exists a blow-up time $0< T=T(s,p) < \infty$, where the solution explodes to infinity, with $$u(x,t) \sim \frac{f(x)}{T-t}, \ \ \ t \to T^-$$

• Can one prove this, and relate $T$ and the amplitude $f(x)$ to the parameters $s,p$ entering the initial condition ?
• It is possible (and how !) to continue the solution after the blow-up time ?
• A colleague of mine suggested that I could exchange the role of the dependent and independent variables, and/or to look at a paper by Clarkson, Fokas, \& Ablowitz, "Hodograph transformations of linearizable partial differential equations", where they study similar equations and relate them to the Harry-Dym equation and "solvable" (by inverse scattering transform or Painlevé reductions) variants. Needless to say, the functional form of the terms in my pde do not fall under the conditions of applicability of Clarkson-fokas-Ablowitz' main theorem, while the hodograph transformation: $$(t,x,u(x,t)) \longrightarrow (t,y=u,x=v(y,t))$$ so that (using $u_x = 1/v_y$, $u_{x x}=-v_y^{-3}v_{y y}$, $u_t= -v_t/v_y$) $$v_t + y (v_y^{-1})_y =0$$ does not seem to give a pde much easier to solve (or to recognize, such as the porous medium equation)

Edit: an alternative perhaps preferable form (closer to the porous medium equation!) is for the function $w=\frac{1}{2}{\rm Log}(u^2)$, which obeys $$w_t = \left(w_x e^w\right)_x$$ The Lie-point symmetries of the latter are recalled in this review (section 5), but none would match the initial condition I have to deal with.

-
I always get nervous when my solutions grow at spatial infinity. One possibility would be to make the transformation v = (x+1)^{-4}u so that v is, at least initially, well behaved at spatial infinity. The resulting equation has a v^2 term on the right hand side. This gives bounds for finite-time blow up at least for initial data with negative minima. For positive initial data, this form of the equation at least highlights the competition between the curvature and amplitude of a positive maximum. –  Aaron Hoffman Feb 6 '12 at 13:57
Thanks Aaron for the comment. Do you think that such a method would also entail a better-behaved numerical scheme ? –  Ivan Dornic Feb 6 '12 at 14:57
The solution (x^2/2+ax+b)/(T-t) blows up at T, but has the wrong initial value. –  Bob Terrell Feb 7 '12 at 21:49
@Bob. Yes, I had seen this, and it also turns out that the separable solution you mention comes out through an ODE reduction after studying the Lie-point symmetries for the related function $w=\frac{1}{2}{\rm Log}(u^2)$ obeying the perhaps somewhat better-looking $w_t=(w_x e^w)_x$. I have edited my question accordingly, and I have given a ref. concerning the Lie symmetry-analysis. –  Ivan Dornic Feb 9 '12 at 11:27
There is a sort of Backlund transformation. If v_x = log(uv) and v_t = U_x+v_x, then u_t = u u_{xx} iff v_t = v v_{xx}. For example if u=1, then v is a traveling wave with profile satisfying v_x=log v. But I don't see any better examples so far. –  Bob Terrell Feb 9 '12 at 20:55