Hardware

The Dry Powder Inhaler Design and Air Flow Modeling

We used both SOLIDWORKS Premium 2022 SP1.0 software in designing the dry powder inhaler (DPI) and SOLIDWORKS Flow Simulation add-in as a computational fluid dynamics (CFD) software to model the air flow. We were capable of designing a (DPI) to optimize our particles delivery and airflow.

Our innovative dry powder inhaler design and models are important to:

1. Enhance both the airflow and the particles delivery to the lung by targeting airways effectively

2. Improve dispersion.

3. Improve particle deaggregation by modeling complex interactions like collisions and turbulence for higher fine particle fractions.

4. Make Patient-specific adaptations, which accounts for variability in age and severity.

5. Reduce airflow resistance, making it easier for patients to inhale and decrease both the bacterial remnants in the inhaler and the initial dose (D1).

6. Optimizing the inspiratory flow rate and maximizing its mean (standard deviation) peak according to the mild deviation in the age.

7. Validate and reduce testing, by correlating predictions with data to minimize human trials and speed up development.

8. Reduce expenses on physical prototypes by CFD simulating, which reduces the required number of prototypes to do experiments.

We started by creating every part individually as shown:

Mouthpiece

The mouthpiece part is designed to give the optimum volume flow rate at the outlet. The length of the mouthpiece was studied, as this parameter determines the level of flow development through the mouthpiece. Undeveloped flow can contain regions of high velocity that can enhance throat impaction upon inhalation. It is believed that more development of the flow through the mouthpiece leads to a more uniform flow profile at the device exit. This uniform profile reduces the regions of high velocity, potentially reducing throat impaction and improving overall inhaler performance. So we changed the length of the mouthpiece many times until getting the optimum velocity and volume flow rate at the outlet.

at the beginning of the mouthpiece the air flows through a full grid, which gives the least percentage of particles impacting on the mouthpiece of about 22%, which is far less than other possible cases of grid design, which give impaction percentage varying from ~57% to ~88%.

The outlet internal diameter we chose is optimum for the mass median aerodynamic diameter (MMAD) of 3.0 µm, which is well-suited for targeting the lower airways and potentially the alveolar region, which is desired for therapeutic effects in asthma and also for probiotic delivery to the lung periphery.

Fig(?):

Part 1. Mouthpiece

Capsule Chamber

The capsule chamber part - as shown above - has 2 inlets tangential to it’s inner diameter to prevent or reduce any possible chance to generate a flow separation or a recirculating flow region as much as possible. It is designed to get the optimum inlet volume flow rate and inlet velocity, in order to reduce impaction on the grid and the mouthpiece afterwards.

The capsule chamber is designed to host a HPMC capsule of size 3, so we created the whole assembly design according to that size. The HPMC capsule requires small puncture force relatively, compared to other capsule materials (e.g. Gelatine, …..etc).

Saleem et al., looked at the influence of pin number by using DPIs with a different number of punctures (2-pins vs. 8-pins). A significant difference was found, demonstrating that the 2-pin DPI device showed significantly lower MMAD. Another study conducted by Torrisi et al., comparing the puncture force needed in the case of (2 sets of 4-pins vs. 2 single pins DPI) with gelatine and HPMC capsules, showing – as expected – that 2 single pins need significantly less penetration force. So The capsule will be pierced by 2 piercing buttons with a 0.5 mm pierced aperture, which requires more E_dispersion, having a larger FPF less than 5 µm in diameter (FPF > 5 µm), while having a lower MMAD and ED.

Fig(?):

Part 2. Capsule Chamber

Piercing Buttons

The piercing buttons are designed to have 1 piercing pin each, with a 0.5 mm pierced aperture and travelling length of 3 mm to pierce the capsule with total number of 2 apertures. Each piercing button has 4 bias springs with length of 10 mm and diameter of 2 mm. The bias springs’ mission is to return the piercing buttons to their initial positions.

Fig(?):

Part 3, 4. Piercing Buttons (2 parts)

Cover

The cover is simply used to cover the mouthpiece to prevent any possible pollution to it, protecting the inhaler from dirt, dust and moisture, keeping it clean and dry when not in use.

Fig(?):

Part 5. Cover

Cap

The cap is responsible of covering the bottom of the capsule chamber to protect it from dust also, as the capsule chamber has lots of internal cuts to reduce it’s weight by removing the unneeded material from it, making it economical for production.

Fig(?):

Part 6. Cap

DPI Assembly

Then we assembled them:

Fig(?):

