Next Article in Journal
Experimental Study on Mix Proportion Parameter Optimization of Cement Anchoring Material
Next Article in Special Issue
An Experimental and Computational Study of the High-Velocity Impact of Low-Density Aluminum Foam
Previous Article in Journal
Reversible and Irreversible Processes in Drying and Wetting of Soil
Previous Article in Special Issue
Compressive Behaviour of Closed-Cell Aluminium Foam at Different Strain Rates
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Mathematical Modeling of Breast Tumor Destruction Using Fast Heating during Radiofrequency Ablation

Department of Computational Mechanics and Engineering, Silesian University of Technology, 44-100 Gliwice, Poland
Materials 2020, 13(1), 136; https://doi.org/10.3390/ma13010136
Submission received: 3 November 2019 / Revised: 20 December 2019 / Accepted: 23 December 2019 / Published: 28 December 2019

Abstract

:
In oncology, hyperthermia is understood as a planned, controlled technique of heating cancerous changes in order to destroy their cells or stop their growth. In clinical practice, hyperthermia is used in combination with radiotherapy, chemotherapy, or immunological therapy. During the hyperthermia, the tissue is typically exposed to a temperature in the range of 40–45 °C, the exception is thermoablation, during which the temperatures reach much higher values. Thermoablation is characterized by the use of high temperatures up to 90 °C. The electrode using the radiofrequency is inserted into the central area of the tumor. Interstitial thermoablation is used to treat, among others, breast and brain cancer. The therapy consists of inducing coagulation necrosis in an area that is heated to very high temperatures. Mathematical modeling is based on the use of a coupled thermo-electric model, in which the electric field is described by means of the Laplace equation, while the temperature field is based on the Pennes equation. Coupling occurs at the level of the additional source function in the Pennes equation. The temperature field obtained in this way makes it possible to calculate the Arrhenius integral as a determinant of the destruction of biological tissue. As a result of numerical calculations regarding the temperature field and the Arrhenius integral, it can be concluded that, with the help of numerical tools and mathematical modeling, one can simulate the process of destroying cancerous tissue.

1. Introduction

Cancer diseases are one of the main causes of death in both the societies of highly developed and low developed countries. In recent years, their intensity has increased particularly, which caused cancer to become not only a serious health problem, but also a social and economic problem [1].
Breast cancer is the second most common malignancy in the world, right after lung cancer. In the female population, it is the most frequently occurring malignant tumor. According to data from the World Health Organization, in 2018, 2.1 million new cases of this disease were diagnosed, which is 11.6% of all cancers in women [2].
Breast cancer is a malignant tumor arising from glandular tissue cells that develops locally in the breast and provides metastases to the lymph nodes and internal organs (e.g., lungs, liver, bones, and the brain). This cancer is the most common malignancy in women. It accounts for about 23% of all cases of malignant tumors in women and about 15% of deaths for this reason. It is estimated that 2.1 million women are diagnosed with breast cancer every year around the world, and about 620,000 die of it [2,3].
Breast cancer is the biggest oncology problem in developed countries; it is also a growing problem in developing countries. An example is the comparison of the number of patients and deaths: in developed countries, there are slightly more than half of the cases (50.1%), while the majority of deaths are recorded in less developed countries (58.7%) [4,5].
Medicine increasingly uses the possibility of performing treatments with methods of limited invasiveness (so-called minimally invasive procedures). One of such techniques is thermoablation, i.e., destruction of pathological tissues by means of high temperature. It has a special application in oncology. Based on the energy sources and their delivery, the thermal ablative techniques can be classified into five major categories, namely, radiofrequency ablation (RFA), laser ablation (LA), microwave ablation (MWA), cryoablation, and high-intensity focused ultrasound (HIFU) ablation [6]. An interesting review of the literature on thermoablation of a breast tumor using the above-mentioned methods can be found, among others, in [7].
The most common thermoablation, using electromagnetic waves, consists of inserting the electrode in the form of a needle into the central site of the tumor. This can be done traditionally by opening the abdominal cavity (e.g., in the case of a liver tumor) or through the skin during laparotomy or laparoscopy. The process is carried out under the control of ultrasonography, which is widely recognized as the best examination in the imaging of the surgery during its performance, or under the control of a computed tomography.
The needle-electrode used in this technique transmits high frequency current (300–500 kHz, which corresponds to the range of long and medium radio waves, usually 460 kHz). In this way, thermal energy is generated. The high electric intensity, resulting in the highest temperature, is at the tip of the electrode. The high temperature obtained in this way allows carbonization of carbohydrates contained in the tissues. The degree of destruction depends on the temperature obtained and the time of current flow. Temperature between 50 and 100 °C causes coagulation necrosis. Temperatures up to 45–50 °C induce necrosis and contribute to the destruction of cancer tissue due to dehydration, thickening, and denaturation of intracellular proteins, as well as destroying cell membranes. Destruction of the lesion with a diameter of approx. 3 cm, and therefore optimal in terms of the size of the area undergoing treatment, requires the influence of current in 10 to 15 min [8].
The beginning of treatment with a high temperature is at the end of the 19th century. The breakthrough came in 1891, when Jacques-Arsène d’Arsonval discovered that the flow of high-frequency current through the tissues leads to heat in them, without causing muscle contraction. In 1990, the destruction of liver tumors in animals was initiated by placing inside the electrode connected to the power generator. In the case of human liver tumors, thermoablation entered clinical practice in 1995, and, five years later, it was officially approved in the United States by the Federal Drug Administration (FDA) for use in cases of unresectable liver tumors [9].

2. Materials and Methods

2.1. Introduction

