Modelling inert gas exchange in tissue and mixed-venous blood return to the lungs

Inert gas exchange in tissue has been almost exclusively modelled by using an ordinary differential equation. The mathematical model that is used to derive this ordinary differential equation assumes that the partial pressure of an inert gas (which is proportional to the content of that gas) is a function only of time. This mathematical model does not allow for spatial variations in inert gas partial pressure. This model is also dependent only on the ratio of blood flow to tissue volume, and so does not take account of the shape of the body compartment or of the density of the capillaries that supply blood to this tissue. The partial pressure of a given inert gas in mixed-venous blood flowing back to the lungs is calculated from this ordinary differential equation. In this study, we write down the partial differential equations that allow for spatial as well as temporal variations in inert gas partial pressure in tissue. We then solve these partial differential equations and compare them to the solution of the ordinary differential equations described above. It is found that the solution of the ordinary differential equation is very different from the solution of the partial differential equation, and so the ordinary differential equation should not be used if an accurate calculation of inert gas transport to tissue is required. Further, the solution of the PDE is dependent on the shape of the body compartment and on the density of the capillaries that supply blood to this tissue. As a result, techniques that are based on the ordinary differential equation to calculate the mixed-venous blood partial pressure may be in error.     

Heat Transfer in a Micropolar Fluid over a Stretching Sheet with Newtonian Heating

This article looks at the steady flow of Micropolar fluid over a stretching surface with heat transfer in the presence of Newtonian heating. The relevant partial differential equations have been reduced to ordinary differential equations. The reduced ordinary differential equation system has been numerically solved by Runge-Kutta-Fehlberg fourth-fifth order method. Influence of different involved parameters on dimensionless velocity, microrotation and temperature is examined. An excellent agreement is found between the present and previous limiting results.     

Complex dynamic behaviour of autonomous microbial food chains

 The dynamic behaviour of food chains under chemostat conditions is studied. The microbial food chain consists of substrate (non-growing resources), bacteria (prey), ciliates (predator) and carnivore (top predator). The governing equations are formulated at the population level. Yet these equations are derived from a dynamic energy budget model formulated at the individual level. The resulting model is an autonomous system of four first-order ordinary differential equations. These food chains resemble those occuring in ecosystems. Then the prey is generally assumed to grow logistically. Therefore the model of these systems is formed by three first-order ordinary differential equations. As with these ecosystems, there is chaotic behaviour of the autonomous microbial food chain under chemostat conditions with biologically relevant parameter values. It appears that the trajectories on the attractors consists of two superimposed oscillatory behaviours, a slow one for predator–top predator and a fast one for the prey–predator on one branch at which the top predator increases slowly. In some regions of the parameter space there are multiple attractors. Received 8 November 1995; received in revised form 7 January 1997     

脉冲-扩散周期竞争系统解的渐近行为

研究了在周期变化环境中具有扩散及种群密度可能发生突变的两竞争种群动力系统的数学模型．模型由反应扩散方程组以及初边值及脉冲条件组成．文章建立了研究模型的上下解方法,获得了一些比较原理．利用脉冲常微分方程的比较定理以及利用相应的脉冲常微分方程的解控制和估计所讨论模型的解,研究了系统模型的解的渐近性质．     

Optimal control of the chemotherapy of HIV

 Using an existing ordinary differential equation model which describes the interaction of the immune system with the human immunodeficiency virus (HIV), we introduce chemotherapy in an early treatment setting through a dynamic treatment and then solve for an optimal chemotherapy strategy. The control represents the percentage of effect the chemotherapy has on the viral production. Using an objective function based on a combination of maximizing benefit based on T cell counts and minimizing the systemic cost of chemotherapy (based on high drug dose/strength), we solve for the optimal control in the optimality system composed of four ordinary differential equations and four adjoint ordinary differential equations. Received 5 July 1995; received in revised form 3 June 1996     

Some mathematical aspects of chemotherapy: I. One-organ models

A model of the processes occuring in the exchange of a drug between capillary plasma, extracellular space and intracellular space is developed. This leads to an interesting set of differential differences equations, one of which is an integrodifferential equation, another a partial differential equation. Under certain conditions, these may be simplified to a set of ordinary differential equations. The application of Laplace transform techniques to the solution of these equations is discussed.     

Transport phenomena in solid state fermentation: Oxygen transport in static tray fermentors

A mathematical model has been developed for describing the oxygen concentration during the exponential growth of microorganisms, in a static solid substrate bed supported on a tray fermentor. The model equations comprise of one partial differential equation for mass transfer and an ordinary differential equation of growth. After nondimensionlisation, analytical solution to the model has been obtained by the method of Laplace transforms. An expression for critical thickness of bed is deduced from the model equation. The significance of the model in the design of tray fermentors is discussed. The validity of the discussion is verified by taking an illustration from the literature.     

