Eigensolutions of nonlinear wave equations in one dimension

A class of nonlinear equations describing the steady propagation of a disturbance on the infinite interval in one dimensional space are shown under certain conditions to admit solution with a unique velocity of propagation. The class of equations describe both initial and final homogeneous steady states which are asymptotically stable with respect to uniform perturbations, in contrast to the Fisher equation, which does not.     

Mechanisms for stabilisation and destabilisation of systems of reaction-diffusion equations

Potential mechanisms for stabilising and destabilising the spatially uniform steady states of systems of reaction-diffusion equations are examined. In the first instance the effect of introducing small periodic perturbations of the diffusion coefficients in a general system of reaction-diffusion equations is studied. Analytical results are proved for the case where the uniform steady state is marginally stable and demonstrate that the effect on the original system of such perturbations is one of stabilisation. Numerical simulations carried out on an ecological model of Levin and Segel (1976) confirm the analysis as well as extending it to the case where the perturbations are no longer small. Spatio-temporal delay is then introduced into the model. Analytical and numerical results are presented which show that the effect of the delay is to destabilise the original system.     

Traveling wave solutions of a nonlinear reaction–advection equation

 We establish the existence of traveling wave solutions for a nonlinear partial differential equation that models a logistically growing population whose movement is governed by an advective process. Conditions are presented for which traveling wave solutions exist and for which they are stable to small semi-finite domain perturbations. The wave is of mathematical interest because its behavior is determined by a singular differential equation and those with small speed of propagation steepen into a shock-like solutions. Finally, we indicate that the smoothing presence of diffusion allows wave persistence when an advective slow moving wave may collapse. Received: 24 November 1997 / Revised version: 13 July 1998     

Global stability in discrete population models with delayed-density dependence

We address the global stability issue for some discrete population models with delayed-density dependence. Applying a new approach based on the concept of the generalized Yorke conditions, we establish several criteria for the convergence of all solutions to the unique positive steady state. Our results support the conjecture stated by Levin and May in 1976 affirming that the local asymptotic stability of the equilibrium of some delay difference equations (including Ricker's and Pielou's equations) implies its global stability. We also discuss the robustness of the obtained results with respect to perturbations of the model.     

Mapping parameter spaces of biological switches

Since the seminal 1961 paper of Monod and Jacob, mathematical models of biomolecular circuits have guided our understanding of cell regulation. Model-based exploration of the functional capabilities of any given circuit requires systematic mapping of multidimensional spaces of model parameters. Despite significant advances in computational dynamical systems approaches, this analysis remains a nontrivial task. Here, we use a nonlinear system of ordinary differential equations to model oocyte selection in Drosophila, a robust symmetry-breaking event that relies on autoregulatory localization of oocyte-specification factors. By applying an algorithmic approach that implements symbolic computation and topological methods, we enumerate all phase portraits of stable steady states in the limit when nonlinear regulatory interactions become discrete switches. Leveraging this initial exact partitioning and further using numerical exploration, we locate parameter regions that are dense in purely asymmetric steady states when the nonlinearities are not infinitely sharp, enabling systematic identification of parameter regions that correspond to robust oocyte selection. This framework can be generalized to map the full parameter spaces in a broad class of models involving biological switches.     

Spatial patterns in population conflicts

We consider a dynamical model for evolutionary games, and enquire how the introduction of diffusion may lead to the formation of stationary spatially inhomogeneous solutions, that is patterns. For the model equations being used it is already known that if there is an evolutionarily stable strategy (ESS), then it is stable. Equilibrium solutions which are not ESS's and which are stable with respect to spatially constant perturbations may be unstable for certain choices of the dispersal rates. We prove by a bifurcation technique that under appropriate conditions the instability leads to patterns. Computations using a curve-following technique show that the bifurcations exhibit a rich structure with loops joined by symmetry-breaking branches.     

The management of pesticide resistance.

We solve the two-species Volterra type non-linear reactive-diffusive growth equations in one spatial dimension which is taken to be infinite. We find periodic solutions which oscillate about the constant (space independent) equilibrium solutions and are stable with respect to the non-linear perturbations for small values of the diffusion constants.     

Numerical bifurcation analysis of delay differential equations arising from physiological modeling

