Purpose: Although theoretical modeling is widely used to study different aspects of radiofrequency ablation (RFA) its utility is directly related to its realism. An important factor in this realism is the use of mathematical functions to model the temperaturedependence of tissue thermal (k) and electrical () conductivities. Our aim was to review the piecewise mathematical functions most commonly used for modeling the temperaturedependence of k and  in RFA computational modeling. Materials and methods:We built a hepatic RFA theoretical model of a cooled electrode and compared lesion dimensions and impedance evolution with combinations of mathematical functions proposed in previous studies We employed the thermal damage contour D63 to compute the lesion dimension contour, which corresponds to Ω= 1, Ω being local thermal damage assessed by the Arrhenius damage model. Results:The results were very similar in all cases in terms of impedance evolution and lesion size after 6 minutes of ablation. Although the relative differences between cases in terms of time to first roll-off (abrupt increase in impedance) were as much as 12%, the maximum relative differences in terms of the short lesion (transverse) diameter were below 3.5%. Conclusions:The findings suggest that the different methods of modeling temperature dependence of k and  reported in the literature do not significantly affect the computed lesion diameter.3

For the same delivered energy, the new protocol significantly increases the magnitude of the applied electric field, which may be the reason why it is clinically more effective in achieving pain relief.

The objective of this study was to compare three different heat transfer models for radiofrequency ablation of in vivo liver tissue using a cooled electrode and three different voltage levels. The comparison was between the simplest but less realistic Pennes’ equation and two porous media-based models, i.e. the Local Thermal Non-Equilibrium (LTNE) equations and Local Thermal Equilibrium (LTE) equation, both modified to take into account two-phase water vaporization (tissue and blood). Different blood volume fractions in liver were considered and the blood velocity was modeled to simulate a vascular network. Governing equations with the appropriate boundary conditions were solved with Comsol Multiphysics finite-element code. The results in terms of coagulation transverse diameters and temperature distributions at the end of the application showed significant differences, especially between Pennes and the modified LTNE and LTE models. The new modified porous media-based models covered the ranges found in the few in vivo experimental studies in the literature and they were closer to the published results with similar in vivo protocol. The outcomes highlight the importance of considering the three models in the future in order to improve thermal ablation protocols and devices and adapt the model to different organs and patient profiles.

Berjano, E. (2012). Relationship between roll-off occurrence and spatial distribution of dehydrated tissue during RF ablation with cooled electrodes. International Journal of Hyperthermia. 28 (1) AbstractPurpose: To study the relationship between roll-off (sudden increase in impedance) and spatial distribution of dehydrated tissue during RF ablation using a cooled electrode (temperatures around 100ºC). Methods:We used a double approach: 1) theoretical modeling based on the Finite Element Method; and 2) 20 ablations using an experimental study on ex-vivo excised bovine liver in which we measured impedance progress and temperature at three points close to the electrode surface: 0.5 (T1), 1.5 (T2) and 2.5 (T3) mm from the tip. T2 was located exactly at the center of the 30 mm long electrode.Results: Temperatures at T1 and T3 quickly rose to 100ºC (at ≈ 20 and 40 s, respectively), while at the rise at T2 was somewhat slower, stabilized around 50 s and reached a maximum value of 99ºC at about 60 s. Impedance reached a minimum of 65  (plateau), began increasing at 50 s and continued rising throughout the procedure, reaching a value equal to the initial value at 70 s. Likewise, computed impedance dropped to ≈73  (plateau), began increasing at 50 s and reached an impedance value equal to the initial value at ≈78 s, which approximately coincided with the time when the entire zone surrounding the electrode was within the 100ºC isotherm. Conclusion:There is a close relationship between the moment at which roll-off occurs and the time when the entire electrode is completely encircled by the dehydrated tissue. The midelectrode zone is the last in which tissue dessication occurs.

Koruth, J.; Berjano, E. (2014). Radiofrequency cardiac ablation with catheters placed on opposing sides of the ventricular wall: Computer modelling comparing bipolar and unipolar modes. International Journal of Hyperthermia. 30 (6) Materials and methods:We built computational models to study the temperature distributions and lesion dimensions created by BM and UM on IVS and VFW during RFA. Two different UM types were considered: sequential (SEUM) and simultaneous (SIUM). The effect of ventricular wall thickness, catheter misalignment, epicardial fat, and presence of air in the epicardial space were also studied.Results: As regards IVS ablation, BM created transmural and symmetrical lesions for wall thicknesses up to 15 mm. SEUM and SIUM were not able to create transmural lesions with IVS thicknesses ≥ 12.5 and 15 mm, respectively. Lesions were asymmetrical only with SEUM. For VFW ablation, BM also created transmural lesions for wall thicknesses up to 15 mm. However, with SEUM and SIUM transmurality was obtained for VFW thicknesses ≤ 7.5 and 12.5 mm, respectively. With the three modes VFW lesions were always asymmetrical. In the scenario with air or a fat tissue layer on the epicardial side, only SIUM was capable of creating transmural lesions. Overall, BM was superior to UM in IVS and VFW ablation when the catheters were not aligned. Conclusions:Our findings suggest that BM is more effective than UM in achieving transmurality across both ventricular sites, except in the situation of the epicardial catheter tip surrounded by air or placed over a fat tissue layer.