Search

In this paper we study semi-discrete and fully discrete evolving surface finite element schemes for the Cahn-Hilliard equation with a logarithmic potential. Specifically we consider linear finite elements discretising space and backward Euler time discretisation. Our analysis relies on a specific geometric assumption on the evolution of the surface. Our main results are $L^2_{H^1}$ error bounds for both the semi-discrete and fully discrete schemes, and we provide some numerical results., Comment: 60 pages, 4 figures, 1 table

We consider the existence of suitable weak solutions to the Cahn-Hilliard equation with a non-constant (degenerate) mobility on a class of evolving surfaces. We also show weak-strong uniqueness for the case of a positive mobility function, and under some further assumptions on the initial data we show uniqueness for a class of strong solutions for a degenerate mobility function., Comment: 33 pages, fixed a few mistakes and added more exposition where suitable

We consider the surface Stokes equation with Lagrange multiplier and approach it numerically. Using a Taylor-Hood surface finite element method, along with an appropriate estimation for the additional Lagrange multiplier, we establish optimal convergence results for the velocity in $H^1$ and $L^2$, and in $L^2$ for the two pressures. Furthermore, we present a new inf-sup condition to help with the stability and convergence results. This formulation offers the advantage of simplified implementation without approximating geometric quantities, although it requires a higher approximation of the parameterized surface. In addition, we provide numerical simulations that confirm the established error bounds and perform a comparative analysis against the penalty approach., Comment: 42 pages, 10 figures

We focus on the derivation and analysis of a model for multi-component phase separation occurring on biological membranes, inspired by observations of lipid raft formation. The model integrates local membrane composition with local membrane curvature, describing the membrane's geometry through a perturbation method represented as a graph over an undeformed Helfrich minimising surface, such as a sphere. The resulting energy consists of a small deformation functional coupled to a Cahn-Hilliard functional. By applying Onsager's variational principle, we obtain a multi-component Cahn-Hilliard equation for the vector $\boldsymbol\varphi$ of protein concentrations coupled to an evolution equation for the small deformation $u$ along the normal direction to the reference membrane. Then, in the case of a constant mobility matrix, we consider the Cauchy problem and we prove that it is (globally) well posed in a weak setting. We also demonstrate that any weak solution regularises in finite time and satisfies the so-called "strict separation property". This property allows usto show that any weak solution converges to a single stationary state by a suitable version of the Lojasiewicz-Simon inequality.

We derive various novel free boundary problems as limits of a coupled bulk-surface reaction-diffusion system modelling ligand-receptor dynamics on evolving domains. These limiting free boundary problems may be formulated as Stefan-type problems on an evolving hypersurface. Our results are new even in the setting where there is no domain evolution. The models are of particular relevance to a number of applications in cell biology. The analysis utilises $L^\infty$-estimates in the manner of De Giorgi iterations and other technical tools, all in an evolving setting. We also report on numerical simulations.

We study two fully discrete evolving surface finite element schemes for the Cahn-Hilliard equation on an evolving surface, given a smooth potential with polynomial growth. In particular we establish optimal order error bounds for a (fully implicit) backward Euler time-discretisation, and an implicit-explicit time-discretisation, with isoparametric surface finite elements discretising space., Comment: 42 pages, 5 figures (added new section on sharp interface numerics)

We derive a system of equations which can be seen as an evolving surface version of the diffuse interface "Model H" of Hohenberg and Halperin (1977). We then consider the well-posedness for the corresponding (tangential) system when one prescribes the evolution of the surface. Well-posedness is proved for smooth potentials in the Cahn-Hilliard equation with polynomial growth, and also for a thermodynamically relevant singular potential., Comment: 59 pages. Fixed even more typos, and corrected the long time existence result

