Source author record

Thomas Ranner

Thomas Ranner 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

8works
4topics
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

8 published item(s)

preprint2020arXiv

A stable finite element method for low inertia undulatory locomotion in three dimensions

We present and analyse a numerical method for understanding the low-inertia dynamics of an open, inextensible viscoelastic rod - a long and thin three dimensional object - representing the body of a long, thin microswimmer. Our model allows for both elastic and viscous, bending and twisting deformations and describes the evolution of the midline curve of the rod as well as an orthonormal frame which fully determines the rod's three dimensional geometry. The numerical method is based on using a combination of piecewise linear and piecewise constant functions based on a novel rewriting of the model equations. We derive a stability estimate for the semi-discrete scheme and show that at the fully discrete level that we have good control over the length element and preserve the frame orthonormality conditions up to machine precision. Numerical experiments demonstrate both the good properties of the method as well as the applicability of the method for simulating undulatory locomotion in the low-inertia regime.

preprint2020arXiv

A unified theory for continuous in time evolving finite element space approximations to partial differential equations in evolving domains

We develop a unified theory for continuous in time finite element discretisations of partial differential equations posed in evolving domains including the consideration of equations posed on evolving surfaces and bulk domains as well coupled surface bulk systems. We use an abstract variational setting with time dependent function spaces and abstract time dependent finite element spaces. Optimal a priori bounds are shown under usual assumptions on perturbations of bilinear forms and approximation properties of the abstract finite element spaces. The abstract theory is applied to evolving finite elements in both flat and curved spaces. Evolving bulk and surface isoparametric finite element spaces defined on evolving triangulations are defined and developed. These spaces are used to define approximations to parabolic equations in general domains for which the abstract theory is shown to apply. Numerical experiments are described which confirm the rates of convergence.

preprint2016arXiv

An unfitted discontinuous Galerkin scheme for conservation laws on evolving surfaces

Motivated by considering partial differential equations arising from conservation laws posed on evolving surfaces, a new numerical method for an advection problem is developed and simple numerical tests are performed. The method is based on an unfitted discontinuous Galerkin approach where the surface is not explicitly tracked by the mesh which means the method is extremely flexible with respect to geometry. Furthermore, the discontinuous Galerkin approach is well-suited to capture the advection driven by the evolution of the surface without the need for a space-time formulation, back-tracking trajectories or streamline diffusion. The method is illustrated by a one-dimensional example and numerical results are presented that show good convergence properties for a simple test problem.

preprint2016arXiv

Coupled bulk-surface free boundary problems arising from a mathematical model of receptor-ligand dynamics

We consider a coupled bulk-surface system of partial differential equations with nonlinear coupling modelling receptor-ligand dynamics. The model arises as a simplification of a mathematical model for the reaction between cell surface resident receptors and ligands present in the extra-cellular medium. We prove the existence and uniqueness of solutions. We also consider a number of biologically relevant asymptotic limits of the model. We prove convergence to limiting problems which take the form of free boundary problems posed on the cell surface. We also report on numerical simulations illustrating convergence to one of the limiting problems as well as the spatio-temporal distributions of the receptors and ligands in a realistic geometry.

preprint2016arXiv

Solving reaction-diffusion equations on evolving surfaces defined by biological image data

We present a computational approach for solving reaction-diffusion equations on evolving surfaces which have been obtained from cell image data. It is based on finite element spaces defined on surface triangulations extracted from time series of 3D images. A model for the transport of material between the subsequent surfaces is required where we postulate a velocity in normal direction. We apply the technique to image data obtained from a spreading neutrophil cell. By simulating FRAP experiments we investigate the impact of the evolving geometry on the recovery. We find that for idealised FRAP conditions, changes in membrane geometry, easily account for differences of $\times 10$ in recovery half-times, which shows that experimentalists must take great care when interpreting membrane photobleaching results. We also numerically solve an activator -- depleted substrate system and report on the effect of the membrane movement on the pattern evolution.

preprint2015arXiv

A computational approach to an optimal partition problem on surfaces

We explore an optimal partition problem on surfaces using a computational approach. The problem is to minimise the sum of the first Dirichlet Laplace--Beltrami operator eigenvalues over a given number of partitions of a surface. We consider a method based on eigenfunction segregation and perform calculations using modern high performance computing techniques. We first test the accuracy of the method in the case of three partitions on the sphere then explore the problem for higher numbers of partitions and on other surfaces.

preprint2014arXiv

Evolving surface finite element method for the Cahn-Hilliard equation

We use the evolving surface finite element method to solve a Cahn- Hilliard equation on an evolving surface with prescribed velocity. We start by deriving the equation using a conservation law and appropriate transport for- mulae and provide the necessary functional analytic setting. The finite element method relies on evolving an initial triangulation by moving the nodes accord- ing to the prescribed velocity. We go on to show a rigorous well-posedness result for the continuous equations by showing convergence, along a subse- quence, of the finite element scheme. We conclude the paper by deriving error estimates and present various numerical examples.

preprint2014arXiv

Unfitted finite element methods using bulk meshes for surface partial differential equations

In this paper, we define new unfitted finite element methods for numerically approximating the solution of surface partial differential equations using bulk finite elements. The key idea is that the $n$-dimensional hypersurface, $Γ\subset \mathbb{R}^{n+1}$, is embedded in a polyhedral domain in $\mathbb R^{n+1}$ consisting of a union, $\mathcal{T}_h$, of $(n+1)$-simplices. The finite element approximating space is based on continuous piece-wise linear finite element functions on $\mathcal{T}_h$. Our first method is a sharp interface method, \emph{SIF}, which uses the bulk finite element space in an approximating weak formulation obtained from integration on a polygonal approximation, $Γ_{h}$, of $Γ$. The full gradient is used rather than the projected tangential gradient and it is this which distinguishes \emph{SIF} from the method of [42]. The second method, \emph{NBM}, is a narrow band method in which the region of integration is a narrow band of width $O(h)$. \emph{NBM} is similar to the method of [13]. but again the full gradient is used in the discrete weak formulation. The a priori error analysis in this paper shows that the methods are of optimal order in the surface $L^{2}$ and $H^{1}$ norms and have the advantage that the normal derivative of the discrete solution is small and converges to zero. Our third method combines bulk finite elements, discrete sharp interfaces and narrow bands in order to give an unfitted finite element method for parabolic equations on evolving surfaces. We show that our method is conservative so that it preserves mass in the case of an advection diffusion conservation law. Numerical results are given which illustrate the rates of convergence.