{'fraction_non_alphanumeric': 0.03228115847705825, 'fraction_numerical': 0.00986007159127888, 'mean_word_length': 4.841284926819996, 'pattern_counts': {'":': 0, '<': 3, '<?xml version=': 0, '>': 0, 'https://': 0, 'lorem ipsum': 0, 'www.': 0, 'xml': 0}, 'pii_count': 0, 'substrings_counts': 0, 'word_list_counts': {'cursed_substrings.json': 1, 'profanity_word_list.json': 0, 'sexual_word_list.json': 0, 'zh_pornsignals.json': 0}}
Optimization of Field-Free Point Position, Gradient Field and Ferromagnetic Polymer Ratio for Enhanced Navigation of Magnetically Controlled Polymer-Based Microrobots in Blood Vessel

Microscale and nanoscale robots, frequently referred to as future cargo systems for targeted drug delivery, can effectively convert magnetic energy into locomotion. However, navigating and imaging them within a complex colloidal vascular system at a clinical scale is exigent. Hence, a more precise and enhanced hybrid control navigation and imaging system is necessary. Magnetic particle imaging (MPI) has been successfully applied to visualize the ensemble of superparamagnetic nanoparticles (MNPs) with high temporal sensitivity. MPI uses the concept of field-free point (FFP) mechanism in the principal magnetic field. The gradient magnetic field (|∇B|) of MPI scanners can generate sufficient magnetic force in MNPs; hence, it has been recently used to navigate nanosized particles and micron-sized swimmers. In this article, we present a simulation analysis of the optimized navigation of an ensemble of microsized polymer MNP-based drug carriers in blood vessels. Initially, an ideal two-dimensional FFP case is employed for the basic optimization of the FFP position to achieve efficient navigation. Thereafter, a nine-coil electromagnetic actuation simulation system is developed to generate and manipulate the FFP position and |∇B|. Under certain vessel and fluid conditions, the particle trajectories of different ferromagnetic polymer ratios and |∇B| were compared to optimize the FFP position.