Bulk-surface systems on evolving domains are studied. Such problems appear typically from modelling receptor-ligand dynamics in biological cells. Our first main result is the global existence and boundedness of solutions in all dimensions. This is achieved by proving $L^p$-maximal regularity of parabolic equations and duality methods in moving surfaces, which are of independent interest. The second main result is the large time dynamics where we show, under the assumption that the volume/area of the moving domain/surface is unchanged and that the material velocities are decaying for large time, that the solution converges to a unique spatially homogeneous equilibrium. The result is proved by extending the entropy method to bulk-surface systems in evolving domains., Comment: 44 pages. Some typos are corrected. Comments are very welcome

The aim of this paper is to develop a numerical scheme to approximate evolving interface problems for parabolic equations based on the abstract evolving finite element framework proposed in (C M Elliott, T Ranner, IMA J Num Anal, 41:3, 2021, doi:10.1093/imanum/draa062). An appropriate weak formulation of the problem is derived for the use of evolving finite elements designed to accommodate for a moving interface. Optimal order error bounds are proved for arbitrary order evolving isoparametric finite elements. The paper concludes with numerical results for a model problem verifying orders of convergence., Comment: 32 pages, 5 figures, 1 table

We consider the Cahn-Hilliard equation with constant mobility and logarithmic potential on a two-dimensional evolving closed surface embedded in $\mathbb R^3$, as well as a related weighted model. The well-posedness of weak solutions for the corresponding initial value problems on a given time interval $[0,T]$ have already been established by the first two authors. Here we first prove some regularisation properties of weak solutions in finite time. Then, we show the validity of the strict separation property for both the problems. This means that the solutions stay uniformly away from the pure phases $\pm1$ from any positive time on. This property plays an essential role to achieve higher-order regularity for the solutions. Also, it is a rigorous validation of the standard double-well approximation. The present results are a twofold extension of the well-known ones for the classical equation in planar domains.

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.

We propose and unify classes of different models for information propagation over graphs. In a first class, propagation is modelled as a wave which emanates from a set of \emph{known} nodes at an initial time, to all other \emph{unknown} nodes at later times with an ordering determined by the arrival time of the information wave front. A second class of models is based on the notion of a travel time along paths between nodes. The time of information propagation from an initial \emph{known} set of nodes to a node is defined as the minimum of a generalised travel time over subsets of all admissible paths. A final class is given by imposing a local equation of an eikonal form at each \emph{unknown} node, with boundary conditions at the \emph{known} nodes. The solution value of the local equation at a node is coupled to those of neighbouring nodes with lower values. We provide precise formulations of the model classes and prove equivalences between them. Finally we apply the front propagation models on graphs to semi-supervised learning via label propagation and information propagation on trust networks.

We develop a functional framework suitable for the treatment of partial differential equations and variational problems on evolving families of Banach spaces. We propose a definition for the weak time derivative that does not rely on the availability of a Hilbertian structure and explore conditions under which spaces of weakly differentiable functions (with values in an evolving Banach space) relate to classical Sobolev--Bochner spaces. An Aubin--Lions compactness result is proved. We analyse concrete examples of function spaces over time-evolving spatial domains and hypersurfaces for which we explicitly provide the definition of the time derivative and verify isomorphism properties with the aforementioned Sobolev--Bochner spaces. We conclude with the proof of well posedness for a class of nonlinear monotone problems on an abstract evolving space (generalising the evolutionary $p$-Laplace equation on a moving domain or surface) and identify some additional problems that can be formulated with the setting developed in this work.

We describe a functional framework suitable to the analysis of the Cahn-Hilliard equation on an evolving surface whose evolution is assumed to be given \textit{a priori}. The model is derived from balance laws for an order parameter with an associated Cahn-Hilliard energy functional and we establish well-posedness for general regular potentials, satisfying some prescribed growth conditions, and for two singular nonlinearities -- the thermodynamically relevant logarithmic potential and a double obstacle potential. We identify, for the singular potentials, necessary conditions on the initial data and the evolution of the surfaces for global-in-time existence of solutions, which arise from the fact that integrals of solutions are preserved over time, and prove well-posedness for initial data on a suitable set of admissible initial conditions. We then briefly describe an alternative derivation leading to a model that instead preserves a weighted integral of the solution, and explain how our arguments can be adapted in order to obtain global-in-time existence without restrictions on the initial conditions. Some illustrative examples and further research directions are given in the final sections.

