A Two-Species Finite Volume Scalar Model for Modeling the Diffusion of Poly(lactic-co-glycolic acid) into a Coronary Arterial Wall from a Single Half-Embedded Drug Eluting Stent Strut

This paper outlines the methodology and results for a two-species finite volume scalar computational drug transport model developed for simulating the mass transport of Poly(lacticco-glycolic acid (PLGA)) from a half-embedded single strut implanted in a coronary arterial vessel wall. The mathematical drug transport model incorporates the convection-diffusion equation in scalar form (dimensionless) with a two-species (free-drug and bound-drug) mass transport setup, including reversible equilibrium reaction source terms for the free and bound-drug states to account for the pharmaco-kinetic reactions in the arterial wall. The relative reaction rates of the added source terms control the interconversion of the drug between the free and bound-drug states. The model is solved by a 2D finite-volume method for discretizing and solving the free and bound drug transport equations with anisotropic vascular drug diffusivities. This model is an improvement over previously developed models using the finite-difference and finite element method. A dimensionless characteristic scaling pre-analysis was conducted a priori to evaluate the significance of implementing the reaction source terms in the transport equations. This paper reports the findings of an investigation of the interstitial flow profile into the arterial wall and the free and bound drug diffusion profiles with a parametric study of varying the polymer drug concentration (low and high), tortuosity, porosity, and Peclet and DamKöhler numbers over the course of 400 h (16.67 days). The results also reveal how a single species drug delivery model that neglects both a reversible binding reaction source term and the porosity and tortuosity of the arterial wall cannot accurately predict the distribution of both the free and bound drug.