A new mathematical model for avascular tumour growth

The early development of solid tumours has been extensively studied, both experimentally via the multicellular spheroid assay, and theoretically using mathematical modelling. The vast majority of previous models apply specifically to multicell spheroids, which have a characteristic structure of a proliferating rim and a necrotic core, separated by a band of quiescent cells. Many previous models represent these as discrete layers, separated by moving boundaries. Here, the authors develop a new model, formulated in terms of continuum densities of proliferating, quiescent and necrotic cells, together with a generic nutrient/growth factor. The model is oriented towards an in vivo rather than in vitro setting, and crucially allows for nutrient supply from underlying tissue, which will arise in the two-dimensional setting of a tumour growing within an epithelium. In addition, the model involves a new representation of cell movement, which reflects contact inhibition of migration. Model solutions are able to reproduce the classic three layer structure familiar from multicellular spheroids, but also show that new behaviour can occur as a result of the nutrient supply from underlying tissue. The authors analyse these different solution types by approximate solution of the travelling wave equations, enabling a detailed classification of wave front solutions.     

Mathematical modelling of drug transport in tumour multicell spheroids and monolayer cultures

In this paper we adapt an avascular tumour growth model to compare the effects of drug application on multicell spheroids and on monolayer cultures. The model for the tumour is based on nutrient driven growth of a continuum of live cells, whose birth and death generates volume changes described by a velocity field. The drug is modelled as an externally applied, diffusible material capable of killing cells, both linear and Michaelis-Menten kinetics for drug action on cells being studied. Numerical solutions of the resulting system of partial differential equations for the multicell spheroid case are compared with closed form solutions of the monolayer case, particularly with respect to the effects on the cell kill of the drug dosage and of the duration of its application. The results show an enhanced survival rate in multicell spheroids compared to monolayer cultures, consistent with experimental observations, and indicate that the key factor determining this is drug penetration. An analysis of the large time tumour spheroid response to a continuously applied drug at fixed concentration reveals up to three stable large time solutions, namely the trivial solution (i.e. a dead tumour), a travelling wave (continuously growing tumour) and a sublinear growth case in which cells reach a pseudo-steady-state in the core. Each of these possibilities is formulated and studied, with the bifurcations between them being discussed. Numerical solutions reveal that the pseudo-steady-state solutions persist to a significantly higher drug dose than travelling wave solutions.     

The uniqueness of wave solutions for the deterministic non-reducible n-type epidemic

In a recent paper, [8], we investigated the existence of wave solutions for a model of the deterministic non-reducible n-type epidemic. In this paper we first prove two properties left as an open question in that paper. The uniqueness of the wave solutions at all speeds for which a wave solution exists is then established. Only an exceptional case is not covered.     

Concentration wave of a solute in an artery: the influence of curvature

Mass transport and diffusion phenomena in the arterial lumen are studied through a mathematical model. Blood flow is described by the unsteady Navier-Stokes equation and solute dynamics by an advection-diffusion equation, the convective field being provided by the fluid velocity. A linearization procedure over the steady state solution is carried out and an asymptotic analysis is used to study the effect of a small curvature with respect to the straight tube. Analytical and numerical solutions are found: the results show the characteristics of the long wave propagation and the role played by the geometry on the solute distribution and demonstrate the strong influence of curvature induced by the fluid dynamics.     

具时滞的人口模型的行波解

研究具时滞的人口模型的行波解存在性问题,利用[5]中的方法,鹕行波解的存在性问题转化为寻找上下解的问题。     

一个具有年龄结构的单种群离散反应扩散模型的波前解

利用上下解方法研究了一个具有年龄结构的单种群离散反应扩散模型波前解的存在性,并证明了存在具有临界波速的波前解.     

The equation of motion for sperm flagella.

The equation of motion for sperm flagella, in which the elastic bending moment and the active contractile moment are balanced by the moment from the viscous resistance of the surrounding fluid, is solved for a wave solution that superimposes partial solutions. Substitution of the expression for the wave solution into the equation leads to an expression for the active contractile moment. This active moment can be decomposed into two parts. The first part describes an active moment that travels over the flagellum with the mechanical flagellar wave, the second part represents a moment in phase over the entire length of the flagellum, which decreases linearly towards the distal tip. The linear synchronous moment, to which an amount of traveling moment has been added as a perturbation, leads to wave solutions that closely resemble flagellar waves. Properties such as wavelength and wave amplitudes and also the shape of the waves in sea urchin sperm flagella at different frequencies are accurately described by the theory. The change in wave shape in sea urchin sperm flagella at raised viscosity is predicted well by the theory. The different wave properties caused in bull sperm flagella by different boundary conditions at the proximal junction are explained. When only a traveling active moment is present in a flagellum, the wave solutions describe waves of a small wave length in a long flagellum. Some properties of the wave motion of sperm flagella are derived from the theory and verified experimentally.     

Rogue Waves: From Nonlinear Schr?dinger Breather Solutions to Sea-Keeping Test

