Expanding the Bioheat Equation to Include Tissue Internal Water Evaporation During Heating

Numerical study of the induction of intratumoral apoptosis under microwave ablation by changing slot length of microwave coaxial antenna

Recent advances in technology have led to an increase in the detection of previously undetected deep-located tumor tissue. As a result, the medical field is using a variety of methods to treat deep-located tumors, and minimally invasive treatment techniques are being explored. In this study, therapeutic effect of microwave ablation (MWA) on tumor generated inside liver tissue was analyzed through numerical analysis. The distribution of electromagnetic fields in biological tissues emitted by microwave coaxial antenna (MCA) was calculated through the wave equation, and the thermal behavior of the tissue was analyzed through the Pennes bioheat equation. Among various treatment conditions constituting MWA, tumor radius and the slot length inside the MCA were changed, and the resulting treatment effect was quantitatively confirmed through three apoptotic variables. As a result, each tumor radius has optimal power condition for MWA, 2.6W, 2.4W, and 3.0W respectively. This study confirmed optimal therapeutic conditions for MWA. Three apoptotic variables were used to quantitatively identify apoptotic temperature maintenance inside tumor tissue and thermal damage to surrounding normal tissue. The findings of this study are expected to serve as a standard for treatment based on actual MWA treatment.

Non-Fourier Bioheat Transfer Analysis in Brain Tissue During Interstitial Laser Ablation: Analysis of Multiple Influential Factors

This work presents the dual-phase lag-based non-Fourier bioheat transfer model of brain tissue subjected to interstitial laser ablation. The finite element method has been utilized to predict the brain tissue's temperature distributions and ablation volumes. A sensitivity analysis has been conducted to quantify the effect of variations in the input laser power, treatment time, laser fiber diameter, laser wavelength, and non-Fourier phase lags. Notably, in this work, the temperature-dependent thermal properties of brain tissue have been considered. The developed model has been validated by comparing the temperature obtained from the numerical and ex vivo brain tissue during interstitial laser ablation. The ex vivo brain model has been further extended to in vivo settings by incorporating the blood perfusion effects. The results of the systematic analysis highlight the importance of considering temperature-dependent thermal properties of the brain tissue, non-Fourier behavior, and microvascular perfusion effects in the computational models for accurate predictions of the treatment outcomes during interstitial laser ablation, thereby minimizing the damage to surrounding healthy tissue. The developed model and parametric analysis reported in this study would assist in a more accurate and precise prediction of the temperature distribution, thus allowing to optimize the thermal dosage during laser therapy in the brain.

Proactive esophageal cooling during laser cardiac ablation: A computer modeling study

BACKGROUND AND OBJECTIVES

Laser ablation is increasingly used to treat atrial fibrillation (AF). However, atrioesophageal injury remains a potentially serious complication. While proactive esophageal cooling (PEC) reduces esophageal injury during radiofrequency ablation, the effects of PEC during laser ablation have not previously been determined. We aimed to evaluate the protective effects of PEC during laser ablation of AF by means of a theoretical study based on computer modeling.

METHODS

Three-dimensional mathematical models were built for 20 different cases including a fragment of atrial wall (myocardium), epicardial fat (adipose tissue), connective tissue, and esophageal wall. The esophagus was considered with and without PEC. Laser-tissue interaction was modeled using Beer-Lambert's law, Pennes' Bioheat equation was used to compute the resultant heating, and the Arrhenius equation was used to estimate the fraction of tissue damage (FOD), assuming a threshold of 63% to assess induced necrosis. We modeled laser irradiation power of 8.5 W over 20 s. Thermal simulations extended up to 250 s to account for thermal latency.

RESULTS

PEC significantly altered the temperature distribution around the cooling device, resulting in lower temperatures (around 22°C less in the esophagus and 9°C in the atrial wall) compared to the case without PEC. This thermal reduction translated into the absence of transmural lesions in the esophagus. The esophagus was thermally damaged only in the cases without PEC and with a distance equal to or shorter than 3.5 mm between the esophagus and endocardium (inner boundary of the atrial wall). Furthermore, PEC demonstrated minimal impact on the lesion created across the atrial wall, either in terms of maximum temperature or FOD.

CONCLUSIONS

PEC reduces the potential for esophageal injury without degrading the intended cardiac lesions for a variety of different tissue thicknesses. Thermal latency may influence lesion formation during laser ablation and may play a part in any collateral damage.

Optimal design of aperiodic tri-slot antennas for the conformal ablation of liver tumors using an experimentally validated MWA computer model

OBJECTIVE

This study aims to demonstrate that the conformal microwave ablation (MWA) of liver tumors could be attained by optimizing the structure of an aperiodic tri-slot coaxial antenna, its insertion depth, and input power.

METHODS

A computational MWA model with an aperiodic tri-slot coaxial antenna operating at the frequency of 2.45 GHz was built and validated by both an ex vivo and a pilot in vivo experiment with porcine healthy livers. The validated in vivo computational MWA model implemented with a liver tumor was then used as a testbed to investigate the conformal ablation of liver tumors. Five liver tumors in different sizes and shapes were investigated. A genetic algorithm optimization method (NSGA-II) was used to optimize the structure of antenna, insertion depth of antenna, and microwave antenna input power for the conformal ablation of liver tumors.

RESULTS

The validation results showed that a good agreement in both the spatiotemporal temperature distribution and ablation zone was found between the computer model and the ex vivo experiments at both 45 W, 5 min and 60 W, 3 min treatments and the in vivo experiment at 45 W, 5 min treatment. The optimized simulation results confirmed that five cases of liver tumors in different sizes and shapes can be conformally ablated by optimizing the aperiodic tri-slot coaxial antenna, antenna insertion depth, and microwave antenna input power.

CONCLUSION

This paper demonstrates that the aperiodic tri-slot coaxial antenna can be optimized with the insertion depth and input power for the conformal ablation of liver tumors, regardless the size and shape of liver tumors.

First validation of a model-based hepatic percutaneous microwave ablation planning on a clinical dataset

A model-based planning tool, integrated in an imaging system, is envisioned for CT-guided percutaneous microwave ablation. This study aims to evaluate the biophysical model performance, by comparing its prediction retrospectively with the actual ablation ground truth from a clinical dataset in liver. The biophysical model uses a simplified formulation of heat deposition on the applicator and a heat sink related to vasculature to solve the bioheat equation. A performance metric is defined to assess how the planned ablation overlaps the actual ground truth. Results demonstrate superiority of this model prediction compared to manufacturer tabulated data and a significant influence of the vasculature cooling effect. Nevertheless, vasculature shortage due to branches occlusion and applicator misalignment due to registration error between scans affects the thermal prediction. With a more accurate vasculature segmentation, occlusion risk can be estimated, whereas branches can be used as liver landmarks to improve the registration accuracy. Overall, this study emphasizes the benefit of a model-based thermal ablation solution in better planning the ablation procedures. Contrast and registration protocols must be adapted to facilitate its integration into the clinical workflow.

