Research Papers

J Biomech Eng. 2012;134(5):051001-051001-8. doi:10.1115/1.4006681.

Patient-specific computational fluid dynamics (CFD) is a powerful tool for researching the role of blood flow in disease processes. Modern clinical imaging technology such as MRI and CT can provide high resolution information about vessel geometry, but in many situations, patient-specific inlet velocity information is not available. In these situations, a simplified velocity profile must be selected. We studied how idealized inlet velocity profiles (blunt, parabolic, and Womersley flow) affect patient-specific CFD results when compared to simulations employing a “reference standard” of the patient’s own measured velocity profile in the carotid bifurcation. To place the magnitude of these effects in context, we also investigated the effect of geometry and the use of subject-specific flow waveform on the CFD results. We quantified these differences by examining the pointwise percent error of the mean wall shear stress (WSS) and the oscillatory shear index (OSI) and by computing the intra-class correlation coefficient (ICC) between axial profiles of the mean WSS and OSI in the internal carotid artery bulb. The parabolic inlet velocity profile produced the most similar mean WSS and OSI to simulations employing the real patient-specific inlet velocity profile. However, anatomic variation in vessel geometry and the use of a nonpatient-specific flow waveform both affected the WSS and OSI results more than did the choice of inlet velocity profile. Although careful selection of boundary conditions is essential for all CFD analysis, accurate patient-specific geometry reconstruction and measurement of vessel flow rate waveform are more important than the choice of velocity profile. A parabolic velocity profile provided results most similar to the patient-specific velocity profile.

Commentary by Dr. Valentin Fuster
J Biomech Eng. 2012;134(5):051002-051002-13. doi:10.1115/1.4006814.

An idealized systemic-to-pulmonary shunt anatomy is parameterized and coupled to a closed loop, lumped parameter network (LPN) in a multidomain model of the Norwood surgical anatomy. The LPN approach is essential for obtaining information on global changes in cardiac output and oxygen delivery resulting from changes in local geometry and physiology. The LPN is fully coupled to a custom 3D finite element solver using a semi-implicit approach to model the heart and downstream circulation. This closed loop multidomain model is then integrated with a fully automated derivative-free optimization algorithm to obtain optimal shunt geometries with variable parameters of shunt diameter, anastomosis location, and angles. Three objective functions: (1) systemic; (2) coronary; and (3) combined systemic and coronary oxygen deliveries are maximized. Results show that a smaller shunt diameter with a distal shunt-brachiocephalic anastomosis is optimal for systemic oxygen delivery, whereas a more proximal anastomosis is optimal for coronary oxygen delivery and a shunt between these two anatomies is optimal for both systemic and coronary oxygen deliveries. Results are used to quantify the origin of blood flow going through the shunt and its relationship with shunt geometry. Results show that coronary artery flow is directly related to shunt position.

Commentary by Dr. Valentin Fuster
J Biomech Eng. 2012;134(5):051003-051003-7. doi:10.1115/1.4006683.

The presence of bacterial biofilms is detrimental in a wide range of healthcare situations especially wound healing. Physical debridement of biofilms is a method widely used to remove them. This study evaluates the use of microfluidic jet impingement to debride biofilms. In this case, a biofilm is treated as a saturated porous medium also having linear elastic properties. A numerical modeling approach is used to calculate the von Mises stress distribution within a porous medium under fluid-structure interaction (FSI) loading to determine the initial rupture of the biofilm structure. The segregated model first simulates the flow field to obtain the FSI interface loading along the fluid-solid interface and body force loading within the porous medium. A stress-strain model is consequently used to calculate the von Mises stress distribution to obtain the biofilm deformation. Under a vertical jet, 60% of the deformation of the porous medium can be accounted for by treating the medium as if it was an impermeable solid. However, the maximum deformation in the porous medium corresponds to the point of maximum shear stress which is a different position in the porous medium than that of the maximum normal stress in an impermeable solid. The study shows that a jet nozzle of 500 μm internal diameter (ID) with flow of Reynolds number (Re) of 200 can remove the majority of biofilm species.

Commentary by Dr. Valentin Fuster
J Biomech Eng. 2012;134(5):051004-051004-10. doi:10.1115/1.4006684.