Introduction
Despite spectacular medical science innovations and studies in various advanced cancer treatments, cancer remains the second leading cause of death [1]. The conventional therapeutic methods include chemotherapy, radiation therapy, and photodynamic therapy (PDT). However, these traditional therapies have certain limitations; for example, in chemotherapy the entire body is exposed to high concentrations of toxic drugs, which cause severe damage to healthy tissue despite killing cancer cells; radiation therapy might induce severe side effects on the human body due to the radiation exposure. Moreover, PDT is expensive but it is also difficult to target cancer tissue for photosensitizing, and there is limited knowledge for clinical cancer treatment [2,3]. Chemotherapy is the most frequently used method for cancer treatment. There are several major classifications of anticancer drugs. Doxorubicin (DOX) is widely used to treat cancer as a broad-spectrum anticancer agent with an optimistically strong effect on solid tumors. DOX is a hydrophilic anthracycline antibiotic that, when inserted into the DNA of cancer cells, destroys the DNA double helix [4]. However, DOX not only acts on cancer cells but also badly affects normal cells and may also cause cardiac toxicity depending on the cumulative dose. Therefore, it is crucial to reduce toxicity to normal tissues without affecting its ability to kill cancer tissues.
Targeted drug delivery (TDD) can precisely target and deliver an anticancer drug to infected tissue, which has emerged as a prominent solution to the toxicity problem [5,6]. A wide variety of both nano-and microscale carriers have been utilized as TDD vehicles. These carriers include both organic and inorganic materials. Some are created from metals, while others are made from polymers [3,7,8]. Different approaches have been suggested for steering these TDD vehicles; the use of acoustic radiation force, thermo-electromagnetically reacting bilayer-structured microrobots, octagram-shaped microgrippers, and macrophagebased microrobots have been proposed [9][10][11]. However, the use of electromagnetic actuators (EMAs) is the most widely utilized technique for navigating both nano-and microsized magnetic drug carriers [12][13][14]. The first human trial was reported in 1996 when cancer drugs were attached to iron core particles with diameters of 100 nm and then steered by an external magnet field to treat tumors [15].
Along with their biocompatibility, these small-scale carriers must satisfy several other requirements, such as the capability to move in different types of liquids and blood vessels. They must also be capable of overcoming obstacles and penetrating complex barriers. Finding carriers that can satisfy all these requisites is exigent, and the solutions formulated to resolve them usually require the combination of multiple materials on the same platform. For example, microrobots made of gelatin/polymer/hydrogel, encapsulating superparamagnetic nanoparticles (MNPs), and carrying DOX particles are potential candidates [16]. The magnetic properties of MNPs can be manipulated for navigation using external magnetic fields. First, the microrobot reaches a predetermined target location by the gradient magnetic field of the electromagnetic actuation system. Next, after near-infrared (NIR) irradiation, the gelatin/PVA of the microrobot is decomposed, and the MNPs and DOX drug particles are left in the target area. The disassembled MNPs are recovered from the target lesion by the magnetic field of the EMA system, only DOX particles remain in the target area to generate a therapeutic effect in the target lesion only. Therefore, MNP-mediated polymer-based carriers have emerged as one of the most promising cargo systems for TDD [3].
For efficient TDD, a real-time tracking of MNPs and microrobots in vivo is required. Several methods for developing a targeting scheme with feedback control, such as using ultrasound to locate solid microsized particles [17] or a microscope to track visible particles [18], have been investigated. However, none of the currently available clinical imaging modalities satisfies all the requirements. In the literature, the imaging aspect of magnetic nanoparticles has been thoroughly discussed, and the advantages and limitations of each imaging technique have been extensively described [19]. A major problem is to identify a specific navigation method that fits a particular imaging technique; a hybrid navigation and imaging method is highly preferable. The magnetic resonance navigation (MRN) system was developed based on magnetic resonance imaging (MRI) to steer and capture images of nanoparticles simultaneously using electromagnetic coils [20]; MRI is already clinically available and can provide relatively improved resolution and contrast of images. In vivo investigations of the MRN system were successfully conducted in the artery of a living animal using a magnetic bead with a diameter of 1.5 mm [21]. However, to date, the MRI of nanoscale and microscale particles has not been fully successful; thus, a new approach capable of both guiding and imaging magnetic nanoparticles should be conceived.
Magnetic particle imaging (MPI) is a tracer-based (i.e., MNP-based) imaging technique that has gained considerable interest from MRI and TDD researchers [22]. Due to the unique magnetic features of MNPs, they can function at both the cellular and molecular levels. They can also produce a unique secondary magnetization signal as a reaction to an external magnetic field. Hence, they can suitably be used as contrast agents in MPI and drug carrier particles in TDD. Moreover, MPI scanners have been observed to achieve higher temporal and spatial resolutions and have high potentials for revolutionizing the field of biomedical imaging and TDD [23].
The fundamental principle of MPI is the detection of nonlinear responses of superparamagnetic iron oxide particles induced by a secondary oscillating magnetic field in a specific region, where the value of the primary magnetic field is zero-i.e., a field-free point (FFP). These FFPs can be generated and scanned by several combinations of permanent magnets and electromagnetic coils [24]. The simplest magnetic arrangement involves using four permanent magnets arranged in the form of a quadrupole to achieve static FFPs and an additional device to drive the FFP along a suitable trajectory [22,25]. The main disadvantage of this method is its low scanning speed. Other methods include the use of supportive drive coils that generate homogenous magnetic fields to drive a static FFP to a designated position [26]. The development of improved high-end MPI scanners has motivated TDD researchers to employ MPI modules for the navigation of magnetic carriers. Nothnagel et al. demonstrated this MPI capability by steering soft magnetic spheres by varying the FFP position using an MPI coil [27]. This somehow resolved one of the major problems encountered in an EMA system with MPI imaging-i.e., the problem of high-gradient fields interfering with the steering and tracking of nanoparticles in real time. However, the author did not consider any fluidic flow, MNPs are moved in stationary fluid. Recently, a hybrid guidance system capable of performing both the actuation and monitoring of nanoparticles was reported [28]. The paper proposed the stepwise navigation and two-dimensional (2D) tracking of nanoparticles in real time using a feedback control method. Simultaneous MPI imaging and MNP navigation in bifurcation flow experiments have also been recently demonstrated [29].
In this article, we present a simulation analysis of the optimized navigation of an ensemble of microsized polymer-MNP-based drug carriers in blood vessels ( Figure 1). Initially, an ideal two-dimensional FFP case is employed for the basic optimization of the FFP position to achieve efficient navigation. Thereafter, a nine-coil electromagnetic actuation simulation system is developed to generate and manipulate the FFP position and |∇B|. Under certain vessel and fluid conditions, the trajectories of microrobots with different ferromagnetic polymer ratios and |∇B| were compared and optimized for efficient navigation.

