Search

Fitting the measured bioimpedance spectroscopy (BIS) data to the Cole model and then extracting the Cole parameters is a common practice in BIS applications. The extracted Cole parameters then can be analysed as descriptors of tissue electrical properties. To have a better evaluation of physiological or pathological properties of biological tissue, accurate extraction of Cole parameters is of great importance. This paper proposes an improved Cole parameter extraction based on bacterial foraging optimization (BFO) algorithm. We employed simulated datasets to test the performance of the BFO fitting method regarding parameter extraction accuracy and noise sensitivity, and we compared the results with those of a least squares (LS) fitting method. The BFO method showed better robustness to the noise and higher accuracy in terms of extracted parameters. In addition, we applied our method to experimental data where bioimpedance measurements were obtained from forearm in three different positions of the arm. The goal of the experiment was to explore how robust Cole parameters are in classifying position of the arm for different people, and measured at different times. The extracted Cole parameters obtained by LS and BFO methods were applied to different classifiers. Two other evolutionary algorithms, GA and PSO were also used for comparison purpose. We showed that when the classifiers are fed with the extracted feature sets by BFO fitting method, higher accuracy is obtained both when applying on training data and test data. [ABSTRACT FROM AUTHOR]

Previous studies showed that all the vascular parameters from both the morphological and topological parameters were affected with the altering of imaging resolutions. However, neither the sensitivity analysis of the vascular parameters at multiple resolutions nor the distinguishability estimation of vascular parameters from different data groups has been discussed. In this paper, we proposed a quantitative analysis method of vascular parameters for vascular networks of multi-resolution, by analyzing the sensitivity of vascular parameters at multiple resolutions and estimating the distinguishability of vascular parameters from different data groups. Combining the sensitivity and distinguishability, we designed a hybrid formulation to estimate the integrated performance of vascular parameters in a multi-resolution framework. Among the vascular parameters, degree of anisotropy and junction degree were two insensitive parameters that were nearly irrelevant with resolution degradation; vascular area, connectivity density, vascular length, vascular junction and segment number were five parameters that could better distinguish the vascular networks from different groups and abide by the ground truth. Vascular area, connectivity density, vascular length and segment number not only were insensitive to multi-resolution but could also better distinguish vascular networks from different groups, which provided guidance for the quantification of the vascular networks in multi-resolution frameworks. [ABSTRACT FROM AUTHOR]

TiO2 scaffolds have previously shown to have promising osteoconductive properties in previous in vivo experiments. Appropriate mechanical stimuli can further promote this osteoconductive behaviour. However, the complex mechanical environment and the mechanical stimuli enhancing bone regeneration for porous bioceramics have not yet been fully elucidated. This paper aims to compare and evaluate mechanical environment of TiO2 scaffold with three commercial CaP biomaterials, i.e. Bio-Oss, Cerabone and Maxresorb under simulated perfusion culture conditions. The solid phase and fluid phase were modelled as linear elastic material and Newtonian fluid, respectively. The mechanical stimulus was analysed within these porous scaffolds quantitatively. The results showed that the TiO2 had nearly heterogeneous stress distributions, however lower effective Young's modulus than Cerabone and Maxresorb. The permeability and wall shear stress (WSS) for the TiO2 scaffold was significantly higher than other commercial bone substitute materials. Maxresorb and Bio-Oss showed lowest permeability and local areas of very high WSS. The detailed description of the mechanical performance of these scaffolds could help researchers to predict cell behaviour and to select the most appropriate scaffold for different in vitro and in vivo performances. Graphical abstract Schematic representation of the establishment procedure. Take the establishment process of Cerabone as an example. Left shows a slice of micro-CT image from Cerabone, and 1.5 mm × 1.5 mm region of interest was shown in the red box. A 1.5-mm3 cube was cut out by Boolean operation in Mimics (Materialise, Belgium), and the cubic model was remeshed in 3-Matic 6.0 (Materialise, Belgium). The cubic model is shown in blue, and the empty space in red. [ABSTRACT FROM AUTHOR]

The paper describes a spherical harmonic-based finite element scheme for solving Poisson-type equations throughout volumes characterised by irregularly shaped inner and outer surfaces. The inner and outer surfaces are defined by spherical harmonics, and the volume in between these surfaces is divided into nested shells that are weighted averages of the inner and outer surfaces. The resulting mesh comprises hexahedral elements, wherein each hexahedral element is defined by inner and outer shells in the radial direction and divisions in the polar and azimuthal directions. The spacing between shells can be set to any desired value. Similarly, the size of the polar and azimuthal divisions can be specified. A test of the scheme on an anisotropic sphere, meshed with 720 nodes, yielded a relative error of 0.78% on the sphere's surface. As a comparison, a previously published combined finite element/boundary element scheme with a 946-node mesh produced a corresponding error of 3.57%. [ABSTRACT FROM AUTHOR]

Electrode designs and strategies for electroneurogram recordings are often tested first by computer simulations and then by animal models, but they are rarely implanted for long-term evaluation in humans. The models show that the amplitude of the potential at the surface of an axon is higher in front of the nodes of Ranvier than at the internodes; however, this has not been investigated through in vivo measurements. An original experimental method is presented to emulate a single fiber action potential in an infinite conductive volume, allowing the potential of an axon to be recorded at both the nodes of Ranvier and the internodes, for a wide range of electrode-to-fiber radial distances. The paper particularly investigates the differences in the action potential amplitude along the longitudinal axis of an axon. At a short radial distance, the action potential amplitude measured in front of a node of Ranvier is two times larger than in the middle of two nodes. Moreover, farther from the phantom axon, the measured action potential amplitude is almost constant along the longitudinal axis. The results of this new method confirm the computer simulations, with a correlation of 97.6 %. [ABSTRACT FROM AUTHOR]

This paper presents a novel hybrid medical stent device. This hybrid stent device formed by fractal mesh structures provides a flow-diverting effect and stent-assisted coil embolization. Flow-diverter stents decrease blood flow into an aneurysm to prevent its rupture. In general, the mesh size of a flow-diverter stent needs to be small enough to prevent blood flow into the aneurysm. Conventional flow-diverter stents are not available for stent-assisted coil embolization, which is an effective method for aneurysm occlusion, because the mesh size is too small to insert a micro-catheter for coil embolization. The proposed hybrid stent device is capable of stent-assisted coil embolization while simultaneously providing a flow-diverting effect. The fractal stent device is composed of mesh structures with fine and rough mesh areas. The rough mesh area can be used to insert a micro-catheter for stent-assisted coil embolization. Flow-diverting effects of two fractal stent designs were composed to three commercially available stent designs. Flow-diverting effects were analyzed using computational fluid dynamics (CFD) analysis and particle image velocimetry (PIV) experiment. Based on the CFD and PIV results, the fractal stent devices reduce the flow velocity inside an aneurism just as much as the commercially available flow-diverting stents while allowing stent-assisted coil embolization. [ABSTRACT FROM AUTHOR]