interactive DPI assembly or isometric view photo

DPI Assembly testing

Afterwards we tested the assembly for any possible errors or interferences between bodies and checked the clearances to make any modifications if needed.

Fig(?):

Fig 1. Checking for interferences using Interference Detection

After that we have made the engineering drawing of the assembly:

Fig(?):

Fig 2. Assembly Engineering Drawing

After creating the DPI assembly, which models the air flow from the inlets of the capsule chamber going through the mouthpiece to the outlet, we set the boundary conditions.

Fig(?):

Fig 3. Inlets Boundary Conditions

We set the thermodynamic parameters at the inlets to be atmospheric pressure (101325 Pa) and temperature of 20.05 °C. For turbulence parameters we used Turbulent Energy and Dissipation turbulence model (k-ε), we set them for 1 J/kg and 1 W/kg, respectively.

Fig(?):

Fig 4. Outlet Boundary Conditions

We set the volume flow rate at the outlet to be 100 L/min = 0.001667 m^3/s. this number is carefully chosen according to a specific age and severity, adult with type 2 severe asthma (uncontrolled). The mean peak inspiratory flow rate = 106 L/min with standard deviation of 16 L/min.

The volume flow rate though a DPI (Q) is proportional to the square root of the pressure drop (ΔP) the patient develops across it, namely:

Fig(?):

Equation

The constant of proportionality is termed: the device resistance (R).

Figure 5 presents volume flow rate resistance curves for a variety of inhalers, covering a range of pressure drops that are typically achievable by patients (2-6 Kpa).

Patients inhale faster throw low-resistance devices and slower through high-resistance devices because the pressure drops they generate tend to be similar.

Fig(?):

Fig 5. Volume Flow Rate vs. Device Resistance for several commercial DPIs

The wall conditions set is adiabatic wall with 0.0 µm wall roughness (smooth wall)

After specifying of the boundary conditions and air flow region, we will conduct the mesh, getting ready to run setup.

The mesh below models the DPI and the air flow region:

Fig(?):

Fig 6. Mesh of the DPI

To start the stage of mesh and simulation afterwards, we removed the unneeded parts (cover, cap, piercing buttons), external fillets and chamfers, …etc, in order to reduce the number of mesh cells by removing unnecessary parts, making it faster to run simulation and reducing the unneeded load on the CPU cores and then we have the ability to increase the number of mesh cells in the flow domain, and then to make the results more accurate, in addition to make the results more visualizable when removing these parts.

For the mesh settings we set the level of initial mesh to be 7 out of 7 which is the most accurate level, with a minimum gap size of 1e-06 m, giving a very accurate simulation of air flow and results. We also activated the advanced channel refinement.

After boundary conditions were specified, and the whole volume was meshed. We are ready to run the simulation’s setup.

Setup

We specified some goal plots to get sufficient information in order to see the quality of our design.

The inlets

Fig(?):

Fig 7. Goal Plot at the Inlets

The outlets

Fig(?):

Fig 8. Goal Plot at the Outlet

We did the simulation for 2 cases:

1. DPI device without a capsule

2. DPI device with a capsule

And we will discuss them in details below:

1. DPI device without a capsule

Fig(?):

Fig(?):

These are the average values of pressure, velocity and volume flow rate at the both inlets and the outlet, showing a pressure drop of 4.8236 KPa. Having the volume flow rate, we got the resistance of the DPI device of about 0.02196269565〖 KPa〗^(1⁄2)/L/min, which is less than ~86% of the commercial DPI devices worldwide while not exhibiting the capsule.

We can clearly see that the flow is a turbulent flow, as Reynolds number at the inlet and the outlet, respectively:

Fig(?):

The following are the contours and plots we got as a result of the simulation:

Velocity Magnitude Plots

Fig(?):

Velocity Magnitude Flow Trajectories

Fig(?):

Velocity Magnitude Flow Trajectories

Showing the air trajectories from the inlet through the DPI device to the outlet. The variety of the velocity magnitude is presented as a variety of colors using contour palette.

Fig(?):

Plan Cut Isometric View Velocity Magnitude Contours

Showing an XZ plane cut plot from the internal base of the capsule chamber through the DPI device’s flow domain to the outlet. The variety of the velocity magnitude is presented as a variety of colors using contour palette.

Pressure Magnitude Plots

Fig(?):

Pressure Magnitude Flow Trajectories

Fig(?):

Pressure Magnitude Flow Trajectories

Showing the air trajectories from the inlet through the DPI device to the outlet. The variety of the pressure magnitude is presented as a variety of colors using contour palette.

