Development of a general method for designing microvascular networks using distribution of wall shear stress

In the present study, theoretical formulations for calculation of optimal bifurcation angle and relationship between the diameters of mother and daughter vessels using the power law model for non-Newtonian fluids are developed. The method is based on the distribution of wall shear stress in the mother and daughter vessels. Also, the effect of distribution of wall shear stress on the minimization of energy loss and flow resistance is considered. It is shown that constant wall shear stress in the mother and daughter vessels provides the minimum flow resistance and energy loss of biological flows. Moreover, the effects of different wall shear stresses in the mother and daughter branches, different lengths of daughter branches in the asymmetric bifurcations and non-Newtonian effect of biological fluid flows on the bifurcation angle and the relationship between the diameters of mother and daughter branches are considered. Using numerical simulations for non-Newtonian models such as power law and Carreau models, the effects of optimal bifurcation angle on the pressure drop and flow resistance of blood flow in the symmetric bifurcation are investigated. Numerical simulations show that optimal bifurcation angle decreases the pressure drop and flow resistance especially for bifurcations at large Reynolds number.     

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

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

The blood flow of the brain

The rigidity of the skull and the inertial characteristics and incompressibility of its contents cause the elastic cerebral arteries and veins to act over brief periods of time like rigid tubes of relatively small diameter. Poiseuille's law is applicable to their behavior. The use of this law, in combination with the fact that, during brief intervals, the total volume of the cerebral arteries and veins remains constant, permits derivation of a mathematical expression for the average arterial flow in terms of an average arterial radius. The differentiated equation has five positive roots which represent maxima and minima of the average flow in terms of the average arterial radius. The theoretical results have physiological implications and potential clinical usefulness, which are discussed.     

On connecting large vessels to small. The meaning of Murray's law

A large part of the branching vasculature of the mammalian circulatory and respiratory systems obeys Murray's law, which states that the cube of the radius of a parent vessel equals the sum of the cubes of the radii of the daughters. Where this law is obeyed, a functional relationship exists between vessel radius and volumetric flow, average linear velocity of flow, velocity profile, vessel-wall shear stress, Reynolds number, and pressure gradient in individual vessels. In homogeneous, full-flow sets of vessels, a relation is also established between vessel radius and the conductance, resistance, and cross- sectional area of a full-flow set.     

Patient-specific initial wall stress in abdominal aortic aneurysms with a backward incremental method

Patient-specific wall stress simulations on abdominal aortic aneurysms may provide a better criterion for surgical intervention than the currently used maximum transverse diameter. In these simulations, it is common practice to compute the peak wall stress by applying the full systolic pressure directly on the aneurysm geometry as it appears in medical images. Since this approach does not account for the fact that the measured geometry is already experiencing a substantial load, it may lead to an incorrect systolic aneurysm shape. We have developed an approach to compute the wall stress on the true diastolic geometry at a given pressure with a backward incremental method. The method has been evaluated with a neo-Hookean material law for several simple test problems. The results show that the method can predict an unloaded configuration if the loaded geometry and the load applied are known. The effect of incorporating the initial diastolic stress has been assessed by using three patient-specific geometries acquired with cardiac triggered MR. The comparison shows that the commonly used approach leads to an unrealistically smooth systolic geometry and therefore provides an underestimation for the peak wall stress. Our backward incremental modelling approach overcomes these issues and provides a more plausible estimate for the systolic aneurysm volume and a significantly different estimate for the peak wall stress. When the approach is applied with a more complex material law which has been proposed specifically for abdominal aortic aneurysm similar effects are observed and the same conclusion can be drawn.     

A numerical assessment of wall shear stress changes after endovascular stenting

This theoretical/numerical study aims at assessing the haemodynamic changes induced by endovascular stenting. By using the classical one-dimensional linear pressure waves theory in elastic vessels, we first show that the modulus of the reflection coefficient induced by an endovascular prosthesis is most likely small since it is proportional to the stent-to-wavelength ratio. As a direct consequence, the wall motion of the elastic (stented) artery can be prescribed a priori and the coupled fluid-structure problem does not have to be solved for assessing the haemodynamic changes due to stenting. Several 2D axisymetric calculations are performed to solve the unsteady incompressible Navier-Stokes equations on moving meshes for different types of (stented) arteries. The numerical results suggest that endovascular stenting increases the systo-diastolic variations of the wall shear stress (by 35% at the middle of the stent, by almost 50% in the proximal transition region). Additional calculations show that over-dilated stents produce less haemodynamic perturbations. Indeed, the increase of the amplitude of the wall shear stress variations over the cardiac cycle is only 10% when the stent radius is equal to the radius of the elastic artery at systole (instead of being equal to the mean artery radius).     