In the case of electric field impact, numerical modeling consists of solving a coupled problem, and, more specifically, analysis of electro-thermal coupling. A tissue model is selected for the mathematical description of heat transfer, namely the classic Pennes model [10] (the most commonly used mathematical model describing the process of heat transfer in living organisms, on a macro scale, which is commonly used in particular to predict temperature changes during hyperthermia, or hypothermia treatments, i.e., broadly understood thermotherapy, e.g., [11,12,13,14,15,16]). The Pennes equation is actually a Fourier heat conduction equation with an additional source component that takes into account the global heat flow effect between blood and the surrounding tissue. In 1948, Pennes [10] proposed a model of the heat transfer in which the presence of a system of small blood vessels and metabolic processes occurring in tissues was taken into account by adding appropriate source components to the Fourier equation. The source function proposed by Pennes is the sum of two components called the perfusion source and the metabolic source. The idea of this equation is that blood at a temperature equal to blood temperature in the aorta flows through the tissue at a certain speed, achieving thermal balance in it, and then accumulates in the veins. The rate of blood perfusion is therefore a key parameter in the calculation of the bioheat transfer. Pennes differential equation is presented in Section 2.3. In turn, the mathematical description of the electric field is based on the Laplace equation. The Laplace differential equation is a second order linear partial differential equation. The name of the equation comes from the name of Pierre Simon de Laplace, who formulated it in the 18th century. The Laplace equation describes many physical processes. In a typical situation, the function contained in it is interpreted as a density of a certain magnitude (e.g., chemical concentration, temperature, potential). This equation expresses the following property of the potential field: divergence of the vector potential field (i.e., potential gradient), in the absence of the source is equal to zero. The Laplace equation is a special case of the Poisson equation expressing an analogous relationship in the case of the field sources. Therefore, it describes many processes occurring in nature, e.g., in electrostatics, electrostatic potential in the absence of electric charge, in thermodynamics, it describes the steady-state heat transfer process, in fluid mechanics, it describes the steady flow of liquids, in electrodynamics, it describes current flow in extended media, and, in mechanics, it describes deformations of elastic material. Simplification involving the use of the above equation to the mathematical description of the electric field is related to the following assumption. During medical procedures based on radiofrequency ablation caused by the electric field, a current with frequencies in the range of tens to several hundred kHz is induced. In such a frequency range, the electromagnetic wave length is much larger than the depth of the human body; therefore, the current flow is via electrical conductivity and can be analyzed as a so-called quasi-static formulation, which allows coupling the electrical and thermal conductivity problem [8,17,18,19,20,21,22]. To determine the internal heat source resulting from the influence of the electric field, it is necessary to know the intensity of the electric field, which is dependent on the electric conductivity coefficient.
The coupling of these two problems takes place by means of the weak formulation, which consists of taking into account the source function resulting from the interaction of the electric field in the Pennes equation, describing the heat transfer in the biological tissue.
Besides determining the temperature field in the biological tissue affected by the tumor and subjected to external thermal action, it is also necessary to estimate the degree of tissue destruction under the influence of high temperature. In the literature, e.g., [15,20,23,24], we can find information that the temperature above 45 degrees of Celsius is the so-called ablation temperature, i.e., causing irreversible destruction of biological tissue (the duration of the heating process is also important). It is obvious that tissue destruction is reflected in the values of its parameters, which may assume nonlinear values depending on the processes involved. Very often, the perfusion coefficient is treated as an indicator of the main tissue damage. A simple solution is the dependence of the perfusion coefficient on the temperature, assuming that exceeding the temperature at the point under consideration above 45 degrees of Celsius causes the disappearance of blood perfusion. In the literature, we can find relations describing the relationship between temperature and perfusion, as well as between perfusion and tissue destruction. One such method is to determine the value of the Arrhenius integral, which will be discussed in more detail in the next section. At this point, however, several facts should be presented. The "Arrhenius Integral" or "Arrhenius Equation" was physical justification and interpretation in 1889 by Svante Arrhenius, a Swedish chemist. Arrhenius performed experiments that correlated chemical reaction rate constants with temperature. After observing that many chemical reaction rates depended on the temperature, Arrhenius developed this equation to characterize the temperature-dependent reactions. The Arrhenius equation is successfully used not only in estimating the degree of destruction of biological tissue, as is the case in this article, but also in issues of physical chemistry, e.g., [25], or issues related to the analysis of the behavior of metals or alloys during high-temperature deformation, e.g., [26], and I suppose in many other problems.

2.2. Breast Modeling

Both men and women have breast glands, but there are many more tissues in the women breast. Each breast is located on the chest muscle, called the pectoral muscle. The feminine breast covers a relatively large area. It stretches from the place just below the collarbone, to the armpit, and across the bridge. They are made up of lobules (milk-producing glands) and ducts (tubes that carry milk to the nipple). These are surrounded by glandular, muscular, and fatty tissue (cf. Figure 1a [27]).
The breast is a very complex and largest human gland. During human life, it undergoes more changes than any other part of the human body—from birth through puberty, pregnancy and breastfeeding, until menopause. It contains from 15 to 20 glands called lobes. The lobes are connected to the nipple with several tubes called ducts that carry milk to the nipple. The breast consists of lymph nodes and vessels that allow the transmission of lymph fluid and white blood cells. The rest is mainly fatty tissue. The breast, like any other part of the body, is made up of billions of microscopic cells. These cells usually reproduce in an orderly way—new healthy cells divide and then replace those that have died. However, sometimes, it happens that cells develop abnormally (mutations appear). This happens when the genes responsible for checking that cells are reproducing correctly do not detect changes or mutations. In this situation, pathological cells divide and multiply, and their growth is sometimes quite fast. At this stage, a growth may not be cancerous. At this point, we can deal with the so-called "non-invasive tumor", which is located in the lobe or the duct. When the tumor can infect the surrounding tissue, it is treated as cancerous. These types of cancers require immediate treatment because, if they continue to grow and propagate, they will be life-threatening. Breast cancer most often develops in glandular tissue (cf. Figure 1b [28]).
For the purpose of numerical calculations, a multi-layered breast geometric model is prepared. Among the various configurations, the following types of models can be distinguished: extremely dense (ED), heterogeneously dense (HD), scattered fibro-glandular (SFG), and predominantly fatty (PF) [29,30]. Figure 2 presents the cross-sectional view of the developed breast models with the percentage of tissue compositions for each layer.
Because the number of breast cancer cases among women increases after the age of 35, the HD model was selected for analysis (see: Figure 2b), whose volume share of individual layers is appropriate for this age. Following [30], the percentage of layers of HD model can be determined at the level: muscle layer—20%, glandular layer—60%, and fat layer—20%.

