A numerical model for the oxygenation of blood in lung capillaries--effect of nth order one-step kinetics of oxygen uptake by haemoglobin

A numerical model is described for the oxygenation of blood in lung capillaries by considering the transport mechanisms of molecular diffusion, convection and the facilitated diffusion due to the presence of haemoglobin. In order to represent the oxygen dissociation curve accurately in the model, the nth order one-step kinetics of oxygen uptake by haemoglobin has been used. The resulting system of coupled, non-linear partial differential equations is solved numerically. It is shown that the blood is required to traverse a larger distance in the capillary before becoming fully oxygenated with nth order one-step kinetics in comparison to first-order one-step kinetics.     

On the Distribution of a Permeable Solute during Poiseuille Flow in Capillary Tubes

Equations are derived describing the dispersion of a permeable solute during Poiseuille flow in a capillary model. It is shown that for the normal range of physiological parameters such as capillary radius, capillary length, blood flow, permeability coefficients, and diffusion constants, the center of mass of a bolus of solute moves at a speed very close to the mean speed of flow and that the solute leaves the capillary with an exponential time course depending on the permeability but not on the diffusion constant. There is no appreciable difference in the dispersion of the solute or in its rate of permeation from the capillary whether one considers piston flow or Poiseuille flow. A bolus of arbitrary radial shape tends to become radially uniform very close to the arterial end of the capillary.     

An electrochemical model of the transport of charged molecules through the capillary glycocalyx

An electrochemical theory of the glycocalyx surface layer on capillary endothelial cells is developed as a model to study the electrochemical dynamics of anionic molecular transport within capillaries. Combining a constitutive relationship for electrochemical transport, derived from Fick's and Ohm's laws, with the conservation of mass and Gauss's law from electrostatics, a system of three nonlinear, coupled, second-order, partial, integro-differential equations is obtained for the concentrations of the diffusing anionic molecules and the cations and anions in the blood. With the exception of small departures from electroneutrality that arise locally near the apical region of the glycocalyx, the model assumes that cations in the blood counterbalance the fixed negative charges bound to the macromolecular matrix of the glycocalyx in equilibrium. In the presence of anionic molecular tracers injected into the capillary lumen, the model predicts the size- and charge-dependent electrophoretic mobility of ions and tracers within the layer. In particular, the model predicts that anionic molecules are excluded from the glycocalyx at equilibrium and that the extent of this exclusion, which increases with increasing tracer and/or glycocalyx electronegativity, is a fundamental determinant of anionic molecular transport through the layer. The model equations were integrated numerically using a Crank-Nicolson finite-difference scheme and Newton-Raphson iteration. When the concentration of the anionic molecular tracer is small compared with the concentration of ions in the blood, a linearized version of the model can be obtained and solved as an eigenvalue problem. The results of the linear and nonlinear models were found to be in good agreement for this physiologically important case. Furthermore, if the fixed-charge density of the glycocalyx is of the order of the concentration of ions in the blood, or larger, or if the magnitude of the anionic molecular valence is large, a closed-form asymptotic solution for the diffusion time can be obtained from the eigenvalue problem that compares favorably with the numerical solution. In either case, if leakage of anionic molecules out of the capillary occurs, diffusion time is seen to vary exponentially with anionic valence and in inverse proportion to the steady-state anionic tracer concentration in the layer relative to the lumen. These findings suggest several methods for obtaining an estimate of the glycocalyx fixed-charge density in vivo.     

血浆层对血管壁剪应力和剪应力梯度的影响

在细小血管中,由于血细胞明显的趋轴效应,管中的血液分为两个不同的区域,即具有血细胞的核心区和邻近管壁和血浆层。应用两相分层流模型,研究在相同的流量和管径下,当核心区中的血液分别为牛顿流体和Casson流体时,不同的血浆层厚度对细小血管壁剪应力和剪应力梯度的影响。结果表明,血浆层的存在对壁剪应力和壁剪应力梯度有较大影响,当血浆层厚度仅为血管半径的1％和3％时,壁剪应力梯度分别下降约10％和20％。     

Transient effects on the initial rate of oxygenation of red blood cells