Building on previous studies in which the transport and targeting of 90 Y microspheres for liver tumor treatment were numerically analyzed based on medical data sets, this two-part paper discusses the influence of an anchored, radially adjustable catheter on local blood flow and microsphere delivery in an idealized hepatic artery system (Part I). In Part II a patient-inspired case study with necessary conditions for optimal targeting of radioactive microspheres (i.e., yttrium 90) onto liver tumors is presented. A new concept of optimal catheter positioning is introduced for selective targeting of two daughter-vessel exits potentially connected to liver tumors. Assuming laminar flow in rigid blood vessels with an anchored catheter in three controlled positions, the transient three-dimensional (3D) transport phenomena were simulated employing user-enhanced engineering software. The catheter position as well as injection speed and delivery function may influence fluid flow and particle transport. Although the local influences of the catheter may not be negligible, unique cross-sectional particle release zones exist, with which selectively the new controlled targeting methodology would allow optimal microsphere delivery. The insight gained from this analysis paves the way for improved design and testing of a smart microcatheter (SMC) system as well as new investigations leading to even more successful treatment with 90 Y microspheres or combined internal radiation and chemotherapy.

Commentary by Dr. Valentin Fuster
J Biomech Eng. 2012;134(5):051005-051005-10. doi:10.1115/1.4006685.

In this second part, the methodology for optimal tumor-targeting is further explored, employing a patient-inspired hepatic artery system which differs significantly from the idealized configuration discussed in Part I. Furthermore, the fluid dynamics of a microsphere supply apparatus is also analyzed. The best radial catheter positions and particle-release intervals for tumor targeting were determined for both the idealized and patient-inspired configurations. This was accomplished by numerically analyzing generated particle release maps (PRMs) for ten equally spaced intervals throughout the pulse. As in Part I, the effects of introducing a catheter were also investigated. In addition to the determination of micro-catheter positioning and, hence, optimal microsphere release, a microsphere-supply apparatus (MSA) was analyzed, which transports the particles to the catheter-nozzle, considering different axial particle injection functions, i.e., step, ramp, and S-curve. A refined targeting methodology was developed which demonstrates how the optimal injection region and interval can be determined with the presence of a catheter for any geometric configuration. Additionally, the less abrupt injection functions (i.e., ramp and S-curve) were shown to provide a more compact particle stream, making them better choices for targeting. The results of this study aid in designing the smart micro-catheter (SMC) in conjunction with the MSA, bringing this innovative treatment procedure one step closer to implementation in clinical practice.

Commentary by Dr. Valentin Fuster
J Biomech Eng. 2012;134(5):051006-051006-6. doi:10.1115/1.4006688.

In most finite element (FE) studies of vertebral bodies, axial compression is the loading mode of choice to investigate structural properties, but this might not adequately reflect the various loads to which the spine is subjected during daily activities or the increased fracture risk associated with shearing or bending loads. This work aims at proposing a patient-specific computer tomography (CT)-based methodology, using the currently most advanced, clinically applicable finite element approach to perform a structural investigation of the vertebral body by calculation of its full six dimensional (6D) stiffness matrix. FE models were created from voxel images after smoothing of the peripheral voxels and extrusion of a cortical shell, with material laws describing heterogeneous, anisotropic elasticity for trabecular bone, isotropic elasticity for the cortex based on experimental data. Validated against experimental axial stiffness, these models were loaded in the six canonical modes and their 6D stiffness matrix calculated. Results show that, on average, the major vertebral rigidities correlated well or excellently with the axial rigidity but that weaker correlations were observed for the minor coupling rigidities and for the image-based density measurements. This suggests that axial rigidity is representative of the overall stiffness of the vertebral body and that finite element analysis brings more insight in vertebral fragility than densitometric approaches. Finally, this extended patient-specific FE methodology provides a more complete quantification of structural properties for clinical studies at the spine.

Commentary by Dr. Valentin Fuster
J Biomech Eng. 2012;134(5):051007-051007-8. doi:10.1115/1.4006680.

Advanced prosthetic knees for transfemoral amputees are currently based on controlled damper mechanisms. Such devices require little energy to operate, but can only produce negative or zero joint power, while normal knee joint function requires alternative phases of positive and negative work. The inability to generate positive work may limit the user’s functional capabilities, may cause undesirable adaptive behavior, and may contribute to excessive metabolic energy cost for locomotion. In order to overcome these problems, we present a novel concept for an energy-storing prosthetic knee, consisting of a rotary hydraulic actuator, two valves, and a spring-loaded hydraulic accumulator. In this paper, performance of the proposed device will be assessed by computational modeling and by simulation of functional activities. A computational model of the hydraulic system was developed, with methods to obtain optimal valve control patterns for any given activity. The objective function for optimal control was based on tracking of joint angles, tracking of joint moments, and the energy cost of operating the valves. Optimal control solutions were obtained, based on data collected from three subjects during walking, running, and a sit-stand-sit cycle. Optimal control simulations showed that the proposed device allows near-normal knee function during all three activities, provided that the accumulator stiffness was tuned to each activity. When the energy storage mechanism was turned off in the simulations, the system functioned as a controlled damper device and optimal control results were similar to literature data on human performance with such devices. When the accumulator stiffness was tuned to walking, simulated performance for the other activities was sub-optimal but still better than with a controlled damper. We conclude that the energy-storing knee concept is valid for the three activities studied, that modeling and optimal control can assist the design process, and that further studies using human subjects are justified.

