BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Charles M. Elliott

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

ResearcherUnclaimed source record
math.APmath.NAmath.OCNumerical AnalysisMethodology

Catalog footprint

What is connected

20works
5topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

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

20 published item(s)

preprint2022arXiv

Numerical analysis for the interaction of mean curvature flow and diffusion on closed surfaces

An evolving surface finite element discretisation is analysed for the evolution of a closed two-dimensional surface governed by a system coupling a generalised forced mean curvature flow and a reaction--diffusion process on the surface, inspired by a gradient flow of a coupled energy. Two algorithms are proposed, both based on a system coupling the diffusion equation to evolution equations for geometric quantities in the velocity law for the surface. One of the numerical methods is proved to be convergent in the $H^1$ norm with optimal-order for finite elements of degree at least two. We present numerical experiments illustrating the convergence behaviour and demonstrating the qualitative properties of the flow: preservation of mean convexity, loss of convexity, weak maximum principles, and the occurrence of self-intersections.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Reconciling Bayesian and perimeter regularization for binary inversion

A central theme in classical algorithms for the reconstruction of discontinuous functions from observational data is perimeter regularization via the use of the total variation. On the other hand, sparse or noisy data often demands a probabilistic approach to the reconstruction of images, to enable uncertainty quantification; the Bayesian approach to inversion, which itself introduces a form of regularization, is a natural framework in which to carry this out. In this paper the link between Bayesian inversion methods and perimeter regularization is explored. In this paper two links are studied: (i) the maximum a posteriori (MAP) objective function of a suitably chosen Bayesian phase-field approach is shown to be closely related to a least squares plus perimeter regularization objective; (ii) sample paths of a suitably chosen Bayesian level set formulation are shown to possess finite perimeter and to have the ability to learn about the true perimeter.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

A Variational Approach to Particles in Lipid Membranes

A variety of models for the membrane-mediated interaction of particles in lipid membranes, mostly well-established in theoretical physics, is reviewed from a mathematical perspective. We provide mathematically consistent formulations in a variational framework, relate apparently different modelling approaches in terms of successive approximation, and investigate existence and uniqueness. Numerical computations illustrate that the new variational formulations are directly accessible to effective numerical methods.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Double obstacle phase field approach to an inverse problem for a discontinuous diffusion coefficient

We propose a double obstacle phase field approach to the recovery of piece-wise constant diffusion coefficients for elliptic partial differential equations. The approach to this inverse problem is that of optimal control in which we have a quadratic fidelity term to which we add a perimeter regularisation weighted by a parameter sigma. This yields a functional which is optimised over a set of diffusion coefficients subject to a state equation which is the underlying elliptic PDE. In order to derive a problem which is amenable to computation the perimeter functional is relaxed using a gradient energy functional together with an obstacle potential in which there is an interface parameter epsilon. This phase field approach is justified by proving Gamma-convergence to the functional with perimeter regularisation as epsilon tends to zero. The computational approach is based on a finite element approximation. This discretisation is shown to converge in an appropriate way to the solution of the phase field problem. We derive an iterative method which is shown to yield an energy decreasing sequence converging to a discrete critical point. The efficacy of the approach is illustrated with numerical experiments.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

On algorithms with good mesh properties for problems with moving boundaries based on the Harmonic Map Heat Flow and the DeTurck trick

In this paper, we present a general approach to obtain numerical schemes with good mesh properties for problems with moving boundaries, that is for evolving submanifolds with boundaries. This includes moving domains and surfaces with boundaries. Our approach is based on a variant of the so-called the DeTurck trick. By reparametrizing the evolution of the submanifold via solutions to the harmonic map heat flow of manifolds with boundary, we obtain a new velocity field for the motion of the submanifold. Moving the vertices of the computational mesh according to this velocity field automatically leads to computational meshes of high quality both for the submanifold and its boundary. Using the ALE-method in [16], this idea can be easily built into algorithms for the computation of physical problems with moving boundaries.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

On Approximations of the Curve Shortening Flow and of the Mean Curvature Flow based on the DeTurck trick

In this paper we discuss novel numerical schemes for the computation of the curve shortening and mean curvature flows that are based on special reparametrizations. The main idea is to use special solutions to the harmonic map heat flow in order to reparametrize the equations of motion. This idea is widely known from the Ricci flow as the DeTurck trick. By introducing a variable time scale for the harmonic map heat flow, we obtain families of numerical schemes for the reparametrized flows. For the curve shortening flow this family unveils a surprising geometric connection between the numerical schemes in [5] and [9]. For the mean curvature flow we obtain families of schemes with good mesh properties similar to those in [3]. We prove error estimates for the semi-discrete scheme of the curve shortening flow. The behaviour of the fully-discrete schemes with respect to the redistribution of mesh points is studied in numerical experiments. We also discuss possible generalizations of our ideas to other extrinsic flows.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Small deformations of Helfrich energy minimising surfaces with applications to biomembranes