Fig(?):

Plan Cut Isometric View Pressure Magnitude Contours

Showing an XZ plane cut plot from the internal base of the capsule chamber through the DPI device’s flow domain to the outlet. The variety of the pressure magnitude is presented as a variety of colors using contour palette.

Air Density Magnitude Plots

Fig(?):

Air Density Magnitude Flow Trajectories

Fig(?):

Air Density Magnitude Flow Trajectories

Showing the air trajectories from the inlet through the DPI device to the outlet. The variety of the air density magnitude is presented as a variety of colors using contour palette.

Temperature Magnitude Plots

Fig(?):

Temperature Magnitude Flow Trajectories

Fig(?):

Temperature Magnitude Flow Trajectories

Showing the air trajectories from the inlet through the DPI device to the outlet. The variety of the temperature magnitude is presented as a variety of colors using contour palette.

At the outlet, there are some parameters that vary according to the position on the outlet perimeter. Hereafter, we mentioned some of them:

Vorticity

Fig(?):

Vorticity variation based on position on the outlet perimeter

Pressure

Fig(?):

Pressure variation based on position on the outlet perimeter

Relative Pressure

Fig(?):

Relative pressure variation based on position on the outlet perimeter

Air Density

Fig(?):

Air density variation based on position on the outlet perimeter

Goal plots

Here are some of the goal plots mentioned before in the setup section:

Fig(?):

Average Static Pressure

Fig(?):

Average Dynamic Pressure

Fig(?):

Average Total Pressure

Fig(?):

Average Inlet Velocity

Fig(?):

Average Outlet Velocity

2. DPI device with a capsule

Fig(?):

Fig(?):

These are the average values of pressure, velocity and volume flow rate at the both inlets and the outlet, showing a pressure difference of 4.90896 KPa while exhibiting the capsule. Having the volume flow rate, we got the resistance of the DPI device of about 0.02215617295〖 KPa〗^(1⁄2)/L/min, which is less than ~78% of commercial DPI devices worldwide and somewhat close to the least resistance commercial DPI devices in the world, whose resistance = 0.0125 〖 KPa〗^(1⁄2)/L/min, such as Rotahaler.

We can clearly see that the flow is a turbulent flow here also, as Reynolds number at the inlet and the outlet, respectively:

Fig(?):

Streamlines and vectors

Fig(?):

Right Plane Cut Front View of Streamlines and vectors

Showing the air streamlines and vectors from the inlet through the DPI device to the outlet.

Velocity Magnitude Plots

Fig(?):

Fig(?):

Velocity Magnitude Flow Trajectories

Showing the air trajectories from the inlet through the DPI device to the outlet. The variety of the velocity magnitude is presented as a variety of colors using contour palette.

Fig(?):

Plan Cut Isometric View Velocity Magnitude Contours

Showing an XZ plane cut plot from the internal base of the capsule chamber through the DPI device’s flow domain to the outlet. The variety of the velocity magnitude is presented as a variety of colors using contour palette.

Fig(?):

Right Plane Cut Front View Velocity Magnitude Contours

Showing an XY plane cut plot from the internal base of the capsule chamber through the DPI device’s flow domain to the outlet. The variety of the velocity magnitude is presented as a variety of colors using contour palette.

Pressure Magnitude Plots

Fig(?):

Fig(?):

Pressure Magnitude Flow Trajectories

Showing the air trajectories from the inlet through the DPI device to the outlet. The variety of the pressure magnitude is presented as a variety of colors using contour palette.

Fig(?):

Plan Cut Isometric View Pressure Magnitude Contours

Showing an XZ plane cut plot from the internal base of the capsule chamber through the DPI device’s flow domain to the outlet. The variety of the pressure magnitude is presented as a variety of colors using contour palette.

Fig(?):

Right Plane Cut Front View Pressure Magnitude Contours

Showing an XY plane cut plot from the internal base of the capsule chamber through the DPI device’s flow domain to the outlet. The variety of the pressure magnitude is presented as a variety of colors using contour palette.

Fig(?):

Isometric View Capsule Pressure Magnitude Contours

Showing the variety of the pressure magnitude on the capsule as a variety of colors using contour palette.

Air Density Magnitude Plots

Fig(?):

Fig(?):

Air Density Magnitude Flow Trajectories

Showing the air trajectories from the inlet through the DPI device to the outlet. The variety of the air density magnitude is presented as a variety of colors using contour palette.

Fig(?):

Plan Cut Isometric View Air Density Magnitude Contours

