A Computational Study on Magnetic Nanoparticles Hyperthermia of Ellipsoidal Tumors

: The modelling of magnetic hyperthermia using nanoparticles of ellipsoid tumor shapes has not been studied adequately. To ﬁll this gap, a computational study has been carried out to determine two key treatment parameters: the therapeutic temperature distribution and the extent of thermal damage. Prolate and oblate spheroidal tumors, of various aspect ratios, surrounded by a large healthy tissue region are assumed. Tissue temperatures are determined from the solution of Pennes’ bio-heat transfer equation. The mortality of the tissues is determined by the Arrhenius kinetic model. The computational model is successfully veriﬁed against a closed-form solution for a perfectly spherical tumor. The therapeutic temperature and the thermal damage in the tumor center decrease as the aspect ratio increases and it is insensitive to whether tumors of the same aspect ratio are oblate or prolate spheroids. The necrotic tumor area is affected by the tumor prolateness and oblateness. Good comparison is obtained of the present model with three sets of experimental measurements taken from the literature, for animal tumors exhibiting ellipsoid-like geometry. The computational model enables the determination of the therapeutic temperature and tissue thermal damage for magnetic hyperthermia of ellipsoidal tumors. It can be easily reproduced for various treatment scenarios and may be useful for an effective treatment planning of ellipsoidal tumor geometries.


Introduction
The most well-known techniques for cancer treatment are chemotherapy, radiation, immunotherapy and targeted therapy [1][2][3]. These treatments suffer from various severe side effects such as hair loss, nausea, infertility, nerve damage and general discomfort to the patient [3].
The hyperthermia of tumors is regarded as a viable treatment provided that the tumor locations are known and have not metastasized or when it is difficult to surgically remove the tumor, such as in liver cancer [16]. The ultimate treatment goal is to raise the tumor temperature above the physiological temperature for a certain time and to damage the malignant cells through several mechanisms such as protein denaturation, apoptosis and necrosis [17], without significantly affecting the surrounding healthy tissue. In mild hyperthermia the tumor temperature usually ranges between 40-45 • C and typically above 50 • C for ablation therapy [18][19][20][21][22].
MNPs have a broad spectrum of applications such as heat exchangers, solar collectors and separation of heavy metals from water [23][24][25]. Hyperthermia using MNPs (magnetic hyperthermia) has received increased attention due to its ability to deliver adequate heating power to several types of tumors [26]. The therapy involves the administration sarcoma tumor both experimentally and computationally. The shape of the tumor was also approximated by an ellipsoidal shape. Kandala et al. [94] proposed a computational model for the utilization of power modulation for magnetic nanoparticle hyperthermia of elliptic (2D) and ellipsoidal (3D) tumors. In the above-mentioned studies, the aspect ratio of the ellipsoid tumors was fixed.
Egolf et al. [95] developed an analytical model for the transient temperature evolution in three tumor shapes of equal volume: a perfect spherical, a prolate spheroid with an aspect ratio of approximately three and an oblate spheroid with an aspect ratio of eight. Spatial temperature distributions in the tumor and the surrounding healthy tissue were neglected. Their results show that the uniform temperature in the spherical tumor was higher than in the prolate spheroid tumor and much higher than the oblate spheroid tumor. Tehrani et al. [96] studied numerically oblate and prolate spheroid tumors of equal volumes for the treatment of microwave ablation using a coaxial antenna. In their work the aspect ratio of the ellipsoids varied from one to five. Their results show that the aspect ratio has a significant effect on the extent of the ablation zone in the tumor.
The objective of the present investigation is to provide a systematic study for magnetic nanoparticles hyperthermia of ellipsoidal tumors of various aspect ratios and to compare the results of the numerical predictions to experimental data. The tumors are modeled as prolate and oblate spheroids of equal volumes.

