ISPUB.com / IJBE/4/1/10671
  • Author/Editor Login
  • Registration
  • Facebook
  • Google Plus

ISPUB.com

Internet
Scientific
Publications

  • Home
  • Journals
  • Latest Articles
  • Disclaimers
  • Article Submissions
  • Contact
  • Help
  • The Internet Journal of Bioengineering
  • Volume 4
  • Number 1

Original Article

Stochastic Flow as An Inherent Part of Microcirculation

V Kislukhin

Keywords

indicator dilution, markov chain, microcirculation, stochastic flow, vasomotion

Citation

V Kislukhin. Stochastic Flow as An Inherent Part of Microcirculation. The Internet Journal of Bioengineering. 2008 Volume 4 Number 1.

Abstract

There are two phenomena: (a) a bolus of an intravascular tracer after short period of damped oscillation starts to monotonically decrease toward the level of complete mixing; (b) the pass of erythrocytes through microvessels exhibits an irregular flow. The connection between these phenomena is established by mathematical model based on assumptions: (a) hemodynamic stability, (b) the independence of the future trajectory of any blood particles from a past trajectory.Analysis reveals that stochastic behavior of flow is a property of microcirculation, meaning (a) only a fraction of microvessels is open for flow; (b) the state of microvessels, open or closed, is being governed by a random cause. The rate of reassignment of the openness (R) is a characteristic of microcirculation. Low R is the cause for some capillaries be closed for a long time resulting in impairment of nutrient delivery to the surrounding tissue.

 

Introduction

Two phenomena do not appear to be connected: (a) fluctuations of the flow in the microcirculation; (b) a monotonic decrease in concentration of an intravascular tracer (such as erythrocytes-Cr 51 ). Fluctuations can be revealed by observing the red blood cells, and are generated by regular, as well as irregular contractions of smooth muscles, and has been termed as vasomotion [1]. A monotonic decrease in concentration of an intravascular tracer is a well established in the blood volume (BV) measurements [2].

However, mathematical analysis presented in the current manuscript reveals that the only cause for the monotonic decrease is a stochastic flow within the microcirculation. Analysis is based on assumptions: (a) an indicator is in the blood circulation till complete mixing, (b) hemodynamic is stable at time of BV measurement (in other words, the monotonic drop is not due to the leakage of the tracer out of intravascular space and/or the extension of BV during time of measurement) and (c) the independence of any blood particles future trajectory from its past trajectory. Also we assume that the transport of any tracer throughout cardio-vascular system (CVS) is a recurrent process. That is, any particle leaving the left ventricle after passing any of finite chains of artery-arteriole-capillary-venules-vein (systemic as well pulmonary) will return, and times of those various returns will generate a distribution of transit time [3]. Thus, despite the huge number of the paths within vascular system, we assume that the CVS is a finite set of elements. These assumptions lead to the description of the transport of a tracer, as system of differential equations with constant coefficients. It is known that any differential system has as a solution the exponentials series [4] and that is the basis for the connection between stochastic flow and monotonic decrease of a concentration.

Stochastic flow means there is a network of microvessels and: (a) each microvessel from the network is either open for flow or closed; (b) the opening/closing is being governed by a random cause. The probabilities to open and to close are under the central (neuro and hormonal) and local (metabolites) influence.

The recruitment of microvessels, as the mechanism of adaptation to the new level of oxygen demand, has been well established since the research of Krogh [5]. Stochasticity of flow reveals opportunity to change delivery of nutrients without recruitment of microvessels. By changing only reassignment of the openness and keeping fraction of open microvessel constant one could significantly change delivery of nutrients [67]. Thus a stochasticity of flow could be a powerful tool for an adjustment of perfusion to the demand of tissue and the aim of the paper is to reveal that stochastic flow is the inherent property of microcirculation.

Mathematical model for the passage of an intravascular tracer

We start with assumptions: (a) the stability of hemodynamic, and (b) the future trajectory of any particle depends only on the current place of the particle. A model based on the given assumptions is a Markov chain [4] if it includes the following three components: (a) a structure of the CVS, (b) a distribution of a tracer throughout the CVS, and (c) an operator of the transition of an indicator throughout the CVS. In detail:

(a) A structure of the CVS is a set of segments, such as the heart chambers, conductive and micro vessels. The segments are numbered {Sk; k=1...N}. Three remarks: (i) S1 is the notation for the right atrium; (ii) subscript m is used for segment Sm corresponding to the aorta, and (iii) SN is the segment corresponding to the vena cava. Numeration of others segments goes under rule: blood flows from the segments with low subscript to the segments with higher subscript.

