Next Article in Journal
Insights into the Functions of LncRNAs in Drosophila
Next Article in Special Issue
Assembling Magnetic Nanoparticles on Nanomechanical Resonators for Torque Magnetometry
Previous Article in Journal / Special Issue
Instrument-Free and Visual Detection of Salmonella Based on Magnetic Nanoparticles and an Antibody Probe Immunosensor
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Numerical Model for Magnetic Fluid Hyperthermia in a Realistic Breast Phantom: Calorimetric Calibration and Treatment Planning

by
Arkadiusz Miaskowski
1 and
Mahendran Subramanian
2,*
1
Department of Applied Mathematics and Computer Science, University of Life Sciences Lublin, 20-950 Lublin, Poland
2
Department of Bioengineering and Department of Computing, Royal School of Mines, Imperial College London, London SW7 2AZ, UK
*
Author to whom correspondence should be addressed.
Int. J. Mol. Sci. 2019, 20(18), 4644; https://doi.org/10.3390/ijms20184644
Submission received: 9 August 2019 / Revised: 13 September 2019 / Accepted: 17 September 2019 / Published: 19 September 2019
(This article belongs to the Special Issue Magnetic Nanoparticles 2019)

Abstract

:
This paper aims to apply a proposed, based on calorimetric measurements, a reliable numerical model for magnetic fluid hyperthermia (MFH) treatment planning of breast cancer. Furthermore, we perform a comparative analysis of magnetic nanoparticles (MNPs) and tumour tissue interactions by means of the magnetic-field-dependent Néel and Brownian relaxation times. The analysis was based on an anatomically correct breast model (developed in-house) and a modified linear response theory, which was applied to investigate the heat dissipation from the magnetic nanoparticles dispersed in the breast tumour. The calculations of the single-domain magnetic power losses were conducted for a case where the magnetic field value and the applied frequency were known, but also for the different concentrations of the MNPs in the tumour. Two scenarios were considered: The MNPs mobilised and immobilised in the tumour. In parallel, the eddy currents effect, together with the related temperature distributions, were considered in order to analyse safety issues. By changing the MNP concentration in the tumour, the corresponding temperature distributions were calculated. The eddy current effect, together with the related temperature distribution, were considered in order to analyse safety issues. Varying the MNP concentration in the tumour, the corresponding temperature distribution was calculated. Moreover, the cumulative equivalent minutes at 43   were analysed. In the anatomically correct breast phantoms, the tissue location can lead to “hot spots” due to the eddy current effect and subsequently to the high gradients of the temperature. That is why the analysis of safety issues related to the overheating side effect should be taken into consideration during the treatment planning of magnetic fluid hyperthermia. The phenomenon of heat dissipation from MNPs is very sophisticated and depends on their concentration, the distribution and the relaxation mechanism in the tumour, together with magnetic field strength and frequency. Furthermore, we inferred that the phenomenon of heat dissipation from MNPs equally depends on MNP-tissue interactions, and it can lead to 30% differences in the power assessment. Nevertheless, the aforementioned factors should be considered in parallel using anatomical, volume-dependent models to enhance the efficiency of in vivo treatment.

1. Introduction

