Laser anemometry measurements of pulsatile flow past aortic valve prostheses

Experimental results are presented on physiological pulsatile flow past caged ball and tilting disc aortic valve prostheses mounted in an axisymmetric chamber incorporated in a mock circulatory system. The measurements of velocity profiles and turbulent normal stresses during several times in a cardiac cycle were obtained using laser-Doppler anemometry. Our results show that with increased angle of opening for the tilting disc valves, a large but locally confined vortex is observed along the wall in the minor flow region throughout most of the cardiac cycle. The turbulent normal stresses measured downstream to the tilting disc in the minor flow region parallel to the tilt axis were found to be larger than those measured downstream to the caged ball valves. Comparison of measurements with steady flow at flow rates comparable to peak pulsatile flow rate show that the turbulent normal stresses are larger by a factor of two in pulsatile flow with a frequency of 1.2 Hz.     

A structural and dynamic molecular model for the sodium channel of Electrophorus electricus

Chemical logic and single group rotation (SGR) theory are applied to the primary structure determined by Noda et al. [(1984) Nature 312, 121-127] to construct a molecular model of the sodium channel of Electrophorus electricus. Both structural and dynamic aspects of the channel are accounted for, including gating current, sensitivity to changes in membrane potential, channel opening, a binding site for sodium, selectivity for sodium over potassium, capacity for rapid sodium flow, sensitivity to batrachotoxin (or other toxins) and inactivation.     

Pulsatile flow past aortic valve bioprostheses in a model human aorta

Pulsatile flow development past tissue valve prostheses in a model human aorta has been studied using qualitative flow visualization and quantitative laser-Doppler techniques. Experiments were conducted both in steady and physiological pulsatile flow situations and the measurements included the pressure drop across the valve, the instantaneous flow rate as well as the velocity profiles and turbulent stresses downstream to the valves. Our study shows that the velocity profiles with pericardial valves are closer to those measured past natural aortic valves. The porcine valves with a smaller valve opening area produce a narrower and stronger jet downstream from the valve with relatively larger turbulent axial stresses in the boundary of the jet. Our study suggests that the pericardial valves with turbulent stresses comparable to those of caged ball and tilting disc valves are preferable from a hemodynamic point of view.     

Direction-dependent constriction flow in a poroelastic solid: the intervertebral disc valve

We hypothesize that a direction-dependent flow resistance exists in the intervertebral disc due to constriction flow in the cartilage endplates. A comparison of the hydrostatic pressure in the nucleus of the healthy intervertebral disc during daily loading with the relatively low osmotic swelling pressure during rest, suggests the necessity of such direction-dependent flow resistance to ensure that all the fluid exuded from the disc during loading is recovered during rest. A physical model demonstrating the direction-dependent resistance of constriction flow in a poroelastic solid is presented. A finite element model was developed and validated against this physical model. The finite element model showed that decrease of the constriction hole area not only increases the resistance to fluid flow, but also causes the direction-dependency of flow resistance to decrease. Through this mechanism, endplate sclerosis could affect normal daily fluid exchange in the intervertebral disc, resulting in decreased mass transport and/or dehydration of the disc.     

Monte Carlo-based assessment of the trade-off between spatial resolution,field-of-view and scattered radiation in the variable resolution X-ray CT scanner

