- Research
- Open Access

# Computational study of the effects of arterial bifurcation on the temperature distribution during cryosurgery

- Yong-Chang Zheng†
^{1}, - Jun-Hong Wu†
^{2}, - Zhi-Zhu He
^{3}Email author and - Shao-Jiong Huang
^{3}Email author

**Received:**25 May 2017**Accepted:**10 January 2018**Published:**16 January 2018

## Abstract

### Background

Thermally significant blood flows into locally cooled diseased tissues and warm them during cryosurgery so that the iceball is often hard to cover the whole diseased volume. This paper is aimed at investigating the effects of large arterial bifurcation on the temperature distribution during cryosurgery through simulation method.

### Methods

A parametric geometry model is introduced to construct a close-to-real arterial bifurcation. The three-dimensional transient conjugate heat transfer between bifurcated artery and solid tissues with phase change during cryosurgery is performed by finite volume method.

### Results

The discussion was then made on the effects of the relative position between cryoprobe and artery bifurcation, the inlet velocity of root artery and the layout of multiple cryoprobes on the temperature distribution and iceball evolution. The results show that the thermal interaction between blood flow and iceball growth near bifurcation is considerable complex. The thermal effects of bifurcation could modulate the iceball morphology, severely weaken its freezing volume and prevent the blood vessel from being frozen.

### Conclusion

The present work is expected to be valuable in optimizing cryosurgery scheme of the situation that the bifurcated artery is embedded into the disease tissue.

## Keywords

- Bioheat transfer
- Cryosurgery
- Numerical simulation
- Arterial bifurcation
- Blood thermal effect

## Background

Cryosurgery has been widely demonstrated as an excellent therapeutic approach to destroy the diseased tissues (such as tumor) due to its minimal invasiveness [1–3]. In order to freeze and kill the target cells, largely decreasing the temperature of target tissues to below 0 °C is necessary. The temperature distribution in target tissue is the major parameter to evaluate the cryosurgical output. Accordingly, it is very important to accurately control the temperature distribution in order to enhance the destruction of diseased tissues and avoid the injury of healthy tissues.

Blood flow could remarkably affect the temperature distributions during freezing, especially with the presence of large blood vessels (larger than 0.5 mm in diameter) [4]. Thermally significant blood flows into locally cooled diseased tissues and warm them during cryosurgery so that the iceball is often hard to cover the whole diseased volume. Recently the study on thermal behavior of large blood vessels has attracted much attention in the cryosurgery area. Deng et al. [5] adopted infrared thermography system to investigate the thermal effects of large vessels during cryosurgery based on simulation and animal experiments. The results showed that the heating nature of the flowing blood in the large vessels could produce steep temperature gradients and inadequate cooling to the frozen tissues. Numerical simulation method was also used to study the thermal effects of the large blood vessel. In the early research [6], the blood velocity in the line-like vessel was simply considered as constant. The cylindrical cryoprobes and blood vessel were both approximated as cubes. Such simplification would induce numerical errors to the temperature distribution [7]. Latterly, a finite element method based on FEM commercial software was introduced to obtain a more accurate numerical solution of temperature field near blood vessel [8]. However, compared with the investigation for the thermal effects of large vessels on the hyperthermia ablation [9–14], study on the similar issues in cryosurgery is still rare.

Compared to single large line-like artery, the bifurcated artery (such as inlet artery in liver) has complicated structure so that the complex blood flow distribution would induce heterogeneous heat transfer surrounding artery bifurcation. In recent years, some investigations focus on the cooling effects of artery bifurcation on the hyperthermia ablation, such as microwave ablation [15, 16] and radio-frequency ablation [17]. However, few investigation contributes to the warm effects of artery bifurcation on the cryosurgery. The details of three-dimensional transient temperature distribution during cryosurgery surrounding artery bifurcation, which is more real and complex, are still unknown. The phase-change heat transfer combined with the convective mechanism of blood flow in arterial bifurcation during cryosurgery would be more complex than that during hyperthermia ablation. Tremendous contributions are needed to probe into such important issues, which are very useful for optimizing cryosurgery scheme of the situation that the bifurcated artery is embedded into the disease tissues. The aim of this paper is to disclose the detailed temperature characteristics of cryosurgery in the vicinity of an arterial bifurcation.

