Source author record

Anna Marciniak-Czochra

Anna Marciniak-Czochra 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

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

20 published item(s)

preprint2022arXiv

$Γ$-limit for a sharp interface model related to pattern formation on biomembranes

We derive a macroscopic limit for a sharp interface version of a model proposed in [29] to investigate pattern formation due to competition of chemical and mechanical forces in biomembranes. We identify sub- and supercrital parameter regimes and show with the introduction of the autocorrelation function that the ground state energy leads to the isoperimetric problem in the subcritical regime, which is interpreted to not form fine scale patterns.

preprint2022arXiv

Asymptotic behaviour of a structured population model on a space of measures

In this paper we consider a physiologically structured population model with distributed states at birth, formulated on the space of non-negative Radon measures. Using a characterisation of the pre-dual space of bounded Lipschitz functions, we show how to apply the theory of strongly continuous positive semigroups to such a model. In particular, we establish the exponential convergence of solutions to a one-dimensional global attractor.

preprint2022arXiv

Discontinuous stationary solutions to certain reaction-diffusion systems

Systems consisting of a single ordinary differential equation coupled with one reaction-diffusion equation in a bounded domain and with the Neumann boundary conditions are studied in the case of particular nonlinearities from the Brusselator model, the Gray-Scott model, the Oregonator model and a certain predator-prey model. It is shown that the considered systems have the both smooth and discontinuous stationary solutions, however, only discontinuous ones can be stable.

preprint2020arXiv

Dynamical properties of feedback signalling in B lymphopoiesis: A mathematical modelling approach

Haematopoiesis is the process of generation of blood cells. Lymphopoiesis generates lymphocytes, the cells in charge of the adaptive immune response. Disruptions of this process are associated with diseases like leukaemia, which is especially incident in children. The characteristics of self-regulation of this process make them suitable for a mathematical study. In this paper we develop mathematical models of lymphopoiesis using currently available data. We do this by drawing inspiration from existing structured models of cell lineage development and integrating them with paediatric bone marrow data, with special focus on regulatory mechanisms. A formal analysis of the models is carried out, giving steady states and their stability conditions. We use this analysis to obtain biologically relevant regions of the parameter space and to understand the dynamical behaviour of B-cell renovation. Finally, we use numerical simulations to obtain further insight into the influence of proliferation and maturation rates on the reconstitution of the cells in the B line. We conclude that a model including feedback regulation of cell proliferation represents a biologically plausible depiction for B-cell reconstitution in bone marrow. Research into haematological disorders could benefit from a precise dynamical description of B lymphopoiesis.

preprint2020arXiv

Local asymptotic stability of a system of integro-differential equations describing clonal evolution of a self-renewing cell population under mutation

In this paper we consider a system of non-linear integro-differential equations (IDEs) describing evolution of a clonally heterogeneous population of malignant white blood cells (leukemic cells) undergoing mutation and clonal selection. We prove existence and uniqueness of non-trivial steady states and study their asymptotic stability. The results are compared to those of the system without mutation. Existence of equilibria is proved by formulating the steady state problem as an eigenvalue problem and applying a version of the Krein-Rutmann theorem for Banach lattices. The stability at equilibrium is analysed using linearisation and the Weinstein-Aronszajn determinant which allows to conclude local asymptotic stability.

preprint2016arXiv

Instability of Turing patterns in reaction-diffusion-ODE systems

The aim of this paper is to contribute to the understanding of the pattern formation phenomenon in reaction-diffusion equations coupled with ordinary differential equations. Such systems of equations arise, for example, from modeling of interactions between cellular processes such as cell growth, differentiation or transformation and diffusing signaling factors. We focus on stability analysis of solutions of a prototype model consisting of a single reaction-diffusion equation coupled to an ordinary differential equation. We show that such systems are very different from classical reaction-diffusion models. They exhibit diffusion-driven instability (Turing instability) under a condition of autocatalysis of non-diffusing component. However, the same mechanism which destabilizes constant solutions of such models, destabilizes also all continuous spatially heterogeneous stationary solutions, and consequently, there exist no stable Turing patterns in such reaction-diffusion-ODE systems. We provide a rigorous result on the nonlinear instability, which involves the analysis of a continuous spectrum of a linear operator induced by the lack of diffusion in the destabilizing equation. These results are extended to discontinuous patterns for a class of nonlinearities.

