Mathematical Modeling of Targeted Drug Delivery Using Magnetic Nanoparticles during Intraperitoneal Chemotherapy

Intraperitoneal (IP) chemotherapy has emerged as a promising method for the treatment of peritoneal malignancies (PMs). However, microenvironmental barriers in the tumor limit the delivery of drug particles and their deep penetration into the tumor, leading to reduced efficiency of treatment. Therefore, new drug delivery systems should be developed to overcome these microenvironmental barriers. One promising technique is magnetically controlled drug targeting (MCDT) in which an external magnetic field is utilized to concentrate drug-coated magnetic nanoparticles (MNPs) to the desired area. In this work, a mathematical model is developed to investigate the efficacy of MCDT in IP chemotherapy. In this model, considering the mechanism of drug binding and internalization into cancer cells, the efficacy of drug delivery using MNPs is evaluated and compared with conventional IP chemotherapy. The results indicate that over 60 min of treatment with MNPs, drug penetration depth increased more than 13 times compared to conventional IPC. Moreover, the drug penetration area (DPA) increased more than 1.4 times compared to the conventional IP injection. The fraction of killed cells in the tumor in magnetic drug delivery was 6.5%, which shows an increase of more than 2.5 times compared to that of the conventional method (2.54%). Furthermore, the effects of magnetic strength, the distance of the magnet to the tumor, and the magnetic nanoparticles’ size were evaluated. The results show that MDT can be used as an effective technique to increase the efficiency of IP chemotherapy.


Introduction
Cancers that occur in organs associated with the peritoneal cavity are prone to metastasize to the intraperitoneal space. Patients with these malignancies experience low quality of life due to the fact of complications such as urinary blockage, ascites, pain, and inability to eat and drink. The five-year survival rate is 49% in these types of cancer. The average survival rate in peritoneal carcinomatosis of gastric origin is 1-3 months [1]. Management of peritoneal metastases (PMs) has always been a challenging issue. Chemotherapy with systemic injections in the 1980s had a palliative approach and predicted a survival time of less than a few months. With the introduction of peritoneal plasma membrane as an opportunity for the regional treatment of PMs by IP injection of chemotherapy drugs,

Materials and Methods
In clinical IP chemotherapy, chemotherapeutic agents are transferred to the desired location within a 1-2 h cycle. A schematic of drug delivery by IP injection of MNPs is shown in Figure 1. Drug-coated MNPs are injected into the peritoneal cavity and guided by a permanent magnet. In the following, the governing equations of a magnetically controlled IP drug targeting system are presented.

Interstitial Fluid Flow
Considering the intercapillary distance (33-98 μm), which is usually 2-3 orders of magnitude smaller than the length scale of drug transfer [26,27], tumor tissue acts as a porous media. Hence, Darcy's law in porous media was used to explain the interstitial fluid flow [28,29]: Here, κ is the interstitium's hydraulic conductivity. In addition, Pi and are the interstitial fluid pressure and velocity, respectively. The steady-state mass conservation equation for the incompressible interstitial fluid is as follows: In this equation, is the net flow rate from blood vessels into the interstitium per unit volume, and is the net lymphatic drainage per unit volume. Starling's law was used to calculate and according to the equation below: where LP and indicate the hydraulic conductivity of the microvascular wall and vascular surface area per unit volume, respectively, and and , are the intravascular pressure and interstitial pressure, respectively.
and πB indicate the average osmotic reflection coefficient of plasma proteins and plasma osmotic pressure. Likewise, πi is the interstitial fluid osmotic pressure. Lymphatic drainage ( ) is related to the difference between interstitial fluid pressure and lymphatic pressure: where LPL is the hydraulic conductivity of the lymphatic vessel wall, is the surface area per unit volume of lymphatic vessels, and is the pressure in the lymphatic vessels. Here, is considered zero because of the lack of an effective lymphatic system in solid tumors.

Interstitial Fluid Flow
Considering the intercapillary distance (33-98 µm), which is usually 2-3 orders of magnitude smaller than the length scale of drug transfer [26,27], tumor tissue acts as a porous media. Hence, Darcy's law in porous media was used to explain the interstitial fluid flow [28,29]: Here, κ is the interstitium's hydraulic conductivity. In addition, P i and v i are the interstitial fluid pressure and velocity, respectively. The steady-state mass conservation equation for the incompressible interstitial fluid is as follows: In this equation, φ B is the net flow rate from blood vessels into the interstitium per unit volume, and φ L is the net lymphatic drainage per unit volume. Starling's law was used to calculate φ B and φ L according to the equation below: where L P and S V indicate the hydraulic conductivity of the microvascular wall and vascular surface area per unit volume, respectively, and P B and P i , are the intravascular pressure and interstitial pressure, respectively. σ s and π B indicate the average osmotic reflection coefficient of plasma proteins and plasma osmotic pressure. Likewise, π i is the interstitial fluid osmotic pressure. Lymphatic drainage (φ L ) is related to the difference between interstitial fluid pressure and lymphatic pressure: where L PL is the hydraulic conductivity of the lymphatic vessel wall, S L V is the surface area per unit volume of lymphatic vessels, and P L is the pressure in the lymphatic vessels. Here, φ L is considered zero because of the lack of an effective lymphatic system in solid tumors.

Mass Transport
Maxwell-Ampere's law was used to correlate the magnetic field and the current intensity based on the magneto-static nature of the problem (Equation (5)). Moreover, Gauss's law (Equation (6)) was used to describe the magnetic flux density [18,30]: Here, H is the magnetic field, J is the current intensity, and B is the magnetic flux density. As we are using a permanent magnet, the current density (J) was considered 0. For the air and tissue domain, the constitutive equation B = µ 0 H was used, while for the magnet domain, the equation B = µ 0 µ r H + B rem was applied [31]. In these equations, the vacuum magnetic permeability was considered µ 0 = 4 π × 10 −7 NA −2 , the relative magnetic permeability µ r = 1000, and the remnant magnetic flux was defined as B rem . In a magnetic field (H), the magnetic force exerted on an MNP is determined as [32]: Here, V MNP = 4 3 πa 3 and χ are the volume and magnetic susceptibility of MNPs, respectively. In the case of a strong magnetic field, which can cause particle saturation [33], M sat indicates the saturated magnetization of the particle [33][34][35]. For well-made particles, → M sat is on the order of 0.5 T and lines up with → H [36,37]. The saturation of the particle does not affect the direction of the force, only its magnitude.
The free drug's concentration in the interstitial fluid is calculated using the convectiondiffusion equation as follows: where C F is the free drug concentration, D F is the diffusion coefficient of free drug in the porous medium, and C rec and ϕ are the concentration of cell surface receptors and the tumor volume fraction available for the drugs, respectively. K ON is the constants of the drug association rate to the cancer cells. In addition, v and D mnp are the velocity and effective diffusion coefficient of MNPs in tissue. The velocity (v) was calculated as the sum of the local velocity of interstitial fluid (v i ) and the equilibrium velocity (v e ). v e corresponds to when the Stokes drag force, Fs = (6πaη)ve, is equal to the magnetic force (Fm) [15,18]: Here, η = 1.12 × 10 −3 Pa.s is the dynamic viscosity of the interstitial fluid. Φ in Equation (8), is the net rate of free drug exchange from blood and lymphatic vessels and is calculated as follows: Using the pore model [38,39] for trans-capillary transfer, Φ B and Φ L are expressed as follows: where P is the permeability of tumor microvessels, Cp is the concentration of MNPs in the blood plasma, and σ f is the osmotic reflection coefficient for the MNPs. Pe is the Peclet number defined as follows: The equation for the bound drug concentration is as follows: where C B is the drug concentration bound to the cancer cell, and K INT is the constant rate of drug internalization into the cancer cells. Moreover, the equation internalized drug concentration is expressed as: where C I is the concentration of the drug that is internalized into cancer cells.

