Timo Heister

Timo Heister, Thomas Wick

We present a scalable, parallel implementation of a solver for the solution of a phase-field model for quasi-static brittle fracture. The code is available as open source. Numerical solutions in 2d and 3d with adaptive mesh refinement show optimal scaling of the linear solver based on algebraic multigrid, and convergence of the phase-field model towards exact values of functionals of interests such as the crack opening displacement or the total crack volume. In contrast to uniform refinement, adaptive mesh refinement allows us to recover optimal convergence rates for the non-smooth solutions encountered in typical test problems. We also present numerical studies of the influence of the finite domain size on functional evaluations used to approximate the infinite domain.

Thomas C. Clevenger, Timo Heister, Guido Kanschat et al.

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.

Sergey Charnyi, Timo Heister, Maxim A. Olshanskii et al.

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.

Timo Heister, Maxim A. Olshanskii, Leo G. Rebholz

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.

Sergey Charnyi, Timo Heister, Maxim A. Olshanskii et al.

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.