A weakly nonlinear analysis of a model of avascular solid tumour growth

 In this paper we study a mathematical model that describes the growth of an avascular solid tumour. Our analysis concentrates on the stability of steady, radially-symmetric model solutions with respect to perturbations taken from the class of spherical harmonics. Using weakly nonlinear analysis, previous results are extended to show how the amplitudes of the asymmetric modes interact. Attention focuses on a special case for which the model equations simplify. Analysis of the simplified model equations leads to the identification of a two-parameter family of asymmetric steady solutions, the dimensions of whose stable and unstable manifolds depend on the system parameters. The asymmetric steady solutions limit the basin of attraction of the radially-symmetric steady state when it is linearly stable. On the basis of these numerical and analytical results we postulate the existence of fully nonlinear steady solutions which are stable with respect to time-dependent perturbations. Received: 25 October 1998 / Revised version: 20 June 1998     

A Biomechanical Triphasic Approach to the Transport of Nondilute Solutions in Articular Cartilage

Biomechanical models for biological tissues such as articular cartilage generally contain an ideal, dilute solution assumption. In this article, a biomechanical triphasic model of cartilage is described that includes nondilute treatment of concentrated solutions such as those applied in vitrification of biological tissues. The chemical potential equations of the triphasic model are modified and the transport equations are adjusted for the volume fraction and frictional coefficients of the solutes that are not negligible in such solutions. Four transport parameters, i.e., water permeability, solute permeability, diffusion coefficient of solute in solvent within the cartilage, and the cartilage stiffness modulus, are defined as four degrees of freedom for the model. Water and solute transport in cartilage were simulated using the model and predictions of average concentration increase and cartilage weight were fit to experimental data to obtain the values of the four transport parameters. As far as we know, this is the first study to formulate the solvent and solute transport equations of nondilute solutions in the cartilage matrix. It is shown that the values obtained for the transport parameters are within the ranges reported in the available literature, which confirms the proposed model approach.     

Evolution of discontinuous solutions to transport equations in stellarators

It was shown earlier that, in the range of rare collisions, transport equations for stellarators allow steady discontinuous solutions for the ambipolar electric field and for the plasma density and temperature gradients. Moreover, such solutions are non-single-valued; that is, their explicit form depends on the initial values of the ambipolar electric field. The time-independent transport equations are derived under the conventional quasineutrality condition; i.e., it is assumed that the electron and ion densities, N _e and N _i , are related by the relationship N _e = ZN _i (where Z is the ion charge number). In other words, the plasma charge density is assumed to be much less than the product e _i N _i . Under typical conditions, the corresponding inequality is satisfied by a large margin. However, if the electric field E has discontinuities, then it can be seen from the equation ▿·E = 4πρ that, at the discontinuity points, the charge density becomes infinite and the relationship N _e = ZN _i fails to hold, so it is necessary to replace it with N _e = ZN _i + ρ/e _e . In the transport equations, this latter replacement produces additional terms, proportional to the second radial derivative of the field E. With these additional terms, the steady solutions are modified substantially. First, the ambipolar field and the derivatives of the density and temperatures all become continuous functions of the coordinates, a result that seems to be quite obvious. The second, not-so-obvious result is that the steady solutions become single-valued, i.e., independent of the initial values of the ambipolar electric field. It turns out that, in this case, two regimes are possible, depending on the values of the plasma parameters. In the first regime, the solution is unique and is independent of the initial conditions. In the second regime, two steady solutions can exist, depending on the initial conditions. One of the solution is similar to that obtained in the first regime, and the other differs from the first one both in the ambipolar field profile and in the dependence of the density and temperatures on the minor plasma radius. It cannot be excluded that different plasma confinement modes revealed in experiments are associated with the existence of such solutions.     

Transport across homoporous and heteroporous membranes in nonideal, nondilute solutions. I. Inequality of reflection coefficients for volume flow and solute flow.

The Kirkwood formulation of the Stefan-Maxwell equations is used to develop the transport equations for a membrane bounded by nonideal, nondilute solutions. The reflection coefficients for volume flow and solute flow are not equal but are related by a simple expression that depends on the concentration of the bounding solutions. The ratio of the two coefficients is independent of heteroporous membrane structure and the thickness of adjacent boundary layers. Experimental measurements of these reflection coefficients for sucrose transport across Cuprophan verify this relationship; this indicates that the Onsager reciprocal relation, which is assumed by the theory, holds for nonideal, nondilute solutions. The two reflection coefficients may be made operationally identical by a simple redefination of the osmotic driving force.     

Isoelectric focusing using non-amphoteric buffers in free solution: I. Determination of stable concentration profiles