Cell Survival Model
The fraction of killed cells was calculated as 1 − S F , where S F is the fraction of surviving cells. S F was obtained from Equation (16) [40], which is based on an in vitro study by Ker et al. [41] for the fraction of non-small surviving lung tumor cells.
In Equation (16), ω = 0.4938 and C I indicates the internalized drug. The values of the model parameters, including tissue and solute transport parameters, are listed in Tables 1 and 2, respectively.

Model Geometry and Boundary Conditions
A solid tumor may have a necrotic area that lacks functional blood or lymph vessels; thus, there is no exchange of fluid with the interstitium in this region. In contrast, the regions of the tumor that are outside of the necrotic core contain rapidly dividing cells and blood vessels. In order to consider this non-uniform perfusion rate in the tumor, a biologically related structure, including an area with a necrotic core and a surrounding with leaky vasculature, is defined as a tumor [49]. A representation of the model's geometry is presented in Figure 2. A magnet with length, l, and width, h, lies at a distance, d, from the tumor. A circular tumor with radius, R, and a necrotic core with radius, R n , were considered. Doxorubicin has been used as one of the most widely used drugs in the model.

Model Geometry and Boundary Conditions
A solid tumor may have a necrotic area that lacks functional blood or lymph vessels; thus, there is no exchange of fluid with the interstitium in this region. In contrast, the regions of the tumor that are outside of the necrotic core contain rapidly dividing cells and blood vessels. In order to consider this non-uniform perfusion rate in the tumor, a biologically related structure, including an area with a necrotic core and a surrounding with leaky vasculature, is defined as a tumor [49]. A representation of the model's geometry is presented in Figure 2. A magnet with length, l, and width, h, lies at a distance, d, from the tumor. A circular tumor with radius, R, and a necrotic core with radius, Rn, were considered. Doxorubicin has been used as one of the most widely used drugs in the model. The boundary conditions for the present study are listed in Table 3. For internal boundaries in the tumor, the continuity of IFP, concentration, and its flux were considered, where Ω and Ω demonstrate the necrotic and viable areas of the tumor at the boundary. Moreover, for the outer boundary, the Dirichlet boundary condition was applied.

Solution Strategy
To create the initial conditions for time-dependent numerical simulations, fluid flow equations were primarily solved to achieve a steady-state solution in the computational domain. Then, the values of velocity and pressure obtained at time zero were used to simulate drug delivery. The solution strategy for this section is also shown in Figure 3.  Table 3. For internal boundaries in the tumor, the continuity of IFP, concentration, and its flux were considered, where Ω − and Ω + demonstrate the necrotic and viable areas of the tumor at the boundary. Moreover, for the outer boundary, the Dirichlet boundary condition was applied.
The inner boundary of the tumor The outer boundary of the tumor

Solution Strategy
To create the initial conditions for time-dependent numerical simulations, fluid flow equations were primarily solved to achieve a steady-state solution in the computational domain. Then, the values of velocity and pressure obtained at time zero were used to simulate drug delivery. The solution strategy for this section is also shown in Figure 3.

Conventional IPC
In this section, we examine the results of drug transfer by conventional IP chemotherapy. The spatiotemporal distribution of the free drug (CF), bound drug (CB), and drugs internalized into cancer cells (CI) were examined in the tumor during one hour of chemotherapy. Two main criteria were considered for the evaluation of the performance of chemotherapy with this method: 1. Fraction of killed cells (FK) was utilized as the major parameter for quantitatively evaluating drug delivery efficiency. Since drug particles do not reach all tumor parts in chemotherapy, this fraction was first studied as the fraction of killed cells in the drug penetration area (FKPA). Then, in order to evaluate the treatment of the whole tumor, this parameter was calculated and examined as the effective fraction of killed cells (FKeff), which represents the fraction of killed cells by considering the whole tumor; 2. To evaluate the performance of the drug delivery system in improving the penetration of the drug into the tumor, W1/2 is defined as the distance from the outer border of the tumor where the concentration of a free drug is equal to 50% of the concentration at the tumor border [50]. This parameter considers the slope of decreasing the drug concentration by moving from the tumor border to its depth. The larger the value of w1/2, the lower the concentration drop when moving to the depth of the tumor and the greater the penetration of the drug into the tumor.

Distribution of IFV and IFP in the Tumor
Tumor microenvironment features have a crucial role in the way a drug is delivered to the tumor. Higher cell density in tumors results in reducing the permeability of tumor tissue compared to the normal one. Figures 4 and 5 show the distribution of IFP and IFV in the tumor, respectively. As indicated in Figure 4a,b, IFP was highest at the tumor center

Conventional IPC
In this section, we examine the results of drug transfer by conventional IP chemotherapy. The spatiotemporal distribution of the free drug (C F ), bound drug (C B ), and drugs internalized into cancer cells (C I ) were examined in the tumor during one hour of chemotherapy. Two main criteria were considered for the evaluation of the performance of chemotherapy with this method:

1.
Fraction of killed cells (FK) was utilized as the major parameter for quantitatively evaluating drug delivery efficiency. Since drug particles do not reach all tumor parts in chemotherapy, this fraction was first studied as the fraction of killed cells in the drug penetration area (FK PA ). Then, in order to evaluate the treatment of the whole tumor, this parameter was calculated and examined as the effective fraction of killed cells (FK eff ), which represents the fraction of killed cells by considering the whole tumor; 2.
To evaluate the performance of the drug delivery system in improving the penetration of the drug into the tumor, w 1/2 is defined as the distance from the outer border of the tumor where the concentration of a free drug is equal to 50% of the concentration at the tumor border [50]. This parameter considers the slope of decreasing the drug concentration by moving from the tumor border to its depth. The larger the value of w 1/2 , the lower the concentration drop when moving to the depth of the tumor and the greater the penetration of the drug into the tumor.