Novel therapy methods without negative impacts are needed to reduce mortality rates in patients with breast cancer. One such method is magnetic fluid hyperthermia (MFH) in which magnetic nanoparticles (MNPs) dissipate heat during radiofrequency (RF) application [1,2,3,4,5], activating biochemical pathways leading to necrosis or apoptosis [6,7,8,9]. In the past decade, attention has been given to MFH treatment for breast cancer either on its own or in combination with chemotherapy and radiotherapy [10,11,12,13,14]. However, dexterous techniques of modelling and validation will help to improve and standardise this potential clinical treatment. Effective MFH is accomplished when maximum cytotoxicity with low or no side effects are attained with a limited dosage of Alternating Magnetic Field (AMF) susceptible MNPs, so suitable heat characterisation of MNPs and treatment planning should be carried out prior to therapy. Besides, ethical factors, cost-effective clinical trials and infrastructure requirements restrict the number of teams performing in vivo studies, hampering significant progress in understanding MFH for the therapy of breast cancer. Furthermore, constraints such as those for the pain threshold of 4.85 × 108 A/m/s for induced eddy currents in patients during MFH treatment [9] and the field limit being sufficient to saturate all the particles within the distribution at a particular frequency, limit MFH treatment. However, it is possible to use virtual breast phantom models to evaluate the efficacy of MFH and the risk assessment of the radiofrequency.
Three dimensional (3-D) models were used in the past to investigate the effects of MNP volume fraction, nanocomposite geometry, and treatment parameters on thermal profiles regarding female breast cancer [15,16,17]. However, the proposed models were simplified to the semi-ellipsoidal block, the semi-sphere or to the anatomically correct shape without naturalistic breast tissues. This study aims to develop an improved numerical breast phantom model derived from human magnetic resonance images [18] to obtain a lifelike model. Moreover, the literature lacks numerical models, which include 3-D temperature distribution due to eddy currents and single-domain magnetic power losses in anatomically correct phantoms. Furthermore, the numerical investigation conducted in this work was based on the results received from calorimetric measurements of the ferrofluid using the magneThermTM system [19,20], i.e., the numerical model developed for heat generation terms. Specific loss powers (SLP) were validated with calorimetric measurements and are now used to analyse breast cancer hyperthermia treatment mediated by magnetic fluid.
It is clear that the reported SLPs values [20] do not reflect the situation one might expect to find in human tissues where, for example, blood perfusion, spatial tissues location and magnetic nanoparticles–tissue interactions affect heat dissipation. That is why, in this work, the reported SLPs values together with the numerical model, have been transposed to realistic breast cancer.
Our previous studies presented mathematical models, i.e., the single-domain magnetic power losses model and the eddy current model. Both of them were used later in the modified Pennes equation to evaluate a temperature distribution. The current work details the anatomically correct breast model used in the simulations and its parameters, followed by several analyses of volumetric power losses, which determine spatial temperature distribution as a function of time for various concentrations of magnetic fluid preparations. Furthermore, mobilised and immobilised MNPs in the tumour were considered.

2. Results and Discussion

The mathematical model described in the previous section was used to investigate the temperature distribution during MFH. Different magnetic fluid concentrations were applied in the tumour to control the temperature in the target tissue, then the tumour position in the breast model was determined by maximum electrical loss. These power density values are referred to as “hot spots”. In addition to the power losses due to the eddy current effect, a power dissipation from magnetic nanoparticles should be considered in parallel when dealing with MFH treatment planning. Moreover, particular interest was paid to the comparative analysis of MNPs–tissue interactions, as in MFH, and the heating rates of MNPs which were determined using Néel and Brownian losses, which can be inappropriate for MFH treatment. For a reliable model of heating, two cases were considered, i.e., the MNPs free to move in the tumour (Néel and Brownian relaxations included in the tumour) and the MNPs immobilised in the tumour (in which only Néel relaxation takes place).

2.1. Breast Phantom Model

In our previous study, the breast model containing only skin, breast fat and tumour tissues was evaluated based on field application by a pancake coil [17]. In the present study, more tissues were added to achieve a very detailed anatomical structure. The breast phantom model was taken from [18], and it was assigned to the Class 3 category—heterogeneously dense (51–75% glandular). The following breast tissues were included in the model: Breast gland (blue), breast fat (yellow), fat (green), skin (red) and muscle (orange) as shown in Figure 1. Moreover, the physical parameters of the tumour were considered as the same as the muscle tissues ones but with perfusion. The physical parameters of tissues were calculated for the frequency f = 171 kHz based on the Information Technologies in Society (IT’IS) database [21]. The parameters are shown in Table 1, where HTR = ρ b c b ρ ω is the heat transfer rate.
In our previous work [17], the tissues included in the breast phantom models were averaged due to their volume fraction together with their physical parameters. Such approximation can be useful while planning the treatment when the applicator position pattern is required. In the current work, the model was improved, and the effect of the eddy current and “hot spots” (the places where high values of power losses may occur due to the anatomical configuration of breast tissues) were examined. The model cross-section is shown in Figure 1.

2.2. Eddy Currents Effect on Temperature Distribution

