In this chapter, based on (Russo, Numerical simulation of magneto-hydro-dynamics inside cardiovascular and respiratory systems. University of Rome Tor Vergata, A.A, 2014/2015), we present the results of the interactions between a magnetic field, generated by a probe and the magnetic nanoparticles transported by the blood in the hepatic artery. The geometry of the hepatic artery is extracted from DICOM images, taken with a CT scan, and the reconstruction is done with VMTK. As far as the blood velocity is concerned, we use the Womersley solution, (Womersley, J Physiol 127:553–563, 1955), to represent a realistic flow rate from literature data, (Tang, Abdominal aortic hemodynamics in young healthy adults at rest and during lower limb exercise: quantification using image-based computer modeling. American Journal of Physiology—Heart and Circulatory Physiology, 2006). Three different meshes are generated in order to evaluate the grid independence of the solution. The optimal grid is evaluated on the basis of the pressure and wall shear stress results. A new solver, coupling the Lagrangian dynamics of the nanoparticles with the Eulerian dynamics of the blood, is implemented in OpenFOAM to perform the simulations. A single rectangular coil is used to generate the external magnetic field. The resistive pressure, Womersley’s profile for the inlet velocity and the magnetic field of a rectangular coil are implemented in the software as boundary conditions. The results show the influence of the position of the probe, as well as the limitations associated with the rectangular coil configuration.