Search

We present the design and implementation details of a geometric multigrid method on adaptively refined meshes for massively parallel computations. The method uses local smoothing on the refined part of the mesh. Partitioning is achieved by using a space filling curve for the leaf mesh and distributing ancestors in the hierarchy based on the leaves. We present a model of the efficiency of mesh hierarchy distribution and compare its predictions to runtime measurements. The algorithm is implemented as part of the deal.II finite-element library and as such available to the public. [ABSTRACT FROM AUTHOR]

We study discretizations of the incompressible Navier–Stokes equations, written in the newly developed energy–momentum–angular momentum conserving (EMAC) formulation. We consider linearizations of the problem, which at each time step will reduce the computational cost, but can alter the conservation properties. We show that a skew-symmetrized linearization delivers the correct balance of (only) energy and that the Newton linearization conserves momentum and angular momentum, but conserves energy only up to the nonlinear residual. Numerical tests show that linearizing with 2 Newton steps at each time step is very effective at preserving all conservation laws at once, and giving accurate answers on long time intervals. The tests also show that the skew-symmetrized linearization is significantly less accurate. The tests also show that the Newton linearization of EMAC finite element formulation compares favorably to other traditionally used finite element formulation of the incompressible Navier–Stokes equations in primitive variables. [ABSTRACT FROM AUTHOR]

This paper provides an overview of the new features of the finite element library deal.II version 9.0. [ABSTRACT FROM AUTHOR]

In this paper, we propose a new numerical method to find the ground state energy of a given physical system within the Kohn–Sham density functional theory. The h‐adaptive finite element method is adopted for spatial discretization and implemented with matrix‐free operator evaluation. The ground state energy is found by performing unconstrained minimization with non‐orthogonal orbitals using the limited memory Broyden–Fletcher–Goldfarb–Shanno (BFGS) method. A geometric multigrid preconditioner is applied to improve the convergence. The clear advantage of the proposed approach is demonstrated on selected examples by comparing the performance to other methods such as preconditioned steepest descent minimization. The proposed method provides a solid framework toward O(N) complexity for the locally adaptive real‐space solution of density functional theory with finite elements. [ABSTRACT FROM AUTHOR]

This paper provides an overview of the new features of the finite element library deal.II version 8.5. [ABSTRACT FROM AUTHOR]

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 $$L^2$$ and $$H^1$$ norms, without any timestep restriction. Moreover, our analysis avoids the use of Gronwall-type estimates, which leads us to stability bounds with only polynomial (instead of exponential) dependence on the Reynolds number. Numerical experiments are given that demonstrate the effectiveness of the method. [ABSTRACT FROM AUTHOR]

Melt generation and migration are an important link between surface processes and the thermal and chemical evolution of the Earth's interior. However, their vastly different timescales make it difficult to study mantle convection and melt migration in a unified framework, especially for 3-D global models. And although experiments suggest an increase in melt volume of up to 20 per cent from the depth of melt generation to the surface, previous computations have neglected the individual compressibilities of the solid and the fluid phase. Here, we describe our extension of the finite element mantle convection code ASPECT that adds melt generation and migration. We use the original compressible formulation of the McKenzie equations, augmented by an equation for the conservation of energy. Applying adaptive mesh refinement to this type of problems is particularly advantageous, as the resolution can be increased in areas where melt is present and viscosity gradients are high, whereas a lower resolution is sufficient in regions without melt. Together with a high-performance, massively parallel implementation, this allows for high-resolution, 3-D, compressible, global mantle convection simulations coupled with melt migration. We evaluate the functionality and potential of this method using a series of benchmarks and model setups, compare results of the compressible and incompressible formulation, and show the effectiveness of adaptive mesh refinement when applied to melt migration. Our model of magma dynamics provides a framework for modelling processes on different scales and investigating links between processes occurring in the deep mantle and melt generation and migration. This approach could prove particularly useful applied to modelling the generation of komatiites or other melts originating in greater depths. The implementation is available in the Open Source ASPECT repository. [ABSTRACT FROM AUTHOR]

We investigate a flux-preserving enforcement of inhomogeneous Dirichlet boundary conditions for velocity, u | ∂ Ω = g , for use with finite element methods for incompressible flow problems that strongly enforce mass conservation. Typical enforcement via nodal interpolation is not flux-preserving in general, and it can create divergence error even when divergence-free elements are used. We show with analysis and numerical tests that by slightly (and locally) changing nodal interpolation, the enforcement recovers flux-preservation, is optimally accurate, and delivers divergence-free solutions when used with divergence-free finite elements. [ABSTRACT FROM AUTHOR]

Quasi-static elastoplastic contact problems are ubiquitous in many industrial processes and other contexts, and their numerical simulation is consequently of great interest in accurately describing and optimizing production processes. The key component in these simulations is the solution of a single load step of a time iteration. From a mathematical perspective, the problems to be solved in each time step are characterized by the difficulties of variational inequalities for both the plastic behavior and the contact problem. Computationally, they also often lead to very large problems. In this paper, we present and evaluate a complete set of methods that are (1) designed to work well together and (2) allow for the efficient solution of such problems. In particular, we use adaptive finite element meshes with linear and quadratic elements, a Newton linearization of the plasticity, active set methods for the contact problem, and multigrid-preconditioned linear solvers. Through a sequence of numerical experiments, we show the performance of these methods. This includes highly accurate solutions of a three-dimensional benchmark problem and scaling our methods in parallel to 1024 cores and more than a billion unknowns. Copyright © 2015 John Wiley & Sons, Ltd. [ABSTRACT FROM AUTHOR]

We derive a new kind of boundary conditions for the vorticity equation with solid wall boundaries for fluid flow problems. The formulation uses a Dirichlet condition for the normal component of vorticity and Neumann type conditions for the tangential components. In a Galerkin (integral) formulation the tangential condition is natural, i.e., it is enforced by a right-hand side functional and does not impose a boundary constraint on trial and test spaces. The functional involves the pressure variable, and we discuss several velocity–vorticity formulations where the proposed condition is appropriate. Several numerical experiments are given that illustrate the validity of the new approach. [ABSTRACT FROM AUTHOR]

Abstract: In this paper the grad–div stabilization for the incompressible Navier–Stokes finite element approximations is considered from two different viewpoints: (i) as a variational multiscale approach for the pressure subgrid modeling and (ii) as a stabilization procedure of least-square type. Some new error estimates for the linearized problem with the grad–div stabilization are proved with the help of norms induced by the pressure Schur complement operator. We discuss the stabilization parameter choice arising in the frameworks of least-square and multiscale methods and consider assumptions which allow to relate both approaches. [Copyright &y& Elsevier]