The Frank-Kamenetskii or "spontaneous combustion" equation is the PDE

$$ {\partial u\over \partial t} = {\partial^2u\over \partial x^2} + A\exp(u). $$

On the interval $[-1,1]$ with zero initial and boundary conditions, solutions to this equation blow up to infinity in finite time if $A$ is bigger than about $0.878$. For smaller $A$, solutions converge to a steady state.

Here we compute some of these steady-state solutions, which are solutions of the ODE boundary value problem

$$ u''+A\exp(u)=0,\qquad u(-1)=u(1)=0. $$

N = chebop([-1 1]);
N.bc = 'dirichlet';
FS = 'fontsize';
for A = [.2 .4 .6 .8 .87]
  N.op = @(u) diff(u,2) + A*exp(u);
  u = N\0;
  plot(u,'color',[.6 0 .5],'linewidth',2), grid on, hold on
  text(-.1,max(u)+.04,['A = ' num2str(A)],FS,14)
end
axis([-1 1 0 1.2])
title('Frank-Kamenetskii blowup equation',FS,14)

References

  1. H. Fujita, On the nonlinear equations $\Delta u + \exp(u) = 0$ and $dv/dt = \Delta v + \exp(v)$, Bulletin of the American Mathematical Society, 75 (1969), 132-135.