Source author record

Andreas Dedner

Andreas Dedner appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

11works
7topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

11 published item(s)

preprint2021arXiv

Robust nonconforming virtual element methods for general fourth order problems with varying coefficients

We present a class of nonconforming virtual element methods for general fourth order partial differential equations in two dimensions. We develop a generic approach for constructing the necessary projection operators and virtual element spaces. Optimal error estimates in the energy norm are provided for general linear fourth order problems with varying coefficients. We also discuss fourth order perturbation problems and present a novel nonconforming scheme which is uniformly convergent with respect to the perturbation parameter without requiring an enlargement of the space. Numerical tests are carried out to verify the theoretical results. We conclude with a brief discussion on how our approach can easily be applied to nonlinear fourth order problems.

preprint2020arXiv

A surface finite element method for computational modelling of cell blebbing

Cell blebs are protrusions of the cell membrane and can be instrumental for cell migration. We derive a continuum model for the mechanical and geometrical aspects of the onset of blebbing in terms of a force balance. It is abstract and flexible in that it allows for amending force contributions related to membrane tension or the presence of linker molecules between membrane and cell cortex. The deforming membrane and all forces are expressed by means of a parametrisation over a stationary reference surface. A variational formulation is presented and analysed for well-posedness. For this purpose, we derive a semi-discrete scheme based on the surface finite element method. We provide a convergence result and estimates of the error due to the spatial discretisation. Furthermore, we present a computational framework where specific models can be implemented and later on conveniently amended if desired, using a domain specific language implemented in Python. While the high level program control can be done within the Python scripting environment, the actual computationally expensive step of evolving the solution over time is carried out by binding to an efficient software backend. Cell membrane geometries given in terms of a parametrisation or obtained from image data can be accounted for. A couple of numerical simulation results illustrate the approach.

preprint2020arXiv

The DUNE Framework: Basic Concepts and Recent Developments

This paper presents the basic concepts and the module structure of the Distributed and Unified Numerics Environment and reflects on recent developments and general changes that happened since the release of the first Dune version in 2007 and the main papers describing that state [1, 2]. This discussion is accompanied with a description of various advanced features, such as coupling of domains and cut cells, grid modifications such as adaptation and moving domains, high order discretizations and node level performance, non-smooth multigrid methods, and multiscale methods. A brief discussion on current and future development directions of the framework concludes the paper.

preprint2016arXiv

Convection in a Single Column -- Modelling, Algorithm and Analysis

The group focused on a model problem of idealised moist air convection in a single column of atmosphere. Height, temperature and moisture variables were chosen to simplify the mathematical representation (along the lines of the Boussinesq approximation in a height variable defined in terms of pressure). This allowed exact simple solutions of the numerical and partial differential equation problems to be found. By examining these, we identify column behaviour, stability issues and explore the feasibility of a more general solution process.

preprint2015arXiv

Discontinuous Galerkin methods for hyperbolic and advection-dominated problems on surfaces

We extend the discontinuous Galerkin (DG) framework to the analysis of first-order hyperbolic and advection-dominated problems posed on implicitly defined surfaces. The focus will be on the hyperbolic part, which is discretised using a "discrete surface" generalisation of the jump-stabilised upwind flux. A key issue arising in the analysis (which does not appear in the planar setting) is the treatment of the discrete velocity field, choices of which play an important role in the stability of the scheme. We then prove optimal error estimates in an appropriate norm given a number of assumptions on the discrete velocity field, which are then investigated and discussed in more detail. The theoretical results are verified numerically for a number of test problems exhibiting advection-dominated behaviour.

preprint2015arXiv

Efficient Multigrid Preconditioners for Atmospheric Flow Simulations at High Aspect Ratio

