Analysis of viscoelastic wall properties in ovine arteries

Constituent-based quasi-linear viscoelasticity: a revised quasi-linear modelling framework to capture nonlinear viscoelasticity in arteries

Arteries exhibit fully nonlinear viscoelastic behaviours (i.e. both elastically and viscously nonlinear). While elastically nonlinear arterial models are well established, effective mathematical descriptions of nonlinear viscoelasticity are lacking. Quasi-linear viscoelasticity (QLV) offers a convenient way to mathematically describe viscoelasticity, but its viscous linearity assumption is unsuitable for whole-wall vascular applications. Conversely, application of fully nonlinear viscoelastic models, involving deformation-dependent viscous parameters, to experimental data is impractical and often reduces to identifying specific solutions for each tested loading condition. The present study aims to address this limitation: By applying QLV theory at the wall constituent rather than at the whole-wall level, the deformation-dependent relative contribution of the constituents allows to capture nonlinear viscoelasticity with a unique set of deformation-independent model parameters. Five murine common carotid arteries were subjected to a protocol of quasi-static and harmonic, pseudo-physiological biaxial loading conditions to characterise their viscoelastic behaviour. The arterial wall was modelled as a constrained mixture of an isotropic elastin matrix and four families of collagen fibres. Constituent-based QLV was implemented by assigning different relaxation functions to collagen- and elastin-borne parts of the wall stress. Nonlinearity in viscoelasticity was assessed via the pressure dependency of the dynamic-to-quasi-static stiffness ratio. The experimentally measured ratio increased with pressure, from 1.03 [Formula: see text] 0.03 (mean [Formula: see text] standard deviation) at 80-40 mmHg to 1.58 [Formula: see text] 0.22 at 160-120 mmHg. Constituent-based QLV captured well this trend by attributing the wall viscosity predominantly to collagen fibres, whose recruitment starts at physiological pressures. In conclusion, constituent-based QLV offers a practical and effective solution to model arterial viscoelasticity.

Viscoelasticity of human descending thoracic aorta in a mock circulatory loop

Healthy human descending thoracic aortas, obtained during organ donation for transplant and research, were tested in a mock circulatory loop to measure the mechanical response to physiological pulsatile pressure and flow. The viscoelastic properties of the aortic segments were investigated at three different pulse rates. The same aortic segments were also subjected to quasi-static pressure tests in order to identify the aortic dynamic stiffness ratio, which is defined as the ratio between the stiffness in case of pulsatile pressure and the stiffness measured for static pressurization, both at the same value of pressure. The loss factor was also identified. The shape of the deformed aorta under static and dynamic pressure was measured by image processing to verify the compatibility of the end supports with the natural deformation of the aorta in the human body. In addition, layer-specific experiments on 10 human descending thoracic aortas allowed to precisely identify the mass density of the aortic tissue, which is an important parameter in cardiovascular dynamic models.

Role of smooth muscle activation in the static and dynamic mechanical characterization of human aortas

The rupture of aortic aneurysms causes around 10,000 deaths each year in the United States. Prosthetic tubes for aortic repair present a large mismatch of mechanical properties with the natural aorta, which has negative consequences for perfusion. This motivates research into the mechanical characterization of human aortas to develop a new generation of mechanically compatible aortic grafts. Experimental data and a suitable material model for human aortas with vascular smooth muscle (VSM) activation are not available. Hence, the present study provides experimental data that are needed. These data made it possible to develop a precise structure-based model of active aortic tissue. The results show the importance of VSM activation on the static and dynamic mechanical response of human aortas.

Experimental data and a suitable material model for human aortas with smooth muscle activation are not available in the literature despite the need for developing advanced grafts; the present study closes this gap. Mechanical characterization of human descending thoracic aortas was performed with and without vascular smooth muscle (VSM) activation. Specimens were taken from 13 heart-beating donors. The aortic segments were cooled in Belzer UW solution during transport and tested within a few hours after explantation. VSM activation was achieved through the use of potassium depolarization and noradrenaline as vasoactive agents. In addition to isometric activation experiments, the quasistatic passive and active stress–strain curves were obtained for circumferential and longitudinal strips of the aortic material. This characterization made it possible to create an original mechanical model of the active aortic material that accurately fits the experimental data. The dynamic mechanical characterization was executed using cyclic strain at different frequencies of physiological interest. An initial prestretch, which corresponded to the physiological conditions, was applied before cyclic loading. Dynamic tests made it possible to identify the differences in the viscoelastic behavior of the passive and active tissue. This work illustrates the importance of VSM activation for the static and dynamic mechanical response of human aortas. Most importantly, this study provides material data and a material model for the development of a future generation of active aortic grafts that mimic natural behavior and help regulate blood pressure.

