Simulations and visualizations for interpretation of brain microdialysis data during deep brain stimulation.

Microdialysis of the basal ganglia was used in parallel to deep brain stimulation (DBS) for patients with Parkinson's disease. The aim of this study was to patient-specifically simulate and visualize the maximum tissue volume of influence (TVI(max)) for each microdialysis catheter and the electric field generated around each DBS electrode. The finite element method (FEM) was used for the simulations. The method allowed mapping of the anatomical origin of the microdialysis data and the electric stimulation for each patient. It was seen that the sampling and stimulation targets differed among the patients, and the results will therefore be used in the future interpretation of the biochemical data.

A model for simulation and patient-specific visualization of the tissue volume of influence during brain microdialysis.

Microdialysis can be used in parallel to deep brain stimulation (DBS) to relate biochemical changes to the clinical outcome. The aim of the study was to use the finite element method to predict the tissue volume of influence (TVI(max)) and its cross-sectional radius (r (TVImax)) when using brain microdialysis, and visualize the TVI(max) in relation to patient anatomy. An equation based on Fick's law was used to simulate the TVI(max). Factorial design and regression analysis were used to investigate the impact of the diffusion coefficient, tortuosity and loss rate on the r (TVImax). A calf brain tissue experiment was performed to further evaluate these parameters. The model was implemented with pre-(MRI) and post-(CT) operative patient images for simulation of the TVI(max) for four patients undergoing microdialysis in parallel to DBS. Using physiologically relevant parameter values, the r (TVImax) for analytes with a diffusion coefficient D = 7.5 × 10⁻6 cm²/s was estimated to 0.85 ± 0.25 mm. The simulations showed agreement with experimental data. Due to an implanted gold thread, the catheter positions were visible in the post-operative images. The TVI(max) was visualized for each catheter. The biochemical changes could thereby be related to their anatomical origin, facilitating interpretation of results.

Impact of cysts during radiofrequency lesioning in deep brain structures--a simulation and in vitro study.

Radiofrequency lesioning of nuclei in the thalamus or the basal ganglia can be used to reduce symptoms caused by e.g. movement disorders such as Parkinson's disease. Enlarged cavities containing cerebrospinal fluid (CSF) are commonly present in the basal ganglia and tend to increase in size and number with age. Since the cavities have different electrical and thermal properties compared with brain tissue, it is likely that they can affect the lesioning process and thereby the treatment outcome. Computer simulations using the finite element method and in vitro experiments have been used to investigate the impact of cysts on lesions' size and shape. Simulations of the electric current and temperature distributions as well as convective movements have been conducted for various sizes, shapes and locations of the cysts as well as different target temperatures. Circulation of the CSF caused by the heating was found to spread heat effectively and the higher electric conductivity of the CSF increased heating of the cyst. These two effects were together able to greatly alter the resulting lesion size and shape when the cyst was in contact with the electrode tip. Similar results were obtained for the experiments.

Thermally induced convective movements in a standard experimental model for characterization of lesions prior to radiofrequency functional neurosurgery.

Experimental exploration of equipment for stereotactic functional neurosurgery based on heating induced by radio-frequency current is most often carried out prior to surgery in order to secure a correct function of the equipment. The experiments are normally conducted in an experimental model including an albumin solution in which the treatment electrode is submerged, followed by a heating session during which a protein clot is generated around the electrode tip. The clot is believed to reflect the lesion generated in the brain during treatment. It is thereby presupposed that both the thermal and electric properties of the model are similar to brain tissue. This study investigates the presence of convective movements in the albumin solution using laser Doppler velocimetry. The result clearly shows that convective movements that depend on the time dependent heating characteristics of the equipment arise in the solution upon heating. The convective movements detected show a clear discrepancy compared with the in vivo situation that the experimental model tries to mimic; both the velocity (maximum velocity of about 5 mms) and mass flux are greater in this experimental setting. Furthermore the flow geometry is completely different since only a small fraction of the tissue surrounding the electrode in vivo consists of moving blood, whereas the entire surrounding given by the albumin solution in the experimental model is moving. Earlier investigations by our group (Eriksson et al., 1999, Med. Biol. Eng. Comput. 37, pp. 737-741; Wren, 2001, Ph.D. thesis; and Wren et al., 2001, Med. Biol. Eng. Comput. 39, pp. 255-262) indicate that the heat flux is an essential parameter for the lesion growth and final size, and that presence of convective movements in the model might substantially increase the heat flux. Thus, convective movements of the magnitude presented here will very likely underestimate the size of the brain lesion, a finding that definitely should be taken into consideration when using the model prior to patient treatment.

Radio-frequency lesioning in brain tissue with coagulation-dependent thermal conductivity: modelling, simulation and analysis of parameter influence and interaction.

Radio-frequency brain lesioning is a method for reducing e.g. symptoms of movement disorders. A small electrode is used to thermally coagulate malfunctioning tissue. Influence on lesion size from thermal and electric conductivity of the tissue, microvascular perfusion and preset electrode temperature was investigated using a finite-element model. Perfusion was modelled as an increased thermal conductivity in non-coagulated tissue. The parameters were analysed using a 2(4)-factorial design (n=16) and quadratic regression analysis (n=47). Increased thermal conductivity of the tissue increased lesion volume, while increased perfusion decreased it since coagulation creates a thermally insulating layer due to the cessation of blood perfusion. These effects were strengthened with increased preset temperature. The electric conductivity had negligible effect. Simulations were found realistic compared to in vivo experimental lesions.

