Locally implicit discontinuous Galerkin time domain method for electromagnetic wave propagation in dispersive media applied to numerical dosimetry in biological tissues

Abstract : We are concerned here with the numerical simulation of electromagnetic wave propagation in biological media. Because of their water content, these media are dispersive i.e. their electromagnetic material characteristics depend of the frequency. In the time-domain, this translates in a time dependency of these parameters that can be taken into account through and additional (auxiliary) differential equation for, e.g, the electric polarization, which is coupled to the system of Maxwell's equations. From the application point of view, the problems at hand most often involve irregularly shaped structures corresponding to biological tissues. Modeling realistically the interfaces between tissues is particularly important if one is interested in evaluating accurately the impact of field discontinuities at these interfaces. In this paper, we propose and study a locally implicit discon-tinuous Galerkin time-domain method formulated on an unstructured tetrahedral mesh for solving the resulting system of differential equations in the case of Debye-type media. Three-dimensional numerical simulations are presented concerning the exposure of head tissues to a localized source radiation. 1. Introduction. This article is concerned with the numerical simulation of electromagnetic wave propagation in dispersive media. These are materials in which either or both of the electromagnetic material parameters ε and µ are functions of frequency. Note that the conductivity σ may also be a function of frequency, but its effect can be rolled into the complex permittivity. In reality, all materials have frequency-dependent ε and µ, but many materials can be approximated as frequency-independent over a frequency band of interest, simplifying their analysis and simulation. Here, we will focus on the much more common case of frequency-dependent permittivity. A lot of practical electromagnetic wave propagation problems involve such propagation media, such as those involving the interaction of an electromagnetic wave with biological tissues. The numerical modeling of the propagation of electromagnetic waves through human tissues is at the heart of many biomedical applications such as the microwave imaging of cancer tumors or the treatment of the latter by hy-perthermia. For example, microwave imaging for breast cancer detection is expected to be safe for the patient and has the potential to detect very small cancerous tumors in the breast [3, 14, 24]. Beside, the definition of microwave-based hyperthermia as an immunotherapy strategy for cancer can also be cited [2, 11]. The electroporation technique can also be an application, which consists of applying nanopulses to the tissues, enabling only intracellular membranes to be affected, and then opening the route to therapeutic strategies such as electrochemotherapy or gene transfer [21, 30, 25, 26, 28]. Because for all these biomedical applications experimental modeling is almost impossible , computer simulation is the approach of choice for understanding the underlying
Complete list of metadatas

Cited literature [41 references]  Display  Hide  Download

https://hal.inria.fr/hal-01133694
Contributor : Stéphane Descombes <>
Submitted on : Friday, September 23, 2016 - 9:43:45 PM
Last modification on : Wednesday, December 4, 2019 - 10:03:07 AM

File

sisc-article-lastversion.pdf
Files produced by the author(s)

Identifiers

Citation

Stéphane Descombes, Stéphane Lanteri, Ludovic Moya. Locally implicit discontinuous Galerkin time domain method for electromagnetic wave propagation in dispersive media applied to numerical dosimetry in biological tissues. SIAM Journal on Scientific Computing, Society for Industrial and Applied Mathematics, 2016, 38 (5), pp.A2611-A2633. ⟨10.1137/15M1010282⟩. ⟨hal-01133694v2⟩

Share

Metrics

Record views

296

Files downloads

290