ObjectiveThe purpose of this work is to evaluate the impact of optimization of magnification on performance parameters of the variable resolution X-ray (VRX) CT scanner.MethodsA realistic model based on an actual VRX CT scanner was implemented in the GATE Monte Carlo simulation platform. To evaluate the influence of system magnification, spatial resolution, field-of-view (FOV) and scatter-to-primary ratio of the scanner were estimated for both fixed and optimum object magnification at each detector rotation angle. Comparison and inference between these performance parameters were performed angle by angle to determine appropriate object position at each opening half angle.ResultsOptimization of magnification resulted in a trade-off between spatial resolution and FOV of the scanner at opening half angles of 90°–12°, where the spatial resolution increased up to 50% and the scatter-to-primary ratio decreased from 4.8% to 3.8% at a detector angle of about 90° for the same FOV and X-ray energy spectrum. The disadvantage of magnification optimization at these angles is the significant reduction of the FOV (up to 50%). Moreover, magnification optimization was definitely beneficial for opening half angles below 12° improving the spatial resolution from 7.5 cy/mm to 20 cy/mm. Meanwhile, the FOV increased by more than 50% at these angles.ConclusionIt can be concluded that optimization of magnification is essential for opening half angles below 12°. For opening half angles between 90° and 12°, the VRX CT scanner magnification should be set according to the desired spatial resolution and FOV.     

Complex distributions of residual stress and strain in the mouse left ventricle: experimental and theoretical models

 Most soft biological tissues, including ventricular myocardium, are not stress free when all external loads are removed. Residual stress has implications for mechanical performance of the heart, and may be an indicator of patterns of regional growth and remodeling. Cross-sectional rings of arrested ventricles opened up when a radial cut was made (initial mean opening angles were 64 ± 17°), but further circumferential cuts revealed the presence of additional residual stresses in the tissue with further opening of the rings. In normal mouse hearts, the inner half of a short-axis ring opened more than the outer half, and this change was dependent on apex–base location. At the apex the inner section vs. outer section opening angles were 226 ± 47° vs. 89 ± 28°, while at the base the same two angles were 160 ± 30° vs. 123 ± 35°. A simple theoretical cylindrical shell model with incompressible hyperelastic material properties was used to model the experimental deformations based on the cutting experiments. The model predicts different residual stress fields depending on the nature of the opening after the circumferential cut (which is done after the conventional radial cut). The observed opening angles were consistent with steep stress gradients near the endocardium compared with those predicted if the first cut was assumed to relieve all residual stresses. These results imply a more complex distribution of residual stress and strain in ventricular myocardium than previously thought. Received: 23 May 2002 / Accepted: 30 September 2002 We would like to acknowledge the surgical skills and data analysis of Zuangjie Li. This work was supported in part by National Heart, Lung, and Blood Institute Grants HL-43026 and HL-64321.     

Effect of collagen fibre orientation on intervertebral disc torsion mechanics

The intervertebral disc is a complex fibro-cartilaginous material, consisting of a pressurized nucleus pulposus surrounded by the annulus fibrosus, which has an angle-ply structure. Disc injury and degeneration are noted by significant changes in tissue structure and function, which significantly alters stress distribution and disc joint stiffness. Differences in fibre orientation are thought to contribute to changes in disc torsion mechanics. Therefore, the objective of this study was to evaluate the effect of collagen fibre orientation on internal disc mechanics under compression combined with axial rotation. We developed and validated a finite element model (FEM) to delineate changes in disc mechanics due to fibre orientation from differences in material properties. FEM simulations were performed with fibres oriented at \(\pm 30^{\circ }\) throughout the disc (uniform by region and fibre layer). The initial model was validated by published experimental results for two load conditions, including \(0.48\,\hbox {MPa}\) axial compression and \(10\,\hbox {Nm}\) axial rotation. Once validated, fibre orientation was rotated by \(4^{\circ }\) or \(8^{\circ }\) towards the horizontal plane, resulting in a decrease in disc joint torsional stiffness. Furthermore, we observed that axial rotation caused a sinusoidal change in disc height and radial bulge, which may be beneficial for nutrient transport. In conclusion, including anatomically relevant fibre angles in disc joint FEMs is important for understanding stress distribution throughout the disc and will be important for understanding potential causes for disc injury. Future models will include regional differences in fibre orientation to better represent the fibre architecture of the native disc.     

Growth and remodeling in a thick-walled artery model: effects of spatial variations in wall constituents