Distribution of IFV and IFP in the Tumor
Tumor microenvironment features have a crucial role in the way a drug is delivered to the tumor. Higher cell density in tumors results in reducing the permeability of tumor tissue compared to the normal one. Figures 4 and 5 show the distribution of IFP and IFV in the tumor, respectively. As indicated in Figure 4a,b, IFP was highest at the tumor center (1533 Pa), except in a small area near the outer boundary of the tumor where the pressure decreased significantly. According to Darcy's law (Equation (1)), since the pressure gradient is zero in a large part of the tumor, IFV had an insignificant value in this area (Figure 5a). Similarly, in the outer boundary of the tumor, because of the high amount of pressure gradient, the IFP increased sharply and reached its maximum level (0.17 µm/s) (Figure 5b). (1533 Pa), except in a small area near the outer boundary of the tumor where the pressure decreased significantly. According to Darcy's law (Equation (1)), since the pressure gradient is zero in a large part of the tumor, IFV had an insignificant value in this area ( Figure  5a). Similarly, in the outer boundary of the tumor, because of the high amount of pressure gradient, the IFP increased sharply and reached its maximum level (0.17 μm/s) ( Figure  5b).
(a) (b)  Figure 6 shows the profiles of the average concentration of doxorubicin in the forms of free, bound, and internalized over 60 min of injection time. As is clear, immediately after injection, the concentration of free doxorubicin in the tumor skyrocketed to its maximum level (0.0013 mol/m 3 ) and then plateaued. The same applied to the concentration of bound drug, but the difference was that the concentration of bound drug reached its maximum more gradually. Unlike free and bound drug concentrations, the concentration of internalized drug was constantly increasing. This increasing internalized drug concentration cell as well as the non-decreasing trend of the free and bound drugs were due to the constant concentration of the drug at the outer boundary of the tumor within one hour of administration. In other words, with the continuous injection of free drug, the concentration of free and bound drug reached a steady-state during treatment. Obviously, with the (1533 Pa), except in a small area near the outer boundary of the tumor where the pressure decreased significantly. According to Darcy's law (Equation (1)), since the pressure gradient is zero in a large part of the tumor, IFV had an insignificant value in this area ( Figure  5a). Similarly, in the outer boundary of the tumor, because of the high amount of pressure gradient, the IFP increased sharply and reached its maximum level (0.17 μm/s) ( Figure  5b).

Concentration Distribution in the Tumor
(a) (b)  Figure 6 shows the profiles of the average concentration of doxorubicin in the forms of free, bound, and internalized over 60 min of injection time. As is clear, immediately after injection, the concentration of free doxorubicin in the tumor skyrocketed to its maximum level (0.0013 mol/m 3 ) and then plateaued. The same applied to the concentration of bound drug, but the difference was that the concentration of bound drug reached its maximum more gradually. Unlike free and bound drug concentrations, the concentration of internalized drug was constantly increasing. This increasing internalized drug concentration cell as well as the non-decreasing trend of the free and bound drugs were due to the constant concentration of the drug at the outer boundary of the tumor within one hour of administration. In other words, with the continuous injection of free drug, the concentration of free and bound drug reached a steady-state during treatment. Obviously, with the  Figure 6 shows the profiles of the average concentration of doxorubicin in the forms of free, bound, and internalized over 60 min of injection time. As is clear, immediately after injection, the concentration of free doxorubicin in the tumor skyrocketed to its maximum level (0.0013 mol/m 3 ) and then plateaued. The same applied to the concentration of bound drug, but the difference was that the concentration of bound drug reached its maximum more gradually. Unlike free and bound drug concentrations, the concentration of internalized drug was constantly increasing. This increasing internalized drug concentration cell as well as the non-decreasing trend of the free and bound drugs were due to the constant concentration of the drug at the outer boundary of the tumor within one hour of administration. In other words, with the continuous injection of free drug, the concentration of free and bound drug reached a steady-state during treatment. Obviously, with the end of the drug injection into the peritoneal cavity, the mentioned concentrations will decrease. However, the main focus in this study was a one-hour time span in which a continuous concentration of the drug was injected into the peritoneal cavity.