Showing an XZ plane cut plot from the internal base of the capsule chamber through the DPI device’s flow domain to the outlet. The variety of the air density magnitude is presented as a variety of colors using contour palette.

Fig(?):

Right Plane Cut Front View Air Density Magnitude Contours

Showing an XY plane cut plot from the internal base of the capsule chamber through the DPI device’s flow domain to the outlet. The variety of the air density magnitude is presented as a variety of colors using contour palette.

Temperature Magnitude Plots

Fig(?):

Fig(?):

Temperature Magnitude Flow Trajectories

Showing the air trajectories from the inlet through the DPI device to the outlet. The variety of the temperature magnitude is presented as a variety of colors using contour palette.

Fig(?):

Plan Cut Isometric View Temperature Magnitude Contours

Showing an XZ plane cut plot from the internal base of the capsule chamber through the DPI device’s flow domain to the outlet. The variety of the temperature magnitude is presented as a variety of colors using contour palette.

Fig(?):

Right Plane Cut Front View Temperature Magnitude Contours

Showing an XY plane cut plot from the internal base of the capsule chamber through the DPI device’s flow domain to the outlet. The variety of the temperature magnitude is presented as a variety of colors using contour palette.

At the outlet, there are some parameters that vary according to the position on the outlet perimeter. Hereafter, we mentioned some of them:

Vorticity

Fig(?):

Vorticity variation based on position on the outlet perimeter

Pressure

Fig(?):

Pressure variation based on position on the outlet perimeter

Relative Pressure

Fig(?):

Relative pressure variation based on position on the outlet perimeter

Air Density

Fig(?):

Air density variation based on position on the outlet perimeter

Goal plots

Here are some of the goal plots mentioned before in the setup section:

Fig(?):

Velocity

Fig(?):

Static Pressure

Fig(?):

Total Pressure

Fig(?):

Volume Flow Rate

Fig(?):

Absolute Total Enthalpy Rate

References

[1] Clark, A. R., Weers, J. G., & Dhand, R. (2020). The confusing world of dry powder inhalers: It is all about inspiratory pressures, not inspiratory flow rates. Journal of Aerosol Medicine and Pulmonary Drug Delivery, 33(1), 1–11.

[2] Elkins, M. R., Robinson, P., Anderson, S. D., Perry, C. P., Daviskas, E., & Charlton, B. (2014). Inspiratory flows and volumes in subjects with cystic fibrosis using a new dry powder inhaler device. The Open Respiratory Medicine Journal, 8, 1-7.

[3] Benque, B., & Khinast, J. G. (2019). Understanding the motion of hard-shell capsules in dry powder inhalers. International Journal of Pharmaceutics, 118481.

[4] Yang, M. Y., Verschuer, J., Shi, Y., Song, Y., Katsifis, A., Eberl, S., Wong, K., Brannan, J. D., Cai, W., Finlay, W. H., & Chan, H.-K. (2016). The effect of device resistance and inhalation flow rate on the lung deposition of orally inhaled mannitol dry powder. International Journal of Pharmaceutics, 513(1–2), 256–264.

[5] Ke, W.-R., Chang, R. Y. K., Kwok, P. C. L., Tang, P., Chen, L., Chen, D., & Chan, H.-K. (2021). Administration of dry powders during respiratory supports. Annals of Translational Medicine, 9(7), 596.

[6] Kopsch, T., Murnane, D., & Symons, D. (2016). Optimizing the entrainment geometry of a dry powder inhaler: Methodology and preliminary results. Pharmaceutical Research, 33(11), 2668–2679.

[7] Huynh, B. K., Chen, Y., Fletcher, D. F., Young, P., Zhu, B., & Traini, D. (2015). An investigation into the powder release behavior from capsule-based dry powder inhalers. Aerosol Science and Technology, 49(10), 902–911.

[8] Laube, B. L. (2020). The confusing world of dry powder inhalers: It is all about inspiratory pressures, not inspiratory flow rates. Journal of Aerosol Medicine and Pulmonary Drug Delivery, 33(1), 1–11.

[9] Lexmond, A. J., Kruizinga, T. J., Hagedoorn, P., & de Boer, A. H. (2014). Effect of inhaler design variables on paediatric use of dry powder inhalers. PLoS ONE, 9(6), e99304.

[10] Demoly, P., Hagedoorn, P., de Boer, A. H., & Frijlink, H. W. (2014). The clinical relevance of dry powder inhaler performance for drug delivery. Respiratory Medicine, 108(8), 1195–1203.

