Unveiling advanced mechanisms of inhalable drug aerosol dynamics using computational fluid dynamics and discrete element method

Abstract

Capsule-based dry powder inhalers (DPIs) are widely used to treat chronic obstructive pulmonary disease (COPD) by delivering active pharmaceutical ingredients (APIs) via inhalation into human respiratory systems. Previous research has shown that the actuation flow rate, aerodynamic particle size distribution (APSD), and particle shape of lactose carriers are factors that can influence the particle deposition patterns in human respiratory systems. Understanding the dynamics of APIs transport in DPIs and airways can provide significant value for the design optimization of DPIs and particle shapes to enhance the delivery of APIs to the designated lung sites, i.e., small airways. Thus, it is necessary to investigate how to modulate the above-mentioned factors to increase the delivery efficacy to small airways and enhance the therapeutic effect to treat COPD. Compared with in vitro and in vivo methods, computational fluid-particle dynamics (CFPD) models allow researchers to study the transport dynamics of inhalable therapeutic dry powders in both DPI and human respiratory systems. However, existing CFPD models neglect particle-particle interactions, and most existing airway models lack peripheral lung airway and neglect the airway deformation kinematics. Such deficiencies can lead to inaccurate predictions of particle transport and deposition. This study developed a one-way coupled computational fluid dynamics (CFD) and discrete element method (DEM) model to simulate the particle-particle and particle-device interactions, and the transport of API-carrier dry powder mixtures with different shapes of carriers in a DPI flow channel. The influence of actuation flow rate (30 to 90 L/min) and particle shape (aspect ratio equals 1, 5, and 10) on lactose carrier dynamics in a representative DPI, i.e., SpirivaTM HandihalerTM, has been investigated. Subsequently, an elastic truncated whole-lung model has also been developed to predict particle transport and deposition from mouth to alveoli, with disease-specific airway deformation kinematics. Numerical results indicate that 90 L/min actuation flow rate generates the highest delivery efficiency of Handihaler, as approximately 26% API reaches the deep lung region. The elastic truncated whole-lung modeling results show that noticeable differences of predictions between static and elastic lung models can be found, which demonstrates the necessity to model airway deformation kinematics in virtual lung models.