Source author record

Nikolaos Sfakianakis

Nikolaos Sfakianakis 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

9works
6topics
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

9 published item(s)

preprint2022arXiv

An adaptive rectangular mesh administration and refinement technique with application in cancer invasion models

We present an administration technique for the bookkeeping of adaptive mesh refinement on (hyper-)rectangular meshes. Our technique is a unified approach for h-refinement on 1-, 2- and 3D domains, which is easy to use and avoids traversing the connectivity graph of the ancestry of mesh cells. Due to the employed rectangular mesh structure, the identification of the siblings and the neighbouring cells is greatly simplified. The administration technique is particularly designed for smooth meshes, where the smoothness is dynamically used in the matrix operations. It has a small memory footprint that makes it affordable for a wide range of mesh resolutions over a large class of problems. We present three applications of this technique, one of which addresses h-refinement and its benefits in a 2D tumour growth and invasion problem.

preprint2020arXiv

Modeling multiple taxis: tumor invasion with phenotypic heterogeneity, haptotaxis, and unilateral interspecies repellence

We provide a short review of existing models with multiple taxis performed by (at least) one species and consider a new mathematical model for tumor invasion featuring two mutually exclusive cell phenotypes (migrating and proliferating). The migrating cells perform nonlinear diffusion and two types of taxis in response to non-diffusing cues: away from proliferating cells and up the gradient of surrounding tissue. Transitions between the two cell subpopulations are influenced by subcellular (receptor binding) dynamics, thus conferring the setting a multiscale character. We prove global existence of weak solutions to a simplified model version and perform numerical simulations for the full setting under several phenotype switching and motility scenarios. We also compare (via simulations) this model with the corresponding haptotaxis-chemotaxis one featuring indirect chemorepellent production and provide a discussion about possible model extensions and mathematical challenges.

preprint2016arXiv

A multiscale approach to the migration of cancer stem cells: mathematical modelling and simulations

We propose a multiscale model of the invasion of the extracellular matrix by two types of cancer cells, the differentiated cancer cells and the cancer stem cells. We assume that the epithelial mesenchymal-like transition between them is driven primarily by the epidermal growth factors. We moreover take into account the transidifferentiation program of the cancer stem cells and the cancer associated fibroblast cells as well as the fibroblast-driven remodelling of the extracellular matrix. The proposed haptotaxis model combines the macroscopic phenomenon of the invasion of the extracellular matrix with the microscopic dynamics of the epidermal growth factors. We analyse our model in a component-wise manner and compare our findings with the literature. We investigate pathological situations regarding the epidermal growth factors and accordingly propose "mathematical-treatment" scenarios to control the aggressiveness of the tumour.

preprint2016arXiv

Energy intake functions of ectotherms and endotherms derived from their body mass growth

How animals allocate energy to different body functions is still not completely understood and a challenging topic until recently. Here, we investigate in more detail the allocation of energy intake to growth, reproduction or heat production by developing energy budget models for ectothermic and endothermic vertebrates using a mathematical approach. We calculated energy intake functions of ectotherms and endotherms derived from their body mass growth. We show that our energy budget model produces energy intake patterns and distributions as observed in ectothermic and endothermic species. Our results comply consistently with some empirical studies that in endothermic species, like birds and mammals, energy is used for heat production instead of growth. Our model additionally offers an explanation on known differences in absolute energy intake between ectothermic fish and reptiles and endothermic birds and mammals. From a mathematical point of view, the model comes in two equivalent formulations, a differential and an integral one. It is derived from a discrete level approach, and it is shown to be well-posed and to attain a unique solution for (almost) every parameter set. Numerically, the integral formulation of the model is considered as an inverse problem with unknown parameters that are estimated using a series of experiments/realistic data.

preprint2016arXiv

Numerical simulation of a contractivity based multiscale cancer invasion model