The rate-controlling process in the oxygenation of red blood cells is investigated using a Roughton-like model for oxygen diffusion and reaction with hemoglobin. The mathematical equations describing the model are solved using two independent techniques, numerical inversions of the Laplace transform of the equations and numerical solutions via an implicit-explicit finite difference form of the equations. The model is used to re-examine previous theoretical models that incorporate either a red cell membrane that is resistive to oxygen diffusion or an unstirred layer of water surrounding the cell. Although both models have been postulated to be equivalent, the results of the computer simulations demonstrate significant differences between the two models in the rate of oxygenation of the red cells, depending upon the values chosen for the diffusion coefficient for O₂ in the membrane and the thickness of the water layer. The difference is apparently due to differences in the induction and transient periods of the water layer model relative to the membrane model.     

Electrokinetic effects on luminal and transmural fluid flow in capillaries

The influence on fluid flow of the fixed charge on the surface of capillaries is calculated using the linearised Poisson-Boltzmann equations. The results depend strongly upon the ratio of the capillary radius to the Debye length. At physiological ionic strength, the Debye length is less than 1 nm and electrostatic effects are negligible. In particular, they can not explain the Copley-Scott Blair phenomenon in artificial capillaries. Electrostatic effects can be significant in smaller channels and it is calculated that in intercellular clefts in the capillary endothelium the apparent viscosity of the fluid may increase more than 50%. These effects can also be important in the flow in the narrow gap between a red cell and the blood capillary wall. Using the Fitzgerald-Lighthill model of this flow and parameters typical of the human microcirculation, the theory predicts that the apparent viscosity in the gap will be increased by about 5%.     

Effect of Dimensional Changes on Plasma Characteristics in Electrothermal Capillary Discharges for Optimized Performance in Fusion Pellet Injection

Geometrical changes in capillary discharges influence the plasma properties and can control exit parameters to certain desired values. For a fixed capillary radius of 2 mm and a 72-μs 43.9-kA peak discharge current, the plasma temperature is about 2.7 eV for different capillary lengths due to the constant input energy, while the number densities tend to saturate for capillary lengths greater than 12 cm. The electrical conductivity reaches 4.02 × 104 Ω^?1 m^?1 and then tends to saturate for 9-cm capillary length. The maximum bulk velocity at the capillary exit slightly increases with the increase in the capillary length from 6.15 to 6.26 km/s for lengths below 18 cm and decreases to 5.88 km/s for longer capillaries due to the higher amount of ablated mass and increased drag forces. For a 9-cm length with the same 72-μs 43.9-kA peak discharge current, the increase in the capillary radius reduces the energy density, which in turn reduces the total ablate mass, plasma density, electrical conductivity, and exit pressure. It is shown that the plasma temperature decreases from 4.6 to 2.1 eV by increasing the capillary radius and radiant heat flux also drops from 463 to 18.1 GW/m². The exit bulk velocity drops from 8.7 to 5.3 km/s as the radius increases from 0.5 to 3.6 mm, respectively. The design features of a capillary discharge can be adjusted for the radius and length, to produce specific plasma parameters for desired applications. Scaling laws relating exit peak plasma parameters to radius and length are obtained to facilitate quick estimate of plasma parameters. The validation of this model has been confirmed by confronting with experimental measurements.     

A semi-empirical model for flow of blood and other particulate suspensions through narrow tubes

A semi-empirical model applicable to the flow of blood and other particulate suspensions through narrow tubes has been developed. It envisages a central core of blood surrounded by a wall layer of reduced hematocrit. With the help of this model the wall layer thickness and extent of plug flow may be calculated using pressure drop, flow rate and hematocrit reduction data. It has been found from the available data in the literature that for a given sample of blood the extent of plug flow increases with decreasing tube diameter. Also for a flow through a given tube it increases with hematocrit. The wall layer thickness is found to decrease with increase in blood hematocrit. A comparison between the results of rigid particulate suspensions and blood reveals that the thicker wall layer and smaller plug flow radius in the case of blood may be attributed to the deformability of the erythrocytes.     

A numerical model for blood oxygenation in the pulmonary capillaries — Effect of pulmonary membrane resistance