Viscoelastic characterization of human descending thoracic aortas under cyclic load

Experiments were carried out on 15 human descending thoracic aortas from heart-beating healthy donors who donated organs for transplant. The aortas were kept refrigerated in organ preservation solution and tested were completed within 48 hours from explant. Donors' age was comprised between 25 and 70 years, with an average of 51.7 ± 12.8 years. Quasi-static and dynamic uniaxial tensile test were carried out in thermally controlled physiological saline solution in order to characterize the viscoelastic behavior. Strips were tested under harmonic deformation of different frequency, between 1 and 11 Hz, at three initial pre-stretches. Cyclic deformations of two different amplitudes were used: a physiological one and a small one, the latter one for comparison purposes to understand the accuracy limits of viscoelastic models. Aortic strips in circumferential and longitudinal directions were cut from each aorta. Some strips were dissected to separate the three layers: intima, media and adventitia. They were tested individually in order to obtain layer-specific data. However, strips of the intact wall were also tested. Therefore, 8 strips per donors were tested. Viscoelastic parameters are accurately evaluated from the hysteresis loops. Results show that small-amplitude cyclic strain over-estimate the storage modulus and under-estimate the loss-factor. Therefore, cyclic deformation of physiological amplitude is necessary to obtain correct viscoelastic data of aortic tissue. The value of the applied pre-stretch is significant on the dynamic stiffness ratio (storage modulus divided by the corresponding quasi-static stiffness), while it is less significant for the loss factor. The median of the dynamic stiffness ratios, in physiological conditions, varies between 1.14 and 1.33 for the different layers and the intact wall; the corresponding median of the loss factors varies between 0.050 and 0.066. The lowest dynamic stiffness ratios and loss factors were obtained from donors of the youngest age group. STATEMENT OF SIGNIFICANCE: There is an increasing interest in replacing traditional Dacron grafts used to repair thoracic aortas after acute dissection and aneurysm, with grafts in innovative biomaterials that mimic the mechanical properties and the dynamic behavior of the aorta. The human aorta is a complex laminated structure with hyperelastic and viscoelastic material properties and residual stresses. This study aims to characterize the nonlinear viscoelastic properties of ex-vivo human descending thoracic aortas by measuring hysteresis loops of physiological amplitude under harmonic strain. Results show the necessity to characterize the viscoelastic material properties of the aorta under physiological conditions, as well as the necessity to introduce improved models that take better into account the influence of the initial pre-stretch and amplitude of the cyclic load.

Identifying the parameters of viscoelastic model for a gel-type material as representative of cardiac muscle in dynamic tests

In this paper, mechanical parameters of a calf heart muscle are identified and a gel-type material as the representative of the cardiac muscle in dynamic tests is introduced. The motivation of this study is to introduce a replacement material of the heart muscle to use in experimental studies of the leadless pacemaker. A particular test setup is developed to capture the experimental data based on the stress relaxation test method where its outputs are time histories of the force and displacement. The standard linear solid model is used for mathematical modeling of the heart muscle sample and a gel-type material specimen namely α-gel. Five tests with different strain history (13.6%,17.1%,20.6%22.4%and,23.8%) are performed by regarding and disregarding the influence of the initial ramp of the loading. The mechanical parameters of the standard linear solid model were identified with precise curve fitting. Consideration of the initial ramp significantly influences the consequences and they are so close to their experimental counterparts. The identified parameters of the standard linear solid model by regarding the influence of the initial ramp for the gel-type material are within an acceptable range for the viscoelastic properties of the calf heart tissue. These results show that the gel-type material has the potential to represent the cardiac muscle in the leadless pacemaker experimental studies. Dynamic mechanical analysis is used to characterize the dynamic viscoelastic properties for the gel by utilizing the identified parameters with taking into account the initial ramp in the frequency domain. Results show that Storage modulus, Loss modulus, and Loss tangent are strongly frequency-dependent especially at low-frequency around the heartbeat frequency range (0-2 Hz).

Stress distribution analysis in healthy and stenosed carotid artery models reconstructed from in vivo ultrasonography

Purpose

This study investigated the accuracy of models reconstructed from ultrasound image processing by comparing the radial displacement waveforms of a subject-specific artery model and evaluated stress changes in the proximal shoulder, throat, and distal shoulder of the plaques depending on the degree of carotid artery stenosis.

Methods

