Transport lattice models of heat transport in skin with spatially heterogeneous, temperature-dependent perfusion

BACKGROUND: Investigation of bioheat transfer problems requires the evaluation of temporal and spatial distributions of temperature. This class of problems has been traditionally addressed using the Pennes bioheat equation. Transport of heat by conduction, and by temperature-dependent, spatially heterogeneous blood perfusion is modeled here using a transport lattice approach. METHODS: We represent heat transport processes by using a lattice that represents the Pennes bioheat equation in perfused tissues, and diffusion in nonperfused regions. The three layer skin model has a nonperfused viable epidermis, and deeper regions of dermis and subcutaneous tissue with perfusion that is constant or temperature-dependent. Two cases are considered: (1) surface contact heating and (2) spatially distributed heating. The model is relevant to the prediction of the transient and steady state temperature rise for different methods of power deposition within the skin. Accumulated thermal damage is estimated by using an Arrhenius type rate equation at locations where viable tissue temperature exceeds 42°C. Prediction of spatial temperature distributions is also illustrated with a two-dimensional model of skin created from a histological image. RESULTS: The transport lattice approach was validated by comparison with an analytical solution for a slab with homogeneous thermal properties and spatially distributed uniform sink held at constant temperatures at the ends. For typical transcutaneous blood gas sensing conditions the estimated damage is small, even with prolonged skin contact to a 45°C surface. Spatial heterogeneity in skin thermal properties leads to a non-uniform temperature distribution during a 10 GHz electromagnetic field exposure. A realistic two-dimensional model of the skin shows that tissue heterogeneity does not lead to a significant local temperature increase when heated by a hot wire tip. CONCLUSIONS: The heat transport system model of the skin was solved by exploiting the mathematical analogy between local thermal models and local electrical (charge transport) models, thereby allowing robust, circuit simulation software to obtain solutions to Kirchhoff's laws for the system model. Transport lattices allow systematic introduction of realistic geometry and spatially heterogeneous heat transport mechanisms. Local representations for both simple, passive functions and more complex local models can be easily and intuitively included into the system model of a tissue

Image-guided focused ultrasound ablation of breast cancer: current status, challenges, and future directions

Image-guided focussed ultrasound (FUS) ablation is a non-invasive procedure that has been used for treatment of benign or malignant breast tumours. Image-guidance during ablation is achieved either by using real-time ultrasound (US) or magnetic resonance imaging (MRI). The past decade phase I studies have proven MRI-guided and US-guided FUS ablation of breast cancer to be technically feasible and safe. We provide an overview of studies assessing the efficacy of FUS for breast tumour ablation as measured by percentages of complete tumour necrosis. Successful ablation ranged from 20% to 100%, depending on FUS system type, imaging technique, ablation protocol, and patient selection. Specific issues related to FUS ablation of breast cancer, such as increased treatment time for larger tumours, size of ablation margins, methods used for margin assessment and residual tumour detection after FUS ablation, and impact of FUS ablation on sentinel node procedure are presented. Finally, potential future applications of FUS for breast cancer treatment such as FUS-induced anti-tumour immune response, FUS-mediated gene transfer, and enhanced drug delivery are discussed. Currently, breast-conserving surgery remains the gold standard for breast cancer treatment

PLANET: An Ellipse Fitting Approach for Simultaneous T1 and T2 Mapping Using Phase-Cycled Balanced Steady-State Free Precession