We present a problem-suited numerical method for a particularly challenging cancer invasion model. This model is a multiscale haptotaxis advection-reaction-diffusion system that describes the macroscopic dynamics of two types of cancer cells coupled with microscopic dynamics of the cells adhesion on the extracellular matrix. The difficulties to overcome arises from the non-constant advection and diffusion coefficients, a time delay term, as well as stiff reaction terms. Our numerical method is a second order finite volume implicit-explicit scheme adjusted to include a) non-constant diffusion coefficients in the implicit part, b) an interpolation technique for the time delay, and c) a restriction on the time increment for the stiff reaction terms.

preprint2015arXiv

A mathematical insight in the epithelial-mesenchymal-like transition in cancer cells and its effect in the invasion of the extracellular matrix

Current biological knowledge supports the existence of a secondary group of cancer cells within the body of the tumour that exhibits stem cell-like properties. These cells are termed Cancer Stem Cells (CSCs}, and as opposed to the more usual Differentiated Cancer Cells (DCCs), they exhibit higher motility, they are more resilient to therapy, and are able to metastasize to secondary locations within the organism and produce new tumours. The origin of the CSCs is not completely clear; they seem to stem from the DCCs via a transition process related to the Epithelial-Mesenchymal Transition (EMT) that can also be found in normal tissue. In the current work we model and numerically study the transition between these two types of cancer cells, and the resulting "ensemble" invasion of the extracellular matrix. This leads to the derivation and numerical simulation of two systems: an algebraic-elliptic system for the transition and an advection-reaction-diffusion system of Keller-Segel taxis type for the invasion.

preprint2015arXiv

An Extended Filament Based Lamellipodium Model Produces Various Moving Cell Shapes in the Presence of Chemotactic Signals

The Filament Based Lamellipodium Model (FBLM) is a two-phase two-dimensional continuum model, describing the dynamcis of two interacting families of locally parallel actin filaments (C.Schmeiser and D.Oelz, How do cells move? Mathematical modeling of cytoskeleton dynamics and cell migration. Cell mechanics: from single scale-based models to multiscale modeling. Chapman and Hall, 2010). It contains accounts of the filaments' bending stiffness, of adhesion to the substrate, and of cross-links connecting the two families. An extension of the model is presented with contributions from nucleation of filaments by branching, from capping, from contraction by actin-myosin interaction, and from a pressure-like repulsion between parallel filaments due to Coulomb interaction. The effect of a chemoattractant is described by a simple signal transduction model influencing the polymerization speed. Simulations with the extended model show its potential for describing various moving cell shapes, depending on the signal transduction procedure, and for predicting transients between nonmoving and moving states as well as changes of direction.

preprint2014arXiv

Numerical study of cancer cell invasion dynamics using adaptive mesh refinement: the urokinase model

In the present work we investigate the chemotactically and proteolytically driven tissue invasion by cancer cells. The model employed is a system of advection-reaction-diffusion equations that features the role of the serine protease urokinase-type plasminogen activator. The analytical and numerical study of this system constitutes a challenge due to the merging, emerging, and travelling concentrations that the solutions exhibit. Classical numerical methods applied to this system necessitate very fine discretization grids to resolve these dynamics in an accurate way. To reduce the computational cost without sacrificing the accuracy of the solution, we apply adaptive mesh refinement techniques, in particular h-refinement. Extended numerical experiments exhibit that this approach provides with a higher order, stable, and robust numerical method for this system. We elaborate on several mesh refinement criteria and compare the results with the ones in the literature. We prove, for a simpler version of this model, $L^\infty$ bounds for the solutions, we study the stability of its conditional steady states, and conclude that it can serve as a test case for further development of mesh refinement techniques for cancer invasion simulations.

preprint2012arXiv

Entropy dissipation of moving mesh adaptation

Non-uniform grids and mesh adaptation have been a growing part of numerical simulation over the past years. It has been experimentally noted that mesh adaptation leads not only to locally improved solution but also to numerical stability of the underlying method. There have been though only few results on the mathematical analysis of these schemes due to the lack of proper tools that incorporate both the time evolution and the mesh adaptation step of the overall algorithm. In this paper we provide a method to perform the analysis of the mesh adaptation method, including both the mesh reconstruction and evolution of the solution. We moreover employ this method to extract sufficient conditions -on the adaptation of the mesh- that stabilize a numerical scheme in the sense of the entropy dissipation.