bluuebunny/crossref_metadata_2025 huggingface dataset:
Errors Induced by Off-Axis Measurement of the Elastic Properties of Bone
Authors: Turner C. H. , Cowin S. C.
Published: June 2009

Misalignment between the axes of measurement and the material symmetry axes of bone causes error in anisotropic elastic property measurements. Measurements of Poisson’s ratio were strongly affected by misalignment errors. The mean errors in the measured Young’s moduli were 9.5 and 1.3 percent for cancellous and cortical bone, respectively, at a misalignment angle of 10 degrees. Mean errors of 1.1 and 5.0 percent in the measured shear moduli for cancellous and cortical bone, respectively, were found at a misalignment angle of 10 degrees. Although, cancellous bone tissue was assumed to have orthotropic elastic symmetry, the possibility of the greater symmetry of transverse isotropy was investigated. When the nine orthotropic elastic constants were forced to approximate the five transverse isotropic elastic constants, errors of over 60 percent were introduced. Therefore, it was concluded that cancellous bone is truly orthotropic and not transversely isotropic. A similar but less strong result for cortical bone tissue was obtained.
Simulation of the Double Limb Support Phase of Human Gait
Authors: Ju Ming-Shaung , Mansour J. M.
Published: June 2009

Dynamic mechanical models of the double limb support phase of human gait were developed for both two-dimensional (sagittal plane) and three-dimensional motion. A “foot” model with a curved plantar surface was also developed such that the model foot motion was kinematically equivalent to that of a walking subject. This foot model was incorporated into the planar model for double limb support. The dynamic formulations were based on Kane’s method and were implemented symbolically using MACSYMA. The development of the formulations for the constrained systems, application of these formulations to the study of normal gait, the sensitivity of the simulation to the frequency content of the input data, the sensitivity of limb displacements to changes in joint moments and the application of a nonlinear feedback controller to correct for perturbations in limb trajectories were investigated.
The Dynamics of Quadrupedal Locomotion
Authors: Pandy M. G. , Kumar V. , Berme N. , Waldron K. J.
Published: June 2009

This paper presents a dynamical analysis of quadrupedal locomotion, with specific reference to an adult Nubian goat. Measurements of ground reaction forces and limb motion are used to assess variations in intersegmental forces, joint moments, and instantaneous power for three discernible gaits: walking, running, and jumping. In each case, inertial effects of the torso are shown to dominate to the extent that lower-extremity contributions may be considered negligible. Footforces generated by the forelimbs exceed those exerted by the hindlimbs; and, in general, ground reactions increase with speed. The shoulder and hip dominate mechanical energy production during walking, while the knee plays a more significant role in running. In both cases, however, the elbow absorbs energy, and by so doing functions primarily as a damping (control) element. As opposed to either walking or running, jumping requires total horizontal retardation of the body’s center of mass. In this instance, generating the necessary vertical thrust amounts to energy absorption at all joints of the lower extremities.
Description and Error Evaluation of an In Vitro Knee Joint Testing System
Authors: Lewis J. L. , Lew W. D. , Schmidt J.
Published: June 2009

An experimental system for the analysis of knee joint biomechanics is presented. The system provides for the simultaneous recording of ligament forces using buckle transducers and three-dimensional joint motion using an instrumented spatial linkage, as in vitro specimens are subjected to a variety of external loads by a pneumatic loading apparatus with associated force transducers. The system components are described, and results of an evaluation of system errors and accuracy are presented. The experimental setup has been successfully used in the analysis of normal knee ligament mechanics, as well as surgical reconstructions of the anterior cruciate ligament. The system can also be adapted to test other human or animal in vitro joints.
Convection-Diffusion Interaction for Oscillatory Flow in a Tapered Tube
Authors: Godleski D. A. , Grotberg J. B.
Published: June 2009

Transport of soluble material is analyzed for volume-cycle oscillatory flow in a tapered tube. The equations of motion are solved using a regular perturbation method for small taper angle and order unity amplitude over a range of the Womersley parameter. The transport equation is also solved by a regular perturbation method where uniform end concentrations and no wall flux are assumed. The time-averaged axial transport of solute is calculated for several tapered tubes. There is substantial modification of transport compared to the straight tube case and the results are interpreted with respect to pulmonary gas exchange.
Numerical Solutions for Steady and Unsteady Flow in a Model of the Pulmonary Airways
Authors: Kimmel E. , Kamm R. D. , Shapiro A. H.
Published: June 2009

A computational model is presented for unsteady flow through a collapsible tube with variable wall stiffness. The one-dimensional flow equations are solved for inlet, outlet and external conditions that vary with time and for a tube with time-dependent, spatially-distributed local properties. In particular, the effects of nonuniformities and local perturbations in stiffness distribution in the tube are studied. By allowing the flow to evolve in time, asymptotically steady flows are calculated. When simulating a quasi-steady reduction in downstream pressure, the model demonstrates critical transitions, the phenomena of wave-speed limitation and the sites of flow limitation. It also exhibits conditions for which viscous flow limitation occurs. Computations of rapid, unsteady changes of the exit pressure illustrate the phenomena occurring at the onset of a cough, and the generation and propagation of elastic jumps.
Shear-Layer Detection in Poststenotic Flow by Spectrum Analysis of Doppler Signals
Authors: Tamura T. , Fronek A.
Published: June 2009