In this paper we introduce a mathematical model for small deformations induced by external forces of closed surfaces that are minimisers of Helfrich-type energies. Our model is suitable for the study of deformations of cell membranes induced by the cytoskeleton. We describe the deformation of the surface as a graph over the undeformed surface. A new Lagrangian and the associated Euler-Lagrange equations for the height function of the graph are derived. This is the natural generalisation of the well known linearisation in the Monge gauge for initially flat surfaces. We discuss energy perturbations of point constraints and point forces acting on the surface. We establish existence and uniqueness results for weak solutions on spheres and on tori. Algorithms for the computation of numerical solutions in the general setting are provided. We present numerical examples which highlight the behaviour of the surface deformations in different settings at the end of the paper.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Well-posedness of a fractional porous medium equation on an evolving surface

We investigate the existence, uniqueness, and $L^1$-contractivity of weak solutions to a porous medium equation with fractional diffusion on an evolving hypersurface. To settle the existence, we reformulate the equation as a local problem on a semi-infinite cylinder, regularise the porous medium nonlinearity and truncate the cylinder. Then we pass to the limit first in the truncation parameter and then in the nonlinearity, and the identification of limits is done using the theory of subdifferentials of convex functionals. In order to facilitate all of this, we begin by studying (in the setting of closed Riemannian manifolds and Sobolev spaces) the fractional Laplace--Beltrami operator which can be seen as the Dirichlet-to-Neumann map of a harmonic extension problem. A truncated harmonic extension problem will also be examined and convergence results to the solution of the harmonic extension will be given. For a technical reason, we will also consider some related extension problems on evolving hypersurfaces which will provide us with the minimal time regularity required on the harmonic extensions in order to properly formulate the moving domain problem. This functional analytic theory is of course independent of the fractional porous medium equation and will be of use generally in the analysis of fractional elliptic and parabolic problems on manifolds.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

A Stefan problem on an evolving surface

We formulate a Stefan problem on an evolving hypersurface and study the well-posedness of weak solutions given $L^1$ data. To do this, we first develop function spaces and results to handle equations on evolving surfaces in order to give a natural treatment of the problem. Then we consider the existence of solutions for $L^\infty$ data; this is done by regularisation of the nonlinearity. The regularised problem is solved by a fixed point theorem and then uniform estimates are obtained in order to pass to the limit. By using a duality method we show continuous dependence which allows us to extend the results to $L^1$ data.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

An abstract framework for parabolic PDEs on evolving spaces

We present an abstract framework for treating the theory of well-posedness of solutions to abstract parabolic partial differential equations on evolving Hilbert spaces. This theory is applicable to variational formulations of PDEs on evolving spatial domains including moving hypersurfaces. We formulate an appropriate time derivative on evolving spaces called the material derivative and define a weak material derivative in analogy with the usual time derivative in fixed domain problems; our setting is abstract and not restricted to evolving domains or surfaces. Then we show well-posedness to a certain class of parabolic PDEs under some assumptions on the parabolic operator and the data.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

On some linear parabolic PDEs on moving hypersurfaces

We consider existence and uniqueness for several examples of linear parabolic equations formulated on moving hypersurfaces. Specifically, we study in turn a surface heat equation, an equation posed on a bulk domain, a novel coupled bulk-surface system and an equation with a dynamic boundary condition. In order to prove the well-posedness, we make use of an abstract framework presented in a recent work by the authors which dealt with the formulation and well-posedness of linear parabolic equations on arbitrary evolving Hilbert spaces. Here, after recalling all of the necessary concepts and theorems, we show that the abstract framework can applied to the case of evolving (or moving) hypersurfaces, and then we demonstrate the utility of the framework to the aforementioned problems.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Phase field methods for binary recovery

We consider the inverse problem of recovering a binary function from blurred and noisy data. Such problems arise in many applications, for example image processing and optimal control of PDEs. Our formulation is based on the Mumford-Shah model, but with a phase field approximation to the perimeter regularisation. We use a double obstacle potential as well as a smooth double well potential. We introduce an iterative method for solving the problem, develop a suitable discretisation of this iterative method, and prove some convergence results. Numerical simulations are presented which illustrate the usefulness of the approach and the relative merits of the phase field models.

0 citations0 reviewsNo code linkedOpen access