An adaptive local mesh refinement method for time-dependent partial differential equations

Abstract We discuss an adaptive local mesh refinement procedure for solving time-dependent initial boundary value problems for vector systems of partial differential equations on rectangular spatial domains. The method identifies and groups regions having large local error indicators into rectangular clusters. The time step and computational cells within clustered rectangles are recursively divided until a prescribed tolerance is satisfied. The refined meshes are properly nested within coarser mesh boundaries; thus, simplifying the prescription of interface conditions at boundaries between fine and coarse meshes. Our method may be used with several numerical solution strategies and error indicators. The meshes may be nonuniform and either stationary or moving. A code based on our refinement procedure is used with a MacCormack finite difference method to solve some examples involving systems of hyperbolic conservation laws.

[1]  I. I. Glass,et al.  Domains and boundaries of non-stationary oblique shock-wave reflexions. 1. Diatomic gas , 1979, Journal of Fluid Mechanics.

[2]  J. Flaherty,et al.  A two-dimensional mesh moving technique for time-dependent partial differential equations , 1986 .

[3]  Ivo Babuška,et al.  Adaptive Methods and Error Estimation for Elliptic Problems of Structural Mechanics. , 1983 .

[4]  William Gropp,et al.  A Test of Moving Mesh Refinement for 2-D Scalar Hyperbolic Problems , 1980 .

[5]  S. F. Davis TVD finite difference schemes and artificial viscosity , 1984 .

[6]  Marsha Berger Data structures for adaptive mesh refinement , 1983 .

[7]  Joseph E. Flaherty,et al.  A moving finite element method with error estimation and refinement for one-dimensional time dependent partial differential equations , 1986 .

[8]  J. Flaherty,et al.  A Local Refinement Finite-Element Method for Two-Dimensional Parabolic Systems , 1988 .

[9]  Alfred V. Aho,et al.  The Design and Analysis of Computer Algorithms , 1974 .

[10]  Joseph E. Flaherty,et al.  A moving-mesh finite element method with local refinement for parabolic partial differential equations , 1986 .

[11]  I. Babuška,et al.  The finite element method for parabolic equations , 1982 .

[12]  P. Woodward,et al.  The numerical simulation of two-dimensional fluid flow with strong shocks , 1984 .

[13]  I. I. Glass,et al.  Nonstationary Oblique Shock-Wave Reflections: Actual Isopycnics and Numerical Experiments , 1978 .

[14]  M. Berger ON CONSERVATION AT GRID INTERFACES. , 1987 .

[15]  Ivo Babuška,et al.  The finite element method for parabolic equations , 1982 .

[16]  R. Maccormack The Effect of Viscosity in Hypervelocity Impact Cratering , 1969 .

[17]  Ellis Horowitz,et al.  Fundamentals of Data Structures , 1984 .

[18]  M. Berger,et al.  Adaptive mesh refinement for hyperbolic partial differential equations , 1982 .

[19]  J. E. Flaherty,et al.  A Local Refinement Finite Element Method for Time Dependent Partial Differential Equations , 1985 .

[20]  William R. Goodin,et al.  Numerical solution of the atmospheric diffusion equation for chemically reacting flows , 1982 .

[21]  Werner C Rheinboldt,et al.  Feedback Systems and Adaptivity for Numerical Computations. , 1983 .

[22]  Richard G. Hindman,et al.  Generalized Coordinate Forms of Governing Fluid Equations and Associated Geometrically Induced Errors , 1982 .

[23]  Joe F. Thompson A survey of dynamically-adaptive grids in the numerical solution of partial differential equations , 1984 .