key: cord-0972324-1skit2oc authors: Darvishi, Vahid; Navidbakhsh, Mahdi; Amanpour, Saeid title: Heat and mass transfer in the hyperthermia cancer treatment by magnetic nanoparticles date: 2021-11-26 journal: Heat Mass Transf DOI: 10.1007/s00231-021-03161-3 sha: cd89a241e0bf8a2c83d69b33fa402ae837589f5c doc_id: 972324 cord_uid: 1skit2oc In this study, a more precise and cost-effective method is used for studying the drug delivery and distribution of magnetic nanoparticles in fluid hyperthermia cancer treatment, and numerical methods are employed to determine the effect of blood circulation on heat transfer and estimate the success of cancer treatment. A combination of numerical, analytical, and experimental researches is being conducted, which illustrates the essential role of numerical methods in medical and biomedical science. Magnetic NanoParticles’ distribution and effects of infusion rate on the treatment are also discussed by considering the real distribution of MNPs. To increase accuracy and reduce costs in the in-vitro section, direct cutting and image processing methods are used instead of MRI. Based on the results of this section, with a tenfold increase in the infusion rate (4 μl/min to 40 μl/min), the penetration depth increases by 1 mm, which represents a nearly 17 percent increase. Concentrations of MNPs also decrease significantly at higher infusion rates. The simulations of heat transfer reveal that maximum temperatures occur at the lowest infusion rate (1.25 μl/min), and blood flow also has a significant effect on heat transfer. With an increase in the infusion rate, necrosis tissue recedes from the tumor center and approaches the border between the tumor and healthy tissue. Results also show that, in lower MNPs’ concentrations, higher infusion rates result in better treatment even though minimum infusion rates are suggested to be the best rates to facilitate distribution and treatment. Cancer is a disease that has high mortality rates, and there is no definitive treatment for it yet. Most of the cancer types appear in the form of a solid tumor. Cancerous cells can lead to cancer finally [1] . There are some prevalent treatments such as chemotherapy, radiotherapy and surgical ablation of cancerous tumors already [2] but all of them have many side effects, and they do not have enough accuracy [3, 4] . In recent years, a new treatment called hyperthermia has been investigated widely and several experiments show its high efficacy while having fewer side effects [5] [6] [7] [8] [9] . Hyperthermia is a type of thermal treatment for cancer, in which tumor cells are heated to a temperature above 42 °C and this leads to necrosis in living cells [10] . Increase in temperature also has been used for the diagnosis of the circulating cancer cells in blood [11] . There are many ways to raise the temperature of a tissue such as a laser, ultrasound and magnetic fluid. Laser irradiation as one of these methods was used successfully to make thermal damage, and the effects of some parameters such as the laser exposure time and the blood perfusion were considered [12] [13] [14] . Among all these methods, using magnetic fluid is the most effective way to raise the temperature in the body while has fewer side effects. In this method, a fluid containing magnetic nanoparticles (MNPs) is delivered into the tumor and generates adequate heating when exposed to the low alternating magnetic field. Killing tumor cells and making no side effects in healthy cells around the tumor is the primary purpose of this method. In magnetic fluid hyperthermia, the temperature can reach 42 °C to 45 °C. Raising to 42 °C is mostly used to boost other methods of treatment, whereas raising over 45 °C has been deployed as a cancer treatment method individually [15] . Heating the tumor cells leads to their necrosis or apoptosis in the magnetic fluid hyperthermia method [16, 17] . Some parameters affect the results of killing cells in this method such as nanoparticle diameter, types of coatings of MNPs, aggregation of nanoparticles in magnetic fluid, and material used in the synthesis of MNPs. Usov et al. worked on the relation between nanoparticle diameter and specific absorption rate. Results proved a substantial dependency [18] . An experimental trial indicated that cell-nanoparticle interaction reduces the cytotoxic effects of MNPs hyperthermia [19] . Dabbagh et al. suggested using a thermosensitive polymeric shell for coating the MNPs for obtaining a synergistic effect of hyperthermia. A 77% decrease in cell viability due to the synergistic effects was observed using a new coating [20] . De Mello et al. worked on a series of zinc-and manganese-co-doped MNPs. The highest hyperthermia rate was reached by a manganeseonly doped and Zn doping decreased hyperthermia property [21] . Fabris et al. suggested a new method to control the heat generation of magnetic colloids by changing the shell composition of bimagnetic nanoparticles [22] . Engelmann et al. worked on the relation between the agglomeration of nanoparticles and heating efficiency. Results showed An increase of up to 23% for agglomerated MNPs in suspension and a decrease up to 28% for mixed states of immobilized and agglomerated MNPs [23] . Mai et al. showed it is possible to prevent the aggregation of the nanoparticle by simple irradiation of copper-mediated polymerization with ultraviolet light [24] . Esmaeili et al. suggested dendrimers as a material for MNPs synthesis because of their small size, well-defined globular shape and monodispersed. As a result of using dendrimers, superior heat generation was achieved [25] . Ota et al. worked on harmonic signals and measured them with respect to different MNPs structures to examine the key effects of particle structures on magnetization harmonics. A dependency between core size and anisotropy of MNPs and harmonics was observed in their work [26] . The distribution of MNPs in tissue after injection is another parameter that has effects on hyperthermia efficacy. As mentioned before, raising the temperature of tumor cells lead them to necrosis and it is highly related to the position of MNPs after the injection, so it is vital to know about MNPs distribution and concentration after injection. For a long time, MNPs were believed to be distributed uniformly after injection, but it is now been established that the concentration is not uniform after injection [2, 27, 28] . There are some parameters that can cause a non-uniform distribution. Chen et al. investigated the effects of changing the rate of the infusion, cannula size, and the concentration of infusate on the fluid dispersion after injection. They found that higher infusion rates can cause more leak back of the infusate than lower infusion rates. It was also claimed that varying the concentration of the infusate does not affect the volume of the distribution [29] . Salloum et al. examined some parameters and their effects on MNPs distribution such as infusion flow rate and agarose gel concentration. It was concluded that the infusion flow rate and concentration could affect the distribution of MNPs. By decreasing the infusion rate, distribution becomes more spherically shaped [28] . Navidbakhsh et al. also investigated the effect of infusion rate on MNPs distribution with the help of computer simulation in an agarose gel. They developed a model for temperature distribution [30] . Kandala et al. used six different MNPs distributions to simulate heat transfer in tissue during hyperthermia by magnetic fluid. It was concluded that in concentrated nanoparticles distribution, modulated power heating makes less damage on surrounding normal tissue compared to constant power heating [31] . Gu et al. examined effects of heating on MNPs concentration distribution in prostatic cancer tumors. It was found that MNPs distribution volume in high heat generation rate was 10 percent smaller in the heating group, while in low heat generation rate it was 95 percent larger in the heating group [32] . In the hyperthermia cancer treatment by MNPs, several parameters may affect the results of treatment such as infusion rate, nanoparticle concentration in the magnetic fluid, nanoparticles diameter, coating type and surface charge of nanoparticles, type of the tissue and type of the needle used for injection and also blood flow. Among these parameters, the infusion rate has more effects, and it is also controllable. Moreover, blood flow can play a vital role in heat transfer during the treatment. Since the thermocouples had errors while a magnetic field was applying, simulation of heat transfer could help scientists to determine. Different investigations showed the effectiveness of using numerical methods for the estimation of the temperature [33] [34] [35] . To investigate the effects of infusion rate on distribution and as a consequence on hyperthermia cancer treatment, a fluid Hamilton syringe with a maximum reservoir volume of 250 l was used to inject fluid into the tissue. A 26-gauge Hamilton needle was also used to avoid the backflow of magnetic fluid. For injecting MNPs into the tissue at a constant rate, a syringe pump was needed. The only available pump was the Te-331 syringe pump with the retaining clip that only was able to handle large syringes with a size between 10 and 50 ml, so a 50 ml syringe was prepared and the Hamilton syringe was fixed inside it to be put into the syringe pump. Subsequently, since input numbers for the pump showed the infusion rate and volume of injections per the 50 ml syringe, a conversion between the 50 ml syringe and the fixed Hamilton syringe must be obtained, which was achieved by taking some tests and determining the time of each test. By performing these experiments, the ratio of 1 to 160 was obtained and formulated as Eq. 1 for the pump infusion rate and the actual rate of Hamilton syringe injection. In Eq. (1), = 9.6 and V H and V P denote the input velocity for the Hamilton syringe (μl/min) and input velocity for the infusion pump (ml/hour), respectively. MoS 2 ∕CoF 2 O 4 that had been examined and used in previous in-vivo tests was prepared and considered as magnetic fluid. MoS 2 had semi-conductor properties in addition to a large surface area, and introduction of CoF 2 O 4 on the surface of MoS 2 led to create hetero-structures. Results of previous investigations showed an improvement in electromagnetic absorptions of hybrids of MoS 2 and CoF 2 O 4 [36] . Moreover, no sign of inflammation or necrosis in mice liver was reported in the injection of MoS 2 ∕CoF 2 O 4 . The size of MNPs was also an important parameter in hyperthermia since it affected the specific loss power (SLP) [37] . MNPs with the size of 17 ± 4nm was considered since it could provide a sufficient SLP for hyperthermia treatment and was dispersed in aqueous very well [36] . The chicken breast tissue was used in this study because of its availability, low cost and brightness. The brightness of tissue was important since this brightness made it possible to determine the distribution of a substance after injection by using image processing techniques. Using chicken tissue for in-vitro investigation of cancer was reported before [38] [39] [40] . In fluid injection, 14 different slices were selected for testing various infusion rates between 1.25 l min and 104.06 l min (Maximum infusion rate produced by the infusion pump). To have comparable results with previous experimental works, a range of 10 l min to 40 l min was considered for infusion rate. Each infusion rate was repeated three times at best because of limitation of the available nanofluid. For cutting the tissue, a standard surgical blade was used. The cutting section was considered perpendicular to the injection site as shown in Fig. 1 . The crucial issue in this part was the probability of human errors in the cutting of the tissue that could be eliminated by using micron cutting devices. Access to these kinds of devices had been confronted with limitations. Therefore, the number of tests was considered 3 in each case, to minimize human errors in the cutting section. A Canon PC1680 camera was used to take photos, under a standard light condition and in the absence of sunlight. In the next step, all images were prepared for the image processing section. The preparation process began with cutting the part of the texture image off. Since the real tumors had an oval shape with a large diameter of 9 mm and a small diameter of 6 mm [36] , in the image processing section, a part of the original image which its dimensions were 2 cm * 2 cm, was cut to a rectangle of 9 mm * 6 mm to make the results of this study comparable with previous in-vivo works. A Schematic of these steps is illustrated in Fig. 2 . The MATLAB software was deployed for image processing in this study. To achieve the distribution of MNPs, all color images were converted into black-and-white images. In the black-and-white pictures, a threshold between 0 and 1 was determined for the MATLAB software to delete points of images that were brighter to separate darker points as nanoparticles position. By extracting the coordinates of all black pixels, the distribution of MNPs was obtained. Thermocouples had errors while a magnetic field was applying. Due to limitations in measuring the temperature in an experimental trial, simulation of heat transfer was deployed to determine the temperature of any point in the tissue which was warmed by nanoparticles. The role of blood flow on heat transfer was also investigated using numerical methods. A two-dimensional geometry was considered to make the simulation results comparable with future laboratory and animal works for the tumor. Based on previous animal works, an ellipse with a large diameter of 9 mm and a small diameter of 6 mm was created in COMSOL software using geometry tools. Around this ellipse, another ellipse with a large diameter of 19 mm and a small diameter of 16 mm was considered to represent the healthy tissue surrounding the tumor. The geometry is schematically illustrated in Fig. 3A . After creating the desired geometries, coordinates of points obtained by the MATLAB software as the distribution of MNPs were imported into the COMSOL software and added to the geometry. An example of the geometry after adding points as nanoparticles is illustrated in Fig. 3B. Properties of materials used in this simulation are presented in Table 1 . Although there was no blood in the chicken breast tissue, blood properties were also added in the model ( Table 2 ). The reason for considering the properties of the blood was to investigate its effect on the heat transfer and also made the numerical model more realistic. The point heat source was considered 0.016 W m . This consideration was based on having a maximum temperature of about 46 °C in the uniform distribution which was a standard value in clinical studies [7] . In the simulation of physics of the model, Pennes' bio-heat equation, as illustrated in Eq. (2), was considered as the governing equation. The reason was the correspondence of this equation with experimental studies [45] . In Eq. (2), , ∁, T, k denote the tumor/healthy tissue density, tumor/healthy tissue specific heat, local tumor/healthy tissue temperature and thermal conductivity of the tumor/ healthy tissue, respectively. t denotes time and b , ∁ b , b , T b denote blood density, specific heat of blood, blood perfusion rate and blood temperature, respectively. Q m denotes the metabolic heat generation rate in the tumor/healthy tissue and Q p denotes the point heat source in the tumor/healthy tissue. The Arrhenius equation was also deployed in the simulation section. This equation indicates the amount of necrotized tissue after a certain period of time and is expressed in Eq. (3) [46] : In Eq. 3, A, E a , R u denote pre-exponential factor, activation energy, and universal gas constant, respectively. For boundary conditions, the continuity of temperature and heat flux in the boundary between tumor and healthy tissue was added to the model (Eq. (4) and Eq. (5)): Which indicated the equivalence of temperature and flow flux in the boundary between the tumor and the healthy tissue. The second boundary condition was related to the outer surface temperature of the healthy tissue, which was constant and equal to body temperature or 37 °C. This boundary condition can be expressed as follows: The only initial condition that had been used in this study was the initial temperature of the healthy tissue and tumor, which was considered 37 °C and can be expressed as follows: In Eq. (4) to Eq. (7), a 1 = 4.5 mm, a 2 = 9.5 mm, b 1 = 3 mm and b 2 = 8 mm. The bio-heat transfer equation could be summarized and solved using finite element methods (FEM). The finite element technique was a powerful method in numerical solutions that could solve Eq. (2) with boundary and initial conditions Eq. (4) to Eq. (7). FEM divided a complex system into smaller parts called finite elements to solve a problem. This could be achieved by implementing space discretization, which was accomplished by building a mesh for the model. An algebraic equation system could be ultimately derived from the finite element method formulation of a boundary value problem. These equations were then grouped into a larger system that fit the whole model. Using the calculus of variations, the FEM approximated a solution by minimizing a function of errors. The overall time of 1200 s was considered in simulations since laboratory animals used in clinical trials could be anesthetized and treated up to 20 min. The Timestep of 1 s was also considered to have the temperature and the percentage of necrosis tissue in each second. A triangular mesh was used as a mesh type and generated by tools provided in COMSOL software. Each triangle had three sides with a specific magnitude. The maximum and minimum magnitude was considered for element's (sides of triangles) size in different mesh types to control the total number of triangles and as a result the resolution of the final results. To ensure the results were not affected by the time step and mesh size and were accurate, the maximum temperature was calculated for different mesh sizes at a constant time step of 1 s. Results represented an error of about 0.005 °C that could be ignored. The maximum temperature was also calculated for different time steps at constant normal mesh size. In this section, no error was observed. Results of this section which are illustrated in Tables 3 and 4 , showed the independence of the results from the mesh size and the time step. Considering the results of this section, the time step of 1 s and normal mesh type were chosen for simulation. Numerical models could be validated by comparing their results with an experimental or analytical investigation. This was a common method to validate a numerical model and authors had used it in other numerical researches such as coronavirus area too [47] . It was vital to have the same properties, initial and boundary conditions and set up to be able to compare the results. The numerical model used in this study was verified against the results of Lin and Liu [48] , for the radial distribution of temperature in the tissue with the same tissue property, initial, and boundary condition. The model was developed then with an elliptical shape instead of the circular but with the same initial and boundary conditions as Lin and Liu's model. A comparison of the simulation with the results of Lin and Liu is illustrated in Fig. 4 . In addition to verifying the numerical model with an analytical investigation, independence of the model from time step and mesh size was also checked in the previous section to be sure about the validation of the results. Increasing the infusion rate leads to an increase in nanoparticle penetration in tissue. Results show congruence with the results of Attaluri's work [49] . However, the notable point is the magnitude of the increase. According to Attaluri's work, with an increase in the infusion rate from 5 l min to 10 l min , an increase of 50% in penetration depth occurs, while in current experiments, with an increase in the infusion rate from 4 l min to 40 l min or in other words, with ten Fig. 4 Radial distribution of temperature in the tissue after 600 s, a comparison between present work and Lin and Liu's [48] length of dispersion(mm) times the infusion rate, there is a 1-mm increase in penetration depth, which represents a nearly 17 percent increase in penetration depth. The reason may be the type of tissue tested and the difference in the tests performed with the 2011 Attaluri's work. Another reason could be the use of different procedures to examine the distribution of MNPs since Attaluri used CT scan imaging techniques, while in the present study, cutting and direct imaging method has been used. As shown in Fig. 5 , the best distribution occurs in the lowest infusion rate of 1.25 l min , which is almost circularly distributed and has similarity with the results of 2008 Salloum's work on the agarose gel. Comparing different samples under the microscope reveals a lower concentration of MNPs could be achieved at higher infusion rates. This decrease in concentration could happen since, in higher infusion rates, nanoparticles are distributed in a bigger area. The effects of infusion rate on MNPs concentration are illustrated in Fig. 6 . By having the results obtained in the simulation of heat transfer section, and considering the gradual increase in penetration depth by increasing the infusion rate and consequently excretion of nanoparticles from the border of the tumor and healthy tissue, it is determined that the tissue necrosis happens in healthy tissue at higher rates. This necrosis of healthy tissue occurs while some parts of the tumor do not feel any increase in temperature and consequently will not be necrosis at all. Therefore, it could be concluded that higher infusion rates should be avoided to have a circular shape distribution and also to prevent the penetration of nanoparticles into healthy tissue. By obtaining the specific distribution of magnetic nanoparticles in the tissue, the maximum temperature in different cases is higher than the maximum temperature in a uniform distribution (less than 46 °C). This difference indicates the impact of the distribution of MNPs in tissues on the maximum produced temperature. Due to the maximum temperatures obtained in each section, it can not be concluded that there is a correlation between the infusion rate and the maximum temperature. As shown in Table 5 , for an infusion rate of 1.25 l min , the maximum temperature is about 51.985 °C, which by increasing the infusion rate up to 10 l min , reaches about 51.033 °C. The reason can be reducing the distribution area, or in other words, concentration of nanoparticles in a smaller area which is consistent with the results of the experimental section (Concentration of MNPs for different infusion rates). The lack of uniform distribution of nanoparticles also affects the magnitude of maximum temperature at a specific point. For instance, in Sample 1 of the fifth test, the maximum temperature is 49.189 °C, while in the similar sample of the same test (Sample 3 of the fifth test), the maximum temperature is about 50.953 °C, which is about 1.8 °C of difference. Also, other reasons such as the impact of the image processing and the type of cutting and photographing may also affect the results in each section, but in general, and for samples after the fourth test, which the distribution width is almost constant, it shows a very slight decrease in maximum temperature that is similar to 2011 Attaluri results. In other words, maximum temperature occurs in minimum infusion rates (1.25 l min ). Another noteworthy point obtained from the results is the role of blood flow in heat transfer. As it is shown in Table 5 , blood has a significant effect on heat transfer in tissue, and as a result, magnetic fluid hyperthermia cancer treatment. Maximum temperature can be increased up to 9 °C (test8) in the absence of blood flow in the tissue. Results of the effects of blood on heat transfer are divided into three sections. Having no blood in both healthy and tumor tissue, having blood just in healthy tissue, and finally having blood in both healthy and tumor tissue. The last point to be taken from the temperature distribution in the tissue is the maximum temperature point, which by increasing the infusion rate, recedes from the tumor center (which is the injection site) and approaches the borders of the tumor and healthy tissue (Fig. 7. ). Cancer is one of the most vital causes of death in recent years, and several methods have been suggested for treating it. Hyperthermia is one of the emerging methods of treatment that is more accurate and has fewer side effects compared to other treatments. Although many investigations have been conducted in this area, there is a lack of knowledge in some parameters such as the role of blood flow on MNPs hyperthermia, diffusion of MNPs, temperature distribution during treatment, and tumor necrosis after treatment yet. In this study, a different technique is used to measure the depth of penetration compared to previous works and a more precise ratio between penetration depth and infusion rate is achieved as a result. The results of experiments on MNPs show that by increasing the infusion rate, the penetration depth increases. With an increase in the infusion from 4 l min to 40 l min , there is a 1-mm increase in penetration depth (nearly 17 percent increase). This increase in infusion rate also has a significant impact on the concentration of MNPs in tissue. Increasing the infusion rates leads to a lower concentration of MNPs. The numerical simulation reveals that blood flow in tumor and healthy tissue plays an important role in heat transfer in the tissue (in addition to mass transfer of MNPs out of tumor area made by blood flow, which has been investigated in previous investigations). The numerical simulation also showed that there is not a clear and precise relationship between increasing the penetration depth and the maximum temperature created in the tissue, but it can be concluded that by increasing the infusion rate and consequently increasing penetration depth, the maximum temperature decreases slightly which means maximum temperature occurs by the infusion rate of 1.25 l min . Finally, the results of the necrosis section reveal that by increasing the infusion rate, necrosis has occurred in healthy tissue rather than cancerous tissue, which should be avoided by reducing the infusion rate. Moreover, the location of the necrosis recedes from the tumor center and approaches the borders of the tumor and healthy tissues by increasing the infusion rate. Funding The authors did not receive support from any organization for the submitted work. Data availability Data are easily accessible in any time. Code availability Codes are available at any time. Ethics approval Not applicable. The authors agree to participate. The authors agree to publication this article. The authors have no conflicts of interest to declare that are relevant to the content of this article. Nanoparticle distribution and temperature elevations in prostatic tumours in mice during magnetic nanoparticle hyperthermia Transport phenomena in biological systems The effect of thermotherapy using magnetic nanoparticles on rat malignant glioma Enhancement of tumor thermal therapy using gold nanoparticle-assisted tumor necrosis factor-α delivery Thermotherapy of prostate cancer using magnetic nanoparticles: feasibility, imaging, and three-dimensional temperature distribution Tumour cell toxicity of intracellular hyperthermia mediated by magnetic nanoparticles Anti-cancer effect of hyperthermia on breast cancer by magnetite nanoparticleloaded anti-HER2 immunoliposomes Hyperthermia in cancer therapy: where are we today and where are we going? Thermophoretic isolation of circulating tumor cells, numerical simulation and design of a microfluidic chip Analytical solutions of fractional bioheat model in a spherical tissue Analytical estimations of temperature in a living tissue generated by laser irradiation using experimental data An analytical study on the fractional transient heating within the skin tissue during the thermal therapy A Topical Review of Magnetic Fluid Hyperthermia Characterization of Magnetite/Iodized Oil Magnetic Fluid for Embolization and Hyperthermia Application Synthesis and magnetic hyperthermia properties of zwitterionic dopamine sulfonate ligated magnesium ferrite and zinc ferrite nanoparticles Heating ability of magnetic nanoparticles with cubic and combined anisotropy Cell-Promoted Nanoparticle Aggregation Decreases Nanoparticle-Induced Hyperthermia under an Alternating Magnetic Field Independently of Nanoparticle Coating, Core Size, and Subcellular Localization Polyethylene glycol-coated porous magnetic nanoparticles for targeted delivery of chemotherapeutics under magnetic hyperthermia condition Coprecipitation synthesis of (Zn-Mn)-co-doped magnetite nanoparticles and their application in magnetic hyperthermia Controlling the dominant magnetic relaxation mechanisms for magnetic hyperthermia in bimagnetic core-shell nanoparticles Magnetic Relaxation of Agglomerated and Immobilized Iron Oxide Nanoparticles for Hyperthermia and Imaging Applications Thermo-Responsive Iron Oxide Nanocubes for an Effective Clinical Translation of Magnetic Hyperthermia and Heat-Mediated Chemotherapy Dendrimer functionalized magnetic nanoparticles as a promising platform for localized hyperthermia and magnetic resonance imaging diagnosis Effects of size and anisotropy of magnetic nanoparticles associated with dynamics of easy axis for magnetic particle imaging An in-vivo experimental study of temperature elevations in animal tissue during magnetic nanoparticle hyperthermia Controlling nanoparticle delivery in magnetic nanoparticle hyperthermia for cancer treatment: Experimental study in agarose gel Variables affecting convection-enhanced delivery to the striatum: a systematic examination of rate of infusion, cannula size, infusate concentration, and tissue-cannula sealing time Evaluation of the effects of injection velocity and different gel concentrations on nanoparticles in hyperthermia therapy Temperaturecontrolled power modulation compensates for heterogeneous nanoparticle distributions: a computational optimization analysis for magnetic hyperthermia Nanoparticle Redistribution in PC3 Tumors Induced by Local Heating in Magnetic Nanoparticle Hyperthermia: In Vivo Experimental Study Thermal response of cylindrical tissue induced by laser irradiation with experimental study Nonlinear analysis of dual-phase lag bio-heat model in living tissues induced by laser irradiation Theoretical analysis of thermal damages in skin tissue induced by intense moving heat source Distribution of "molybdenum disulfide/cobalt ferrite" nanocomposite in animal model of breast cancer, following injection via differential infusion flow rates Sizedependant heating rates of iron oxide nanoparticles for magnetic fluid hyperthermia Thermoacoustic and photoacoustic tomography of thick biological tissues toward breast imaging The chicken model of spontaneous ovarian cancer Mechanical characterization of human brain tumors from patients and comparison to potential surgical phantoms Effects of temperature distribution in the tissue around the tumor on the quality of hyperthermia Thermoregulation in the canine prostate during transurethral microwave hyperthermia, part II: blood flow response Optimizing heat shock protein expression induced by prostate cancer laser therapy through predictive computational models A comprehensive tissue properties database provided for the thermal assessment of a human at rest Analysis of tissue and arterial blood temperatures in the resting human forearm Thermal modeling of lesion growth with radiofrequency ablation devices Centrifugal isolation of SARS-CoV-2: numerical simulation for purification of hospitals' air Estimation for the heating effect of magnetic nanoparticles in perfused tissues Quantification of Nanostructure Distribution in Tissue Using MicroCT Imaging