Spectrum analysis of the Doppler signals was performed 0.5 tube diameters downstream from an axisymmetric constriction with an area reduction of 80 percent in steady flow at a jet Reynolds number of 2840. Both pulsed and continuous wave (CW) Doppler spectra showed significant reverse flow components in the separated flow. The pulsed Doppler spectra exhibited sudden changes when the sample volume crossed the shear layer between the center jet and the separated flow. A power spectrum equation was theoretically derived from continuity of flow to define the Doppler shift frequency for the shear layer velocity. The CW Doppler spectrum showed a minimum spectrum density at a frequency which equalled the shear layer Doppler shift frequency derived from the equation. The pulsed spectra exhibited the sudden changes at the same frequency as well.
Hot-Film Wall Shear Probe Measurements Inside a Ventricular Assist Device
Authors: Baldwin J. T. , Tarbell J. M. , Deutsch S. , Geselowitz D. B. , Rosenberg G.
Published: June 2009

Wall shear rates at eleven sites within the Penn State Electric Ventricular Assist Device (EVAD) were determined with the pump operating under conditions of 30 and 50 percent systolic duration and a mean flow rate of 5.8 L/min using a flush-mounted hot-film probe. Probe calibrations were performed with the hot-film in two orientations relative to the flow direction: a standard orientation and an orientation in which the hot-film was rotated by 90 deg from the standard orientation. The magnitude and direction of the wall shear stress at each site within the EVAD were estimated from ensemble averaged voltage data recorded for similar standard and rotated film orientations. The results indicate that, during diastole the wall shear stress direction around the pump’s periphery for both operating conditions is predominantly perpendicular to the inflow-outflow plane (in the direction of the pusher plate motion) and reaches a peak value of approximately 350 dynes/cm2. The highest wall shear stresses were found near the prosthetic aortic valve (inside the EVAD) under the 30 percent systolic duration condition and are estimated to be as high as 2700 dynes/cm2. Peak shear stress values of 1400 dynes/cm2 were observed in the vicinity of the prosthetic mitral valve under both operating conditions. The results suggested that the valve regions are substantially more hemolytic than other wall regions of the EVAD; the magnitudes of the wall shear stresses are sensitive to operating conditions; and that wall shear in the direction of pusher plate motion can be significant.
Model Studies at Mechanical Aortic Heart Valve Prostheses—Part I: Steady-State Flow Fields and Pressure Loss Coefficients
Authors: Knoch Martin , Reul Helmut , Kro¨ger Ralf , Rau Gu¨nter
Published: June 2009

In a 3:1 scaled model of the human aorta models of the Omniscience, Bjo¨rk-Shiley Convexo-Concave, Bjo¨rk-Shiley Monostrut, Medtronic-Hall, Duromedics (Hemex) and the Saint Jude Medical heart valve prostheses are studied in steady flow representing the systolic peak flow phase. Detailed flow visualization experiments show flow separations at all inner ring surfaces as well as at most of the occluders. The resulting stagnation areas increase the risk of thrombus accumulation. Flow separations also stimulate vortex formation and turbulent mixing at the downstream jet boundaries and thus may intensify blood damage by turbulent shear stresses. The different influences of struts and occluder guides on the flow around the occluders are discussed. The effects of the individual valve components on the flow fields are analyzed and correlated with the resulting pressure losses.
Finite Elasticity Modeling of the Biaxial and Uniaxial Properties of Compliant Vascular Grafts
Authors: Stewart S. F. C. , Lyman D. J.
Published: June 2009

Compliant vascular grafts were modeled by finite elasticity theory. A linear, biaxial model satisfactorily described the stress-strain behavior in inflation tests, where the sample length was fixed longitudinally and inflated. The model was then used to predict the behavior in a longitudinal test (where the longitudinal stretch was varied while keeping the pressure zero), and in a uniaxial test of a circumferential strip. The model satisfactorily predicted the longitudinal Young’s modulus measured in the longitudinal test. The model was less successful in predicting the circumferential Young’s modulus measured in the uniaxial test, possibly because the state of stress was not purely uniaxial.
The Micromechanics Versus the Macromechanics of Cortical Bone—A Comprehensive Presentation
Authors: Ascenzi A.
Published: June 2009

Secondary cortical bone is a complicated patchwork of structures which can be viewed as a hierarchy of four different orders. As far as the biomechanical properties of cortical bone are concerned, the lamella is the most important of the four. The relative distribution of longitudinal lamellae (whose fiber bundles and crystallites have a longitudinal course and withstand loading by tension) with respect to transverse lamellae (whose fiber bundles and crystallites have a transverse course and withstand loading by compression) governs the mechanical properties of bone at macroscopic level both in normal and pathological conditions.
The Three-Dimensional Kinematics and Flexibility Characteristics of the Human Ankle and Subtalar Joints—Part I: Kinematics
Authors: Siegler Sorin , Chen Jie , Schneck C. D.
Published: June 2009