Commentary by Dr. Valentin Fuster

Technical Briefs

J Biomech Eng. 2012;134(5):054501-054501-7. doi:10.1115/1.4006686.

Recent studies have shown that capacitance measurements of large arteries provide better prognosis and diagnosis than tests of resistance alone in pulmonary hypertension (Mahapatra , 2006, “Relationship of Pulmonary Arterial Capacitance and Mortality in Idiopathic Pulmonary Arterial Hypertension,” J. Am. Coll. Cardiol., 47 (4), pp. 799–803; Reuben, 1971, “Compliance of the Human Pulmonary Arterial System in Disease,” Circ. Res., 29 , pp. 40–50]. Decreased arterial capacitance causes increased load to the heart and is the direct result of increased stiffness and elastic modulus of the arterial wall. Here, we validate a pressure-diameter (PD) method for comparing the elastic modulus and collagen engagement for post-hilar pulmonary arteries with a large range of arterial diameter. The tissue mechanics of the post-hilar arteries are not well-characterized in pulmonary hypertension. It is believed that future studies with this method will provide useful insight into the role of passive tissue mechanics of these arteries in the pathophysiology of pulmonary hypertension, eventually improving clinical diagnosis, prognosis, and treatment. Post-hilar pulmonary arteries, excised from healthy and hypertensive calves and healthy cows, were inflated over a range of 0 [mm Hg] to 110 [mm Hg] in an isolated tissue bath. Internal pressure was recorded with an electric pressure catheter. Artery diameter and longitudinal stretch were recorded photographically. Stress-strain data curves were extracted using Lame’s law of thick-walled tubes. Radial strips were removed from each section and tested in a uniaxial (MTS) tester for validation. Both the elastic modulus and collagen engagement strain were similar to results obtained by more traditional means. The average difference between measured values of the two methods for collagen engagement strain was 3.3% of the average value of the engagement strain. The average difference between the measured values of the two methods for modulus of elasticity was 7.4% of the average value of the modulus. The maximum, theoretical, relative error for the stress determined with the PD method was calculated at 20.3%. The PD method proved to be a suitable replacement for uniaxial strain tests in comparing collagen engagement strains. The method allowed faster testing of tissues of multiple diameters, while removing the effect of end conditions. The PD method will be of further utility in continued study of tissue mechanics in pulmonary hypertension studies.

Commentary by Dr. Valentin Fuster
J Biomech Eng. 2012;134(5):054502-054502-5. doi:10.1115/1.4006687.

Scoliosis is defined as a spinal pathology characterized as a three-dimensional deformity of the spine combined with vertebral rotation. Treatment for severe scoliosis is achieved when the scoliotic spine is surgically corrected and fixed using implanted rods and screws. Several studies performed biomechanical modeling and corrective forces measurements of scoliosis correction. These studies were able to predict the clinical outcome and measured the corrective forces acting on screws, however, they were not able to measure the intraoperative three-dimensional geometry of the spinal rod. In effect, the results of biomechanical modeling might not be so realistic and the corrective forces during the surgical correction procedure were intra-operatively difficult to measure. Projective geometry has been shown to be successful in the reconstruction of a three-dimensional structure using a series of images obtained from different views. In this study, we propose a new method to measure the three-dimensional geometry of an implant rod using two cameras. The reconstruction method requires only a few parameters, the included angle θ between the two cameras, the actual length of the rod in mm, and the location of points for curve fitting. The implant rod utilized in spine surgery was used to evaluate the accuracy of the current method. The three-dimensional geometry of the rod was measured from the image obtained by a scanner and compared to the proposed method using two cameras. The mean error in the reconstruction measurements ranged from 0.32 to 0.45 mm. The method presented here demonstrated the possibility of intra-operatively measuring the three-dimensional geometry of spinal rod. The proposed method could be used in surgical procedures to better understand the biomechanics of scoliosis correction through real-time measurement of three-dimensional implant rod geometry in vivo.

Commentary by Dr. Valentin Fuster

Sorry! You do not have access to this content. For assistance or to subscribe, please contact us:

  • TELEPHONE: 1-800-843-2763 (Toll-free in the USA)
  • EMAIL: asmedigitalcollection@asme.org
Sign In