An investigation of pulmonary surfactant physicochemical behavior under airway reopening conditions.

Airway reopening mechanics depend on surfactant physicochemical properties. During reopening, the progression of a finger of air down an airway creates an interface that is continually expanding into the bulk fluid. Conventional surfactometers are not capable of evaluating physicochemical behavior under these conditions. To study these aspects, we investigated the pressure required to push a semi-infinite bubble of air down a fluid-filled cylindrical capillary of radius R. The ionic surfactant SDS and pulmonary surfactant analogs L-alpha-dipalmitoylphosphatidylcholine and Infasurf were investigated. We found that the nonequilibrium adsorption of surfactant can create a large nonequilibrium normal stress and a surface shear stress (Marangoni stress) that increase the bubble pressure. The nonphysiological surfactant SDS is capable of eliminating the normal stress and partially reducing the Marangoni stress. The main component of pulmonary surfactant, L-alpha-dipalmitoylphosphatidylcholine, is not capable of reducing either stress, demonstrating slow adsorption properties. The clinically relevant surfactant Infasurf is shown to have intermediate adsorption properties, such that the nonequilibrium normal stress is reduced but the Marangoni stress remains large. Infasurf's behavior suggests that an optimal surfactant solution will have sorption properties that are fast enough to reduce the reopening pressure that may damage airway wall epithelial cells but slow enough to maintain the Marangoni stress that enhances airway stability.     

Heterogeneous mechanics of the mouse pulmonary arterial network

Individualized modeling and simulation of blood flow mechanics find applications in both animal research and patient care. Individual animal or patient models for blood vessel mechanics are based on combining measured vascular geometry with a fluid structure model coupling formulations describing dynamics of the fluid and mechanics of the wall. For example, one-dimensional fluid flow modeling requires a constitutive law relating vessel cross-sectional deformation to pressure in the lumen. To investigate means of identifying appropriate constitutive relationships, an automated segmentation algorithm was applied to micro-computerized tomography images from a mouse lung obtained at four different static pressures to identify the static pressure–radius relationship for four generations of vessels in the pulmonary arterial network. A shape-fitting function was parameterized for each vessel in the network to characterize the nonlinear and heterogeneous nature of vessel distensibility in the pulmonary arteries. These data on morphometric and mechanical properties were used to simulate pressure and flow velocity propagation in the network using one-dimensional representations of fluid and vessel wall mechanics. Moreover, wave intensity analysis was used to study effects of wall mechanics on generation and propagation of pressure wave reflections. Simulations were conducted to investigate the role of linear versus nonlinear formulations of wall elasticity and homogeneous versus heterogeneous treatments of vessel wall properties. Accounting for heterogeneity, by parameterizing the pressure/distention equation of state individually for each vessel segment, was found to have little effect on the predicted pressure profiles and wave propagation compared to a homogeneous parameterization based on average behavior. However, substantially different results were obtained using a linear elastic thin-shell model than were obtained using a nonlinear model that has a more physiologically realistic pressure versus radius relationship.     

Evolution of acoustically vaporized microdroplets in gas embolotherapy

Acoustic vaporization dynamics of a superheated dodecafluoropentane (DDFP) microdroplet inside a microtube and the resulting bubble evolution is investigated in the present work. This work is motivated by a developmental gas embolotherapy technique that is intended to treat cancers by infarcting tumors using gas bubbles. A combined theoretical and computational approach is utilized and compared with the experiments to understand the evolution process and to estimate the resulting stress distribution associated with vaporization event. The transient bubble growth is first studied by ultra-high speed imaging and then theoretical and computational modeling is used to predict the entire bubble evolution process. The evolution process consists of three regimes: an initial linear rapid spherical growth followed by a linear compressed oval shaped growth and finally a slow asymptotic nonlinear spherical bubble growth. Although the droplets are small compared to the tube diameter, the bubble evolution is influenced by the tube wall. The final bubble radius is found to scale linearly with the initial droplet radius and is approximately five times the initial droplet radius. A short pressure pulse with amplitude almost twice as that of ambient conditions is observed. The width of this pressure pulse increases with increasing droplet size whereas the amplitude is weakly dependent. Although the rise in shear stress along the tube wall is found to be under peak physiological limits, the shear stress amplitude is found to be more prominently influenced by the initial droplet size. The role of viscous dissipation along the tube wall and ambient bulk fluid pressure is found to be significant in bubble evolution dynamics.     

