ABSTRACT
In this work the dynamic flow as well as the particle motion and deposition in a commercial dry powder inhaler, DPI (i.e., Turbuhaler) is described using computational fluid dynamics, CFD. The dynamic flow model presented here is an extension of a steady flow model previously described in Milenkovic et al. (2013). The model integrates CFD simulations for dynamic flow, an Eulerian-fluid/Lagrangian-particle description of particle motion as well as a particle/wall interaction model providing the sticking efficiency of particles colliding with the DPI walls. The dynamic flow is imposed by a time varying outlet pressure and the particle injections into the DPI are assumed to occur instantaneously and follow a prescribed particle size distribution, PSD. The total particle deposition and the production of fine particles in the DPI are determined for different peak inspiratory flow rates, PIFR, flow increase rates, FIR, and particle injection times. The simulation results for particle deposition are found to agree well with available experimental data for different values of PIFR and FIR. The predicted values of fine particle fraction are in agreement with available experimental results when the mean size of the injected PSD is taken to depend on the PIFR.
Subject(s)
Computer Simulation , Dry Powder Inhalers , Hydrodynamics , Models, Theoretical , Administration, Inhalation , Aerosols , Equipment Design , Particle Size , Powders , Time FactorsABSTRACT
In this work the steady-state flow in a commercial dry powder inhaler device, DPI (i.e., Turbuhaler) is described using computational fluid dynamics. The Navier-Stokes equations are solved using commercial CFD software considering different flow models, i.e., laminar, k-ε, k-ε RNG, and k-ω SST as well as large Eddy simulation. Particle motion and deposition are described using a Eulerian-fluid/Lagrangian-particle approach. Particle collisions with the DPI walls are taken to result in deposition when the normal collision velocity is less than a critical capture velocity. Flow and particle deposition, for a range of mouthpiece pressure drops (i.e., 800-8800 Pa), as well as particle sizes corresponding to single particles and aggregates (i.e., 0.5-20 µm), are examined. The total volumetric outflow rate, the overall particle deposition as well as the spatial distribution of deposition sites in the DPI are determined. The transitional k-ω SST model for turbulent flow was found to produce results most similar to a reference solution obtained with LES, as well as experimental results for the pressure drop in the DPI. Overall, the simulation results are found to be in agreement with the available experimental data for local and total particle deposition.