This paper investigates the effect of developing arterial blockage at the abdominal aorta on the blood pressure waves at an externally accessible location suitable for invasive measurements such as the brachial and the femoral arteries. Arterial blockages are created surgically within the abdominal aorta of healthy Wistar rats to create narrowing resemblance conditions. Blood pressure is measured using a catheter inserted into the right femoral artery. Measurements are taken at the baseline healthy condition as well as at four different severities (20, 50, 80 and 100 %) of arterial blockage. In vivo and in vitro measurements of the lumen diameter and wall thickness are taken using magnetic resonance imaging and microscopic techniques, respectively. These data are used to validate a 3D computational fluid dynamics model which is developed to generalize the outcomes of this work and to determine the arterial stress and strain under the blockage conditions. This work indicates that an arterial blockage in excess of 20 % of the lumen diameter significantly influences the pressure wave and reduces the systolic blood pressure at the right femoral artery. High wall shear stresses and low circumferential strains are also generated at the blockage site. [ABSTRACT FROM AUTHOR]

With the advent of powerful computers and parallel processing including Grid technology, the use of Monte Carlo (MC) techniques for radiation transport simulation has become the most popular method for modeling radiological imaging systems and particularly X-ray computed tomography (CT). The stochastic nature of involved processes such as X-ray photons generation, interaction with matter and detection makes MC the ideal tool for accurate modeling. MC calculations can be used to assess the impact of different physical design parameters on overall scanner performance, clinical image quality and absorbed dose assessment in CT examinations, which can be difficult or even impossible to estimate by experimental measurements and theoretical analysis. Simulations can also be used to develop and assess correction methods and reconstruction algorithms aiming at improving image quality and quantitative procedures. This paper focuses mainly on recent developments and future trends in X-ray CT MC modeling tools and their areas of application. An overview of existing programs and their useful features will be given together with recent developments in the design of computational anthropomorphic models of the human anatomy. It should be noted that due to limited space, the references contained herein are for illustrative purposes and are not inclusive; no implication that those chosen are better than others not mentioned is intended. [ABSTRACT FROM AUTHOR]

This paper presents a three-dimensional (3D) reconstruction system of the human spine for the routine evaluation of musculoskeletal pathologies like idiopathic scoliosis. The main objective of this 3D reconstruction system is to offer a versatile and robust tool for the 3D analysis of spines in any healthcare centre with standard clinical setup using standard uncalibrated radiographic images. The novel system uses a self-calibration algorithm and a weak-perspective method to reconstruct the 3D coordinates of anatomical landmarks from bi-planar radiographic images of a patient's trunk. Additionally, a small planar object of known dimensions is proposed to warrant an accurately scaled model of the spine. In order to assess the validity of the 3D reconstructions yielded by the proposed system, a clinical study using 60 pairs of digitized X-rays of adolescents was conducted. The subject cohort in the study group was composed of 51 scoliotic and 9 non-scoliotic patients, with an average Cobb angle on the frontal plane of 25 degrees. For each case, a 3D reconstruction of the spine and pelvis was obtained with the previous system used at our hospital (which requires a positioning apparatus and a calibration jacket), and with the proposed method. Results show that 3D reconstructions obtained with the new system using uncalibrated X-ray images yield geometrically accurate models with insignificant differences for 2D and 3D clinical indexes commonly used in the evaluation of spinal deformities. This demonstrates the system to be a viable and accurate tool for clinical studies and biomechanical analysis purposes, with the added advantage of versatility to any clinical setup for routine follow-ups and surgical planning. [ABSTRACT FROM AUTHOR]

Analysis of three-dimensional (3D) images of human torsos for torso deformities such as scoliosis requires classifying torso distortion. Assessing torso distortion from 3D images is not trivial as actual torsos are non-symmetric and show an outstanding range of variations leading to high classification errors. As the degree of spinal deformity (and classification of torso shape) influences scoliosis treatment options, the development of more accurate classification procedures is desirable. This paper presents a technique for assessing torso shape and classifying scoliosis into mild, moderate and severe categories using two indices, 'twist' and 'bend', obtained from orthogonally transformed images of the complete torso surface called orthogonal maps. Four transforms (axial line, unfolded cylinder, enclosing cylinder and subtracting cylinder) were used. Blind tests on 361 computer models with known deformation parameter values show 100% classification accuracy. Tests on eight volunteers without scoliosis validated the system and tests on 22 torso images of volunteers with scoliosis showed up to 95.5% classification accuracy. In addition to classifying scoliosis, orthogonal maps present the entire torso in one view and are viable for use in scoliosis clinics for monitoring the progression of scoliosis. [ABSTRACT FROM AUTHOR]

In this paper, we describe a computational model dedicated to building an apnoea monitoring system for newborn babies. The proposed model is based on whole body plethysmography, which involves non-invasive measurement of lung ventilation indirectly from the pressure deflections generated when a subject breathes inside a chamber of fixed volume (Bert in C R Soc Biol Paris 20:22-23, 1868). The computational model simulates thermal and environmental flow conditions occurring in the neonate chamber, especially steady state flow with heat transfer and carbon dioxide (CO2) transport during the exhalation phase. This permits the variance of all critical parameters and the analysis of their effects on the distributions of interest. The main objective is to study thermal and air quality comfort conditions under which infants can be monitored for long-term periods. The method deploys computational fluid dynamics techniques and parametric modelling which, by allowing input parameters to be modulated, represent a more efficient and flexible analytical tool than previous experimental techniques. Simulation data reveal that the largest flow rates occur in areas near the openings with slight formation of air recirculation zones; temperature distribution shows signs of stratification, with higher temperatures than the supplied air, CO2 distribution presents acceptable air quality level and predicted mean vote index affords a relatively acceptable thermal comfort level. This analytical approach can be considered as innovative, and can find a new application in clinical infant apnoea monitoring in a way that allows determination of the optimal location for placing a sensor to detect respiration activity without any contact with the infant's body, and without any risk, in contrast to available whole body plethysmography techniques previously tested in infants (Fleming et al. in J Appl Physiol 55:1924-1931, 1983). [ABSTRACT FROM AUTHOR]

Uniaxial stress-strain data were obtained from in vitro experiments on 20 porcine livers for compressions, elongations and cycles of compression and then elongation. There were about 70 cylindrical samples, with diameter 7mm and varying height (4-11 mm). The combined compression and elongation test provide a unified framework for both compression and elongation for applications such as computer-aided surgical simulation. It enable the zero stress state of the experimental liver sample to be precisely determined. A new equation that combined both logarithmic and polynomial strain energy forms was proposed in modelling these experimental data. The assumption of incompressibility was justified from a preliminary Poisson's ratio for elongation and compression at 0.43+/-0.16 and 0.47+/-0.15, respectively. This equation provided a good fit for the observed mechanical properties of liver during compression-elongation cycles and for separate compressions or elongations. The root mean square errors were 91.92+/-17.43 Pa, 57.55+/-13.23 Pa and 29.78+/-17.67 Pa, respectively. In comparison with existing strain energy functions, this combined model was the better constitutive equation. Application of this theoretical model to small liver samples and other tissues demonstrated its suitability as the material model of choice for soft tissue. [ABSTRACT FROM AUTHOR]