The Rheology of Liquid Pentane Isomers by Non-Equilibrium Molecular Dynamics Simulations

Abstract

In this paper, non-equilibrium molecular dynamics (NEMD) simulations of planar Couette flow are reported for an expanded collapsed atom model for liquid pentane isomers at 273.15 K. The strain rate dependent viscosity for liquid pentane isomers exhibits shear-thinning and a linear dependence on γ^1/2. Newtonian viscosities for liquid pentane isomers obtained by a linear extrapolation to zero strain rate are: 0.256cP for normal pentane, 0.219cP for isopentane, and 0.168cP for neopentane. The strain rate dependent pressure difference and normal stress difference vary nearly linearly with the γ^3/2 law and the γ law, respectively, for all three liquid pentane isomers. The overall trend of the square of radius of gyration and end-to-end distance for normal pentane is a linear increase with strain rate. For isopentane, the trend hardly changes for the range of shear rate in this study. The alignment angle decreases with increasing strain rate and the alignment angle of the straight chain alkane is less than that of the branched chain alkane. The average percentage of C?C?C?C trans for normal pentane as a function of strain rate is in excellent correlation with the square of the radius of gyration and the average end-to-end distance. Applying the strain rate in the x-direction, the alignment angle is forced to decrease and the percentage of C?C?C?C trans increases with increasing strain rate.     

Numerical simulation of sinusoidal flow in a straight elastic tube: effects of phase angles

Numerical simulations of flow in straight elastic (moving wall) tubes subjected to a sinusoidal pressure gradient were performed for conditions prevailing in large and medium sized arteries. The effects of varying the phase angle between the pressure gradient and the tube radius, the amplitude of wall motion, and the unsteadiness parameter (alpha) on flow rate and wall shear stress were investigated. Mean and peak flow rates and shear stresses were found to be strongly affected by the phase angle between the pressure gradient and the tube radius with greater sensitivity at higher diameter variation and higher alpha. In large artery simulations (alpha = 12), means flow rate was found to be 60% higher and peak flow rate to be 73% higher than corresponding rigid tube values for certain phase angles, while a threefold increase in mean wall shear stress and sevenfold increase in peak wall shear stress were observed in a sensitive phase angle range. Significant reversal in the wall shear stress direction occurred in the sensitive phase angle range even when there was negligible flow rate reversal. All effects were greatly diminished in simulations of medium sized vessels (alpha = 4). Some experimental evidence to support the predictions of a strong effect of phase angle on wall shear stress in large vessels is presented. Finally, physiological implications of the present work are discussed from a basis of aortic input impedance data, and a physical explanation for the extreme sensitivity of the flow field to small amplitude wall motion at high alpha is given.     

Adaptation of passive rat left ventricle in diastolic dysfunction.

This article deals with providing a theoretical explanation for quantitative changes in the geometry, the opening angle and the deformation parameters of the rat ventricular wall during adaptation of the passive left ventricle in diastolic dysfunction. A large deformation theory is applied to analyse transmural stress and strain distribution in the left ventricular wall considering it to be made of homogeneous, incompressible, transversely isotropic, non-linear elastic material. The basic assumptions made for computing stress distributions are that the average circumferential stress and strain for the adaptive ventricle is equal to the average circumferential stress and strain in the normotensive ventricle, respectively.All the relevant parameters, such as opening angle, twist per unit length, axial extension, internal and external radii and others, in the stress-free, unloaded and loaded states of normotensive, hypertensive and adaptive left ventricle are determined. The circumferential stress and strain distribution through the ventricular wall are also computed. Our analysis predicts that during adaptation, wall thickness and wall mass of the ventricle increase. These results are consistent with experimental findings and are the indications of initiation of congestive heart failure.     

Long time evolution of atherosclerotic plaques

