Source author record

Chandrasekhar Venkataraman

Chandrasekhar Venkataraman 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

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

12 published item(s)

preprint2020arXiv

Multiscale analysis and simulation of a signalling process with surface diffusion

We present and analyse a model for cell signalling processes in biological tissues. The model includes diffusion and nonlinear reactions on the cell surfaces, and both inter- and intracellular signalling. Using techniques from the theory of two-scale convergence as well the unfolding method, we show convergence of the solutions to the model to solutions of a two-scale macroscopic problem. We also present a two-scale bulk-surface finite element method for the approximation of the macroscopic model. We report on some benchmarking results as well as numerical simulations in a biologically relevant regime that illustrate the influence of cell-scale heterogeneities on macroscopic concentrations.

preprint2016arXiv

A Bayesian approach to parameter identification with an application to Turing systems

We present a Bayesian methodology for infinite as well as finite dimensional parameter identification for partial differential equation models. The Bayesian framework provides a rigorous mathematical framework for incorporating prior knowledge on uncertainty in the observations and the parameters themselves, resulting in an approximation of the full probability distribution for the parameters, given the data. Although the numerical approximation of the full probability distribution is computationally expensive, parallelised algorithms can make many practically relevant problems computationally feasible. The probability distribution not only provides estimates for the values of the parameters, but also provides information about the inferability of parameters and the sensitivity of the model. This information is crucial when a mathematical model is used to study the outcome of real-world experiments. Keeping in mind the applicability of our approach to tackle real-world practical problems with data from experiments, in this initial proof of concept work, we apply this theoretical and computational framework to parameter identification for a well studied semilinear reaction-diffusion system with activator-depleted reaction kinetics, posed on evolving and stationary domains.

preprint2016arXiv

A computational approach for mode isolation for reaction-diffusion systems on arbitrary geometries

In this article we present a computational framework for isolating spatial patterns arising in the steady states of reaction-diffusion systems. Such systems have been used to model many different phenomena in areas such as developmental and cancer biology, cell motility and material science. Often one is interested in identifying parameters which will lead to a particular pattern. To attempt to answer this, we compute eigenpairs of the Laplacian on a variety of domains and use linear stability analysis to determine parameter values for the system that will lead to spatially inhomogeneous steady states whose patterns correspond to particular eigenfunctions. This method has previously been used on domains and surfaces where the eigenvalues and eigenfunctions are found analytically in closed form. Our contribution to this methodology is that we numerically compute eigenpairs on arbitrary domains and surfaces. Here we present various examples and demonstrate that mode isolation is straightforward especially for low eigenvalues. Additionally we see that if two or more eigenvalues are in a permissible range then the inhomogeneous steady state can be a linear combination of the respective eigenfunctions. Finally we show an example which suggests that pattern formation is robust on similar surfaces in cases that the surface either has or does not have a boundary.

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

Lumped finite element method for reaction-diffusion systems on compact surfaces

We propose and analyse a novel surface finite element method that preserves the invariant regions of systems of semilinear parabolic equations on closed compact surfaces in $\mathbb{R}^3$ under discretisation. We also provide a fully-discrete scheme by applying the implicit-explicit (IMEX) Euler method in time. We prove the preservation of the invariant rectangles of the continuous problem under spatial and full discretizations. For scalar equations, these results reduce to the well-known discrete maximum principle. Furthermore, we prove optimal error bounds for the semi- and fully-discrete methods, that is the convergence rates are quadratic in the meshsize and linear in the timestep. Numerical experiments are provided to support the theoretical findings. In particular we provide examples in which, in the absence of lumping, the numerical solution violates the invariant region leading to blow-up due to the nature of the kinetics.

preprint2015arXiv

Stability analysis and simulations of coupled bulk-surface reaction-diffusion systems

In this article we formulate new models for coupled systems of bulk-surface reaction-diffusion equations on stationary volumes. The bulk reaction-diffusion equations are coupled to the surface reaction-diffusion equations through linear Robin-type boundary conditions. We then state and prove the necessary conditions for diffusion-driven instability for the coupled system. Due to the nature of the coupling between bulk and surface dynamics, we are able to decouple the stability analysis of the bulk and surface dynamics. Under a suitable choice of model parameter values, the bulk reaction-diffusion system can induce patterning on the surface independent of whether the surface reaction-diffusion system produces or not, patterning. On the other hand, the surface reaction-diffusion system can not generate patterns everywhere in the bulk in the absence of patterning from the bulk reaction-diffusion system. For this case, patterns can only be induced in regions close to the surface membrane. Various numerical experiments are presented to support our theoretical findings. Our most revealing numerical result is that, Robin-type boundary conditions seem to introduce a boundary layer coupling the bulk and surface dynamics.

preprint2015arXiv

Whole cell tracking through the optimal control of geometric evolution laws

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.

preprint2013arXiv

Implicit-explicit timestepping with finite element approximation of reaction-diffusion systems on evolving domains

We present and analyse an implicit-explicit timestepping procedure with finite element spatial approximation for a semilinear reaction-diffusion systems on evolving domains arising from biological models, such as Schnakenberg's (1979). We employ a Lagrangian formulation of the model equations which permits the error analysis for parabolic equations on a fixed domain but introduces technical difficulties, foremost the space-time dependent conductivity and diffusion. We prove optimal-order error estimates in the $\Lp{\infty}(0,T;\Lp{2}(\W))$ and $\Lp{2}(0,T;\Hil{1}(\W))$ norms, and a pointwise stability result. We remark that these apply to Eulerian solutions. Details on the implementation of the Lagrangian and the Eulerian scheme are provided. We also report on a numerical experiment for an application to pattern formation on an evolving domain.

preprint2013arXiv

Parameter identification problems in the modelling of cell motility

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.

preprint2010arXiv

Global existence for semilinear reaction-diffusion systems on evolving domains

We present global existence results for solutions of reaction-diffusion systems on evolving domains. Global existence results for a class of reaction-diffusion systems on fixed domains are extended to the same systems posed on spatially linear isotropically evolving domains. The results hold without any assumptions on the sign of the growth rate. The analysis is valid for many systems that commonly arise in the theory of pattern formation. We present numerical results illustrating our theoretical findings.