Many problems in fluid modelling require the efficient solution of highly anisotropic elliptic partial differential equations (PDEs) in "flat" domains. For example, in numerical weather- and climate-prediction an elliptic PDE for the pressure correction has to be solved at every time step in a thin spherical shell representing the global atmosphere. This elliptic solve can be one of the computationally most demanding components in semi-implicit semi-Lagrangian time stepping methods which are very popular as they allow for larger model time steps and better overall performance. With increasing model resolution, algorithmically efficient and scalable algorithms are essential to run the code under tight operational time constraints. We discuss the theory and practical application of bespoke geometric multigrid preconditioners for equations of this type. The algorithms deal with the strong anisotropy in the vertical direction by using the tensor-product approach originally analysed by Börm and Hiptmair [Numer. Algorithms, 26/3 (2001), pp. 219-234]. We extend the analysis to three dimensions under slightly weakened assumptions, and numerically demonstrate its efficiency for the solution of the elliptic PDE for the global pressure correction in atmospheric forecast models. For this we compare the performance of different multigrid preconditioners on a tensor-product grid with a semi-structured and quasi-uniform horizontal mesh and a one dimensional vertical grid. The code is implemented in the Distributed and Unified Numerics Environment (DUNE), which provides an easy-to-use and scalable environment for algorithms operating on tensor-product grids. Parallel scalability of our solvers on up to 20,480 cores is demonstrated on the HECToR supercomputer.

preprint2015arXiv

The DUNE-ALUGrid Module

In this paper we present the new DUNE-ALUGrid module. This module contains a major overhaul of the sources from the ALUgrid library and the binding to the DUNE software framework. The main changes include user defined load balancing, parallel grid construction, and an redesign of the 2d grid which can now also be used for parallel computations. In addition many improvements have been introduced into the code to increase the parallel efficiency and to decrease the memory footprint. The original ALUGrid library is widely used within the DUNE community due to its good parallel performance for problems requiring local adaptivity and dynamic load balancing. Therefore, this new model will benefit a number of DUNE users. In addition we have added features to increase the range of problems for which the grid manager can be used, for example, introducing a 3d tetrahedral grid using a parallel newest vertex bisection algorithm for conforming grid refinement. In this paper we will discuss the new features, extensions to the DUNE interface, and explain for various examples how the code is used in parallel environments.

preprint2014arXiv

Adaptive discontinuous Galerkin methods on surfaces

We present a dual weighted residual-based a posteriori error estimate for a discontinuous Galerkin (DG) approximation of a linear second-order elliptic problem on compact smooth connected and oriented surfaces in $\mathbb{R}^{3}$ which are implicitly represented as level sets of a smooth function. We show that the error in the energy norm may be split into a "residual part" and a higher order "geometric part". Upper and lower bounds for the resulting a posteriori error estimator are proven and we consider a number of challenging test problems to demonstrate the reliability and efficiency of the estimator. We also present a novel "geometric" driven refinement strategy for PDEs on surfaces which considerably improves the performance of the method on complex surfaces.

preprint2014arXiv

High order discontinuous Galerkin methods on surfaces

We derive and analyze high order discontinuous Galerkin methods for second-order elliptic problems on implicitely defined surfaces in $\mathbb{R}^{3}$. This is done by carefully adapting the unified discontinuous Galerkin framework of Arnold et al. [2002] on a triangulated surface approximating the smooth surface. We prove optimal error estimates in both a (mesh dependent) energy norm and the $L^2$ norm.

preprint2013arXiv

Analysis of the discontinuous Galerkin method for elliptic problems on surfaces

We extend the discontinuous Galerkin (DG) framework to a linear second-order elliptic problem on a compact smooth connected and oriented surface. An interior penalty (IP) method is introduced on a discrete surface and we derive a-priori error estimates by relating the latter to the original surface via the lift introduced in Dziuk (1988). The estimates suggest that the geometric error terms arising from the surface discretisation do not affect the overall convergence rate of the IP method when using linear ansatz functions. This is then verified numerically for a number of test problems. An intricate issue is the approximation of the surface conormal required in the IP formulation, choices of which are investigated numerically. Furthermore, we present a generic implementation of test problems on surfaces.

preprint2013arXiv

Discontinuous Galerkin methods for nonvariational problems

We extend the finite element method introduced by Lakkis and Pryer [2011] to approximate the solution of second order elliptic problems in nonvariational form to incorporate the discontinuous Galerkin (DG) framework. This is done by viewing the NVFEM as a mixed method whereby the finite element Hessian is an auxiliary variable in the formulation. Representing the finite element Hessian in a discontinuous setting yields a linear system of the same size and having the same sparsity pattern of the compact DG methods for variational elliptic problems. Furthermore, the system matrix is very easy to assemble, Thus this approach greatly reduces the computational complexity of the discretisation compared to the continuous approach. We conduct a stability and consistency analysis making use of the unified framework set out in Arnold et. al. [2001]. We also give an apriori analysis of the method. The analysis applies to any consistent representation of the finite element Hessian, thus is applicable to the previous works making use of continuous Galerkin approximations.