The evolution of atherosclerosis in general, and the influence of wall shear stress on the growth of atherosclerotic plaques in particular, is an intricate phenomenon which is still only partly understood. We therefore propose a qualitative mathematical model which consists of a number of ordinary differential equations for the concentrations of the most relevant constituents of the atherosclerotic plaque. These equations were studied both for the case that the wall shear stress is a parameter (model A), and for the case in which the plaque evolution is coupled to the blood flow (model B) which results in a time dependent wall shear stress. We find that both models exhibit a class of marginally stable equilibria, all reflecting states in which the plaque only grows for a short period of time after a perturbation. The uncoupled model A, however, shows bi-stability between this class of equilibria and another equilibrium state in which the plaque experiences unlimited growth in time, if the LDL cholesterol intake exceeds a threshold value. In model B the bi-stability vanishes, but we find that there is still a critical value of the LDL cholesterol intake beyond which the lumen radius drastically decreases. We show that this decrease is quite sensitive to the value of the wall shear stress.     

Left ventricular wall stress compendium

Left ventricular (LV) wall stress has intrigued scientists and cardiologists since the time of Lame and Laplace in 1800s. The left ventricle is an intriguing organ structure, whose intrinsic design enables it to fill and contract. The development of wall stress is intriguing to cardiologists and biomedical engineers. The role of left ventricle wall stress in cardiac perfusion and pumping as well as in cardiac pathophysiology is a relatively unexplored phenomenon. But even for us to assess this role, we first need accurate determination of in vivo wall stress. However, at this point, 150 years after Lame estimated left ventricle wall stress using the elasticity theory, we are still in the exploratory stage of (i) developing left ventricle models that properly represent left ventricle anatomy and physiology and (ii) obtaining data on left ventricle dynamics. In this paper, we are responding to the need for a comprehensive survey of left ventricle wall stress models, their mechanics, stress computation and results. We have provided herein a compendium of major type of wall stress models: thin-wall models based on the Laplace law, thick-wall shell models, elasticity theory model, thick-wall large deformation models and finite element models. We have compared the mean stress values of these models as well as the variation of stress across the wall. All of the thin-wall and thick-wall shell models are based on idealised ellipsoidal and spherical geometries. However, the elasticity model's shape can vary through the cycle, to simulate the more ellipsoidal shape of the left ventricle in the systolic phase. The finite element models have more representative geometries, but are generally based on animal data, which limits their medical relevance. This paper can enable readers to obtain a comprehensive perspective of left ventricle wall stress models, of how to employ them to determine wall stresses, and be cognizant of the assumptions involved in the use of specific models.     

Peristaltic biofluids flow through vertical porous human vessels using third-grade non-Newtonian fluids model

In this paper, the heat and flow characteristic of third-grade non-Newtonian biofluids flow through a vertical porous human vessel due to peristaltic wall motion are studied. The third-grade model can describe shear thinning (or shear thickening) and normal stress differences, which is acceptable for biofluids modeling. In order to solve the governing equations, the assumption of long-wavelength approximation is utilized. This hypothesis emphasizes that the wavelength of the peristaltic wall motion is large in comparison with the radius of the human vessel, which is widely acceptable in biological investigations. The analytical perturbation method is employed to solve the governing equations. Consequently, analytical expressions for the velocity profile, shear stress, temperature field, and biofluid flow rate are obtained. In addition, the effects of the governing parameters such as the third-grade non-Newtonian parameter, Grashof Number, Eckert number, and porosity, on the results are examined.     

A framework for incorporating 3D hyperelastic vascular wall models in 1D blood flow simulations

We present a novel framework for investigating the role of vascular structure on arterial haemodynamics in large vessels, with a special focus on the human common carotid artery (CCA). The analysis is carried out by adopting a three-dimensional (3D) derived, fibre-reinforced, hyperelastic structural model, which is coupled with an axisymmetric, reduced order model describing blood flow. The vessel transmural pressure and lumen area are related via a Holzapfel–Ogden type of law, and the residual stresses along the thickness and length of the vessel are also accounted for. After a structural characterization of the adopted hyperelastic model, we investigate the link underlying the vascular wall response and blood-flow dynamics by comparing the proposed framework results against a popular tube law. The comparison shows that the behaviour of the model can be captured by the simpler linear surrogate only if a representative value of compliance is applied. Sobol’s multi-variable sensitivity analysis is then carried out in order to identify the extent to which the structural parameters have an impact on the CCA haemodynamics. In this case, the local pulse wave velocity (PWV) is used as index for representing the arterial transmission capacity of blood pressure waveforms. The sensitivity analysis suggests that some geometrical factors, such as the stress-free inner radius and opening angle, play a major role on the system’s haemodynamics. Subsequently, we quantified the differences in haemodynamic variables obtained from different virtual CCAs, tube laws and flow conditions. Although each artery presents a distinct vascular response, the differences obtained across different flow regimes are not significant. As expected, the linear tube law is unable to accurately capture all the haemodynamic features characterizing the current model. The findings from the sensitivity analysis are further confirmed by investigating the axial stretching effect on the CCA fluid dynamics. This factor does not seem to alter the pressure and flow waveforms. On the contrary, it is shown that, for an axially stretched vessel, the vascular wall exhibits an attenuation in absolute distension and an increase in circumferential stress, corroborating the findings of previous studies. This analysis shows that the new model offers a good balance between computational complexity and physics captured, making it an ideal framework for studies aiming to investigate the profound link between vascular mechanobiology and blood flow.

     