The in-vitro, three dimensional kinematic characteristics of the human ankle and subtalar joint were investigated in this study. The main goals of this investigation were: 1) To determine the range of motion of the foot-shank complex and the associated range of motion of the ankle and subtalar joints; 2) To determine the kinematic coupling characteristics of the foot-shank complex, and 3) To identify the relationship between movements at the ankle and subtalar joints and the resulting motion produced between the foot and the shank. The tests were conducted on fifteen fresh amputated lower limbs and consisted of incrementally displacing the foot with respect to the shank while the motion of the articulating bones was measured through a three dimensional position data acquisition system. The kinematic analysis was based on the helical axis parameters describing the incremental displacements between any two of the three articulating bones and on a joint coordinate system used to describe the relative position between the bones. From the results of this investigation it was concluded that: 1) The range of motion of the foot-shank complex in any direction (dorsiflexion/plantarflexion, inversion/eversion and internal rotation/external rotation) is larger than that of either the ankle joint or the subtalar joint.; 2) Large kinematic coupling values are present at the foot-shank complex in inversion/eversion and in internal rotation/external rotation. However, only a slight amount of coupling was observed to occur in dorsiflexion/plantarflexion.; 3) Neither the ankle joint nor the subtalar joint are acting as ideal hinge joints with a fixed axis of rotation.; 4) Motion of the foot-shank complex in any direction is the result of rotations at both the ankle and the subtalar joints. However, the contribution of the ankle joint to dorsiflexion/plantarflexion of the foot-shank complex is larger than that of the subtalar joint and the contribution of the subtalar joint to inversion/eversion is larger than that of the ankle joint.; 5) The ankle and the subtalar joints have an approximately equal contribution to internal rotation/external rotation movements of the foot-shank complex.
The Three-Dimensional Kinematics and Flexibility Characteristics of the Human Ankle and Subtalar Joint—Part II: Flexibility Characteristics
Authors: Chen Jie , Siegler Sorin , Schneck Carson D.
Published: June 2009

The objective of the present study was to investigate the in-vitro, coupled, three-dimensional load-displacement and flexibility characteristics of the human ankle joint complex consisting of the talocrural and the talocalcaneal joints and to determine the effects that sectioning of the anterior talofibular ligament has on these characteristics. Similar to other anatomical joints such as the knee and the intervertebral joint, the ankle joint complex was found to exhibit highly nonlinear load-displacement characteristics with the angular displacement approaching asymptotic values as the external load was increased. Therefore, a procedure of incremental linearization was used to derive the flexibility characteristics of this structure. According to this procedure, external loads were applied to the calcaneus in small increments and its resulting three dimensional displacements were recorded. The incremental flexibility coefficients were then derived by assuming linear load-displacement relationship for each increment. From the results obtained from fifteen human ankle specimens, it was evident that the ankle joint complex exhibit highly coupled flexibility and load-displacement characteristics. It was further concluded that the ankle joint complex is the most flexible in the neighborhood of the unloaded, neutral position and that all the flexibility coefficients of the structure decrease rapidly toward the extremes of the range of motion. Rupture of the anterior talofibular ligament was found to have a significant effect on the load-displacement and flexibility characteristics of the ankle joint complex. This effect was manifested as a change in the load-displacement characteristics and a large increase in the flexibility coefficients primarily in those corresponding to rotations in the transverse and the coronal plane. The results of the present study can provide the necessary data base for the development of quantitative diagnostic technique for identifying the site and the extent of injury to the collateral ligaments of the ankle.
Continuous Monitoring of Erythrocyte Sedimentation Process: A New Possible Mechanism of Erythrocyte Sedimentation
Authors: Kuo Cheng-Deng , Bai Jing-Jang , Chang I-Te , Wang Jia-Horng , Chien Shu
Published: June 2009

An automated system is constructed to record the complete course of erythrocyte sedimentation process. In this system a light source and a paired photodetector are employed to monitor the change of light transmittance at the junction of plasma and the sedimenting red blood cell column, thus providing a continuous record of erythrocyte sedimentation as a function of time. Differentiation of this sedimentation—time curve yields a velocity—time curve of erythrocyte sedimentation. Frequently recorded “spikes” on top of the velocity—time curve imply the episodes of very rapid fall of erythrocytes in the sedimentation tube that cannot be explained by the currently accepted theory of erythrocyte sedimentation based mainly on Stokes’ law, and a new mechanism of rouleau coalescing and fracturing is proposed to account for them.
Measurement of Local Mass Transfer Coefficients in a Cast Model of the Human Upper Respiratory Tract
Authors: Hanna L. M. , Scherer P. W.
Published: June 2009

Local mass transfer coefficients measured using the naphthalene sublimation technique in an acrylic cast model of the human upper respiratory tract are reported as the Sherwood numbers for the corresponding regions. A steady air flow rate of 12 L per min was used for all measurements. Values of the Sherwood number are seen to be highest in the nasal cavity and proximal nasopharynx while a minimum value occurs just downstream from the larynx. Local values of the Nusselt number obtained in the trachea and proximal nasal cavity assuming a complete heat and mass transfer analogy agree well with in-vivo physiological measurements. The mass transfer coefficients found can be incorporated into an analytical model of respiratory heat and water vapor transfer or into a model of pollutant gas uptake in the respiratory tract.
Material Constants for a Finite Element Model of the Intervertebral Disk With a Fiber Composite Annulus
Authors: Spilker R. L. , Jakobs D. M. , Schultz A. B.
Published: June 2009