Analytical Model
Consider a spherical ferromagnetic microrobot of radius r immersed in blood flowing in a vessel (Figure 1c). The microrobot is governed by several forces: magnetic force ( ), blood hydrodynamic drag force ( ), gravitational force ( ), and contact force between MNPs and the vessel wall ( ). The robot's translational motion is given by where is the translational velocity and is the robot mass.

Hydrodynamic Drag Force
The low-velocity approximation of nonturbulent creeping flow was considered, and blood was assumed to be incompressible. Using the Stokes flow model, the hydrodynamic drag force ( ) acting on a spherical body of radius r in the fluid is expressed as

Analytical Model
Consider a spherical ferromagnetic microrobot of radius r immersed in blood flowing in a vessel ( Figure 1c). The microrobot is governed by several forces: magnetic force (F m ), blood hydrodynamic drag force (F d ), gravitational force (F g ), and contact force between MNPs and the vessel wall (F c ). The robot's translational motion is given by where v is the translational velocity and m p is the robot mass.

Hydrodynamic Drag Force
The low-velocity approximation of nonturbulent creeping flow was considered, and blood was assumed to be incompressible. Using the Stokes flow model, the hydrodynamic drag force (F d ) acting on a spherical body of radius r in the fluid is expressed as where m p is the particle mass; τ p is the particle velocity response time (SI unit: s); and u and v are the velocities of the fluid and robot, respectively. The particle velocity response time, however, depends on the flow conditions and parameters. The Reynolds number for laminar (Stokes) flow must satisfy Re 1 to be valid; this number is given by where η is the dynamic viscosity of the fluid, and d p is the diameter of microparticles. The particle velocity response time of spherical particles in a laminar flow is defined as follows: For a relatively moderate Re (i.e., 1 < Re < 800), the drag force can be set according to the Schiller-Naumann drag law [30]. Accordingly, the particle velocity response time is redefined as follows: where C D = 24 Re 1 + 0.15 × Re 0.687 .