This study focused on the application of real-time Kalman filters to biomechanical data and, in particular, the simulation environment used to compare the performance of modified and standard two-state Kalman filters when estimating displacement and velocity from noisy displacement data. The modification proposed in this paper was the numerical tachometer, augmented by a median smoother. The numerical tachometer integrated the derivative estimates from finite differences of noisy sampled data into the Kalman filter structure; the median smoother acted before differentiation, to protect from grossly erroneous measurements. The numerical tachometer allowed better fits to the simulated data than can be achieved without it: the root mean square errors decreased by 10% in the displacement domain and by 54% in the velocity domain, for sampling frequencies and signal contamination levels that were typical in human movement sciences. The sensitivity to errors in the modelling of the signal and noise characteristics was less than in the standard filter implementation. The use of the median smoother improved the robustness of the filtering algorithm against additive white Gaussian measurement noise and allowed the cancellation of isolated noise spikes. [ABSTRACT FROM AUTHOR]

In electrical impedance tomography (EIT), current patterns are injected into a subject and boundary voltages are measured to reconstruct a cross-sectional image of resistivity distribution. Static EIT image reconstruction requires a computer model of a subject, an efficient data-collection method and robust and fast reconstruction algorithms. The finite-element method is used as the computer model. The paper describes the finite-element analysis software package developed, including an interactive graphical mesh generator and fast algorithms for solving linear systems of equations using sparse-matrix and vector techniques. Various models of irregularly shaped subjects are developed using mesh-design tools, including automatic mesh generation and optimisation using the Delaunay algorithm. Even though the software package is customised for use in electrical impedance tomography, it can be used for other biomedical research areas, such as impedance cardiography, cardiac defibrillation and impedance pneumography. [ABSTRACT FROM AUTHOR]

The fetal electrocardiogram (fECG) contains important information regarding the health of the fetus. However, the fECG obtained noninvasively from the abdominal surface electrical recordings of a pregnant woman are dominated by strong interference from the maternal electrocardiogram (mECG). In this paper, based on the H(infinity) principle, two adaptive algorithms are proposed for the extraction of fECG from the trans-abdominal recordings of pregnant women. The motivation behind the application of H(infinity) techniques is the fact that they are robust with respect to model uncertainties and lack of statistical information regarding noise. The proposed algorithms are applied to simulated as well as real multichannel ECG recordings and their performances are compared to that of the well-known least-mean-square (LMS) adaptive algorithm. It is found that the proposed H(infinity) based algorithms perform superior to the LMS algorithm in extracting the fECG signal. [ABSTRACT FROM AUTHOR]

Surface-recorded somatosensory-evoked potentials (SEPs) are neural signals elicited by an external stimulus. In the case of electrically induced SEPs, the artifact generated by the stimulation process can severely distort the signal. In some cases, the artifact tail often lasts well into the initiation of the SEP making the determination of absolute latency very difficult. In this work, a new approach was taken to identify factors that affect the tail of the artifact. The methodology adopted was the development of a lumped electrical circuit model of the artifact generation process. While the modeling of the instrumentation hardware is relatively simple, this is not the case with tissue and electrode/skin interface effects. Consequently, this paper describes a novel tissue modeling approach that uses an autoregressive moving average (ARMA) parametric technique and an artificial neural network (ANN) to estimate tissue parameters from experimental data. This coupled with an estimation of the stimulation electrode-skin impedance completes the lumped circuit model. Simulink (The Mathworks Inc.) was used to evaluate the model under several different conditions. These results show that both the stimulation electrode-skin interface impedance and nature of the body tissue directly under the recording electrodes have a profound effect on the appearance of the stimulus artifact tail. This was verified by experimentally recorded data obtained from the median nerve using surface electrodes. Conclusions drawn from this work include that stimulation electrodes with low series capacitance should be used whenever possible to minimize the duration of the artifact tail. [ABSTRACT FROM AUTHOR]

The electrogastrogram (EGG), a cutaneous measurement of gastric electrical activity, is a mixture of gastric slow waves and noise. To detect the propagation of gastric slow waves, it is desired to obtain gastric slow waves in each of multichannel EGGs. Recently, independent component analysis (ICA) has shown its efficiency in separating the gastric slow wave from noisy multichannel EGGs. However, this method is not able to recover gastric slow waves in each of the multichannel EGGs. In this paper, a two-stage combined method was proposed for extracting gastric slow waves. First, ICA was performed to separate the gastric slow wave component from noisy multichannel EGGs. Second, adaptive signal enhancement with a reference input derived by the ICA in the first stage was employed to extract gastric slow waves in each channel. Quantitative analysis showed that, with the proposed method, the maximum root-mean-square error between the estimated time lag and its theoretical value in the simulations was only 0.65. The results from real EGG data demonstrated that the combined method was able to extract gastric slow waves from individual channels of EGGs which are important to identify the slow wave propagation. Therefore, the proposed method can be used to detect propagation of gastric slow waves from multichannel EGGs. [ABSTRACT FROM AUTHOR]

To achieve an accurate estimation of the instantaneous turbulent velocity fluctuations downstream of prosthetic heart valves in vivo, the variability of the spectral method used to measure the mean frequency shift of the Doppler signal (i.e. the Doppler velocity) should be minimised. This paper investigates the performance of various short-time spectral parametric methods such as the short-time Fourier transform, autoregressive modelling based on two different approaches, autoregressive moving average modelling based on the Steiglitz-McBride method, and Prony's spectral method. A simulated Doppler signal was used to evaluate the performance of the above mentioned spectral methods and Gaussian noise was added to obtain a set of signals with various signal-to-noise ratios. Two different parameters were used to evaluate the performance of each method in terms of variability and accurate matching of the theoretical Doppler mean instantaneous frequency variation within the cardiac cycle. Results show that autoregressive modelling outperforms the other investigated spectral techniques for window lengths varying between 1 and 10 ms. Among the autoregressive algorithms implemented, it is shown that the maximum entropy method based on a block data processing technique gives the best results for a signal-to-noise ratio of 20 dB. However, at 10 and 0 dB, the Levinson-Durbin algorithm surpasses the performance of the maximum entropy method. It is expected that the intrinsic variance of the spectral methods can be an important source of error for the estimation of the turbulence intensity. The range of this error varies from 0.38% to 24% depending on the parameters of the spectral method and the signal-to-noise ratio. [ABSTRACT FROM AUTHOR]

A computer model is used to calculate the optimum geometry of an epidural electrode, consisting of a longitudinal contact array, for spinal cord stimulation in the management of chronic, intractable pain. 3D models of the spinal area are used for the computation of stimulation induced fields, and a cable model of myelinated nerve fibre is used for the calculation of the threshold stimulus to excite large dorsal column and dorsal root fibres. The criteria for the geometry of the longitudinal contact array are: a low threshold for the stimulation of dorsal column fibres compared with dorsal root fibres; and a low stimulation voltage (and current). For both percutaneous and laminectomy electrodes, the contact length should be approximately 1.5 mm, and the optimum contact separation, as determined by the computer model, is 2-2.5 mm. The contacts for a laminectomy electrode should be approximately 4 mm wide. This electrode geometry is applicable to all spinal levels where the dorsal columns can be stimulated (C1-2 down to L1). The stimulating electrode should preferably be used as a tripole with one (central) cathode. [ABSTRACT FROM AUTHOR]

