As a part of ongoing efforts to develop computerized planning tools for cryosurgery, the current study focuses on developing an efficient numerical technique for bioheat transfer simulations. Our long-term goal is to develop a planning tool for cryosurgery that takes a 3D reconstruction of a target region, and suggests the best cryoprobe layout. Toward that goal, a planning algorithm, termed "force-field analogy," has been recently presented, based on a sequence of bioheat transfer simulations, which are by far the most computationally expensive part of the planning method. The objective in the current study is to develop a finite difference numerical scheme for bioheat transfer simulations, which reduces the overall run time of computerized planning, thereby making it clinically relevant. While the general concept of variable grid size and time intervals is not new, its application to the phase change problem of cryosurgery is the unique contribution of the current study.

Background: Cryosurgery is the destruction of undesired tissues by freezing, as in prostate cryosurgery, for example. Minimally-invasive cryosurgery is currently performed by means of an array of cryoprobes, each in the shape of a long hypodermic needle. The optimal arrangement of the cryoprobes, which is known to have a dramatic effect on the quality of the cryoprocedure, remains an art held by the cryosurgeon, based on the cryosurgeon's experience and "rules of thumb." An automated computerized technique for cryosurgery planning is the subject matter of the current report, in an effort to improve the quality of cryosurgery. Method of Approach:A two-phase optimization method is proposed for this purpose, based on two previous and independent developments by this research team. Phase I is based on a bubblepacking method, previously used as an efficient method for finite elements meshing. Phase II is based on a force-field analogy method, which has proven to be robust at the expense of a typically long runtime. Results:As a proof-of-concept, results are demonstrated on a 2D case of a prostate cross-section. The major contribution of this study is to affirm that in many instances cryosurgery planning can be performed without extremely expensive simulations of bioheat transfer, achieved in Phase I. Conclusions:This new method of planning has proven to reduce planning runtime from hours to minutes, making automated planning practical in a clinical time frame.

As part of an ongoing project to develop automated tools for cryosurgery planning, the current study focuses on the development of a 3D bubble packing method. A proof-of-concept for the new method is demonstrated on five prostate models, reconstructed from ultrasound images. The new method is a modification of an established method in 2D. Ellipsoidal bubbles are packed in the volume of the prostate in the current study; such bubbles can be viewed as a first-order approximation of a frozen region around a single cryoprobe. When all cryoprobes are inserted to the same depth, optimum planning was found to occur at about 60% of the length of the prostate (measured from its apex), which leads to cooling of approximately 75% of the prostate volume below a specific temperature threshold of −22°C. Bubble packing has the potential to dramatically reduce the run time for automated planning.

Due to time constraints in a clinical setup, bubble packing alone may be considered adequate for computerized planning. Furthermore, only bubble packing is demonstrated to be adequate for intra-operative planning.

The current study presents a computerized planning scheme for prostate cryosurgery using a variable insertion-depth strategy. This study is a part of an ongoing effort to develop computerized tools for cryosurgery. Based on typical clinical practices, previous automated planning schemes have required that all cryoprobes be aligned at a single insertion depth. The current study investigates the benefit of removing this constraint, in comparison with results based on uniform insertion-depth planning as well as the so-called "pullback procedure". Planning is based on the so-called "bubble-packing method", and its quality is evaluated with bioheat transfer simulations. This study is based on five 3D prostate models, reconstructed from ultrasound imaging, and cryoprobe active length in the range of 15 mm to 35 mm. The variable insertion depth technique is found to consistently provide superior results when compared to the other placement methods. Furthermore, it is shown that both the optimal active length and the optimal number of cryoprobes vary among prostate models, based on the size and shape of the target region. Due to its low computational cost, the new scheme can be used to determine the optimal cryoprobe layout for a given prostate model in real time.