A study of the blood oxygenation in pulmonary capillaries is made by considering the transport mechanisms of molecular diffusion, convection and the facilitated diffusion due to the presence of haemoglobin. The resistance offered by the pulmonary membrane on the transport of gases has been incorporated. The resulting system of coupled, non-linear partial differential equations is solved numerically.

It is found that, in the immediate neighbourhood of the entry, the amount of dissolved O₂ decreases. This decreases further as the resistance offered by the pulmonary membrane increases. The rate of oxygenation of blood increases as the permeability coefficient for O₂(P_o) increases. It is shown that the ideally permeable case for both O₂ and CO₂ can be approximated by taking P_o ˜ 10 cm/s. Further, it is shown that the oxygen takes longest and CO₂ is the fastest to attain equilibration. The equilibration length increases as the resistance offered by the membrane increases. Finally, some of the pulmonary diseases such as pulmonary oedema and fibrosis have been analyzed.     

A mathematical model of tumor oxygen and glucose mass transport and metabolism with complex reaction kinetics

Hypoxia imparts radioresistance to tumors, and various approaches have been developed to enhance oxygenation, thereby improving radiosensitivity. This study explores the influence of kinetic and physical factors on substrate metabolism in a tumor model, based on a Krogh cylinder. In tissue, aerobic metabolism is assumed to depend on glucose and oxygen, represented by the product of Michaelis-Menten expressions. For the base case, an inlet pO(2) of 40 mmHg, a hypoxic limit of 5 mmHg, and a tissue/capillary radius ratio of 10 are used. For purely aerobic metabolism, a hypoxic fraction of 0.16 and volume-average pO(2) of 8 mmHg are calculated. Reducing the maximum oxygen rate constant by 9%, decreasing the tissue cylinder radius by 5%, or increasing the capillary radius by 8% abolishes the hypoxic fraction. When a glycolytic term is added, concentration profiles are similar to the base case. Using a distribution of tissue/capillary radius ratios increases the hypoxic fraction and reduces sensitivity to the oxygen consumption rate, compared to the case with a single tissue/capillary radius ratio. This model demonstrates that hypoxia is quite sensitive to metabolic rate and geometric factors. It also predicts quantitatively the effects of inhibited oxygen metabolism and enhanced mass transfer on tumor oxygenation.     

Two-layered blood flow in stenosed tubes for different diseases

A two-fluid model for blood flow through a stenosed tube has been developed. The model consists of a core (suspension of RBCs) and peripheral plasma layer. The core is assumed to be represented by a polar fluid and the plasma layer by a Newtonian fluid. The flow is assumed to be steady and laminar, and the fluids incompressible. The flow variables are computed for normal blood and for the cases of polycythemia, plasma cell dyscrasias and for Hb SS diseases. Resistance to flow has been computed for different stenosis length and for different stenosis height. Shear stress distribution along the axial distance has been computed for different stenosis height. The impact of size effects (particle size to tube diameter) on blood diseases is discussed.     

Some factors affecting oxygen uptake by red blood cells in the pulmonary capillaries

In this study we investigate the equations governing the transport of oxygen in pulmonary capillaries. We use a mathematical model consisting of a red blood cell completely surrounded by plasma within a cylindrical pulmonary capillary. This model takes account of convection and diffusion of oxygen through plasma, diffusion of oxygen through the red blood cell, and the reaction between oxygen and haemoglobin molecules. The velocity field within the plasma is calculated by solving the slow flow equations. We investigate the effect on the solution of the governing equations of: (i) mixed-venous blood oxygen partial pressure (the initial conditions); (ii) alveolar gas oxygen partial pressure (the boundary conditions); (iii) neglecting the convection term; and (iv) assuming an instantaneous reaction between the oxygen and haemoglobin molecules. It is found that: (a) equilibrium is reached much more rapidly for high values of mixed-venous blood and alveolar gas oxygen partial pressure; (b) the convection term has a negligible effect on the time taken to reach a prescribed degree of equilibrium; and (c) an instantaneous reaction may be assumed. Explanations are given for each of these results.     