(b) Distribution of a tracer throughout the CVS is a vector z(t) = {zk(t), k=1,..,N}, where the k th component of z(t) is the fraction of a tracer within the Sk at time t. In particular, the zm(t) is the tracer recorded in aorta. An initial distribution of an indicator, z(0), is z1(0)=1, and for all k>1 zk(0)=0, meaning that a tracer is injected into the right atrium at time t=0,

(c) An operator A= {aij} provides the transition of a tracer during one cardio cycle, where aij is the fraction of a tracer within Si that passes during one cardio-cycle into Sj : as result the distribution of the tracer at time t, z(t), transforms to the distribution at time t+1: z(t+1) = z(t)A, and, recursively: z(t) = z(0)A t . (2.1)

The standard approach to deal with (2.1) is based on the decomposition of matrices [4] and the dilution curve in the aorta, zm(t), is given by (2.2);

Figure 1

where {si} are the roots of the equation Det(sA-E)=0; bm1, {bmi}, and {bmj} are constants and can be obtained from a combination of eigenvectors of A [4].

The right side of (2.2) is the sum of three sequential terms: the constant, the steadily decreasing term, and the damped oscillating term.

A stochastic flow necessary to produce the monotone decrease of a tracer is generated by vasomotion presented on Figure 1. It is a network of microvessels in which every microvessel is either open or closed. Blood passes through open microvessels and is sequestered in closed microvessels. Greek letters, Figure 1, represent the probabilities for microvessels to change their state during one cardio-cycle: α is the probability to be open and remain open, β is the probability to be open and become closed, v is the probability to be closed and remain closed, and µ is the probability to be closed and become open. There are two characteristics of flow: the fraction of open microvessels, no=µ/(µ+β), and the rate of reassignment of the openness (or the rate of switching between open and closed microvessels [8]), R=β+µ. To insert such a scheme into the mathematical model, let introduce two segments of the microcirculation: Sj-1 (the segment of open microvessels), Sj, (the segment of closed microvessels). The exchange between Sj-1 and Sj is as on Figure 1. Additionally, let introduce segment Sj+1, this is the segment receiving blood from microcirculation given by the segments Sj-1 and Sj. This implies the following equalities for the elements of matrix A ajj = v; aj-1j = β. aj-1j+1 = α and ajj+1 = µ.

Figure 2
Figure 1. Microcirculation as stochastic system White circles are for open vessels; dark circles are for closed vessels

Results

We start with the analysis of diagonal elements of matrix A, these are the elements {aii i=1,…,N}. The aii is the probability of a tracer in segment Si to spend next cardio-cycle in the same segment. There are four segments, heart chambers, with non-zero aii, since during cardio-cycle part of blood remains in each chamber. As established in the following statements only the presence of non heart and non zero {aii} leads to the existence of a steadily decreasing term in (2.2).

First statement: If all non-heart aii are zero, then s1=1 is the only real characteristic number, in other words the dilution curve (2.2) contains only damped oscillation around the level of complete mixing. The frequencies of oscillations are:

Figure 3

with T as the mean time to pass the whole CVS.

Second statement: If some of a non-heart ajj is non-zero (ajj=v=1-µ, and aj-1j = β.) then real positive characteristic number, s2, exists, and 1+ µ< s2<1+β+µ. In other words, random flow leads to the monotone decrease in a tracer concentration.

The proof for these statements is given in Appendix.

Discussion

Having proof of the existence of a stochastic flow, two questions remain: (1) is the given proof acceptable, (2) what is the significance of a stochastic flow.

(1) Acceptance of the proof. Validity of conclusions of mathematic model is based on used assumptions. We made three main assumptions: (i) hemodynamic is stable, (ii) a tracer is intravascular, and (iii) the number of the segments of CVS is finite. Let analyze the assumptions. (i) Most elusive is suggestion about stability of hemodynamic. However, for time of 10 min, with some precautions at time of BV measurement, hypothesis of stability is plausible. (ii) The CVS contains enormous number of segments, however, the huge amount of elements is due to the parallelness of the paths from heart to heart. Mathematically it means that the range of A is relatively low (about 50 at most), and there are only few singularities of the function F(s), (A.3). (iii) Any real indicator is leaking from vascular system. However, erythrocytes 51 Cr remains in vascular system for days and could be considered as strictly intravascular, even albumin 131 I with possible leakage in hours is considered as intravascular if correction for the leakage is done [9]