In present study, a parametric geometry model [18] is introduced to construct a close-to-real bifurcated artery. The steady state blood flow is considered here. The three-dimensional transient conjugate heat transfer between bifurcated artery and solid tissues during cryosurgery is performed by finite volume method. Then the relative position between cryoprobe and artery bifurcation, the changes of root vessel inlet velocity and the layout of multiple cryoprobes would be considered to investigate the detailed temperature distribution between bifurcation and cryoprobe. The mechanism of heat transfer between artery bifurcation and solid tissues, the iceball edge evolution would be then revealed to evaluate the thermal effects of large artery bifurcation on the tissue temperature distribution and iceball growth during cryosurgery.

## Geometric and mathematical model

### Geometric model

_{1}= 20 mm is the distance from root vessel inlet to bifurcation, L

_{2}= 80 mm is the height of daughter vessel along vertical direction, D

_{0}= 10 mm is the diameter of root vessel [17], \({\text{D}}_{1} = {\text{D}}_{2} = \sqrt[3]{{{1 \mathord{\left/ {\vphantom {1 2}} \right. \kern-0pt} 2}}}{\text{D}}_{0}\) 7.94 mm is the diameter of daughter vessels, Ф = 60° is the angle of bifurcation. The arterial bifurcation is embedded in cylindrical solid tissues with diameter D

_{ t }= 120 mm and length L

_{ t }= L

_{1}+ L

_{2}= 100 mm (see Fig. 1b). Single cylindrical cryoprobe with diameter D

_{p}= 4 mm inserts the solid tissue near the bifurcation site along horizontal direction. The distance from cryoprobe center to bifurcation point denotes as L

_{d}(see Fig. 1a). The probe shape is illustrated in Fig. 1c. The cryoprobe has two parts, adiabatic shaft and active tip with length L

_{p}= 22 mm, which could lead to an extremely low temperature (such as − 196 °C of liquid nitrogen).

### Governing equation

*T*>

*T*

_{ u }), transition field (

*T*

_{ l }≤

*T*≤

*T*

_{ u }

*)*and frozen area (

*T*<

*T*

_{ l }

*)*according to temperature distribution, where

*T*

_{ l }and

*T*

_{ u }denote respectively the lower and upper phase transition temperature of solid tissues. Then the Pennes bioheat transfer equation for whole solid tissues based on the effective heat capacity method [6] could be written as

*T*

_{ t }is the temperature of solid tissues and

*T*

_{ cb }denotes the temperature of blood perfusion from capillary vessel.

*C*

_{ b }is the heat capacity of blood. In addition, the detailed expression of \(\hat{C}\) the effective tissue heat capacity, \(\hat{\kappa }\) the effective tissue thermal conductivity, \(\hat{\omega }_{cb}\) the effective blood perfusion and \(\hat{Q}_{m}\) the effective tissue metabolic heat generation are determined by solid tissue state, which is given by

The thermo-physical parameters involved in the above equation could be found in [6]:*C*_{
t
} = *C*_{
b
} = 3.6 MJ/m^{3}K, *C*_{
f
} = 1.8 MJ/m^{3}K, *ω*_{
cb
} = 5 × 10^{−4}/s, *Q*_{
m
} = 420 J/m^{3}, *κ*_{
f
} = 2 W/mK, *κ*_{
t
} = *κ*_{
b
} = 0.5 W/mK, *L*_{
f
} = 250 MJ/m^{3}, T_{cb} = 37 °C, *T*_{
u
} = − 1 °C and *T*_{
l
} = − 8 °C. Blood viscosity is *μ* = 2.5 × 10^{−3} Ns/m^{2}.

*T*

_{ b }denotes the temperature in bifurcated artery. The heat transfer in solid tissues described by Eq. (1) and bifurcated artery presented in Eq. (5) could be combined by the interface between them, which is a typical conjugate heat transfer problem.

### Boundary condition