Three groups of subjects (healthy and with less than 50% or more carotid stenosis) were evaluated with ultrasonography. Two-dimensional transverse imaging of the common carotid artery was performed to reconstruct the geometry. A longitudinal view of the same region was recorded to extract the Kelvin viscoelastic model parameters. The pulse pressure waveform and the effective pressure of perivascular tissue were loaded onto the internal and external walls of the model. Effective, circumferential, and principal stresses applied to the plaque throat, proximal shoulder, and distal shoulder in the transverse planes were extracted.

Results

The radial displacement waveforms of the model were closely correlated with those of image processing in all three groups. The mean of the effective, circumferential, and principal stresses of the healthy arteries were 15.01±4.93, 12.97±5.07, and 12.39±2.86 kPa, respectively. As stenosis increased from mild to significant, the mean values of the effective, circumferential, and first principal stresses increased significantly (97%, 74%, and 103% at the plaque throat, respectively) (P<0.05). The minimum effective stress was at the lipid pool. The effective stress in calcified areas was higher than in other parts of the artery wall.

Conclusion

This model can discriminate differences in stresses applied to mildly and severely stenotic plaques.

Computational hemodynamics in arteries with the one-dimensional augmented fluid-structure interaction system: viscoelastic parameters estimation and comparison with in-vivo data

Mathematical models are widely recognized as a valuable tool for cardiovascular diagnosis and the study of circulatory diseases, especially to obtain data that require otherwise invasive measurements. To correctly simulate body hemodynamics, the viscoelastic properties of vessels walls are a key aspect to be taken into account as they play an essential role in cardiovascular behavior. The present work aims to apply the augmented fluid-structure interaction system of blood flow to real case studies to assess the validity of the model as a valuable resource to improve cardiovascular diagnostics and the treatment of pathologies. Main contributions of the paper include the evaluation of viscoelastic tube laws, estimation of viscoelastic parameters and comparison of models with literature results and in-vivo experiments. The ability of the model to correctly simulate pulse waveforms in single arterial segments is verified using literature benchmark test cases, designed taking into account a simple elastic behavior of the wall in the upper thoracic aorta and in the common carotid artery. Furthermore, in-vivo pressure waveforms, extracted from tonometric measurements performed on four human common carotid arteries and two common femoral arteries, are compared to numerical solutions. It is highlighted that the viscoelastic damping effect of arterial walls is required to avoid an overestimation of pressure peaks. Finally, an effective procedure to estimate the viscoelastic parameters of the model is herein proposed, which returns hysteresis curves of the common carotid arteries dissipating energy fractions in line with values calculated from literature hysteresis loops in the same vessel.

Numerical modelling of a peripheral arterial stenosis using dimensionally reduced models and kernel methods

In this work, we consider 2 kinds of model reduction techniques to simulate blood flow through the largest systemic arteries, where a stenosis is located in a peripheral artery, i.e., in an artery that is located far away from the heart. For our simulations, we place the stenosis in one of the tibial arteries belonging to the right lower leg (right posterior tibial artery). The model reduction techniques that are used are on the one hand dimensionally reduced models (1-D and 0-D models, the so-called mixed-dimension model) and on the other hand surrogate models produced by kernel methods. Both methods are combined in such a way that the mixed-dimension models yield training data for the surrogate model, where the surrogate model is parametrised by the degree of narrowing of the peripheral stenosis. By means of a well-trained surrogate model, we show that simulation data can be reproduced with a satisfactory accuracy and that parameter optimisation or state estimation problems can be solved in a very efficient way. Furthermore, it is demonstrated that a surrogate model enables us to present after a very short simulation time the impact of a varying degree of stenosis on blood flow, obtaining a speedup of several orders over the full model.

Quantitative Vascular Evaluation: From Laboratory Experiments to Point-of-Care Patient (Experimental Approach)

This paper illustrates the evolution of our knowledge of arterial mechanics from our initial research works up to the present time. Several techniques focusing on this topic in terms of our experience are dis-cussed. An interdisciplinary team composed by different institutions from Argentina, Uruguay, France and Spain was created to conduct research, to train human resources and to fulfill the inevitable social role of gaining access to technological inno-vation to improve public health.

Lumped-parameter models of the pulmonary vasculature during the progression of pulmonary arterial hypertension