For large-scale separations of proteins, the use of simple non-amphoteric buffers in free solution and in multicompartment electrolyzers seems promising for industrial applications. The stabilization of a pH profile with this type of buffer requires the strict observation of two conditions: choice of an adequate buffer; stationary profiles of concentrations. During electrolysis in free solution, the ions of the buffer are displaced across the compartments by migration and by diffusion. To keep a stationary composition, the inflow and outflow of all individual ionic species through each compartment must be identical. At high current, diffusion may be neglected against migration and the ionic flows will be identical if the transport number of each ion is constant at each location within the cell. In these conditions, stationary compositions will be independent of the electric current. This condition of constant transport numbers implies the use of profiles of buffer concentrations different from those published up to now. The new equations for these profiles of concentrations are given in the present paper. The constant migration of the ions must be compensated in the end compartments of the isoelectric focusing cell to provide a stable steady state. Two methods are proposed in the literature: the buffer renewal method and the external recycling method (rheoelectrolysis). Here modified buffer renewal method is proposed. Using stationary mass balances, analytical equations are given to calculate the flows and the composition of the solutions to be recycled or added. Using these equations and the profiles of concentrations to keep constant transport numbers, it is demonstrated that only a renewal of the buffers in the end compartments may lead to stable pH profiles and thus to valid conditions of separation.     

Solution of a set of Maxwell-Lorentz equations for a ring relativistic electron beam

Nonlinear solutions to a set of Maxwell’s equations and the relativistic equations of electron motion are obtained that describe the equilibrium of a high-power ring relativistic electron beam against the background of immobile ions. By transforming the basic equations, a set of equations for a three-component vortex vector field is derived that describes ring beam configurations for plasma confinement. An example of a numerical calculation of the steady state of a compact beam torus of immobile ions and relativistic electrons is presented.     

Contribution of Na/Ca transport to the resting membrane potential

Relations are derived that describe the combined effects of electrodiffusion, the Na/K pump, and Na/Ca transport by carrier on the resting membrane potential. Equations are derived that apply to both steady-state and non-steady-state conditions. Some example calculations from the equations are plotted at different permeability coefficient ratios, PK:PCa:PNa. The equations predict a depolarizing action of Na/Ca transport when more than two Na ions per Ca ion are transported by the carrier. For all permeability ratios examined, a steady state for Ca ions is achieved with at most a few millivolts of depolarization.     

Periodicities of open linear systems with positive steady states

Systems of two, three, and four linear non-homogeneous differential equations are examined with a view toward determining whether they can possibly serve as mathematical models to describe periodicities in the concentrations of substances which enhance or inhibit each other's rate of production (or dissipation). The nature of the model demands that the solutions of the differential equations be non-negative at all times, i.e., that all the steady states be positive. Conditions for periodicity and for positive steady states are derived, and it is shown that these conditions are not always compatible with each other. In particular it is shown that certain three- and four-hormone models proposed to account for the periodicities observed in the menstrual cycle cannot satisfy the above conditions for any values of the parameters and hence are inadequate.     

A Mathematical Model of the Human Metabolic System and Metabolic Flexibility

In healthy subjects some tissues in the human body display metabolic flexibility, by this we mean the ability for the tissue to switch its fuel source between predominantly carbohydrates in the postprandial state and predominantly fats in the fasted state. Many of the pathways involved with human metabolism are controlled by insulin and insulin-resistant states such as obesity and type-2 diabetes are characterised by a loss or impairment of metabolic flexibility. In this paper we derive a system of 12 first-order coupled differential equations that describe the transport between and storage in different tissues of the human body. We find steady state solutions to these equations and use these results to nondimensionalise the model. We then solve the model numerically to simulate a healthy balanced meal and a high fat meal and we discuss and compare these results. Our numerical results show good agreement with experimental data where we have data available to us and the results show behaviour that agrees with intuition where we currently have no data with which to compare.     

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.     

Bifurcation analysis of nonlinear reaction-diffusion equations—I. Evolution equations and the steady state solutions

A model nonlinear network involving chemical reactions and diffusion is studied. The time evolution and bounds on the steady state solutions are analyzed. Spatially ordered solutions of the equations of the dissipative structure type are found by bifurcation theory. These solutions are calculated analytically and their qualitative properties are discussed.     

Nonequilibrium thermodynamic model of the rat proximal tubule epithelium.

The rat proximal tubule epithelium is represented as well-stirred, compliant cellular and paracellular compartments bounded by mucosal and serosal bathing solutions. With a uniform pCO2 throughout the epithelium, the model variables include the concentrations of Na, K, Cl, HCO3, H2PO4, HPO4, and H, as well as hydrostatic pressure and electrical potential. Except for a metabolically driven Na-K exchanger at the basolateral cell membrane, all membrane transport within the epithelium is passive and is represented by the linear equations of nonequilibrium thermodynamics. In particular, this includes the cotransport of Na-Cl and Na-H2PO4 and countertransport of Na-H at the apical cell membrane. Experimental constraints on the choice of ionic conductivities are satisfied by allowing K-Cl cotransport at the basolateral membrane. The model equations include those for mass balance of the nonreacting species, as well as chemical equilibrium for the acidification reactions. Time-dependent terms are retained to permit the study of transient phenomena. In the steady state the energy dissipation is computed and verified equal to the sum of input from the Na-K exchanger plus the Gibbs free energy of mass addition to the system. The parameter dependence of coupled water transport is studied and shown to be consistent with the predictions of previous analytical models of the lateral intercellular space. Water transport in the presence of an end-proximal (HCO3-depleted) luminal solution is investigated. Here the lower permeability and higher reflection coefficient of HCO3 enhance net sodium and water transport. Due to enhanced flux across the tight junction, this process may permit proximal tubule Na transport to proceed with diminished energy dissipation.     