The authors would like to emphasise that the aim of this part of the study is not the design a breast cancer applicator, but to fulfil the condition for the magnetic field strength ( H max = 10.3 kA / m in the middle of the applicator) and the frequency, as was performed during the calorimetric experiments using the magneThermTM system. The location of the spherical tumour (2 mL volume) in the breast model was determined by the maximum volumetric eddy power density ( P eddy ) to consider the worst-case scenario.
The volumetric power density due to the eddy currents was evaluated for frequency   f 1 = 171   kHz , and is shown in Figure 2. It is obvious that the highest value of eddy power density occurred on the skin tissue, reaching 70 kW/m3. That is why the data was rescaled to the maximum value (35.5 kW/m3), i.e., to one of the “hot spots”, which were placed in the model excluding the skin layer to improve visualisation. It is important to note that the presence of “hot spots” should be considered during treatment planning as they are hard to predict and they can limit the amount of energy that can be deposited in a patient’s body. Moreover, one should remember that maximal electrical power losses will occur on the surface of the tumour as the eddy currents have to circulate around the tumour centre where they are the smallest. This sometimes can be used as an advantage because high-conductivity tumours surrounded by low-conductivity tissues will have a local eddy current flowing around the approximate centre of the tumour leading to increased heating of deep-seated tumours.
In our case, to investigate the worst-case scenario, the tumour was positioned in one of the “hot spots”, and the power density was calculated (see Figure 2—right). One can see that the tumour position has changed the power distribution pattern and the “hot spot” shifted due to the physical properties of the tumour. However, taking the temperature distribution into account (see Figure 3), the area where the temperature rises above 39 °C (which is relatively high), means that overheating of tissues can occur primarily on the skin.

2.3. Single-Domain Magnetic Power Loss Analysis

As the first step, the single-domain magnetic power losses for both the mobilised and the immobilised MNPs were calculated. The ferrofluid concentration was the same as what is was during the calorimetry experiments ( φ = 5.0   mg / mL ) , and it was used as a reference. The parameters which have been used in Formula (1) are provided in Table 2.
Next, the ferrofluid concentration was changed ( ϕ r ) to control the temperature in the tumour. The values of maximum magnetic power losses for the different relative concentrations ϕ r are given in Table 3 and also presented in Figure 4.
Taking into account the single-domain magnetic power losses, one can see that the powers had decreased by 30% when MNPs were immobilised in the tumour. It is also easy to linearly approximate both cases as follows: P mobilzed = 244.4 ϕ r 0.7826 , P immoblized = 170.3 ϕ r 0.2539 in order to predict the magnetic power losses for other magnetic fluid concentrations.
For both cases (the mobilised MNPs on the left and the immobilised MNPs on the right) in Figure 5, the temperature distribution patterns are shown in the breast phantom for different concentrations of the ferrofluid after 1800 secs of heating. The black area in Figure 5A, 5B, 5C indicates the temperature of 40 °C for different ϕ r   = [1, 0.5, 0.35]. Moreover, the temperatures over time in the middle of the tumour are presented in Figure 6 for the mobilised case.
In the case of the immobilised MNPs in the tumour, it can be seen that both the maximum power dissipation (see Table 3), as well as the temperature, has reduced (see Figure 6). For the reference concentration, ϕ r = 1.0 , the power decreased by 30%, while the temperature difference, Δ T , was approximately 5.4   when comparing the two cases of MNPs in the tumour. Nevertheless, the relationship between MNPs concentration and the temperature increase is not straightforward, for example, for ϕ r = 0.5 the temperature difference is Δ T = 2.1   and for ϕ r = 0.35 it is Δ T = 1.2   . However, in order to compare the aforementioned cases, the concentration, ϕ r , was changed for the corresponding values of the power dissipation (for example, to receive P moblil = 121.0   kW / m 3 , the concentration should be ϕ r = 0.5 for the mobilised MNPs, while for immobilised MNPs concentration should be ϕ r = 0.71 , see Table 3). Finally, Figure 7 presents the relative concentration as the function of the temperature difference between the two cases together with its linear approximation.
According to the data in Table 3 and Figure 5, it is possible to control the temperature in the tumour and its vicinity by changing the magnetic fluid concentration. The volumetric power density, required to raise the temperature to 43   in the tumour was at least 84.2   kW / m 3 , as depicted in Figure 8, where the isosurface 43   can be seen inside the tumour volume.
However, in order to reach 43   in the whole tumour volume 121.0   kW / m 3 was required but, in this case, the temperature inside the tumour was 45   . On the other hand, it should be underlined that the heating power may have increased due to a possible agglomeration of MNPs in the tumour. Moreover, the heterogeneous temperature distribution in the tumour vicinity, which can be observed in Figure 4, Figure 6 and Figure 9 may lead to high-temperature gradients even if the eddy currents power losses are minimal when compared with the magnetic power ones. This phenomenon should be considered in patient treatment planning.