3D modeling of vector/edge finite element method for multi-ablation technique for large tumor-computational approach

Microwave ablation (MWA) is a cancer thermal ablation treatment that uses electromagnetic waves to generate heat within the tissue. The goal of this treatment is to eliminate tumor cells while leaving healthy cells unharmed. During MWA, excess heat generation can kill healthy cells. Hence, mathematical models and numerical techniques are required to analyze the heat distribution in the tissue before the treatment. The aim of this research is to explain the implementation of the 3D vector finite element method in a wave propagation model that simulates the specific absorption rate in the liver. The 3D Nedelec elements from H(curl; Ω) space are used to discretize the wave propagation model, and this implementation is helpful in solving many real-world problems that involve electromagnetic propagation with perfect conducting and absorbing boundary conditions. One of the difficulties in ablation treatment is creating a large ablation zone for a large tumor (diameter greater than 3 cm) in a short period of time with minimum damage to the surrounding tissue. This article addresses the aforementioned issue by introducing four antennas into the different places of the tumor sequentially and producing heat uniformly over the tumor. The results demonstrated that 95.5% of the tumor cells were killed with minimal damage to the healthy cells when the heating time was increased to 4 minutes at each position. Subsequently, we studied the temperature distribution and localised tissue contraction in the tissue using the three-dimensional bio-heat equation and temperature-time dependent model, respectively. The local tissue contraction is measured at arbitrary points in the domain and is more noticeable at temperatures higher than 102°C. The thermal damage in the liver during MWA treatment is investigated using the three-state cell death model. The system of partial differential equations is solved numerically due to the complex geometry of the domain, and the results are compared with experimental data to validate the models and parameters.

Optimization of laser dosimetry based on patient-specific anatomical models for the ablation of pancreatic ductal adenocarcinoma tumor

Laser-induced thermotherapy has shown promising potential for the treatment of unresectable primary pancreatic ductal adenocarcinoma tumors. Nevertheless, heterogeneous tumor environment and complex thermal interaction phenomena that are established under hyperthermic conditions can lead to under/over estimation of laser thermotherapy efficacy. Using numerical modeling, this paper presents an optimized laser setting for Nd:YAG laser delivered by a bare optical fiber (300 µm in diameter) at 1064 nm working in continuous mode within a power range of 2-10 W. For the thermal analysis, patient-specific 3D models were used, consisting of tumors in different portions of the pancreas. The optimized laser power and time for ablating the tumor completely and producing thermal toxic effects on the possible residual tumor cells beyond the tumor margins were found to be 5 W for 550 s, 7 W for 550 s, and 8 W for 550 s for the pancreatic tail, body, and head tumors, respectively. Based on the results, during the laser irradiation at the optimized doses, thermal injury was not evident either in the 15 mm lateral distances from the optical fiber or in the nearby healthy organs. The present computational-based predictions are also in line with the previous ex vivo and in vivo studies, hence, they can assist in the estimation of the therapeutic outcome of laser ablation for pancreatic neoplasms prior to clinical trials.

Theoretical Estimation of Tissue Thermal Response and Associated Thermal Damage During Gold Nanorod-enhanced Photothermal Therapy of Tumors

In the present work, we implemented a computational framework of in vivo gold nanorod (GNR)-enhanced photothermal therapy (PTT) for tumor treatment. The temperature-dependent thermophysical properties of biological tissue and the optical properties of both GNRs and the biological media were included. The latter were modulated during the treatment simulation to account for their variation, from the native to the coagulated state. The contribution of tissue injury-dependent blood perfusion was also considered. The developed model allowed for the estimation of temperature distribution during the photothermal procedure at different procedural settings and amounts of GNRs embedded in the tumor region (i.e., 12.5 μg, 25 μg, and 50 μg). Furthermore, the influence of GNRs on thermal injury, estimated with different damage models, was assessed. The inclusion of GNRs in the tumor entailed an increment of maximum tissue temperature, and faster heating kinetics, as witnessed by the lower time needed to reach complete thermal damage at the tumor center. The percentage of tumor thermal damage evaluated at the end of the simulated treatment was 48%, 69%, and 90%, for PTT in the presence of 12.5 μg, 25 μg, and 50 μg of GNRs, respectively.Clinical Relevance-This establishes that simulation-based tools, modeling the tissue properties variation during the photothermal treatment, can serve as promising preplanning platforms for nanoparticle-assisted light therapies.

Multiphysics Modeling of Low-Intensity Pulsed Ultrasound Induced Chemotherapeutic Drug Release from the Surface of Gold Nanoparticles

Currently, no numerical model for low-intensity pulsed ultrasound (LIPUS)-triggered anticancer drug release from gold nanoparticle (GNP) drug carriers exists in the literature. In this work, LIPUS-induced doxorubicin (DOX) release from GNPs was achieved in an ex vivo tissue model. Transmission electronic microscopy (TEM) imaging was performed before and after LIPUS exposure, and significant aggregation of the GNPs was observed upon DOX release. Subsequently, GNP surface potential was determined before and after LIPUS-induced DOX release, using a Zetasizer. A numerical model was then created to predict GNP aggregation, and the subsequent DOX release, via combining a thermal field simulation by solving the bioheat transfer equation (in COMSOL) and the Derjaguin, Landau, Verwey, and Overbeek (DLVO) total interaction potential (in MATLAB). The DLVO model was applied to the colloidal DOX-loaded GNPs by summing the attractive van der Waals and electrostatic repulsion interaction potentials for any given GNP pair. DLVO total interaction potential was found before and after LIPUS exposure, and an energy barrier for aggregation was determined. The DLVO interaction potential peak amplitude was found to drop from 1.36 k_BT to 0.24 k_BT after LIPUS exposure, translating to an 82.4% decrease in peak amplitude value. It was concluded that the interaction potential energy threshold for GNP aggregation (and, as a result, DOX release) was equal to 0.24 k_BT.

A vector finite element approach to temperature dependent parameters of microwave ablation for liver cancer

Microwave ablation (MWA) is a minimally invasive treatment for cancer that uses electromagnetic waves to kill the tumor cells without significantly damaging the surrounding healthy cells. A three-state cell death model calculates the thermal damage around the Hepatocellular carcinoma (HCC) tumor in the liver tissue. The temperature profile is simulated for a single-slot co-axial antenna with a 1 mm air slot located near the tip of the antenna to produce an adequate amount of heat. The aims of this study are (1) to use the vector/edge finite element method (VFEM) to simulate the electromagnetic wave propagation to obtain the specific absorption rate, which is an input for the bio-heat equation that predicts the heat distribution in the liver tissue during MWA treatment, and (2) to compare the computational costs of VFEM and the finite element method (FEM) when different types of input powers and dielectric properties are used in the wave propagation equation. This study claims that the accuracy level increases marginally with less computation cost while using VFEM for temperature-dependent wave propagation equation.

