STABILITY AND ACCURACY OF DIFFERENCING METHODS FOR VISCOPLASTIC MODELS IN WAVECODES

The numerical stability and truncation error of a family of differencing schemes for viscoplastic constitutive relations in wavecodes is investigated. A von Neumann stability analysis is performed for a one-dimensional model problem. This analysis identifies two differencing methods that have no restriction on the time step size beyond the usual Courant-Friedrichs-Lewy condition. One of these methods is first-order accurate, and the other is second-order accurate. Implementation of one of these methods in the three-dimensional wavecode CTH is discussed. © 1993 Academic Press. All rights reserved.