Alcohol consumption may result in electrocardiographic changes and arrhythmias, at least partly due to effects of ethanol on cardiac ionic currents. Contractility and intracellular Ca(2+) dynamics seem to be altered as well. In this study, we integrated the available (mostly animal) experimental data into previously published models of the rat and human ventricular myocytes to assess the share of ionic current components in ethanol-induced changes in AP configuration and cytosolic Ca(2+) transient in ventricular cardiomyocytes. The rat model reproduced well the experimentally observed changes in AP duration (APD) under ethanol (slight prolongation at 0.8 mM and shortening at ≥8 mM). These changes were almost exclusively caused by the ethanol-induced alterations of I K1. The cytosolic Ca(2+) transient decreased gradually with the increasing ethanol concentration as a result of the ethanol-induced inhibition of I Ca. In the human model, ethanol produced a dose-dependent APD lengthening, dominated by ethanol effect on I Kr, the key repolarising current in human ventricles. This effect might contribute to the clinically observed proarrhythmic effects of ethanol in predisposed individuals. [ABSTRACT FROM AUTHOR]

We estimated the dynamic cell metabolic activity and the distribution of the pH value and oxygen concentration in tissue samples cultured in vitro by using real-time sensor records and a numerical simulation of the underlying reaction-diffusion processes. As an experimental tissue model, we used chicken spleen slices. A finite element method model representing the biochemical processes and including the relevant sensor data was set up. By fitting the calculated results to the measured data, we derived the spatiotemporal values of the pH value, the oxygen concentration and the absolute metabolic activity (extracellular acidification and oxygen uptake rate) of the samples. Notably, the location of the samples in relation to the sensors has a great influence on the detectable metabolic rates. The long-term vitality of the tissue samples strongly depends on their size. We further discuss the benefits and limitations of the model. [ABSTRACT FROM AUTHOR]

The clinical advantages of reducing spondylolisthesis over fusion in situ have several intuitive reasons such as restore the spinal column into a more anatomic relationship and alignment. However, there is only little evidence in the literature supporting the theoretical advantages of reduction, and its effect on spinopelvic alignment remains poorly defined. In this study, a comprehensive finite element model was developed to analyze the biomechanics of the spine after spinal fusion at L5-S1 in both types of high-grade spondylolisthesis (balanced and unbalanced pelvis). The relevant clinical indices (i.e. spondylolisthesis grade and Dubousset lumbosacral angle), the displacement of L4-L5, pressure within the annulus and nucleus, and stress at L4-L5 were evaluated and compared. The model can well predict the changes of the important clinical indices during the surgery. For a balanced pelvis, the reduction has a minimal effect on the biomechanical conditions at the adjacent level during postsurgical activities. In the unbalanced case, reduction induced larger deformation in the lumbosacral region and a higher stress concentration at adjacent level. Whether such a stress concentration can lead to long-term disc degeneration is not known. The results provide additional information for the clinician considering reduction of high-grade spondylolisthesis. [ABSTRACT FROM AUTHOR]

Since the advent of electromyogram recording, precise measures of tremor and gait have been used to study movement disorders such as Parkinson's disease. Now, a wide range of accelerometers and other motion-tracking technologies exist to better inform researchers and clinicians, yet such systems are rarely tested for accuracy or suitability before use. Our inexpensive test-rig can produce sinusoidal displacements using a simple cantilever system driven by a subwoofer. Controlled sinusoids were generated using computer software, and the displacement amplitudes of the test-rig were verified with fiducial marker tracking. To illustrate the use of the test-rig, we evaluated an accelerometer and an electromagnetic motion tracker. Accelerometry recordings were accurate to within ±0.09 g of actual peak-to-peak amplitude with a frequency response close to unity gain between 1 and 20 Hz. The electromagnetic sensor underestimated peak displacement by 2.68 mm, which was largely due to a diminishing gain with increasing frequency. Both sensors had low distortion. Overall sensitivity was limited by noise for the accelerometer and quantisation resolution for the electromagnetic sensor. Our simple and low-cost test-rig can be used to bench-test sensors used in movement disorders research. It was able to produce reliable sinusoidal displacements and worked across the 1- to 20-Hz frequency range. [ABSTRACT FROM AUTHOR]

In coronary artery bypass graft (CABG) surgery the involved tissues are overstretched, which may lead to intimal hyperplasia and graft failure. We propose a computational methodology for the simulation of traditional CABG surgery, and analyze the effect of two clinically relevant parameters on the artery and graft responses, i.e., incision length and insertion angle for a given graft diameter. The computational structural analyses are based on actual three-dimensional vessel dimensions of a human coronary artery and a human saphenous vein. The analyses consider the structure of the end-to-side anastomosis, the residual stresses and the typical anisotropic and nonlinear vessel behaviors. The coronary artery is modeled as a three-layer thick-walled tube. The finite element method is employed to predict deformation and stress distribution at various stages of CABG surgery. Small variations of the arterial incision have relatively big effects on the size of the arterial opening, which depends solely on the residual stress state. The incision length has a critical influence on the graft shape and the stress in the graft wall. Stresses at the heel region are higher than those at the toe region. The changes in the mechanical environment are severe along all transitions between the venous tissue and the host artery. Particular stress concentrations occur at the incision ends. The proposed computational methodology may be useful in designing a coronary anastomotic device for reducing surgical trauma. It may improve the quantitative knowledge of vessel diseases and serve as a tool for virtual planning of vascular surgery. [ABSTRACT FROM AUTHOR]