2.3. Mathematical Modeling of Electric and Temperature Fields

Because the wavelength of the RF current in tissues is much greater that depth of a human body, the quasi-static electric field approximation is applied. Using the quasi-static formulation, the electric field intensity E (V/m) inside the biological tissue for 2D problem can be calculated as follows:
E e ( x 1 , x 2 ) = φ e ( x 1 , x 2 ) = [ φ e ( x 1 , x 2 ) x 1 φ e ( x 1 , x 2 ) x 2 ] T
where subscript e denotes the sub-regions of breast, electrode and trocar (i.e., 1—muscle, 2—gland, 3—fat, 4—tumor, 5—electrode, 6—trocar), {x1, x2} is a geometrical co-ordinates and φe (V) is an electric potential, while the external heat generation Qeext(x1, x2) due to the electric heating for healthy tissue and tumor is defined as follows:
Q e e x t ( x 1 , x 2 ) = σ e 2 | E e ( x 1 , x 2 ) | 2 = σ e 2 ( φ e ( x 1 , x 2 ) x 1 + φ e ( x 1 , x 2 ) x 2 ) 2 ,
where σe (S/m) is an electrical conductivity.
The electric potential φ e ( x 1 , x 2 ) inside the breast, cancer, and internal applicator is described by the system of Laplace equations [16,19,21,31]:
x 1 , x 2 Ω e : σ e 2 φ e ( x 1 , x 2 ) = 0 .
For a heat transfer process in biological tissue, the Pennes model has been proposed e.g., [10,12,14,16,17,19,21,32,33,34,35]
x 1 , x 2 Ω e : c e ρ e T e ( x 1 , x 2 , t ) t = λ e ( T ) T e ( x 1 , x 2 , t ) + k e [ T B T e ( x 1 , x 2 , t ) ] + Q m e t e + Q e e x t ( x 1 , x 2 ) ,
where t (s) denotes time, ρe (kg/m3) is a density, ce (J/(kgK)) is a specific heat, λe (W/(mK)) is a thermal conductivity, Te (K) is a temperature, ke = GBecBρB (W/(m3K)) is a perfusion rate (GBe (1/s) is a blood perfusion coefficient, cB (J/(m3K)) is a specific heat of blood, ρB (kg/m3) is a blood density), TB is a supplying arterial blood temperature (baseline physiological temperature), and Qmet e (W/m3) is a metabolic heat source. The material properties used in calculations are collected in Table 1, based on [8]. At this point, it should be mentioned that the blood perfusion coefficient GB is a key parameter in the calculation of the bioheat transfer problem. The value of this parameter is different for different types of tissues (e.g., fat, skin, liver, heart, kidney, etc.) and increases by an order of magnitude in conditions different from the natural state (e.g., tumor).
The electrical and thermal conductivity are assumed as linear functions of temperature, approximated by Equations (5) and (6), respectively, for both healthy and tumor tissues [8,36]
σ ( T ) = σ 0 [ 1 + 0.02 ( T T B ) ]
λ ( T ) = λ 0 + 0.0013 ( T T B )
where σ0 and λ0 are the baseline electrical and thermal conductivities, respectively, at baseline physiological temperature, TB = 37 °C (see: Table 1).
Differential equations which describe the electric (cf. Euqation (3)) and temperature (cf. Equation (4)) fields are supplemented by appropriate boundary conditions. Constant voltage boundary condition φ = U (U (V) is an electric potential of the electrode) is applied on the active part of the electrode, on the boundaries Γ1 (see: Figure 3) (simulating a ground pad) the zero volts (U = 0) is considered. On the remaining boundaries, an electric insulation is assumed. In the case of temperature field, on the boundaries Γ1, the body core temperature Tb = 37 °C is assumed. The outer surface of breast is assumed to be subjected to natural convection
λ T ( x 1 , x 2 , t ) n = α e f f [ T ( x 1 , x 2 , t ) T a m b ]
where αeff = 13.5 (W/(m2K)) represents the combined effective heat transfer coefficient due to convection, radiation, and evaporation [8,30,37,38], while the Tamb = 25 °C is an ambient temperature. On the contact surfaces between breast sub-regions and tumor, the ideal electric and thermal contacts are assumed. Equation (4) is also supplemented by the initial conditions:
  • the initial temperature of the electrode tine and trocar domains (indicated as Ω5 and Ω6) have been set to 25 °C, simulating the ambient room temperature condition
    x 1 , x 2 Ω 5 Ω 6 : T ( x 1 , x 2 , t ) = T 5 6 | t = 0 = 25 ° C
  • the initial temperature of the breast and tumor domains (indicated as Ω1, Ω2, Ω3 and Ω4) have been considered to be uniform and the same as the body core temperature of the human body
    x 1 , x 2 Ω 1 Ω 2 Ω 3 Ω 4 : T ( x 1 , x 2 , t ) = T 1 4 | t = 0 = 37 ° C

2.4. Arrhenius Scheme—A Model of Tissue Destruction