3. Materials and Methods

3.1. Single-Domain Magnetic Power Losses Model

When polydisperse magnetic nanoparticles in the superparamagnetic regime are exposed to an AMF with the given magnetic field strength, H max , and frequency, f, the magnetisation lags behind the external field. According to linear response theory (LRT), the heat dissipation in W/kg can be expressed as the specific loss power (SLP) [23]:
S L P L R T = 0 V π μ 0 χ H max 2 f ρ V M g ( μ , s , V M ) d V     ( W kg )
where g ( μ ,   s , V M ) is the volume-weighted distribution (log-normal distribution), ρ and V M are the density and the magnetic volume of MNP, respectively, and χ is the average out-of-phase component of the susceptibility given by:
χ = μ 0 m 3 k B T   2 π f τ ( 1 + ( 2 π f τ ) 2 )
where k B is the Boltzmann constant, m = M s V M is the magnetic moment where M s is the magnetic saturation and τ is the Néel-Brown relaxation time. The Néel ( τ N ) and Brownian ( τ B ) metrics, as dependent on an external magnetic field (Hmax), are given by [23]:
[ τ N ( H max     ) ] 2 = f 0 ( 1 h 2 ) { ( 1 + h ) exp [ ( K V k B T ) ( 1 + h ) 2 ] + + ( 1 h ) exp [ ( K V k B T ) ( 1 h ) 2   ] }
[ τ B ( H max ) ] 1 = τ B 1 ( 1 + 0.07 ( μ 0 μ H max k B T ) 2 ) 0.5
where τ B = 3 η V H k B T and η is the viscosity, h = H max H k   and μ 0 H k = 2 K M s represents the anisotropy field and f 0 = 10 9   s 1 is the attempt frequency. Finally, the existence of two relaxation times leads to the effective relaxation time given by τ ( H max ) 1 = τ N ( H max ) 1 + τ B ( H max ) 1 . In this case, it is worth mentioning, that the above model contains few simplifications such as follows: All the MNPs have spherical shapes, MNPs do not create clusters or chains, and therefore there are not any inter-particle interactions and finally—the hydrodynamic particle size distribution is uniform. One should realise that if the concentration of the MNPs in the ferrofluid is relatively small and, for example, does not exceed 5   mg / mL , the effects caused by dipole–dipole interactions between the nanoparticles are considered insignificant at room temperature. Thus, in our case, the nanoparticles were considered to be magnetically isolated from each other.

3.2. Eddy Current Effect Model

When human tissues are exposed to an alternating electromagnetic field, the eddy current effect is observed due to non-zero conductivity of the tissues, which finally leads to their heating.
The magnetic field strength ( H ) generated by an applicator can be solved using the magnetic vector potential ( A ):
2 A = μ 0 J s
where J s is the current density vector (A/m2) and H = μ 0 1 × A . In order to include a volumetric power density ( P eddy ) produced by an eddy current ( J eddy ) effect, a magneto-quasi static algorithm was applied. Considering the scalar potential (u) together with a magnetic vector potential (A), one can formulate the issue as:
J e d d y = σ u σ ω A
where σ is the conductivity (S/m) and ω = 2πf is the angular frequency. Finally, P eddy can be expressed as:
P eddy = J eddy 2 σ
Discussing the eddy current or Foucault current is essential due to their effect on different types of tissues when exposed to a time-varying magnetic field. The eddy current effect can lead to unwanted non-specific heating of healthy tissues, so it must be considered in the applicator design and during treatment planning.

3.3. Pennes Equation

