科研之友
微信
新浪微博
Facebook
分享链接
摘要
Evidence is presented to show that self-sustained oscillations of purely hemodynamical origin are possible in some arcade-type microvascular networks supplied with steady boundary conditions, but that in others the oscillations disappear with sufficient reduction of the time step Delta t, showing them to be numerical artefacts. In an attempt to elucidate the mechanisms involved in the onset of fluctuations, we proceed to perform a linear stability analysis for the convective model of Kiani et al. (Microvasc. Res. 45:219-232, 1993; Am. J. Physiol. 266(35):H1822-H1828, 1994), and show that this leads via a system of delay differential equations to a nonlinear eigenvalue problem. This result generalises the characteristic equation obtained by Carr et al. (Ann. Biomed. Eng. 33:764-771, 2005) and Geddes et al. (SIAM J. Appl. Dyn. Syst. 6(4):694-727, 2007) who solved a special case in a two node network. An implicit numerical method is proposed for the computation of blood flows in networks using the convective model. In a moderate size subnetwork of one of the networks chosen by Kiani et al. (Am. J. Physiol. 266(35):H1822-H1828, 1994), the topology, vessel lengths, and diameters of which were based on microvascular networks in the rat mesentery, we compare results generated using the original explicit numerical method of Kiani et al. (Am. J. Physiol. 266(35):H1822-H1828, 1994) with those from our implicit scheme. From the linear stability theory, a critical value D (RBC,crit) of a red blood cell diameter parameter D (RBC) in the plasma skimming model of Fenton et al. (Pflugers Arch. 403:396-401, 1985b) is identified for the onset of oscillations about steady state and both the explicit and implicit methods are used to calculate the inflow hematocrit solutions in all vessels of the subnetwork at the critical parameter value, subject to perturbed initial conditions. The results of the implicit method are demonstrated to be in excellent and superior agreement with the predictions of the linear analysis in this case. For values of D (RBC) slightly larger than D (RBC,crit) the bifurcating periodic solutions calculated using either the explicit or implicit schemes are characteristic of those of a supercritical Hopf bifurcation and the graphs of D (RBC) vs. oscillation amplitude would seem to converge as Delta t -> 0.
- 出版日期2013-4