To estimate the degree of tissue destruction the Arrhenius integral, which describes the relationship between temperature and tissue damage, is used [39,40,41,42,43,44], where the relation between the tissue death with the temperature and the treatment time is given by
A r r ( x 1 , x 2 , t ) = 0 t f A exp [ Δ E R T ( x 1 , x 2 , t ) ] d t ,
where R (J/(molK)) is a universal gas constant (R = 8.3143), ΔE (J/mol) is an activation energy, A (1/s) is a pre-exponential factor, T(x1, x2, t) denotes tissue temperature at the point {x1, x2}, and time t f, while [0, tf] is a considered time interval. In the Arrhenius equation, the activation energy ΔE is defined as the minimum energy required to initiate the chemical reaction. The parameters A and ΔE are determined experimentally. For example, Table 2 summarizes the values of the above parameters for different types of tissues [45,46].
The main assumption of the Arrhenius scheme is that the thermal tissue damage is irreversible and total. It is assumed that irreversible tissue damage occurred, if
A r r ( x 1 , x 2 , t f ) 1
and the degree of tissue destruction as a result of thermal interaction corresponds to the share of damaged cells in the total volume of tissue VD, which is expressed by
V D ( x 1 , x 2 , t ) = κ 0 κ ( t ) κ 0 = 1 exp [ A r r ( x 1 , x 2 , t ) ] ,
where κ0 is the initial concentration of healthy cells, whereas κ(t) the concentration of healthy cells after time t (t—thermal exposure time). Based on the knowledge of cell concentration, the relationship that defines the degree of tissue damage is written as follows:
A r r ( x 1 , x 2 , t ) = ln ( κ ( t ) κ 0 ) .
In computational practice. it is more convenient to use the integral relationship (10). According to the relationship (12), for the threshold value Arr = 1, 63% of the cells are permanently and irreversibly damaged (thermally damaged). It should be noted that, for the value Arr = 4.6, the degree of tissue damage is equal to 99%.
In addition, the value of perfusion coefficient GBe (cf. Equation (4)) depends on the value of the Arrhenius integral. The equation describing the relationship between tissue damage should reproduce the initial increase in perfusion caused by hyperemia (vasolidation), and then its fall due to collapse of vasculature (closing the vessels). Therefore, according to the literature [8,33,39], the perfusion coefficient can be taken as the following function:
G B e ( A r r ) = { G B e 0 , for A r r = 0 ( 1 + 25 A r r 260 A r r 2 ) G B e 0 , for 0 < A r r 0.1 ( 1 A r r ) G B e 0 , for 0.1 < A r r 1 0 , for A r r > 1 ,
where GBe0 is the baseline blood perfusion coefficient (cf. Table 1).

3. Results

A heterogeneous (multi-layered) two-dimensional model of woman breast (see: Figure 2b) has been considered, in which a tumor of 2 cm diameter has been embedded. Dimensions of the domain considered, position of tumor center, and the position of internal applicator are presented in Figure 3a, while Figure 3b shows a fragment of the analyzed geometry with the marked control points. Table 3 contains the coordinates of the control points. The applicator consists of a stainless steel trocar and the nickel-titanium electrode tine. The applicator is located in the tumor area at an angle of 45° and its tip is in the center of the tumor. The electrical and thermophysical parameters of breast sub-regions, cancer, blood, electrode tip, and trocar are collected in Table 1, while the other parameters used for the calculations, discussed in the text, are collected in Table 4.
In order to solve the equations describing the potential of electric and temperature fields, the finite element method is applied using commercial software MSC MARC/MENTAT 2019 (64-bit). Bearing in mind a good representation of the geometry, and due to the expected large electrical and thermal gradients, the following mesh parameters were adopted: mesh density is equal to 17,704 quadrilateral, four-node elements. The computations have been performed using the personal workstation with a processor Intel® Core i7 CPU 950 @3.07 GHz with 12 GB RAM.
Figure 4 illustrates contour plots of electric potential for the exemplary voltage U = 10V, while Figure 5 shows the temperature distribution after 10 min during ablation heating. In Figure 6, the temperature history at the control points Pi (i = 1, 2, …, 5) (cf. Figure 3b), obtained for the electric potential U = 10V and exposure time equal to 10 min, are shown, while, in Figure 7, the Arrhenius integral courses at the control points are presented. The contour maps presented in Figure 4 and Figure 5 and heating curves (see: Figure 6) were obtained on the basis of numerical calculations. In Figure 8, the temperature history for different values of electric potential U = {10, 12.5, 15, 17.5, 20}V at a control point P1 are presented, while, in Figure 9, the Arrhenius integral courses at a control point P1 are shown. The exposure time is also equal to 10 min.
As can be seen in Figure 9, the voltage adopted at 12.5 V causes the Arrhenius integral to assume a value exceeding 4.6 at control point P4 (probability of tissue destruction equal to 99%). Therefore, Figure 10 shows the value of the Arrhenius integral for all control points for an assumed voltage of 12.5 V.
In Figure 11, the temperature history at the control points Pi (i = 6, 7, …, 11) (cf. Figure 3b), obtained for the electric potential U = 12.5 V and exposure time equal to 10 min, are shown, while, in Figure 12, the Arrhenius integral courses at these control points are presented.
Table 5 presents the value of the Arrhenius integral (cf. Figure 10), the probability of tissue destruction based on the formula (12), and temperature values for all control points after 10 min.

4. Discussion

Numerical modeling of the thermoablation process taking into account the influence of the electric field gives very good and promising results. Computer software using FEM allows simulation of the electro-thermal coupling process, and thus a simulation of breast tumor destruction using an internal applicator generating an electric field and causing the effect of thermal ablation. The mathematical model supplemented by the Arrhenius integral makes it possible to estimate the degree of destruction of cancerous tissue under the influence of high temperature.
As shown in the calculation results, the increase of electrical potential caused an increase of temperature within the analyzed area (cf. Figure 8). At the same time, an increase of temperature caused an increase of the value of the Arrhenius integral (cf. Figure 9).
In Figure 7, which shows the course of the Arrhenius integral for selected control points (see: Figure 3b), it can be observed that, only for point P4, the value of the Arr parameter approaches the value of the damage criterion (cf. Equation (11)). For the other control points, i.e., P1, P2, P3 and P5, the value resulting from Equation (10) is definitely less than 1.
The effect caused by the increase of the value of the electric potential on the active part of the electrode is shown in Figure 8 and Figure 9. The increase of electrical potential had a direct impact on the temperature rise and exposure time needed to obtain an Arrhenius integral value equal to 1.0 or 4.6. Analyzing Figure 9, it can be concluded that the optimal solution seems to be the choice of voltage equal to 12.5 V or 15 V, taking into account both the exposure time and the value of the Arrhenius integral. After about 10 min, tissue destruction determined by the location of the P4 control point can be seen with a 99% probability for U = 12.5 V, and after about 6.5 min for U = 15 V—(Arr > 4.6). If the potential value increased further, the Arrhenius integral increased very rapidly, which could also cause the destruction of healthy tissue surrounding the tumor.
Subsequent numerical simulations confirmed the above conclusions. Figure 10 shows the course of the Arrhenius integral at control points in the tumor for a voltage of 12.5 V. For points P1, P2, and, P4, the probability of destruction exceeds 80% (cf. Table 5), which corresponds to the temperature values obtained during the simulation. Of course, the exposure time to high temperature is also important.
The introduction of additional control points in healthy tissue, i.e., P6 (muscle), P7P9 (gland), P10 (on the skin surface), and P11 (fat) allowed for numerical calculation of heating/cooling curves at these points (cf. Figure 11) and estimating the degree of tissue damage based on the Arrhenius integral (cf. Figure 12). Based on numerical simulations, it can be stated that the voltage selection on the active part of the electrode, equal to 12.5 V, was carried out correctly. The criterion of destruction (cf. Equation (11)) was not exceeded at any of the additional points analyzed. Only at point P9 located in the glandular tissue above the tumor can it be observed that the Arrhenius integral value is close to 1.0 (probability of destruction about 60%). However, this point is located about 5 mm from the tumor tissue, and it can be assumed that it is an acceptable region.
The concentration of thermal energy in a specific area can be achieved by introducing paramagnetic nanoparticles into the tumor, e.g., [12,21,33], although such a procedure will cause an increase of temperature at the nanoparticle injection site but will not prevent possible damage to healthy tissue. An interesting solution seems to be the introduction into the tissue surrounding the tumor, nanoparticles with such thermophysical parameters that could constitute a kind of barrier against temperature increase, and thus protect healthy tissue from overheating. This issue will be analyzed in the future.
From the point of view of the possibility of supporting thermotherapeutic techniques, it is important to examine the impact of individual parameters occurring in the mathematical model on the degree of tissue destruction. Therefore, it is necessary to estimate, based on the sensitivity analysis methods, which parameter has the greatest effect on the computations’ results [47,48]. The results of calculations related to the sensitivity analysis were presented at the 4th Polish Congress of Mechanics and 23rd International Conference on Computer Methods in Mechanics [49].
One of the important parameters of the mathematical model of bioheat transfer is undoubtedly the perfusion coefficient because the flowing blood acts as a heat sink. The level of tumor perfusion varies drastically depending on the type of tumor, and the value GB = 0.0053 1/s used for calculations (cf. Table 1) is only one of the possibilities [8,50]. Therefore, it is important to examine the impact of this parameter on the obtained results of calculations. For this purpose, different values of perfusion coefficient appropriate for normoperfusion (GB = 0.0053 1/s), hypoperfusion (GB = 0.00265 1/s), and hyperperfusion (GB = 0.0106 1/s) were adopted. The above values for hypo- and hyperperfusion were obtained by twice decreasing and increasing the reference value of the GB parameter, respectively. Figure 13 and Figure 14 present the temperature history and the course of the Arrhenius integral at point P4 (cf. Figure 3b), in which, due to its location, the changes are best observed.
As could be expected, an increase of the perfusion coefficient caused a decrease of temperature during the numerical simulation of the thermoablation phenomenon (see Figure 13). The opposite situation occurred when the GB coefficient decreased and then the temperature increased. Changes in temperature distribution (cf. Figure 13) are also reflected in the course of the Arrhenius integral (cf. Figure 14). The above simulations indicate the importance of proper selection of parameters values occurring in the mathematical model.

Funding

Publication is partially supported under the Rector’s grant in the area of scientific research and development works, SILESIAN UNIVERSITY OF TECHNOLOGY, grant number 10/040/RGJ19/0082, and the statutory grant, FACULTY OF MECHANICAL ENGINEERING, SILESIAN UNIVERSITY OF TECHNOLOGY, grant number 10/990/BK_19/0139.

Conflicts of Interest

The author declares no conflict of interest.

References

  1. Bojakowska, U.; Kalinowski, P.; Kowalska, M.E. Epidemiology and prophylaxis of breast cancer. J. Edu. Health Sport 2016, 6, 701–710. [Google Scholar]
  2. Lokuhetty, D.; White, V.A.; Watanabe, R.; Cree, I.A.; International Agency for Research on Cancer (IARC). WHO Classification of Tumours. Breast Tumours, 5th ed.; WHO Classification of Tumours Editorial Board; World Health Organization: Genève, Switzerland, 2019. [Google Scholar]
  3. Wojciechowska, U.; Didkowska, J.; Morbidity and Deaths from Malignant Tumors in Poland. National Cancer Registry, Maria Skłodowska-Curie Institute of Oncology. Available online: http://onkologia.org.pl/rak-piersi-kobiet/ (accessed on 17 June 2019). (In Polish).
  4. Didkowska, J.; Wojciechowska, U. Breast cancer in Poland and Europe—Population and statistics. J. Oncol. 2013, 63, 111–118. [Google Scholar]
  5. Ferlay, J.; Shin, H.; Bray, F.; Forman, D.; Mathers, C.; Parkin, D.M. Estimates of worldwide burden of cancer in 2008: GLOBOCAN 2008. Int. J. Cancer 2010, 127, 2893–2917. [Google Scholar] [CrossRef] [PubMed]
  6. Chu, K.F.; Dupuy, D.E. Thermal ablation of tumours: Biological mechanisms and advances in therapy. Nat. Rev. Cancer 2014, 14, 199–208. [Google Scholar] [CrossRef]
  7. Zhao, Z.; Wu, F. Minimally-invasive thermal ablation of early-stage breast cancer: A systemic review. EJSO-Eur. J. Surg. Oncol. 2010, 36, 1149–1155. [Google Scholar] [CrossRef]
  8. Sundeep, S.; Repaka, R. Parametric sensitivity analysis of critical factors affecting the thermal damage during RFA of breast tumor. Int. J. Therm. Sci. 2018, 124, 366–374. [Google Scholar]
  9. Rzepka, E.; Püsküllüoglu, M. The role of hyperthermia in oncological treatment. Oncol. Clin. Pract. 2012, 8, 178–188. [Google Scholar]
  10. Pennes, H.H. Analysis of tissue and arterial blood temperature in the resting human forearm. J. Appl. Physiol. 1948, 1, 93–122. [Google Scholar] [CrossRef]
  11. Jamil, M.; Ng, E.Y.K. To optimize the efficacy of bioheat transfer in capacitive hyperthermia: A physical perspective. J. Therm. Biol. 2013, 38, 272–279. [Google Scholar] [CrossRef]
  12. Lv, Y.G.; Deng, Z.S.; Liu, J. 3D numerical study on the induced heating effects of embedded micro/nanoparticles on human body subject to external medical electromagnetic field. IEEE Trans. Nanobiosci. 2005, 4, 284–294. [Google Scholar]
  13. Schutt, D.; Berjano, E.J.; Haemmerich, D. Effect of electrode thermal conductivity in cardiac radiofrequency catheter ablation: A computational modeling study. Int. J. Hyperther. 2009, 25, 99–107. [Google Scholar] [CrossRef]
  14. Majchrzak, E.; Mochnacki, B.; Dziewonski, M.; Jasinski, M. Numerical modelling of hyperthermia and hypothermia processes. In Advanced Materials Research, Proceedings of the International Conference on Computational Materials Science (CMS 2011), Guangzhou, China, 17–18 April 2011; Xiong, F., Ed.; Trans Tech Publications: Guangzhou, China, 2011; pp. 268–270. [Google Scholar]
  15. Gasselhuber, A.; Dreher, M.R.; Negussie, A.; Wood, B.J.; Rattay, F.; Haemmerich, D. Mathematical spatio-temporal model of drug delivery from low temperature sensitive liposomes during radiofrequency tumour ablation. Int. J. Hyperther. 2010, 26, 499–513. [Google Scholar] [CrossRef] [Green Version]
  16. Paruch, M. Cancer ablation during RF hyperthermia using internal electrode. In Advances in Mechanics: Theoretical, Computational and Interdisciplinary Issues, 1st ed.; Kleiber, M., Burczynski, T., Wilde, K., Gorski, J., Winkelmann, K., Smakosz, L., Eds.; Taylor&Francis Group, CRC Press: Boca Raton, FL, USA, 2016; pp. 455–458. [Google Scholar]
  17. Andreuccetti, D.; Zoppetti, N. Quasi-static electromagnetic dosimetry: From basic principles to examples of applications. Int. J. Occup. Saf. Ergon. 2006, 12, 201–215. [Google Scholar] [CrossRef] [Green Version]
  18. Varon, L.A.B.; Orlande, H.R.B.; Elicabe, G.E. Estimation of state variables in the hyperthermia therapy of cancer with heating imposed by radiofrequency electromagnetic waves. Int. J. Therm. Sci. 2015, 98, 228–236. [Google Scholar] [CrossRef]
  19. Lamien, B.; Varon, L.A.B.; Orlande, H.R.B.; Elicabe, G.E. State estimation in bioheat transfer: A comparison of particle filter algorithms. Int. J. Numer. Methods Heat Fluid Flow 2017, 27, 615–638. [Google Scholar] [CrossRef]
  20. Ng, E.Y.K.; Jamil, M. Parametric sensitivity analysis of radiofrequency ablation with efficient experimental design. Int. J. Therm. Sci. 2014, 80, 41–47. [Google Scholar] [CrossRef]
  21. Paruch, M. Identification of the cancer ablation parameters during RF hyperthermia using gradient, evolutionary and hybrid algorithms. Int. J. Numer. Methods Heat Fluid Flow 2017, 27, 674–697. [Google Scholar] [CrossRef]
  22. Doss, J.D. Calculation of electric fields in conductive media. Med. Phys. 1982, 9, 566–573. [Google Scholar] [CrossRef]
  23. Van Sonnenberg, E.; McNullen, W.N.; Solbiati, L. Tumor Ablation. Principles and Practice; Springer: New York, NY, USA, 2005; pp. 3–16. [Google Scholar]
  24. Goldberg, S.N.; Saldinger, P.F.; Gazelle, G.S.; Huertas, J.C.; Stuart, K.E.; Jacobs, T.; Kruskal, J.B. Percutaneous tumor ablation: Increased necrosis with combined radio-frequency ablation and intratumoral doxorubicin injection in a rat breast tumor model. Radiology 2001, 220, 420–427. [Google Scholar] [CrossRef]
  25. Hung, K.Y.; Lin, Y.C.; Feng, H.P. The effects of acid etching on the nanomorphological surface characteristics and activation energy of titanium medical materials. Materials 2017, 10, 1164–1177. [Google Scholar] [CrossRef]
  26. Yiming, L.; Hongchao, J.; Zhongman, C.; Xuefeng, T.; Yaogang, L.; Jinping, L. Comparative study on constitutive models for 21-4N heat resistant steel during high temperature deformation. Materials 2019, 12, 1893–1915. [Google Scholar]
  27. Siteman Cancer Center. Breast Cancer Treatment. General Information about Breast Cancer. Available online: https://siteman.wustl.edu/ncipdq/cdr0000062955/ (accessed on 16 June 2019).
  28. Siteman Cancer Center. Breast Cancer Treatment. Stages of Breast Cancer. Available online: https://siteman.wustl.edu/ncipdq/cdr0000062955/ (accessed on 16 June 2019).
  29. Ng, E.Y.K.; Sudharsan, N.M. An improved three-dimensional direct numerical modelling and thermal analysis of a female breast with tumour. J. Eng. Med. 2001, 215, 25–37. [Google Scholar] [CrossRef]
  30. Wahab, A.A.; Salim, M.I.M.; Ahamat, M.A.; Manaf, N.A.; Yunus, J.; Lai, K.W. Thermal distribution analysis of three-dimensional tumor-embedded breast models with different breast density compositions. Med. Biol. Eng. Comput. 2016, 54, 1363–1373. [Google Scholar] [CrossRef]
  31. Cheng, D.K. Fundamentals of Engineering Electromagnetics; Addison-Wesley Publishing Company: Boston, MA, USA, 1993. [Google Scholar]
  32. Paruch, M.; Turchan, L. Mathematical modelling of the destruction degree of cancer under the influence of a RF hyperthermia. AIP Conf. Proc. 2018, 1922, 1–10. [Google Scholar]
  33. Paruch, M. Identification of the degree of tumor destruction on the basis of the Arrhenius integral using the evolutionary algorithm. Int. J. Therm. Sci. 2018, 130, 507–517. [Google Scholar] [CrossRef]
  34. Jasinski, M. Numerical analysis of soft tissue damage process caused by laser action. AIP Conf. Proc. 2018, 1922, 1–10. [Google Scholar]
  35. Gas, P. Optimization of multi-slot coaxial antennas for microwave thermotherapy based on the S11-parameter analysis. Biocybern. Biomed. Eng. 2017, 37, 78–93. [Google Scholar] [CrossRef]
  36. Zhang, B.; Moser, M.A.; Zhang, E.M.; Luo, Y.; Zhang, H.; Zhang, W. Study of the relationship between the target tissue necrosis volume and the target tissue size in liver tumours using two-compartment finite element RFA modelling. Int. J. Hyperth. 2014, 30, 593–602. [Google Scholar] [CrossRef]
  37. Nilsson, A.L. Blood flow, temperature, and heat loss of skin exposed to local radiative and convective cooling. J. Investg. Dermatol. 1987, 88, 586–593. [Google Scholar] [CrossRef] [Green Version]
  38. Osman, M.M.; Afify, E.M. Thermal modelling of the normal woman’s breast. J. Biomech. Eng. 1984, 106, 123–130. [Google Scholar] [CrossRef]
  39. Abraham, J.P.; Sparrow, E.M. A thermal-ablation bioheat model including liquid-to-vapor phase change, pressure- and necrosis-dependent perfusion, and moisture-dependent properties. Int. J. Heat Mass Transf. 2007, 50, 2537–2544. [Google Scholar] [CrossRef]
  40. Henriques, F.C. Studies of thermal injuries, V. The predictability and the significance of thermally induced rate process leading to irreversible epidermal injury. Arch. Pathol. 1947, 43, 489–502. [Google Scholar]
  41. Korczak, A.; Jasinski, M. Modelling of biological tissue damage process with application of interval arithmetic. J. Theor. Appl. Mech. 2019, 57, 249–261. [Google Scholar] [CrossRef]
  42. Gas, P.; Wyszkowska, J. Influence of multi-tine electrode configuration in realistic hepatic RF ablative heating. Arch. Electr. Eng. 2019, 68, 521–533. [Google Scholar]
  43. Singh, S.; Repaka, R. Numerical study to establish relationship between coagulation volume and target tip temperature during temperature-controlled radiofrequency ablation. Electromagn. Biol. Med. 2018, 37, 13–22. [Google Scholar] [CrossRef] [PubMed]
  44. Majchrzak, E.; Turchan, L.; Jasinski, M. Identification of laser intensity assuring the destruction of target region of biological tissue using the gradient method and generalized dual-phase lag equation. Iran. J. Sci. Technol. Trans. Mech. Eng. 2019, 43, 539–548. [Google Scholar] [CrossRef] [Green Version]
  45. Chang, I.A.; Nguyen, U.D. Thermal modeling of lesion growth with radiofrequency ablation devices. Biomed. Eng. Online 2004, 3, 1–19. [Google Scholar] [CrossRef] [Green Version]
  46. Jasinski, M. Mathematical Modelling of Tissue Damage Process Caused by External Heat Sources; Silesian University of Technology: Gliwice, Poland, 2016; pp. 23–25. (In Polish) [Google Scholar]
  47. Majchrzak, E.; Paruch, M.; Dziewonski, M.; Freus, S.; Freus, K. Sensitivity analysis of temperature field and parameter identification in burned and healthy skin tissue. In Book Series: Advanced Structured Materials: Computational Modeling, Optimization and Manufacturing Simulation of Advanced Engineering Materials; Munoz-Rojas, P.A., Ed.; Springer International Publishing: Cham, Switzerland, 2016; Volume 49, pp. 89–112. [Google Scholar]
  48. Kałuża, G.; Majchrzak, E.; Turchan, Ł. Sensitivity analysis of temperature field in the heated soft tissue with respect to the perturbations of porosity. Appl. Math. Model. 2017, 49, 498–513. [Google Scholar] [CrossRef]
  49. Paruch, M. Sensitivity analysis and the inverse problem in the mathematical modelling of tumor ablation using the interstitial hyperthermia. In Proceedings of the 4th Polish Congress of Mechanics and 23rd International Conference on Computer Methods in Mechanics, Cracow, Poland, 8–12 September 2019. [Google Scholar]
  50. Jain, R.K.; Ward-Hartley, K. Tumor blood flow-characterization, modifications and role in hyperthermia. IEEE Trans. Sonics Ultrason. 1984, 31, 504–526. [Google Scholar] [CrossRef]
Figure 1. Anatomy of the female breast: (a) healthy breast; (b) breast with cancer. Reproduced by permission from Terese Winslow; Adapted from [27,28].
Figure 1. Anatomy of the female breast: (a) healthy breast; (b) breast with cancer. Reproduced by permission from Terese Winslow; Adapted from [27,28].
Materials 13 00136 g001aMaterials 13 00136 g001b
Figure 2. Cross-sectional view of the developed breast models (a) extremely dense (ED); (b) heterogeneously dense (HD); (c) scattered fibro-glandular (SFG); (d) predominantly fatty (PF).
Figure 2. Cross-sectional view of the developed breast models (a) extremely dense (ED); (b) heterogeneously dense (HD); (c) scattered fibro-glandular (SFG); (d) predominantly fatty (PF).
Materials 13 00136 g002
Figure 3. (a) Dimensions of the domain considered, position of tumor center and location of the internal applicator; (b) fragment of the analyzed geometry with the marked control points.
Figure 3. (a) Dimensions of the domain considered, position of tumor center and location of the internal applicator; (b) fragment of the analyzed geometry with the marked control points.
Materials 13 00136 g003
Figure 4. Electric potential distribution (U = 10V) (numerical calculations).
Figure 4. Electric potential distribution (U = 10V) (numerical calculations).
Materials 13 00136 g004
Figure 5. Temperature distribution after 10 min (U = 10 V) (numerical calculations).
Figure 5. Temperature distribution after 10 min (U = 10 V) (numerical calculations).
Materials 13 00136 g005
Figure 6. Temperature history at the control points in tumor (U = 10 V) (numerical calculations).
Figure 6. Temperature history at the control points in tumor (U = 10 V) (numerical calculations).
Materials 13 00136 g006
Figure 7. Arrhenius integral courses at the control points in tumor (U = 10 V).
Figure 7. Arrhenius integral courses at the control points in tumor (U = 10 V).
Materials 13 00136 g007
Figure 8. Temperature history at the control point P4 for different values of electric potential U (numerical calculations).
Figure 8. Temperature history at the control point P4 for different values of electric potential U (numerical calculations).
Materials 13 00136 g008
Figure 9. Arrhenius integral courses at the control point P4 for different values of electric potential U.
Figure 9. Arrhenius integral courses at the control point P4 for different values of electric potential U.
Materials 13 00136 g009
Figure 10. Arrhenius integral courses at the control points in tumors for U = 12.5 V.
Figure 10. Arrhenius integral courses at the control points in tumors for U = 12.5 V.
Materials 13 00136 g010
Figure 11. Temperature history at the control points (U = 12.5 V) (numerical calculations).
Figure 11. Temperature history at the control points (U = 12.5 V) (numerical calculations).
Materials 13 00136 g011
Figure 12. Arrhenius integral courses at the control points for U = 12.5 V.
Figure 12. Arrhenius integral courses at the control points for U = 12.5 V.
Materials 13 00136 g012
Figure 13. Temperature history at the control point P4 for different values of tumor perfusion coefficient.
Figure 13. Temperature history at the control point P4 for different values of tumor perfusion coefficient.
Materials 13 00136 g013
Figure 14. Arrhenius integral courses at the control point P4 for different values of tumor perfusion coefficient.
Figure 14. Arrhenius integral courses at the control point P4 for different values of tumor perfusion coefficient.
Materials 13 00136 g014
Table 1. Electrical and thermo-physical properties.
Table 1. Electrical and thermo-physical properties.
TissueElectrical Conductivity σ (S/m)Specific Heat c (J/(kg·K))Thermal Conductivity λ (W/(mK))Density ρ (kg/m3)Metabolic Heat Source Qmet (W/m3)Blood Perfusion GB (1/s)
Gland0.56329600.3310417000.0005
Fat0.025423480.219114000.0002
Muscle0.43934210.4910907000.0008
Tumor0.7937700.48105077920.0053
Blood
Electrode
Trocar
-
108
105
3617
840
1045
-
18
0.026
1050
6450
70
-
-
-
-
-
-
Table 2. The parameters A and ΔE for different tissue types.
Table 2. The parameters A and ΔE for different tissue types.
Tissue TypeA (1/s)ΔE (J/mol)
Breast1.18 × 10443.02 × 105
Liver7.39 × 10392.58 × 105
Skin1.80 × 10513.27 × 105
Tissue with the capillaries1.98 × 101066.67 × 105
Aorta5.60 × 10634.30 × 105
Table 3. Coordinates of control points.
Table 3. Coordinates of control points.
Control Pointx (m)y (m)
P1 (in tumor)0.0220.0655
P2 (in tumor)0.0320.0755
P3 (in tumor)0.0320.0555
P4 (in tumor)0.0320.0655
P5 (in tumor)0.0420.0655
P6 (in muscle)0.010.0655
P7 (in gland)0.040.05
P8 (in gland)0.050.0655
P9 (in gland)0.0320.08
P10 (on the skin)0.05750.086
P11 (in fat)0.05650.08
Table 4. Values of parameters used during the numerical simulation.
Table 4. Values of parameters used during the numerical simulation.
ParameterValue
Baseline physiological temperature TB
Body core temperature Tb
37 °C
37 °C
Initial temperature of the breast and tumor T1–4
Initial temperature of the electrode and trocar T5–6
37 °C
25 °C
Ambient temperature Tamb25 °C
Effective heat transfer coefficient αeff13.5 W/(m2K)
Time step Δt
Length of the active part of the applicator
Width of the applicator
1 s
10 mm
1.41 mm
Table 5. Arrhenius integral value and probability of tissue destruction.
Table 5. Arrhenius integral value and probability of tissue destruction.
Control PointArrT °CProbability of Destruction %
P11.738850.682.43
P23.084952.595.43
P30.631346.946.81
P44.972853.999.31
P50.539346.441.68
P60.353344.429.76
P70.164241.515.14
P80.110939.710.5
P90.915548.359.97
P100.02132.72.08
P110.039434.63.86

Share and Cite

MDPI and ACS Style

Paruch, M. Mathematical Modeling of Breast Tumor Destruction Using Fast Heating during Radiofrequency Ablation. Materials 2020, 13, 136. https://doi.org/10.3390/ma13010136

AMA Style

Paruch M. Mathematical Modeling of Breast Tumor Destruction Using Fast Heating during Radiofrequency Ablation. Materials. 2020; 13(1):136. https://doi.org/10.3390/ma13010136

Chicago/Turabian Style

Paruch, Marek. 2020. "Mathematical Modeling of Breast Tumor Destruction Using Fast Heating during Radiofrequency Ablation" Materials 13, no. 1: 136. https://doi.org/10.3390/ma13010136

APA Style

Paruch, M. (2020). Mathematical Modeling of Breast Tumor Destruction Using Fast Heating during Radiofrequency Ablation. Materials, 13(1), 136. https://doi.org/10.3390/ma13010136

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