preprint2016arXiv

Mass concentration in a nonlocal model of clonal selection

Self-renewal is a constitutive property of stem cells. Testing the cancer stem cell hypothesis requires investigation of the impact of self-renewal on cancer expansion. To understand better this impact, we propose a mathematical model describing dynamics of a continuum of cell clones structured by the self-renewal potential. The model is an extension of the finite multi-compartment models of interactions between normal and cancer cells in acute leukemias. It takes a form of a system of integro-differential equations with a nonlinear and nonlocal coupling, which describes regulatory feedback loops in cell proliferation and differentiation process. We show that such coupling leads to mass concentration in points corresponding to maximum of the self-renewal potential and the model solutions tend asymptotically to a linear combination of Dirac measures. Furthermore, using a Lyapunov function constructed for a finite dimensional counterpart of the model, we prove that the total mass of the solution converges to a globally stable equilibrium. Additionally, we show stability of the model in space of positive Radon measures equipped with flat metric. The analytical results are illustrated by numerical simulations.

preprint2016arXiv

Stable patterns with jump discontinuity in systems with Turing instability and hysteresis

Classical models of pattern formation are based on diffusion-driven instability (DDI) of constant stationary solutions of reaction-diffusion equations, which leads to emergence of stable, regular Turing patterns formed around that equilibrium. In this paper we show that coupling reaction-diffusion equations with ordinary differential equations (ODE) may lead to a novel pattern formation phenomenon in which DDI causes destabilization of both constant solutions and Turing patterns. Bistability and hysteresis effects in the null sets of model nonlinearities yield formation of far from the equilibrium patterns with jump discontinuity. We derive conditions for stability of stationary solutions with jump discontinuity in a suitable topology which allows disclosing the discontinuity points and leads to the definition of (ε0 , A)-stability. Additionally, we provide conditions on stability of patterns in a quasi-stationary model reduction. The analysis is illustrated on the example of three-component model of receptor-ligand binding. The proposed model provides an example of a mechanism of de novo formation of far from the equilibrium patterns in reaction-diffusion-ODE models involving co-existence of DDI and hysteresis.

preprint2015arXiv

Diffusion-driven blowup of nonnegative solutions to reaction-diffusion-ODE systems

In this paper we provide an example of a class of two reaction-diffusion-ODE equations with homogeneous Neumann boundary conditions, in which Turing-type instability not only destabilizes constant steady states but also induces blow-up of nonnegative spatially heterogeneous solutions. Solutions of this problem preserve nonnegativity and uniform boundedness of the total mass. Moreover, for the corresponding system with two non-zero diffusion coefficients, all nonnegative solutions are global in time. We prove that a removal of diffusion in one of the equations leads to a finite-time blow-up of some nonnegative spatially heterogeneous solutions.

preprint2015arXiv

The Escalator Boxcar Train method for a system of aged-structured equations

The Escalator Boxcar Train method (EBT) is a numerical method for structured population models of McKendrick-von Foerster type. Those models consist of a certain class of hyperbolic partial differential equations and describe time evolution of the distribution density of the structure variable describing a feature of individuals in the population. The method was introduced in late eighties and widely used in theoretical biology, but its convergence was proven only in recent years using the framework of measure-valued solutions. Till now the EBT method was developed only for scalar equation models. In this paper we derive a full numerical EBT scheme for age-structured, two-sex population model (Fredrickson-Hoppensteadt model), which consists of three coupled hyperbolic partial differential equations with nonlocal boundary conditions. It is the first step towards extending the EBT method to systems of structured population equations.