A simple axisymmetric finite element model of a human spine segment containing two adjacent vertebrae and the intervening intervertebral disk was constructed. The model incorporated four substructures: one to represent each of the vertebral bodies, the annulus fibrosus, and the nucleus pulposus. A semi-analytic technique was used to maintain the computational economies of a two-dimensional analysis when nonaxisymmetric loads were imposed on the model. The annulus material was represented as a layered fiber-reinforced composite. This paper describes the selection of material constants to represent the anisotropic layers of the annulus. It shows that a single set of material constants can be chosen so that model predictions of gross disk behavior under compression, torsion, shear, and moment loading are in reasonable agreement with the mean and range of experimentally measured disk behaviors. It also examines the effects of varying annular material properties.
O2 Transport in Cerebral Microregions (Mathematical Simulation)
Authors: Kislyakov Yu. Ya. , Ivanov K. P.
Published: June 2009

The effects of the circulation rate in capillaries, the intensity of O2 consumption by nerve cells and the capillary network density on the O2 tension distribution in the cerebral cortex have been studied, utilizing a mathematical model simulating actual neuron-capillary relationships. The model has been written as a system of equations in partial derivatives, its solution obtained by the net-point method. Regulatory variations of the capillary circulation rate in certain cerebral microregions have been shown to ensure similar changes in oxygen supply throughout the region. A drop of the pO2 level in a cerebral microregion with a rising O2 consumption by nerve cells is shown to be due, by 75 percent, to the increase of O2 consumption and by 25 percent, to the lower pO2 in the capillaries. Conversely, an increase in pO2 in microregions resulting from a lower O2 consumption by neurons is due by 75 percent, to a pO2 rise in capillaries and by 25 percent, at the expense of an O2 consumption decrease. In cerebral regions differing in capillary network density by 20 percent, changes in the conditions for oxygen supply to tissue are due by 1/3 to pO2 variations in the capillaries and by 2/3 to alterations in the diffusion distances.
Investigation of Temperature Fields Around Embedded Cryoprobes
Authors: Budman H. , Shitzer A. , Del Giudice S.
Published: June 2009

The temperature fields around cryoprobes were investigated analytically and experimentally. Two cryoprobes were employed: a spherically shaped general purpose probe utilizing liquid nitrogen and a cylindrical “glaucoma” probe utilizing the Joule-Thomson effect in gaseous CO2. Both probes were operated by commercial cryostats. The analytical solutions included a one-dimensional integral solution for the general purpose cryoprobe, and finite element solutions for both cryoprobes. Both solutions were based on the enthalpy method. Analytical and experimental results compared reasonably well. Deviations of these results are believed to be due, mainly, to the incomplete specification of the boundary conditions on the surface of the cryoprobe.
Blood Platelet Surface Interactions on Fibrinogen Under Flow as Viewed With Fluorescent Video-Microscopy
Authors: Feuerstein I. A. , Kush J.
Published: June 2009

The interaction of fluorescently labeled blood platelets with fibrinogen-coated glass was studied in Poiseuille flow at 3 wall shear rates, 40, 80 and 944 s−1. Observations were made via video-microscopy at a distance of 0.5 cm from a tube’s entrance over a 1370 μm2 portion of luminal area. The rates of arrival and detachment, and the net rate of adhesion of cells increased nonlinearly with flow rate. The fraction of arriving cells, first contacts, which adhered without subsequent movement and the fraction of arriving cells which adhered, moved to new positions and then remained adherent, were maximal at 80 s−1. For platelets which adhere and then move to a number of new positions, the likelihood of permanent adhesion is greater than 85 percent. The adhesion process is one in which 40–60 percent of cells permanently adhere on first contact with an additional 30 percent adhering after several moves along the surface. Cells contacting where a platelet was previously adherent had a greater chance of adhering than they would on an unaltered fibrinogen surface. The efficiency of platelet adhesion is greater for second contacts than for first contacts on unaltered fibrinogen coated surface.
A Sensitivity Analysis of the Thermal Pulse Decay Method for Measurement of Local Tissue Conductivity and Blood Perfusion
Authors: Arkin H. , Holmes K. R. , Chen M. M.
Published: June 2009

The Thermal Pulse Decay (TPD) method for the determination of local tissue thermal conductivity and blood perfusion rate is based on a comparison of measured with theoretically calculated temperatures. A sensitivity analysis of the theoretical model is performed. This analysis supports the establishment of an experimental protocol which reduces the measurement errors: An “optimal” measurement time inverval for typical perfusion rates (up to 6 mL/mL/min) was found to be between 3 and 11 s after the heat pulse is turned off. Within this interval, the maximum error in determination of tissue conductivity and blood perfusion caused by experimental measurement errors is expected not to exceed 5 percent. The presently chosen pulse duration of 3 s is in agreement with the analysis as a good compromise between accuracy and excessive tissue heating.
Two-Component Laser Velocimeter Measurements Downstream of Heart Valve Prostheses in Pulsatile Flow
Authors: Tiederman W. G. , Steinle M. J. , Phillips W. M.
Published: June 2009