Under suitable assumptions, the nonlinear dynamics of surface gravity waves can be modeled by the one-dimensional nonlinear Schrödinger equation. Besides traveling wave solutions like solitons, this model admits also breather solutions that are now considered as prototypes of rogue waves in ocean. We propose a novel technique to study the interaction between waves and ships/structures during extreme ocean conditions using such breather solutions. In particular, we discuss a state of the art sea-keeping test in a 90-meter long wave tank by creating a Peregrine breather solution hitting a scaled chemical tanker and we discuss its potential devastating effects on the ship.     

Neural fields with fast learning dynamic kernel

We introduce a modified-firing-rate model based on Hebbian-type changing synaptic connections. The existence and stability of solutions such as rest state, bumps, and traveling waves are shown for this type of model. Three types of kernels, namely exponential, Mexican hat, and periodic synaptic connections, are considered. In the former two cases, the existence of a rest state solution is proved and the conditions for their stability are found. Bump solutions are shown for two kinds of synaptic kernels, and their stability is investigated by constructing a corresponding Evans function that holds for a specific range of values of the kernel coefficient strength (KCS). Applying a similar method, we consider exponential synaptic connections, where traveling wave solutions are shown to exist. Simulation and numerical analysis are presented for all these cases to illustrate the resulting solutions and their stability.     

Dispersal and competition models for plants

New models for seed dispersal and competition between plant species are formulated and analyzed. The models are integrodifference equations, discrete in time and continuous in space, and have applications to annual and perennial species. The spread or invasion of a single plant species into a geographic region is investigated by studying the travelling wave solutions of these equations. Travelling wave solutions are shown to exist in the single-species models and are compared numerically. The asymptotic wave speed is calculated for various parameter values. The single-species integrodifference equations are extended to a model for two competing annual plants. Competition in the two-species model is based on a difference equation model developed by Pakes and Maller [26]. The two-species model with competition and dispersal yields a system of integrodifference equations. The effects of competition on the travelling wave solutions of invading plant species is investigated numerically.     

A Traveling Wave Model for Invasion by Precursor and Differentiated Cells

We develop and investigate a continuum model for invasion of a domain by cells that migrate, proliferate and differentiate. The model is applicable to neural crest cell invasion in the developing enteric nervous system, but is presented in general terms and is of broader applicability. Two cell populations are identified and modeled explicitly; a population of precursor cells that migrate and proliferate, and a population of differentiated cells derived from the precursors which have impaired migration and proliferation. The equation describing the precursor cells is based on Fisher’s equation with the addition of a carrying-capacity limited differentiation term. Two variations of the proliferation term are considered and compared. For most parameter values, the model admits a traveling wave solution for each population, both traveling at the same speed. The traveling wave solutions are investigated using perturbation analysis, phase plane methods, and numerical techniques. Analytical and numerical results suggest the existence of two wavespeed selection regimes. Regions of the parameter space are characterized according to existence, shape, and speed of traveling wave solutions. Our observations may be used in conjunction with experimental results to identify key parameters determining the invasion speed for a particular biological system. Furthermore, our results may assist experimentalists in identifying the resource that is limiting proliferation of precursor cells.     

Steady states and outbreaks of two-phase nonlinear age-structured model of population dynamics with discrete time delay

In this paper we study the nonlinear age-structured model of a polycyclic two-phase population dynamics including delayed effect of population density growth on the mortality. Both phases are modelled as a system of initial boundary values problem for semi-linear transport equation with delay and initial problem for nonlinear delay ODE. The obtained system is studied both theoretically and numerically. Three different regimes of population dynamics for asymptotically stable states of autonomous systems are obtained in numerical experiments for the different initial values of population density. The quasi-periodical travelling wave solutions are studied numerically for the autonomous system with the different values of time delays and for the system with oscillating death rate and birth modulus. In both cases it is observed three types of travelling wave solutions: harmonic oscillations, pulse sequence and single pulse.     

一类互惠生态模型解的周期性和爆破

研究了变系数两种群互惠模型解的周期性和爆破性.首先利用上下解和单调迭代方法研究了互惠模型周期解的存在性、稳定性和吸引性,结果表明种群之间的竞争强于互惠时周期解存在.接着利用比较原理得到了解在有限时刻爆破的充分条件,结果表明种群之间的互惠强于之间的竞争时爆破发生.最后通过数值模拟验证了定理的正确性.     

三种群互惠模型抛物系统解的整体存在与爆破

研究了三种群互惠模型的抛物系统,用上下解方法研究解的整体存在性与爆破问题,并给出了相应的条件.结果表明当种群自身的竞争强时解整体存在,反之则有可能爆破.     

Analysis of a time-delayed mathematical model for tumour growth with an almost periodic supply of external nutrients

In this paper, the existence, uniqueness and exponential stability of almost periodic solutions for a mathematical model of tumour growth are studied. The establishment of the model is based on the reaction–diffusion dynamics and mass conservation law and is considered with a delay in the cell proliferation process. Using a fixed-point theorem in cones, the existence and uniqueness of almost periodic solutions for different parameter values of the model is proved. Moreover, by the Gronwall inequality, sufficient conditions are established for the exponential stability of the unique almost periodic solution. Results are illustrated by computer simulations.     