Purpose: To demonstrate the feasibility of a novel, ellipse fitting approach, named PLANET, for simultaneous estimation of relaxation times T1 and T2 from a single 3D phase-cycled balanced steady-state free precession (bSSFP) sequence. Methods: A method is presented in which the elliptical signal model is used to describe the phase-cycled bSSFP steadystate signal. The fitting of the model to the acquired data is reformulated into a linear convex problem, which is solved directly by a linear least squares method, specific to ellipses. Subsequently, the relaxation times T1 and T2, the banding free magnitude, and the off-resonance are calculated from the fitting results. Results: Maps of T1 and T2, as well as an off-resonance and a banding free magnitude can be simultaneously, quickly, and robustly estimated from a single 3D phase-cycled bSSFP sequence. The feasibility of the method was demonstrated in a phantom and in the brain of healthy volunteers on a clinical MR scanner. The results were in good agreement for the phantom, but a systematic underestimation of T1 was observed in the brain. Conclusion: The presented method allows for accurate mapping of relaxation times and off-resonance, and for the reconstruction of banding free magnitude images at realistic signalto- noise ratios. Magn Reson Med 000:000–000, 2017. VC 2017 The Authors Magnetic Resonance in Medicine published by Wiley Periodicals, Inc. on behalf of International Society for Magnetic Resonance in Medicine. This is an open access article under the terms of the Creative Commons Attribution-NonCommercial-NoDerivs License, which permits use and distribution in any medium, provided the original work is properly cited, the use is non-commercial and no modifications or adaptations are made

PLANET: An Ellipse Fitting Approach for Simultaneous T1 and T2 Mapping Using Phase-Cycled Balanced Steady-State Free Precession

Purpose: To demonstrate the feasibility of a novel, ellipse fitting approach, named PLANET, for simultaneous estimation of relaxation times T1 and T2 from a single 3D phase-cycled balanced steady-state free precession (bSSFP) sequence. Methods: A method is presented in which the elliptical signal model is used to describe the phase-cycled bSSFP steadystate signal. The fitting of the model to the acquired data is reformulated into a linear convex problem, which is solved directly by a linear least squares method, specific to ellipses. Subsequently, the relaxation times T1 and T2, the banding free magnitude, and the off-resonance are calculated from the fitting results. Results: Maps of T1 and T2, as well as an off-resonance and a banding free magnitude can be simultaneously, quickly, and robustly estimated from a single 3D phase-cycled bSSFP sequence. The feasibility of the method was demonstrated in a phantom and in the brain of healthy volunteers on a clinical MR scanner. The results were in good agreement for the phantom, but a systematic underestimation of T1 was observed in the brain. Conclusion: The presented method allows for accurate mapping of relaxation times and off-resonance, and for the reconstruction of banding free magnitude images at realistic signalto- noise ratios. Magn Reson Med 000:000–000, 2017. VC 2017 The Authors Magnetic Resonance in Medicine published by Wiley Periodicals, Inc. on behalf of International Society for Magnetic Resonance in Medicine. This is an open access article under the terms of the Creative Commons Attribution-NonCommercial-NoDerivs License, which permits use and distribution in any medium, provided the original work is properly cited, the use is non-commercial and no modifications or adaptations are made

MRI methods for the evaluation of high intensity focused ultrasound tumor treatment:current status and future needs

\u3cp\u3eThermal ablation with high intensity focused ultrasound (HIFU) is an emerging noninvasive technique for the treatment of solid tumors. HIFU treatment of malignant tumors requires accurate treatment planning, monitoring and evaluation, which can be facilitated by performing the procedure in an MR-guided HIFU system. The MR-based evaluation of HIFU treatment is most often restricted to contrast-enhanced T\u3csub\u3e1\u3c/sub\u3e-weighted imaging, while it has been shown that the non-perfused volume may not reflect the extent of nonviable tumor tissue after HIFU treatment. There are multiple studies in which more advanced MRI methods were assessed for their suitability for the evaluation of HIFU treatment. While several of these methods seem promising regarding their sensitivity to HIFU-induced tissue changes, there is still ample room for improvement of MRI protocols for HIFU treatment evaluation. In this review article, we describe the major acute and delayed effects of HIFU treatment. For each effect, the MRI methods that have been - or could be - used to detect the associated tissue changes are described. In addition, the potential value of multiparametric MRI for the evaluation of HIFU treatment is discussed. The review ends with a discussion on future directions for the MRI-based evaluation of HIFU treatment.\u3c/p\u3