Rate of excitation propagation in a reduced Hodgkins-Huxley model. III. Integrodifferential equations]

The Hodgkin - Huxley system of equations is reduced to single integral-differential equation in neglection of slow variables dynamics. Two limiting cases of fast and slow sodium activation processes are considered. The first case leads to a nonlinear differential equation for the potential, the second one - to an ordinary differential equation with a known source as a function of coordinate. Such a simplification is due to approximation of steady-state sodium activation variable with the help of Heviside function. The validity of this approximation is discussed; the corresponding error is estimated by calculation of the second approximation for the source function.     

Mathematical analysis of a proposed mechanism for oscillatory insulin secretion in perifused HIT-15 cells

Oscillatory secretion of insulin has been observed in many different experimental preparations ranging from pancreatic islets to the whole pancreas. Here we examine the mathematical features underlying a possible model for oscillatory secretion from the perifused, insulin-secreting cell line, HIT-15. The model includes the kinetics of uptake of glucose by GLUT transporters, the rate of glucose metabolism within the cell, and the effect of glucose on the rate of insulin secretion. Putative feedback by insulin on the rate of glucose transport into the cells is treated phenomenologically and leads to insulin oscillations similar to those observed experimentally in HIT cells. The resulting set of ordinary differential equations is simplified by time-scale analysis to a two-variable set of ordinary differential equations. Because of this simplification we can explore, in great detail, the characteristics of the oscillations and their sensitivity to parameter variation using phase plane analysis.     

Top predator persistence in differential equation models of food chains: The effects of omnivory and external forcing of lower trophic levels

Persistence criteria are given for the highest trophic level predator in ordinary differential equation models of food chains exhibiting arbitrary omnivory and external supplementation of food source or an intermediate predator. The results are expressed in terms of inequalities involving the bounds on the intrinsic growth and interaction rates. Whether omnivory or external forcing enhances persistence is discussed, particularly for the examples of three-, four-, and five-link Lotka-Volterra food chains.     

Exact fundamental solutions of linear parabolic equations with spatially varying coefficients

The exact general solution is obtained to a linear second order ordinary differential equation which has quite general coefficients depending on an arbitrary function of the independent variable. From this, the exact fundamental solution is derived for the corresponding linear parabolic partial differential equation with coefficients depending on the single space coordinate. In a special case this latter equation reduces to one of the Fokker-Planck type. These coefficients are then generalised and the appropriate fundamental solution is obtained. Extensions are given to linear parabolic equations in two andn space dimensions. The paper provides a collection of basic examples which illustrate and develop the theory for the generation of the exact fundamental solutions. Reduction to, and the corresponding fundamental solutions of the Fokker-Planck equations is presented, where appropriate.     

A random-walk model of human mortality and aging

Consideration is made of the roles of certain types of state space and time scales for a random-walk model of individual physiological status change and death. Because the actual measurement of physiological variables omits many variables relevant to survival, we are forced to view this model as operating in a stochastic state space for a population of individuals where only the frequency distributions are deterministic. In this stochastic state space, under the assumption that the “history” of prior movement contains no additional information, the forward partial differential equation is obtained for the distribution of a population whose movement in the selected space is determined by the randomwalk equations. If the initial distribution of the population in the state space is normal, then certain assumptions about movement and mortality will operate to preserve normality thereafter. Under the assumption of normality, simultaneous ordinary differential equations can be derived from the forward partial differential equation defining the distribution function. Examination of the ordinary simultaneous differential equations shows how parameters for certain models of aging and mortality can be obtained.     

Curvature-Driven Pore Growth in Charged Membranes during Charge-Pulse and Voltage-Clamp Experiments

We find that curvature-driven growth of pores in electrically charged membranes correctly reproduces charge-pulse experiments. Our model, consisting of a Langevin equation for the time dependence of the pore radius coupled to an ordinary differential equation for the number of pores, captures the statistics of the pore population and its effect on the membrane conductance. The calculated pore radius is a linear, and not an exponential, function of time, as observed experimentally. Two other important features of charge-pulse experiments are recovered: pores reseal for low and high voltages but grow irreversibly for intermediate values of the voltage. Our set of coupled ordinary differential equations is equivalent to the partial differential equation used previously to study pore dynamics, but permits the study of longer timescales necessary for the simulations of voltage-clamp experiments. An effective phase diagram for such experiments is obtained.     

Numerical methods and parameter estimation of a structured population model with discrete events in the life history