Early development and quorum sensing in bacterial biofilms

 We develop mathematical models to examine the formation, growth and quorum sensing activity of bacterial biofilms. The growth aspects of the model are based on the assumption of a continuum of bacterial cells whose growth generates movement, within the developing biofilm, described by a velocity field. A model proposed in Ward et al. (2001) to describe quorum sensing, a process by which bacteria monitor their own population density by the use of quorum sensing molecules (QSMs), is coupled with the growth model. The resulting system of nonlinear partial differential equations is solved numerically, revealing results which are qualitatively consistent with experimental ones. Analytical solutions derived by assuming uniform initial conditions demonstrate that, for large time, a biofilm grows algebraically with time; criteria for linear growth of the biofilm biomass, consistent with experimental data, are established. The analysis reveals, for a biologically realistic limit, the existence of a bifurcation between non-active and active quorum sensing in the biofilm. The model also predicts that travelling waves of quorum sensing behaviour can occur within a certain time frame; while the travelling wave analysis reveals a range of possible travelling wave speeds, numerical solutions suggest that the minimum wave speed, determined by linearisation, is realised for a wide class of initial conditions. Received: 10 February 2002 / Revised version: 29 October 2002 / Published online: 19 March 2003 Key words or phrases: Bacterial biofilm – Quorum sensing – Mathematical modelling – Numerical solution – Asymptotic analysis – Travelling wave analysis     

Travelling wave solutions of diffusive Lotka-Volterra equations

We establish the existence of travelling wave solutions for two reaction diffusion systems based on the Lotka-Volterra model for predator and prey interactions. For simplicity, we consider only 1 space dimension. The waves are of transition front type, analogous to the travelling wave solutions discussed by Fisher and Kolmogorov et al. for a scalar reaction diffusion equation. The waves discussed here are not necessarily monotone. For any speed c there is a travelling wave solution of transition front type. For one of the systems discussed here, there is a distinguished speed c^* dividing the waves into two types, waves of speed c < c^* being one type, waves of speed c ? c^* being of the other type. We present numerical evidence that for this system the wave of speed c^* is stable, and that c^* is an asymptotic speed of propagation in some sense. For the other system, waves of all speeds are in some sense stable. The proof of existence uses a shooting argument and a Lyapunov function. We also discuss some possible biological implications of the existence of these waves.     

Mathematical analysis of a model describing evolution of an asexual population in a changing environment

We investigate a mathematical model for an asexual population with non-overlapping (discrete) generations, that exists in a changing environment. Sexual populations are also briefly discussed at the end of the paper. It is assumed that selection occurs on the value of a single polygenic trait, which is controlled by a finite number of loci with discrete-effect alleles. The environmental change results in a moving fitness optimum, causing the trait to be subject to a combination of stabilising and directional selection.This model is different from that investigated by Waxman and Peck [Genetics 153 (1999) 1041] where overlapping generations and continuous effect alleles were considered. In this paper, we consider non-overlapping generations and discrete effect alleles. However in [Genetics 153 (1999) 1041] and the present work, there is the same pattern of environmental change, namely a constant rate of change of the optimum.From [Genetics 153 (1999) 1041], no rigorous theoretical conclusion can be drawn about the form of the solutions as t grows large. Numerical work carried out in [Genetics 153 (1999) 1041] suggests that the solution is a lagged travelling wave solution, but no mathematical proof exists for the continuous model. Only partial results, regarding existence of travelling wave solutions and perturbed solutions, have been established (see [Nonlin. Anal. 53 (2003) 683; An integral equation describing an asexual population in a changing environment, Preprint]). For the discrete case of this paper, under the assumption that the ratio between the unit of genotypic value and the speed of environment change is a rational number, we are able to give rigorous proof of the following conclusion: the population follows the environmental change with a small lag behind, moreover, the lag is represented using a calculable quantity.     

Accurate computation of the stable solitary wave for the FitzHugh-Nagumo equations

Comparisons are made between three different methods for computing the stable solitary wave solution for the FitzHugh-Nagumo equations which consist of a nonlinear diffusion equation coupled to an ordinary differential equation in time. They model the Hodgkin-Huxley equations which describe the propagation of the nerve impulse down the axon. Two of the methods involve the travelling wave equations. Previous accurate numerical computations of these equations as an initial-value problem using a shooting method lead to inaccurate values for the wave speed; however, nonlinear corrections to the initial values are shown to yield accurate values. A boundary-value method applies asymptotic boundary conditions and uses a spline-collocation code called COLSYS for numerical solution of boundary-value problems which leads to accurate wave profiles and speeds. The third method is to solve an initial-boundary-value problem with an adaptive outgoing wave condition for the partial differential equations where the solitary wave emerges as the stable long time solution. The concept of a wave integral is introduced and they are derived to determine the wave speed used in the adaptive boundary condition and to measure the closeness of the computed solutions to the exact solitary wave solution.This work was supported in part by the Natural Sciences and Engineering Research Council Canada under Grant A4559 and by the John Simon Guggenheim Memorial Foundation