Geometrical Description
The general equation of an ellipsoid is given by [97]: where a, b and c are the lengths of the principal semi-axes. For the case of all lengths equal a = s = c = R, Equation (1) describes a perfect sphere with radius R. In the present work we are interested for ellipsoids with a = c (symmetric around the y axis), while perfect spherical tumors constitute only a limit-case scenario. Such shapes are usually referred to as ellipsoids by revolution. Here, the y-axis is set as the axis of revolution. Two basic cases may be distinguished: (i) oblate spheroids with semi-axis a > b (ii) prolate spheroids with semi-axis a < b as shown in Figure 1.
In addition, we define the aspect ratio AR for the generated ellipsoids using the following notation [96]: AR = major axis length minor axis length (2) Increasing AR leads to ellipsoidal tumors with more elongated shapes. The surface S of the ellipsoids is expressed via the following formulation [98]: where e is the eccentricity of the ellipsoid. The volume of the ellipsoids is given by [95]: All the generated ellipsoidal tumors are set to have equal volumes.
The dimensions of the ellipsoid tumors used in this work are shown in Table 1. The tumor geometries are taken to have the same volume, as calculated from Equation (4). The range of the selected tumor dimensions are within the range of earlier works [80,86,95,96]. It is also assumed that the ellipsoidal tumors are surrounded by healthy tissue of spherical geometry, as shown in Figure 2. The region of the healthy tissue is assumed to be significantly larger than the tumor. In particular, the radius of the healthy tissue R h is taken approximately eight times larger than the semi-major axis of the tumor with the highest aspect ratio. Due to the rotational symmetry of the geometries, the present thermal problem can be solved as an axisymmetric problem instead of a 3D one, which substantially decreases the computational cost of the numerical simulations [99].  For the discretization of the computational domains, we used a combination of regular and unstructured meshes consisting of triangular cells. All meshes were constructed using GMSH software [100]. The unstructured mesh is used to discretize the tumor region as well as a healthy tissue layer around the tumor. We followed this approach to better capture the surface geometry of the tumors with high aspect ratios (e.g., AR = 8). Two sample meshes for AR = 2 are shown in Figure 3.  . Two representative computational meshes used in the study focused at the tumor region and the close area around it. Magnified views close to the tumor/healthy tissue boundary are also shown. Both meshes correspond to tumors with aspect ratio AR = 2.

Bio-Heat Transfer Analysis
Bio-heat transfer between the ellipsoidal tumor and the surrounding healthy tissue is expressed by the thermal energy balance for perfused tissues described by the Pennes bio-heat equation [93,94]: where the subscript n stands for the tissue under consideration (n = 1 for tumor and n = 2 for healthy tissue) and the subscript b corresponds to blood properties. Also, ρ n and ρ b denote the densities of the tissues and the blood respectively, c n and c b are the corresponding heat capacities, T(x,y,t) is the local tissue temperature, k n is the tissue thermal conductivity, w b is the blood perfusion rate, and T b = 37 • C is the blood temperature. The left-hand side term in Equation (5) expresses the time rate of change of internal energy per unit volume. The first term on the right-hand side of Equation (5) represents the heat conduction in the tissue. The second term represents an additional change in the internal energy per unit volume associated with blood perfusion in tissue, assuming that the rate of heat transfer between tissue and blood is proportional to the blood perfusion rate and the difference between the local tissue temperature and the blood temperature, as suggested in [65]. In addition, Q met,n is the internal heat generation rate per unit volume associated with the metabolic heat production. Finally, Q s is the power dissipation density by the MNPs. It is assumed no leakage of MNPs to the surrounding healthy tissue. Therefore, Q s is only applied to the cancerous region filled with the MNPs, similarly to several earlier investigations, e.g., [15,67,69,92]. Following the assumptions by [42,44,92,95,101], a uniform distribution of nanoparticles is assumed. The temperature of the surface of the healthy tissue region (see Figure 2), which is far away the tumor surface as explained in the previous section, is set to T s = T b [69,75]. For the initial temperature condition we set T(x,y,t = 0) = T b . A symmetry (zero heat flux) boundary condition is applied on the xand y-axes (line OA and line OB respectively in Figure 2). Also, heat flux and temperature continuity are considered on the tumor-healthy tissue interface as follows [12,93]: where n is the direction vector perpendicular to the tumor surface. To carry out numerical simulations for the different AR cases, typical tissue thermophysical values (for example prostate tissue) are chosen from earlier works [36,63,67,71,94] and shown in Table 2. Note that we have assumed that the perfusion rate is independent from region and temperature [31,58]. We have also assumed that the thermophysical properties of the tumor and the healthy tissue are the same, based on earlier works [36,67,71,94] For a complete mathematical formulation of heat transfer in a solid body and the implementation of "internal heat generation" terms within a numerical setting, the interested reader can be referred to [102].  [99,103] which is based on the finite volume method (FVM). The equations are solved utilizing a preconditioned bi-conjugate gradient (PBiCG) with diagonal incomplete-LU (DILU) preconditioning [104]. The tolerance was set to 10 −8 . To speed up the simulations we ran them in parallel mode on a system of distributed processors by geometric domain and associated field decomposition, as described in [105].