A model dedicated to the numerical analysis of heat transfer in human tissues was proposed by Pennes [24]. In this model, blood perfusion is assumed to be uniform throughout the tissue and all heat leaving the artery is absorbed by the local tissue, with no venous rewarming. Moreover, the Pennes model is limited to the perfusion source term, i.e., the arterial temperature is assumed to be equal to the body core temperature.
To account for the strong temperature dependence due to bio-regulatory processes [25], the breast model tissues perfusion rates were assumed to be constant ( ω ( T ) = c o n s t . ) while the tumour perfusion was assumed nonlinear according to the following formula [26]:
ω tumour ( T ) = 0.416 + 0.416 exp ( ( T 37   ) 4 220 )         ( kg m 3 s )
which indicates the mass flow rate as presented in Figure 9.
Hence, the temperature distribution in the breast model was investigated using the modified bio-heat transfer equation:
ρ c T t = ( κ T ) + ρ b c b ρ ω ( T ) ( T b T ) + ρ HGR + ρ S L P LRT + P eddy     ( W m 3 )
where ρ is the tissue density, c is the tissue-specific heat, ρ b is the density of blood, c b is blood specific heat, κ is the tissue thermal conductivity, ω(T) is the perfusion rate, T b is the arterial blood temperature [C], T is the local temperature, HGR is the metabolic heat generated rate and S L P LRT   is the external power losses due to magnetic nanoparticles. The convection boundary condition (Robins type) was used:
κ T n = h ( T ext T )
where h is the heat transfer coefficient, and κ is the thermal conductivity of the skin layer.

4. Conclusions

The modified Pennes heat equation with Robin’s boundary condition to compute the temperature distribution within the realistic breast phantom was investigated. The single-domain magnetic and eddy current mediated power loss densities were applied as the parallel external heat sources. It was shown that due to the anatomical structure of the breast tissues, the “hot spots” could occur together with high-temperature gradients. That is why this phenomenon should be considered during treatment planning.
More importantly, the distribution of the magnetic power losses along with the temperature was not homogeneous in the cancer area. It was also demonstrated that the viscosity played a crucial role, and the magnetic fluid concentration was critical for effective MFH treatment. The behaviour of MNPs within the tumour volume was also analysed as their lack of mobilisation leads to lower power dissipation. Consequently, it is vital to know if the MNPs are free to move in the interstitial fluid or whether they are immobilised in the tumour tissue. Thus, magnetic nanoparticles–tissue interactions hold the equal potential to affect treatment.
To conclude, from the MFH-treatment planning point of view, it is crucial to prevent healthy surrounding tissues from overheating, so particular attention should be given to the temperature distribution in the tumour vicinity. The accuracy of magnetic fluid thermal therapy determines its role in prospective treatment planning for the individual patient because inadequate thermal doses to the tumour can cause a failure of treatment. Hence, the proposed numerical model will help to calibrate magnetic fluid temperature measurements when exposed to radiofrequency, and realistic models like the one presented here will enable in silico evaluation of MFH efficiency prior to clinical treatment.

Author Contributions

Conceptualisation and methodology were provided by A.M. Simulations were conducted by A.M. Data analysis and their visualisation, writing and editing were prepared by M.S.

Funding