Elevated turbulent shear stresses resulting from disturbed blood flow through prosthetic heart valves can cause damage to red blood cells and platelets. The purpose of this study was to measure the turbulent shear stresses occurring downstream of aortic prosthetic valves during in-vitro pulsatile flow. By matching the indices of refraction of the blood analog fluid and model aorta, correlated, simultaneous two-component laser velocimeter measurements of the axial and radial velocity components were made immediately downstream of two aortic prosthetic valves. Velocity data were ensemble averaged over 200 or more cycles for a 15-ms window opened at peak systolic flow. The systolic duration for cardiac flows of 8.4 L/min was 200 ms. Ensemble-averaged total shear stress levels of 2820 dynes/cm2 and 2070 dynes/cm2 were found downstream of a trileaflet valve and a tilting disk valve, respectively. These shear stress levels decreased with axial distance downstream much faster for the tilting disk valve than for the trileaflet valve.
Swing Phase Simulation and Design of Above Knee Prostheses
Authors: Tsai C.-S. , Mansour J. M.
Published: June 2009

A detailed dynamic model of the stump-prosthesis system for an above knee amputee was developed. The model was used to examine the influence of controls and design parameters on the limb system performance during the swing phase of gait. The model duplicated the clinically known fact that hydraulic knee controllers allow the amputee to change walking speed while mechanical knee controllers limit the amputee to a single walking speed. Contrary to current practice, the simulations suggest that light weight prosthesis designs do not perform as well as heavier designs. A simple design based on a constant friction knee is shown to yield good overall performance.
Wolff’s Law of Trabecular Architecture at Remodeling Equilibrium
Authors: Cowin S. C.
Published: June 2009

An elastic constitutive relation for cancellous bone tissue is developed. This relationship involves the stress tensor T, the strain tensor E and the fabric tensor H for cancellous bone. The fabric tensor is a symmetric second rank tensor that is a quantitative stereological measure of the microstructural arrangement of trabeculae and pores in the cancellous bone tissue. The constitutive relation obtained is part of an algebraic formulation of Wolff’s law of trabecular architecture in remodeling equilibrium. In particular, with the general constitutive relationship between T, H and E, the statement of Wolff’s law at remodeling equilibrium is simply the requirement of the commutativity of the matrix multiplication of the stress tensor and the fabric tensor at remodeling equilibrium, T* H* = H* T*. The asterisk on the stress and fabric tensor indicates their values in remodeling equilibrium. It is shown that the constitutive relation also requires that E* H* = H* E*. Thus, the principal axes of the stress, strain and fabric tensors all coincide at remodeling equilibrium.
Mechanical Effects of Compression of Peripheral Nerves
Authors: Dahlin L. B. , Danielsen N. , Ehira T. , Lundborg G. , Rydevik B.
Published: June 2009

Effects of graded compression on nerve function were analyzed in order to evaluate the relative importance of pressure level and duration of compression for functional deterioration. The pressure was applied by means of a small inflatable cuff. The effects of two pressure levels, i.e., 80 mm Hg applied for 2 hr or 400 mm Hg applied for 15 min, were studied in rabbit tibial nerves. The lower pressure tested, which is known to induce ischemia of the compressed nerve segment, also causes some degree of mechanical deformation of the nerve trunk, which leads to incomplete recovery following pressure release. The duration of compression is of importance for the degree of nerve injury even at the higher pressure level tested.
The Apparent Viscoelastic Behavior of Articular Cartilage—The Contributions From the Intrinsic Matrix Viscoelasticity and Interstitial Fluid Flows
Authors: Mak A. F.
Published: June 2009

Articular cartilage was modeled rheologically as a biphasic poroviscoelastic material. A specific integral-type linear viscoelastic model was used to describe the constitutive relation of the collagen-proteoglycan matrix in shear. For bulk deformation, the matrix was assumed either to be linearly elastic, or viscoelastic with an identical reduced relaxation spectrum as in shear. The interstitial fluid was considered to be incompressible and inviscid. The creep and the rate-controlled stressrelaxation experiments on articular cartilage under confined compression were analyzed using this model. Using the material data available in the literature, it was concluded that both the interstitial fluid flow and the intrinsic matrix viscoelasticity contribute significantly to the apparent viscoelastic behavior of this tissue under confined compression.
Fatigue Life Estimation Procedures for the Endurance of a Cardiac Valve Prosthesis: Stress/Life and Damage-Tolerant Analyses
Authors: Ritchie R. O. , Lubock P.
Published: June 2009

Projected fatigue life analyses are performed to estimate the endurance of a cardiac valve prosthesis under physiological environmental and mechanical conditions. The analyses are conducted using both the classical stress-strain/life and the fracture mechanics-based damage-tolerant approaches, and provide estimates of expected life in terms of initial flaw sizes which may pre-exist in the metal prior to the valve entering service. The damage-tolerant analysis further is supplemented by consideration of the question of “short cracks,” which represents a developing area in metal fatigue research, not commonly applied to data in standard engineering design practice.
A Finite-Element Model of Blood Flow in Arteries Including Taper, Branches, and Obstructions
Authors: Porenta G. , Young D. F. , Rogge T. R.
Published: June 2009