We consider two numerical methods for the solution of a physiologically structured population (PSP) model with multiple life stages and discrete event reproduction. The model describes the dynamic behaviour of a predator-prey system consisting of rotifers predating on algae. The nitrate limited algal prey population is modelled unstructured and described by an ordinary differential equation (ODE). The formulation of the rotifer dynamics is based on a simple physiological model for their two life stages, the egg and the adult stage. An egg is produced when an energy buffer reaches a threshold value. The governing equations are coupled partial differential equations (PDE) with initial and boundary conditions. The population models together with the equation for the dynamics of the nutrient result in a chemostat model. Experimental data are used to estimate the model parameters. The results obtained with the explicit finite difference (FD) technique compare well with those of the Escalator Boxcar Train (EBT) method. This justifies the use of the fast FD method for the parameter estimation, a procedure which involves repeated solution of the model equations.     

Current status and challenges in connecting models of erythrocyte metabolism to experimental reality

Detailed kinetic models of human erythrocyte metabolism have served to summarize the vast literature and to predict outcomes from laboratory and “Nature's” experiments on this simple cell. Mathematical methods for handling the large array of nonlinear ordinary differential equations that describe the time dependence of this system are well developed, but experimental methods that can guide the evolution of the models are in short supply. NMR spectroscopy is one method that is non-selective with respect to analyte detection but is highly specific with respect to their identification and quantification. Thus time courses of metabolism are readily recorded for easily changed experimental conditions. While the data can be simulated, the systems of equations are too complex to allow solutions of the inverse problem, namely parameter-value estimation for the large number of enzyme and membrane-transport reactions operating in situ as opposed to in vitro. Other complications with the modelling include the dependence of cell volume on time, and the rates of membrane transport processes are often dependent on the membrane potential. These matters are discussed in the light of new modelling strategies.     

Kolmogorov’s Differential Equations and Positive Semigroups on First Moment Sequence Spaces

Spatially implicit metapopulation models with discrete patch-size structure and host-macroparasite models which distinguish hosts by their parasite loads lead to infinite systems of ordinary differential equations. In several papers, a this-related theory will be developed in sufficient generality to cover these applications. In this paper the linear foundations are laid. They are of own interest as they apply to continuous-time population growth processes (Markov chains). Conditions are derived that the solutions of an infinite linear system of differential equations, known as Kolmogorov’s differential equations, induce a C ₀-semigroup on an appropriate sequence space allowing for first moments. We derive estimates for the growth bound and the essential growth bound and study the asymptotic behavior. Our results will be illustrated for birth and death processes with immigration and catastrophes. An erratum to this article can be found at     

State-dependent neutral delay equations from population dynamics

A novel class of state-dependent delay equations is derived from the balance laws of age-structured population dynamics, assuming that birth rates and death rates, as functions of age, are piece-wise constant and that the length of the juvenile phase depends on the total adult population size. The resulting class of equations includes also neutral delay equations. All these equations are very different from the standard delay equations with state-dependent delay since the balance laws require non-linear correction factors. These equations can be written as systems for two variables consisting of an ordinary differential equation (ODE) and a generalized shift, a form suitable for numerical calculations. It is shown that the neutral equation (and the corresponding ODE—shift system) is a limiting case of a system of two standard delay equations.     

Cancer chemotherapy: Optimal control using the Verhulst-Pearl equation

General (deterministic) ordinary differential equations for the representation of cancer growth are presented when the growth is perturbed due to the action of a chemotherapeutic agent. The Verhulst-Pearl equation is introduced as a particular example of a growth equation applicable to human tumors. An optimal control problem with general performance criterion and state equation is formulated and shown to possess a novel feedback control relationship. This relationship is used in two continuous drug delivery problems involving the Verhulst-Pearl equation.     

Birfurcation theory applied to a simple model of a biochemical oscillator

Some results are presented relating to the question whether self-sustained oscillations are possible in a sequence of biochemical reactions with end- point inhibition. The model used has a single nonlinear ordinary differential equation coupled to a set of linear equations, with all coefficients in the linear terms equal. The explicit algebraic form of the Hopf-Friedrich bifurcation theory is used to show that when the number of coupled equations is large enough this model has a stable periodic solution when the equilibrium point of the equations has just become unstable.     

一类分数阶微分方程的最优控制

主要研究了一类分数阶微分方程的最优控制问题.通过Oustaloup迭代逼近,可以将分数阶微分算子在频率域范围内进行近似.那么原先的分数阶微分方程就转换为一般的常微分方程.利用这个关系,可以得到关于分数阶微分方程的两个定理和一个引理.最后给出一个例子说明该方法的有效性.