We consider sharp interface asymptotics for a phase field model of two phase near spherical biomembranes involving a coupling between the local mean curvature and the local composition proposed by the first and second authors. The model is motivated by lipid raft formation. We introduce a reduced diffuse interface energy depending only on the membrane composition and derive the $\Gamma-$limit. We demonstrate that the Euler-Lagrange equations for the limiting functional and the sharp interface energy coincide. Finally, we consider a system of gradient flow equations with conserved Allen-Cahn dynamics for the phase field model. Performing a formal asymptotic analysis we obtain a system of gradient flow equations for the sharp interface energy coupling geodesic curvature flow for the phase interface to a fourth order PDE free boundary problem for the surface deformation.

We consider a model of a biomembrane with attached proteins. The membrane is represented by a near spherical continuous surface and attached proteins are described as discrete rigid structures which attach to the membrane at a finite number of points. The resulting surface minimises a quadratic elastic energy (obtained by a perturbation of the Canham-Helfrich energy) subject to the point constraints which are imposed by the attachment of the proteins. We calculate the derivative of the energy with respect to protein configurations. The proteins are constrained to move tangentially by translation and by rotation in the axis normal to a reference point. Previous studies have typically restricted themselves to a nearly flat membrane and circular inclusions. A numerically accessible representation of this derivative is derived and employed in some numerical experiments.

We derive and analyse an energy to model lipid raft formation on biological membranes involving a coupling between the local mean curvature and the local composition. We apply a perturbation method recently introduced by Fritz, Hobbs and the first author to describe the geometry of the surface as a graph over an undeformed Helfrich energy minimising surface. The result is a surface Cahn-Hilliard functional coupled with a small deformation energy We show that suitable minimisers of this energy exist and consider a gradient flow with conserved Allen-Cahn dynamics, for which existence and uniqueness results are proven. Finally, numerical simulations show that for the long time behaviour raft-like structures can emerge and stablise, and their parameter dependence is further explored., Comment: 27 pages, 11 figures

We formulate a well-posedness and approximation theory for a class of generalised saddle point problems with a specific form of constraints. In this way we develop an approach to a class of fourth order elliptic partial differential equations with point constraints using the idea of splitting into coupled second order equations. An approach is formulated using a penalty method to impose the constraints. Our main motivation is to treat certain fourth order equations involving the biharmonic operator and point Dirichlet constraints for example arising in the modelling of biomembranes on curved and flat surfaces but the approach may be applied more generally. The theory for well-posedness and approximation is presented in an abstract setting. Several examples are described together with some numerical experiments.

In this contribution to the proceedings of the 11th Mathematical Society of Japan (MSJ) Seasonal Institute (July 2018) we give an overview of some recent work on a mathematical model for small deformations of a spherical membrane. The idea is to consider perturbations to minimisers of a surface geometric energy. The model is obtained from consideration of second order approximations to a perturbed energy. In particular, the considered problems involve particle constraints and surface phase field energies., Comment: Submission to the proceedings of the 11th Mathematical Society of Japan (MSJ) Seasonal Institute (July 2018)

In this article, we analyse the domain mapping method approach to approximate statistical moments of solutions to linear elliptic partial differential equations posed over random geometries including smooth surfaces and bulk-surface systems. In particular, we present the necessary geometric analysis required by the domain mapping method to reformulate elliptic equations on random surfaces onto a fix deterministic surface using a prescribed stochastic parametrisation of the random domain. An abstract analysis of a finite element discretisation coupled with a Monte-Carlo sampling is presented for the resulting elliptic equations with random coefficients posed over the fixed curved reference domain and optimal error estimates are derived. The results from the abstract framework are applied to a model elliptic problem on a random surface and a coupled elliptic bulk-surface system and the theoretical convergence rates are confirmed by numerical experiments.