The research was funded by the University of Life Sciences Lublin [AM], nanoTherics Ltd. [MS]. Currently, [MS] is funded by the Engineering and Physical Sciences Research Council (EPSRC, UK) and the Imperial College London [EP/N509486/1: 1979819].

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Subramanian, M.; Pearce, G.; Guldu, O.K.; Tekin, V.; Miaskowski, A.; Aras, O.; Unak, P. A pilot study into the use of FDG-mNP as an alternative approach in neuroblastoma cell hyperthermia. IEEE Trans. Nanobiosci. 2016, 15, 517–525. [Google Scholar] [CrossRef] [PubMed]
  2. Subramanian, M.; Miaskowski, A.; Pearce, G.; Dobson, J. A coil system for real-time magnetic fluid hyperthermia microscopy studies. Int. J. Hyperth. 2016, 32, 112–120. [Google Scholar] [CrossRef] [PubMed]
  3. Miaskowski, A.; Sawicki, B.; Subramanian, M. Identification of diffusion coefficients in heat equation on the base of non-adiabatic measurements of ferrofluids. In Proceedings of the 17th International Conference in Computational Problems of Electrical Engineering (CPEE), Sandomierz, Poland, 14 September 2016; IEEE: Piscataway, NJ, USA, 2016; pp. 1–4. [Google Scholar]
  4. Subramanian, M.; Miaskowski, A.; Mahapatro, A.K.; Dobson, J. Practical bioinstrumentation developments for AC magnetic field-mediated magnetic nanoparticle heating applications. Appl. Phys. A 2019, 125, 194. [Google Scholar] [CrossRef] [Green Version]
  5. Miaskowski, M.; Balakrishnan, P.; Subramanian, M.; Hovorka, O. An in vivo coil setup for AC magnetic field-mediated magnetic nanoparticle heating experiments. In Proceedings of the 2019 41st Annual International Conference of the IEEE Engineering in Medicine and Biology Society, Berlin, Germany, 2019. in press. [Google Scholar]
  6. Creixell, M.; Bohorquez, A.C.; Torres-Lugo, M.; Rinaldi, C. EGFR-targeted magnetic nanoparticle heaters kill cancer cells without a perceptible temperature rise. ACS Nano 2011, 5, 7124–7129. [Google Scholar] [CrossRef] [PubMed]
  7. Ivkov, R. Magnetic nanoparticle hyperthermia: A new frontier in biomedicine? Int. J. Hyperth. 2013, 29, 703–705. [Google Scholar] [CrossRef] [PubMed]
  8. Jordan, A.; Scholz, R.; Wust, P.; Fähling, H.; Felix, R. Magnetic fluid hyperthermia (MFH): Cancer treatment with AC magnetic field induced excitation of biocompatible superparamagnetic nanoparticles. J. Magn. Magn. Mater. 1999, 201, 413–419. [Google Scholar] [CrossRef]
  9. Kozissnik, B.; Bohorquez, A.C.; Dobson, J.; Rinaldi, C. Magnetic fluid hyperthermia: Advances, challenges, and opportunity. Int. J. Hyperth. 2013, 29, 706–714. [Google Scholar] [CrossRef] [PubMed]
  10. Hilger, I.; Hergt, R.; Kaiser, W.A. Use of magnetic nanoparticle heating in the treatment of breast cancer. In IEE Proceedings-Nanobiotechnology; IET Digital Library: Stevenage, UK, 2005. [Google Scholar]
  11. Hilger, I.; Hergt, R.; Kaiser, W.A. Towards breast cancer treatment by magnetic heating. J. Magn. Magn. Mater. 2005, 293, 314–319. [Google Scholar] [CrossRef]
  12. Wang, G.; Xu, D.; Chai, Q.; Tan, X.; Zhang, Y.; Gu, N.; Tang, J. Magnetic fluid hyperthermia inhibits the growth of breast carcinoma and downregulates vascular endothelial growth factor expression. Oncol. Lett. 2014, 7, 1370–1374. [Google Scholar] [CrossRef] [PubMed]
  13. Kossatz, S.; Grandke, J.; Couleaud, P.; Latorre, A.; Aires, A.; Crosbie-Staunton, K.; Ludwig, R.; Dähring, H.; Ettelt, V.; Lazaro-Carrillo, A.; et al. Efficient treatment of breast cancer xenografts with multifunctionalized iron oxide nanoparticles combining magnetic hyperthermia and anti-cancer drug delivery. Breast Cancer Res. 2015, 17, 66. [Google Scholar] [CrossRef] [PubMed]
  14. Petryk, A.A.; Giustini, A.J.; Gottesman, R.E.; Kaufman, P.A.; Hoopes, P.J. Magnetic nanoparticle hyperthermia enhancement of cisplatin chemotherapy cancer treatment. Int. J. Hyperth. 2013, 29, 845–851. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  15. Kan-Dapaah, K.; Rahbar, N.; Soboyejo, W. Implantable magnetic nanocomposites for the localized treatment of breast cancer. J. Appl. Phys. 2014, 116. [Google Scholar] [CrossRef]
  16. Reynoso, F.J.; Lee, C.D.; Cheong, S.K.; Cho, S.H. Implementation of a multisource model for gold nanoparticle-mediated plasmonic heating with near-infrared laser by the finite element method. Med. Phys. 2013, 40, 073301. [Google Scholar] [CrossRef] [PubMed]
  17. Miaskowski, A.; Sawicki, B. Magnetic fluid hyperthermia modeling based on phantom measurements and realistic breast model. IEEE Trans. Biomed. Eng. 2013, 60, 1806–1813. [Google Scholar] [CrossRef] [PubMed]
  18. Burfeindt, M.J.; Colgan, T.J.; Mays, R.O.; Shea, J.D.; Behdad, N.; Van Veen, B.D.; Hagness, S.C. MRI-derived 3-D-printed breast phantom for microwave breast imaging validation. IEEE Antennas Wirel. Propag. Lett. 2012, 11, 1610–1613. [Google Scholar] [CrossRef] [PubMed]
  19. Miaskowski, A.; Subramanian, M. Specific absorption rate parameter model in magnetic hyperthermia. In Proceedings of the 2017 18th International Conference on Computational Problems of Electrical Engineering (CPEE), Kutna Hora, Czech Republic, 11–13 September 2017; IEEE: Piscataway, NJ, USA, 2017; pp. 1–4. [Google Scholar]
  20. Miaskowski, A.; Sawicki, B.; Subramanian, M. Single-domain nanoparticle magnetic power losses calibrated with calorimetric measurements. Bulletin of the Polish Academy of Sciences. Tech. Sci. 2018, 66, 509–516. [Google Scholar]
  21. Hasgall, P.A.; Di Gennaro, F.; Baumgartner, C.; Neufeld, E.; Gosselin, M.C.; Payne, D.; Klingenböck, A.; Kuster, N. IT’IS Database for Thermal and Electromagnetic Parameters of Biological Tissues. Version 2.6, 13 January 2015. Available online: www.itis.ethz.ch/database (accessed on 10 May 2016).
  22. Miaskowski, A. Magnetic Fluid Hyperthermia Treatment Planning Correlated with Calorimetic Measurements under Non-Adiabatic Conditions; Science Publishing Society Libropolis: Lublin, Poland, 2018; ISBN 978-83-63761-34-9. [Google Scholar]
  23. Mamiya, H. Recent Advances in Understanding Magnetic Nanoparticles in AC Magnetic Fields and Optimal Design for Targeted Hyperthermia. J. Nanomater. 2013. [Google Scholar] [CrossRef]
  24. Pennes, H.H. Analysis of tissue and arterial blood temperatures in the resting human forearm. J. Appl. Physiol. 1948, 1, 93–122. [Google Scholar] [CrossRef] [PubMed]
  25. Rossmanna, C.; Haemmerich, D. Review of temperature dependence of thermal properties, dielectric properties, and perfusion of biological tissues at hyperthermic and ablation temperatures. Crit. Rev. Biomed. Eng. 2014, 42, 467–492. [Google Scholar] [CrossRef] [PubMed]
  26. Sawicki, B.; Miaskowski, A. Nonlinear higher-order transient solver for magnetic fluid hyperthermia. J. Comput. Appl. Math. 2014, 270, 143–151. [Google Scholar] [CrossRef]