Magnetic Force
The magnetic gradient induces a magnetic motive force on the nanoparticles and this process is known as magnetophoresis [31]. Assuming that the nanoparticles are homogeneous spheres with radius R and magnetization → M in a magnetically linear fluid of permeability µ 1 . When subjected to a magnetic intensity → H, then the effective magnetic moment m e f f is defined as: where µ 0 = 4π × 10 −7 N/A 2 is the permeability of free space. The magnetophoretic force on a single MNP can be expressed in the gradient of the magnetic intensity (∇ → H) and the effective magnetic dipole moment (m e f f ) of the MNPs.
The magnetization and magnetic intensity are linked by: where χ = µ 2 /µ 0 − 1, is the magnetic susceptibility of the particles. In this case, Equation The magnetophoretic force can be expressed in the gradient of the magnetic intensity (∇ → H) and the effective magnetic dipole moment (m e f f ) of the MNPs Combining Equations (8) and (9) for effective moment, the magnetic force on single MNP in a gradient magnetic field may be written as The nonlinear behavior in magnetic materials is very dominant and it cannot be ignored when modeling the electromechanics of magnetic particles. These materials can be classified as hard and soft magnetic materials. In both cases, the magnetization is dependent on magnetic field intensity. The most important phenomena in the nonlinear behavior of magnetic materials is the saturation, which limits the magnitude of the magnetization vector to a finite value, M sat . Considering that the MNPs are magnetically soft, and the applied magnetic field is strong, the particles will saturate so that Similarly, let us consider a spherical microrobot which is made of nonmagnetic binding polymer + MNPs. The magnetic force on this microparticle is defined by the quantity of MNPs in the microparticle. Therefore, if we define a volumetric ratio parameter, τ m = V m V , then the effective force on the microrobot will be modified as: where V and τ m M sat are the volume and saturation magnetization of the microrobot, respectively; → B is the external magnetic field; and τ m is the ferromagnetic ratio [32].

Gravitational Force
The gravitational force ( → F g ) can be interpreted as the apparent weight-i.e., combined action of weight and fluid buoyancy: In the foregoing, ρ = τ m ρ m + (1τ m )ρ poly , where ρ m and ρ poly are the densities of the magnetic material and polymer, respectively; ρ f is the carrier fluid density; and g is the gravitational acceleration.

Contact Force
The contact force results from particle-vessel collisions. This effect can be introduced in different ways, such as general reflection, elastic bouncing, and electrostatic attraction. Only the general reflection and bouncing effects have been considered in this current work. The robot velocity after collision is given by where v c is the velocity of the robot before collision, and n is the normal vector on the wall surface. In order to reduce computational time, we have neglected of electrostatic interactions.

Ideal 2D Case
The COMSOL Multiphysics ® software (5.4, COMSOL Inc, Burlington, MA, USA) was applied to perform all numerical modeling and animated simulation results with the table of animation list (Table S1) are in Supplementary Materials. The authors started with the simplest case, i.e., an ideal 2D FFP generation using a set of equations in the magnetic field, no current (mfnc) module. Then, the microrobots were navigated and tracked in 2D blood vessel laminar flow (spf) and particle tracing in fluid flow (fpt) modules consecutively. The vessel geometry was a 2D bifurcation with a diameter, total length, and bifurcation angle of 5 mm, 30 mm, and 40 • , respectively (Figure 2b). A fluid initial velocity of 1 cm/s was set at the inlet, and the velocities at the outlets were controlled under suitable pressure conditions. The Reynolds number in a typical microfluidic channel is relatively small (Re < 1); hence, the flow can be considered to be laminar. The velocity profile of the developed flow is shown in Figure 2a. Each microrobot has a diameter of 30 µ m. Their magnetic properties can be controlled using the ferromagnetic ratio ( ), which actually controls the volume concentration of MNPs in the binding polymer. In the 2D simulation, = 0.2 for most cases. A detailed comparison between the magnetic force and is explained in Section 3.2.  Figure 3a shows that, without external forces, the robots simply follow the fluid flow with the drag force; hence, they are uniformly dispersed in the vessel. Using the particle counter feature of COMSOL, the exact number of microrobots reaching both outlets was monitored. In the absence of FFP, an equal number of microrobots, i.e., 75, reached each outlet with the same temporal pattern. When the FFP was applied, the microrobots reacted to the gradient field and their trajectories deviated. However, the trajectories are considerably dependent on the FFP position. Figure 2b shows the FFP positions used in the simulations. When the FFP was placed at position ⑥, 80% of the microrobots were navigated to the left bifurcation (Figure 3b). We also simulated the trajectories of microrobots creeping towards right bifurcation when FFP is positioned at ① and ③ (Figure 4). It is evident that forces microrobots toward right bifurcation while acts downward, decelerating and allowing microrobots to linger with . The FFP strength and position were manually controlled using the mfnc module; for each FFP position, particle trajectories were traced and compared. Considering clinical limitations, we restricted our maximum gradient field to |∇B| = 1 T/m ∼ 0.8 × 10 6 A/m 2 . An ensemble of 150 microrobots were released at the inlet with an initial velocity of 0.5 cm/s. These were released in three steps with 50 microrobots in each temporal step. Each microrobot has a diameter of 30 µm. Their magnetic properties can be controlled using the ferromagnetic ratio (τ m ), which actually controls the volume concentration of MNPs in the binding polymer. In the 2D simulation, τ m = 0.2 for most cases. A detailed comparison between the magnetic force and τ m is explained in Section 3.2. Figure 3a shows that, without external forces, the robots simply follow the fluid flow with the drag force; hence, they are uniformly dispersed in the vessel. Using the particle counter feature of COMSOL, the exact number of microrobots reaching both outlets was monitored. In the absence of FFP, an equal number of microrobots, i.e., 75, reached each outlet with the same temporal pattern. When the FFP was applied, the microrobots reacted to the gradient field and their trajectories deviated. However, the trajectories are considerably dependent on the FFP position. Figure 2b shows the FFP positions used in the simulations. When the FFP was placed at position 6 , 80% of the microrobots were navigated to the left bifurcation (Figure 3b). We also simulated the trajectories of microrobots creeping towards right bifurcation when FFP is positioned at 1 and 3 (Figure 4). It is evident that ∇B x forces microrobots toward right bifurcation while ∇B y acts downward, decelerating and allowing microrobots to linger with ∇B x .

