For irreversible-electroporation (IRE)-based therapies, the underlying electric field distribution in the target tissue is influenced by the electroporation-induced conductivity changes and is important for predicting the treatment zone.
Objective: In this study, we characterized the liver tissue conductivity changes during high-frequency irreversible electroporation (H-FIRE) treatments of widths 5 and 10 μs and proposed a method for predicting the ablation zones.
Methods: To achieve this, we created a finite-element model of the tissue treated with H-FIRE and IRE pulses based on experiments conducted in an in-vivo rabbit liver study. We performed a parametric sweep on a Heaviside function that captured the tissue conductivity versus electric field behavior to yield a model current close to the experimental current during the first burst/pulse. A temperature module was added to account for the current increase in subsequent bursts/pulses. The evolution of the electric field at the end of the treatment was overlaid on the experimental ablation zones determined from hematoxylin and eosin staining to find the field thresholds of ablation.
Results: Dynamic conductivity curves that provided a statistically significant relation between the model and experimental results were determined for H-FIRE. In addition, the field thresholds of ablation were obtained for the tested H-FIRE parameters.
Conclusion: The proposed numerical model can simulate the electroporation process during H-FIRE.
Significance: The treatment planning method developed in this study can be translated to H-FIRE treatments of different widths and for different tissue types.