Numerical simulation of radio frequency ablation with state dependent material parameters in three space dimensions

We present a model for the numerical simulation of radio frequency (RF) ablation of tumors with mono- or bipolar probes. This model includes the electrostatic equation and a variant of the well-known bio-heat transfer equation for the distribution of the electric potential and the induced heat. The equations are nonlinearly coupled by material parameters that change with temperature, dehydration and damage of the tissue. A fixed point iteration scheme for the nonlinear model and the spatial discretization with finite elements are presented. Moreover, we incorporate the effect of evaporation of water from the cells at high temperatures using a predictor-corrector like approach. The comparison of the approach to a real ablation concludes the paper.