Adaptive Method with Mesh Moving and Local Mesh Refinement for Time-Dependent Partial Differential Equations

摘要

The authors discuss mesh moving, static mesh regeneration, and local mesh refinement algorithms that can be used with a finite difference or finite element scheme to solve initial-boundary value problems for vector systems of time-dependent partial differential equations in two-space dimensions and time. A coarse base mesh of quadrilateral cells is moved by an algebraic mesh movement function so that it may follow and isolate spatially distinct phenomena. The local mesh refinement method recursively divides the time step and spatial cells of the moving base mesh in regions were error indicators are high until a prescribed tolerance is satisfied. The static mesh regeneration procedure is used to create a new base mesh when the existing ones become too distorted. In order to test our adaptive algorithms, the authors implemented them in a system code with an initial mesh generator, a MacCormack finite difference scheme for hyperbolic systems, and an error indicator based upon estimates of the local discretization error obtained by Richardson extrapolation. Results are presented for several computational examples. (kr)