We propose a double obstacle phase field methodology for binary recovery of the slowness function of an Eikonal equation found in first traveltime tomography. We treat the inverse problem as an optimization problem with quadratic misfit functional added to a phase field relaxation of the perimeter penalization functional. Our approach yields solutions as we account for well posedness of the forward problem by choosing regular priors. We obtain a convergent finite difference and mixed finite element based discretization and a well defined descent scheme by accounting for the non-differentiability of the forward problem. We validate the phase field technique with a $\Gamma$ - convergence result and numerically by conducting parameter studies for the scheme, and by applying it to a variety of test problems with different geometries, boundary conditions, and source - receiver locations.

We consider the well-posedness and numerical approximation of a Hamilton--Jacobi equation on an evolving hypersurface in $\mathbb R^3$. Definitions of viscosity sub- and supersolutions are extended in a natural way to evolving hypersurfaces and provide uniqueness by comparison. An explicit in time monotone numerical approximation is derived on evolving interpolating triangulated surfaces. The scheme relies on a finite volume discretisation which does not require acute triangles. The scheme is shown to be stable and consistent leading to an existence proof via the proof of convergence. Finally an error bound is proved of the same order as in the flat stationary case.

We formulate a well-posedness and approximation theory for a class of generalised saddle point problems. In this way we develop an approach to a class of fourth order elliptic partial differential equations using the idea of splitting into coupled second order equations. Our main motivation is to treat certain fourth order surface equations arising in the modelling of biomembranes but the approach may be applied more generally. In particular, we are interested in equations with non-smooth right hand sides and operators which have non-trivial kernels.The theory for well posedness and approximation is presented in an abstract setting. Several examples are described together with some numerical experiments.

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., Comment: 30 pages, 15 figures

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.

In this paper, we introduce and analyse a surface finite element discretization of advection-diffusion equations with uncertain coefficients on evolving hypersurfaces. After stating unique solvability of the resulting semi-discrete problem, we prove optimal error bounds for the semi-discrete solution and Monte Carlo samplings of its expectation in appropriate Bochner spaces. Our theoretical findings are illustrated by numerical experiments in two and three space dimensions.

We prove existence, uniqueness, and regularity for a reaction-diffusion system of coupled bulk-surface equations on a moving domain modelling receptor-ligand dynamics in cells. The nonlinear coupling between the three unknowns is through the Robin boundary condition for the bulk quantity and the right hand sides of the two surface equations. Our results are new even in the non-moving setting, and in this case we also show exponential convergence to a steady state. The primary complications in the analysis are indeed the nonlinear coupling and the Robin boundary condition. For the well posedness and essential boundedness of solutions we use several De Giorgi-type arguments, and we also develop some useful estimates to allow us to apply a Steklov averaging technique for time-dependent operators to prove that solutions are strong. Some of these auxiliary results presented in this paper are of independent interest by themselves., Comment: 34 pages, 1 figure

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.

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.

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.

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.

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., Comment: Accepted for publication in SIAM Journal of Mathematical Analysis (SIMA)

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.

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., Comment: 38 pages

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.

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., Comment: 21 pages

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., Comment: 32 pages. Section 3 recalls the main results and notation of our earlier work arXiv:1403.4500

We consider an elliptic optimal control problem where the objective functional contains an integral along a surface of codimension 1, also known as a hypersurface. In particular, we use a fidelity term that encourages the state to take certain values along a curve in 2D or a surface in 3D. In the discretisation of this problem, which uses piecewise linear finite elements, we allow the hypersurface to be approximated e.g. by a polyhedral hypersurface. This can lead to simpler numerical methods, however it complicates the numerical analysis. We prove a priori $L^2$ error estimates for the control and present numerical results that agree with these. A comparison is also made to point control problems.

We consider an elliptic optimal control problem where the objective functional contains evaluations of the state at a finite number of points. In particular, we use a fidelity term that encourages the state to take certain values at these points, which means our problem is related to ones with state constraints at points. The analysis and numerical analysis differs from when the fidelity is in the $L^2$ norm because we need the state space to embed into the space of continuous functions. In this paper we discretise the problem using two different piecewise linear finite element methods. For each discretisation we use two different approaches to prove a priori $L^2$ error estimates for the control. We discuss the differences between these methods and approaches and present numerical results that agree with our analytical results.

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.

In this paper we study time-periodic solutions to advection-diffusion equations of a scalar quantity $u$ on a periodically moving $n$-dimensional hypersurface $\Gamma(t) \subset \mathbb{R}^{n+1}$. We prove existence and uniqueness of solutions in suitable H\"older 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., Comment: 38 pages. Minor typos corrected

We consider an arbitrary-Lagrangian-Eulerian evolving surface finite element method for the numerical approximation of advection and diffusion of a conserved scalar quantity on a moving surface. We describe the method, prove optimal order error bounds and present numerical simulations that agree with the theoretical results.

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, $\Gamma \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, $\Gamma_{h}$, of $\Gamma$. 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.

Reverberation mapping (RM) is an important technique in studies of active galactic nuclei (AGN). The key idea of RM is to measure the time lag $\tau$ between variations in the continuum emission from the accretion disc and subsequent response of the broad line region (BLR). The measurement of $\tau$ is typically used to estimate the physical size of the BLR and is combined with other measurements to estimate the black hole mass $M_{\rm BH}$. A major difficulty with RM campaigns is the large amount of data needed to measure $\tau$. Recently, Fine et al (2012) introduced a new approach to RM where the BLR light curve is sparsely sampled, but this is counteracted by observing a large sample of AGN, rather than a single system. The results are combined to infer properties of the sample of AGN. In this letter we implement this method using a hierarchical Bayesian model and contrast this with the results from the previous stacked cross-correlation technique. We find that our inferences are more precise and allow for more straightforward interpretation than the stacked cross-correlation results., Comment: Accepted for publication in MNRAS Letters. 5 pages, 5 figures. Source code for this paper is available at https://github.com/eggplantbren/RMHB

We present a novel parameter identification algorithm for the estimation of parameters in models of cell motility using imaging data of migrating cells. Two alternative formulations of the objective functional that measures the difference between the computed and observed data are proposed and the parameter identification problem is formulated as a minimisation problem of nonlinear least squares type. A Levenberg-Marquardt based optimisation method is applied to the solution of the minimisation problem and the details of the implementation are discussed. A number of numerical experiments are presented which illustrate the robustness of the algorithm to parameter identification in the presence of large deformations and noisy data and parameter identification in three dimensional models of cell motility. An application to experimental data is also presented in which we seek to identify parameters in a model for the monopolar growth of fission yeast cells using experimental imaging data., Comment: 31 Pages, 13 Figures

Lateral diffusion of molecules on surfaces plays a very important role in various biological processes, including lipid transport across the cell membrane, synaptic transmission and other phenomena such as exo- and endocytosis, signal transduction, chemotaxis and cell growth. In many cases, the surfaces can possess spatial inhomogeneities and/or be rapidly changing shape. Using a generalisation of the model for a thermally excited Helfrich elastic membrane, we consider the problem of lateral diffusion on quasi-planar surfaces, possessing both spatial and temporal fluctuations. Using results from homogenisation theory, we show that, under the assumption of scale separation between the characteristic length and time scales of the membrane fluctuations and the characteristic scale of the diffusing particle, the lateral diffusion process can be well approximated by a Brownian motion on the plane with constant diffusion tensor $D$ which depends in a highly nonlinear way on the detailed properties of the surface. The effective diffusion tensor will depend on the relative scales of the spatial and temporal fluctuations and, for different scaling regimes, we prove the existence of a macroscopic limit in each case., Comment: 56 pages, 9 figures, submitted to J. Nonlin. Sci

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.

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.

We formulate stochastic partial differential equations on Riemannian manifolds, moving surfaces, general evolving Riemannian manifolds (with appropriate assumptions) and Riemannian manifolds with random metrics, in the variational setting of the analysis to stochastic partial differential equations. Considering mainly linear stochastic partial differential equations, we establish various existence and uniqueness theorems., Comment: 42 pages, 1 figure