A fixed grid method for capturing the motion of self-intersecting wavefronts and related PDEs

Moving surfaces that self-intersect arise naturally in the geometric optics model of wavefront motion. Standard ray tracing techniques can be used to compute these motions, but they lose resolution as rays diverge. In this pa per we develop numerical methods that maintain uniform spatial resolution o f the front at all times. Our approach is a fixed grid, wavefront capturing formulation based on the Dynamic Surface Extension method of Steinhoff and Fan (Technical report, University of Tennessee Space Institute). The new m ethods can treat arbitrarily complicated self intersecting fronts, as well as refraction, reflection, and focusing. We also develop methods fur curvat ure-dependent front motions and the motion of filaments. We validate our me thods with numerical experiments. (C) 2000 Academic Press.