Resources for and from Research, Teaching and Studying
Recent Submissions
This paper is concerned with the application of Finite Element Methods (FEM) and Newton-Multigrid solvers to simulate thixotropic flows using quasi-Newtonian modeling. The thixotropy phenomena are introduced to yield stress material by taking into consideration the in-ternal material microstructure using a structure parameter. Firstly, the viscoplastic stress is modified to include the thixotropy throughout the structure parameter. Secondly, an evolution equation for the struc-ture parameter...
We discuss how ‘parallel-in-space & simultaneous-in-time’ Newton-multigrid approaches can be designed which improve the scaling behavior of the spatial parallelism by reducing the latency costs. The idea is to solve many time steps at once and therefore solving fewer but larger systems. These large systems are reordered and interpreted as a space-only problem leading to multigrid algorithm with semi-coarsening in space and line smoothing in time direction. The smoother is further improved by ...
We consider a Gierer-Meinhardt system on a surface coupled with aparabolic PDE in the bulk, the domain confined by this surface. Such a model was recently proposed and analyzed for two-dimensional bulk domains by Gomez, Ward and Wei (SIAM J. Appl. Dyn. Syst. 18, 2019).We prove the well-posedness of the bulk-surface system in arbitrary space dimensions and show that solutions remain uniformly bounded in parabolic Hölder spaces for all times. The proof uses Schauders fixed po...
This work presents an enriched Galerkin (EG) discretization for the two-dimensional shallow-water equations. The EG finite element spaces are obtained by extending the approximation spaces of the classical finite elements by discontinuous functions supported on elements. The simplest EG space is constructed by enriching the piecewise linear continuous Galerkin space with discontinuous, element-wise constant functions. Similarly to discontinuous Galerkin (DG) discretizations, the EG scheme is ...
In this work we present a new approach for coupled CFD-Optics problems that consists of a combination of a Finite Element Method (FEM) based flow solver with a ray tracing based tool for optic forces that are induced by a laser. This is a setup that is mainly encountered in the field of optical traps. We combined the open-source computational fluid dynamics (CFD) package FEATFLOW with the raytracing software of the LAT-RUB with this task in mind. We benchmark and analyz...
Building on recent advances in the analysis and design of algebraic flux correction (AFC) schemes, new differentiable limiter functions are constructed for efficient nonlinear solution strategies. The proposed scaling parameters are used to limit artificial diffusion operators incorporated into the residual of a high order target scheme to produce accurate and bound-preserving finite element approximations to hyperbolic problems. Due to this stabilization procedure, the occurring system becom...
It is an open question if the threshold condition θ < θ_* for the Dörfler marking parameter is necessary to obtain optimal algebraic rates of adaptive finite element methods. We present a (non-PDE) example fitting into the common abstract convergence framework (axioms of adaptivity) and which is potentially converging with exponential rates. However, for Dörfler marking θ > θ_* the algebraic converges rate can be made arbitrarily small.
The paper is concerned with an optimal control problem governed by the rate-independent system of quasi-static perfect elasto-plasticity. The objective is optimize the displacement field in the domain occupied by the body by means of prescribed Dirichlet boundary data, which serve as control variables. The arising optimization problem is nonsmooth for several reasons, in particular, since the control-to-state mapping is not single-valued. We therefore apply a Yosida regularization to obtain a ...
This paper is concerned with necessary optimality conditions for optimal control problems governed by variational inequalities of the second kind. So-called strong stationarity conditions are derived in an abstract framework. Strong stationarity conditions are regarded as the most rigorous ones, since they imply all other types of stationarity concepts and are equivalent to purely primal optimality conditions. The abstract framework is afterwards applied to four application-driven exa...
We approximate the solution of the Stokes equations by a new quasi-optimal and pressure robust discontinuous Galerkin discretization of arbitrary order. This means quasi-optimality of the velocity error independent of the pressure. Moreover, the discretization is well-defined for any load which is admissible for the continuous problem and it also provides classical quasioptimal estimates for the sum of velocity and pressure errors. The key design principle is a careful discretization of the l...
In this work, we introduce algebraic flux correction schemes for enriched (P1 ⊕ P0 and Q1 ⊕ P0) Galerkin discretizations of the linear advection equation. The piecewise-constant component stabilizes the continuous Galerkin approximation without introducing free parameters. However, violations of discrete maximum principles are possible in the neighborhood of discontinuities and steep fronts. To keep the cell averages and the degrees of freedom of the continuous P1/Q1 component in the admissibl...
The paper is concerned with an optimal control problem governed by the rate-independent system of quasi-static perfect elasto-plasticity. The objective is to optimize the stress field by controlling the displacement at prescribed parts of the boundary. The control thus enters the system in the Dirichlet boundary conditions. Therefore, the safe load condition is automatically fulfilled so that the system admits a solution, whose stress field is unique. This gives rise to a well defined control-to-...
To benefit from current trends in HPC hardware, such as increasing avail-ability of low precision hardware, we present the concept of prehandling as a direct way of preconditioning and the hierarchical finite element method which is exceptionally well-suited to apply prehandling to Poisson-like problems, at least in 1D and 2D. Such problems are known to cause ill-conditioned stiffness matrices and therefore high computational errors due to round-off. We show by means of numerical results that...
This paper is focused on the aspects of limiting in residual distribution (RD) schemes for high-order finite element approximations to advection problems. Both continuous and discontinuous Galerkin methods are considered in this work. Discrete maximum principles are enforced using algebraic manipula-tions of element contributions to the global nonlinear system. The required modifications can be carried out without calculating the element matrices and assembling their global counterparts. The co...
We present a time-simultaneous multigrid scheme for parabolic equations that is motivated by blocking multiple time steps together. The resulting method is closely related to multigrid waveform relaxation and is robust with respect to the spatial and temporal grid size and the number of simultaneously computed time steps. We give an intuitive understanding of the convergence behavior and briefly discuss how the theory for multigrid waveform relaxation can be applied in some special cases. Fin...
In this paper we present a holistic software approach based on the FEAT3 software for solving multidimensional PDEs with the Finite Element Method that is built for a maximum of performance, scalability, maintainability and extensibilty. We introduce basic paradigms how modern computational hardware architectures such as GPUs are exploited in a numerically scalable fashion. We show, how the framework is extended to make even the most recent advances on the hardware market accessible to the fr...
In this work, the novel “Tensor Diffusion” approach for simulating viscoelastic fluids is proposed, which is based on the idea, that the extra-stress tensor in the momentum equation of the flow model is replaced by a product of the strain-rate tensor and a tensor-valued viscosity. At least for simple flows, this approach offers the possibility to reduce the full nonlinear viscoelastic model to a generalized “Tensor Stokes” problem, avoiding the need of considering a separate stress tensor in ...
In this work, the novel "Tensor Diffusion" approach for simulating viscoelastic fluids is proposed, which is based on the idea, that the extra-stress tensor in the momentum equation of the flow model is replaced by a product of the strain-rate tensor and a tensor-valued viscosity. At least for simple flows, this approach offers the possibility to reduce the full nonlinear viscoelastic model to a generalized "Tensor Stokes" problem, avoiding the need of considering a separate stress tensor in ...