Modeling of soft tissue thermal damage based on GPU acceleration

Jinao Zhang, Jeremy Hills, Yongmin Zhong, Bijan Shirinzadeh, Julian A Smith, Chengfan Gu

Research output: Contribution to journalArticleResearchpeer-review

1 Citation (Scopus)

Abstract

Hyperthermia treatments require precise control of thermal energy to form the coagulation zones which sufficiently cover the tumor without affecting surrounding healthy tissues. This has led modeling of soft tissue thermal damage to become important in hyperthermia treatments to completely eradicate tumors without inducing tissue damage to surrounding healthy tissues. This paper presents a methodology based on GPU acceleration for modeling and analysis of bio-heat conduction and associated thermal-induced tissue damage for prediction of soft tissue damage in thermal ablation, which is a typical hyperthermia therapy. The proposed methodology combines the Arrhenius Burn integration with Pennes' bio-heat transfer for prediction of temperature field and thermal damage in soft tissues. The problem domain is spatially discretized on 3-D linear tetrahedral meshes by the Galerkin finite element method and temporally discretized by the explicit forward finite difference method. To address the expensive computation load involved in the finite element method, GPU acceleration is implemented using the High-Level Shader Language and achieved via a sequential execution of compute shaders in the GPU rendering pipeline. Simulations on a cube-shape specimen and comparison analysis with standalone CPU execution were conducted, demonstrating the proposed GPU-accelerated finite element method can effectively predict the temperature distribution and associated thermal damage in real time. Results show that the peak temperature is achieved at the heat source point and the variation of temperature is mainly dominated in its direct neighbourhood. It is also found that by the continuous application of point-source heat energy, the tissue at the heat source point is quickly necrotized in a matter of seconds, while the entire neighbouring tissues are fully necrotized in several minutes. Further, the proposed GPU acceleration significantly improves the computational performance for soft tissue thermal damage prediction, leading to a maximum reduction of 55.3 times in computation time comparing to standalone CPU execution.

Original languageEnglish
Pages (from-to)5-12
Number of pages8
JournalComputer Assisted Surgery
Volume24
Issue numberS1
DOIs
Publication statusPublished - 1 Oct 2019

Keywords

  • Arrhenius Burn integration
  • finite element method
  • Graphics Processing Unit
  • Pennes’ bio-heat equation
  • thermal ablation
  • Thermal damage

Cite this

@article{2746c19dcf474a658e5657a488186aad,
title = "Modeling of soft tissue thermal damage based on GPU acceleration",
abstract = "Hyperthermia treatments require precise control of thermal energy to form the coagulation zones which sufficiently cover the tumor without affecting surrounding healthy tissues. This has led modeling of soft tissue thermal damage to become important in hyperthermia treatments to completely eradicate tumors without inducing tissue damage to surrounding healthy tissues. This paper presents a methodology based on GPU acceleration for modeling and analysis of bio-heat conduction and associated thermal-induced tissue damage for prediction of soft tissue damage in thermal ablation, which is a typical hyperthermia therapy. The proposed methodology combines the Arrhenius Burn integration with Pennes' bio-heat transfer for prediction of temperature field and thermal damage in soft tissues. The problem domain is spatially discretized on 3-D linear tetrahedral meshes by the Galerkin finite element method and temporally discretized by the explicit forward finite difference method. To address the expensive computation load involved in the finite element method, GPU acceleration is implemented using the High-Level Shader Language and achieved via a sequential execution of compute shaders in the GPU rendering pipeline. Simulations on a cube-shape specimen and comparison analysis with standalone CPU execution were conducted, demonstrating the proposed GPU-accelerated finite element method can effectively predict the temperature distribution and associated thermal damage in real time. Results show that the peak temperature is achieved at the heat source point and the variation of temperature is mainly dominated in its direct neighbourhood. It is also found that by the continuous application of point-source heat energy, the tissue at the heat source point is quickly necrotized in a matter of seconds, while the entire neighbouring tissues are fully necrotized in several minutes. Further, the proposed GPU acceleration significantly improves the computational performance for soft tissue thermal damage prediction, leading to a maximum reduction of 55.3 times in computation time comparing to standalone CPU execution.",
keywords = "Arrhenius Burn integration, finite element method, Graphics Processing Unit, Pennes’ bio-heat equation, thermal ablation, Thermal damage",
author = "Jinao Zhang and Jeremy Hills and Yongmin Zhong and Bijan Shirinzadeh and Smith, {Julian A} and Chengfan Gu",
year = "2019",
month = "10",
day = "1",
doi = "10.1080/24699322.2018.1557891",
language = "English",
volume = "24",
pages = "5--12",
journal = "Computer Assisted Surgery",
issn = "2469-9322",
publisher = "Taylor & Francis",
number = "S1",

}

