Mathematical modeling of cryogenic processes in biotissues and

Mathematical modeling of cryogenic processes in biotissues and optimization of the cryosurgery operations N. A. Kudryashov, K. E. Shilnikov National Research Nuclear University “MEPh. I” Department of Applied Mathematics VIII-th conference Working group on mathematical models and numerical methods in biomathematics Moscow, October 31 – November 3, 2016

Introduction 1. 130% growth of prostate cancer incidence during the decade 2. Wide range of applicability of the cryosurgery 3. Difficulties of cryosurgery planning on the real tumors of an arbitrary shape 4. Finding of the ways to improve the effectiveness and safety of carrying out the cryosurgery 2

Aims and problems solved Aim: The development of effective numerical algorithms and program complexes for the planning and optimization of the cryosurgery operations Problems solved: 1. 2. 3. Multiscale and macroscopically averaged models for bioheat transfer processes taking into account the peculiarities of cell-scale cryogenic processes Numerical algorithms and computerized tools for three dimensional macroscopic modeling of cryotip-based and cutaneous cryosurgery Numerical algorithms for cryosurgery optimization based on the direct operation modeling 3

Introduction. Cell scale. Cryogenic processes in biotissues 1. Intracellular crystallization 2. Cell dehydratation 3. Extracellular crystallization Influence of cooling rate on the cell processes: a) Low cooling rate: extracellular crystallization is dominant (cryoconservation) b) High cooling rate: intracellular crystallization is dominant (cryusurgery) ►►Dependence of freezing interval on the cooling rate 4

Cryogenic processes in biotissues Dehydratation as the mechanism of osmotic regulation Mazur (1965) Cell membrane permeability Levin et. al. (1976) 5

Cryogenic processes in biotissues. Crystallization Intracellular crystallization Toner et al. (1990) Extracellular crystallization Pazhayannur et. al. (1996) 6

Differential model for cell scale crystallization Fraction of unfrozen water Cell phase distribution: + initial conditions 7

Dependence of phase distribution on the cooling rate. Numerical results Cooling rate 80 K/min Cooling rate 5 K/min

Multiscale model for cryogenic bioheat transfer + initial and boundary conditions 9

Mathematical modeling of the cryosurgery with the spherical cryotip. Governing equations 10

Mathematical modeling of the cryosurgery with the spherical cryotip. Initial and boundary conditions Conditions on the cryotip’s working surface : i. e. Consistence condition 11

Mathematical modeling of the cryosurgery with the spherical cryotip Numerical algorithm 1. Explicit finite volume approximation for energy conservation law 2. Semi-implicit Euler-Cauchy scheme for the dynamic system on phase distribution parameters in each grid cell 12

Numerical results Macroscopic mesh: Microscopic time discetization: 13

Method of averaged phase distribution Averaging parameters: 1. - characteristic temperature of full freezing 2. , 14

Mathematical modeling of the cryosurgery with the spherical cryotip Macroscopic statement + initial and boundary conditions 15

Comparison of numerical results obtained by both models proposed with the experimental data Time dependence of freezing front propagation depth Experiment: Junkun et. al. (1999) 16

Three dimensional numerical modeling of the cryosurgery Main assumptions 1. Tumor tissue sizes are much less then the organ’s sizes 2. 3. 4. Computation area allows the construction of indexed mesh Cryotips’ sizes are much less compared to the characteristic computation area sizes Phase averaging parameters are close to apriori estimations 17

Three dimensional numerical modeling of the cryosurgery Problem statement 18

Three dimensional numerical modeling of the cryosurgery Numerical algorithm (1) Shilnikov E. V. (2014) 19

Three dimensional numerical modeling of the cryosurgery Numerical algorithm (2) Flux relaxation method for the improving the scheme stability conditions: Chetverushkin B. N. , Shilnikov E. V. (2010) 20

Three dimensional numerical modeling of the cryosurgery Grid convergence Step, м 5. 44 e-7 3. 41 e-6 0. 00125 5. 27 e-7 3. 28 e-6 0. 001 5. 17 e-7 3. 35 e-6 0. 00083 5. 21 e-7 3. 25 e-6 0. 00071 5. 24 e-7 3. 3 e-6 0. 00063 5. 15 e-7 3. 31 e-6 0. 0005 grid cells Tumor tissue undestructed, м^3 0. 00167 grid cells Healthy tissue destructed, м^3 5. 2 e-7 3. 29 e-6 grid cells 21

Three dimensional numerical modeling of the cryosurgery Numerical results(1) Geometry of the problem Resulting critical iso-surfaces localization relatively to the tumor tissue Propagation of critical iso-surfaces t=30 s t=90 s t=180 s t=285 s 22

Three dimensional numerical modeling of the cryosurgery Numerical results(2) t=15 sec t=90 sec t=150 sec t=300 sec t=450 sec t=600 sec 23

Pareto optimization of the cryosurgery operations Objective variables: Problem statemanent: Objective functions: volume of injured healthy tissue (VH) and volume of undestructed tumor tissue(VT) Quasi-gradient scheme: 24

Optimization of probe based cryosurgery Compromise optimum. Numerical results. t=90 sec t=300 sec t=15 sec t=450 sec t=150 sec t=600 sec 1. Increasing of tumor tissue destruction effectiveness by 23% 2. 3 -fold decrease of healthy tissue injury 25

Optimization of the cutaneous cryosurgery Suppressing heaters method(1) Boundary conditions: Free propagation of the tissue necrosis front 26

Optimization of the cutaneous cryosurgery Suppressing heaters method(2) Resulting tissue necrosis front localization H ealthy tissue injury decreased by 30% against the increasing of duration of the operation by 25% 27

Conclusion 1. The multiscale and phase averaged macroscopic models for the cryogenic bioheat transfer are proposed as the tools for predicting the results of cryosurgery with respect to peculiarities of phase change processes in biological solutes. 2. The quasi-gradient method is adopted for the optimization of the cryosurgery operations 3. Numerical experiments on the model problems showed that the developed planning tools can cause a valuable increasing of effectiveness and safety of cryosurgery operations performing 28

Publications 1. N. A. Kudryashov, K. E. Shilnikov. Numerical modeling and optimization of the cryosurgery operations // Journal of Computational and Applied Mathematics. 2015, No. 290, pp. 259– 267. 2. N. A. Kudryashov, K. E. Shilnikov. Numerical simulation of cryosurgeries and optimization of probe placement // Computational Mathematics and Mathematical Physics. Volume 55, Issue 9, pp. 1579– 1589. 3. N. A. Kudryashov, K. E. Shilnikov. Numerical Simulation of the Effect of Localization of the Front of Cellular Necrosis during Cutaneous Cryosurgery // Mathematical Models and Computer Simulations. 2016, Vol. 8, No. 6, pp. 680– 688. 4. N. A. Kudryashov, K. E. Shilnikov. Nonlinear bioheat transfer models and multiobjective numerical optimization of the cryosurgery operations // Proceedings of 13 -th International Conference of Numerical Analysis and Applied Mathematics, Greece, 2015. AIP Conf. Proc. 1738, 230009 (2016); http: //dx. doi. org/10. 1063/1. 4952018 29

Thank you for your attention 30