Heat Generation by the Magnetic Nanoparticles
As explained in the work by Miaskowski and Sawicki [49], the superparamagnetic heating phenomenon is the base of magnetic nanoparticles hyperthermia. Rosensweig [33] formulated a theoretical model for the heating of colloidal magnetic fluid (ferrofluid) due to application of an AFM. In his work, the power dissipation density (in W/m 3 ) from a monodispersed magnetic fluid filled with superparamagnetic particles subjected to an alternating magnetic field [15] is given by: where µ 0 = 4π × 10 −7 H·m −1 is the vacuum magnetic permeability, f is the frequency of the magnetic field with magnetic field intensity H 0 (A·m −1 ), χ 0 corresponds to the magnetic susceptibility of the magnetic fluid defined by where χ i the initial susceptibility given by and with M d the domain magnetization of a suspended particle, V m = 4/3R 3 is the magnetic volume for a nanoparticle of radius R, k B = 1.381 × 10 −23 J·K −1 is the Boltzmann constant and T is the absolute temperature. In Equation (10), ϕ is the solid volume fraction of the nanoparticles. The power dissipation is accomplished via two physical processes, the Néel and Brownian relaxation. Because these processes take place in parallel, the effective relaxation time in Equation (5) is given by where τ B the Brownian characteristic relaxation time with η the viscosity of the matrix fluid and V H is taken as the hydrodynamic volume of the nanoparticle related to V m as V H = (1 + δ/R) 3 V m where δ is the thickness of a sorbed surfactant layer (δ = 2 nm according to Rosensweig [33]). In Equation (13) τ N is the Néel characteristic relaxation time given by [15]: where τ 0 ≈ 10 −9 s is an attempt time [15,49] and K is the anisotropy constant (J/m).

Tissue Thermal Damage
In the present work, the extent of the tissue thermal damage is determined with the Arrhenius kinetic model, which has been used in several studies, e.g., [21,76,106]. This model was originally proposed by Henriques and Moritz [107,108], where the tissue damage is expressed via a dimensionless damage parameter Ω, given by: where τ is treatment duration, C(0) is the original concentration of the tissue constituent, C(τ) the undamaged tissue constituent at the end of treatment heating, A the frequency factor (s −1 ), E a the activation energy (J·mol −1 ) and R the gas constant. The temperature T(x,y,t) in Equation (15) is in Kelvin. Ω = 1 means that the damage process is 63.2% complete [21,54] and the tissue may be assumed to be irreversibly damaged [54,106]. The values of the frequency factor and activation energy depend upon the cell line. For the computational results of the present investigation, the constituent cells of the tissue are assumed to be the AT1 subline of Dunning R3327 rat prostate cells with the corresponding values obtained from earlier works [76,92], namely: A = 2.99 × 10 37 s −1 and E a = 244.8 kJ·mol −1 .

