BZPEERPapers, people, evidence and research pathsDiscover
Research tools
Network
Opportunities
Account
Source author record
Vanessa Styles appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.
Catalog footprint
Research graph
Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.
BZPEER is loading the nearby papers, people, topics and institutions for this page.
Published work
We consider a fully practical finite element approximation of a diffuse interface model for tumour growth that takes the form of a degenerate parabolic system. In addition to showing stability bounds for the approximation, we prove convergence, and hence existence of a solution to this system in two space dimensions. Several numerical experiments demonstrate the practicality and robustness of the proposed method.
We consider a diffuse interface approach for solving an elliptic PDE on a given closed hypersurface. The method is based on a (bulk) finite element scheme employing numerical quadrature for the phase field function and hence is very easy to implement compared to other approaches. We estimate the error in natural norms in terms of the spatial grid size, the interface width and the order of the underlying quadrature rule. Numerical test calculations are presented which confirm the form of the error bounds.
We consider a finite element approximation for a system consisting of the evolution of a curve evolving by forced curve shortening flow coupled to a reaction-diffusion equation on the evolving curve. The curve evolves inside a given domain $Ω\subset \mathbb{R}^2$ and meets $\partial Ω$ orthogonally. The scheme for the coupled system is based on the schemes derived in [BDS17] and [DE98]. We present numerical experiments and show the experimental order of convergence of the approximation.
We consider a semi-discrete finite element approximation for a system consisting of the evolution of a planar curve evolving by forced curve shortening flow inside a given bounded domain $Ω\subset \mathbb{R}^2$, such that the curve meets the boundary $\partial Ω$ orthogonally, and the forcing is a function of the solution of a reaction-diffusion equation that holds on the evolving curve. We prove optimal error bounds for the resulting approximation and present numerical experiments.
Using basic thermodynamic principles we derive a Cahn--Hilliard--Darcy model for tumour growth including nutrient diffusion, chemotaxis, active transport, adhesion, apoptosis and proliferation. The model generalises earlier models and in particular includes active transport mechanisms which ensure thermodynamic consistency. We perform a formally matched asymptotic expansion and develop several sharp interface models. Some of them are classical and some are new which for example include a jump in the nutrient density at the interface. A linear stability analysis for a growing nucleus is performed and in particular the role of the new active transport term is analysed. Numerical computations are performed to study the influence of the active transport term for specific growth scenarios.
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 consider a finite element approximation for a system consisting of the evolution of a closed planar curve by forced curve shortening flow coupled to a reaction-diffusion equation on the evolving curve. The scheme for the curve evolution is based on a parametric description allowing for tangential motion, whereas the discretisation for the PDE on the curve uses an idea from [6]. We prove optimal error bounds for the resulting fully discrete approximation and present numerical experiments. These confirm our estimates and also illustrate the advantage of the tangential motion of the mesh points in practice.
In this paper we analyze a fully discrete numerical scheme for solving a parabolic PDE on a moving surface. The method is based on a diffuse interface approach that involves a level set description of the moving surface. Under suitable conditions on the spatial grid size, the time step and the interface width we obtain stability and error bounds with respect to natural norms. Furthermore, we present test calculations that confirm our analysis.
Cell tracking algorithms which automate and systematise the analysis of time lapse image data sets of cells are an indispensable tool in the modelling and understanding of cellular phenomena. In this study we present a theoretical framework and an algorithm for whole cell tracking. Within this work we consider that "tracking" is equivalent to a dynamic reconstruction of the whole cell data (morphologies) from static image datasets. The novelty of our work is that the tracking algorithm is driven by a model for the motion of the cell. This model may be regarded as a simplification of a recently developed physically meaningful model for cell motility. The resulting problem is the optimal control of a geometric evolution law and we discuss the formulation and numerical approximation of the optimal control problem. The overall goal of this work is to design a framework for cell tracking within which the recovered data reflects the physics of the forward model. A number of numerical simulations are presented that illustrate the applicability of our approach.