Timo Heister

The deal.II Library, Version 8.4

Abstract This paper provides an overview of the new features of the finite element library deal.II version 8.5.

Algorithms and data structures for massively parallel generic adaptive finite element codes

Todays largest supercomputers have 100,000s of processor cores and offer the potential to solve partial differential equations discretized by billions of unknowns. However, the complexity of scaling to such large machines and problem sizes has so far prevented the emergence of generic software libraries that support such computations, although these would lower the threshold of entry and enable many more applications to benefit from large-scale computing. We are concerned with providing this functionality for mesh-adaptive finite element computations. We assume the existence of an “oracle” that implements the generation and modification of an adaptive mesh distributed across many processors, and that responds to queries about its structure. Based on querying the oracle, we develop scalable algorithms and data structures for generic finite element methods. Specifically, we consider the parallel distribution of mesh data, global enumeration of degrees of freedom, constraints, and postprocessing. Our algorithms remove the bottlenecks that typically limit large-scale adaptive finite element analyses. We demonstrate scalability of complete finite element workflows on up to 16,384 processors. An implementation of the proposed algorithms, based on the open source software p4est as mesh oracle, is provided under an open source license through the widely used deal.II finite element software library.

BurnMan: A lower mantle mineral physics toolkit

We present BurnMan, an open-source mineral physics toolbox to determine elastic properties for specified compositions in the lower mantle by solving an Equation of State (EoS). The toolbox, written in Python, can be used to evaluate seismic velocities of new mineral physics data or geodynamic models, and as the forward model in inversions for mantle composition. The user can define the composition from a list of minerals provided for the lower mantle or easily include their own. BurnMan provides choices in methodology, both for the EoS and for the multiphase averaging scheme. The results can be visually or quantitatively compared to observed seismic models. Example user scripts show how to go through these steps. This paper includes several examples realized with BurnMan: First, we benchmark the computations to check for correctness. Second, we exemplify two pitfalls in EoS modeling: using a different EoS than the one used to derive the mineral physical parameters or using an incorrect averaging scheme. Both pitfalls have led to incorrect conclusions on lower mantle composition and temperature in the literature. We further illustrate that fitting elastic velocities separately or jointly leads to different Mg/Si ratios for the lower mantle. However, we find that, within mineral physical uncertainties, a pyrolitic composition can match PREM very well. Finally, we find that uncertainties on specific input parameters result in a considerable amount of variation in both magnitude and gradient of the seismic velocities.

On conservation laws of Navier–Stokes Galerkin discretizations

Abstract We study conservation properties of Galerkin methods for the incompressible Navier–Stokes equations, without the divergence constraint strongly enforced. In typical discretizations such as the mixed finite element method, the conservation of mass is enforced only weakly, and this leads to discrete solutions which may not conserve energy, momentum, angular momentum, helicity, or vorticity, even though the physics of the Navier–Stokes equations dictate that they should. We aim in this work to construct discrete formulations that conserve as many physical laws as possible without utilizing a strong enforcement of the divergence constraint, and doing so leads us to a new formulation that conserves each of energy, momentum, angular momentum, enstrophy in 2D, helicity and vorticity (for reference, the usual convective formulation does not conserve most of these quantities). Several numerical experiments are performed, which verify the theory and test the new formulation.

High accuracy mantle convection simulation through modern numerical methods – II: realistic models and problems

Computations have helped elucidate the dynamics of Earths mantle for several decades already. The numerical methods that underlie these simulations have greatly evolved within this time span, and today include dynamically changing and adaptively refined meshes, sophisticated and efficient solvers, and parallelization to large clusters of computers. At the same time, many of these methods -- discussed in detail in a previous paper in this series -- were developed and tested primarily using model problems that lack many of the complexities that are common to the realistic models our community wants to solve today. With several years of experience solving complex and realistic models, we here revisit some of the algorithm designs of the earlier paper and discuss the incorporation of more complex physics. In particular, we re-consider time stepping and mesh refinement algorithms, evaluate approaches to incorporate compressibility, and discuss dealing with strongly varying material coefficients, latent heat, and how to track chemical compositions and heterogeneities. Taken together and implemented in a high-performance, massively parallel code, the techniques discussed in this paper then allow for high resolution, 3d, compressible, global mantle convection simulations with phase transitions, strongly temperature dependent viscosity and realistic material properties based on mineral physics data.

Unconditional long-time stability of a velocity---vorticity method for the 2D Navier---Stokes equations

We prove unconditional long-time stability for a particular velocity–vorticity discretization of the 2D Navier–Stokes equations. The scheme begins with a formulation that uses the Lamb vector to couple the usual velocity–pressure system to the vorticity dynamics equation, and then discretizes with the finite element method in space and implicit–explicit BDF2 in time, with the vorticity equation decoupling at each time step. We prove the method’s vorticity and velocity are both long-time stable in the

Finitely Presented Groups Related to Kaplansky’s Direct Finiteness Conjecture

L^2

Massively parallel finite element programming

L2 and