Mesh and Timestep Sensitivity Analysis
A mesh sensitivity analysis was carried out to determine the size of the mesh. The computational sample meshes are shown in Table 3. The mesh sensitivity was performed on an oblate spheroidal tumor with AR = 8. The quantity for which the analysis was performed is the tumor temperature at a distance 2 mm above the tumor geometric center that lies on the y-axis (see Figure 2) after 30 min of treatment. The simulation results in Table 3 show that increasing the mesh size and the temperature on the above-mentioned location generally increases. However, a closer look at the values shows that from mesh 3 to mesh 4 the temperature values change only on the third decimal, which means that temperature change between these two meshes is about 0.01%. Since this change is very small, mesh 3 is selected for the numerical simulations. Furthermore, the timestep in the present work is set to 1 s. Simulation runs with a smaller time step were also performed, namely 0.1 s, which resulted in no significant difference (<0.001%) in the solution. Moreover, the treatment temperature behavior of the computational model is verified with the closed-form transient solution proposed by Liangruksa et al. [67] for a tumor with AR = 1 (perfect sphere). In their work the solution is given in a dimensionless form (Equations (16) and (17) in [67]). Our computational results are in excellent agreement with the closed-form solution, as shown in Figure 4.  [67]. T* is the dimensionless temperature, r* the dimensionless distance from the tumor center and t* is the dimensionless time, as defined in Liangruksa et al. [67].

Computational Results and Discussion
Magnetite (Fe 3 O 4 ) nanoparticles are selected as heat mediators assuming typical magnetic properties, as shown in Table 4. The magnetic field properties are also presented in Table 4. Note that for the selected H 0 and f values we find H 0 f = 1.496 × 10 9 A·m −1 ·s −1 , which falls between the limits of Atkinson-Brezovich (4.85 × 10 8 A·m −1 ·s −1 ) and Dutz-Hergt (5 × 10 9 A·m −1 ·s −1 ) criterions [29,30]. Also, the nanoparticles volume fraction we used is low. Therefore, the effective tumor parameters of MNPs-saturated tissue are nearly identical to tumor parameters without nanoparticles that are used in the model. By substituting these parameters in Equation (8) we find Q s = 1.91 × 10 5 W/m 3 , which is within the range of earlier publications [63,65,68]. Table 4. Magnetic nanoparticles and magnetic field parameters [33,36,47,49,63].

Parameter
Value The computational results are carried out for a 30 min treatment, since in magnetic hyperthermia it is desirable to have a treatment duration as short as possible for safety purposes [76,109]. The duration the AMF is switched on and heats the nanoparticles is assumed to be 22 min [76]. After that time and for the remaining eight minutes of the treatment, the magnetic field is off and stops heating.
To gain an initial understanding of the tissue temperature distribution, in Figure 5 the temperature field is presented for a region near oblate tumors after 22 min of treatment at various AR values. Note that all tumor shapes have the same volume. In all cases the maximum treatment temperature is observed at the tumor center. As the aspect ratio AR increases, this maximum temperature decreases. This is also the case for the temperature on other regions inside the tumor. A similar behavior is observed for the prolate spheroidal tumors as shown in Figure 6. This behavior is consistent with the results of our earlier preliminary work [99]. Note that the tissue and nanoparticle parameters used in [99] are different than the ones used in the present work. Moreover, the bio-heat equation in [99] was solved under a steady state condition. In the current investigation we have used the more realistic temperature time dependent approach which further allows us to determine the extent of the tissue thermal damage with the Arrhenius thermal damage model.   Figure 7 shows time-dependent temperature profiles at the tumor center for all the considered shapes. The AR value appears to have a significant effect on the tumor temperature evolution. At very initial times the temperature in the center increases rapidly and it is relatively independent of the aspect ratio and whether the tumor is an oblate or prolate spheroid. However, at intermediate times, the temperature becomes significantly dependent on the aspect ratio and it remains nearly insensitive to prolate or oblate spheroid shapes of the same AR. For instance, for AR = 1 after 22 min of treatment, the temperature is 46.62 • C, while for AR = 8 the maximum temperature is approximately 3.5 • C less for both prolate and oblate spheroidal tumors. The above treatment temperature trends could be explained in terms of the surface area of the tumors. From Equation (3), increasing the AR-value of an ellipsoid tumor leads to an increase of the surface area. For instance, for an oblate tumor with AR = 8 the surface area is approximately twice the surface area of a perfectly spherical tumor of the same volume and for a prolate tumor of AR = 8 it is approximately 1.5 times. Since the available surface area increases, the heat transfer from the nanoparticle-saturated tumor to the cooler healthy tissue increases, leading to a drop of the tumor temperature. Figure 8 shows the temperature values at two representative tumor surface probe locations as shown in Figure 8a. Observation point P1 is fixed on the tumor surface at the semi-minor axis and P2 fixed on the tumor surface at the semi-major axis. For a tumor with AR = 1, both probes give exactly the same temperature, as shown in Figure 8b. Of course, this is due to the perfect spherical symmetry for this tumor and has been shown by several other works. However, as shown in Figure 8c for a prolate and an oblate tumor with AR = 2, the treatment temperature as a function of time, at P1 and P2, differs significantly. While at the start of the treatment the temperature rises rapidly in a similar manner at both probe locations, after a relatively short time the temperature history at P1 is substantially higher than the temperature at P2. For the oblate tumor, after 22 min of treatment, the temperature at P1 is approximately 2 • C higher than the temperature at P2. A similar temperature difference is observed at the corresponding locations for a prolate spheroidal tumor of the same aspect ratio, but with lower treatment temperature values than the oblate. It should be also noted that for the oblate spheroid the temperature at P1 is higher than the temperature at the corresponding location of a perfectly spherical tumor (AR = 1), as compared to Figure 8b. For tumors with AR = 4 (see Figure 8d) the temperature difference between the two probe locations, after 22 min of treatment, increases to approximately 3 • C (for both oblate and prolate shapes), which is also in the case of tumors with an AR of 8 (see Figure 8e). As expected, due to the larger available surface for heat transfer, the temperatures for the case of AR = 8 are lower than for tumors with an AR of 4. Once the treatment temperature of the proposed model is determined, the degree of thermal damage in the tumor and healthy tissue can be calculated using the Arrhenius kinetic model [21]. The maximum treatment temperature profiles as a function of treatment time presented in the previous figures can be used to evaluate the degree of thermal damage in the tumor center.
The thermal damage at the tumor center as a function of treatment time is shown in Figure 9. At the beginning of the treatment, the thermal damage of the malignant tissue at the tumor center is low and relatively unaffected by the tumor aspect ratio and whether it is a prolate or oblate spheroid. As the treatment time progresses, the highest thermal damage is observed in the tumors with AR = 1 and the lowest for tumors with AR = 8. For example, at 20 min of treatment the Arrhenius damage index is Ω ≈ 2.75 for AR = 1 and Ω ≈ 1 for AR = 8. The thermal damage evolution in the tumor center is generally independent of the tumor prolateness or oblateness, but it is strongly affected by the aspect ratio. Also, after approximately 22 min of procedure, Ω appears to be relatively stable for each aspect ratio. If we select the frequently quoted Ω = 1 threshold value above, which has a relatively high chance of causing irreversible damage the malignant tissue, we find that for AR = 1 this value is reached after approximately 10 min. of treatment. However, for AR = 8, more than 20 min of treatment are required to reach the Ω = 1 threshold value, and shorter times are required for tumors with smaller AR values.  Figure 10 shows the Arrhenius integral Ω in the tumor for three representative geometries. Figure 10a is for AR = 1, Figure 10b for a prolate tumor with AR = 2 and in Figure 10c for an oblate one with AR = 2. For each case, the thermal damage in the central tumor region is higher than regions near the tumor surface, which is the direct result of the temperature distribution (see Figures 5 and 6). To determine the extent of the damaged region in the ellipsoidal tumors, the frequently quoted Ω = 1 criterion is used as carried out by Eltejaei et al. [106] and Andreozzi et al. [110]. In these studies, the boundary between the necrotic tumor area (Ω > 1) and the undamaged tumor region (Ω < 1) was determined using the Ω = 1 iso-contour. If we follow this approach in the sample tumor cases shown in Figure 10, we find that the location of this boundary is significantly affected by the aspect ratio and by whether the tumor is a prolate or oblate spheroid. Interestingly, for an oblate tumor of AR = 2, a small portion of the surrounding healthy tissue may necrose, which is not the case for a prolate tumor of the same aspect ratio.

Comparison with Experiments
In the present section, the computational model is compared with three sets of experimental measurements by different researchers.
The first comparison is with the experimental measurements by Hamaguchi et al. [86]. In their work, the tumor was a squamous cell carcinoma on the cervical lymph node of female Japanese white rabbits. According to [86], the ratio of the volume of the cancerous lymph node to that of normal lymph nodes was approximately 11 (cancerous lymph node volume 1161 ± 276.4 mm 3 , normal lymph node volume 105.6 ± 43.37 mm 3 ). Twenty days after tumor transplantation, magnetite nanoparticles, of average core magnetite size D = 10 nm, were injected from the rabbit tongue. The average nanoparticle uptake from the cancerous lymph was approximately 4 mg ± 1 mg. For the hyperthermia treatment, a transistor inverter was used with frequency 118 kHz and magnetic field strength of 384 Oe or 30.6 kA/m. From a histological section of the swollen lymph, we approximated the tumor shape with a prolate spheroid that we fitted on top of the tumor. Two tumor-shaped approximations are considered as shown in Figure 11a,b. In Case A we find AR ≈ 1.8, and for case B, AR ≈ 2.2. Inserting the tumor volume value in Equation (4) we calculate a ≈ 5.1 mm and from Equation (2) we find b ≈ 9.18 mm for Case A. In Case B we find a ≈ 4.78 mm and b ≈ 10.44 mm. From the values reported by Hamaguchi et al. [86] and using Rosensweig's theory (Equations (8)- (14)) we find the heat dissipated by the nanoparticles equal to 2.1 × 10 5 W/m 3 . For the blood perfusion we use 1.3 × 10 −3 s −1 within the range of earlier works [63,[92][93][94]. The treatment temperature simulation results, for Case A and Case B, are shown in Figure 11c,d, respectively. For the 4 mg dosage, the predictions are in qualitative agreement with the temperature measurements by Hamaguchi et al. [86]. Some small differences are observed between the numerical result of Case A and Case B, with Case A being slightly closer to the measurements. It should be pointed it out that Hamaguchi et al. [86] report that the 4 mg nanoparticle uptake from the cancerous lymph has approximately ±1 mg uncertainty in the measurement. Interestingly, if we use a 5 mg dosage for Case A and Case B our results are in better agreement with the experimental temperature measurements by Hamaguchi et al. [86]. Figure 11. Two cases approximating the tumor shape from a histological cross-section by Hamaguchi et al. [86], with a prolate spheroid. Note that the tumor histological cross-section has been redrawn from the original: (a) prolate spheroid shape, case A with AR ≈ 1.8, on top of the redrawn tumor and (b) prolate spheroid shape, case B with AR ≈ 2.2, on top of the redrawn tumor. Plots (c,d) show parametric comparison of the numerically determined temperature at the tumor center with the measured temperature by [86]. Temperature data points and bars are mean values and standard deviation respectively of 5 independent experiments.
Subsequently, the computational model predictions are compared with experimental measurements and with 3D computational results by Pearce et al. [92] for murine mammary adenocarcinoma tumors. The tumor volume was 329 mm 3 and was heated for 600 s. In their work, iron oxide nanoparticles (IONP) of 10-20 nm in diameter were. The IONPs were exposed to magnetic field strengths between 20 and 50 kA/m (rms) at 162 kHz. Pearce et al. [92] report that the transient temperature was recorded at a location called "center" and another location separated by 3 mm, called "tip". They also mention that the center probe location was placed as close as possible to the approximate center of the tumor. A redrawn histologic section of the tumor in Pearce et al. [92] is shown in Figure 12. As in the previous experimental comparison, we approximated the tumor shape with a prolate spheroid that we fitted on top of the tumor. Two tumor shape approximations were considered, as shown in Figure 12a,b. For Case A we found AR ≈ 1.29 and for case B, AR ≈ 1.6. We then found a ≈ 3.9 mm and b ≈ 5.1 mm for Case A and for Case B we find a ≈ 3.6 mm and b ≈ 5.8 mm. The experimental temperature measurements close to the tumor center (probe location center) and about 3 mm from the tumor center (probe location tip), are shown in Figure 12c,f. According to Pearce et al. [92], the value of heat generated by the nanoparticles was ". . . adjusted somewhat until the experiment maximum transient temperature (or steady state) temperature record from the embedded probes was closely approximated by the numerical model result.". They also report that the same approach was followed for the blood perfusion: ". . . adjusted to improve match to the measurements. . . ". The numerical results given by [92] are shown in Figure 12 with broken lines. The adjusted by Pearce et al. [92] value for the generated heat by the nanoparticles was 1.1 × 10 6 W/m 3 . For the adjusted perfusion, according to Pearce et al. [92], the initial tumor perfusion, 3 × 10 −3 s −1 was increased to as much as 7 × 10 −3 s −1 , as required to match experimental results. If we follow the Pearce et al. [92] approach of adjusting the heat generated and the perfusion rate we find good agreement with the measurements for the probe location center, as shown in Figure 12c (Case A), using the values of 1.75 × 10 6 W/m 3 and 2.5 × 10 −3 s −1 . It should be pointed out that at t = 0 we have used the experimentally measured temperature (32 • C), while in the numerical model in [92] a higher temperature of approximately 36 • C was assumed by Pearce et al. [92], without providing an explanation for this choice. This perhaps explains the differences between our adjusted values with the ones by Pearce et al. [92]. Good agreement with the measured temperature and our model is also observed for the tip location, seen in Figure 12e, while in the prediction by Pearce et al. [92], the computational model gives higher temperatures than the experiment at this location. For the tumor geometry of Case B, we use the adjusted heat generated and blood perfusion values from Case A and compare our predictions with the experiments in Figure 12d (center location) and Figure 12f (tip location). Of course, due to the larger AR of the tumor than in Case A, the maximum temperatures are somewhat lower but reasonably close to the measurements. Unfortunately, due to the large range of two simultaneous parameters, namely, the nanoparticle diameter (10 to 20 nm) and the applied magnetic field (20 to 50 kA/m) reported in Pearce et al. [92], we could not apply Rosensweig's theory as we did for Hamaguchi et al. [86]. Subsequently, we compared the cumulative equivalent minutes at 43 • C (CEM43) of our model with the CEM43 measurements and model predictions reported by Pearce et al. [92]. According to Pearce et al. [92], the CEM43 in discrete interval form is written as where R CEM is the time scaling ratio, 43 • C is the reference temperature and ∆t i (min) is spent at temperature T i ( • C). In their work R CEM = 0.45 was chosen. Using Equation (16) for our model predictions in Figure 12 we obtain CEM43 values close to the calculated by Pearce et al. [92], as shown in Table 5.  In Figure 13 the present model is compared with experimental measurements by Ling et al. [91]. In their study, cancerous human breast cells were injected in nine mice. Hyperthermia was applied when the tumor reached a volume of 0.453 ± 0.109 cm 3 . The tumor volume in Ling et al. [91] was calculated from Equation (4). Then 0.1 mL of a mixture that was injected in the tumor tissue containing 10 wt% Fe 3 O 4 powder, 36 wt% polymethylmethacrylate (PMMA) powder and 54 wt% MMA liquid. A magnetic field of frequency 626 kHz and strength 28.6 kA/m was applied for 180 s. Ling et al. [91] report that the ferrous powder diameter was in the range of D = 20-50 nm. An ultrasound image section of the tumor from Ling et al. [91] is shown at the top of Figure 13. The red dashed line corresponds to the tumor boundary as marked by Ling et al. [91]. For the comparison, we recognize two cases. Given the tumor volume mentioned above for Case A (Figure 13a) we find AR ≈ 2.5 with a ≈ 3.5 mm and b ≈ 8.8 mm, while for Case B, AR ≈ 2.82 with a ≈ 3.37 mm and b ≈ 9.5 mm. Subsequently, we applied Rosensweig's theory to estimate the amount of heat dissipated by the nanoparticles. We initially calculated the solid volume fraction of the nanoparticles in the tumor volume. The density of the injected mixture can be calculated from the following equation where ρ Fe3O4 = 5180 kg/m 3 [33], ρ MMA = 940 kg/m 3 [111], ρ PMMA = 1180 kg/m 3 [112] and the corresponding weight fractions are mentioned above. Substituting the values in the above equation we find ρ mixture = 1450 kg/m 3 . The 0.1 mL mixture injected in the 453 mm 3 tumor volume has a mass of 0.1 mL × 1450 mg/mL = 145 mg. Since 10 wt% of the mixture is the ferrous mass, the nanoparticles weight in the tumor is 14.5 mg. This corresponds a solid volume fraction in the tumor of ϕ ≈ 6.16 × 10 −3 . Also, Ling et al. [91] measured the magnetization saturation of the nanoparticles M s = 1560 A/m (or mass magnetization saturation σ s = 0.3 emu·g −1 ). According to Rosensweig [33], the domain magnetization (M d in Equations (10) and (11)) is given by M d = M s /ϕ, where substitution in the present case gives M d ≈ 253 kA/m. If we select the value of 35 nm for the ferrous powder average size in the range 20-50 nm reported by Ling et al. [91], and use the above-mentioned values, we find (from Rosensweig's theory) 1.7 × 10 6 W/m 3 of dissipated heat. The comparison of the present model tumor surface temperature, using this value, with the measurements by Ling et al. [91] is shown in Figure 13c for Case A and Figure 13d for Case B. The blood perfusion value used for the Hamaguchi et al. [86] comparison is also used here. Note that in [91] the temperature is recorded using thermal imaging. The exact location of the tumor surface at which temperature was measured is not reported. The thermal images presented in Ling et al. [91] suggest that it is measured at a central surface location. For the comparison we have selected in our model a temperature at the tumor surface along the semi-minor axis (see Figure 13a or Figure 13b). Both Case A and B model predictions are close to the low measured temperature limit, with Case A being marginally better. If we select an average value of 30 nm for the powder size, which is also within the 20-50 nm range reported by Ling et al. [91], our model predictions are in good agreement with the measured temperatures as shown in Figure 13c,d.

Concluding Remarks
A computational study for magnetic hyperthermia using nanoparticles of ellipsoidal tumors has been presented. The tumors were approximated as equal volume prolate and oblate spheroids of various aspect ratios, surrounded by a large spherical healthy tissue region. The nanoparticles are assumed to be uniformly distributed in the whole tumor. The bio-heat transfer analysis is carried out using the Pennes bio-heat equation. The results indicate that the highest temperature is achieved in the ellipsoidal tumor center, the value of which decreases by increasing the aspect ratio of the tumor. This value appears to be insensitive to whether the ellipsoidal tumor is a prolate or oblate spheroid. Probing the temperature at the tumor surface at two locations, one along the major and one along the minor axis, reveals that oblate tumors have generally higher surface temperatures than oblate ones, the values of which strongly depend on the aspect ratio. Using the Arrhenius kinetic model for thermal damage, we find that the thermal damage in the tumor center is unaffected by whether the tumor is oblate or prolate and decreases for increasing aspect ratio. Also, the computational model produces results for the extent of the tumor necrotic area, which is affected by the aspect ratio as well as the prolateness and oblateness of the ellipsoid tumors.
The numerical model was compared with three different sets of experimental measurements involving nanoparticle hyperthermia in animal tumors which are available in the literature. In all comparisons, we have approximated each tumor shape with two prolate spheroid geometries of slightly different aspect ratios to describe as best as possible the tumor shape. Both case geometries produced results reasonably close to the measured ones. Model predictions were generally in satisfactory or perhaps good agreement with the experiments when uncertainties in the measured properties of the nanoparticles are taken into account. Also, although the parameters of the tissue used in the model are derived from different tissues (muscle [86], liver [91], prostate [92]), the comparisons show good agreement with the experimental measurements presented by other authors with the proposed numerical method.
It should be pointed out that according to Giustini et al. [113], available technologies that convey heat to tumors, such as RF, microwave, ultrasound and conductive, have not been able to target heat specifically to tumors in an effective manner, especially to metastatic ones. Hyperthermia using magnetic nanoparticles is a minimally invasive treatment that appears to offer more localized heating of the tumors as reported by many authors cited in the present work. Nonetheless, the treatment needs to be precisely assessed by taking into account several factors, such as temperature dependence over time, treatment time, tumor shape, tissue damage and nanoparticle dosage. Computational modeling is absolutely necessary for an in-depth understanding on the effect of the above parameters on the treatment outcome. The present computational model has been focused on the determination of such parameters, specifically for ellipsoidal (prolate and oblate) tumors with several aspect ratios. The model can be easily implemented and reproduced for various treatment scenarios and may be useful for further treatment planning.