FFP Generation
The COMSOL EMA simulation model the adoption of our experimentally validated EMA system; it has the same physical conditions as the coils (e.g., number of turns and scales) [33]. It can generate a magnetic field gradient of up to |∇B| ∼ 3 T/m, and each coil can support a maximum current of 20 A. The region of interest (ROI) in our system is 40 × 40 mm 2 . This nine-coil EMA setup not only generates higher gradient fields but also provides the gradient field in the z-direction, thus minimizing the use of a mechanical microstage in most conventional four-coil EMA setups and can control the microrobots in 5-DOFs [34]. The detailed parameters of the EMA coil system are summarized in Table 1.   To measure the current (I) combination for a certain FFP at any given point P(x, y, z) in the workspace, the magnetic field created by the nth single electromagnetic coil can be expressed by vector B n (P). This is computed by the product of the input current (i) and magnetic field per unit current (b n (P)), as follows: B n (P) = [ B x,n (P) B y,n (P) B z,n (P) ] = [b x,n (P)b y,n (P)b z,n (P) ]i. (15) Figure 3. (a) Particle trajectories without any external magnetic field influence; (b) FFP at position ⑥.  Using the superposition properties of electromagnetic coils, the resultant magnetic field at point P for the n-coil system can be expressed by a combination of linearly independent magnetic field vectors: whereb(P) ∈ R 3×9 is the mapping matrix in tesla per ampere from the current input matrix (i) to the magnetic field. The partial derivative of the magnetic field (B) in the x, y, and z directions can similarly be derived to compute the magnetic force in Equation (14): Equation (14) can be rewritten in the following matrix form: By combining Equations (16)- (18), the resultant magnetic field and magnetic force acing on the unit volume object can be calculated as where D = [B F] T is the desired matrix and X u ∈ R 12×9 represents the conversion matrix corresponding to the nine coils in the ROI. To obtain the value of the input current in Equation (19), a matrix inversion was applied to X u to obtain the current matrix via the following equation, which is analogous to the methodology reported in [35]: With B x = B y = B z = 0 at P(x, y, z), Maxwell's equations, and G xx + G yy + G zz = 0, the suitable solution for i was obtained using the COMSOL built-in coefficient-based partial differential equation solver [25]. The simulation results of FFP at two different positions based on the proposed i solution are depicted in Figure 5b.
where = [ ] is the desired matrix and ∈ 12×9 represents the conversion matrix corresponding to the nine coils in the ROI. To obtain the value of the input current in Equation (19), a matrix inversion was applied to to obtain the current matrix via the following equation, which is analogous to the methodology reported in [36]: With = = = 0 ( , , ), Maxwell's equations, and Gxx + Gyy + Gzz = 0, the suitable solution for i was obtained using the COMSOL built-in coefficient-based partial differential equation solver [25]. The simulation results of FFP at two different positions based on the proposed i solution are depicted in Figure 5b.

Microrobot Navigation
The vessel geometry is simply a modification of the 2D case with the same Y-shaped bifurcation (Figure 6a); the diameter, width, and height are 6, 30, and 40 mm, respectively. A maximum fluid velocity of 1.5 cm/s was set by controlling the pressure at the inlet and outlets. A nonpulsating steady creeping flow of incompressible fluid was considered. The cross-sectional velocity profile is shown in Figure 6b. In the case of blood, the viscosity (η) depends on both vessel diameter and hematocrit rate; however, in small vessels, its effect is virtually negligible [37].

Microrobot Navigation
The vessel geometry is simply a modification of the 2D case with the same Y-shaped bifurcation (Figure 6a); the diameter, width, and height are 6, 30, and 40 mm, respectively. A maximum fluid velocity of 1.5 cm/s was set by controlling the pressure at the inlet and outlets. A nonpulsating steady creeping flow of incompressible fluid was considered. The cross-sectional velocity profile is shown in Figure 6b. In the case of blood, the viscosity (η) depends on both vessel diameter and hematocrit rate; however, in small vessels, its effect is virtually negligible [36].
For a better comparison with 2D results, similar conditions were employed. All the animated simulation results of microrobot navigation in bifurcation flow with FFP field and the table of animation list (Table S1) are in Supplementary Materials. An ensemble of 150 microrobots were released at the inlet, each with an initial velocity of 0.5 cm/s. The microrobots were released in three steps with 50 microrobots in each temporal step. Each microrobot has a diameter of 30 µm, and the magnetic properties were controlled using the ferromagnetic ratio (τ m ). Figure 7 shows the comparison of microrobot trajectories with and without FFP at certain positions for τ m = 0.5 and ∇B = 1 T/m. Different from the 2D simulations where the magnetic fields at points 1 and 6 are mirror images of each other, the ∇B generated by the EMA system is slightly different. The difference in ∇B was observed after checking the particle trajectories at both points in Figure 7. There were six and eight particles on the right bifurcation at points 1 and 6 , respectively, with a slight difference in velocities. This difference is caused by the five coils on the lower side that are not symmetrically arranged relative to these points, resulting in slightly different FFP patterns. For a better comparison with 2D results, similar conditions were employed. An e semble of 150 microrobots were released at the inlet, each with an initial velocity 0.5 cm/s. The microrobots were released in three steps with 50 microrobots in each tem poral step. Each microrobot has a diameter of 30 μm, and the magnetic properties we controlled using the ferromagnetic ratio ( ). Figure 7 shows the comparison of micr robot trajectories with and without FFP at certain positions for = 0.5 and = 1 T/m Different from the 2D simulations where the magnetic fields at points ① and ⑥ are mi ror images of each other, the generated by the EMA system is slightly different. Th difference in was observed after checking the particle trajectories at both points in Fi ure 7. There were six and eight particles on the right bifurcation at points ① and ⑥ respectively, with a slight difference in velocities. This difference is caused by the five coi on the lower side that are not symmetrically arranged relative to these points, resulting slightly different FFP patterns.  Figure 8 shows the dependence of navigation efficiency on τ m . With increasing τ m , the microrobots were observed to move slower as they dragged across the inner wall of the blood vessel. Although no significant particle-wall interactions were considered apart from general reflection, the accumulation of static particles on the wall cannot be ignored. Therefore, to achieve efficient navigation without blood vessel blockage, it is important to choose a combination of τ m and ∇B suitable for the fluid flow conditions. In the literature, the reported average value of τ m is between 0.1 and 0.2. The authors also attempted to find a suitable ∇B for τ m = 0.2. The comparison of particle trajectories with ∇B = 1.5 and 2 T/m is shown in Figure 9. forces microrobots toward left bifurcation and decelerates microrobots, allowing them to linger with . Figure 8 shows the dependence of navigation efficiency on . With increasing , the microrobots were observed to move slower as they dragged across the inner wall of the blood vessel. Although no significant particle-wall interactions were considered apart from general reflection, the accumulation of static particles on the wall cannot be ignored. Therefore, to achieve efficient navigation without blood vessel blockage, it is important to choose a combination of and suitable for the fluid flow conditions. In the literature, the reported average value of is between 0.1 and 0.2. The authors also attempted to find a suitable for = 0.2 . The comparison of particle trajectories with = 1.5 and 2 T/m is shown in Figure 9.

