Computing small solutions of Burgers' equation backwards in time

AbstractWe construct and analyze an algorithm for the numerical computation of Burgers' equation for preceding times, given an a priori bound for the solution and an approximation to the terminal data. The method is based on the “backward beam equation” coupled with an iterative procedure for the solution of the nonlinear problem via a sequence of linear problems. We also present the results of several numerical experiments. It turns out that the procedure converges “asymptotically,” i.e., in the same manner in which an asymptotic expansion converges. This phenomenon seems related to the “destruction of information,” at t = 0, which is typical in backwards dissipative equations. We derive a priori stability estimates for the analytic backwards problem, and we observe that in many numerical experiments, the distance backwards in time where significant accuracy can be attained is much larger than would be expected on the basis of such estimates. The method is useful for small solutions. Problems where steep gradients occur require considerably more precision in measurement. The algorithm is applicable to other semilinear problems