STABILITY OF EXPLICIT SCHEMES IN THE PHYSICAL AND FREQUENCY DOMAINS

We look at some of the issues involved in designing stable explicit nu merical schemes for linear advection equations from two perspectives: (a) in the physical domain, where each scheme represents a particular interpolation of discrete data, and (b) in the frequency domain, where the behavior of each scheme is determined by the spectral characteris tics of the operator that is acting on discrete data. We show that (1) the fully discrete form is equivalent to choosing a value for the dep endent variable from an interpolation of the data in the spatial domai n at the previous time level, (2) interpolation generates a continuous function (polynomial) in the physical space, (3) size of the time ste p used in updating the solution determines the location from where the interpolated value is obtained, and (4) if a choice of step size show s amplification in the spectral domain, interpolation in the physical domain exceeds the bounds set by the discrete data at a spatial locati on corresponding to this step size. Comparisons are made between the b ehavior of the operator in the frequency and physical domains; and the amplification in the frequency domain matches the value of extrema ge nerated by the interpolation. Examples to illustrate both perspectives include first and second difference operators, spatial averaging, and various central and upwind schemes for the linear advection equation. (C) 1996 Academic Press, Inc.