A longitudinal study of monocrotaline-induced pulmonary arterial hypertension (PAH) was carried out in Sprague-Dawley rats to investigate the changes in impedance (comprising resistance and compliance) produced by elevated blood pressure. Using invasively measured blood flow as an input, blood pressure was predicted using 3- and 4-element Windkessel (3WK, 4WK) type lumped-parameter models. Resistance, compliance, and inductance model parameters were obtained for the five different treatment groups via least-squares errors. The treated animals reached levels of hypertension, where blood pressure increased two folds from control to chronic stage of PAH (mean pressure went from 24 ± 5 to 44 ± 6 mmHg, P < 0.0001) but blood flow remained overall unaffected. Like blood pressure, the wave-reflection coefficient significantly increased at the advanced stage of PAH (0.26 ± 0.09 to 0.52 ± 0.09, P < 0.0002). Our modeling efforts revealed that resistances and compliance changed during the disease progression, where changes in compliance occur before the changes in resistance. However, resistance and compliance are not directly inversely related. As PAH develops, resistances increase nonlinearly (R_d exponentially and R at a slower rate) while compliance linearly decreases. And while 3WK and 4WK models capture the pressure-flow relation in the pulmonary vasculature during PAH, results from Akaike Information Criterion and sensitivity analysis allow us to conclude that the 3WK is the most robust and accurate model for this system. Ninety-five percent confidence intervals of the predicted model parameters are included for the population studied. This work establishes insight into the complex remodeling process occurring in PAH.

Cardiovascular mechanics in the early stages of pulmonary hypertension: a computational study

We formulate and study a new mathematical model of pulmonary hypertension. Based on principles of fluid and elastic dynamics, we introduce a model that quantifies the stiffening of pulmonary vasculature (arteries and arterioles) to reproduce the hemodynamics of the pulmonary system, including physiologically consistent dependence between compliance and resistance. This pulmonary model is embedded in a closed-loop network of the major vessels in the body, approximated as one-dimensional elastic tubes, and zero-dimensional models for the heart and other organs. Increasingly severe pulmonary hypertension is modeled in the context of two extreme scenarios: (1) no cardiac compensation and (2) compensation to achieve constant cardiac output. Simulations from the computational model are used to estimate cardiac workload, as well as pressure and flow traces at several locations. We also quantify the sensitivity of several diagnostic indicators to the progression of pulmonary arterial stiffening. Simulation results indicate that pulmonary pulse pressure, pulmonary vascular compliance, pulmonary RC time, luminal distensibility of the pulmonary artery, and pulmonary vascular impedance are much better suited to detect the early stages of pulmonary hypertension than mean pulmonary arterial pressure and pulmonary vascular resistance, which are conventionally employed as diagnostic indicators for this disease.

Comparative in vivo analysis of the role of the adventitia and the endothelium on arterial mechanical function: relevance for aortic counterpulsation

PURPOSE

The comparative effect of the intimal and adventitial layers on arterial biomechanics control, in basal and altered conditions, remains to be elucidated. This study aimed (1) to characterize the arterial conduit (CF) and buffering (distensibility) function of the iliac arteries in in vivo animals, in which the intimal and adventitial layers were removed; (2) to determine the effects of intra-aortic ballon pumping (IABP) on simultaneously de-adventitialized (DA) and de-endothelialized (DE) iliac arteries before and after induced heart failure.

METHODS

Pressure and diameter signals were measured in the iliac arteries of sheep (n = 7) in which the adventitial and intima layer were removed. Intra-aortic balloon pump (IABP) assistance was used in a control state and after heart failure induction.

RESULTS

Both DE and DA determined significant changes in arterial diameter, distensibility and CF. Changes were higher after DA than after DE in terms of distensibility and CF (p<0.05). DA followed by DE (DA + DE) showed significant increases in arterial diameter and CF, accompanied by a decrease in distensibility (p<0.05) with respect to intact arteries. Heart failure induction caused significant hemodynamic changes without modifying the already impaired local biomechanical parameters. Nonsignificant improvements in the biomechanical parameters of DA + DE iliac arteries were observed during IABP before and after heart failure induction.

CONCLUSIONS

Biomechanical changes caused by DA of iliac arteries were more important than those observed after DE. The DA + DE arteries showed significant differences with respect to intact arteries and with DA or DE arteries. IABP-related effects on arterial mechanics were absent in DA + DE arteries.

Uncertainty Quantification in a Patient-Specific One-Dimensional Arterial Network Model: EnKF-Based Inflow Estimator