For solid tissues boundary, the thermal condition is considered as adiabatic wall boundary. The inlet velocity of the root artery adopts parabolic velocity profile *u*(*r*) = 2 V[1 − (2*r*/*D*_{0})^{2}], where V is average velocity and *r* the radial position. The inlet temperature of the root artery assumes as constant *T*_{0} = 37 °C. The outlet boundary of the daughter bifurcation artery adopts the pressure boundary, where the reference pressure is set as zero. The interface between solid tissues and cryoprobes tip is set as *T*_{
p
} = − 196 °C.

### Finite volume analysis tool

Firstly, the geometrical model of bifurcated artery constructed from software Solidworks imports into mesh generation software Gambit. Mesh size in artery domains is 1 mm so that about 10,000 tetrahedral elements and 20,000 nodes are obtained. The nonuniform mesh is used to map the solid tissues. The mesh size 0.6 mm is used for the field close to active cryoprobe, where the large temperature gradient happens during cryosurgery. The mesh size 2 mm is used for the field far away from the active cryoprobe. Thus there are about 1,500,000 tetrahedral elements and 270,000 nodes for the whole solid tissues. After creating the mesh, Gambit generates an input file to the finite volume package Fluent 6.3, which has been extensively utilized to address a variety of practical engineering problems nowadays. The blood flow and temperature analysis are solved according to SIMPLE and second order up wind algorithms.

## Results and discussion

### Steady blood velocity distribution

### Thermal interaction between iceball evolution and arterial bifurcation

_{d}= 20 mm between cryoprobe and bifurcation point is considered here. Figure 3 shows the temperature distribution in the plane x = 0, y = 0 and z = − 30 mm of bifurcated artery (Fig. 3a) and the iceball pattern view from x axis direction and temperature contours in the plane x = 0 (Fig. 3b) at freezing time t = 20 min. The iceball boundary in present study is defined as iso-surface of T = 0 °C. In order to focus on the temperature variation in the whole bifurcated artery, the contour interval is here set as (32.5, 37.5 °C) (Fig. 3a). In fact, the lowest temperature for artery domains, which happens on the inner wall of branching artery closest to cryoprobe, could reach 5 °C. This temperature could not decrease with freezing time evolution. It indicates that the warm blood flow prevents the artery from being frozen. However, one should note that the blood in the artery may be frozen for a shorter distance L

_{d}or a smaller inlet velocity of the root artery. The temperature in the vicinity of bifurcation point keeps at a relatively high level about 35 °C. From Fig. 3a, one can see that the bifurcated artery domains close to cryoprobe has a lower temperature distribution so that the following catchment area would be cooling. From Fig. 3b, we can see that the iceball stops evolution when encountering the arterial bifurcation. The large temperature gradient between iceball and artery wall is helpful to convective heat transfer which would enhance the thermal effects of blood flow. Different from line-like artery, the arterial bifurcation would lead to a more irregular iceball, which may make the target tissue frozen incompletely.

_{p}= 4 mm) for its different positions, L

_{d}= 20 mm (Fig. 4a), L

_{d}= 30 mm (Fig. 4b) and Ld = 40 mm (Fig. 4c) at freezing time t = 20 min. It is easy to find from Fig. 4 that the intensity of thermal interaction between iceball and artery bifurcation becomes weaker with increase of the relative distance between cryoprobe and bifurcation. The shorter distance between cryoprobe and artery would advance their thermal interaction during freezing. The iceball formation induced by cryoprobe has a limited volume. Thus the thermal effects from warm blood flow on the iceball would disappear when their distance is beyond available range. Clinical practice for cryosurgery has demonstrated that the most serious cryoinjury is achieved when the target tissue undergo below so-called lethal temperature, which is considered currently as in the range of − 50 to − 40 °C [3]. Table 1 presents the volume of iceball and lethal domains (where the temperature is below − 40 °C) for different cryoprobe positions. Both the volume of iceball and lethal domain increase when the cryoprobe is far away from arterial bifurcation.

The iceball volume (T = 0 °C) and the lethal volume (T = − 40 °C) for different cryoprobe positions at freezing time t = 20 min

Cryoprobe position (mm) | Ld = 20 | Ld = 30 | Ld = 40 |
---|---|---|---|

Iceball volume (cm | 39.13 | 41.90 | 51.90 |

Lethal volume (cm | 6.81 | 8.33 | 9.48 |

_{p}= 4 mm) positions L

_{d}= 20 mm, L

_{d}= 30 mm and L

_{d}= 40 mm. The temperature remains at a high level for L

_{d}= 20 mm due to strong thermal effects of warm blood flow. At the beginning of the freezing process, the temperature of the recorded point has the consistent cooling rate for each position. Such consistence would keep longer time with increase of L

_{d}. The reason is that at freezing initial stage a large distance between the iceball and artery leads to a weak thermal interaction. Thus the temperature of the recorded point close to cryoprobe is mainly determined by iceball at initial freezing stage. Figure 6 shows the total heat flux evolution of bifurcated artery surface for different cryoprobe (D

_{p}= 4 mm) positions. With increase of Ld, the total heat flux of bifurcated artery surface significantly decrease. This is mainly because the longer distance weaken the thermal effects of the blood flow.

The iceball volume (T = 0 °C) and the lethal volume (T = − 40 °C) for different inlet velocity with the same cryoprobe position L_{d} = 30 mm and freezing time t = 20 min

Inlet velocity (m/s) | V = 0.01 | V = 0.05 | V = 0.10 | V = 0.20 | V = 0.30 |
---|---|---|---|---|---|

Iceball volume (cm | 44.32 | 42.59 | 42.00 | 41.90 | 41.83 |

Lethal volume (cm | 8.56 | 8.41 | 8.35 | 8.33 | 8.32 |

The treatment time of cryosurgery when the minimum temperature of the bifurcated artery surface approaches to the freezing temperature (T = 0 °C) according to different cryoprobe positions and inlet velocities of artery root

Ld = 10 (mm) | Ld = 15 (mm) | Ld = 20 (mm) | Ld = 30 (mm) | Ld = 40 (mm) | |
---|---|---|---|---|---|

V = 0.01 (m/s) | 9.93 s | 40.50 s | 109.15 s | 731.50 s | ∞ |

V = 0.05 (m/s) | 11.10 s | 51.80 s | 296.40 s | ∞ | ∞ |

V = 0.10 (m/s) | 11.40 s | 58.05 s | ∞ | ∞ | ∞ |

V = 0.20 (m/s) | 11.50 s | 61.90 s | ∞ | ∞ | ∞ |

V = 0.30 (m/s) | 11.60 s | 65.30 s | ∞ | ∞ | ∞ |

_{p}= 2 mm and active length L

_{p}= 20 mm are distributed uniformly on the circumference with radius 20 mm. Its center locates on the center-line of cylinder with L

_{d}= 30 mm. Figure 9 shows that the iceball view from x direction and y direction for three cryoprobes at different freezing time. Three single iceballs are formed at initial cooling stage and combined fast to a large ingle iceball under a high-strength freezing. It can be seen from Fig. 9 that the three iceball induced by different cryoprobes have different growth rates due to varied cryoprobe positions and the inhomogeneous convective heat transfer between solid tissue and arterial bifurcation.

## Conclusions

In summary, the present paper has adopted three dimensional numerical simulation method to investigate the thermal effects of arterial bifurcation on temperature responses during cryosurgery based on single and multiple cryoprobe system. We have investigated in detail that the blood velocity distribution in arterial bifurcation and its effects on the iceball growth. The results indicate that complex blood velocity distribution could induce the inhomogeneous convective heat transfer between solid tissue and arterial bifurcation. Thus the iceball near arterial bifurcation presents strong irregular geometry. The blood flow of bifurcated artery has significant heating effects on the target freezing domains. It is also noteworthy that the artery wall is easily suffering from cold injury, which should be paid by special attention. In order to protect the artery wall, nanoparticle and external fields [19, 20] could be applied to enhance heat transfer near artery wall with low temperature.

## Notes

## Declarations

### Authors’ contributions

YCZ and JHW performed the computations, analyzed the data, and wrote the draft paper. SJH analyzed the data and wrote the draft paper. ZZH conceived and designed this study. All authors read and approved the final manuscript.

### Acknowledgements

This work is supported by the National Natural Science Foundation of China under Grant Nos. 51476181 and 61675236.

### Competing interests

The authors declare that they have no competing interests.

### Consent for publication

Not applicable.

### Ethics approval and consent to participate

The authors declare that this study does not involve human subjects, human material and human data.

### Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

**Open Access**This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

## Authors’ Affiliations

## References

- Rubinsky B. Cryosurgery. Annu Rev Biomed Eng. 2000;2:157–87.View ArticleGoogle Scholar
- Gage AA, Baust JG. Cryosurgery for tumors—a clinical overview. Technol Cancer Res Treat. 2004;3:187–99.View ArticleGoogle Scholar
- Gage AA, Baust J. Mechanisms of tissue injury in cryosurgery. Cryobiology. 1998;37:171–86.View ArticleGoogle Scholar
- Kotte A, Leeuwen GV, Bree JD, et al. A description of discrete vessel segments in thermal modeling of tissues. Phys Med Biol. 1996;41:865–84.View ArticleGoogle Scholar
- Deng ZS, Liu J, Wang HW. Disclosure of the significant thermal effects of large blood vessels during cryosurgery through infrared temperature mapping. Int J Therm Sci. 2008;47:530–45.View ArticleGoogle Scholar
- Deng ZS, Liu J. Numerical study of the effects of large blood vessels on three-dimensional tissue temperature profiles during cryosurgery. Numer Heat Transf Part A Appl. 2006;49:47–67.View ArticleGoogle Scholar
- He ZZ, Liu J. The effects of blood flow on the iceball evolution during a multiple probes cryosurgery. In: ASME 2011 International Mechanical Engineering Congress & Exposition. 2011.Google Scholar
- Kim C, O’Rourke AP, Will JA, et al. Finite-element analysis of hepatic cryoablation around a large blood vessel. IEEE Trans Biomed Eng. 2008;55:2087–93.View ArticleGoogle Scholar
- Horng TL, Lin WL, Liauh CT, et al. Effects of pulsatile blood flow in large vessels on thermal dose distribution during thermal therapy. Med Phys. 2007;34:1312–20.View ArticleGoogle Scholar
- Kolios MC, Sherar MD, Hunt JW. Large blood-vessel cooling in heated tissues—a numerical study. Phys Med Biol. 1995;40:477–94.View ArticleGoogle Scholar
- Payne S, Flanagan R, Pollari M, et al. Image-based multi-scale modelling and validation of radio-frequency ablation in liver tumours. Phil Trans R Soc A Math Phys Eng Sci. 2011;369:4233–54.MathSciNetView ArticleMATHGoogle Scholar
- Huang HW, Shih TC, Liauh CT, et al. Computer simulation of 3-d temperature and power distributions in tissue with a countercurrent blood vessels network during hyperthermia. J Med Biol Eng. 2009;29:252–8.Google Scholar
- Lagendijk JJW. Hyperthermia treatment planning. Phys Med Biol. 2000;45:61–76.View ArticleGoogle Scholar
- Chen ZP, Roemer RB. The effects of large blood-vessels on temperature distributions during simulated hyperthermia. J Biomech Eng Trans ASME. 1992;114:473–81.View ArticleGoogle Scholar
- Liu YJ, Qiao AK, Nan Q, et al. Thermal characteristics of microwave ablation in the vicinity of an arterial bifurcation. Int J Hyperth. 2006;22:491–506.View ArticleGoogle Scholar
- Lu YL, Nan Q, Du JJ, et al. Experimental study on thermal field in the vicinity of arterial bifurcation in microwave ablation therapy. Int J Hyperth. 2010;26:316–26.View ArticleGoogle Scholar
- Hariharan P, Chang I, Myers MR, et al. Radio-frequency ablation in a realistic reconstructed hepatic tissue. J Biomech Eng Trans ASME. 2007;129:354–64.Google Scholar
- Zakaria H, Robertson AM, Kerber CW. A parametric model for studies of flow in arterial bifurcations. Ann Biomed Eng. 2008;36:1515–30.View ArticleGoogle Scholar
- Bhatti MM, Zeeshan A, Ellahi R. Endoscope analysis on peristaltic blood flow of Sisko fluid with titanium magneto-nanoparticles. Comput Biol Med. 2016;78:29–41.View ArticleGoogle Scholar
- Abbas MA, Bai YQ, Rashidi MM, Bhatti MM. Application of drug delivery in Magnetohydrodynamics peristaltic blood flow of nanofluid in a non-uniform channel. J Mech Med Biol. 2016;16:1650052–66.View ArticleGoogle Scholar