Computational Assessment of Unsteady Flow Effects on Magnetic Nanoparticle Targeting Efficiency in a Magnetic Stented Carotid Bifurcation Artery.

PURPOSE: Worldwide, cardiovascular disease is the leading cause of hospitalization and death. Recently, the use of magnetizable nanoparticles for medical drug delivery has received much attention for potential treatment of both cancer and cardiovascular disease. However, proper understanding of the interacting magnetic field forces and the hydrodynamics of blood flow is needed for effective implementation. This paper presents the computational results of simulated implant assisted medical drug targeting (IA-MDT) via induced magnetism intended for administering patient specific doses of therapeutic agents to specific sites in the cardiovascular system. The drug delivery scheme presented in this paper functions via placement of a faintly magnetizable stent at a diseased location in the carotid artery, followed by delivery of magnetically susceptible drug carriers guided by the local magnetic field. Using this method, the magnetic stent can apply high localized magnetic field gradients within the diseased artery, while only exposing the neighboring tissues, arteries, and organs to a modest magnetic field. The localized field gradients also produce the forces needed to attract and hold drug-containing magnetic nanoparticles at the implant site for delivering therapeutic agents to treat in-stent restenosis. METHODS: The multi-physics computational model used in this work is from our previous work and has been slightly modified for the case scenario presented in this paper. The computational model is used to analyze pulsatile blood flow, particle motion, and particle capture efficiency in a magnetic stented region using the magnetic properties of magnetite (Fe3O4) and equations describing the magnetic forces acting on particles produced by an external cylindrical electromagnetic coil. The electromagnetic coil produces a uniform magnetic field in the computational arterial flow model domain, while both the particles and the implanted stent are paramagnetic. A Eulerian-Lagrangian technique is adopted to resolve the hemodynamic flow and the motion of particles under the influence of a range of magnetic field strengths (Br = 2T, 4T, 6T, and 8T). Particle diameter sizes of 10 nm-4 µm in diameter were evaluated. Two dimensionless numbers were evaluated in this work to characterize relative effects of Brownian motion (BM), magnetic force induced particle motion, and convective blood flow on particle motion. RESULTS: The computational simulations demonstrate that the greatest particle capture efficiency results for particle diameters within the micron range of 0.7-4 µm, specifically in regions where flow separation and vortices are at a minimum. Similar to our previous work (which did not involve the use of a magnetic stent), it was also observed that the capture efficiency of particles decreases substantially with particle diameter, especially in the superparamagnetic regime. Contrary to our previous work, using a magnetic stent tripled the capture efficiency of superparamagnetic particles. The highest capture efficiency observed for superparamagnetic particles was 78% with an 8 T magnetic field strength and 65% with a 2 T magnetic field strength when analyzing 100 nm particles. For 10 nm particles and an 8 T magnetic field strength, the particle capture efficiency was 55% and for a 2 T magnetic field strength the particle capture efficiency was observed to be 43%. Furthermore, it was found that larger magnetic field strengths, large particle diameter sizes (1 µm and above), and slower blood flow velocity improves the particle capture efficiency. The distribution of captured particles on the vessel wall along the axial and azimuthal directions is also discussed. Results for captured particles on the vessel wall along the axial flow direction showed that the particle density decreased along the axial direction, especially after the stented region. For the entrance section of the stented region, the captured particle density distribution along the axial direction is large, corresponding to the center-symmetrical distribution of the magnetic force in that section. CONCLUSION: The simulation results presented in this work have shown to yield favorable capture efficiencies for micron range particles and superparamagnetic particles using magnetized implants such as the stent discussed in this work. The results presented in this work justify further investigation of MDT as a treatment technique for cardiovascular disease.

Development of a Two-Way Coupled Eulerian-Lagrangian Computational Magnetic Nanoparticle Targeting Model for Pulsatile Flow in a Patient-Specific Diseased Left Carotid Bifurcation Artery.

PURPOSE: The aim of the present work is to present the development of a computational two-way coupled (fluid and particle coupled) magnetic nanoparticle targeting model to investigate the efficacy of magnetic drug targeting (MDT) in a patient-specific diseased left carotid bifurcation artery. MDT of therapeutic agents using multifunctional carrier particles has the potential to provide effective treatment of both cancer and cardiovascular disease by enabling a variety of localized treatment and diagnostic modalities while minimizing side effects. METHODS: A computational model is developed to analyze pulsatile blood flow, particle motion, and particle capture efficiency in a diseased left carotid bifurcation artery using the magnetic properties of magnetite (Fe3O4) and equations describing the magnetic forces acting on particles produced by an external cylindrical electromagnetic coil. A Eulerian-Lagrangian technique is adopted to resolve the hemodynamic flow and the motion of particles under the influence of a magnetic field (Br= 2T). Particle diameter sizes of 20 nm-4 µm in diameter were considered. RESULTS: The computational simulations demonstrate that the greatest particle capture efficiency results for particle diameters within the micron range, specifically 4 µm in regions where flow separation and vortices are at a minimum. It was also determined that the capture efficiency of particles decreases substantially with particle diameter, especially in the superparamagnetic regime. Particles larger than 2 µm were targeted and captured at the desired location by the external magnetic field, and the largest capture efficiency observed was approximately 98%. CONCLUSION: The simulation results presented in the present work have shown to yield favorable capture efficiencies for micron range particles and a potential for enhancing capture efficiency of superparamagnetic particles in smaller arteries and/or using magnetized implants such as cardiovascular stents. The present work presents results for justifying further investigation of MDT as a treatment technique for cardiovascular disease.