Transient changes during microwave ablation simulation : a comparative shape analysis

Microwave ablation therapy is a hyperthermic treatment for killing cancerous tumours whereby microwave energy is dispersed into a target tissue region. Modelling can provide a prediction for the outcome of ablation, this paper explores changes in size and shape of temperature and Specific absorption rate fields throughout the course of simulated treatment with different probe concepts. Here, an axisymmetric geometry of a probe embedded within a tissue material is created, solving coupled electromagnetic and bioheat equations using the finite element method, utilizing hp discretisation with the NGSolve library. Results show dynamic changes across all metrics, with different responses from different probe concepts. The sleeve probe yielded the most circular specific absorption rate pattern with circularity of 0.81 initially but suffered the largest reduction throughout ablation. Similarly, reflection coefficients differ drastically from their initial values, with the sleeve probe again experiencing the largest change, suggesting that it is the most sensitive the changes in the tissue dielectric properties in these select probe designs. These collective characteristic observations highlight the need to consider dielectric property changes and probe specific responses during the design cycle.

MWA Performed at 5.8 GHz through 'Side Firing' Approach: An Exploratory Study

Recent studies have shown that ablation techniques have the potential to eradicate adrenal adenomas while preserving the functionalities of the adrenal gland and the surrounding anatomical structures. This study explores a new microwave ablation (MWA) approach operating at 5.8 GHz and using anatomical and dielectric characteristics of the target tissue to create directional heating patterns. Numerical simulations are executed in planar and 3D adrenal models, considering two energy doses. The numerical study is refined accounting for the vaporization of the tissue water content. Ex vivo experimental evaluations on porcine adrenal models complete the study. The numerical and experimental results show that spherical ablation zones are able to cover the target for both energy doses considered. Nonetheless, most of the non-targeted tissues can be preserved from excessive heating when low energy level is used. Numerical models accounting for water vaporization are capable to foresee the experimental temperature values. This study shows that the proposed MWA directional approach operating at 5.8 GHz can be considered for creating effective and selective ablation zones.

An Analysis of Microwave Ablation Parameters for Treatment of Liver Tumors from the 3D-IRCADb-01 Database

Simulation techniques are powerful tools for determining the optimal conditions necessary for microwave ablation to be efficient and safe for treating liver tumors. Owing to the complexity and computational resource consumption, most of the existing numerical models are two-dimensional axisymmetric models that emulate actual three-dimensional cancers and the surrounding tissue, which is often far from reality. Different tumor shapes and sizes require different input powers and ablation times to ensure the preservation of healthy tissues that can be determined only by the full three-dimensional simulations. This study aimed to tailor microwave ablation therapeutic conditions for complete tumor ablation with an adequate safety margin, while avoiding injury to the surrounding healthy tissue. Three-dimensional simulations were performed for a multi-slot microwave antenna immersed in two tumors obtained from the 3D-IRCADb-01 liver tumors database. The temperature dependence of the dielectric and thermal properties of healthy and tumoral liver tissues, blood perfusion, and water content are crucial for calculating the correct ablation time and, thereby, the correct ablation process. The developed three-dimensional simulation model may help practitioners in planning patient-individual procedures by determining the optimal input power and duration of the ablation process for the actual shape of the tumor. With proper input power, necrotic tissue is placed mainly in the tumor, and only a small amount of surrounding tissue is damaged.

Modeling and ex vivo experimental validation of liver tissue carbonization with laser ablation

OBJECTIVE

The purpose of this study was to model the process of liver tissue carbonization with laser ablation (LA).

METHODS

A dynamic heat source model was proposed and combined with the light distribution model as well as bioheat transfer model to predict the development of tissue carbonization with laser ablation (LA) using an ex vivo porcine liver tissue model. An ex vivo laser ablation experiment with porcine liver tissues using a custom-made 1064 nm bare fiber was then used to verify the simulation results at 3, 5, and 7 W laser administrations for 5 min. The spatiotemporal temperature distribution was monitored by measuring the temperature changes at three points close the fiber during LA. Both the experiment and simulation of the temperature, tissue carbonization zone, and ablation zone were then compared.

RESULTS

Four stages were recognized in the development of liver tissue carbonization during LA. The growth of the carbonization zone along the fiber axial and radial directions were different in the four stages. The carbonization zone along the fiber axial direction (L2) grew in the four stages with a sharp increase in the initial period and a minor increase in Stage 4. However, the change in the carbonization zone along the fiber radial direction (D2) increased dramatically (Stage 1) to a long-time plateau (Stages 2 and 3) followed by a slow growth in Stage 4. An acceptable agreement between the computer simulation and ex vivo experiment in the temperature changes at the three points was found at all three testing laser administrations. A similar result was also obtained for the dimensions of coagulation zone and ablation zone between the computer simulation and ex vivo experiment (carbonization zone: 2.99± 0.10 vs. 2.78 mm², 67.39± 0.09 vs. 63.53 mm², and 90.53± 0.11 vs. 85.15 mm²; ablation zone: 68.95± 0.28 vs. 65.29 mm², 182.11± 0.24 vs. 213.81 mm², and 244.80± 0.06 vs. 251.79 mm² at 3, 5, and 7 W, respectively).

CONCLUSION

This study demonstrates that the proposed dynamic heat source model combined with the light distribution model as well as bioheat transfer model can predict the development of liver tissue carbonization with an acceptable accuracy. This study contributes to an improved understanding of the LA process in the treatment of liver tumors.

Resistance matching materials nanoarchitectonics for better performances in water evaporation-driven generators

The improvement of electricity production for water evaporation-driven generators (WEGs) remains a challenge. Herein, two types of WEGs were designed to study the resistance matching for improving the electricity production using the method of nanoarchitectonics. One type of reduced graphene oxide/carbon nanotube (RGO/CNT) WEG was constructed using RGO with adjustable resistances as working material and CNTs with fixed resistance as electrode material. The other type of graphene oxide (GO)/RGO WEG was constructed using RGO with adjustable resistance as electrode material and GO with fixed resistance was used as working material. The open circuit voltage of RGO/CNT increased from 15 to 56 mV and then decreased to 22 mV with increasing RGO resistance. The short circuit current of RGO/CNT also first increased and then decreased. The performance of GO/RGO was similar with that of RGO/CNT. Typically, the RGO/CNT and GO/RGO WEG showed the highest performance when the working material to electrode material resistance ratio was 2272 and 2365, respectively. It showed that the best resistance ratio of working material to electrode material was in the range of 2000-2500, which helped to improve about 2-5 times of electricity efficiency in the WEG. The present work provides a new direction for optimizing performance of WEGs.