This paper has a dual purpose. First, we describe numerical methods for continuation and bifurcation analysis of steady state solutions and periodic solutions of systems of delay differential equations with an arbitrary number of fixed, discrete delays. Second, we demonstrate how these methods can be used to obtain insight into complex biological regulatory systems in which interactions occur with time delays: for this, we consider a system of two equations for the plasma glucose and insulin concentrations in a diabetic patient subject to a system of external assistance. The model has two delays: the technological delay of the external system, and the physiological delay of the patient's liver. We compute stability of the steady state solution as a function of two parameters, compare with analytical results and compute several branches of periodic solutions and their stability. These numerical results allow to infer two categories of diabetic patients for which the external system has different efficiency.     

Harmless delays in model systems

An autonomous system of two differential equations with an asymptotically stable steady state is considered. It is shown that there exists a class of time delays which preserves the stability of such a steady state. The charactristic return time (or decay rate) for the solutions to approach the steady state is estimated using the variational system for the perturbations.     

Phase resetting and bifurcation in the ventricular myocardium.

With the dynamic differential equations of Beeler, G. W., and H. Reuter (1977, J. Physiol. [Lond.]. 268:177-210), we have studied the oscillatory behavior of the ventricular muscle fiber stimulated by a depolarizing applied current I app. The dynamic solutions of BR equations revealed that as I app increases, a periodic repetitive spiking mode appears above the subthreshold I app, which transforms to a periodic spiking-bursting mode of oscillations, and finally to chaos near the suprathreshold I app (i.e., near the termination of the periodic state). Phase resetting and annihilation of repetitive firing in the ventricular myocardium were demonstrated by a brief current pulse of the proper magnitude applied at the proper phase. These phenomena were further examined by a bifurcation analysis. A bifurcation diagram constructed as a function of I app revealed the existence of a stable periodic solution for a certain range of current values. Two Hopf bifurcation points exist in the solution, one just above the lower periodic limit point and the other substantially below the upper periodic limit point. Between each periodic limit point and the Hopf bifurcation, the cell exhibited the coexistence of two different stable modes of operation; the oscillatory repetitive firing state and the time-independent steady state. As in the Hodgkin-Huxley case, there was a low amplitude unstable periodic state, which separates the domain of the stable periodic state from the stable steady state. Thus, in support of the dynamic perturbation methods, the bifurcation diagram of the BR equation predicts the region where instantaneous perturbations, such as brief current pulses, can send the stable repetitive rhythmic state into the stable steady state.     

Travelling waves of a delayed SIR epidemic model with nonlinear incidence rate and spatial diffusion

This paper is concerned with the existence of travelling waves to a SIR epidemic model with nonlinear incidence rate, spatial diffusion and time delay. By analyzing the corresponding characteristic equations, the local stability of a disease-free steady state and an endemic steady state to this system under homogeneous Neumann boundary conditions is discussed. By using the cross iteration method and the Schauder's fixed point theorem, we reduce the existence of travelling waves to the existence of a pair of upper-lower solutions. By constructing a pair of upper-lower solutions, we derive the existence of a travelling wave connecting the disease-free steady state and the endemic steady state. Numerical simulations are carried out to illustrate the main results.     

Pattern sensitivity to boundary and initial conditions in reaction-diffusion models

We consider Turing-type reaction-diffusion equations and study (via computer simulations) how the relationship between initial conditions and the asymptotic steady state solutions varies as a function of the boundary conditions. The results indicate that boundary conditions which are non-homogeneous with respect to the kinetic steady state give rise to spatial patterns which are much less sensitive to variations in the initial conditions than those obtained with homogeneous boundary conditions, such as zero flux conditions. We also compare linear pattern predictions with the numerical solutions of the full nonlinear problem.This work supported in part by U.S. Army Grant DAJA 37-81-C-0220 and the Science and Engineering Research Council of Great Britain Grant GR/c/63595     

Bifurcation analysis of a nonlinear field model of regeneration

The analysis of bifurcating solutions in the Totafurno and Trainor [23] model of supernumerary limb production in salamanders is re-examined using the symmetry analysis developed by Totafurno [22]. In particular, we show analytically that the appearance of field solutions possessing 2 and 4 singularities (the 2- and 4-centered solutions, respectively) also correspond to true bifurcations with reduced symmetries, just as had been previously found for a solution to the field equations not possessing such singularities (the twist solution). While the results have significance primarily for the biological problem, this work serves as an instructive example of the application of symmetry groups to the bifurcation analysis of nonlinear field equations arising from a variational principle. The relationship between the solutions of the nonlinear equations and the corresponding linear equations is discussed.Supported by the National Sciences and Engineering Research Council and the Medical Research Council of CanadaTo whom correspondence should be sent     