Concentration Distribution in the Tumor
Pharmaceutics 2022, 14, x FOR PEER REVIEW 9 of 24 end of the drug injection into the peritoneal cavity, the mentioned concentrations will decrease. However, the main focus in this study was a one-hour time span in which a continuous concentration of the drug was injected into the peritoneal cavity.  Figure 7 show the contours of the concentration of free, bound, and internalized drugs with the IP injection of doxorubicin at one hour after injection. As seen in Figure 7, the penetration of the drug into the tumor by IP chemotherapy was curbed to a limited area of the outer boundary of the tumor and, as a result, a large part of the tumor was out of reach for drug delivery. This can be examined more precisely in Figure 8a, where the mean drug concentration profiles in the tumor for the free, bound, and internalized forms along the tumor radius after 60 min of treatment are presented. As shown in this figure, the concentration decreased rapidly to zero as they moved away from the tumor boundary. Since the internalized drug was supplied by the drug that was bound to the cancer cells and the bound drug was also a part of the free drug available in the tumor, the concentration of the internalized drug was therefore always much lower than the free drug concentration available in the tumor. In order to acquire a higher concentration of the drug internalized into the cancer cells, more extracellular drugs should be available. Moreover, the drug must penetrate the tumor at the proper depth so that the tumor is homogeneously exposed to the drug. For this reason, the W1/2 was examined to evaluate the drug's penetration into the tumor during treatment. The value obtained for W1/2 in this section was equal to 60 μm and, given that the tumor's radius was 10 mm, the relative penetration depth of the tumor W1/2% was obtained as 0.6%. In addition, calculating the relative drug penetration area (PArel) re-  Figure 7 show the contours of the concentration of free, bound, and internalized drugs with the IP injection of doxorubicin at one hour after injection. As seen in Figure 7, the penetration of the drug into the tumor by IP chemotherapy was curbed to a limited area of the outer boundary of the tumor and, as a result, a large part of the tumor was out of reach for drug delivery. This can be examined more precisely in Figure 8a, where the mean drug concentration profiles in the tumor for the free, bound, and internalized forms along the tumor radius after 60 min of treatment are presented. As shown in this figure, the concentration decreased rapidly to zero as they moved away from the tumor boundary. Since the internalized drug was supplied by the drug that was bound to the cancer cells and the bound drug was also a part of the free drug available in the tumor, the concentration of the internalized drug was therefore always much lower than the free drug concentration available in the tumor. end of the drug injection into the peritoneal cavity, the mentioned concentrations will decrease. However, the main focus in this study was a one-hour time span in which a continuous concentration of the drug was injected into the peritoneal cavity.  Figure 7 show the contours of the concentration of free, bound, and internalized drugs with the IP injection of doxorubicin at one hour after injection. As seen in Figure 7, the penetration of the drug into the tumor by IP chemotherapy was curbed to a limited area of the outer boundary of the tumor and, as a result, a large part of the tumor was out of reach for drug delivery. This can be examined more precisely in Figure 8a, where the mean drug concentration profiles in the tumor for the free, bound, and internalized forms along the tumor radius after 60 min of treatment are presented. As shown in this figure, the concentration decreased rapidly to zero as they moved away from the tumor boundary. Since the internalized drug was supplied by the drug that was bound to the cancer cells and the bound drug was also a part of the free drug available in the tumor, the concentration of the internalized drug was therefore always much lower than the free drug concentration available in the tumor. In order to acquire a higher concentration of the drug internalized into the cancer cells, more extracellular drugs should be available. Moreover, the drug must penetrate the tumor at the proper depth so that the tumor is homogeneously exposed to the drug. For this reason, the W1/2 was examined to evaluate the drug's penetration into the tumor during treatment. The value obtained for W1/2 in this section was equal to 60 μm and, given that the tumor's radius was 10 mm, the relative penetration depth of the tumor W1/2% was obtained as 0.6%. In addition, calculating the relative drug penetration area (PArel) re- In order to acquire a higher concentration of the drug internalized into the cancer cells, more extracellular drugs should be available. Moreover, the drug must penetrate the tumor at the proper depth so that the tumor is homogeneously exposed to the drug. For this reason, the w 1/2 was examined to evaluate the drug's penetration into the tumor during treatment. The value obtained for w 1/2 in this section was equal to 60 µm and, given that the tumor's radius was 10 mm, the relative penetration depth of the tumor w 1/2 % was obtained as 0.6%. In addition, calculating the relative drug penetration area (PA rel ) revealed that 87.3% of the tumor area remained untreated (Figure 8c). These outcomes accurately show one of the main challenges of the IPC, which is the very low penetration depth of the drug. The poor drug penetration depth was a result of the adverse pressure gradient at the tumor boundary followed by the outward convection flow, which exists at the tumor boundary as described in Figures 4 and 5. In order to evaluate the effectiveness of treatment during 60 min of chemotherapy by IP injection, the time profile of the fraction of killed cells (FK) in the drug penetration area is presented in Figure 8c. The FK PA values increased over 60 min of treatment and reached 28.3% at the end of treatment.
depth of the drug. The poor drug penetration depth was a result of the adverse pressure gradient at the tumor boundary followed by the outward convection flow, which exists at the tumor boundary as described in Figures 4 and 5. In order to evaluate the effectiveness of treatment during 60 min of chemotherapy by IP injection, the time profile of the fraction of killed cells (FK) in the drug penetration area is presented in Figure 8c. The FKPA values increased over 60 min of treatment and reached 28.3% at the end of treatment.  Table 4 summarized the calculated values for the treatment evaluation parameters, including drug penetration depth in tumor W1⁄2, the relative drug penetration area (PArel), fraction of killed cells in the area of drug penetration (FKPA), and the effective fraction of killed cells (FKeff) after 60 min of treatment. According to this table, PArel and FKeff are 8.97% and 2.54%, respectively.   Table 4 summarized the calculated values for the treatment evaluation parameters, including drug penetration depth in tumor w 1/2 , the relative drug penetration area (PA rel ), fraction of killed cells in the area of drug penetration (FK PA ), and the effective fraction of killed cells (FK eff ) after 60 min of treatment. According to this table, PA rel and FK eff are 8.97% and 2.54%, respectively. The results of this section generally indicate the low efficacy of the conventional IPC. Even though the IP injection has been conspicuous because of its potential in the regional treatment and direct transfer of drugs to the tumor, it is necessary to provide strategies that can improve treatment using IP injection. The results presented in this section provide a general insight into drug delivery by conventional IP injection and review it from the perspective of fluid and mass transfer, which can help design a treatment strategy that can overcome obstacles of drug delivery to the tumor.

Magnetically Controlled IP Drug Targeting
In this section, the results of the MNP delivery are presented. To evaluate the improvement of treatment using this method. The results of this section were compared to the results of conventional IPC. Firstly, a baseline model is presented in order to compare magnetic IP drug delivery with the conventional IPC. Next, the impact of the MNPs' size, the magnetic strength, and the distance of the magnet from the tumor were investigated.

Baseline Model
In the baseline model, a permanent magnet of 20 by 10 cm in dimension, which was placed at a distance of 5 cm from the tumor, was used to assist in drug delivery to the tumor. The radius of the MNPs used in the baseline model was 100 nm, and the magnetic strength was considered as 1.5 T. As larger tumors are more dangerous due to the possibility of disease recurrence after treatment [1,2], in this study a large tumor with a radius of 10 mm was selected.
In Figure 9, the mean concentrations of the free, bound, and internalized drugs are presented over 60 min of treatment. At the end of the treatment, C F , C B , and C I were calculated to be 0.0054, 0.0022, and 0.00025, respectively. Compared to the corresponding values for conventional IPC (i.e., 0.0012, 0.0005, and 0.009), an increase of 4.5, 4.4, and 2.77 times in the mean concentration of the free, bound, and internalized drug were achieved, respectively. Contours of C F , C B , and C I are shown in Figure 10a-c. Comparison of these contours with the contours of conventional IPC in Figure 7 shows that the drug distribution became much better in all of its three forms using MNPs. In traditional intraperitoneal drug delivery, the drug penetrates the outer border of the tumor uniformly from all radial directions. But in drug delivery using a magnet, the magnetic force determines the area where the drug enters the tumor. In the case we simulated here, the direction of the magnetic force was upward. Thus, in the lower area of the tumor, the magnetic force pushed the drug into the tumor tissue, while in the upper area of the tumor, it the opposite effect and impeded the drug from penetrating the tumor.
In Figure 11a, the C F , C B , and C I distributions along the vertical line passing through the center of the tumor are shown 60 min after the start of injection. The drug had a relatively good penetration into the tumor at the lower area of the tumor (r = 2 cm), the concentration of the drug decreased in all of its three forms as we moved to the center of the tumor. Figure 11b presents the fraction of killed cells in the area with high drug penetration (FK PA ) over time. FK PA gradually increased from zero at the starting point to 51.2% at the end of the treatment. For more accurate evaluation and comparing the efficacy of treatment using MNPs, w 1/2 , PA rel , FK PA , and FK eff were calculated and presented in Table 5. w 1/2 was equal to 0.08 cm using magnetic nanoparticles with a radius of 100 nm, which shows the significant effect of the proposed technique in increasing the penetration depth of the drug by more than 13 times. As mentioned above, the penetration of the drug in the magnetic nanoparticles' drug delivery occurred only in the lower area of the tumor. Moreover, the value of the drug penetration area in the magnetic drug delivery was equal to 0.40 cm 2 , which, compared to conventional IPC, showed an increase of more than 1.4 times. The value of FK eff in the tumor using the magnetic drug delivery was also calculated to be 6.5%, which shows a more than 2.5 times increase than the obtained value for the conventional method (2.54%). Therefore, considering all these criteria for the treatment's evaluation, it can be said that magnetically controlled IP drug targeting has improved drug delivery to the tumor compared to the conventional IPC using doxorubicin. The performance of the proposed method was also compared with the results of two other drugs including paclitaxel and cisplatin (Supplementary Materials Figure S1). The average concentration of paclitaxel and cisplatin 60 min after the treatment were 0.0013 and 0.0006 mol/m 3 , respectively, all of which were less than the amount obtained for magnetically controlled drug delivery. Therefore, using this method can also be useful for the transfer of other types of drugs. In the following, by changing the controllable parameters, such as the size of the MNPs, magnetic strength, and the magnet's distance from the tumor, drug delivery was studied to achieve higher treatment efficacy using this technique. times in the mean concentration of the free, bound, and internalized drug were achieved, respectively. Contours of CF, CB, and CI are shown in Figure 10a-c. Comparison of these contours with the contours of conventional IPC in Figure 7 shows that the drug distribution became much better in all of its three forms using MNPs. In traditional intraperitoneal drug delivery, the drug penetrates the outer border of the tumor uniformly from all radial directions. But in drug delivery using a magnet, the magnetic force determines the area where the drug enters the tumor. In the case we simulated here, the direction of the magnetic force was upward. Thus, in the lower area of the tumor, the magnetic force pushed the drug into the tumor tissue, while in the upper area of the tumor, it the opposite effect and impeded the drug from penetrating the tumor. In Figure 11a, the CF, CB, and CI distributions along the vertical line passing through the center of the tumor are shown 60 min after the start of injection. The drug had a relatively good penetration into the tumor at the lower area of the tumor (r = 2 cm), the concentration of the drug decreased in all of its three forms as we moved to the center of the tumor. Figure 11b presents the fraction of killed cells in the area with high drug penetration (FKPA) over time. FKPA gradually increased from zero at the starting point to 51.2% at the end of the treatment. For more accurate evaluation and comparing the efficacy of treatment using MNPs, w1/2, PArel, FKPA, and FKeff were calculated and presented in Table 5. w1/2 was equal to 0.08 cm using magnetic nanoparticles with a radius of 100 nm, which shows the significant effect of the proposed technique in increasing the penetration depth of the drug by more than 13 times. As mentioned above, the penetration of the drug in the magnetic nanoparticles' drug delivery occurred only in the lower area of the tumor. Moreover, the value of the drug penetration area in the magnetic drug delivery was equal to 0.40 cm 2 , which, compared to conventional IPC, showed an increase of more than 1.4 times. The value of FKeff in the tumor using the magnetic drug delivery was also calculated to be 6.5%, which shows a more than 2.5 times increase than the obtained value for the conventional method (2.54%). Therefore, considering all these criteria for the treatment's evaluation, it can be said that magnetically controlled IP drug targeting has improved drug delivery to the tumor compared to the conventional IPC using doxorubicin. The performance of the proposed method was also compared with the results of two other drugs including paclitaxel and cisplatin (Supplementary Materials Figure S1). The average concentration of paclitaxel and cisplatin 60 min after the treatment were 0.0013 and 0.0006 mol/m 3 , respectively, all of which were less than the amount obtained for magnetically controlled drug delivery. Therefore, using this method can also be useful for the transfer

Effect of the Magnetic Nanoparticles' Size
In this section, the impact of the size of the MNPs used in magnetically controlled IP chemotherapy was studied. Eight different nanoparticles with radii in the range of 25-600 nm were studied. In Figure 12a-h, contours of the free drug concentration distribution are

Effect of the Magnetic Nanoparticles' Size
In this section, the impact of the size of the MNPs used in magnetically controlled IP chemotherapy was studied. Eight different nanoparticles with radii in the range of 25-600 nm were studied. In Figure 12a-h, contours of the free drug concentration distribution are shown for radius ranges of 25, 50, 100, 200, 300, 400, 500, and 600 nm. These contours indicate that magnetic nanoparticles with radii of 25 and 50 nm have poor penetration into the tumor, while increasing the size of nanoparticles to 100 nm significantly increased the penetration of the drug from the lower part of the tumor. Increasing the radius of the nanoparticles to a value greater than 100 nm will dramatically improve the penetration of the drug into the tumor. tion.

Effect of the Magnetic Nanoparticles' Size
In this section, the impact of the size of the MNPs used in magnetically controlled IP chemotherapy was studied. Eight different nanoparticles with radii in the range of 25-600 nm were studied. In Figure 12a-h, contours of the free drug concentration distribution are shown for radius ranges of 25, 50, 100, 200, 300, 400, 500, and 600 nm. These contours indicate that magnetic nanoparticles with radii of 25 and 50 nm have poor penetration into the tumor, while increasing the size of nanoparticles to 100 nm significantly increased the penetration of the drug from the lower part of the tumor. Increasing the radius of the nanoparticles to a value greater than 100 nm will dramatically improve the penetration of the drug into the tumor. Figure 13 shows the mean concentration of free drugs in the tumor for different sizes of MNPs. Comparing the diagrams, we can see that by increasing the size of the MNPs with a radius in the range of 25-300 nm, CF increased. The change in the size from 300 to 400 nm and to more than 400 nm had different effects on CF. In this stage of the treatment, the increase in the size of MNPs from 300 to 400 nm resulted in a negligible decrease in the concentration; in contrast, increasing the nanoparticles' size to more than 400 nm led to a considerable reduction in the concentration of the free drug in the tumor. This was due to the competition of two factors in the penetration of nanoparticles: the magnetic force was greater for larger MNPs, while the microenvironmental barriers of the tumor that impeded the nanoparticles' transport in the tissue had much more of an effect on the larger particles. Thus, the competition between these two factors, namely, the magnetic force exerted on the particle and the microenvironmental barriers, determined the distribution of MNPs with different sizes.   Figure 13 shows the mean concentration of free drugs in the tumor for different sizes of MNPs. Comparing the diagrams, we can see that by increasing the size of the MNPs with a radius in the range of 25-300 nm, C F increased. The change in the size from 300 to 400 nm and to more than 400 nm had different effects on C F . In this stage of the treatment, the increase in the size of MNPs from 300 to 400 nm resulted in a negligible decrease in the concentration; in contrast, increasing the nanoparticles' size to more than 400 nm led to a considerable reduction in the concentration of the free drug in the tumor. This was due to the competition of two factors in the penetration of nanoparticles: the magnetic force was greater for larger MNPs, while the microenvironmental barriers of the tumor that impeded the nanoparticles' transport in the tissue had much more of an effect on the larger particles. Thus, the competition between these two factors, namely, the magnetic force exerted on the particle and the microenvironmental barriers, determined the distribution of MNPs with different sizes.  In Figure 14, CF distribution along the vertical line of the tumor is presented. Drug delivery using nanoparticles with a radius of 300 and 400 nm resulted in deeper penetration of the drug into the tumor and, subsequently, it increased CF in the tumor. However, 25 and 50 nm particles showed poor drug penetration into the tumor.
Efficacy parameters w1/2, PArel, FKPA, and FKeff are reported in Table 6. There was an increase in w1/2 values from 0.003 cm for 25 nm particles to 0.235 cm for 300 nm particles. As the nanoparticle size further increased, w1/2 decreased, until it reached 0.143 for 600 nm particles. The smallest PArel among the compared nanoparticles was related to 25 nm nanoparticles (0.20%), and 300 nm nanoparticles with 55.02% had the largest PArel among the studied MNPs. FKeff values for different sizes of MNPs after 60 min of treatment are presented in Figure 15. FKeff was calculated as the lowest value for 25 nm nanoparticles, which was equal to 0.11%, whereas it was the highest for 300 nm particles with 24.1%. In Figure 14, C F distribution along the vertical line of the tumor is presented. Drug delivery using nanoparticles with a radius of 300 and 400 nm resulted in deeper penetration of the drug into the tumor and, subsequently, it increased C F in the tumor. However, 25 and 50 nm particles showed poor drug penetration into the tumor.   Efficacy parameters w 1/2 , PA rel , FK PA , and FK eff are reported in Table 6. There was an increase in w 1/2 values from 0.003 cm for 25 nm particles to 0.235 cm for 300 nm particles. As the nanoparticle size further increased, w 1/2 decreased, until it reached 0.143 for 600 nm particles. The smallest PA rel among the compared nanoparticles was related to 25 nm nanoparticles (0.20%), and 300 nm nanoparticles with 55.02% had the largest PA rel among the studied MNPs. FK eff values for different sizes of MNPs after 60 min of treatment are presented in Figure 15. FK eff was calculated as the lowest value for 25 nm nanoparticles, which was equal to 0.11%, whereas it was the highest for 300 nm particles with 24.1%.

The Influence of Magnetic Strength (Brem)
In this section, the effect of magnetic strength by changing the remnant magnetic flux (Brem) in the range of 0.5-2.5 T on the treatment output is investigated. In Figure 16, the contours of CF distribution in the tumor are presented for the three magnetic strengths of 0.5, 1.5, and 2.5 T. As can be seen, an increase in the magnetic strength will increase the drug's penetration into the tumor. At the lowest magnetic strength (0.5 T), the drug penetration into the tumor was negligible, whereas at the highest magnetic strength (2.5 T), the drug particles significantly penetrated the tumor. 0.5 T 1.5 T 2.5 T Figure 16. Comparison of the mean concentration of free drug (CF) in the tumor after 60 min of treatment for different applied magnetic strengths. The results of this section show that the size of nanoparticles can affect the efficiency of treatment using this method. By selecting an appropriate MNP size, the effective killed cells fraction can increase by 9.5 times compared with conventional IPC. Whereas an improper selection of nanoparticle size (in the range of 25-100 nm) leads to a reduction in the treatment efficacy even lower than the conventional IPC. Additionally, it should be noted that experimental results [51] demonstrate that nanoparticles with a larger size have a longer residence in the peritoneal cavity and will be available to the tumor for a longer time.

The Influence of Magnetic Strength (B rem )
In this section, the effect of magnetic strength by changing the remnant magnetic flux (B rem ) in the range of 0.5-2.5 T on the treatment output is investigated. In Figure 16, the contours of C F distribution in the tumor are presented for the three magnetic strengths of 0.5, 1.5, and 2.5 T. As can be seen, an increase in the magnetic strength will increase the drug's penetration into the tumor. At the lowest magnetic strength (0.5 T), the drug penetration into the tumor was negligible, whereas at the highest magnetic strength (2.5 T), the drug particles significantly penetrated the tumor.
(Brem) in the range of 0.5-2.5 T on the treatment output is investigated. In Figure 16, the contours of CF distribution in the tumor are presented for the three magnetic strengths of 0.5, 1.5, and 2.5 T. As can be seen, an increase in the magnetic strength will increase the drug's penetration into the tumor. At the lowest magnetic strength (0.5 T), the drug penetration into the tumor was negligible, whereas at the highest magnetic strength (2.5 T), the drug particles significantly penetrated the tumor. 0.5 T 1.5 T 2.5 T Figure 16. Comparison of the mean concentration of free drug (CF) in the tumor after 60 min of treatment for different applied magnetic strengths. Figure 17 presents the mean CF profiles over time for different values of magnetic strengths. The mean free drug concentration 60 min after the injection for the three magnetic strengths of 0.5, 1.5, and 2.5 tesla is calculated as 0.00035, 0.00540, and 0.01464 mol/m 3 , respectively. For magnetic strengths of 1.5 and 2.5 tesla, there was an increase in the mean CF by 4.5 and 12.2 times, respectively. While for the magnetic strength of 0.5 tesla, CF was 70% less than the corresponding value for conventional IPC.  Figure 17 presents the mean C F profiles over time for different values of magnetic strengths. The mean free drug concentration 60 min after the injection for the three magnetic strengths of 0.5, 1.5, and 2.5 tesla is calculated as 0.00035, 0.00540, and 0.01464 mol/m 3 , respectively. For magnetic strengths of 1.5 and 2.5 tesla, there was an increase in the mean C F by 4.5 and 12.2 times, respectively. While for the magnetic strength of 0.5 tesla, C F was 70% less than the corresponding value for conventional IPC.  Table 7 presents the efficacy parameters w1/2, PArel, FKPA, and FKeff. According to the results, the penetration depth was 0.006 cm for the magnetic strength of 0.5 tesla, which was less than the corresponding value for the conventional method. For the magnetic strengths of 1.5 and 2.5 tesla, the values obtained for the penetration depth of the drug into the tumor were equal to 0.071 and 0.192 cm, respectively. This indicates an increase of 11.8 and 32 times in penetration depth for the magnetic strengths of 1.5 and 2.5, respectively. In addition, the PArel value for the magnetic strength of 0.5 was 0.95%, which is less than the calculated PArel for the conventional IPC. Meanwhile, the calculated values for larger magnetic strengths were 12.82% and 38.32% for the magnetic strengths of 1.5 and 2.5 tesla, which showed a 1.43 and 4.2 times increase compared to the conventional IPC, respectively.  Figure 18 shows the FKeff values for different magnetic strengths after 60 min of treatment. The FKeff values obtained for the three magnetic strengths of 0.5 tesla, 1.5 tesla, and 2.5 tesla were 0.43%, 6.52%, and 18.1%, respectively. Similar to other efficacy parameters,  Table 7 presents the efficacy parameters w 1/2 , PA rel , FK PA , and FK eff . According to the results, the penetration depth was 0.006 cm for the magnetic strength of 0.5 tesla, which was less than the corresponding value for the conventional method. For the magnetic strengths of 1.5 and 2.5 tesla, the values obtained for the penetration depth of the drug into the tumor were equal to 0.071 and 0.192 cm, respectively. This indicates an increase of 11.8 and 32 times in penetration depth for the magnetic strengths of 1.5 and 2.5, respectively. In addition, the PA rel value for the magnetic strength of 0.5 was 0.95%, which is less than the calculated PA rel for the conventional IPC. Meanwhile, the calculated values for larger magnetic strengths were 12.82% and 38.32% for the magnetic strengths of 1.5 and 2.5 tesla, which showed a 1.43 and 4.2 times increase compared to the conventional IPC, respectively.  Figure 18 shows the FK eff values for different magnetic strengths after 60 min of treatment. The FK eff values obtained for the three magnetic strengths of 0.5 tesla, 1.5 tesla, and 2.5 tesla were 0.43%, 6.52%, and 18.1%, respectively. Similar to other efficacy parameters, FK eff for the magnetic strength of 0.5 tesla was lower than the conventional method, while for the two larger magnetic strengths, it showed an increase by more than 2.5 and 7 times, respectively. This was due to the competition between the magnetic force exerted on the nanoparticles and the microenvironmental barriers, which means that the magnetic force was not large enough to overcome the microenvironmental barriers of the tumor at lower magnetic strengths. The results of this section show that increasing the magnetic strength improves the penetration of the MNPs into the tumor by increasing the magnetic force exerted on the nanoparticles.

The Effect of Tumor Distance from the Magnet
The distance of the permanent magnet from the tumor is another parameter that can influence the performance of the magnetically controlled IP drug targeting system. In reality, the tumor may be located in different sites of the peritoneum area. Thus, the distance between the magnet and the tumor may vary depending on the tumor's location. Therefore, it is essential to examine the effect of distance and the extent of its impact on the treatment output. In this section, by considering three different distances between the tumor and the permanent magnet (L), this parameter's effect on the treatment using MNPs was studied.
According to Equation (7), the magnetic force exerted on the MNPs was commensurate with ∇ , revealing the importance of the strength and gradient of the magnetic field. To better understand the effect of distance, the spatial distribution of the magnetic field around the permanent magnet is presented in Figure 19. As evident in this figure, by increasing the distance between the tumor and the magnet, the strength and the gradient of the magnetic field decreased significantly, which indicates that the distance of the magnet from the tumor can affect the efficiency of the magnetically controlled IP.

The Effect of Tumor Distance from the Magnet
The distance of the permanent magnet from the tumor is another parameter that can influence the performance of the magnetically controlled IP drug targeting system. In reality, the tumor may be located in different sites of the peritoneum area. Thus, the distance between the magnet and the tumor may vary depending on the tumor's location. Therefore, it is essential to examine the effect of distance and the extent of its impact on the treatment output. In this section, by considering three different distances between the tumor and the permanent magnet (L), this parameter's effect on the treatment using MNPs was studied.
According to Equation (7), the magnetic force exerted on the MNPs was commensurate with ∇ 2 H, revealing the importance of the strength and gradient of the magnetic field. To better understand the effect of distance, the spatial distribution of the magnetic field around the permanent magnet is presented in Figure 19. As evident in this figure, by increasing the distance between the tumor and the magnet, the strength and the gradient of the magnetic field decreased significantly, which indicates that the distance of the magnet from the tumor can affect the efficiency of the magnetically controlled IP.
According to Equation (7), the magnetic force exerted on the MNPs was commensurate with ∇ , revealing the importance of the strength and gradient of the magnetic field. To better understand the effect of distance, the spatial distribution of the magnetic field around the permanent magnet is presented in Figure 19. As evident in this figure, by increasing the distance between the tumor and the magnet, the strength and the gradient of the magnetic field decreased significantly, which indicates that the distance of the magnet from the tumor can affect the efficiency of the magnetically controlled IP.  Figure 20 shows the contours of C F in the tumor for three different distances of the permanent magnet from the tumor including 3, 5, and 10 cm. The comparison of the contours shows that by increasing the distance, a smaller part of the tumor was affected by the drug. This reduction in penetration occurred dramatically, especially by increasing the distance of the magnet from the tumor to 10 cm. According to Figure 21, C F values 60 min after the injection at a distance of 3 cm was 0.0069 mol/m 3 , while increasing the distance to 10 cm decreased the C F value to 0.0019 mol/m 3 .  Figure 20 shows the contours of CF in the tumor for three different distances of the permanent magnet from the tumor including 3, 5, and 10 cm. The comparison of the contours shows that by increasing the distance, a smaller part of the tumor was affected by the drug. This reduction in penetration occurred dramatically, especially by increasing the distance of the magnet from the tumor to 10 cm. According to Figure 21, CF values 60 min after the injection at a distance of 3 cm was 0.0069 mol/m 3 , while increasing the distance to 10 cm decreased the CF value to 0.0019 mol/m 3 .       Table 8 presents the values obtained for w1/2, PArel, FKPA, and FKeff for three distances between the magnet and the tumor after 60 min of treatment. Reducing the distance of the magnet from the tumor from 10 to 3 cm increased the w1/2 and PArel values by 4 and 4.4 times, respectively. Figure 22 compares the FKeff values for different magnetic distances after 60 min of treatment. FKeff increased by 3.88 times by reducing the distance of the magnet from the tumor from 10 to 3 cm.
According to the results of this section, changing the distance between the tumor and the magnet significantly affected the efficiency of the magnetically controlled IP chemotherapy. Considering that tumors may be located in different sites of the peritoneal cavity, treatment using MNPs will be more difficult for tumors in deeper parts of the body. Although, using optimal values for other parameters, such as magnetic nanoparticle size and magnetic strength, can help to improve drug delivery to these tumors.  Table 8 presents the values obtained for w 1/2 , PA rel , FK PA , and FK eff for three distances between the magnet and the tumor after 60 min of treatment. Reducing the distance of the magnet from the tumor from 10 to 3 cm increased the w 1/2 and PA rel values by 4 and 4.4 times, respectively. Figure 22 compares the FK eff values for different magnetic distances after 60 min of treatment. FK eff increased by 3.88 times by reducing the distance of the magnet from the tumor from 10 to 3 cm.