A mathematical model is presented for growth and remodeling of arteries. The model is a thick-walled tube composed of a constrained mixture of smooth muscle cells, elastin and collagen. Material properties and radial and axial distributions of each constituent are prescribed according to previously published data. The analysis includes stress-dependent growth and contractility of the muscle and turnover of collagen fibers. Simulations were conducted for homeostatic conditions and for the temporal response following sudden hypertension. Numerical pressure-radius relations and opening angles (residual stress) show reasonable agreement with published experimental results. In particular, for realistic material and structural properties, the model predicts measured variations in opening angles along the length of the aorta with reasonable accuracy. These results provide a better understanding of the determinants of residual stress in arteries and could lend insight into the importance of constituent distributions in both natural and tissue-engineered blood vessels.     

Estimation of the rotational undamped natural frequency of bileaflet cardiac valve prostheses

The angular momentum balance is solved numerically for a size 29 mm CarboMedics prosthetic heart valve. The lift force is estimated from potential flow theory, while the drag force is estimated from the lift force and a blunt body empiricism. Buoyancy and gravitational effects are calculated based on the assumption of homogeneous leaflets. Other assumptions include uniform flow, negligible friction at the pivot axis, negligible viscous damping and fluid inertance, and a symmetry flow condition. Oscillations are predicted in the opening dynamics in the range of 2-25 Hz, for flow rates through one-half of the orifice in the range of 0.1-10.0 l/min. The frequency of these oscillations is dependent upon the orientation of the leaflet in relation to the gravitational field and the magnitude of the flow rate. In vivo and in vitro measurements by other investigators demonstrate similar effects of gravity and flow rate, with flutter frequencies of the order of 10-100 Hz. Excitation frequencies, based on vortex shedding, are estimated to be of the order 2-200 Hz, for the range of flow rates of the theoretical model. These results suggest that the natural frequency of this rotational second order system may, in theory, be a contributing factor to the flutter observed in bileaflet cardiac valve prostheses. The clinical significance of this finding is yet to be established.     

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.     

Effect of residual stress and heterogeneity on circumferential stress in the arterial wall

Quantifying the stress distribution through the arterial wall is essential to studies of arterial growth and disease. Previous studies have shown that both residual stress, as measured by opening angle, and differing material properties for the media-intima and the adventitial layers affect the transmural circumferential stress (sigma theta) distribution. Because a lack of comprehensive data on a single species and artery has led to combinations from multiple sources, this study determined the sensitivity of sigma theta to published variations in both opening angle and layer thickness data. We fit material properties to previously published experimental data for pressure-diameter relations and opening angles of rabbit carotid artery, and predicted sigma theta through the arterial wall at physiologic conditions. Using a one-layer model, the ratio of sigma theta at the internal wall to the mean sigma theta decreased from 2.34 to 0.98 as the opening angle increased from 60 to 130 deg. In a two-layer model using a 95 deg opening angle, mean sigma theta in the adventitia increased (112 percent for 25 percent adventitia) and mean sigma theta in the media decreased (47 percent for 25 percent adventitia). These results suggest that both residual stress and wall layers have important effects on transmural stress distribution. Thus, experimental measurements of loading curves, opening angles, and wall composition from the same species and artery are needed to accurately predict the transmural stress distribution in the arterial wall.     

Fluid flow and convective transport of solutes within the intervertebral disc

Previous experimental and analytical studies of solute transport in the intervertebral disc have demonstrated that for small molecules diffusive transport alone fulfils the nutritional needs of disc cells. It has been often suggested that fluid flow into and within the disc may enhance the transport of larger molecules. The goal of the study was to predict the influence of load-induced interstitial fluid flow on mass transport in the intervertebral disc.An iterative procedure was used to predict the convective transport of physiologically relevant molecules within the disc. An axisymmetric, poroelastic finite-element structural model of the disc was developed. The diurnal loading was divided into discrete time steps. At each time step, the fluid flow within the disc due to compression or swelling was calculated. A sequentially coupled diffusion/convection model was then employed to calculate solute transport, with a constant concentration of solute being provided at the vascularised endplates and outer annulus. Loading was simulated for a complete diurnal cycle, and the relative convective and diffusive transport was compared for solutes with molecular weights ranging from 400 Da to 40 kDa.Consistent with previous studies, fluid flow did not enhance the transport of low-weight solutes. During swelling, interstitial fluid flow increased the unidirectional penetration of large solutes by approximately 100%. Due to the bi-directional temporal nature of disc loading, however, the net effect of convective transport over a full diurnal cycle was more limited (30% increase). Further study is required to determine the significance of large solutes and the timing of their delivery for disc physiology.     