A nonlinear mathematical model of arterial blood flow, which can account for tapering, branching, and the presence of stenosed segments, is presented. With the finite-element method, the model equations are transformed into a system of algebraic equations that can be solved on a high-speed digital computer to yield values of pressure and volume rate of flow as functions of time and arterial position. A model of the human femoral artery is used to compare the effects of linear and nonlinear modeling. During periods of rapid alterations in pressure or flow, the nonlinear model shows significantly different results than the linear model. The effect of a stenosis on pressure and flow waveforms is also simulated, and the results indicate that these waveforms are significantly altered by moderate and severe stenoses.
Everting (Toposcopic) Catheter for Broad Clinical Application
Authors: Shook D. R. , Doppman J. L. , Cattau E. L. , Goldstein S. R.
Published: June 2009

The advanced development of the clinical everting (toposcopic) catheter is described. A detailed discussion of the design and outline of the fabrication techniques are followed by a thorough performance evaluation and summary of the first two clinical applications. The everting element is a low-durometer thermoplastic polyurethane elastomer. Surface treatments include the bonding of a hydrophilic polymeric coating, optimized for lubricity, to the sliding internal surfaces of the catheter. Eversion pressures and infusion/aspiration flow rates have been measured under various conditions and the infusate-in-blood mixing potential investigated. A preliminary assessment is given of the clinical performance of the catheter in the vascular delivery of chemotherapy and standard endoscopic retrograde cholangiopancreatography.
An Analysis of the Sluicing Gate in Pulmonary Blood Flow
Authors: Fung Y. C. , Zhuang F. Y.
Published: June 2009

For pulmonary blood flow in zone 2 condition, in which the blood pressure in the venule (pven) is lower than the alveolar gas pressure (pA), the blood exiting from the capillary sheet and entering a venule must go through a sluicing gate. The sluicing gate exists because the venule remains patent [7] while the capillaries will collapse when the static pressure of blood falls below the alveolar gas pressure. In the original theory of sheet flow [4–6] the effect of the tension in the interalveolar septa on the flow through the sluicing gate was ignored. Since the tension multiplied by the curvature of the membrane is equivalent to a lateral pressure tending to open the gate, and since the curvature of the capillary wall is high in the gate region, this effect may be important. The present analysis improves the original theory and demonstrates that the effect of membrane tension is to cause flow to increase when the venous pressure continues to decrease. The shape of the sluicing gate resembles that of a venturi tube, and can be determined by an iterative integration of the differential equations. The result forms an important link in the theory of pulmonary blood flow in zone 2 condition.
Modeling Transient Heat Transfer Through the Skin and Superficial Tissues—1: Surface Insulation
Authors: Hodson D. A. , Eason G. , Barbenel J. C.
Published: June 2009

Two models of transient heat transfer through the skin and superficial tissues are presented. One model comprises a finite slab and semi-infinite slab, representing the epidermis and subdermal tissues, respectively, and a heat-generating interface representing the thermal effect of blood flow through the dermis. A model is also considered where the three tissue regions are represented more conventionally by three finite slabs. A transient problem arising from surface insulation is examined and analytical solutions derived from the first model are compared with numerical solutions derived from the second.
Pulsed Ultrasonic Doppler Velocity Measurements Inside a Left Ventricular Assist Device
Authors: Tarbell J. M. , Gunshinan J. P. , Geselowitz D. B. , Rosenberg G. , Shung K. K. , Pierce W. S.
Published: June 2009

In this study we have employed a single channel, pulsed ultrasonic Doppler velocimeter to measure instantaneous velocity distributions within the pumping chamber of a ventricular assist device. Instantaneous velocities have been decomposed into periodic mean and turbulent fluctuating components from which estimates of Reynolds stresses within the chamber and mean shear stresses along the wall of the chamber have been obtained. A review of the complete data set indicates a maximum value of the mean wall shear stress of 25 dynes/cm2 and a maximum Reynolds stress of 212 dynes/cm2. These values are lower than those measured distal to aortic valve prostheses in vitro and are well below levels known to damage blood components. Core flow patterns, wall washing patterns and flow stagnation points are also revealed.
Perfused Phantom Models of Microwave Irradiated Tissue
Authors: Baish J. W. , Foster K. R. , Ayyaswamy P. S.
Published: June 2009

The theoretical basis, practical design considerations, and prototype testing of a perfused model suitable for simulation studies of microwave heated tissue are presented. A parallel tube heat exchanger configuration is used to simulate the internal convection effects of blood flow. The global thermal response of the phantom, on a scale of several tube spacings, is shown theoretically to be nearly identical to that predicted by Pennes’ bioheat equation, which is known to give a reasonable representation of tissue under many conditions. A parametric study is provided for the relationships between the tube size, spacing and material properties and the simulated perfusion rate. A prototype with a physiologically reasonable perfusion rate was tested using a typical hyperthermia applicator. The measured thermal response of the phantom compares favorably with the numerical solution of the bioheat equation under the same irradiation conditions. This similarity sheds light on the unexpected success of the bioheat equation for modeling the thermal response of real tissue.
Phasic and Spatial Pressure Measurements in a Femoral Artery Branch Model for Pulsatile Flow
Authors: Back L. H. , Cho Y. I. , Crawford D. W.
Published: June 2009