A Mathematical Study on Three Layered Oscillatory Blood Flow Through Stenosed Arteries

A mathematical model is constructed to examine the characteristics of three layered blood flow through the oscillatory cylindrical tube (stenosed arteries).The proposed model basically consists three layers of blood (viscous fluids with different viscosities) named as core layer (red blood cells),intermediate layer (platelets/white blood cells) and peripheral layer (plasma).The analysis was restricted to propagation of small-amplitude harmonic waves,generated due to blood flow whose wave length is larger compared to the radius of the arterial segment.The impacts of viscosity of fluid in peripheral layer and intermediate layer on the interfaces,average flow rate,mechanical efficiency,trapping and reflux are discussed with the help of numerical and computational results.This model is the generalized form of the preceding models.On the basis of present discussion,it is found that the size of intermediate and peripheral layers reduces in expanded region and enhances in contracted region with the increasing viscosity of fluid in peripheral layer,whereas,opposite effect is observed for viscosity of fluid in intermediate layer.Final conclusion is that the average flow rate and mechanical efficiency increase with the increasing viscosity of fluid in both layers,however,the effects of the viscosity of fluid in both layers on trapping and reflux are opposite to each other.     

Plug Effect of Erythrocytes in Capillary Blood Vessels

As an idealized problem of the motion of blood in small capillary blood vessels, the low Reynolds number flow of plasma (a newtonian fluid) in a circular cylindrical tube involving a series of circular disks is studied. It is assumed in this study that the suspended disks are equally spaced along the axis of the tube, and that their centers remain on the axis of the tube and that their faces are perpendicular to the tube axis. The inertial force of the fluid due to the convective acceleration is neglected on the basis of the smallness of the Reynolds number. The solution of the problem is derived for a quasi-steady flow involving infinitesimally thin disks. The numerical calculation is carried out for a set of different combinations of the interdisk distance and the ratio of the disk radius to the tube radius. The ratio of the velocity of the disk to the average velocity of the fluid is calculated. The different rates of transport of red blood cells and of plasma in capillary blood vessels are discussed. The average pressure gradient along the axis of the tube is computed, and the dependence of the effective viscosity of the blood on the hematocrit and the diameter of the capillary vessel is discussed.     

Mathematical modelling of oxygen transport to tissue

 The equations governing oxygen transport from blood to tissue are presented for a cylindrical tissue compartment, with blood flowing along a co–axial cylindrical capillary inside the tissue. These governing equations take account of: (i) the non–linear reactions between oxygen and haemoglobin in blood and between oxygen and myoglobin in tissue; (ii) diffusion of oxygen in both the axial and radial directions; and (iii) convection of haemoglobin and plasma in the capillary. A non–dimensional analysis is carried out to assess some assumptions made in previous studies. It is predicted that: (i) there is a boundary layer for oxygen partial pressure but not for haemoglobin or myoglobin oxygen saturation close to the inflow boundary in the capillary; (ii) axial diffusion may not be neglected everywhere in the model; (iii) the reaction between oxygen and both haemoglobin and myoglobin may be assumed to be instantaneous in nearly all cases; and (iv) the effect of myoglobin is only significant for tissue with a low oxygen partial pressure. These predictions are validated by solving the full equations numerically and are then interpreted physically. Received: 13 October 2000 / Revised version: 12 June 2001 / Published online: 17 May 2002     

A model for sieving of macromolecules by the glomerular membrane of the kidney.

The transport of water and of macromolecules across the glomerular membrane of the kidney depends on the membrane parameters (radius, length and number of pores) as well as on the hydrostatic and oncotic pressures on either side of the membrane. The filtration pressure decreases along the capillary loops from afferent to efferent end. Water and solute flows are thus given by a system of two differential equations. The sieving coefficient of the macromolecules is the ratio of solute to water flow. In the program described the differential equations are solved by the Runge-Kutta method (fourth order). Rosenbrock's method of minimization is used to adjust the theoretical to the experimental sieving coefficients. The pore radius, total pore area per unit of path length and conductance of the membrane, as well as the intracapillary hydrostatic pressure and its gradient can thus be determined.     

Blood flow in the capillary bed