Accurate prediction of wall shear stress in a stented artery: newtonian versus non-newtonian models

A significant amount of evidence linking wall shear stress to neointimal hyperplasia has been reported in the literature. As a result, numerical and experimental models have been created to study the influence of stent design on wall shear stress. Traditionally, blood has been assumed to behave as a Newtonian fluid, but recently that assumption has been challenged. The use of a linear model; however, can reduce computational cost, and allow the use of Newtonian fluids (e.g., glycerine and water) instead of a blood analog fluid in an experimental setup. Therefore, it is of interest whether a linear model can be used to accurately predict the wall shear stress caused by a non-Newtonian fluid such as blood within a stented arterial segment. The present work compares the resulting wall shear stress obtained using two linear and one nonlinear model under the same flow waveform. All numerical models are fully three-dimensional, transient, and incorporate a realistic stent geometry. It is shown that traditional linear models (based on blood's lowest viscosity limit, 3.5 Pa s) underestimate the wall shear stress within a stented arterial segment, which can lead to an overestimation of the risk of restenosis. The second linear model, which uses a characteristic viscosity (based on an average strain rate, 4.7 Pa s), results in higher wall shear stress levels, but which are still substantially below those of the nonlinear model. It is therefore shown that nonlinear models result in more accurate predictions of wall shear stress within a stented arterial segment.     

Peristaltic flow in tubes

The study is concerned with the analysis of two flow domains of peristaltic motion in tubes. In the first analysis the wall disturbance wavelength is much larger than the average tube radius. There is a simple algebraic relation between the average flow rate and pressure differential across a wavelength. In the second analysis the disturbance wavelength may be as small as the average radius. A numerical technique may be used to determine the relation between average flow rate and pressure differential across a wavelength.     

A mechanism for arteriolar remodeling based on maintenance of smooth muscle cell activation

Structural adaptation in arterioles is part of normal vascular physiology but is also seen in disease states such as hypertension. Smooth muscle cell (SMC) activation has been shown to be central to microvascular remodeling. We hypothesize that, in a remodeling process driven by SMC activation, stress sensitivity of the vascular wall is a key element in the process of achieving a stable vascular structure. We address whether the adaptive changes in arterioles under different conditions can arise through a common mechanism: remodeling in a stress-sensitive wall driven by a shift in SMC activation. We present a simple dynamic model and show that structural remodeling of the vessel radius by rearrangement of the wall material around a lumen of a different diameter and driven by differences in SMC activation can lead to vascular structures similar to those observed experimentally under various conditions. The change in structure simultaneously leads to uniform levels of circumferential wall stress and wall strain, despite differences in transmural pressure. A simulated vasoconstriction caused by increased SMC activation leads to inward remodeling, whereas outward remodeling follows relaxation of the vascular wall. The results are independent of the specific myogenic properties of the vessel. The simulated results are robust in the face of parameter changes and, hence, may be generalized to vessels from different vascular beds.     

Stress-modulated growth, residual stress, and vascular heterogeneity.

A simple phenomenological model is used to study interrelations between material properties, growth-induced residual stresses, and opening angles in arteries. The artery is assumed to be a thick-walled tube composed of an orthotropic pseudoelastic material. In addition, the normal mature vessel is assumed to have uniform circumferential wall stress, which is achieved here via a mechanical growth law. Residual stresses are computed for three configurations: the unloaded intact artery, the artery after a single transmural cut, and the inner and outer rings of the artery created by combined radial and circumferential cuts. The results show that the magnitudes of the opening angles depend strongly on the heterogeneity of the material properties of the vessel wall and that multiple radial and circumferential cuts may be needed to relieve all residual stress. In addition, comparing computed opening angles with published experimental data for the bovine carotid artery suggests that the material properties change continuously across the vessel wall and that stress, not strain, correlates well with growth in arteries.