Phasic and spatial time-averaged pressure distributions were measured in a 60-deg femoral artery branch model over a large range of branch flow ratios and at physiological Reynolds numbers of about 120 and 700. The results obtained with an in-vivolike flow wave form indicated spatial adverse time average pressure gradients in the branch vicinity which increased in magnitude with branch flow ratio, and the importance of the larger inertial effects at the higher Reynolds numbers. Pressure losses in the branch entrance region were relatively large, and corresponding flow resistances may limit branch flow, particularly at higher Reynolds numbers. The effect of branch flow was to reduce the pressure loss in the main lumen.
Poststenotic Flow Disturbance in the Dog Aorta as Measured With Pulsed Doppler Ultrasound
Authors: Talukder N. , Fulenwider J. T. , Mabon R. F. , Giddens D. P.
Published: June 2009

Blood flow velocity was measured in the dog aorta distal to mechanically induced constrictions of various degrees of severity employing an 8-MHz pulsed Doppler ultrasound velocimeter and a phase-lock loop frequency tracking method for extracting velocity from the Doppler quadrature signals. The data were analyzed to construct ensemble average velocity waveforms and random velocity disturbances. In any individual animal the effect of increasing the degree of stenosis beyond approximately 25 percent area reduction was to produce increasing levels of random velocity disturbance. However, variability among animals was such that the sensitivity of random behavior to the degree of stenosis was degraded to the point that it appears difficult to employ Doppler ultrasound measurements of random disturbances to discriminate among stenoses with area reductions less than approximately 75 percent. On the other hand, coherent vortex structures in velocity waveforms consistently occurred distal to mild constrictions (25–50 percent area reduction). Comparison of the phase-lock loop Doppler ultrasound data with simultaneous measurements using invasive hot-film anemometry, which possesses excellent frequency response, demonstrates that the ultrasound method can reliably detect those flow phenomena in such cases. Thus, the identification of coherent, rather than random, flow disturbances may offer improved diagnostic capability for noninvasively detecting arteriosclerotic plaques at relatively early stages of development.
The Effect of a Turbulent Jet on Gas Transport During Oscillatory Flow
Authors: Kamm R. D. , Bullister E. T. , Keramidas C.
Published: June 2009

Axial mass transport due to the combined effects of flow oscillation and a turbulent jet was studied both experimentally and with a simple theoretical model. The experiments show that the distance over which turbulence enhances transport is greatly increased by flow oscillation, and is particularly sensitive to tidal volume. The jet flow rate and jet configuration are relatively less important. To analyze the results, the region influenced by the jet is divided into two zones: a near field in which the time-mean flow velocities are larger than the turbulent fluctuations, and a far field where the time-mean flow is essentially zero. In the far field, axial mass transport is increased due to the turbulence which decays in strength away from the jet. When oscillatory flow is superimposed upon the steady jet flow, the turbulence in the far field interacts with the flow oscillations to augment the transport of turbulence energy and of mass. This transport enhancement is modeled by introducing an effective axial diffusivity analogous to that used in laminar oscillatory flow.
Inflation Waves Induced by Axial Acceleration of the Aorta
Authors: Elad D. , Foux A. , Lanir Y. , Kivity Y.
Published: June 2009

The phenomenon of high-amplitude inflation waves resulting from a sharp axial acceleration of the aorta, as may ocur in road accidents, is investigated theoretically. The aorta is modeled as an axisymmetric tapered membranic shell (tube) made of an incompressible, nonlinear viscoelastic material with cylindrical orthotropy. It is filled with an inviscid, incompressible fluid whose flow is considered as quasi-one dimensional along the tube axis. The equations of motion of the tube and of the fluid are solved numerically, by using a two-step explicit scheme, for several axial acceleration profiles. The solutions shows that an inflation wave is generated and it propagates in opposite direction to that of the acceloeration. The wall stresses, deformations and their time derivatives as well as fluid velocity and pressure are determined along the tube at different time intervals. Peak axial and circumferential stresses are high, with the latter far exceeding the former. These stresses may cause rupture of the aorta.
Branching Pattern of Pulmonary Arterial Tree in Anesthetized Dogs
Authors: Liu Y. H. , Ritman E. L.
Published: June 2009

The geometry of the pulmonary arterial tree of six adult dogs was measured by a high-speed, volume-scanning, X-ray tomographic technique. After the dogs were anesthetized a catheter was advanced to the right ventricular outlfow tract and 2 mL/kg Renovist contrast agent injected rapidly. During the subsequent pulmonary arterial phase of the angiogram the dogs were scanned. Three-dimensional geometry of the pulmonary arterial tree was measured in terms of vessel segment cross-sectional area, branching angles and interbranch segment lengths along axial pathways. The effect of lung inflation and phase of the cardiac cycle on geometry was shown to be most marked on vessel cross-sectional area. The geometric branching patterns in all dogs were similar. The observed, in-vivo branching pattern behaved somewhat like the branching pattern predicted from optimized models proposed by Murray [4, 5], Zamir [10, 11] and Uylings [7].
An Application of the Micropipette Technique to the Measurement of the Mechanical Properties of Cultured Bovine Aortic Endothelial Cells
Authors: Sato M. , Levesque M. J. , Nerem R. M.
Published: June 2009