Pressure-related deep tissue injury is the term recommended by the United States National Pressure Ulcer Advisory Panel to describe a potentially life-threatening form of pressure ulcers, characterized by the presence of necrotic tissue under intact skin, and associated with prolonged compression of muscle tissue under bony prominences. In this study, a theoretical model was used to determine the relative contributions of the backrest inclination angle during prolonged wheelchair sitting, the muscle tissue stiffness and curvature of the ischial tuberosities (ITs) to the risk for injury in the gluteus muscles that pad the IT bones during sitting. The model is based on Hertz’s theory for analysis of contact pressures between a rigid half-sphere (bone) and an elastic half-space (muscle). Hertz’s theory is coupled with an injury threshold and damage law for muscle—both obtained in previous studies in rats. The simulation outputs the time-dependent bone–muscle contact pressures and the injured area in the gluteus. We calculated the full-size (asymptotic) injured area in the gluteus and the time for injury onset for different sitting angles α( (90–150°), muscle tissue long-term shear moduli G (250–1,200 Pa) and bone diameters D (8–18 mm). We then evaluated the sensitivity of model results to variations in these parameters, in order to determine how injury predictions are affected. In reclined sitting ( α = 150°) the full-size injured area was ∼2.1-fold smaller and the time for injury onset was ∼1.3-fold longer compared with erect sitting ( α = 90°). For greater G the full-size injured area was smaller but the time for injury onset was shorter, e.g., increasing G from 250 to 1200 Pa decreased the full-size injured area ∼2.5-fold, but shortened the time for injury onset 6.2-fold. For smaller D the time for injury onset dropped, e.g., decreased ∼1.5-fold when D decreased from 18 to 8 mm. Interestingly, the full-size injured area maximized at D of about 12 mm but decreased for smaller or larger D. The susceptibility to sitting-acquired deep tissue injury strongly depends on the geometrical and biomechanical characteristics of the bone–muscle interface, and, particularly, on the radius of curvature of the IT which mostly influenced the size of the wound, and on the muscle stiffness which dominantly affected the time for injury onset. [ABSTRACT FROM AUTHOR]

To minimize time-consuming and expensive in vitro and in vivo testing, information regarding the effects of implantation and the implants on measurements should be available during the designing of active implantable devices measuring bioelectric signals such as electrocardiograms (ECG). Modeling offers a fairly inexpensive and effective means of studying and demonstrating the effects of implantation on ECG measurements prior to any in vivo tests, and can thus provide the designer with valuable information. Finite difference model (FDM) and lead field approaches offer straightforward and effective modeling methods supporting the designing of active implantable ECG devices. The present study demonstrates such methods in developing and studying ECG implants. They were applied in demonstrating the effects of implant dimensions and of electrode implantation on the measurement sensitivity of the ECG device. The results of the simulations indicated that the interelectrode distance is the factor of the implant design determining the lead sensitivity. Other parameters related implant dimensions and shape have minor effect on the morphology of the ECG or on the average sensitivity of the measurement. This is shown for example when the interelectrode distance was reduced to 1/3 of original the average lead sensitivity decreased by 69.1% while larger relative changes in other dimensions produced clearly smaller changes. It was also observed here that implanting the electrodes deeper under the skin has major effects on the local sensitivities in heart muscle and thus affect to the morphology of the ECG. The study indicated also that non-conducting medium (i.e. implant insulated body) between the electrodes increases the sensitivity on heart muscle compared to cases where only electrodes are implanted. [ABSTRACT FROM AUTHOR]

A general framework for designing an optimum control strategy for the Hemopump is described. An objective function was defined that includes four membership functions, each constructed based on the desired values of one of the four members: stroke volume, mean left atrial pressure, aortic diastolic pressure and mean pump rotation speed. The Hemopump was allowed to operate either at a constant speed or at two different speeds during a cardiac cycle. The goal was to maximise the objective function by varying the magnitude and timing of the pump speed. Using a canine circulatory model, it was demonstrated that, in general, different cardiac conditions or different clinical objectives require different operation parameters. For example, when a left ventricle with minor ischaemia was simulated, and the main objective was to increase stoke volume, the objective function was maximised, from a value of 0.877 when the pump was off, to 0.946 when the pump was operated at speed 2 (18 500 revolutions min(-1)). On the other hand, for a severely ischaemic heart, the optimum pump speed became speed 3 (20 000 revolutions min(-1)), which maximized the objective function to 0.943 (from 0.707 when the pump was off). The results also suggest that it is more beneficial to operate the Hemopump at two different speeds during a cardiac cycle (a higher speed during systole and early diastole, and a lower speed during late diastole) than to maintain a constant speed throughout the cardiac cycle. [ABSTRACT FROM AUTHOR]

The proposed strap-down integration method exploits the cyclical nature of human gait: during the gait swing phase, the quaternion-based attitude representation is integrated using a gyroscope from initial conditions that are determined during stance by an accelerometer. Positioning requires double time integration of the gravity-compensated accelerometer signals during swing. An interpolation technique applied to attitude quaternions was developed to improve the accuracy of orientation and positioning estimates by accounting for the effect of sensor bias and scale factor drifts. A simulation environment was developed for the analysis and testing of the proposed algorithm on a synthetic movement trajectory. The aim was to define the true attitude and positioning used in the computation of estimation errors. By thermal modelling, the changes of bias and scale factor of the inertial sensors, calibrated at a single reference temperature, were analysed over a range of +/- 10 degrees C, for measurement noise standard deviations up to sigma(g) = 2.5 degrees s(-1) (gyroscope) and sigma(a) = 0.05 m s(-2) (accelerometer). The compensation technique reduced the maximum root mean square errors (RMSEs) to: RMSE(theta) = 14.6 degrees (orientation) and RMSE(d) = 17.7cm (positioning) for an integration interval of one gait cycle (an improvement of 3 degrees and 7 cm); RMSE(theta) = 14.8 degrees and RMSE(d) = 30.0 cm for an integration interval of two gait cycles (an improvement of 11 degrees and 262 cm). [ABSTRACT FROM AUTHOR]

Neuroprosthestic systems can be used to restore hand grasp and wrist control in individuals with C5/C6 spinal cord injury. A computer-based system was developed for the implementation, tuning and clinical assessment of neuroprosthetic controllers, using off-the-shelf hardware and software. The computer system turned a Pentium III PC running Windows NT into a non-dedicated, real-time system for the control of neuroprostheses. Software execution (written using the high-level programming languages LabVIEW and MATLAB) was divided into two phases: training and real-time control. During the training phase, the computer system collected input/output data by stimulating the muscles and measuring the muscle outputs in real-time, analysed the recorded data, generated a set of training data and trained an artificial neural network (ANN)-based controller. During real-time control, the computer system stimulated the muscles using stimulus pulsewidths predicted by the ANN controller in response to a sampled input from an external command source, to provide independent control of hand grasp and wrist posture. System timing was stable, reliable and capable of providing muscle stimulation at frequencies up to 24Hz. To demonstrate the application of the test-bed, an ANN-based controller was implemented with three inputs and two independent channels of stimulation. The ANN controller's ability to control hand grasp and wrist angle independently was assessed by quantitative comparison of the outputs of the stimulated muscles with a set of desired grasp or wrist postures determined by the command signal. Controller performance results were mixed, but the platform provided the tools to implement and assess future controller designs. [ABSTRACT FROM AUTHOR]

The objective of this research was to develop a purely biomechanical model, intended to predict the long-term secondary stability of the implant starting from the biomechanical stability immediately after the operation. A continuous rule-based adaptation scheme was formulated as a dynamic system, and the work verified if such a model produced unique and clinically meaningful solutions. It also investigated whether this continuous model provided results comparable with those of a simpler, discrete-states model used in a previous study. The proposed model showed stable convergence behaviour with all investigated initial conditions, with oscillatory behaviour limited to the first steps of the simulation. The results obtained with the wide range of initial conditions support the hypothesis of the existence and uniqueness of the solution for all initial conditions. The differences between the continuous model and the simpler and more efficient finite-states model were found to be extremely modest (less than 4% over the predicted bonded area). Because of these minimal differences, the use of the much faster finite-states model is recommended to investigate asymptotic conditions, and the continuous model described should be used to investigate the evolution over time of the adaptive process. [ABSTRACT FROM AUTHOR]

A computational model of the fluid dynamics of intraventricular flow was used to investigate the importance of the effects of flow disturbances existing within the left ventricle (LV) at the onset of diastole on a diastolic flow field. The simulation started with a quiescent flow state; it continued for a number of cardiac cycles to obtain a cyclically repeatable flow. After the flow became periodic, the initial diastolic flow was not quiescent: flow disturbances, remnants of a systolic flow, were present within the LV. Nevertheless, they faded away during an acceleration phase of diastole and almost ceased by the end of this phase. Consequently, a flow field during a deceleration phase of diastole, characterised by the formation of a vortex ring, was hardly affected by the initial flow disturbances. The propagation velocity of a colour M-mode Doppler echocardiogram obtained by scanning velocity along the LV long axis was 0.58 m s(-1) in the case where diastolic flow was initially quiescent and 0.56 m s(-1) in the case where flow disturbances existed at the beginning of diastole. These results indicated that the colour M-mode Doppler echocardiographic technique captures flow dynamics produced purely by ventricular expansion, with little influence from initial diastolic flow disturbances. [ABSTRACT FROM AUTHOR]

High-frequency alternating current (AC) waveforms have been shown to produce a quickly reversible nerve block in animal models, but the parameters and mechanism of this block are not well understood. A frog sciatic nerve/gastrocnemius muscle preparation was used to examine the parameters for nerve conduction block in vivo, and a computer simulation of the nerve membrane was used to identify the mechanism for block. The results indicated that a 100% block of motor activity can be accomplished with a variety of waveform parameters, including sinusoidal and rectangular waveforms at frequencies from 2 kHz to 20 kHz. A complete and reversible block was achieved in 34 out of 34 nerve preparations tested. The most efficient waveform for conduction block was a 3-5 kHz constant-current biphasic sinusoid, where block could be achieved with stimulus levels as low as 0.01 microCphase(-1). It was demonstrated that the block was not produced indirectly through fatigue. Computer simulation of high-frequency AC demonstrated a steady-state depolarisation of the nerve membrane, and it is hypothesised that the conduction block was due to this tonic depolarisation. The precise relationship between the steady-state depolarisation and the conduction block requires further analysis. The results of this study demonstrated that high-frequency AC can be used to produce a fast-acting, and quickly reversible nerve conduction block that may have multiple applications in the treatment of unwanted neural activity. [ABSTRACT FROM AUTHOR]

A head-neck computer model was comprehensively validated over a range of rear-impact velocities using experiments conducted by the same group of authors in the same laboratory. Validations were based on mean +/- 1 standard deviation response curves, i.e. corridors. Global head-neck angle, segmental angle and local facet joint regional kinematic responses from the model fell within experimental corridors. This was true for all impact velocities (1.3, 1.8 and 2.6 m s(-1)). The non-physiological S-curvature lasted approximately 100 ms. The present, comprehensively validated model can be used to conduct parametric studies and investigate the effects of factors such as active sequential and parallel muscle contractions, thoracic ramping and local tissue strain responses, as a function of cervical level, joint region and impact velocity in whiplash injury assessment. [ABSTRACT FROM AUTHOR]

As part of the development of new modelling tools for the simulation and design of brace treatment of scoliosis, a finite element model of a brace and its interface with the torso was proposed. The model was adapted to represent one scoliotic adolescent girl treated with a Boston brace. The 3D geometry was acquired using multiview radiographs. The model included the osseo-ligamentous structures, thoracic and abdominal soft tissues, brace foam and shell, and brace-torso interface. The simulations consisted of brace opening to include the patient's trunk followed by brace closing. To validate the model, the resulting geometry was compared with the real in-brace geometry, and the resulting contact reaction forces at the brace-torso interface were compared with the equivalent forces calculated from pressure measurements made on the in-brace patient. Differences between coronal equivalent and reaction forces were less than 7N. However, sagittal reaction forces (47N) were computed on the abdomen, whereas negligible equivalent forces were measured. The simulated geometry presented partially reduced coronal Cobb angles (1-4 degrees), over-corrected sagittal Cobb angles and maximum deformation plane (5 degrees), completely corrected coronal shift, and sagittal shift and rib humps that were not corrected. This study demonstrated the feasibility of a new approach that represents the load transfer from the brace to the spine more realistically than does the direct application of forces. [ABSTRACT FROM AUTHOR]

In the study, a new myocardial infarction (MI) estimation method was developed for estimating MI in the three-dimensional myocardium by means of a heart-model-based inverse approach. The site and size of MI are estimated from body surface electrocardiograms by minimising multiple objective functions of the measured body surface potential maps (BSPMs) and the heart-model-generated BSPMs. Computer simulations were conducted to evaluate the performance of the developed method, using a single-site MI and dual-site MI protocols. The simulation results show that, for the single-site MI, the averaged spatial distance (SD) between the weighting centres of the "true" and estimated MIs, and the averaged relative error (RE) between the numbers of the "true" and estimated infarcted units are 3.0 ± 0.6/3.6 ± 0.6mm and 0.11 ± 0.02/0.14 ± 0.02, respectively, when 5µV/10µV Gaussian white noise was added to the body surface potentials. For the dual-site MI, the averaged SD between the weighting centres of the 'true' and estimated MIs, and the averaged RE between the numbers of the 'true' and estimated infarcted units are 3.8 ± 0.7/3.9 ± 0.7mm and 0.12 ± 0.02/0.14 ± 0.03, respectively, when 5µV/10 µV Gaussian white noise was added to the body surface potentials. The simulation results suggest the feasibility of applying the heart-model-based imaging approach to the estimation of myocardial infarction from body surface potentials. [ABSTRACT FROM AUTHOR]

Radio-frequency thermokeratoplasty (RF-TKP) is a technique used to reshape the cornea curvature by means of thermal lesions using radio-frequency currents. This curvature change allows refractive disorders such as hyperopia to be corrected. A new electrode with ring geometry is proposed for RF-TKP. It was designed to create a single thermal lesion with a full-circle shape. Finite element models were developed, and the temperature distributions in the cornea were analysed for different ring electrode characteristics. The computer results indicated that the maximum temperature in the cornea was located in the vicinity of the ring electrode outer perimeter, and that the lesions had a semi-torus shape. The results also indicated that the electrode thickness, electrode radius and electrode thermal conductivity had a significant influence on the temperature distributions. In addition, in vitro experiments were performed on rabbit eyes. At 5 W power, the lesions were fully circular. Some lesions showed non-uniform characteristics along their circular path. Lesion depth depended on heating duration (60% of corneal thickness for 20 s, and 30% for 10 s). The results suggest that the critical shrinkage temperature (55-63 degrees C) was reached at the central stroma and along the entire circular path in all the cases. [ABSTRACT FROM AUTHOR]

Various investigators have used the monodomain model to study cardiac propagation behaviour. In many cases, the governing non-linear parabolic equation is solved using the finite-difference method. An adequate discretisation of cardiac tissue with realistic dimensions, however, often leads to a large model size that is computationally demanding. Recently, it has been demonstrated, for a two-dimensional homogeneous monodomain, that the Chebyshev pseudospectral method can offer higher computational efficiency than the finite-difference technique. Here, an extension of the pseudospectral approach to a three-dimensional inhomogeneous case with fibre rotation is presented. The unknown transmembrane potential is expanded in terms of Chebyshev polynomial trial functions, and the monodomain equation is enforced at the Gauss-Lobatto node points. The forward Euler technique is used to advance the solution in time. Numerical results are presented that demonstrate that the Chebyshev pseudospectral method offered an even larger improvement in computational performance over the finite-difference method in the three-dimensional case. Specifically, the pseudospectral method allowed the number of nodes to be reduced by approximately 85 times, while the same solution accuracy was maintained. Depending on the model size, simulations were performed with approximately 18-41 times less memory and approximately 99-169 times less CPU time. [ABSTRACT FROM AUTHOR]

Investigations were carried out into whether an equivalent generator consisting of two dipoles could be used to detect dual sites of ventricular activity. A computer model of the human ventricular myocardium was used to simulate activation sequences initiated at eight different pairs of sites positioned on the epicardial surface of the atrio-ventricular ring. From these sequences, 117-lead body surface potentials (covering the anterior and posterior torso), 64-lead magnetic field maps (above the anterior chest) and 128-lead magnetic field maps (above the anterior and posterior chest) were simulated and were then used to localise dual accessory pathways employing pairs of equivalent dipoles. Average localisation errors were 12 mm, 12 mm and 9 mm, respectively, when body surface potentials, 64-lead and 128-lead magnetic fields were used. The results of the study suggest that solving the inverse problem for two dipoles could provide additional information on dual accessory pathways prior to electrophysiological study. [ABSTRACT FROM AUTHOR]

In the present study, we report body surface Laplacian mapping of atrial depolarization under sinus rhythm in 8 healthy male subjects. For each subject, 95 unipolar disk electrodes with inter-electrode distance of 2 cm were used to record simultaneously potential ECGs over the anterior chest. The Laplacian ECG was then estimated during the P wave using a novel spline Laplacian technique. The body surface potential map (BSPM) and body surface Laplacian map (BSLM) at different time instants or time intervals of the P wave were constructed and compared. The present results showed that the BSPMs during the P wave were characterized by the rotation of a pair of positive/negative potential distribution from right to left around the anterior torso. On the other hand, the corresponding BSLMs revealed more spatial details, including two positive activities (denoted as P1 and P2, appeared in all 8 subjects), and three negative activities (denoted as N1, N2, and N3, appeared in 7, 7, and 4 subjects, respectively). The separation of these activities and their evolving patterns were also compared and confirmed by computer simulation using a realistic geometry heart-torso model. The above findings may be directly related to the underlying activation sequence during atrial depolarization in healthy subjects, suggesting the potential clinical applications of the Laplacian ECG technique. [ABSTRACT FROM AUTHOR]

Transcranial electrical stimulation of the brain causes considerable discomfort to the patient. The purpose of the study was to find out whether this could be affected by the choice of stimulation parameters. A spherical volume conductor model of the head and active compartmental models of a pyramidal motor nerve and scalp nociceptor were used in combination to simulate the scalp nociception to transcranial electrical stimulation. Scalp nociceptors were excited at distances of several centimetres from the electrodes. The size of the excited scalp area correlated with the length of the stimulation pulse. The area was 12.3, 20.4 and 26.0 cm2, for a 10 micros, 100 micros and 1 ms constant current pulse, respectively. With a 100 micros constant current pulse, the threshold for motor excitation was 205mA and, for nociception, it was 51 mA. There was no significant difference between constant current and capacitor discharge pulses or between electrodes of different sizes. The results imply that the use of very short stimulation pulses can reduce the pain. If a topical anaesthesia is used to reduce the pain, it has to be applied on a large area around the electrodes. [ABSTRACT FROM AUTHOR]

Failure of an orthopaedic fixation due to stress shielding and consequent screw loosening is a major concern among surgeons: the loosened screws could not only interfere with the healing process but also endanger adjacent anatomical structures. In this study, the effect of the screw's engineering design (dimensions, profile shape and material properties) on the load sharing with adjacent bone and consequent bone resorption was tested, using a set of two-dimensional computational (finite element) models. An algorithm simulating local bone adaptation to strain energy density (SED) mechanical stimuli was developed and used to evaluate the biomechanical performances of different commercial screws. Two new designs, a 'graded-stiffness' composite screw, with a reduced-stiffness titanium core and outer polymeric threads, and an active-compression hollow screw that generates compressive stresses on the surrounding bone, were also evaluated. A dimensionless set of stress transfer parameters (STPs) were utilised for ranking the performances of the different screws according to the expected screw-bone load sharing and its evolution with adaptation of the surrounding tissue. The results indicated that commercial wide (6 mm thread diameter) trapezoidal and rectangular screw profiles have superior biomechanical compatibility with bone (i.e. predicted to be stable after 2 years). The graded-stiffness and active-compression screws provided the best biomechanical performances: bone loading around them was predicted to decrease by no more than 15% after 3 years, compared with a decrease of 55-70% in bone loading around commercially available screws. Computer simulations of bone adaptation around orthopaedic screws are demonstrated to be effective means for objective and quantitative evaluation of the biomechanical aspects of implant-tissue compatibility. [ABSTRACT FROM AUTHOR]

Muscle fibre conduction velocity (MFCV) can be used as an index of the structural and/or functional modifications that can occur during fatigue or pathological processes. Current evaluation of MFCV from surface electromyography (SEMG) classically produces an average value. However, a single mean value is not sufficient when modifications affect only a small part of the conduction velocity distribution. In such a case, an estimation of the whole motor unit conduction velocity distribution (MUCV) would be advantageous. The aim of this study was the evaluation of the quality of two short-term methods based on cross-correlation (CC) and peak-to-peak (PP) estimation. A comprehensive simulation program was used to generate signals with known MUCV distributions. The Dmax statistic of Kolmogorov-Smirnov was used as an error criterion to quantify the estimation error and to optimise the MUCV distribution computation algorithms. The minimum error was observed for an analysing window of 10ms for PP and 15ms for CC. Dmax was significantly lower for PP (0.195+/-0.054) than for CC (0.343+/-0.073). Various simulations showed the strong effect of the variance of the true distribution on the features of the estimated ones. Clinical data measured on the abductor pollicis brevis were studied. MUCV was estimated on a healthy subject (3.63+/-0.87ms(-1)), a patient suffering from a myopathy (2.73+/-0.51ms(-1)) and one suffering from a neuropathy (4.38+/-0.23ms(-1)). The results demonstrate the overall superiority of a peak-to-peak approach. [ABSTRACT FROM AUTHOR]

Although the single equivalent point dipole model has been used to represent well-localised bio-electrical sources, in realistic situations the source is distributed. Consequently, position estimates of point dipoles determined by inverse algorithms suffer from systematic error due to the non-exact applicability of the inverse model. In realistic situations, this systematic error cannot be avoided, a limitation that is independent of the complexity of the torso model used. This study quantitatively investigates the intrinsic limitations in the assignment of a location to the equivalent dipole due to distributed electrical source. To simulate arrhythmic activity in the heart, a model of a wave of depolarisation spreading from a focal source over the surface of a spherical shell is used. The activity is represented by a sequence of concentric belt sources (obtained by slicing the shell with a sequence of parallel plane pairs), with constant dipole moment per unit length (circumferentially) directed parallel to the propagation direction. The distributed source is represented by N dipoles at equal arc lengths along the belt. The sum of the dipole potentials is calculated at predefined electrode locations. The inverse problem involves finding a single equivalent point dipole that best reproduces the electrode potentials due to the distributed source. The inverse problem is implemented by minimising the χ2 per degree of freedom. It is found that the trajectory traced by the equivalent dipole is sensitive to the location of the spherical shell relative to the fixed electrodes. It is shown that this trajectory does not coincide with the sequence of geometrical centres of the consecutive belt sources. For distributed sources within a bounded spherical medium, displaced from the sphere's centre by 40% of the sphere's radius, it is found that the error in the equivalent dipole location varies from 3 to 20% for sources with size between 5 and 50% of the sphere's radius. Finally, a method is devised to obtain the size of the distributed source during the cardiac cycle. [ABSTRACT FROM AUTHOR]

A baroregulation model, based on a previous pulsatile non-linear multielement cardiovascular model, is extended and used to study short-term regulation mechanisms. Using this model, the responses of several cardiovascular variables to different exercise levels are simulated and compared with the experimental data reported in the literature. The impact of physiological or pathological changes on the short-term regulation of arterial pressure under the stimulus of moderate dynamic exercise is then studied. The simulation results indicate that baroreflex feedback plays a critical role in the short-term regulation of arterial pressure. When the baroreflex gain decreases to one-third of the normal value, the response of the mean arterial pressure to moderate dynamic exercise and post-exercise recovery time increases by factors of 1.7 and 2.3, respectively. Clinical data from 36 subjects (two groups: normal and hypertensive) are collected to validate the model. Computer simulations for the hypertensive group show that the elastic modulus of the arterial vessel wall is increased by 1.5 times, and peripheral resistance is increased by 1.3 times the normal value, and the baroreflex gain decreases from 0.55 (for the normal group) to 0.40. The simulation results for normal and hypertensive groups agree well with the clinical data. [ABSTRACT FROM AUTHOR]

The absolute quantified measurement of haemoglobin skin blood saturation from collected reflectance spectra of the skin is complicated by the fact that the blood content of tissues can vary both in the spatial distribution and in the amount. These measurements require an understanding of which vascular bed is primarily responsible for the detected signal. Knowing the spatial detector depth sensitivity makes it possible to find the best range of different probe geometries for the measurements of signal from the required zones and group of vessels inside the skin. To facilitate this, a Monte Carlo simulation has been developed to estimate the sampling volume offered by fibre-optic probes with a small source-detector spacing (in the current report 250 microm, 400 microm and 800 microm). The optical properties of the modelled medium are taken to be the optical properties of the Caucasian type of skin tissue in the visible range of the spectrum. It is shown that, for a small source-detector separation (800 microm and smaller), rough boundaries between layers of different refractive index can play a significant role in skin optics. Wavy layer interfaces produce a deeper and more homogeneous distribution of photons within the skin and tend to suppress the direct channelling of photons from source to detector. The model predicts that a probe spacing of 250 microm samples primarily epidermal layers and papillary dermis, whereas spacings of 400-800 microm sample upper blood net dermis and dermis. [ABSTRACT FROM AUTHOR]

A modified adaptive Takagi-Sugeno (TS) fuzzy logic controller (FLC) is proposed that allows a simulated elbow-like biomechanical system to accurately track sigmoidal and sinusoidal trajectories in the sagittal plane. The work is a first effort towards the implementation of a system to restore elbow movements in quadriplegics using functional neuromuscular stimulation. The single-joint musculo-skeletal system is composed of a co-contractable pair of electrically stimulated muscles; the muscle model accounts for the increase in fatigue during the tracking exercise. In the proposed controller structure, a reinforcement learning scheme is used to accomplish the parameter tuning, and the parameter projection algorithm guarantees the system stability during the adaptation process. The controller performance is evaluated using computer simulation experiments and compared with the performance achievable when a standard proportional-integrative-derivative (PID) controller is employed for the same application. The modified adaptive TSFLC outperforms the PID controller in all tested situations, with a clear-cut advantage in the case of high-frequency sinusoidal trajectories (angular frequencies spanning the interval 8-12 rad s-1). The standard controller suffers from a dramatic increase in root mean square (RMS) tracking error above the value at 8 rad s-1, e.g. ERMS > or = 0.013, whereas the correlation coefficient between the actual and desired trajectory falls almost to zero, starting from the value rho approximately equal to 0.97 at 8 rad s-1. On the other hand, the adaptive TSFLC yields ERMS < or = 0.015, with rho > or = 0.78, over the whole range of tested angular frequencies. [ABSTRACT FROM AUTHOR]

As the craniospinal space is a pressure loaded system it is difficult to conceptualize and understand the flow dynamics through the ventricular system. Aqueduct stenosis compromises flow, increasing the pressure required to move cerebrospinal fluid (CSF) through the ventricles. Under normal circumstances, less than one pascal (1 Pa) of pressure is required to move a physiological flow of CSF through the aqueduct. This is too small to measure using clinical pressure transducers. A computational fluid dynamics (CFD) program, CFX, has been used to model two forms of aqueduct stenosis: simple narrowing and forking of the aqueduct. This study shows that with mild stenoses, the increase in pressure required to drive flow becomes significant (86-125 Pa), which may result in an increased transmantle pressure difference but not necessarily an increased intraventricular pressure. Severe stenoses will result in both. Wall shear stresses increase concomitantly and may contribute to local damage of the aqueduct wall and further gliosis with narrowing. [ABSTRACT FROM AUTHOR]

The effectiveness of substance delivery through catheters is an important issue in interventional radiology, especially for infusion therapies where the pharmacokinetic advantage of local intra-arterial drug administration has been firmly established. In principle, the procedure is used to provide appropriate local concentrations while maintaining low systemic values so as to minimise the global effect and toxicity of the intervention. However, poor drug mixing may produce excessive local concentrations potentially damaging for the surrounding tissues and may lead to unsuccessful therapies. These phenomena have been observed in the infusion therapies of liver cancers through the hepatic artery and with brain tumour therapies through the carotid artery. Many aspects of the drug delivery methodology have been explored in order to determine the infusion conditions that would provide optimal mixing: the catheter tip design is considered one of the most important characteristics to be investigated for this purpose. Interestingly, it turns out that angiographic procedures could also benefit from this, because better mixing properties are associated with designs that provide potentially less harmful flow conditions such as jets, whipping and recoil of the catheter on the vascular wall. A 2D steady numerical model is proposed, to simulate the main physical processes occurring during catheter substance infusion: blood dynamics is taken into account with the Navier-Stokes equations and substance dispersion by the flowing blood with the advection-diffusion equation. The model is used to evaluate mixing properties of certain catheter designs in different flow conditions. In particular, two types of side hole catheter are compared in the context of water bath injection and in the context of vessel injection. The simulations suggest that the improved mixing reported with water bath experiments would not be maintained in the clinical context of arterial circulation. [ABSTRACT FROM AUTHOR]

A three-layer Newtonian model is investigated using a combined Eulerian-Lagrangian computational method to describe the dynamic behaviour of leukocytes. The model, composed of a cell membrane (outer layer), cytoplasm (middle layer) and nucleus (inner layer), can better describe the recovery characteristics because large viscosity and capillarity differences between layers are considered, and both Newtonian and seemingly non-Newtonian behaviours reported in the literature can be reproduced. It is found that, to describe adequately the various rheological characteristics of leukocytes, the presence of the highly viscous nucleus and its deformation/recovery, as well as the surface energy stored in the fluid interfaces, are critical. Photographs from pipette experiments using a fluorescent technique confirm the theoretical finding of the important role played by the nucleus in cell deformation. [ABSTRACT FROM AUTHOR]

A model study introduces the use of statistical signal processing to analyse the signals from an intracavitary probe. A complete derivation is given for the detection of one type of arrhythmogenic substrate, myocardial infarctions (MIs). Both the use of statistical signal processing and the detection of VT substrates, as opposed to activation maps, are unique. A quasi-stationary electromagnetic model with simplified geometry is presented. The model is used to simulate ventricular pacing in the presence of MI. The likelihood ratio is used for detection. A tabulation of the results from this model shows that an intracavitary probe can be used to detect MIs as small as 400 mm2 in 1 mV of noise with a detectability index of 0.495, where 0.5 indicates perfect detection. Sensitivity to noise can be reduced by analysing multiple heart beats. The results are only slightly affected by changing the probe from a cage frame design, which mechanically supports the electrodes on thin spokes, to a balloon design, which supports the electrodes on the surface of an insulating balloon. [ABSTRACT FROM AUTHOR]