Zero-stress state of intra- and extraparenchymal airways from human, pig, rabbit, and sheep lung.

Alterations in airway wall anatomic properties and the consequential effects on airway narrowing have been assessed by use of computational models. In these models, it is generally assumed that at zero transmural pressure the airway wall exists in a zero-stress state. Many studies have shown that this is often not the case, as evidenced by a nonzero opening angle. In this study, we measured the opening angle of airway rings at zero transmural pressure to test this assumption. The airway tree was dissected from human, pig, sheep, and rabbit lungs. Airways were excised from the tree, and the opening angle was measured. There were obvious species and regional differences in opening angle. Rabbit airways from both extraparenchymal and intraparenchymal sites exhibited marked opening angles (7-82 degrees). Extraparenchymal airways from sheep had large opening angles (up to 50 degrees), but ovine intraparenchymal airways had small opening angles. Measurable opening angles were rarely observed in human and porcine airways of any size. The assumption of a stable zero-stress state at zero transmural pressure is therefore valid for human and porcine, but not rabbit and sheep, airways.     

Steady flow development past valve prostheses in a model human aorta—II. Tilting disc valves

The flow development in the model human aorta with uniform entry as well as with centrally occuluding valves mounted at the root of the aorta was described in Part I of this two-paper sequence. Part II deals with the flow development in the model aorta with tilting disc valves mounted at the root of the aorta. Bjork-Shiley and Hall-Kaster tilting disc valves were mounted in three different orientations with respect to the root of the aorta. The velocity profiles and turbulent stresses were measured with laser-Doppler anemometry. Our results under steady flow conditions in the model human aorta show quantitatively that the flow development in the ascending aorta as well as in the brachio-cephalic artery are strongly dependent on the orientation of the tilting disc valves. With the valves tilting towards the outer wall of curvature, our results suggest a tendency for flow separation at the flow divider region of the brachio-cephalic artery.     

Numerical Models for the Simulation of Flexible Artificial Heart Valves: Part II - Valve Studies

Abstract

Using a coupled Lagrangian dynamic leaflet model and an unsteady potential flow solver the motion of a polyurethane type heart valve is simulated in the aortic position. The simulations incorporate two flow domains; the first comprises only the leaflets which are embedded within an unsteady flow of infinite expanse, and the second incorporates the influence of the aortic geometry via a conformal mapping. Simulations are performed for a cardiac output of 51itres/min and a beat period of 72 b.p.m. corresponding to a typical aortic pulse. Resulting valve motions are computed for various leaflet bending stiffnesses in both flow domains. In addition both the bending stress and strain and their time rate of change are evaluated. Valve motion displays the characteristic rapid opening, stable opening and slow closing phases as detailed in the literature. The computed stress values along the leaflet surface are of the order of those found experimentally.     

A mechanism for morphogen-controlled domain growth

Many developmental systems are organised via the action of graded distributions of morphogens. In the Drosophila wing disc, for example, recent experimental evidence has shown that graded expression of the morphogen Dpp controls cell proliferation and hence disc growth. Our goal is to explore a simple model for regulation of wing growth via the Dpp gradient: we use a system of reaction-diffusion equations to model the dynamics of Dpp and its receptor Tkv, with advection arising as a result of the flow generated by cell proliferation. We analyse the model both numerically and analytically, showing that uniform domain growth across the disc produces an exponentially growing wing disc.     

