Computational algorithm of high-intensity focused ultrasound beam of cancer tissue model for hyperthermia therapy


  • Kitsakorn Locharoenrat Dept. of Physics, Biomedical Physics Research Unit, Faculty of Science, King Mongkut’s Institute of Technology Ladkrabang, Bangkok 10520, Thailand
  • Kittiphot Songkaitiwong Dept. of Physics, Biomedical Physics Research Unit, Faculty of Science, King Mongkut’s Institute of Technology Ladkrabang, Bangkok 10520, Thailand


Pennes bioheat, heat transfer, hyperthermia, tumor, ultrasound


Ultrasound source mask is applied to generate high intensity focused ultrasound radiation force in 2D model of human breast carcinoma for hyperthermia therapy. Computational operation is based on Pennes bioheat’s concept. This method provides precise heat transfer values relied on thermal conduction in soft tissue and thermal convection in blood vessel domain. Ultrasound beam at 1 MHz is firstly lateral focused in tumor size of 15.0 mm x 28.5 mm at different focal depths without elevational focalization. Length of each focusing point of ultrasound beam is 6.8 – 45.4 mm in vertical axis, whereas full width at half maximum is 1.1 – 2.5 mm in horizontal axis. Simulated results have been shown that a discrepancy of ultrasound pressure around focus area rises with focal depth. Like pressure profile, when focal depth closes to ultrasound source, thermal homogeneity around focus area is attained, whereas thermal uniformity around focus area becomes worse with increasing focal depth. Using data visualization arrangement, acquired temperature profile corresponding to obtained pressure profile is converted to attain image of 2D model of human breast carcinoma. It is seen that ablated tumors are achieved, whilst normal tissues surrounding them are safe.


Arshad S., Baleanu, D., Haund, F. & Tang, Y. (2016) Dynamical analysis of fractional order model of immunogenic tumors. Advances in Mechanical Engineering 8(7): Article ID 1687814016656704.

Arshad, S., Baleanu, D., Defterli, O. & Shumaila, A. (2019) Numerical framework for the approximate solution of fractional tumor-obesity model. International Journal of Modeling Simulation Science and Computing; 10(1): Article ID 1941008.

Becker, S. & Kuznetsov, A. (2014) Heat transfer and fluid flow in biological processes. Academic Press, New York. Pp.55.

Chapra, S.C. & Canale, R.P. (2015) Numerical methods for engineers. Mc Graw Hill, New York. Pp.11.

Gomez-Aguilar, J.F., Lopez-Lopez, M.G. & Alvarado-Martinez, V.M. (2017) Chaos in a cancer model via fractional derivatives with exponential decay and Mittag-Leffler Law. Entropy, 19(12): Article ID 681.

Gonzalez, F.J. (2007) Thermal simulation of breast tumor. Revistra Maxixana De Fesica 53(4):323-326.

Grisey, A., Yon, S., Letort, V. & Lafitte, P. (2016) Simulation of high-intensity focused ultrasound lesions in presence of boiling. Journal of Therapeutic Ultrasound, 4(11).

Guzman-Cabrera, R., Guzman-Sepulveda, J.R., Parada, A.G. & Garcia, J.R. (2016) Digital processing of thermographic images for medical applications. Revista De Chimie, 67(1): 53-56.

IT’IS Foundation (2019) Database at a Glance, (accessed 01 Jan 2019).

Li, C., Littrup, P. & Huang, L. (2009) In vivo breast sound-speed imaging with ultrasound tomography. Ultrasound in Medical Biology, 35(10):1615-1628.

Newman, L.A. & Bensenhaver, J.M. (2015) Ductal carcinoma in-situ and microinvasive/borderline breast cancer. Springer, Heidelberg. Pp.22.

Treeby, B., Cox, B. & Jaros, J. (2019) A Matlab toolbox for the time-domain

simulation of acoustic wave fields, (accessed 01 Jan 2019).

Treeby, B.E., Jaros, J., Rendell, A.P. & Cox, B.T. (2012) Modeling nonlinear ultrasound propagation in heterogeneous media with power law absorption using a k-space pseudospectral method. Journal of Acoustic Society of America, 131(6):4324-4336.

Treeby, B.E., Tumen, M. & Cox, B.T. (2011) Medical image computing and computer-assisted intervention. Springer, Heidelberg. Pp.67.

Veronesi, U., Goldhirsch, A., Veronesi, P., Gentilini, O.D. & Leonardi, M.C. (2017) Breast cancer: Innovations in research and management. Springer, Heidelberg. Pp.37.

Yıldız, T.A., Arshad, S. & Baleanu, D. (2018) New observations on optimal cancer treatments for a fractional tumor growth model with and without singular kernel. Chaos, Solitons & Fractals, 117: 226-239.

Yıldız, T.A., Arshad, S. & Baleanu, D. (2018) Optimal chemotherapy and immunotherapy schedules for a cancer-obesity model with Caputo time fractional derivative. Mathematical Method Applied Science, 41(18): 9390-9407.

Zhou, Y. (2015) Principles and applications of therapeutic ultrasound in healthcare. CRC Press, New York. Pp.102.