On Efficacy of Microwave Ablation in the Thermal Treatment of an Early-Stage Hepatocellular Carcinoma

Microwave ablation at 2.45 GHz is gaining popularity as an alternative therapy to hepatic resection with a higher overall survival rate than external beam radiation therapy and proton beam therapy. It also offers better long-term recurrence-free overall survival when compared with radiofrequency ablation. To improve the design and optimization of microwave ablation procedures, numerical models can provide crucial information. A three-dimensional model of the antenna and targeted tissue without homogeneity assumptions are the most realistic representation of the physical problem. Due to complexity and computational resources consumption, most of the existing numerical studies are based on using two-dimensional axisymmetric models to emulate actual three-dimensional cancers and surrounding tissue, which is often far from reality. The main goal of this study is to develop a fully three-dimensional model of a multislot microwave antenna immersed into liver tissue affected by early-stage hepatocellular carcinoma. The geometry of the tumor is taken from the 3D-IRCADb-01 liver tumors database. Simulations were performed involving the temperature dependence of the blood perfusion, dielectric and thermal properties of both healthy and tumoral liver tissues. The water content changes during the ablation process are also included. The optimal values of the input power and the ablation time are determined to ensure complete treatment of the tumor with minimal damage to the healthy tissue. It was found that a multislot antenna is designed to create predictable, large, spherical zones of the ablation that are not influenced by varying tissue environments. The obtained results may be useful for determining optimal conditions necessary for microwave ablation to be as effective as possible for treating early-stage hepatocellular carcinoma, with minimized invasiveness and collateral damages.

Multiple cylindrical interstitial laser ablations (CILAs) of porcine pancreas in ex vivo and in vivo models

BACKGROUND

The therapeutic capacity of multiple cylindrical interstitial laser ablations (CILAs) of pancreatic tissue was evaluated with 1064 nm laser light in ex vivo and in vivo porcine pancreatic models.

METHODS

A diffusing applicator was sequentially employed to deliver 1064 nm laser light in a cylindrical distribution to ablate a large volume of pancreatic tissue. Ex vivo tissue was tested at various power levels (5, 7, and 10 W) under US imaging. An in vivo porcine model was used to evaluate the clinical feasibility of multiple CILAs on pancreatic tissue at 5 W via laparotomy (N = 3).

RESULTS

Multiple CILAs symmetrically ablated a range of ex vivo tissue volumes (2.4-6.0 cm³) at various power levels. Multiple CILAs warranted a therapeutic capacity of symmetrically ablating in vivo pancreatic tissue. Both ex vivo and in vivo pancreatic tissues after multiple CILAs at 5 W confirmed the absence of or minimal thermal injury to the peripheral tissue and carbonization.

CONCLUSIONS

The current findings suggest that the collective thermal effects from multiple CILAs can help widely ablate pancreatic tissue with minimal thermal injury. Further in vivo studies will investigate the safety of the proposed CILA treatment as well as acute/chronic responses of pancreatic tissue for clinical translations.

Short pulsed microwave ablation: computer modeling and ex vivo experiments

PURPOSE

To study the differences between continuous and short-pulse mode microwave ablation (MWA).

METHODS

We built a computational model for MWA including a 200 mm long and 14 G antenna from Amica-Gen and solved an electromagnetic-thermal coupled problem using COMSOL Multiphysics. We compared the coagulation zone (CZ) sizes created with pulsed and continuous modes under ex vivo and in vivo conditions. The model was used to compare long vs. short pulses, and 1000 W high-powered short pulses. Ex vivo experiments were conducted to validate the model.

RESULTS

The computational models predicted the axial diameter of the CZ with an error of 2-3% and overestimated the transverse diameter by 9-11%. For short pulses, the ex vivo computer modeling results showed a trend toward larger CZ when duty cycles decreases. In general, short pulsed mode yielded higher CZ diameters and volumes than continuous mode, but the differences were not significant (<5%), as in terms of CZ sphericity. The same trends were observed in the simulations mimicking in vivo conditions. Both CZ diameter and sphericity were similar with short and long pulses. Short 1000 W pulses produced smaller sphericity and similar CZ sizes under in vivo and ex vivo conditions.

CONCLUSIONS

The characteristics of the CZ created by continuous and pulsed MWA show no significant differences from ex vivo experiments and computer simulations. The proposed idea of enlarging coagulation zones and improving their sphericity in pulsed mode was not evident in this study.

Endoscopic ultrasound (EUS)-guided cylindrical interstitial laser ablation (CILA) on in vivo porcine pancreas

This study aims to demonstrate the feasibility of cylindrical interstitial laser ablation (CILA) in porcine pancreatic tissue to develop a EUS-guided PC ablation technique with enhanced safety. A diffusing applicator created a uniformly symmetrical laser ablation in pancreatic tissue. Ex vivo tests presented that both ablation thickness and volume increased linearly with the applied power (R² = 0.96 and 0.90, respectively) without carbonization and fiber degradation. The numerical simulations matched well with the experimental results in terms of temperature development and thermal damage (deviation of ≤ 15%). In vivo tests with EUS confirmed easy insertion and high durability of the diffusing applicator. EUS-guided CILA warranted a feasible therapeutic capacity of ablating in vivo pancreatic tissue. The proposed EUS-guided CILA can be a feasible therapeutic approach to treat PC with predictable thermal ablation and enhanced safety.

Thermomechanical Modeling of Laser Ablation Therapy of Tumors: Sensitivity Analysis and Optimization of Influential Variables

In cancer treatment, laser ablation is a promising technique used to induce localized thermal damage. Different variables influence the temperature distribution in the tissue and the resulting therapy efficacy; thus, the optimal therapy settings are required for obtaining the desired clinical outcome. In this work, thermomechanical modeling of contactless laser ablation was implemented to analyze the sensitivity of independent variables on the optimal treatment conditions. The Finite Element Method was utilized to solve the governing equations, i.e., the bioheat, mechanical deformation, and the Navier-Stokes equations. Validation of the model was evaluated by comparing experimental and simulated temperatures, which indicated high accuracy for estimating temperature. In particular, the results showed that the model is capable of estimating temperature with a good correlation factor (R=0.98) and low Mean Absolute Error (3.9 C). A sensitivity analysis based on laser irradiation time, power, beam distribution, and the blood vessel depth on temperature distribution and fraction of necrotic tissue was performed. Based on the most significant variables i.e., laser irradiation time and power, an optimization process was performed. This resulted into an indication of the optimal therapy settings for achieving maximum procedure efficiency i.e., the required fraction of necrotic tissue within the target volume, constituted by tumor and safety margins around it.

Measurement of Ex Vivo Liver, Brain and Pancreas Thermal Properties as Function of Temperature