Figure 1. Breast phantom model—breast gland (blue), beast fat (yellow), fat (green), skin (red), muscle (orange) and tumour (grey).
Figure 1. Breast phantom model—breast gland (blue), beast fat (yellow), fat (green), skin (red), muscle (orange) and tumour (grey).
Ijms 20 04644 g001
Figure 2. Distribution of the power density due to eddy currents in the breast model without the tumour (left) and with the tumour emended in the “hot spot” (right). The arrow (white) indicates the “hot spot”.
Figure 2. Distribution of the power density due to eddy currents in the breast model without the tumour (left) and with the tumour emended in the “hot spot” (right). The arrow (white) indicates the “hot spot”.
Ijms 20 04644 g002
Figure 3. Temperature distribution (cross-section) in the breast model after 1800 s due to the eddy current effect together with the tumour position (xz- and zy-cross-section).
Figure 3. Temperature distribution (cross-section) in the breast model after 1800 s due to the eddy current effect together with the tumour position (xz- and zy-cross-section).
Ijms 20 04644 g003
Figure 4. Single domain magnetic power losses—magnetic nanoparticles (MNPs) mobilised vs MNPs immobilised in the tumour.
Figure 4. Single domain magnetic power losses—magnetic nanoparticles (MNPs) mobilised vs MNPs immobilised in the tumour.
Ijms 20 04644 g004
Figure 5. XY-cross section temperature distribution passing through the middle of the tumour after 1800 s for different concentrations: ϕ r (A)   ϕ r = 1.0 , (B)   ϕ r = 0.5 and (C)   ϕ r = 0.35 . The mobilised MNPs are on the left, and the immobilised MNPs is on the right. The black colour indicates the 40 °C area.
Figure 5. XY-cross section temperature distribution passing through the middle of the tumour after 1800 s for different concentrations: ϕ r (A)   ϕ r = 1.0 , (B)   ϕ r = 0.5 and (C)   ϕ r = 0.35 . The mobilised MNPs are on the left, and the immobilised MNPs is on the right. The black colour indicates the 40 °C area.
Ijms 20 04644 g005aIjms 20 04644 g005b
Figure 6. Temperatures over time in the middle of the tumour for different magnetic fluid concentrations (mobilised case).
Figure 6. Temperatures over time in the middle of the tumour for different magnetic fluid concentrations (mobilised case).
Ijms 20 04644 g006
Figure 7. The temperature difference Δ T between mobilised and immobilised cases over ferrofluid relative concentration in the tumour together with linear approximation.
Figure 7. The temperature difference Δ T between mobilised and immobilised cases over ferrofluid relative concentration in the tumour together with linear approximation.
Ijms 20 04644 g007
Figure 8. Isosurfaces 40   and 43   for concentration: ϕ r = 0.35 —mobilised. The tumour is green.
Figure 8. Isosurfaces 40   and 43   for concentration: ϕ r = 0.35 —mobilised. The tumour is green.
Ijms 20 04644 g008
Figure 9. A nonlinear model of temperature-dependent blood perfusion for tumour.
Figure 9. A nonlinear model of temperature-dependent blood perfusion for tumour.
Ijms 20 04644 g009
Table 1. Physical parameters of breast tissues for f = 171 kHz [21]. Muscle tissue parameters were set for the tumour.
Table 1. Physical parameters of breast tissues for f = 171 kHz [21]. Muscle tissue parameters were set for the tumour.
Tissue ρ σ c κ HGR*HGR**
[kg/m3][S/m][J/kg/C][W/m/C][W/kg][W/m3/C]
Skin1109.00.173390.50.37221.6477969.16
Muscle (tumour)1090.40.3553421.20.49490.90612705.95
Breast gland1040.50.54122960.00.33452.32310542.6
Fat911.00.0572348.30.211450.52012.82
Breast fat911.00.02222348.330.2090.72782892.22
HGR*—Heat Generation Rate, HTR**—Heat Transfer Rate.
Table 2. Data used to calculate single-domain magnetic power losses [20].
Table 2. Data used to calculate single-domain magnetic power losses [20].
HmaxKkBTηδLog-normal
kA/mkJ/m3J/KKkg/m/snmg(µ, s, VM)
92.030.01.38−23298.08.94−42.0g(15.2, 0.19, VM)
Table 3. Single-domain (S-D) magnetic power losses with regard to the relative concentration ( ϕ r ) of ferrofluid in the tumour [22].
Table 3. Single-domain (S-D) magnetic power losses with regard to the relative concentration ( ϕ r ) of ferrofluid in the tumour [22].
ϕ r S-D Magnetic Power
(Mobilised)
S-D Magnetic Power
(Immobilised)
kW/m3kW/m3
1.0243.0170.0
0.71174.0121.0
0.5121.084.2
0.3584.259.3
0.372.851.2