Successful clinical use of patient-specific models for cardiovascular dynamics depends on the reliability of the model output in the presence of input uncertainties. For 1D fluid dynamics models of arterial networks, input uncertainties associated with the model output are related to the specification of vessel and network geometry, parameters within the fluid and wall equations, and parameters used to specify inlet and outlet boundary conditions. This study investigates how uncertainty in the flow profile applied at the inlet boundary of a 1D model affects area and pressure predictions at the center of a single vessel. More specifically, this study develops an iterative scheme based on the ensemble Kalman filter (EnKF) to estimate the temporal inflow profile from a prior distribution of curves. The EnKF-based inflow estimator provides a measure of uncertainty in the size and shape of the estimated inflow, which is propagated through the model to determine the corresponding uncertainty in model predictions of area and pressure. Model predictions are compared to ex vivo area and blood pressure measurements in the ascending aorta, the carotid artery, and the femoral artery of a healthy male Merino sheep. Results discuss dynamics obtained using a linear and a nonlinear viscoelastic wall model.

Dynamic Viscoelasticity and Surface Properties of Porcine Left Anterior Descending Coronary Arteries

The aim of this study was, for the first time, to measure and compare quantitatively the viscoelastic properties and surface roughness of coronary arteries. Porcine left anterior descending coronary arteries were dissected ex vivo. Viscoelastic properties were measured longitudinally using dynamic mechanical analysis, for a range of frequencies from 0.5 to 10 Hz. Surface roughness was calculated following three-dimensional reconstructed of surface images obtained using an optical microscope. Storage modulus ranged from 14.47 to 25.82 MPa, and was found to be frequency-dependent, decreasing as the frequency increased. Storage was greater than the loss modulus, with the latter found to be frequency-independent with a mean value of 2.10 ± 0.33 MPa. The circumferential surface roughness was significantly greater (p < 0.05) than the longitudinal surface roughness, ranging from 0.73 to 2.83 and 0.35 to 0.92 µm, respectively. However, if surface roughness values were corrected for shrinkage during processing, circumferential and longitudinal surface roughness were not significantly different (1.04 ± 0.47, 0.89 ± 0.27 µm, respectively; p > 0.05). No correlation was found between the viscoelastic properties and surface roughness. It is feasible to quantitatively measure the viscoelastic properties of coronary arteries and the roughness of their endothelial surface.

Linear and Nonlinear Viscoelastic Arterial Wall Models: Application on Animals

This work deals with the viscoelasticity of the arterial wall and its influence on the pulse waves. We describe the viscoelasticity by a nonlinear Kelvin–Voigt model in which the coefficients are fitted using experimental time series of pressure and radius measured on a sheep's arterial network. We obtained a good agreement between the results of the nonlinear Kelvin–Voigt model and the experimental measurements. We found that the viscoelastic relaxation time—defined by the ratio between the viscoelastic coefficient and the Young's modulus—is nearly constant throughout the network. Therefore, as it is well known that smaller arteries are stiffer, the viscoelastic coefficient rises when approaching the peripheral sites to compensate the rise of the Young's modulus, resulting in a higher damping effect. We incorporated the fitted viscoelastic coefficients in a nonlinear 1D fluid model to compute the pulse waves in the network. The damping effect of viscoelasticity on the high-frequency waves is clear especially at the peripheral sites.

Modeling the differentiation of A- and C-type baroreceptor firing patterns

The baroreceptor neurons serve as the primary transducers of blood pressure for the autonomic nervous system and are thus critical in enabling the body to respond effectively to changes in blood pressure. These neurons can be separated into two types (A and C) based on the myelination of their axons and their distinct firing patterns elicited in response to specific pressure stimuli. This study has developed a comprehensive model of the afferent baroreceptor discharge built on physiological knowledge of arterial wall mechanics, firing rate responses to controlled pressure stimuli, and ion channel dynamics within the baroreceptor neurons. With this model, we were able to predict firing rates observed in previously published experiments in both A- and C-type neurons. These results were obtained by adjusting model parameters determining the maximal ion-channel conductances. The observed variation in the model parameters are hypothesized to correspond to physiological differences between A- and C-type neurons. In agreement with published experimental observations, our simulations suggest that a twofold lower potassium conductance in C-type neurons is responsible for the observed sustained basal firing, where as a tenfold higher mechanosensitive conductance is responsible for the greater firing rate observed in A-type neurons. A better understanding of the difference between the two neuron types can potentially be used to gain more insight about pathophysiology and treatment of diseases related to baroreflex function, e.g. in patients with autonomic failure, a syndrome that is difficult to diagnose in terms of its pathophysiology.

Comparison of the Windkessel model and structured-tree model applied to prescribe outflow boundary conditions for a one-dimensional arterial tree model

One-dimensional (1D) modeling is a widely adopted approach for studying wave propagation phenomena in the arterial system. Despite the frequent use of the Windkessel (WK) model to prescribe outflow boundary conditions for 1D arterial tree models, it remains unclear to what extent the inherent limitation of the WK model in describing wave propagation in distal vasculatures affect hemodynamic variables simulated at the arterial level. In the present study, a 1D model of the arterial tree was coupled respectively with a WK boundary model and a structured-tree (ST) boundary model, yielding two types of arterial tree models. The effective resistances, compliances and inductances of the WK and ST boundary models were matched to facilitate quantitative comparisons. Obtained results showed that pressure/flow waves simulated by the two models were comparable in the aorta, whereas, their discrepancies increased towards the periphery. Wave analysis revealed that the differences in reflected waves generated by the boundary models were the major sources of pressure wave discrepancies observed in large arteries. Additional simulations performed under aging conditions demonstrated that arterial stiffening with age enlarged the discrepancies, but with the effects being partly counteracted by physiological aortic dilatation with age. These findings suggest that the method adopted for modeling the outflow boundary conditions has considerable influence on the performance of a 1D arterial tree model, with the extent of influence varying with the properties of the arterial system.

WAVE PROPAGATION IN A 1D FLUID DYNAMICS MODEL USING PRESSURE-AREA MEASUREMENTS FROM OVINE ARTERIES

This study considers a 1D fluid dynamics arterial network model with 14 vessels developed to assimilate ex vivo 0D temporal data for pressure-area dynamics in individual vessel segments from 11 male Merino sheep. A 0D model was used to estimate vessel wall parameters in a two-parameter elastic model and a four-parameter Kelvin viscoelastic model. This was done using nonlinear optimization minimizing the least squares error between model predictions and measured cross-sectional areas. Subsequently, estimated values for elastic stiffness and unstressed area were related to construct a nonlinear relationship. This relation was used in the network model. A 1D single vessel model of the aorta was then developed and used to estimate the inflow profile and parameters for total resistance and compliance for the downstream network and to demonstrate effects of incorporating viscoelasticity in the arterial wall. Lastly, the extent to which vessel wall parameters estimated from ex vivo data can be used to realistically simulate pressure and area in a vessel network was evaluated. Elastic wall parameters in the network simulations were found to yield pressure-area relationships across all vessel locations and sheep that were in ranges comparable to those in the ex vivo data.

A viscoelastic fluid-structure interaction model for carotid arteries under pulsatile flow

In this study, a fluid-structure interaction model (FSI) incorporating viscoelastic wall behaviour is developed and applied to an idealized model of the carotid artery under pulsatile flow. The shear and bulk moduli of the arterial wall are described by Prony series, where the parameters can be derived from in vivo measurements. The aim is to develop a fully coupled FSI model that can be applied to realistic arterial geometries with normal or pathological viscoelastic wall behaviour. Comparisons between the numerical and analytical solutions for wall displacements demonstrate that the coupled model is capable of predicting the viscoelastic behaviour of carotid arteries. Comparisons are also made between the solid only and FSI viscoelastic models, and the results suggest that the difference in radial displacement between the two models is negligible.

Scalability and in vivo validation of a multiscale numerical model of the left coronary circulation

Multiscale modeling is a promising tool for the study of coronary hemodynamics. A key strength of this approach is that it accounts for microvascular properties and extravascular forces that differ regionally and transmurally, as well as wave propagation effects in the conduit arteries. However, little validation of such models has been reported and no models of the newborn coronary circulation have been described. We therefore validated a multiscale model of the left coronary circulation using high-fidelity data from nine adult sheep and nine newborn lambs and investigated whether wave propagation effects are more prominent in adults, whose body size (and hence wave transit distance) is greater. The model consisted of a one-dimensional (1D) network of the major conduit arteries and a lumped parameter model of microvascular beds. Intramyocardial pressure was considered to arise via contraction-related myocyte thickening and transmission of ventricular cavity pressure into the heart wall. 1D network geometry from published human anatomical data was scaled using myocardial weights, while subject-specific aortic pressure/flow and ventricular pressure formed model inputs. Total vascular resistance was determined iteratively from measured mean circumflex coronary flow (CxQ), but no fitting of phasic aspects of the waveform was performed. Excellent agreement was obtained between simulated and measured CxQ waveforms in most cases. Detailed flow waveform analysis did not clearly reveal a greater prominence of wave propagation effects in adults compared with newborns. This multiscale model is likely to be useful for investigating wave phenomena and phasic aspects of coronary flow in adults and during development.

A practical approach to parameter estimation applied to model predicting heart rate regulation

Mathematical models have long been used for prediction of dynamics in biological systems. Recently, several efforts have been made to render these models patient specific. One way to do so is to employ techniques to estimate parameters that enable model based prediction of observed quantities. Knowledge of variation in parameters within and between groups of subjects have potential to provide insight into biological function. Often it is not possible to estimate all parameters in a given model, in particular if the model is complex and the data is sparse. However, it may be possible to estimate a subset of model parameters reducing the complexity of the problem. In this study, we compare three methods that allow identification of parameter subsets that can be estimated given a model and a set of data. These methods will be used to estimate patient specific parameters in a model predicting baroreceptor feedback regulation of heart rate during head-up tilt. The three methods include: structured analysis of the correlation matrix, analysis via singular value decomposition followed by QR factorization, and identification of the subspace closest to the one spanned by eigenvectors of the model Hessian. Results showed that all three methods facilitate identification of a parameter subset. The "best" subset was obtained using the structured correlation method, though this method was also the most computationally intensive. Subsets obtained using the other two methods were easier to compute, but analysis revealed that the final subsets contained correlated parameters. In conclusion, to avoid lengthy computations, these three methods may be combined for efficient identification of parameter subsets.

A quasi-lumped model for the peripheral distortion of the arterial pulse

As blood circulates through the arterial tree, the flow and pressure pulse distort. Principal factors to this distortion are reflections form arterial bifurcations and the viscous character of the flow of the blood. Both of them are expounded in the literature and included in our analysis. The nonlinearities of inertial effects are usually taken into account in numerical simulations, based on Navier-Stokes like equations. Nevertheless, there isn't any qualitative, analytical formula, which examines the role of blood's inertia on the distortion of the pulse. We derive such an analytical nonlinear formula. It emanates from a generalized Bernoulli's equation for an an-harmonic, linear, viscoelastic, Maxwell fluid flow in a linear, viscoelastic, Kelvin-Voigt, thin, cylindrical vessel. We report that close to the heart, convection effects related to the change in the magnitude of the velocity of blood dominate the alteration of the shape of the pressure pulse, while at remote sites of the vascular tree, convection of vorticity, related to the change in the direction of the velocity of blood with respect to a mean axial flow, prevails. A quantitative comparison between the an-harmonic theory and related pressure measurements is also performed.

Pulse wave propagation in a model human arterial network: Assessment of 1-D visco-elastic simulations against in vitro measurements

The accuracy of the nonlinear one-dimensional (1-D) equations of pressure and flow wave propagation in Voigt-type visco-elastic arteries was tested against measurements in a well-defined experimental 1:1 replica of the 37 largest conduit arteries in the human systemic circulation. The parameters required by the numerical algorithm were directly measured in the in vitro setup and no data fitting was involved. The inclusion of wall visco-elasticity in the numerical model reduced the underdamped high-frequency oscillations obtained using a purely elastic tube law, especially in peripheral vessels, which was previously reported in this paper [Matthys et al., 2007. Pulse wave propagation in a model human arterial network: Assessment of 1-D numerical simulations against in vitro measurements. J. Biomech. 40, 3476–3486]. In comparison to the purely elastic model, visco-elasticity significantly reduced the average relative root-mean-square errors between numerical and experimental waveforms over the 70 locations measured in the in vitro model: from 3.0% to 2.5% (p<0.012) for pressure and from 15.7% to 10.8% (p<0.002) for the flow rate. In the frequency domain, average relative errors between numerical and experimental amplitudes from the 5th to the 20th harmonic decreased from 0.7% to 0.5% (p<0.107) for pressure and from 7.0% to 3.3% (p<10−6) for the flow rate. These results provide additional support for the use of 1-D reduced modelling to accurately simulate clinically relevant problems at a reasonable computational cost.

Linear and nonlinear viscoelastic modeling of aorta and carotid pressure-area dynamics under in vivo and ex vivo conditions

A better understanding of the biomechanical properties of the arterial wall provides important insight into arterial vascular biology under normal (healthy) and pathological conditions. This insight has potential to improve tracking of disease progression and to aid in vascular graft design and implementation. In this study, we use linear and nonlinear viscoelastic models to predict biomechanical properties of the thoracic descending aorta and the carotid artery under ex vivo and in vivo conditions in ovine and human arteries. Models analyzed include a four-parameter (linear) Kelvin viscoelastic model and two five-parameter nonlinear viscoelastic models (an arctangent and a sigmoid model) that relate changes in arterial blood pressure to the vessel cross-sectional area (via estimation of vessel strain). These models were developed using the framework of Quasilinear Viscoelasticity (QLV) theory and were validated using measurements from the thoracic descending aorta and the carotid artery obtained from human and ovine arteries. In vivo measurements were obtained from 10 ovine aortas and 10 human carotid arteries. Ex vivo measurements (from both locations) were made in 11 male Merino sheep. Biomechanical properties were obtained through constrained estimation of model parameters. To further investigate the parameter estimates, we computed standard errors and confidence intervals and we used analysis of variance to compare results within and between groups. Overall, our results indicate that optimal model selection depends on the artery type. Results showed that for the thoracic descending aorta (under both experimental conditions), the best predictions were obtained with the nonlinear sigmoid model, while under healthy physiological pressure loading the carotid arteries nonlinear stiffening with increasing pressure is negligible, and consequently, the linear (Kelvin) viscoelastic model better describes the pressure-area dynamics in this vessel. Results comparing biomechanical properties show that the Kelvin and sigmoid models were able to predict the zero-pressure vessel radius; that under ex vivo conditions vessels are more rigid, and comparatively, that the carotid artery is stiffer than the thoracic descending aorta; and that the viscoelastic gain and relaxation parameters do not differ significantly between vessels or experimental conditions. In conclusion, our study demonstrates that the proposed models can predict pressure-area dynamics and that model parameters can be extracted for further interpretation of biomechanical properties.

Numerical assessment of time-domain methods for the estimation of local arterial pulse wave speed

A local estimation of pulse wave speed c, an important predictor of cardiovascular events, can be obtained at arterial locations where simultaneous measurements of blood pressure (P) and velocity (U), arterial diameter (D) and U, flow rate (Q) and cross-sectional area (A), or P and D are available, using the PU-loop, sum-of-squares (∑(2)), lnDU-loop, QA-loop or new D(2)P-loop methods. Here, these methods were applied to estimate c from numerically generated P, U, D, Q and A waveforms using a visco-elastic one-dimensional model of the 55 larger human systemic arteries in normal conditions. Theoretical c were calculated from the parameters of the model. Estimates of c given by the loop methods were closer to theoretical values and more uniform within each arterial segment than those obtained using the ∑(2). The smaller differences between estimates and theoretical values were obtained using the D(2)P-loop method, with root-mean-square errors (RMSE) smaller than 0.18 ms(-1), followed by averaging the two c given by the PU- and lnDU-loops (RMSE <2.99 ms(-1)). In general, the errors of the PU-, lnDU- and QA-loops decreased at locations where visco-elastic effects were small and nearby junctions were well-matched for forward-travelling waves. The ∑(2) performed better at proximal locations.

Physiological interpretation of inductance and low-resistance terms in four-element windkessel models: assessment by generalized sensitivity function analysis

Physiological relevance of parameters of three arterial models, denominated W4P, W4S and IVW, was assessed by computation of parameter-related generalized sensitivity functions (GSFs), which allow the definition of heart-cycle time intervals where the information content of experimental data, useful for estimation of each model parameter, is concentrated. The W4P and W4S are derived from the three-element windkessel by connecting an inductance, L, in parallel or in series, respectively, with aortic characteristic impedance, R(c). In the IVW, L is placed in series at the input of a viscoelastic windkessel, incorporating a Voigt cell (a resistor, R(d), in series with a capacitor, C). Pressure and flow measured in the ascending aorta of five ferrets and five dogs were used to estimate all model parameters, by fitting to pressure. For each model structure, parameter-related GSFs were generated. Focusing on controversial L, R(c) and R(d) physical meaning, our GSF analysis yielded the conclusion that, in both the W4S and the IVW, but not in the W4P, the L-term is suitable to represent the inertial properties of blood motion. Moreover, the meaning of aortic characteristic impedance ascribed to R(c) is questionable; while R(d) is likely to account for viscous losses of arterial wall motion.

Linear and nonlinear viscoelastic modeling of ovine aortic biomechanical properties under in vivo and ex vivo conditions

This study uses linear and nonlinear viscoelastic models to describe the dynamic distention of the aorta induced by time-varying arterial blood pressure. We employ an inverse mathematical modeling approach on a four-parameter (linear) Kelvin viscoelastic model and two five-parameter nonlinear viscoelastic models (arctangent and sigmoid) to infer vascular biomechanical properties under in vivo and ex vivo experimental conditions in ten and eleven male Merino sheep, respectively. We used the Akaike Information Criterion (AIC) as a goodness-of-fit measure. Results show that under both experimental conditions, the nonlinear models generally outperform the linear Kelvin model, as judged by the AIC. Furthermore, the sigmoid nonlinear viscoelastic model consistently achieves the lowest AIC and also matches the zero-stress vessel radii measured ex vivo. Based on these observations, we conclude that the sigmoid nonlinear viscoelastic model best describes the biomechanical properties of ovine large arteries under both experimental conditions considered in this study.