Conclusions
Our simulation model was mainly designed for polymer-based microparticles, and we investigated the effect of FFP position, gradient field, and ferromagnetic ratio on navigation efficiency. We considered a relatively realistic approach to model bloodparticle interactions in a 3D bifurcation vessel. The vessel diameter, blood velocity, and pressure were all realistic. The microcarriers were navigated in both directions using the FFP, and the navigation efficiency was improved for optimal gradient field and ferromagnetic ratio.
For realistic microparticle navigation in vivo experiments, all possible particleparticle and particle-wall interactions should be included in the simulation model. The role of frictional, Van der Waals, electrostatic, and steric contact forces will be vital and may result in the sticking of microparticle to vessel walls. The exact 3D blood map of each patient should be extracted using an imaging modality such as MPI/MRI. Subsequently,

Conclusions
Our simulation model was mainly designed for polymer-based microparticles, and we investigated the effect of FFP position, gradient field, and ferromagnetic ratio on navigation efficiency. We considered a relatively realistic approach to model bloodparticle interactions in a 3D bifurcation vessel. The vessel diameter, blood velocity, and pressure were all realistic. The microcarriers were navigated in both directions using the FFP, and the navigation efficiency was improved for optimal gradient field and ferromagnetic ratio.
For realistic microparticle navigation in vivo experiments, all possible particleparticle and particle-wall interactions should be included in the simulation model. The role of frictional, Van der Waals, electrostatic, and steric contact forces will be vital and may result in the sticking of microparticle to vessel walls. The exact 3D blood map of each

Conclusions
Our simulation model was mainly designed for polymer-based microparticles, and we investigated the effect of FFP position, gradient field, and ferromagnetic ratio on navigation efficiency. We considered a relatively realistic approach to model blood-particle interactions in a 3D bifurcation vessel. The vessel diameter, blood velocity, and pressure were all realistic. The microcarriers were navigated in both directions using the FFP, and the navigation efficiency was improved for optimal gradient field and ferromagnetic ratio.
For realistic microparticle navigation in vivo experiments, all possible particle-particle and particle-wall interactions should be included in the simulation model. The role of frictional, Van der Waals, electrostatic, and steric contact forces will be vital and may result in the sticking of microparticle to vessel walls. The exact 3D blood map of each patient should be extracted using an imaging modality such as MPI/MRI. Subsequently, multiple FFPs at different bifurcations of blood vessels should be considered and investigated. Our future work will include all these interactions and multiple FFPs. We will then validate our model experimentally both in vitro and as well as in vivo.