From arteries to veins, the blood has to go through the ‘capillary’ blood vessels. These blood vessels are so small that often their diameter is smaller than that of the red blood cells. Intimate interactions occur, therefore, between the blood cells and the blood vessels.

A general survey of recent works on capillary blood flow is given in this article. Some details are presented for two problems: the problem of deformation of the flexible red blood cells, their motion in the capillary blood vessels, and the pressure drop due to the red cell blood vessel interaction; and the problem of the flow of plasma ‘bolus’ between neighboring red cells.

The solution supplies many details about the microcirculation phenomenon. Taken together, a method is offered for the calculation of pressure drop in the capillary as a function of various physical parameters: the red cell volume per unit blood volume, (hematocrit), the ratio of the cell diameter to the blood vessel diameter, the ratio of the length of the blood vessel to its length, the volume of individual red cells, and a parameter relating the cell membrane elasticity, plasma viscosity and the cell velocity.     

Three-layered Couette flow of polar fluid with non-zero particle spin boundary condition at the interfaces with applications to blood flow

In this paper, Couette flow of blood is modelled as a three-layered flow. The model basically consists of a core (red-cell suspension) and plasma (a Newtonian fluid) in the top (near the moving plate) and bottom (near the stationary plate) layers. Flow is assumed to be steady and laminar and fluids are incompressible. A spin boundary condition at the interfaces is used by introducing two parameters. Analytic expressions for velocity, total angular velocity and effective viscosity have been obtained and their variations with spin parameters S and s, layer thickness, coupling number N and characteristic length ratio L are computed and shown graphically. One of the important observations of the analysis is the permissible values of the coupling number N is between 0 and 1/square root2 (in the existing literature, the range of N is 0 to 1). The present model includes Couette flow of one and three-layered Newtonian fluids and one-layered polar fluid models as its special cases. Applications of the proposed model to blood flow have been briefly discussed.     

Dispersion of protein solutions in short, open capillary tubes as a method for rapid molecular weight determination

A method has been developed by which the molecular weight of proteins and other freely diffusing species can be estimated on the basis of chromatographic peak shapes developed by injection of a sample into an open capillary tube in a liquid chromatography system. In chromatographic peaks obtained from such a system, there are contributions from both convection and diffusion. Thus, peak shape is dependent upon the diffusion coefficient of the molecular species, the flow rate, and the length of the capillary tube. In the work reported here it has been found that for samples of different proteins ranging from 2000 to 14,000 molecular weight, each injected at the same mobile phase flow rate, the ratio (R) of h1, the height of the peak primarily due to convection, to h2, the height of the "makeup" peak, primarily due to diffusion from the capillary wall, is a direct measure of protein molecular weight. Linear plots of R vs molecular weight are obtained under certain conditions.     

Scaling of brain metabolism and blood flow in relation to capillary and neural scaling

Brain is one of the most energy demanding organs in mammals, and its total metabolic rate scales with brain volume raised to a power of around 5/6. This value is significantly higher than the more common exponent 3/4 relating whole body resting metabolism with body mass and several other physiological variables in animals and plants. This article investigates the reasons for brain allometric distinction on a level of its microvessels. Based on collected empirical data it is found that regional cerebral blood flow CBF across gray matter scales with cortical volume V as CBF ~ V(-1/6), brain capillary diameter increases as V(1/12), and density of capillary length decreases as V(-1/6). It is predicted that velocity of capillary blood is almost invariant (~V(ε)), capillary transit time scales as V(1/6), capillary length increases as V(1/6+ε), and capillary number as V(2/3-ε), where ε is typically a small correction for medium and large brains, due to blood viscosity dependence on capillary radius. It is shown that the amount of capillary length and blood flow per cortical neuron are essentially conserved across mammals. These results indicate that geometry and dynamics of global neuro-vascular coupling have a proportionate character. Moreover, cerebral metabolic, hemodynamic, and microvascular variables scale with allometric exponents that are simple multiples of 1/6, rather than 1/4, which suggests that brain metabolism is more similar to the metabolism of aerobic than resting body. Relation of these findings to brain functional imaging studies involving the link between cerebral metabolism and blood flow is also discussed.