Share and Cite

MDPI and ACS Style

Miaskowski, A.; Subramanian, M. Numerical Model for Magnetic Fluid Hyperthermia in a Realistic Breast Phantom: Calorimetric Calibration and Treatment Planning. Int. J. Mol. Sci. 2019, 20, 4644. https://doi.org/10.3390/ijms20184644

AMA Style

Miaskowski A, Subramanian M. Numerical Model for Magnetic Fluid Hyperthermia in a Realistic Breast Phantom: Calorimetric Calibration and Treatment Planning. International Journal of Molecular Sciences. 2019; 20(18):4644. https://doi.org/10.3390/ijms20184644

Chicago/Turabian Style

Miaskowski, Arkadiusz, and Mahendran Subramanian. 2019. "Numerical Model for Magnetic Fluid Hyperthermia in a Realistic Breast Phantom: Calorimetric Calibration and Treatment Planning" International Journal of Molecular Sciences 20, no. 18: 4644. https://doi.org/10.3390/ijms20184644

APA Style

Miaskowski, A., & Subramanian, M. (2019). Numerical Model for Magnetic Fluid Hyperthermia in a Realistic Breast Phantom: Calorimetric Calibration and Treatment Planning. International Journal of Molecular Sciences, 20(18), 4644. https://doi.org/10.3390/ijms20184644

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop