0

I'm interested in testing some algorithms on the heat equation, and I'd like to assess their accuracy. When evolving a Hamiltonian system, one has the energy to check the validity/correctness of the algorithm, which turns out to be (arguably) the most important quantity to check. Is there a similar quantity for the heat equation - i.e. a quantity which would very well characterize the error of the algorithm?

I could compute the time-evolution of the "energy" of the heat equation easily by the following

$$u_t(x,t)=\hat{H}u(x,t),~~~\hat{H}=\frac{\partial^2}{\partial x^2}$$

The "energy" would then be the expectation value of the $\hat{H}$ operator, which is simple to evaluate in Fourier space. However, I don't know if this "quantity" would even be useful.

Arturo don Juan
  • 489
  • 2
  • 8
  • The standard approach is the method of manufactured solutions; in fact, it's best to use that method even when there is a conserved quantity. – Kirill Sep 15 '16 at 04:40
  • 1
    I don't think that quantity is conserved, it is my understanding that it should be a non-increasing quantity. – nicoguaro Sep 15 '16 at 14:17

0 Answers0