Ambipolarity condition and possibility of multivalued steady-state solutions to transport equations in stellarators

The possibility of steady-state multivalued solutions to transport equations in stellarators is considered. It is shown that the ambipolarity condition is necessary but not sufficient to find the ambipolar electric field, because the functions entering into it (the plasma density and temperature, as well as their spatial derivatives) depend on the ambipolar field. To do this correctly, it is necessary to solve the full set of time-independent transport equations (including diffusion and heat conduction equations). The possible existence of multivalued solutions to this set of equations is analyzed numerically. It is shown that, under certain conditions that depend on the form and magnitude of particle and heat sources, such solutions can exist. Their form is determined by the initial value of the ambipolar field, the source magnitudes, and the boundary conditions. Discontinuous solutions in which the radial profile of the ambipolar field undergoes jumps are found. In this case, however, the particle and energy fluxes remain continuous, because the discontinuities of the electric field are balanced by the discontinuities of the density and temperature gradients.     

A continuum model for coupled cells

A continuum model of diffusion-coupled cells that more accurately reflects the presence of low-permeability gap junctions between cells is analyzed. It is shown by a multi-scale analysis that to lowest order the slow evolution of the mean concentration is described by the usual ordinary differential equations for a discrete model. Furthermore, stable non-uniform steady solutions are shown to exist in the continuum model of a one component system, whereas this is impossible for the standard reaction-diffusion model of this system. It is also shown how to average the equations in this continuum model to obtain a system of reaction-diffusion equations with constant coefficients.     

Effectiveness factor calculations for immobilized enzyme catalysts

The steady state, nonlinear diffusion equations which describe reactions in constrained enzyme solutions are of great interest in many biological and engineering applications. As in other types of nonlinear differential equations, exact analytical solutions do not exist except in some simplified cases. In this paper, a general procedure is presented for solving numerically for the substrate concentration profile and effectiveness factor utilizing the transformation method suggested by Na and Na. Design correlations for enzyme solutions constrained within spherical membranes are included. The use of a unique definition of the Thiele Modulus in these charts permits the clear illustration of the effects of substrate concentration and external mass transfer resistances on the overall effectiveness factor for the catalyst particle.     

Mixed Convection Boundary Layer Flow over a Moving Vertical Flat Plate in an External Fluid Flow with Viscous Dissipation Effect

The steady boundary layer flow of a viscous and incompressible fluid over a moving vertical flat plate in an external moving fluid with viscous dissipation is theoretically investigated. Using appropriate similarity variables, the governing system of partial differential equations is transformed into a system of ordinary (similarity) differential equations, which is then solved numerically using a Maple software. Results for the skin friction or shear stress coefficient, local Nusselt number, velocity and temperature profiles are presented for different values of the governing parameters. It is found that the set of the similarity equations has unique solutions, dual solutions or no solutions, depending on the values of the mixed convection parameter, the velocity ratio parameter and the Eckert number. The Eckert number significantly affects the surface shear stress as well as the heat transfer rate at the surface.     

The effect of automatic ph maintenance on steady state cultivation regimes and their stability. A theoretical study of the case of a perfect ph controller and perfect stirring

The dynamics of pH-controlled chemostat and pH-auxostat is considered. The medium flow governed by pH-controller is calculated using the ionic (mass-charge) balance equations. The generalized equations of pH-controlled continuous fermentation are developed for the case of two variables: cell biomass and residual substrate concentrations. The study of steady state solutions has shown a substantial dependence of the variable patterns on a concentration of alkali (or acid) in titrant medium flow. The generalized analytical expressions for eigenvalues (reciprocals of transient characteristic times) of the culture dynamics have been obtained. The rates of transients during the steady state establishment are considered using these eigenvalues.     

Numerical bifurcation analysis of the pattern formation in a cell based auxin transport model

Transport models of growth hormones can be used to reproduce the hormone accumulations that occur in plant organs. Mostly, these accumulation patterns are calculated using time step methods, even though only the resulting steady state patterns of the model are of interest. We examine the steady state solutions of the hormone transport model of Smith et al. (Proc Natl Acad Sci USA 103(5):1301–1306, 2006) for a one-dimensional row of plant cells. We search for the steady state solutions as a function of three of the model parameters by using numerical continuation methods and bifurcation analysis. These methods are more adequate for solving steady state problems than time step methods. We discuss a trivial solution where the concentrations of hormones are equal in all cells and examine its stability region. We identify two generic bifurcation scenarios through which the trivial solution loses its stability. The trivial solution becomes either a steady state pattern with regular spaced peaks or a pattern where the concentration is periodic in time.     

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.