preprint2014arXiv

Clonal selection and therapy resistance in acute leukemias: Mathematical modelling explains different proliferation patterns at diagnosis and relapse

Recent experimental evidence suggests that acute myeloid leukemias may originate from multiple clones of malignant cells. Nevertheless it is not known how the observed clones may differ with respect to cell properties such as proliferation and self-renewal. There are scarcely any data on how these cell properties change due to chemotherapy and relapse. We propose a new mathematical model to investigate the impact of cell properties on multi-clonal composition of leukemias. Model results imply that enhanced self-renewal may be a key mechanism in the clonal selection process. Simulations suggest that fast proliferating and highly self-renewing cells dominate at primary diagnosis while relapse following therapy-induced remission is triggered mostly by highly self-renewing but slowly proliferating cells. Comparison of simulation results to patient data demonstrates that the proposed model is consistent with clinically observed dynamics based on a clonal selection process.

preprint2013arXiv

A nonlinear effective slip interface law for transport phenomena between a fracture flow and a porous medium

We present modeling of an incompressible viscous flow through a fracture adjacent to a porous medium. We consider a fast stationary flow, predominantly tangential to the porous medium. Slow flow in such setting can be described by the Beavers-Joseph-Saffman slip. For fast flows, a nonlinear filtration law in the porous medium and a non- linear interface law are expected. In this paper we rigorously derive a quadratic effective slip interface law which holds for a range of Reynolds numbers and fracture widths. The porous medium flow is described by the Darcys law. The result shows that the interface slip law can be nonlinear, independently of the regime for the bulk flow. Since most of the interface and boundary slip laws are obtained via upscaling of complex systems, the result indicates that studying the inviscid limits for the Navier-Stokes equations with linear slip law at the boundary should be rethought.

preprint2013arXiv

Efficient algorithms computing distances between Radon measures on R

In this paper numerical methods of computing distances between two Radon measures on R are discussed. Efficient algorithms for Wasserstein-type metrics are provided. In particular, we propose a novel algorithm to compute the flat metric (bounded Lipschitz distance) with a computational cost O(nlogn). The flat distance has recently proven to be adequate for the Escalator Boxcar Train (EBT) method for solving transport equations with growth terms. Therefore, finding efficient numerical algorithms to compute the flat distance between two measures is important for finding the residual error and validating empirical convergence of different methods.

preprint2013arXiv

Pattern formation in a diffusion-ODE model with hysteresis

Coupling diffusion process of signaling molecules with nonlinear interactions of intracellular processes and cellular growth/transformation leads to a system of reaction-diffusion equations coupled with ordinary differential equations (diffusion-ODE models), which differ from the usual reaction-diffusion systems. One of the mechanisms of pattern formation in such systems is based on the existence of multiple steady states and hysteresis in the ODE subsystem. Diffusion tries to average different states and is the cause of spatio-temporal patterns. In this paper we provide a systematic description of stationary solutions of such systems, having the form of transition or boundary layers. The solutions are discontinuous in the case of non-diffusing variables whose quasi-stationary dynamics exhibit hysteresis. The considered model is motivated by biological applications and elucidates a possible mechanism of formation of patterns with sharp transitions.

preprint2013arXiv

Pressure jump interface law for the Stokes-Darcy coupling: Confirmation by direct numerical simulations

It is generally accepted that the effective velocity of a viscous flow over a porous bed satisfies the Beavers-Joseph slip law. To the contrary, interface law for the effective stress has been a subject of controversy. Recently, a pressure jump interface law has been rigorously derived by Marciniak-Czochra and Mikelić. In this paper, we provide a confirmation of the analytical result using direct numerical simulation of the flow at the microscopic level.

preprint2013arXiv

Spike patterns in a reaction-diffusion-ode model with Turing instability