Modeling of soft tissue thermal damage based on GPU acceleration. / Zhang, Jinao; Hills, Jeremy; Zhong, Yongmin; Shirinzadeh, Bijan; Smith, Julian A; Gu, Chengfan.

In: Computer Assisted Surgery, Vol. 24, No. S1, 01.10.2019, p. 5-12.

Research output: Contribution to journalArticleResearchpeer-review

TY - JOUR

T1 - Modeling of soft tissue thermal damage based on GPU acceleration

AU - Zhang, Jinao

AU - Hills, Jeremy

AU - Zhong, Yongmin

AU - Shirinzadeh, Bijan

AU - Smith, Julian A

AU - Gu, Chengfan

PY - 2019/10/1

Y1 - 2019/10/1

N2 - Hyperthermia treatments require precise control of thermal energy to form the coagulation zones which sufficiently cover the tumor without affecting surrounding healthy tissues. This has led modeling of soft tissue thermal damage to become important in hyperthermia treatments to completely eradicate tumors without inducing tissue damage to surrounding healthy tissues. This paper presents a methodology based on GPU acceleration for modeling and analysis of bio-heat conduction and associated thermal-induced tissue damage for prediction of soft tissue damage in thermal ablation, which is a typical hyperthermia therapy. The proposed methodology combines the Arrhenius Burn integration with Pennes' bio-heat transfer for prediction of temperature field and thermal damage in soft tissues. The problem domain is spatially discretized on 3-D linear tetrahedral meshes by the Galerkin finite element method and temporally discretized by the explicit forward finite difference method. To address the expensive computation load involved in the finite element method, GPU acceleration is implemented using the High-Level Shader Language and achieved via a sequential execution of compute shaders in the GPU rendering pipeline. Simulations on a cube-shape specimen and comparison analysis with standalone CPU execution were conducted, demonstrating the proposed GPU-accelerated finite element method can effectively predict the temperature distribution and associated thermal damage in real time. Results show that the peak temperature is achieved at the heat source point and the variation of temperature is mainly dominated in its direct neighbourhood. It is also found that by the continuous application of point-source heat energy, the tissue at the heat source point is quickly necrotized in a matter of seconds, while the entire neighbouring tissues are fully necrotized in several minutes. Further, the proposed GPU acceleration significantly improves the computational performance for soft tissue thermal damage prediction, leading to a maximum reduction of 55.3 times in computation time comparing to standalone CPU execution.

AB - Hyperthermia treatments require precise control of thermal energy to form the coagulation zones which sufficiently cover the tumor without affecting surrounding healthy tissues. This has led modeling of soft tissue thermal damage to become important in hyperthermia treatments to completely eradicate tumors without inducing tissue damage to surrounding healthy tissues. This paper presents a methodology based on GPU acceleration for modeling and analysis of bio-heat conduction and associated thermal-induced tissue damage for prediction of soft tissue damage in thermal ablation, which is a typical hyperthermia therapy. The proposed methodology combines the Arrhenius Burn integration with Pennes' bio-heat transfer for prediction of temperature field and thermal damage in soft tissues. The problem domain is spatially discretized on 3-D linear tetrahedral meshes by the Galerkin finite element method and temporally discretized by the explicit forward finite difference method. To address the expensive computation load involved in the finite element method, GPU acceleration is implemented using the High-Level Shader Language and achieved via a sequential execution of compute shaders in the GPU rendering pipeline. Simulations on a cube-shape specimen and comparison analysis with standalone CPU execution were conducted, demonstrating the proposed GPU-accelerated finite element method can effectively predict the temperature distribution and associated thermal damage in real time. Results show that the peak temperature is achieved at the heat source point and the variation of temperature is mainly dominated in its direct neighbourhood. It is also found that by the continuous application of point-source heat energy, the tissue at the heat source point is quickly necrotized in a matter of seconds, while the entire neighbouring tissues are fully necrotized in several minutes. Further, the proposed GPU acceleration significantly improves the computational performance for soft tissue thermal damage prediction, leading to a maximum reduction of 55.3 times in computation time comparing to standalone CPU execution.

KW - Arrhenius Burn integration

KW - finite element method

KW - Graphics Processing Unit

KW - Pennes’ bio-heat equation

KW - thermal ablation

KW - Thermal damage

UR - http://www.scopus.com/inward/record.url?scp=85072234444&partnerID=8YFLogxK

U2 - 10.1080/24699322.2018.1557891

DO - 10.1080/24699322.2018.1557891

M3 - Article

C2 - 31340685

AN - SCOPUS:85072234444

VL - 24

SP - 5

EP - 12

JO - Computer Assisted Surgery

JF - Computer Assisted Surgery

SN - 2469-9322

IS - S1

ER -