DeepAI AI Chat
Log In Sign Up

Non-Conforming Mesh Refinement for High-Order Finite Elements

by   Jakub Červený, et al.
Lawrence Livermore National Laboratory

We propose a general algorithm for non-conforming adaptive mesh refinement (AMR) of unstructured meshes in high-order finite element codes. Our focus is on h-refinement with a fixed polynomial order. The algorithm handles triangular, quadrilateral, hexahedral and prismatic meshes of arbitrarily high order curvature, for any order finite element space in the de Rham sequence. We present a flexible data structure for meshes with hanging nodes and a general procedure to construct the conforming interpolation operator, both in serial and in parallel. The algorithm and data structure allow anisotropic refinement of tensor product elements in 2D and 3D, and support unlimited refinement ratios of adjacent elements. We report numerical experiments verifying the correctness of the algorithms, and perform a parallel scaling study to show that we can adapt meshes containing billions of elements and run efficiently on 393,000 parallel tasks. Finally, we illustrate the integration of dynamic AMR into a high-order Lagrangian hydrodynamics solver.


page 19

page 20

page 23

page 24


A simple technique for unstructured mesh generation via adaptive finite elements

This work describes a concise algorithm for the generation of triangular...

High order transition elements: The xNy-element concept – Part I: Statics

Advanced transition elements are of utmost importance in many applicatio...

Tailored meshing for parallel 3D electromagnetic modeling using high-order edge elements

We present numerical experiments for geophysics electromagnetic (EM) mod...

Fully Parallel Mesh I/O using PETSc DMPlex with an Application to Waveform Modeling

Large-scale PDE simulations using high-order finite-element methods on u...

Automatic penalty and degree continuation for parallel pre-conditioned mesh curving on virtual geometry

We present a distributed parallel mesh curving method for virtual geomet...

An Optimized, Parallel Computation of the Ghost Layer for Adaptive Hybrid Forest Meshes

We discuss parallel algorithms to gather topological information about o...

Automatic, high-order, and adaptive algorithms for Brillouin zone integration

We present efficient methods for Brillouin zone integration with a non-z...