PARALLEL STRUCTURES AND DYNAMIC LOAD BALANCING FOR ADAPTIVE FINITE-ELEMENT COMPUTATION

An adaptive technique for a partial differential system automatically adjusts a computational mesh or varies the order of a numerical proced ure to obtain a solution satisfying prescribed accuracy criteria in an optimal fashion. We describe data structures for distributed storage of finite element mesh data as well as software for mesh adaptation, l oad balancing, and solving compressible flow problems. Processor load imbalances are introduced at adaptive enrichment steps during the cour se of a parallel computation. To correct this, we have developed three dynamic load balancing procedures based, respectively, on load imbala nce trees, moment of inertia, and octree traversal. Computational resu lts on an IBM SP2 computer are presented for steady and transient solu tions of the three-dimensional Euler equations of compressible flow. ( C) 1998 Elsevier Science B.V.