Validation
In this section, the validation of the numerical solution is presented. The IFP and IFV distributions were calculated using the Darcy law and, subsequently, the concentration distribution of the free, bound, and internalized drugs could be obtained by solving the mass transfer equations. Hence, each of these physics must be validated separately.  According to the results of this section, changing the distance between the tumor and the magnet significantly affected the efficiency of the magnetically controlled IP chemotherapy. Considering that tumors may be located in different sites of the peritoneal cavity, treatment using MNPs will be more difficult for tumors in deeper parts of the body. Although, using optimal values for other parameters, such as magnetic nanoparticle size and magnetic strength, can help to improve drug delivery to these tumors.

Validation
In this section, the validation of the numerical solution is presented. The IFP and IFV distributions were calculated using the Darcy law and, subsequently, the concentration distribution of the free, bound, and internalized drugs could be obtained by solving the mass transfer equations. Hence, each of these physics must be validated separately.

Validation of the IFP and IFV Distribution
Finding the IFP and IFV is one of the exigent steps of this modeling, which was obtained by solving Equation (1). To validate the distribution of IFP in the tumor, a comparison between the radial distribution of IFP with an experimental work by Boucher et al. [52] was conducted under the same conditions, which is presented in Figure 23a. As can be seen in this figure, there was a good agreement between the values obtained from the modeling and the experimental results. In addition, the values of IFV distribution along the tumor radius were compared with the theoretical values presented by Soltani and Chen [42] (Figure 23b), and there was good agreement between the values of these two graphs.

Validation of the Concentration Distribution
Au et al. [50] evaluated the performance of chemotherapy by IP injection on a mouse. In their study, the radius of the tumor was 2 mm, and 45 mM of paclitaxel was injected intraperitoneally. A comparison of the concentration profiles in terms of penetration depth within 6 h after injection is presented in Figure 24. As can be seen, there was a relatively good agreement between the two profiles. The difference between the concentration values in the two graphs resulted from different properties of the tissue and the drug in the two studies.

Validation of the Concentration Distribution
Au et al. [50] evaluated the performance of chemotherapy by IP injection on a mouse. In their study, the radius of the tumor was 2 mm, and 45 mM of paclitaxel was injected intraperitoneally. A comparison of the concentration profiles in terms of penetration depth within 6 h after injection is presented in Figure 24. As can be seen, there was a relatively good agreement between the two profiles. The difference between the concentration values in the two graphs resulted from different properties of the tissue and the drug in the two studies. Au et al. [50] evaluated the performance of chemotherapy by IP injection on a mouse. In their study, the radius of the tumor was 2 mm, and 45 mM of paclitaxel was injected intraperitoneally. A comparison of the concentration profiles in terms of penetration depth within 6 h after injection is presented in Figure 24. As can be seen, there was a relatively good agreement between the two profiles. The difference between the concentration values in the two graphs resulted from different properties of the tissue and the drug in the two studies.

Conclusions
In the present study, a mathematical model was developed to study a magnetically controlled IP drug targeting system as a solution to improve the drug penetration into the tumor. Considering the binding and internalization of the drug into the cancer cells, the mathematical model analyzed the drug distribution inside the tumor in three forms: free, bound, and internalized. The drug penetration area and the fraction of killed cells in the tumor were used to estimate the efficiency of the proposed drug delivery system. The main conclusions of the present study are as follows:

Conclusions
In the present study, a mathematical model was developed to study a magnetically controlled IP drug targeting system as a solution to improve the drug penetration into the tumor. Considering the binding and internalization of the drug into the cancer cells, the mathematical model analyzed the drug distribution inside the tumor in three forms: free, bound, and internalized. The drug penetration area and the fraction of killed cells in the tumor were used to estimate the efficiency of the proposed drug delivery system. The main conclusions of the present study are as follows: 1.
Using a magnetically controlled IP drug targeting system can overcome the microenvironmental barriers of the tumor against the transfer of drug particles and significantly increase the treatment efficiency compared to conventional IP chemotherapy; 2.
There was an optimal size for MNPs: larger nanoparticles exerted a stronger magnetic force. On the other hand, tumor microenvironmental barriers further hinder the movement of large nanoparticles in the tissue compared to smaller ones. Therefore, the competition between these two factors (i.e., the magnetic force to the particle and the hindering effects of the microenvironmental barriers) determines the extent of penetration of different sizes of nanoparticles; 3.
Using a permanent magnet with more magnetic strength and reducing the magnet's distance from the tumor by increasing the magnetic force exerted on the MNPs improved IP delivery of doxorubicin-coated MNPs.
There are a number of assumptions and limitations in our model. With the lack of appropriate experimental data associated with the drug concentration at the outer edge of the tumor, this value was assumed constant during the simulation. By sampling the intraperitoneal fluid during treatment, a more accurate boundary condition for concentration outside the tumor can be obtained. In addition, the tissue poroelasticity effect was not considered in this study. This effect should be considered in the future developments of the model, given that studies have shown the effect of these phenomena on drug binding and internalization [53]. The parameters used in our model were derived from prior experimental investigations that were consistent with the numerical studies too. We adopted the spherical symmetric approximation assuming that the infusion occurs within a circular domain. Although the present study does not claim to be a replacement for experimental studies, the mathematical model of magnetically controlled IP drug targeting is a powerful tool that can help us gain insights into the impacts of various parameters.
Supplementary Materials: The following supporting information can be downloaded at: https: //www.mdpi.com/article/10.3390/pharmaceutics14020324/s1, Figure S1: The average drug concen-tration in the tumor for conventional IP chemotherapy using Doxorubicin, Paclitaxel, and Cisplatin, compared with the result of magnetically controlled IP drug delivery during 60 minutes of treatment.