The mechanical properties of endothelial cells were measured using the micropipette technique. The cells employed were collected from bovine aortic endothelium and cultured in our laboratory. Endothelial cells from confluent monolayers under noflow conditions were detached from their substrate by trypsin or by a mechanical method and suspended in modified Dulbecco medium (MDM). In the micropipette technique, a part of the cell is aspirated into the tip of the micropipette under a microscope, and the deformation measured from a photograph. In this study, the data obtained were analyzed using a model where the cytoskeletal elements, which are considered to be the primary stress bearing components, are assumed to reside in a submembranous, cortical layer. Detached cells were found to have almost homogeneous mechanical properties based on measurements from different regions of the surface of a single cell. However, a hysteresis loop was observed in the relation between pressure and cell deformation during the loading and unloading processes. The calculated elastic shear moduli obtained for the trypsin-detached cells were as much as 10–20 times larger than those of a red blood cell. Mechanically-detached cells had moduli approximately twice that of the trypsin detached cells. Passage time, i.e., cell culture age, had no influence on the mechanical properties of the trypsin-detached cells, but did have an effect on the mechanically-detached cells, with both the younger and older cells being somewhat stiffer.
Variation of Lumbar Spine Stiffness With Load
Authors: Edwards W. T. , Hayes W. C. , Posner I. , White A. A. , Mann R. W.
Published: June 2009

Mechanical studies of the Functional Spinal Unit (FSU) in-vitro have shown that the slopes of the load-displacement curves increase with load. This nonlinearity implies that the stiffness of the FSU is not constant over the range of physiologic loads, and that measurements obtained for FSU specimens through the application of individual loads cannot be summed to predict the response of the specimens to combined loads. Both experimental and analytical methods were developed in the present study to better quantify the nonlinear FSU load-displacement response and to calculate the coupled stiffness of FSU specimens at combined states of load reflecting in-vivo conditions. Results referenced to the center of the vertebral body indicate that lumbar FSU specimens are stiffer in flexion than in extension, and that FSU specimens loaded in flexion are stiffer at high loads than at low loads. The importance of combined load testing and a nonlinear interpretation of load-displacement data is demonstrated.
Cracks Emanating From a Fluid Filled Void Loaded in Compression: Application to the Bone-Implant Interface
Authors: Santare M. H. , Keer L. M. , Lewis J. L.
Published: June 2009

Loosening of orthopedic implants is believed to be caused, in part, by fracture at the bone-cement interface. This loosening occurs even in regions where the interfacial load is primarily compressive. A model is developed whereby cracks can radiate from an elliptical fluid filled void. The incompressible fluid is allowed to penetrate into the cracks when the system is loaded compressively. The mode I stress intensity factor is calculated to test the feasibility of crack growth, and a numerical scheme which uses piecewise quadratic polynomials is used to solve the resulting singular integral equations. The results show the combinations of parameters for which cracks are likely to grow.
The Frequency Response of Electrochemical Wall Shear Probes in Pulsatile Flow
Authors: Talbot L. , Steinert J. J.
Published: June 2009

The frequency response of surface-mounted electrochemical mass transfer probes used to deduce wall shear rates has been investigated experimentally for the case of fully developed laminar pulsatile flow in a straight tube. Generally good agreement is found with the asymptotic results obtained by Lighthill’s methods. The significance of the results with regard to the investigation of models of pulsatile flows of physiological interest is discussed. It is concluded that the frequency-dependent phase and amplitude corrections required to obtain accurate wall shear measurements are of such magnitudes as to render impractical the use of electrochemical probes to determine wall shear rates in these flows.
A Theoretical Model for Peripheral Tissue Heat Transfer Using the Bioheat Equation of Weinbaum and Jiji
Authors: Song Wei Jie , Weinbaum Sheldon , Jiji Latif M.
Published: June 2009

In this paper the new bioheat equation derived in Weinbaum and Jiji [7] is applied to the three layer conceptual model of microvascular surface tissue organization proposed in [1]. A simplified one-dimensional quantitative model of peripheral tissue energy exchange is then developed for application in limb and whole body heat transfer studies. A representative vasculature is constructed for each layer and the enhancement in the local tensor conductivity of the tissue as a function of vascular geometry and blood flow is examined. Numerical solutions for the boundary value problem coupling the three layers are presented and these results used to study the thermal behavior of peripheral tissue for a wide variety of physiological conditions from supine resting state to maximum exercise.
An Investigation of Particle Flow Through Capillary Models With the Resistive Pulse Technique
Authors: Frank R. S. , Hochmuth R. M.
Published: June 2009

The use of the resistive pulse technique for the measurement of microsphere and red cell transit times through single-pore “Nuclepore” membranes (with pore diameters of 3.5 to 7.0 μm and pore length of approximately 11 μm) is described. The investigation of the fluid mechanics and electrical characteristics of the experimental system provides methods for the determination of particle and cell size, and entrance and transit times. Experimental measurement of the position dependent velocity of spherical particles through the pore shows close agreement with theoretical models. Red cell size and transit time through different sized pores at physiological shear stresses is also measured.