The ability to predict heat transfer during hyperthermal and ablative techniques for cancer treatment relies on understanding the thermal properties of biological tissue. In this work, the thermal properties of ex vivo liver, pancreas and brain tissues are reported as a function of temperature. The thermal diffusivity, thermal conductivity and volumetric heat capacity of these tissues were measured in the temperature range from 22 to around 97 °C. Concerning the pancreas, a phase change occurred around 45 °C; therefore, its thermal properties were investigated only until this temperature. Results indicate that the thermal properties of the liver and brain have a non-linear relationship with temperature in the investigated range. In these tissues, the thermal properties were almost constant until 60 to 70 °C and then gradually changed until 92 °C. In particular, the thermal conductivity increased by 100% for the brain and 60% for the liver up to 92 °C, while thermal diffusivity increased by 90% and 40%, respectively. However, the heat capacity did not significantly change in this temperature range. The thermal conductivity and thermal diffusivity were dramatically increased from 92 to 97 °C, which seems to be due to water vaporization and state transition in the tissues. Moreover, the measurement uncertainty, determined at each temperature, increased after 92 °C. In the temperature range of 22 to 45 °C, the thermal properties of pancreatic tissue did not change significantly, in accordance with the results for the brain and liver. For the three tissues, the best fit curves are provided with regression analysis based on measured data to predict the tissue thermal behavior. These curves describe the temperature dependency of tissue thermal properties in a temperature range relevant for hyperthermia and ablation treatments and may help in constructing more accurate models of bioheat transfer for optimization and pre-planning of thermal procedures.

Laser-induced optothermal response of gold nanoparticles: From a physical viewpoint to cancer treatment application

Gold nanoparticles (GNPs)-based photothermal therapy (PTT) is a promising minimally invasive thermal therapy for the treatment of focal malignancies. Although GNPs-based PTT has been known for over two decades and GNPs possess unique properties as therapeutic agents, the delivery of a safe and effective therapy is still an open question. This review aims at providing relevant and recent information on the usage of GNPs in combination with the laser to treat cancers, pointing out the practical aspects that bear on the therapy outcome. Emphasis is given to the assessment of the GNPs' properties and the physical mechanisms underlying the laser-induced heat generation in GNPs-loaded tissues. The main techniques available for temperature measurement and the current theoretical simulation approaches predicting the therapeutic outcome are reviewed. Topical challenges in delivering safe thermal dosage are also presented with the aim to discuss the state-of-the-art and the future perspective in the field of GNPs-mediated PTT.

Role of saline concentration during saline-infused radiofrequency ablation: Observation of secondary Joule heating along the saline-tissue interface

Infusion of saline prior to radiofrequency ablation (RFA) is known to enlarge the thermal coagulation zone. The abundance of ions in saline elevate the electrical conductivity of the saline-saturated region. This promotes greater electric current flow inside the tissue, which increases the amount of RF energy deposition and subsequently enlarges the coagulation zone. In theory, infusion of higher concentration of saline should lead to larger coagulation zone due to the greater number of ions. Nevertheless, existing studies on the effects of concentration on saline-infused RFA have been conflicting, with the exact role of saline concentration yet to be fully elucidated. In this paper, computational models of saline-infused RFA were developed to investigate the role of saline concentration on the outcome of saline-infused RFA. The elevation in tissue electrical conductivity was modelled using the microscopic mixture model, while RFA was modelled using the coupled dual porosity-Joule heating model. Results obtained indicated that the presence of a concentration threshold to which no further elevation in tissue electrical conductivity and enlargement in thermal coagulation can occur. This threshold was determined to be at 15% NaCl. Analysis of the Joule heating distribution revealed the presence of a secondary Joule heating site located along the interface between wet and dry tissue. This secondary Joule heating was responsible for the enlargement in coagulation volume and its rapid growth phase during ablation.

A continuum thermomechanical model for the electrosurgery of soft hydrated tissues using a moving electrode

Electrosurgical radio-frequency heating of tissue is widely applied in minimally invasive surgical procedures to dissect tissue with simultaneous coagulation to obtain hemostasis. The tissue effect depends on the cumulative heating that occurs in the vicinity of the moving blade electrode. In this work, a continuum thermomechanical model based on mixture theory, which accounts for the multiphase nature of soft hydrated tissues and includes transport and evaporation losses, is used to capture the transient heating effect of a moving electrode. The model takes into account the dependence of electrical conductivity and the evaporation rate on the water content in the tissue, as it changes in response to heating. Temperature prediction is validated with mean experimental temperature measured during in situ experiments performed on porcine liver tissue at different power settings of the electrosurgical unit. The model is shown to closely capture the temperature variation in the tissue for three distinct scenarios; with no visible cutting or coagulation damage at a low 10 W power setting, with coagulation damage but no tissue cutting at an intermediate power setting of 25 W, and with both coagulation and tissue cutting at a higher power setting of 50 W. Furthermore, an Arrhenius model is shown to capture tissue damage observed in the experiments. Increase in applied power was found to correlate with tissue cutting and concentrated damage near the electrode, but had little effect on the observed coagulation damage width. The proposed model provides, for the first time, an accurate tool for predicting temperature rise and evolving damage resulting from a moving electrode in pure-cut electrosurgery.

Experimental assessment of microwave ablation computational modeling with MR thermometry

PURPOSE

Computational models are widely used during the design and characterization of microwave ablation (MWA) devices, and have been proposed for pretreatment planning. Our objective was to assess three-dimensional (3D) transient temperature and ablation profiles predicted by MWA computational models with temperature profiles measured experimentally using magnetic resonance (MR) thermometry in ex vivo bovine liver.

MATERIALS AND METHODS

We performed MWA in ex vivo tissue under MR guidance using a custom, 2.45 GHz water-cooled applicator. MR thermometry data were acquired for 2 min prior to heating, during 5-10 min microwave exposures, and for 3 min following heating. Fiber-optic temperature sensors were used to validate the accuracy of MR temperature measurements. A total of 13 ablation experiments were conducted using 30-50 W applied power at the applicator input. MWA computational models were implemented using the finite element method, and incorporated temperature-dependent changes in tissue physical properties. Model-predicted ablation zone extents were compared against MRI-derived Arrhenius thermal damage maps using the Dice similarity coefficient (DSC).

RESULTS

Prior to heating, the observed standard deviation of MR temperature data was in the range of 0.3-0.7°C. Mean absolute error between MR temperature measurements and fiber-optic temperature probes during heating was in the range of 0.5-2.8°C. The mean DSC between model-predicted ablation zones and MRI-derived Arrhenius thermal damage maps for 13 experimental set-ups was 0.95. When comparing simulated and experimentally (i.e. using MRI) measured temperatures, the mean absolute error (MAE %) relative to maximum temperature change was in the range 5%-8.5%.