[11] Mahler, D. A. (2020). The role of inspiratory flow in selection and use of inhaled therapy for patients with chronic obstructive pulmonary disease. Respiratory Medicine, 161, 105857.

[12] 12 / 12 Elkins, M. R., Anderson, S. D., Perry, C. P., Daviskas, E., & Charlton, B. (2014). Inspiratory flows and volumes in subjects with non-CF bronchiectasis using a new dry powder inhaler device. The Open Respiratory Medicine Journal, 8, 8–13.

[13] McEvoy, C., Argula, R., Sahay, S., Shapiro, S., Eagan, C., Hickey, A. J., Smutney, C., Dillon, C., Winkler, T., Davis, B. N., Broderick, M., & Burger, C. (2023). Tyvaso DPI: Drug-device characteristics and patient clinical considerations. Pulmonary Pharmacology & Therapeutics, 102266.

[14] Weers, J. (2022). Suboptimal inspiratory flow rates with passive dry powder inhalers: Big issue or overstated problem? Frontiers in Drug Delivery, 2, 855234.

[15] Cai, X., Dong, J., Milton-McGurk, L., Lee, A., Shen, Z., Chan, H.-K., Kourmatzis, A., & Cheng, S. (2023). Understanding the effects of inhaler resistance on particle deposition behaviour – A computational modelling study. Computers in Biology and Medicine, 166, 107673.

[16] Altman, P., Wehbe, L., Dederichs, J., Guerin, T., Ament, B., Moronta, M. C., Pino, A. V., & Goyal, P. (2018). Comparison of peak inspiratory flow rate via the Breezhaler®, Ellipta® and HandiHaler® dry powder inhalers in patients with moderate to very severe COPD: A randomized cross-over trial. PMC, PMCID: PMC6001060, PMID: 29898702.

[17] Abadelah, M., Chrystyn, H., Bagherisadeghi, G., Abdalla, G., & Larhrib, H. (2018). Study of the emitted dose after two separate inhalations at different inhalation flow rates and volumes and an assessment of aerodynamic characteristics of Indacaterol Onbrez Breezhaler® 150 and 300 μg. AAPS PharmSciTech, 19, 251–261.

[18] Colthorpe, P., Voshaar, T., Kieckbusch, T., Cuoghi, E., & Jauernig, J. (2013). Delivery characteristics of a low-resistance dry-powder inhaler used to deliver the long-acting muscarinic antagonist glycopyrronium. Journal of Drug Assessment, 2(1), 11–16.

[19] Vartiainen, V. A., Lavorini, F., Murphy, A. C., & Rabe, K. F. (2023). High inhaler resistance does not limit successful inspiratory maneuver among patients with asthma or COPD. Expert Review of Clinical Pharmacology, 16(4), 385–393.

[20] Capstick, T., & Clifton, I. (2012). Inhaler technique and training in people with chronic obstructive pulmonary disease and asthma. Expert Review of Respiratory Medicine, 6(1), 91–101.

[21] Capstick, T. G. D., Gudimetla, S., Harris, D. S., Malone, R., & Usmani, O. S. (2024). Demystifying dry powder inhaler resistance with relevance to optimal patient care. Clinical Drug Investigation, 44(2), 109–114.

[22] Leving, M., Wouters, H., de la Hoz, A., & Kocks, J. H. (2021). Impact of PIF, inhalation technique and medication adherence on health status and exacerbations in COPD: Protocol of a real-world observational study (PIFotal COPD Study). Respiratory Research.

[23] Mahler, D. A. (n.d.). Peak inspiratory flow rate as a criterion for dry powder inhaler use in chronic obstructive pulmonary disease. Geisel School of Medicine, Dartmouth College, Hanover, NH; Respiratory Services, Valley Regional Hospital, Claremont, NH.

[24] Haidl, P., Heindl, S., Siemon, K., Bernacka, M., & Cloes, R. M. (2016). Inhalation device requirements for patients' inhalation maneuvers. Respiratory Medicine.

[25] Srichana, T., Martin, G. P., & Marriott, C. (1998). Dry powder inhalers: The influence of device resistance and powder formulation on drug and lactose deposition in vitro. European Journal of Pharmaceutical Sciences, 7(1), 73–80.

[26] de Boer, A. H., Hagedoorn, P., Hoppentocht, M., Buttini, F., Grasmeijer, F., & Frijlink, H. W. (2016). Dry powder inhalation: Past, present and future. Expert Opinion on Drug Delivery, 13(4), 499–512.

Back to top