Logo image
Back
A Flux-Conservative finite difference scheme for the numerical solution of the nonlinear bioheat equation
Book chapter

A Flux-Conservative finite difference scheme for the numerical solution of the nonlinear bioheat equation

G.C. Bourantas, G.R. Joldes, A. Wittek and K. Miller
Computational Biomechanics for Medicine, pp.69-81
Springer
2018
DOI: https://doi.org/10.1007/978-3-319-75589-2_7
url
Link to Published Version *Subscription may be requiredView

Abstract

We present a flux-conservative finite difference (FCFD) scheme for solving the nonlinear (bio)heat transfer in living tissue. The proposed scheme deals with steep gradients in the material properties for malignant and healthy tissues. The method applies directly on the raw medical image data without the need for sophisticated image analysis algorithms to define the interface between tumor and healthy tissues. We extend the classical finite difference (FD) method to cases with high discontinuities in the material properties. We apply meshless kernels, widely used in Smoothed Particle Hydrodynamics (SPH) method, to approximate properties in the off-grid points introduced by the flux-conservative differential operators. The meshless kernels can accurately capture the steep gradients and provide accurate approximations. We solve the governing equations by using an explicit solver. The relatively small time-step applied is counterbalanced by the small computation effort required at each time-step of the proposed scheme. The FCFD method can accurately compute the numerical solution of the bioheat equation even when noise from the image acquisition is present. Results highlight the applicability of the method and its ability to solve tumor ablation simulations directly on the raw image data, without the need to define the interface between malignant and healthy tissues (segmentation) or meshing.

Details

Metrics

54 Record Views
Logo image