CONCLUSION

We developed a system for characterizing 3D transient temperature and ablation profiles with MR thermometry during MWA in ex vivo liver tissue, and applied the system for experimental validation of MWA computational models.

Use of microwave ablation for thermal treatment of solid tumors with different shapes and sizes-A computational approach

Microwave Ablation (MWA) is one of the most recent developments in the field of thermal therapy. This approach is an effective method for thermal tumor ablation by increasing the temperature above the normal physiological threshold to kill cancer cells with minimum side effects to surrounding organs due to rapid heat dispersive tissues. In the present study, the effects of the shape and size of the tumor on MWA are investigated. To obtain the temperature gradient, coupled bio-heat and electromagnetic equations are solved using a three-dimensional finite element method (FEM). To extract cellular response at different temperatures and times, the three-state mathematical model was employed to achieve the ablation zone size. Results show that treatment of larger tumors is more difficult than that of smaller ones. By doubling the diameter of the tumor, the percentage of dead cancer cells is reduced by 20%. For a spherical tumor smaller than 2 cm, applying 50 W input power compared to 25 W has no significant effects on treatment efficiency and only increases the risk of damage to adjacent tissues. However, for tumors larger than 2 cm, it can increase the ablation zone up to 21%. In the spherical and oblate tumors, the mean percentage of dead cells at 6 GHz is nearly 30% higher than that at 2.45GHz, but for prolate tumors, treatment efficacy is reduced by 10% at a higher frequency. Moreover, the distance between two slots in the coaxial double slot antenna is modified based on the best treatment of prolate tumors. The findings of this study can be used to choose the optimum frequency and the best antenna design according to the shape and size of the tumor.

A numerical study on the no-touch bipolar radiofrequency ablation

The commonly used radiofrequency ablation (RFA) technique for treating liver cancer is in the monopolar mode. This requires the insertion of the RF electrode directly into the tumor tissue, which increases the risks of tumor track seeding (TTS). One way to overcome TTS is by employing the bipolar RFA, implemented in the no-touch mode. In the no-touch mode, two RF electrodes are inserted into the healthy tissue that surrounds the tumor. The distance between the electrodes and the tumor is defined as the no-touch gap. The ability of the no-touch bipolar RFA to overcome TTS has been demonstrated in laboratory studies; however, little is known about the thermo-physiological responses of the tissue during the ablation process of the no-touch procedure. This will be investigated numerically in the present study. A 3D model of the liver tissue is developed and the no-touch bipolar RFA implemented using a pair of RF electrodes is simulated using the finite element method. In particular, the effects of the no-touch gap on the treatment outcome of the RFA procedure are investigated. Results show that a larger no-touch gap may result incomplete tumor destruction due to the central region of the tumor not being directly affected by the Joule heating phenomenon that is more prominent around the electrodes. This suggests that an improperly selected no-touch gap may result in a reduced efficiency of the no-touch bipolar RFA.

Dynamic Impedance Monitoring for Large Diameter Vessel Sealing Using Bipolar Electrosurgery

Vessel sealing using bipolar electrosurgery is becoming a common practice in modern operating rooms. Despite all the advantages such as faster operation, less bleeding, and shorter postsurgery recovery time, side effects including sticking, charring, and rebleeding still occur, leading to increased surgery time and sometimes fatal complications. Tissue impedance during the electrosurgical process has been used to determine the electrical power of the process. However, little has been done to understand the dynamic tissue impedance and its effectiveness in monitoring the vessel sealing process. Moreover, the samples used in previous studies all had small diameters of 2–5 mm. In this study, an experimental setup was developed to perform vessel sealing tests using large-diameter blood vessel samples with mimicking blood flow. The tissue impedance during the heating process was obtained. Burst pressures after sealing were measured. A finite element simulation model was developed to understand the dynamic impedance behavior. It is seen that the tissue impedance increases rapidly in the beginning of the heating process and remains at a level that is several orders of magnitude higher than the initial value. This rapid impedance increase indicates protein denaturing, thus can be used to monitor the electrosurgical vessel sealing process. An impedance-based monitoring algorithm was developed, with which a burst pressure at least twice the normal human systolic blood pressure was achieved.

Thermal ablation of biological tissues in disease treatment: A review of computational models and future directions

Percutaneous thermal ablation has proven to be an effective modality for treating both benign and malignant tumours in various tissues. Among these modalities, radiofrequency ablation (RFA) is the most promising and widely adopted approach that has been extensively studied in the past decades. Microwave ablation (MWA) is a newly emerging modality that is gaining rapid momentum due to its capability of inducing rapid heating and attaining larger ablation volumes, and its lesser susceptibility to the heat sink effects as compared to RFA. Although the goal of both these therapies is to attain cell death in the target tissue by virtue of heating above 50°C, their underlying mechanism of action and principles greatly differs. Computational modelling is a powerful tool for studying the effect of electromagnetic interactions within the biological tissues and predicting the treatment outcomes during thermal ablative therapies. Such a priori estimation can assist the clinical practitioners during treatment planning with the goal of attaining successful tumour destruction and preservation of the surrounding healthy tissue and critical structures. This review provides current state-of-the-art developments and associated challenges in the computational modelling of thermal ablative techniques, viz., RFA and MWA, as well as touch upon several promising avenues in the modelling of laser ablation, nanoparticles assisted magnetic hyperthermia and non-invasive RFA. The application of RFA in pain relief has been extensively reviewed from modelling point of view. Additionally, future directions have also been provided to improve these models for their successful translation and integration into the hospital work flow.

Determination of Tissue Thermal Conductivity as a Function of Thermal Dose and Its Application in Finite Element Modeling of Electrosurgical Vessel Sealing

OBJECTIVE

Electrosurgical vessel sealing is a process commonly used to control bleeding during surgical procedures. Finite element (FE) modeling is often performed to obtain a better understanding of thermal spread during this process. The accuracy of the FE model depends on the implemented material properties. Thermal conductivity is one of the most important properties that affect temperature distribution. The goal of this study is to determine the tissue thermal conductivity as a function of thermal dose. Methods: We developed an iterative approach to correlating tissue thermal conductivity to more accurately calculated thermal dose, which cannot be experimentally measured. The resulting regression model was then implemented into an electrosurgical vessel sealing FE model to examine the accuracy of this FE model. Results: The results show that with the regression model, more reasonable temperature and thermal dose prediction can be achieved at the center of the sealed vessel tissue. The resulting electrical current and impedance from the FE model match with the experimental results. Conclusion: The developed approach can be used to determine the correlation between thermal dose and thermal conductivity. Describing the thermal conductivity as a function of thermal dose allows modeling of irreversible changes in tissue properties. Significance: By having a more accurate temperature estimation at the center of the sealed vessel, more insight is provided into how the tissue reacts during the vessel sealing process.

Numerical Analysis of Human Cancer Therapy Using Microwave Ablation

Microwave ablation is one type of hyperthermia treatment of cancer that involves heating tumor cells. This technique uses electromagnetic wave effects to kill cancer cells. A micro-coaxial antenna is introduced into the biological tissue. The radiation emitted by the antenna is absorbed by the tissue and leads to the heating of cancer cells. The diffuse increase in temperature should reach a certain value to achieve the treatment of cancer cells but it should be less than a certain other value to avoid damaging normal cells. This is why hyperthermia treatment should be carefully monitored. A numerical simulation is useful and may provide valuable information. The bio-heat equation and Maxwell’s equations are solved using the finite element method. Electro-thermal effects, temperature distribution profile, specific absorption rate (SAR), and fraction of necrotic tissue within cancer cells are analyzed. The results show that SAR and temperature distribution are strongly affected by input microwave power. High microwave power causes a high SAR value and raises the temperature above 50 °C, which may destroy healthy cells. It is revealed that with a power of 10 W, the tumor cells will be killed without damaging the surrounding tissue.

Waveform Dependent Electrosurgical Effects on Soft Hydrated Tissues

Electrosurgical procedures are ubiquitously used in surgery. The commonly used power modes, including the coagulation and blend modes, utilize non-sinusoidal or modulated current waveforms. For the same power setting, the coagulation, blend and pure cutting modes have different heating and thermal damage outcomes due to the frequency dependence of electrical conductivity of soft hydrated tissues. In this paper, we propose a multi-physics model of soft tissues to account for the effects of multi-frequency electrosurgical power modes within the framework of a continuum thermomechanical model based on mixture theory. Electrical and frequency spectrum results from different power modes at low and high power settings are presented. Model predictions are compared with in vivo electrosurgical heating experiments on porcine liver tissue. The accuracy of the model in predicting experimentally observed temperature profiles is found to be overall greater when frequency-dependence is included. An Arrhenius type model indicates that more tissue damage is correlated with larger duty cycles in multi-frequency modes.

Bio-Integrated Wearable Systems: A Comprehensive Review

Bio-integrated wearable systems can measure a broad range of biophysical, biochemical, and environmental signals to provide critical insights into overall health status and to quantify human performance. Recent advances in material science, chemical analysis techniques, device designs, and assembly methods form the foundations for a uniquely differentiated type of wearable technology, characterized by noninvasive, intimate integration with the soft, curved, time-dynamic surfaces of the body. This review summarizes the latest advances in this emerging field of "bio-integrated" technologies in a comprehensive manner that connects fundamental developments in chemistry, material science, and engineering with sensing technologies that have the potential for widespread deployment and societal benefit in human health care. An introduction to the chemistries and materials for the active components of these systems contextualizes essential design considerations for sensors and associated platforms that appear in following sections. The subsequent content highlights the most advanced biosensors, classified according to their ability to capture biophysical, biochemical, and environmental information. Additional sections feature schemes for electrically powering these sensors and strategies for achieving fully integrated, wireless systems. The review concludes with an overview of key remaining challenges and a summary of opportunities where advances in materials chemistry will be critically important for continued progress.

Modeling Heat Transfer in Tumors: A Review of Thermal Therapies

It is quite challenging to describe heat transfer phenomena in living systems because of the involved phenomena complexity. Indeed, thermal conduction and convection in tissues, blood perfusion, heat generation due to metabolism, complex vascular structure, changing of tissue properties depending on various conditions, are some of the features that make hard to obtain an accurate knowledge of heat transfer in living systems for all the clinical situations. This theme has a key role to predict accurately the temperature distribution in tissues, especially during biomedical applications, such as hyperthermia treatment of cancer, in which tumoral cells have to be destroyed and at the same time the surrounding healthy tissue has to be preserved. Moreover, the lack of experimentation in this field, due to ethical reasons, makes bioheat models even more significant. The first simple bioheat model was developed in 1948 by Pennes (J Appl Physiol 1:93-122, 1948) but it has some shortcomings that make the equation not so accurate. For this reason, over the years it has been modified and more complex models have been developed. The purpose of this review is to give a clear overview of how the bioheat models have been modified when applied in various hyperthermia treatments of cancer.

A continuum thermomechanical model of in vivo electrosurgical heating of hydrated soft biological tissues

Radio-frequency (RF) heating of soft biological tissues during electrosurgical procedures is a fast process that involves phase change through evaporation and transport of intra- and extra-cellular water, and where variations in physical properties with temperature and water content play significant role. Accurately predicting and capturing these effects would improve the modeling of temperature change in the tissue allowing the development of improved instrument design and better understanding of tissue damage and necrosis. Previous models based on the Pennes' bioheat model neglect both evaporation and transport or consider evaporation through numerical correlations, however, do not account for changes in physical properties due to mass transport or phase change, nor capture the pressure increase due to evaporation within the tissue. While a porous media approach can capture the effects of evaporation, transport, pressure and changes in physical properties, the model assumes free diffusion of liquid and gas without a careful examination of assumptions on transport parameters in intact tissue resulting in significant under prediction of temperature. These different approaches have therefore been associated with errors in temperature prediction exceeding 20% when compared to experiments due to inaccuracies in capturing the effects of evaporation losses and transport. Here, we present a model of RF heating of hydrated soft tissue based on mixture theory where the multiphase nature of tissue is captured within a continuum thermomechanics framework, simultaneously considering the transport, deformation and phase change losses due to evaporation that occur during electrosurgical heating. The model predictions are validated against data obtained for in vivo ablation of porcine liver tissue at various power settings of the electrosurgical unit. The model is able to match the mean experimental temperature data with sharp gradients in the vicinity of the electrode during rapid low and high power ablation procedures with errors less than 7.9%. Additionally, the model is able to capture fast vaporization losses and the corresponding increase in pressure due to vapor buildup which have a significant effect on temperature prediction beyond 100 °C.

Wireless, Battery-Free Epidermal Electronics for Continuous, Quantitative, Multimodal Thermal Characterization of Skin

Precise, quantitative measurements of the thermal properties of human skin can yield insights into thermoregulatory function, hydration, blood perfusion, wound healing, and other parameters of clinical interest. The need for wired power supply systems and data communication hardware limits, however, practical applicability of existing devices designed for measurements of this type. Here, a set of advanced materials, mechanics designs, integration schemes, and wireless circuits is reported as the basis for wireless, battery-free sensors that softly interface to the skin to enable precise measurements of its temperature and thermal transport properties. Calibration processes connect these parameters to the hydration state of the skin, the dynamics of near-surface flow through blood vessels and implanted catheters, and to recovery processes following trauma. Systematic engineering studies yield quantitative metrics in precision and reliability in real-world conditions. Evaluations on five human subjects demonstrate the capabilities in measurements of skin hydration and injury, including examples of continuous wear and monitoring over a period of 1 week, without disrupting natural daily activities.

The Influence of Dynamic Tissue Properties on HIFU Hyperthermia: A Numerical Simulation Study

Accurate temperature and thermal dose prediction are crucial to high-intensity focused ultrasound (HIFU) hyperthermia, which has been used successfully for the non-invasive treatment of solid tumors. For the conventional method of prediction, the tissue properties are usually set as constants. However, the temperature rise induced by HIFU irradiation in tissues will cause changes in the tissue properties that in turn affect the acoustic and temperature field. Herein, an acoustic–thermal coupling model is presented to predict the temperature and thermal damage zone in tissue in terms of the Westervelt equation and Pennes bioheat transfer equation, and the individual influence of each dynamic tissue property and the joint effect of all of the dynamic tissue properties are studied. The simulation results show that the dynamic acoustic absorption coefficient has the greatest influence on the temperature and thermal damage zone among all of the individual dynamic tissue properties. In addition, compared with the conventional method, the dynamic acoustic absorption coefficient leads to a higher focal temperature and a larger thermal damage zone; on the contrary, the dynamic blood perfusion leads to a lower focal temperature and a smaller thermal damage zone. Moreover, the conventional method underestimates the focal temperature and the thermal damage zone, compared with the simulation that was performed using all of the dynamic tissue properties. The results of this study will be helpful to guide the doctors to develop more accurate clinical protocols for HIFU treatment planning.

Lesion modeling, characterization, and visualization for image-guided cardiac ablation therapy monitoring

In spite of significant efforts to improve image-guided ablation therapy, a large number of patients undergoing ablation therapy to treat cardiac arrhythmic conditions require repeat procedures. The delivery of insufficient thermal dose is a significant contributor to incomplete tissue ablation, in turn leading to the arrhythmia recurrence. Ongoing research efforts aim to better characterize and visualize RF delivery to monitor the induced tissue damage during therapy. Here, we propose a method that entails modeling and visualization of the lesions in real-time. The described image-based ablation model relies on classical heat transfer principles to estimate tissue temperature in response to the ablation parameters, tissue properties, and duration. The ablation lesion quality, geometry, and overall progression are quantified on a voxel-by-voxel basis according to each voxel's cumulative temperature and time exposure. The model was evaluated both numerically under different parameter conditions, as well as experimentally, using ex vivo bovine tissue samples undergoing ex vivo clinically relevant ablation protocols. The studies demonstrated less than 5°C difference between the model-predicted and experimentally measured end-ablation temperatures. The model predicted lesion patterns were within 0.5 to 1 mm from the observed lesion patterns, suggesting sufficiently accurate modeling of the ablation lesions. Lastly, our proposed method enables therapy delivery feedback with no significant workflow latency. This study suggests that the proposed technique provides reasonably accurate and sufficiently fast visualizations of the delivered ablation lesions.

Technical Note: On Cardiac Ablation Lesion Visualization for Image-guided Therapy Monitoring

The delivery of insufficient thermal dose is a significant contributor to incomplete tissue ablation and leads to arrhythmia recurrence and a large number of patients requiring repeat procedures. In concert with ongoing research efforts aimed at better characterizing the RF energy delivery, here we propose a method that entails modeling and visualization of the lesions in real time. The described image-based ablation model relies on classical heat transfer principles to estimate tissue temperature in response to the ablation parameters, tissue properties, and duration. The ablation lesion quality, geometry, and overall progression is quantified on a voxel-by-voxel basis according to each voxel's cumulative temperature and time exposure. The model was evaluated both numerically under different parameter conditions, as well as experimentally, using ex vivo bovine tissue samples. This study suggests that the proposed technique provides reasonably accurate and sufficiently fast visualizations of the delivered ablation lesions.

Thermal mapping on male genital and skin tissues of laptop thermal sources and electromagnetic interaction

Since the development of communication devices and expansion of their applications, there have been concerns about their harmful health effects. The main aim of this study was to investigate laptop thermal effects caused by exposure to electromagnetic fields and thermal sources simultaneously; propose a nondestructive, replicable process that is less expensive than clinical measurements; and to study the effects of positioning any new device near the human body in steady state conditions to ensure safety by U.S. and European standard thresholds. A computer simulation was designed to obtain laptop heat flux from SolidWorks flow simulation. Increase in body temperature due to heat flux was calculated, and antenna radiation was calculated using Computer Simulation Technology (CST) Microwave Studio software. Steady state temperature and specific absorption rate (SAR) distribution in user's body, and heat flux beneath the laptop, were obtained from simulations. The laptop in its high performance mode caused 420 (W/m² ) peak two-dimensional heat flux beneath it. The cumulative effect of laptop in high performance mode and 1 W antenna radiation resulted in temperatures of 42.9, 38.1, and 37.2 °C in lap skin, scrotum, and testis, that is, 5.6, 2.1, and 1.4 °C increase in temperature, respectively. Also, 1 W antenna radiation caused 0.37 × 10^-3 and 0.13 × 10^-1 (W/kg) peak three-dimensional SAR at 2.4 and 5 GHz, respectively, which could be ignored in reference to standards and temperature rise due to laptop use. Bioelectromagnetics. 38:550-558, 2017. © 2017 Wiley Periodicals, Inc.

An In Vitro and Numerical Study of Moxibustion Therapy on Biological Tissue

OBJECTIVE

Moxibustion therapy achieves satisfactory therapeutic effects largely depending on the heat stimulation of burning moxa. Understanding the thermal characteristics of heating process is an effective way to reveal the underlying mechanisms of moxibustion therapy.

METHODS

This paper performs experimental study on temperature distributions of burning moxa sticks and fresh in vitro porcine abdominal tissue using an infrared camera and thermocouples. Meanwhile, a moxibustion model incorporating moxa stick burning model and tissue heat transfer model was established with consideration of radiation propagation and water evaporation.

RESULTS

The burning features of moxa sticks were acquired and the radiation energy generated by the burning moxa stick was absorbed and scattered in biological tissue, resulting in a large temperature gradient in the skin layer. And the water evaporation led to a mass loss and reduced skin surface temperature. The numerical model was verified by experimental results and the effects of moxibustion treatment distance and duration can be quantified based on model calculation.

CONCLUSION

The detailed heat transfer process of moxibustion was obtained experimentally and numerically. During moxibustion, the radiation attenuation and water evaporation have a significant influence on the energy transport in biological tissue which cannot be ignored. The treatment distance of 3 cm is the recommended value to achieve the treatment efficacy without thermal damage and pain.

SIGNIFICANCE

This research would reveal the underlying mechanisms of moxibustion therapy. Besides, the developed models are expected to establish a guideline for moxibustion clinical treatment.