Development of wall surface tangent DPIV measurement techniques for arterial branch models.

Experimental techniques for measuring unsteady flow in a glass arterial bifurcation model have been developed to aid in quantifying three-dimensional wall shear fluctuations associated with arterial disease. The unique feature of the current technique is the use of a "curved" laser sheet, which was everywhere tangent to the inner wall of a daughter tube in an arterial bifurcation model. Surface tangent velocity vector field measurements were made to demonstrate the potential of this technique. Ensemble-averaged data showing weak secondary flows as well as statistical distributions of flow angles are presented. Measurements of this type may be used to estimate mean and instantaneous wall shear magnitude and direction, data that are necessary for understanding the importance of circumferential motions on arterial disease.     

Ferrowax microvalves for fully automated serial dilution on centrifugal microfluidic platforms

We herein describe a centrifugal microfluidic system to accomplish a fully automated serial dilution. The liquid flow on the disc was regulated by utilizing ferrowax microvalves systematically integrated into the channels within specially designed metering structures. By opening the differently positioned microvalves through irradiation of IR laser to allow metering, the same amount of diluent was serially eluted to the dilution chamber from the same diluent chamber. After dilution, the diluted samples were automatically delivered to the respective final product chambers by appropriately opening or closing the microvalves in the connecting channels, followed by rotating the disc. Based on this unique design principle, six consecutive two-fold and 10-fold dilutions were successfully achieved, yielding excellent accuracy in a wide dynamic range up to six orders of magnitude. Very importantly, the overall serial dilution process, including the diluent addition, mixing, and product transfer steps, was completed very rapidly within 5 min, due to the minimized procedures enabled by the automated actuation of the ferrowax microvalves at the rationally designed positions. We expect our centrifugal microfluidic system would serve as a powerful elemental tool to realize fully automated diagnostic microsystems involving the serial dilution process.     

Geometry strongly influences the response of numerical models of the lumbar spine--a probabilistic finite element analysis

Typical FE models of the human lumbar spine consider a single, fixed geometry. Such models cannot account for potential effects of the natural variability of the spine's geometry. In this study, we performed a probabilistic uncertainty and sensitivity analysis of a fully parameterized, geometrically simplified model of the L3-L4 segment. We examined the impact of the uncertainty in all 40 geometry parameters, estimated lower and upper bounds for the required sample size and determined the most important geometry parameters. The natural variability of the spine's geometry indeed strongly affects intradiscal pressure, range of motion and facet joint contact forces. Deriving generalized statements from fixed-geometry models as well as transferring those results to different cases thus can easily lead to wrong conclusions and should only be performed with extreme caution. We recommend a sample size of ≈ 100 to obtain reasonable accurate point estimates and a sufficient overview of the remaining uncertainties. Yet, only few parameters, especially those determining the disc geometry (disc height, end-plate width and depth) and the facets' position (intra-articular space, pedicle length, facet angles), proved to be truly important. Accurate measurement and modeling of those structures should therefore be prioritized.     

Effect of prosthetic mitral valve geometry and orientation on flow dynamics in a model human left ventricle

Pulsatile flow dynamics through bileaflet (St Jude and Duromedics), tilting disc (Bjork-Shiley and Omniscience), caged ball (Starr-Edwards), pericardial (Edwards) and porcine (Carpentier-Edwards) mitral valves in a model human left ventricle (LV) were studied. The model human ventricle, obtained from an in situ diastolic casting, was incorporated into a mock circulatory system. Measurements were made at various heart rates and flow rates. These included the transvalvular pressure drop and regurgitation in percent and cm3 beat-1. The effect of valve geometry and the orientation of the valve with respect to the valve annulus was analyzed using a flow visualization technique. Qualitative flow visualization study indicates certain preferred orientations for the tilting disc and bileaflet valve prostheses in order to obtain a smooth washout of flow in the LV chamber.