Development of an Experimental and Digital Cardiovascular Arterial Model for Transient Hemodynamic and Postural Change Studies: "A Preliminary Framework Analysis".

The ultimate goal of the present work is to aid in the development of tools to assist in the treatment of cardiovascular disease. Gaining an understanding of hemodynamic parameters for medical implants allow clinicians to have some patient-specific proposals for intervention planning. In the present work an experimental and digital computational fluid dynamics (CFD) arterial model consisting of a number of major arteries (aorta, carotid bifurcation, cranial, femoral, jejunal, and subclavian arteries) were fabricated to study: (1) the effects of local hemodynamics (flow parameters) on global hemodynamics (2) the effects of transition from bedrest to upright position (postural change) on hemodynamics, and (3) diffusion of dye (medical drug diffusion simulation) in the arterial system via experimental and numerical techniques. The experimental and digital arterial models used in the present study are the first 3-D systems reported in literature to incorporate the major arterial vessels that deliver blood from the heart to the cranial and femoral arteries. These models are also the first reported in literature to be used for flow parameter assessment via medical drug delivery and orthostatic postural change studies. The present work addresses the design of the experimental and digital arterial model in addition to the design of measuring tools used to measure hemodynamic parameters. The experimental and digital arterial model analyzed in the present study was developed from patient specific computed tomography angiography (CTA) scans and simplified geometric data. Segments such as the aorta (ascending and descending) and carotid bifurcation arteries of the experimental and digital arterial model was created from online available patient-specific CTA scan data provided by Charite' Clinical and Research Hospital. The cranial and coronary arteries were simplified arterial geometries developed from dimensional specification data used in previous work. For the patient specific geometries, a MATLAB code was written to upload the CTA scans of each artery, calculate the centroids, and produce surface splines at each discrete cross section along the lumen centerline to create the patient specific arterial geometries. The MATLAB code worked in conjunction with computer aided software (CAD) Solidworks to produce solid models of the patient specific geometries and united them with the simplified geometries to produce the full arterial model (CAD model). The CAD model was also used as a blueprint to fabricate the experimental model which was used for flow visualization via particle imaging velocimetry (PIV) and postural change studies. A custom pulse duplicator (pulsatile pump) was also designed and developed for the present work. The pulse duplicator is capable of producing patient-specific volumetric waveforms for inlet flow to the experimental arterial model. A simple fluid structure interaction (FSI) study was also conducted via optical techniques to establish the magnitude of vessel diameter change due to the pulsatile flow. A medical drug delivery (dye dispersion and tracing) case was simulated via a dye being dispersed into the pulsatile flow stream to measure the transit time of the dye front. Pressure waveforms for diseased cases (hypertension & stenotic cases) were also obtained from the experimental arterial model during postural changes from bedrest (0°) to upright position (90°). The postural changes were simulated via attaching the experimental model to a tile table the can transition from 0° to 90°. The PIV results obtained from the experimental model provided parametric data such as velocity and wall shear stress data. The medical drug delivery simulations (experimental and numerical) studies produce time dependent data which is useful for predicting flow trajectory and transit time of medical drug dispersion. In the case of postural change studies, pressure waveforms were obtained from the common carotid artery and the femoral sections to yield pressure difference data useful for orthostatic hypotension analysis. Flow parametric data such as vorticity (flow reversal), wall shear stress, normal stress, and medical drug transit data was also obtained from the digital arterial model CFD simulations. Although the present work is preliminary work, the experimental and digital models proves to be useful in providing flow parametric data of interest such as: (1) normal stress which is useful for predicting the magnitude of forces which could promote arterial rupture or dislodging of medical implants, (2) wall shear stress which is useful for analyzing the magnitude of drug transport at the arterial wall, (3) vorticity which is useful for predicting the magnitude of flow reversal, and (4) arterial compliance in the case of the experimental model which could be useful in the efforts of developing FSI numerical simulations that incorporates compliance which realistically models the flow in the arterial system.

Microscale surface modifications for heat transfer enhancement.

In this experimental study, two surface modification techniques were investigated for their effect on heat transfer enhancement. One of the methods employed the particle (grit) blasting to create microscale indentations, while the other used plasma spray coating to create microscale protrusions on Al 6061 (aluminum alloy 6061) samples. The test surfaces were characterized using scanning electron microscopy (SEM) and confocal scanning laser microscopy. Because of the surface modifications, the actual surface area was increased up to 2.8× compared to the projected base area, and the arithmetic mean roughness value (Ra) was determined to vary from 0.3 µm for the reference smooth surface to 19.5 µm for the modified surfaces. Selected samples with modified surfaces along with the reference smooth surface were then evaluated for their heat transfer performance in spray cooling tests. The cooling system had vapor-atomizing nozzles and used anhydrous ammonia as the coolant in order to achieve heat fluxes up to 500 W/cm(2) representing a thermal management setting for high power systems. Experimental results showed that the microscale surface modifications enhanced heat transfer coefficients up to 76% at 500 W/cm(2) compared to the smooth surface and demonstrated the benefits of these practical surface modification techniques to enhance two-phase heat transfer process.