We explore a mechanism of pattern formation arising in processes described by a system of a single reaction-diffusion equation coupled with ordinary differential equations. Such systems of equations arise from the modeling of interactions between cellular processes and diffusing growth factors. We focused on the model of early carcinogenesis proposed by Marciniak-Czochra and Kimmel, which is an example of a wider class of pattern formation models with an autocatalytic non-diffusing component. We present a numerical study showing emergence of periodic and irregular spike patterns due to diffusion-driven instability. To control the accuracy of simulations, we develop a numerical code based on the finite element method and adaptive mesh. Simulations, supplemented by numerical analysis, indicate a novel pattern formation phenomenon based on the emergence of nonstationary structures tending asymptotically to the sum of Dirac deltas.

preprint2012arXiv

A Rigorous Derivation of the Equations for the Clamped Biot-Kirchhoff-Love Poroelastic plate

In this paper we investigate the limit behavior of the solution to quasi-static Biot's equations in thin poroelastic plates as the thickness tends to zero. We choose Terzaghi's time corresponding to the plate thickness and obtain the strong convergence of the three-dimensional solid displacement, fluid pressure and total poroelastic stress to the solution of the new class of plate equations. In the new equations the in-plane stretching is described by the 2D Navier's linear elasticity equations, with elastic moduli depending on Gassmann's and Biot's coefficients. The bending equation is coupled with the pressure equation and it contains the bending moment due to the variation in pore pressure across the plate thickness. The pressure equation is parabolic only in the vertical direction. As additional terms it contains the time derivative of the in-plane Laplacean of the vertical deflection of the plate and of the the elastic in-plane compression term.

preprint2012arXiv

Effective pressure interface law for transport phenomena between an unconfined fluid and a porous medium using homogenization

We present modeling of the incompressible viscous flows in the domain containing an unconfined fluid and a porous medium. For such setting a rigorous derivation of the Beavers-Joseph-Saffman interface condition was undertaken by Jäger and Mikelić [SIAM J. Appl. Math. \rm 60 (2000), p. 1111-1127] using the homogenization method. So far the interface law for the pressure was conceived and confirmed only numerically. In this article we justify rigorously the pressure jump condition using the corresponding boundary layer.

preprint2011arXiv

Unstable patterns in reaction-diffusion model of early carcinogenesis

Motivated by numerical simulations showing the emergence of either periodic or irregular patterns, we explore a mechanism of pattern formation arising in the processes described by a system of a single reaction-diffusion equation coupled with ordinary differential equations. We focus on a basic model of early cancerogenesis proposed by Marciniak-Czochra and Kimmel [Comput. Math. Methods Med. {\bf 7} (2006), 189--213], [Math. Models Methods Appl. Sci. {\bf 17} (2007), suppl., 1693--1719], but the theory we develop applies to a wider class of pattern formation models with an autocatalytic non-diffusing component. The model exhibits diffusion-driven instability (Turing-type instability). However, we prove that all Turing-type patterns, {\it i.e.,} regular stationary solutions, are unstable in the Lyapunov sense. Furthermore, we show existence of discontinuous stationary solutions, which are also unstable.

preprint2010arXiv

A Structured Population Model of Cell Differentiation

We introduce and analyze several aspects of a new model for cell differentiation. It assumes that differentiation of progenitor cells is a continuous process. From the mathematical point of view, it is based on partial differential equations of transport type. Specifically, it consists of a structured population equation with a nonlinear feedback loop. This models the signaling process due to cytokines, which regulate the differentiation and proliferation process. We compare the continuous model to its discrete counterpart, a multi-compartmental model of a discrete collection of cell subpopulations recently proposed by Marciniak-Czochra et al. in 2009 to investigate the dynamics of the hematopoietic system. We obtain uniform bounds for the solutions, characterize steady state solutions, and analyze their linearized stability. We show how persistence or extinction might occur according to values of parameters that characterize the stem cells self-renewal. We also perform numerical simulations and discuss the qualitative behavior of the continuous model vis a vis the discrete one.