Solving two-dimensional bioheat optimal control problem in solid and vessel domain by pseudospectral discretization

Abstract

In this paper, an optimal control problem subject to an energy transfer equation in thermally significant blood vessel coupled with the Pennes governing equation is solved. The proposed dynamical system is modeled by a system of conduction-reaction and convection-conduction equations. Chebyshev–Gauss–Lobatto (CGL) and Legendre–Gauss–Lobatto (LGL) points are applied to discretize optimal control problem. Then, the problem transferred into a quadratic programming problem using domain decomposition and matching the solution and its first derivative across the related interfaces. The Tikhonov regularization method is used to compute the optimal regularization parameter for the objective function. Some theories have been discussed to prove that the discrete problem has a solution and it is also consistent with the continuous problem. Numerical results are compared and numerical convergent is shown. Optimal and non-optimal external heat source is compared and it is shown that the mathematical modeling is effective in saving the external heat energy. Numerical results show that by this process of heating at the final time, the temperature is very close to the desired temperature without damaging the vessel domain. Numerical results and corresponding graphs are depicted.

Keywords

Pennes’ bioheat equation optimal control blood vessel pseudospectral method quadratic programming Tikhonov regularization

Get full access to this article

View all access options for this article.

References

Abbasi

Malek

(2019) Hyperthermia cancer therapy by domain decomposition methods using strongly continuous semigroups. Mathematics and Computers in Simulation 165: 1–12.

Bagaria

Johnson

(2005) Transient solution to the bioheat equation and optimization for magnetic fluid hyperthermia treatment. International Journal of Hyperthermia 21(1): 57–75.

Bazaraa

Sherali

Shetty

(2013) Nonlinear Programming: Theory and Algorithms. Hoboken, NJ: John Wiley & Sons.

Betts

Campbell

Thompson

(2016) Solving optimal control problems with control delays using direct transcription. Applied Numerical Mathematics 108: 185–203.

Canuto

Quarteroni

(1982) Approximation results for orthogonal polynomials in Sobolev spaces. Mathematics of Computation 38(157): 67–86.

Canuto

Hussaini

Quarteroni

, et al. (2007) Spectral Methods: Fundamentals in Single Domains. Berlin: Springer Science & Business Media.

Cao

Lesnic

(2019) Simultaneous reconstruction of the perfusion coefficient and initial temperature from time-average integral temperature measurements. Applied Mathematical Modelling 68: 523–539.

Dehghan

Taleei

(2011) A Chebyshev pseudospectral multidomain method for the soliton solution of coupled nonlinear Schrödinger equations. Computer Physics Communications 182(12): 2519–2529.

Gänzler

Volkwein

Weiser

(2006) SQP methods for parameter identification problems arising in hyperthermia. Optimization Methods and Software 21(6): 869–887.

10.

Hansen

(1994) Regularization tools: A Matlab package for analysis and solution of discrete ill-posed problems. Numerical Algorithms 6(1): 1–35.

11.

Horng

Lin

Liauh

, et al. (2007) Effects of pulsatile blood flow in large vessels on thermal dose distribution during thermal therapy. Medical Physics 34(4): 1312–1320.

12.

Huang

Horng

(2015) Bioheat transfer and thermal heating for tumor treatment. In: Becker

Kuznetsov

(eds) Heat Transfer and Fluid Flow in Biological Processes. London: Academic Press, pp. 1–42.

13.

Kumar

Rai

(2017) Numerical solution of non-linear dual-phase-lag bioheat transfer equation within skin tissues. Mathematical Biosciences 293: 56–63.

14.

Kumar

Rai

(2016) Numerical simulation of dual-phase-lag bioheat transfer model during thermal therapy. Mathematical Biosciences 281: 82–91.

15.

Qin

(2017) Multidomain Legendre–Galerkin Chebyshev-collocation method for one-dimensional evolution equations with discontinuity. Applied Numerical Mathematics 111: 246–259.

16.

Malek

Abbasi

(2014) Optimal control solution for Pennes’ equation using strongly continuous semigroup. Kybernetika 50(4): 530–543.

17.

Malek

Abbasi

(2015) Heat treatment modelling using strongly continuous semigroups. Computers in Biology and Medicine 62: 65–75.

18.

Malek

Abbasi

(2016) Optimal control for Pennes’ bioheat equation. Asian Journal of Control 18(2): 674–685.

19.

Malek

Momeni-Masuleh

(2008) A mixed collocation–finite difference method for 3d microscopic heat transport problems. Journal of Computational and Applied Mathematics 217(1): 137–147.

20.

Malinen

Duncan

Huttunen

, et al. (2006) Feedforward and feedback control of ultrasound surgery. Applied Numerical Mathematics 56(1): 55–79.

21.

Mehrali-Varjani

Shamsi

Malek

(2018) Solving a class of Hamilton–Jacobi–Bellman equations using pseudospectral methods. Kybernetika 54(4): 629–647.

22.

Meyer

Philip

Tröltzsch

(2006) Optimal control of a semilinear PDE with nonlocal radiation interface conditions. SIAM Journal on Control and Optimization (SICON) 45(2): 699–721.

23.

Prakash

(2010) Theoretical modeling for hepatic microwave ablation. The Open Biomedical Engineering Journal 4: 27–38.

24.

Shih

Horng

Huang

, et al. (2012) Numerical analysis of coupled effects of pulsatile blood flow and thermal relaxation time during thermal therapy. International Journal of Heat and Mass Transfer 55(13–14): 3763–3773.

25.

Shih

Kou

Lin

(2003) The impact of thermally significant blood vessels in perfused tumor tissue on thermal dose distributions during thermal therapies. International Communications in Heat and Mass Transfer 30(7): 975–985.

26.

Shih

Liu

Horng

ATL

(2006) Cooling effect of thermally significant blood vessels in perfused tumor tissue during thermal therapy. International Communications in Heat and Mass Transfer 33(2): 135–141.