Adaptive High-Order Finite-Difference Method for Nonlinear Wave Problems

We discuss a scheme for the numerical solution of one-dimensional initial value problems exhibiting strongly localized solutions or finite-time singularities. To accurately and efficiently model such phenomena we present a full space-time adaptive scheme, based on a variable order spatial finite-difference scheme and a 4th order temporal integration with adaptively chosen time step. A wavelet analysis is utilized at regular intervals to adaptively select the order and the grid in accordance with the local behavior of the solution. Through several examples, taken from gasdynamics and nonlinear optics, we illustrate the performance of the scheme, the use of which results in several orders of magnitude reduction in the required degrees of freedom to solve a problem to a particular fidelity.