Stability analysis for a peri-implant osseointegration model

We investigate stability of the solution of a set of partial differential equations, which is used to model a peri-implant osseointegration process. For certain parameter values, the solution has a ‘wave-like’ profile, which appears in the distribution of osteogenic cells, osteoblasts, growth factor and bone matrix. This ‘wave-like’ profile contradicts experimental observations. In our study we investigate the conditions, under which such profile appears in the solution. Those conditions are determined in terms of model parameters, by means of linear stability analysis, carried out at one of the constant solutions of the simplified system. The stability analysis was carried out for the reduced system of PDE’s, of which we prove, that it is equivalent to the original system of equations, with respect to the stability properties of constant solutions. The conclusions, derived from the linear stability analysis, are extended for the case of large perturbations. If the constant solution is unstable, then the solution of the system never converges to this constant solution. The analytical results are validated with finite element simulations. The simulations show, that stability of the constant solution could determine the behavior of the solution of the whole system, if certain initial conditions are considered.     

Body forces and pressures in elastic models of the myocardium.

Tension strands are introduced to represent active myocardial fibers. They create one body force proportional to the divergence of the tension-direction vector, and a second equal to the tension divided by the radius of curvature. Explicit solutions to isotropic linearly elastic tensor equations with these body forces are found for the radially-symmetric, linearly-isotropic, elastic spherical heart with arbitrary radial body force. They confirm experiments showing supraluminal intramural pressures. Such pressures may affect coronary perfusion. A tension strand model which is a reasonable compromise between actual myofibrillar geometry and analytical simplicity is the iso-oblique, terminating, nonintersecting model. The body force from that or any other axially symmetric body force can be the forcing term for equations in which the heart is modeled as a thin, ellipsoidal, elastic membrane.     

Feedback-mediated dynamics in a model of a compliant thick ascending limb

The tubuloglomerular feedback (TGF) system in the kidney, which is a key regulator of filtration rate, has been shown in physiologic experiments in rats to mediate oscillations in tubular fluid pressure and flow, and in NaCl concentration in the tubular fluid of the thick ascending limb (TAL). In this study, we developed a mathematical model of the TGF system that represents NaCl transport along a TAL with compliant walls. The model was used to investigate the dynamic behaviors of the TGF system. A bifurcation analysis of the TGF model equations was performed by deriving and finding roots of the characteristic equation, which arises from a linearization of the model equations. Numerical simulations of the full model equations were conducted to assist in the interpretation of the bifurcation analysis. These techniques revealed a complex parameter region that allows a variety of qualitatively different model solutions: a regime having one stable, time-independent steady-state solution; regimes having one stable oscillatory solution only; and regimes having multiple possible stable oscillatory solutions. Model results suggest that the compliance of the TAL walls increases the tendency of the model TGF system to oscillate.     

Multistable Dynamics Mediated by Tubuloglomerular Feedback in a Model of Coupled Nephrons

To help elucidate the causes of irregular tubular flow oscillations found in the nephrons of spontaneously hypertensive rats (SHR), we have conducted a bifurcation analysis of a mathematical model of two nephrons that are coupled through their tubuloglomerular feedback (TGF) systems. This analysis was motivated by a previous modeling study which predicts that NaCl backleak from a nephron’s thick ascending limb permits multiple stable oscillatory states that are mediated by TGF (Layton et al. in Am. J. Physiol. Renal Physiol. 291:F79–F97, 2006); that prediction served as the basis for a comprehensive, multifaceted hypothesis for the emergence of irregular flow oscillations in SHR. However, in that study, we used a characteristic equation obtained via linearization from a single-nephron model, in conjunction with numerical solutions of the full, nonlinear model equations for two and three coupled nephrons. In the present study, we have derived a characteristic equation for a model of any finite number of mutually coupled nephrons having NaCl backleak. Analysis of that characteristic equation for the case of two coupled nephrons has revealed a number of parameter regions having the potential for differing stable dynamic states. Numerical solutions of the full equations for two model nephrons exhibit a variety of behaviors in these regions. Some behaviors exhibit a degree of complexity that is consistent with our hypothesis for the emergence of irregular oscillations in SHR.