(2) Significance of stochastic flow. Let assume that stochastic flow is the inherent part of microcirculation. (a) Since the recruitment of microvessels is the main mechanism for increasing blood flow through any tissue [510], then the variation of flow is due to variation of the fraction of open capillaries, no (no=µ/(µ+β)). This can be done by two different ways: (i) by variation of probability for close microvessels become open (it is µ), (ii) by variation of probability of open microvessels become close (it is β). Both probabilities are under neural control, but metabolic influences are different. Open microvessels are a target of delivered metabolites and hormones; close microvessels are a target of local metabolites, thus the regulation of flow could be a complicated process. (b) The tissue around closed capillaries could experience a lack of delivery/removal of substrates/products of metabolism, if the R (rate of vasomotion) is low. The increase of the R could be a tool to improve microcirculation having low fraction of open microvessels. It has been shown [6] that without changes in blood flow through the tissue, just changes in R, one can change oxygen consumption 3-8 fold.

Thus stochastic flow could be a powerful tool in the regulation of a quality of microcirculation.

Conclusion

There are regions of microcirculation with stochastic flow as an immanent property.

Appendix

In sequel aii of the heart chambers will be denoted as aj, j=1, 2, 3, 4.

Statement 1. If non-heart aii=0, then the equation Det(sA-E)=0 can be written as

Figure 4

where aj stays for the residual fraction of the j-heart chamber and the factor ali1 · K · a (ikN) is the probability to pass through the vascular system in (k+1) cardio-cycles.

Equation (A.1) can be rewritten as the equation for generating function [4] to pass whole CVS:

Figure 5

where H(s) and VS(s) are the generating functions for distribution of transit time to pass heart and vascular system, respectively. Since VS(s) is a polynomial with positive coefficients, the direct calculations show that minimum of F(s) outside singular points of H(s) is higher then 1 under condition that the heart volume is less then the volume of vascular system (formally H’(1) < VS’(1)). That means the only damped oscillations around level of complete mixing could be observed. The estimation of frequencies follows from a Taylor decomposition of lnF(s):

Figure 6

Statement 2. If some of the non-heart ajj is non-zero then (A.2) becomes:

Figure 7

Pj1(s) and Pj2(s) are polynomials. The decomposion of (A.3) in vicinity of 1/v reveals the presence of a real characteristic number:

Figure 8

where Tm is the mean time for the tracer to pass the microcirculation, and Topen is the mean time for the tracer to pass the heart, conductive vessels, and open for circulation microvessels.

References

1. Guyton A.C. Textbook of Medical Physiology. Philadelphia, W.B. Saunders Company, 1986
2. Shoemaker WC Measurement of rapidly and slowly circulating red cell volumes in hemorrhagic shock. Am J Physiol. 1962.; 202:1179-82.
3. Meier P, Zierler KL On the theory of the indicator-dilution method for measurement of blood flow and volume. Appl Physiol. 1954; 6(12):731-44.
4. Feller W An Introduction to Probability Theory and Its Applications. New York, London, Sydney, Toronto,John Wiley & Sons, Inc, 1959
5. Krogh A The Anatomy and Physiology of Capillaries, New York, Hafner Publishing CO, , pp. 270-290; 1959
6. Kislukhin V Regulation of oxygen consumption by vasomotion. .Math Biosci. 2004; 191(1): 101-8,.
7. Ranjan K. Pradhan, Chakravarthy V.S., Anil Prabhakar. Effect of chaotic vasomotion in skeletal muscle on tissue oxygenation. Microvasc Res. 2007; 74:51-64,
8. Kislukhin V. Vasomotion model explanation for urea rebound ASAIO J. 2002; 48(3):296-9.
9. Lawson, H.C. The volume of blood - a critical examination of methods for its measurement. In: The Handbook of Physiology: Section 2, Circulation , Volume I, edtr, W.F.Hamilton and P.Dow. Baltimore: Waverly Press; 1962, p.23-49
10. Zierler K. Indicator dilution methods for measuring blood flow, volume, and other properties of biological systems: a brief history and memoir Ann Biomed Eng. 28(8): 836-48, 2000

Author Information

Victor V. Kislukhin
R&D Scientist, Transonic System

Download PDF

Your free access to ISPUB is funded by the following advertisements:

 

 

BACK TO TOP
  • Facebook
  • Google Plus

© 2013 Internet Scientific Publications, LLC. All rights reserved.    UBM Medica Network Privacy Policy

Close

Enter the site

Login

Password

Remember me

Forgot password?

Login

SIGN IN AS A USER

Use your account on the social network Facebook, to create a profile on BusinessPress