Feasibility of patient specific aortic blood flow CFD simulation.

Patient specific modelling of the blood flow through the human aorta is performed using computational fluid dynamics (CFD) and magnetic resonance imaging (MRI). Velocity patterns are compared between computer simulations and measurements. The workflow includes several steps: MRI measurement to obtain both geometry and velocity, an automatic levelset segmentation followed by meshing of the geometrical model and CFD setup to perform the simulations follwed by the actual simulations. The computational results agree well with the measured data.

Errors in body temperature assessment related to individual variation, measuring technique and equipment.

Errors in body temperature measurement might seriously influence the evaluation of an individual's health condition. We studied individual variation, measurement technique and the equipment used when assessing body temperature. In the first part of the study, three volunteers performed repeated measurements for five mornings. In the second part, the morning rectal, oral, ear and axillary temperatures were measured once in 84 men and women (19-59 years). The repeated measurements showed a daily temperature difference of 0.1-0.4 degrees C in rectal and oral temperatures, 0.2 degrees C-1.7 degrees C in the ear and 0.1-0.9 degrees C in the axillary temperatures. In the sample of 84 subjects, men and postmenopausal women had a lower mean body temperature compared to premenopausal women. The mean deviation between rectal temperature, and oral, ear and axillary temperatures, respectively, was > 0.5 degrees C, with a large individual variation. In conclusion, in order to improve the evaluation of body temperature, the assessment should be based on the individual variation, the same site of measurement and no adjustment of oral, ear or axillary temperatures to the rectal site.

Investigation of medical thermal treatment using a hybrid bio-heat model.

Bio-heat transfer, - heat transfer affecting living organism under the influence of blood perfusion -, is given great and increasing attention in medicine today. One reason is the increasing use of thermal treatment methods in for example heart- and neuro-surgery. Analysis and modelling of the thermal aspects is frequently carried out at every stage of device and method development, as it exhibits unique possibilities to understand the complex interactions present. This work investigates the use of a hybrid bio-heat model/equation, which is subsequently used to analyse temperature measurement during thermal treatment of the prostate.

Comparison between a detailed and a simplified finite element model of radio-frequency lesioning in the brain.

A detailed and a simplified model of a lesioning electrode was made using the finite element method. 15 simulations of the lesioning procedure were performed for each model and the resulting lesion volumes were compared in order to investigate if the simplified model is adequate. The simplified model resulted in a very slight overestimation of the volume compared to the detailed model. It was thus concluded that the simplified model is adequate for simulations.

Net pressure gradients in aortic prosthetic valves can be estimated by Doppler.

BACKGROUND: In aortic prosthetic valves, both the Doppler-estimated gradients and orifice areas are misleading in the assessment of hemodynamic performance. The parameter of major interest is the net pressure gradient after pressure recovery (PR). We, therefore, investigated, in vitro, our ability to predict the net pressure gradient and applied the formulas in a representative patient population with 2 different valve designs. METHODS: We studied the St Jude Medical (SJM) standard valve (size 19-27) and SJM Biocor (size 21-27) in an in vitro steady-flow model with simultaneous Doppler-estimated pressure and catheter pressure measurements. Using echocardiography, we also studied patients who received the SJM (n = 66) and SJM Biocor (n = 45). RESULTS: In the SJM, we observed PR both within the prosthesis and aorta, whereas in the SJM Biocor, PR was only present in the aorta. We estimated the PR within the valve and within the aorta separately from echocardiographic in vitro data, combining a regression equation (valve) with an equation on the basis of fluid mechanics theory (aorta). The difference between estimated and catheter-obtained net gradients (mean +/- SD) was 0.6 +/- 1.6 mm Hg in the SJM and -0.2 +/- 1.9 mm Hg in the SJM Biocor. When these equations were applied in vivo, we found that PR had an overall value of 57 +/- 7% of the peak Doppler gradient in the SJM and 33 +/- 9% in the SJM Biocor. CONCLUSIONS: The in vitro results indicate that it is possible to predict the net pressure gradient by Doppler in bileaflet and stented biologic valves. Our data indicate that important PR is also present in stented biologic valves.

Cross-sectional changes in an asymmetric tube with bearing on non-invasive pressure measurements.

This study aims at investigating the radial dimensional changes, as a result of an applied intraluminal pressure for an elastic tube with non-uniform cross section. The study is related to a method for non-invasive pressure measurement using the extracorporeal tube as part of the sensor. The intended application is for monitoring of pressure in the blood and dialysate tubes during haemodialysis. The intention is to find a tube cross-section geometry that results in an expansion of the tube so that it is suitable to use as a component in a pressure sensor. The tube should have high radial compliance and expand in a well-defined manner to be able to transfer the intraluminal pressure to a transducer element sensing the radial force. Radial expansion was studied experimentally for tubes with different cross-section geometries. For small tube expansions the resolution in the experimental measurements was not sufficient to study the radial expansion. In this case, numerical simulation was performed. We conclude that a tube with essentially elliptic outer surface and circular inner surface, with a relation of 1:2 between the size of the thin and thick wall, results in a radial expansion upon application of pressure indicating that this tube is suitable for use as part of a sensor.