1
|
Jing H, Ge H, Tang H, Farnoud A, Saidul Islam M, Wang L, Wang C, Cui X. Assessing airflow unsteadiness in the human respiratory tract under different expiration conditions. J Biomech 2024; 162:111910. [PMID: 38154261 DOI: 10.1016/j.jbiomech.2023.111910] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/11/2023] [Revised: 11/20/2023] [Accepted: 12/14/2023] [Indexed: 12/30/2023]
Abstract
To enhance the understanding of airflow characteristics in the human respiratory system, the expiratory airflow in a human respiratory tract model was simulated using large eddy simulation and dynamic mesh under different expiration conditions aligned with clinically measured data. The airflow unsteadiness was quantitatively assessed using power spectral density (PSD) and spectral entropy (SE). The following findings were obtained: (1) The airflow is highly turbulent in the mouth-pharynx region during expiration, with its dynamic characteristics being influenced by both the transient expiration flow pattern at mouth piece and the glottis motion. (2) PSD analysis reveals that the expiratory airflow is very unsteady, exhibiting a broad-band attenuation spectrum in the pharynx-trachea region. When only transient expiration or glottis motion is considered, the PSD spectrum changes slightly. When both are ignored, however, the change is significant, with the peak frequency reduced to 10% of the real expiration condition. (3) SE analysis indicates that the airflow transitions into turbulence in the trachea, and there may be multiple transitions in the region of soft palate. The transient expiration or glottis motion alone increases turbulence intensity by 2%-15%, while ignoring both reduces turbulence intensity by 10%-20%. This study implies that turbulence characteristics can be significantly different under different expiratory conditions, and therefore it is necessary to determine the expiratory flow characteristics using clinically measured expiratory data.
Collapse
Affiliation(s)
- Hao Jing
- School of Aerospace Engineering, Huazhong University of Science and Technology, Wuhan, China
| | - Haiwen Ge
- Research Center for Intelligent Supercomputing, Zhejiang Laboratory, Hangzhou, China
| | - Hui Tang
- Department of Mechanical Engineering, The Hong Kong Polytechnic University, Hong Kong, China
| | - Ali Farnoud
- Helmholtz Zentrum München-German Research Center for Environmental Health, Institute of Computational Biology, Neuherberg, Germany
| | - Mohammad Saidul Islam
- School of Mechanical and Mechatronic Engineering, University of Technology Sydney, Sydney, Australia
| | - Li Wang
- School of Environmental Science and Engineering, Tianjin University, Tianjin, China
| | - Chenglei Wang
- Department of Mechanical Engineering, The Hong Kong Polytechnic University, Hong Kong, China
| | - Xinguang Cui
- School of Aerospace Engineering, Huazhong University of Science and Technology, Wuhan, China.
| |
Collapse
|
2
|
Takahashi K, Toyama H, Ejima Y, Yang J, Kikuchi K, Ishikawa T, Yamauchi M. Endotracheal tube, by the venturi effect, reduces the efficacy of increasing inlet pressure in improving pendelluft. PLoS One 2023; 18:e0291319. [PMID: 37708106 PMCID: PMC10501657 DOI: 10.1371/journal.pone.0291319] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/16/2023] [Accepted: 08/27/2023] [Indexed: 09/16/2023] Open
Abstract
In mechanically ventilated severe acute respiratory distress syndrome patients, spontaneous inspiratory effort generates more negative pressure in the dorsal lung than in the ventral lung. The airflow caused by this pressure difference is called pendelluft, which is a possible mechanisms of patient self-inflicted lung injury. This study aimed to use computer simulation to understand how the endotracheal tube and insufficient ventilatory support contribute to pendelluft. We established two models. In the invasive model, an endotracheal tube was connected to the tracheobronchial tree with 34 outlets grouped into six locations: the right and left upper, lower, and middle lobes. In the non-invasive model, the upper airway, including the glottis, was connected to the tracheobronchial tree. To recreate the inspiratory effort of acute respiratory distress syndrome patients, the lower lobe pressure was set at -13 cmH2O, while the upper and middle lobe pressure was set at -6.4 cmH2O. The inlet pressure was set from 10 to 30 cmH2O to recreate ventilatory support. Using the finite volume method, the total flow rates through each model and toward each lobe were calculated. The invasive model had half the total flow rate of the non-invasive model (1.92 L/s versus 3.73 L/s under 10 cmH2O, respectively). More pendelluft (gas flow into the model from the outlets) was observed in the invasive model than in the non-invasive model. The inlet pressure increase from 10 to 30 cmH2O decreased pendelluft by 11% and 29% in the invasive and non-invasive models, respectively. In the invasive model, a faster jet flowed from the tip of the endotracheal tube toward the lower lobes, consequently entraining gas from the upper and middle lobes. Increasing ventilatory support intensifies the jet from the endotracheal tube, causing a venturi effect at the bifurcation in the tracheobronchial tree. Clinically acceptable ventilatory support cannot completely prevent pendelluft.
Collapse
Affiliation(s)
- Kazuhiro Takahashi
- Anesthesiology and Perioperative Medicine, Tohoku University Graduate School of Medicine, Sendai, Japan
| | - Hiroaki Toyama
- Anesthesiology and Perioperative Medicine, Tohoku University Graduate School of Medicine, Sendai, Japan
| | - Yutaka Ejima
- Division of Surgical Center and Supply, Sterilization, Tohoku University Hospital, Sendai, Japan
| | - Jinyou Yang
- Department of Biophysics, School of Intelligent Medicine, China Medical University, Shenyang, China
| | - Kenji Kikuchi
- Department of Finemechanics, Graduate School of Engineering, Tohoku University, Sendai, Japan
| | - Takuji Ishikawa
- Graduate School of Biomedical Engineering, Tohoku University, Sendai, Japan
| | - Masanori Yamauchi
- Anesthesiology and Perioperative Medicine, Tohoku University Graduate School of Medicine, Sendai, Japan
| |
Collapse
|
3
|
Jiang L, Han H, Yang J, Fang R, Xin Y, Chen Q, Yao J, Li Z. Do the stent blind insertion into the main portal vein (MPV) and stent diameter influence the surgical outcome of the transjugular intrahepatic portosystemic shunt (TIPS)? Comput Biol Med 2023; 164:107306. [PMID: 37542920 DOI: 10.1016/j.compbiomed.2023.107306] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/20/2023] [Revised: 06/24/2023] [Accepted: 07/28/2023] [Indexed: 08/07/2023]
Abstract
Transjugular intrahepatic portosystemic shunt (TIPS) surgery is a clinical intervention to treat portal hypertension (PH) by deploying a covered stent to establish a shunt path for the portal vein (PV) system, and proper surgical strategy is of great importance to balance the shunt effect and the risk of complications. To understand the clinical strategies of the stent blind insertion and stent selection in clinic, this study investigated the effects of varying stent insertion positions and diameters on the PV hemodynamics and the shunt effect by computational fluid dynamics (CFD) analysis of five post-TIPS subjects. The results showed that the successful TIPS surgeries of the five PH subjects were confirmed by quantifying their pressure drops. The stent insertion positions at the main portal vein (MPV) slightly affected the clinically concerned hemodynamic indexes (i.e., MPV pressure, stent-outlet velocity) and the shunt index (SI). This indicated that the position of the stent going into the MPV may not need to be deliberately selected. Moreover, the covered stents with 6 mm and 8 mm diameters slightly influenced the hemodynamics as well, but the large-diameter stent better improved the shunt effect compared to the small-diameter one. Despite this, the 6 mm stent was suggested thanks to the higher risk of the hepatic encephalopathy (HE) observed in clinic, which indicated the excessive shunt of the 8 mm stent. The current work revealed the effects of different TIPS strategies on the surgical outcome, and could be useful for potential clinical practices.
Collapse
Affiliation(s)
- Lingyu Jiang
- Biomechanics Laboratory, School of Biological Science & Medical Engineering, Southeast University, Nanjing, 210096, China
| | - Hao Han
- Department of Ultrasound Medicine, Affiliated Drum Tower Hospital of Nanjing University Medical School, Nanjing, 210008, Jiangsu, People's Republic of China; Medical Imaging Center, Affiliated Drum Tower Hospital of Nanjing University Medical School, Nanjing, 210008, Jiangsu, People's Republic of China
| | - Jian Yang
- Department of Ultrasound Medicine, Affiliated Drum Tower Hospital of Nanjing University Medical School, Nanjing, 210008, Jiangsu, People's Republic of China; Medical Imaging Center, Affiliated Drum Tower Hospital of Nanjing University Medical School, Nanjing, 210008, Jiangsu, People's Republic of China
| | - Runxin Fang
- Biomechanics Laboratory, School of Biological Science & Medical Engineering, Southeast University, Nanjing, 210096, China
| | - Ying Xin
- Department of Ultrasound Medicine, Affiliated Drum Tower Hospital of Nanjing University Medical School, Nanjing, 210008, Jiangsu, People's Republic of China; Medical Imaging Center, Affiliated Drum Tower Hospital of Nanjing University Medical School, Nanjing, 210008, Jiangsu, People's Republic of China
| | - Qiang Chen
- Biomechanics Laboratory, School of Biological Science & Medical Engineering, Southeast University, Nanjing, 210096, China.
| | - Jing Yao
- Department of Ultrasound Medicine, Affiliated Drum Tower Hospital of Nanjing University Medical School, Nanjing, 210008, Jiangsu, People's Republic of China; Medical Imaging Center, Affiliated Drum Tower Hospital of Nanjing University Medical School, Nanjing, 210008, Jiangsu, People's Republic of China.
| | - Zhiyong Li
- Biomechanics Laboratory, School of Biological Science & Medical Engineering, Southeast University, Nanjing, 210096, China; School of Mechanical, Medical and Process Engineering, Queensland University of Technology (QUT), Brisbane, QLD4001, Australia; Faculty of Sports Science, Ningbo University, Ningbo, 315211, China
| |
Collapse
|
4
|
Ou C, Hang J, Hua J, Li Y, Deng Q, Zhao B, Ling H. Particle Deposition in Large-Scale Human Tracheobronchial Airways Predicted by Single-Path Modelling. INTERNATIONAL JOURNAL OF ENVIRONMENTAL RESEARCH AND PUBLIC HEALTH 2023; 20:4583. [PMID: 36901592 PMCID: PMC10002109 DOI: 10.3390/ijerph20054583] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 01/29/2023] [Revised: 02/22/2023] [Accepted: 02/28/2023] [Indexed: 06/18/2023]
Abstract
The health effects of particles are directly related to their deposition patterns (deposition site and amount) in human airways. However, estimating the particle trajectory in a large-scale human lung airway model is still a challenge. In this work, a truncated single-path, large-scale human airway model (G3-G10) with a stochastically coupled boundary method were employed to investigate the particle trajectory and the roles of their deposition mechanisms. The deposition patterns of particles with diameters (dp) of 1-10 μm are investigated under various inlet Reynolds numbers (Re = 100-2000). Inertial impaction, gravitational sedimentation, and combined mechanism were considered. With the increasing airway generations, the deposition of smaller particles (dp < 4 μm) increased due to gravitational sedimentation, while that of larger particles decreased due to inertial impaction. The obtained formulas of Stokes number and Re can predict the deposition efficiency due to the combined mechanism in the present model, and the prediction can be used to assess the dose-effect of atmospheric aerosols on the human body. Diseases in deeper generations are mainly attributed to the deposition of smaller particles under lower inhalation rates, while diseases at the proximal generations mainly result from the deposition of larger particles under higher inhalation rates.
Collapse
Affiliation(s)
- Cuiyun Ou
- School of Atmospheric Sciences, Sun Yat-sen University, and Southern Marine Science and Engineering Guangdong Laboratory (Zhuhai), Zhuhai 519082, China
| | - Jian Hang
- School of Atmospheric Sciences, Sun Yat-sen University, and Southern Marine Science and Engineering Guangdong Laboratory (Zhuhai), Zhuhai 519082, China
| | - Jiajia Hua
- China Meteorological Administration Xiong’an Atmospheric Boundary Layer Key Laboratory, Baoding 071800, China
| | - Yuguo Li
- Department of Mechanical Engineering, The University of Hong Kong, Hong Kong, China
| | - Qihong Deng
- College of Public Health, Zhengzhou University, Zhengzhou 450001, China
| | - Bo Zhao
- China Meteorological Administration Xiong’an Atmospheric Boundary Layer Key Laboratory, Baoding 071800, China
| | - Hong Ling
- School of Atmospheric Sciences, Sun Yat-sen University, and Southern Marine Science and Engineering Guangdong Laboratory (Zhuhai), Zhuhai 519082, China
| |
Collapse
|
5
|
Spasov G, Rossi R, Vanossi A, Cottini C, Benassi A. A critical analysis of the CFD-DEM simulation of pharmaceutical aerosols deposition in extra-thoracic airways. Int J Pharm 2022; 629:122331. [DOI: 10.1016/j.ijpharm.2022.122331] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/27/2022] [Revised: 10/17/2022] [Accepted: 10/18/2022] [Indexed: 11/07/2022]
|
6
|
CFD Simulations of Respiratory Airflow in Human Upper Airways Response to Walking and Running for Oral Breathing Condition. Heliyon 2022; 8:e10039. [PMID: 35982840 PMCID: PMC9379579 DOI: 10.1016/j.heliyon.2022.e10039] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/07/2021] [Revised: 01/08/2022] [Accepted: 07/18/2022] [Indexed: 11/22/2022] Open
Abstract
Walking and running are common types of physical activities that people do in day to day living, to improve health and physical fitness or for recreation. During a physical activity, rate and depth of breathing increase because working muscles need extra oxygen in order to produce energy. In this study, computational fluid dynamics (CFD) simulations were used to investigate respiratory airflow dynamics in human upper airways response to walking and running for oral breathing. The numerical simulations were done in a realistic CT-scan airway model using ANAYS Fluent 19.0 software. Flow fields were analysed numerically and flow patterns were investigated in the airway model during inspiration and expiration response to walking and running. The axial velocity distributions and secondary flow patterns for the two respiratory phases were analysed response to the two physical activities at different cross-sections of the airway model. The maximum velocity, wall pressure and wall shear stress values for running were respectively 3.2, 9.4 and 5.9 times higher than that of walking during inspiration. The mixing of flow streamlines was observed to be higher during running than walking because of more significant turbulence. More skewed flows at airway curvatures were observed at inspiration than expiration. The results of this study supported the fact that running is a more intense activity than walking from a respiratory dynamics point of view.
Collapse
|
7
|
Lower Inspiratory Breathing Depth Enhances Pulmonary Delivery Efficiency of ProAir Sprays. Pharmaceuticals (Basel) 2022; 15:ph15060706. [PMID: 35745624 PMCID: PMC9227885 DOI: 10.3390/ph15060706] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/27/2022] [Revised: 05/31/2022] [Accepted: 06/01/2022] [Indexed: 12/17/2022] Open
Abstract
Effective pulmonary drug delivery using a metered-dose inhaler (MDI) requires a match between the MDI sprays, the patient’s breathing, and respiratory physiology. Different inhalers generate aerosols with distinct aerosol sizes and speeds, which require specific breathing coordination to achieve optimized delivery efficiency. Inability to perform the instructed breathing maneuver is one of the frequently reported issues during MDI applications; however, their effects on MDI dosimetry are unclear. The objective of this study is to systemically evaluate the effects of breathing depths on regional deposition in the respiratory tract using a ProAir-HFA inhaler. An integrated inhaler mouth-throat-lung geometry model was developed that extends to the ninth bifurcation (G9). Large-eddy simulation (LES) was used to compute the airflow dynamics due to concurrent inhalation and orifice flows. The discrete-phase Lagrangian model was used to track droplet motions. Experimental measurements of ProAir spray droplet sizes and speeds were used as initial and boundary conditions to develop the computational model for ProAir-pulmonary drug delivery. The time-varying spray plume from a ProAir-HFA inhaler into the open air was visualized using a high-speed imaging system and was further used to validate the computational model. The inhalation dosimetry of ProAir spray droplets in the respiratory tract was compared among five breathing depths on a regional, sub-regional, and local basis. The results show remarkable differences in airflow dynamics within the MDI mouthpiece and the droplet deposition distribution in the oral cavity. The inhalation depth had a positive relationship with the deposition in the mouth and a negative relationship with the deposition in the five lobes beyond G9 (small airways). The highest delivery efficiency to small airways was highest at 15 L/min and declined with an increasing inhalation depth. The drug loss inside the MDI was maximal at 45–60 L/min. Comparisons to previous experimental and numerical studies revealed a high dosimetry sensitivity to the inhaler type and patient breathing condition. Considering the appropriate inhalation waveform, spray actuation time, and spray properties (size and velocity) is essential to accurately predict inhalation dosimetry from MDIs. The results highlight the importance of personalized inhalation therapy to match the patient’s breathing patterns for optimal delivery efficiencies. Further complimentary in vitro or in vivo experiments are needed to validate the enhanced pulmonary delivery at 15 L/min.
Collapse
|
8
|
Thomas ML, Longest PW. Evaluation of the Polyhedral Mesh Style for Predicting Aerosol Deposition in Representative Models of the Conducting Airways. JOURNAL OF AEROSOL SCIENCE 2022; 159:105851. [PMID: 34658403 PMCID: PMC8513711 DOI: 10.1016/j.jaerosci.2021.105851] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/13/2023]
Abstract
A critical factor affecting the accuracy of Computational Fluid Dynamic (CFD) simulations and the time required to conduct them is construction of the computational mesh. This study aimed to evaluate the relatively new polyhedral mesh style for simulating aerosol deposition in the upper conducting airways compared with established meshing techniques and experimental data. Hexahedral and polyhedral mesh solutions were compared in two benchmark geometries: 1) a 90°-bend with flow characteristics similar to the extrathoracic airways of an adolescent child, and 2) a double bifurcation representing bifurcations B3-B5 in an adult. Both 4-block and 5-block hexahedral meshes were used in the 90°-bend to capture the potential of fully-structured hexahedral meshes. In the 90°-bend, polyhedral elements matched polydisperse in vitro deposition data with 20% relative error (RE; averaged across the particle sizes considered), which is an improvement on the accuracy of the 4-block hexahedral mesh (35% RE) and is similar to the accuracy of the 5-block hexahedral mesh (19% RE). In the double bifurcation, deposition fraction relative differences evaluated between polyhedral and hexahedral meshes ranged from 0.3% to 28.6% for the different particle sizes assessed, which is an order of magnitude improvement compared with previous studies that considered hexahedral vs. hybrid tetrahedral-prism meshes for the same flow field. Solution convergence time with polyhedral elements was found to be 50% to 140% higher than with hexahedral meshes of comparable size. While application dependent, the increase in simulation time observed with polyhedral meshes will likely be outweighed by the ease and convenience of polyhedral mesh construction. It was concluded that the polyhedral mesh style, with sufficient resolution especially near the walls, is an excellent alternative to the highly regarded hexahedral mesh style for predicting upper airway aerosol transport and deposition and provides a powerful new tool in the assessment of respiratory aerosol dosimetry.
Collapse
Affiliation(s)
- Morgan L. Thomas
- Department of Mechanical and Nuclear Engineering Virginia Commonwealth University, Richmond, VA
| | - P. Worth Longest
- Department of Mechanical and Nuclear Engineering Virginia Commonwealth University, Richmond, VA
- Department of Pharmaceutics Virginia Commonwealth University, Richmond, VA
| |
Collapse
|
9
|
Nguyen VL, Obara H. Investigation of vessel occlusion during cell seeding process. Biomech Model Mechanobiol 2021; 20:2437-2450. [PMID: 34480225 DOI: 10.1007/s10237-021-01517-6] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/27/2021] [Accepted: 08/25/2021] [Indexed: 11/26/2022]
Abstract
The seeding of cells into an organ is an important step in cell therapy because the final functional properties of the organ are related to the initial cell distribution throughout the organ. However, vessel occlusion is a serious problem that prevents uniform distribution of the cells in the entire organ. Understanding the mechanism of vessel occlusion can help optimize the seeding process. In this study, the vessel occlusion phenomenon under perfusion conditions during cell seeding was investigated. First, we applied a microfluidic system that enabled the observation of the occlusion events during injection. Second, we applied a multiphase numerical model that can describe the cell-cell interactions and cell-fluid interactions to investigate the vessel occlusion phenomenon during the seeding process. In particular, the effects of cell concentration and flow rate were investigated. The results indicate the importance of cell-cell interactions and cell-vessel interactions for the occurrence of vessel occlusion. In addition, it is found that the probability of occurrence of vessel occlusion increases with the increase in cell concentration and decrease in flow rate. The simulation model can help determine the optimum parameters to enhance cell seeding efficiency.
Collapse
Affiliation(s)
- Van Lap Nguyen
- Department of Mechanical Systems Engineering, Tokyo Metropolitan University, 1-1 Minami-Osawa, Hachioji, Tokyo, 192-0397, Japan.
- Faculty of Mechanical Engineering, Thuyloi University, 175 Tay Son, Dong Da, Hanoi, Vietnam.
| | - Hiromichi Obara
- Department of Mechanical Systems Engineering, Tokyo Metropolitan University, 1-1 Minami-Osawa, Hachioji, Tokyo, 192-0397, Japan
| |
Collapse
|
10
|
Huang F, Zhu Q, Zhou X, Gou D, Yu J, Li R, Tong Z, Yang R. Role of CFD based in silico modelling in establishing an in vitro-in vivo correlation of aerosol deposition in the respiratory tract. Adv Drug Deliv Rev 2021; 170:369-385. [PMID: 32971228 DOI: 10.1016/j.addr.2020.09.007] [Citation(s) in RCA: 31] [Impact Index Per Article: 10.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/08/2020] [Revised: 08/02/2020] [Accepted: 09/17/2020] [Indexed: 12/14/2022]
Abstract
Effective evaluation and prediction of aerosol transport deposition in the human respiratory tracts are critical to aerosol drug delivery and evaluation of inhalation products. Establishment of an in vitro-in vivo correlation (IVIVC) requires the understanding of flow and aerosol behaviour and underlying mechanisms at the microscopic scale. The achievement of the aim can be facilitated via computational fluid dynamics (CFD) based in silico modelling which treats the aerosol delivery as a two-phase flow. CFD modelling research, in particular coupling with discrete phase model (DPM) and discrete element method (DEM) approaches, has been rapidly developed in the past two decades. This paper reviews the recent development in this area. The paper covers the following aspects: geometric models of the respiratory tract, CFD turbulence models for gas phase and its coupling with DPM/DEM for aerosols, and CFD investigation of the effects of key factors associated with geometric variations, flow and powder characteristics. The review showed that in silico study based on CFD models can effectively evaluate and predict aerosol deposition pattern in human respiratory tracts. The review concludes with recommendations on future research to improve in silico prediction to achieve better IVIVC.
Collapse
Affiliation(s)
- Fen Huang
- School of Energy and Environment, Southeast University, Nanjing 210096, China; Department of Chemical Engineering, Monash University, Clayton, Vic 3800, Australia
| | - Qixuan Zhu
- School of Materials Science and Engineering, University of New South Wales, Sydney, NSW 2052, Australia
| | - Xudong Zhou
- Department of Chemical Engineering, Monash University, Clayton, Vic 3800, Australia
| | - Dazhao Gou
- School of Materials Science and Engineering, University of New South Wales, Sydney, NSW 2052, Australia
| | - Jiaqi Yu
- Institute for Process Modelling and Optimization, JITRI, Suzhou 215000, China
| | - Renjie Li
- Institute for Process Modelling and Optimization, JITRI, Suzhou 215000, China
| | - Zhenbo Tong
- School of Energy and Environment, Southeast University, Nanjing 210096, China.
| | - Runyu Yang
- School of Materials Science and Engineering, University of New South Wales, Sydney, NSW 2052, Australia.
| |
Collapse
|
11
|
Talebizadehsardari P, Rahimzadeh H, Ahmadi G, Inthavong K, Keshtkar MM, Moghimi M. Nano-particle deposition in laminar annular pipe flows. ADV POWDER TECHNOL 2020. [DOI: 10.1016/j.apt.2020.06.005] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
|
12
|
Bass K, Longest W. Development of Dry Powder Inhaler Patient Interfaces for Improved Aerosol Delivery to Children. AAPS PharmSciTech 2020; 21:157. [PMID: 32451773 DOI: 10.1208/s12249-020-01667-3] [Citation(s) in RCA: 12] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/06/2020] [Accepted: 03/22/2020] [Indexed: 02/06/2023] Open
Abstract
The objective of this study was to explore different internal flow passages in the patient interface region of a new air-jet-based dry powder inhaler (DPI) in order to minimize device and extrathoracic aerosol depositional losses using computational fluid dynamics (CFD) simulations. The best-performing flow passages were used for oral and nose-to-lung (N2L) aerosol delivery in pediatric extrathoracic airway geometries consistent with a 5-year-old child. Aerosol delivery conditions were based on a previously developed and tested air-jet DPI device and included a base flow rate of 13.3 LPM (delivered from a small ventilation bag) and an inhaled air volume of 750 mL. Initial CFD models of the system clearly established that deposition on either the back of the throat or nasal cannula bifurcation was strongly correlated with the maximum velocity exiting the flow passage. Of all designs tested, the combination of a 3D rod array and rapid expansion of the flow passage side walls was found to dramatically reduce interface and device deposition and improve lung delivery of the aerosol. For oral aerosol administration, the optimal flow passage compared with a base case reduced device, mouthpiece, and mouth-throat deposition efficiencies by factors of 8-, 3-, and 2-fold, respectively. For N2L aerosol administration, the optimal flow pathway compared with a base case reduced device, nasal cannula, and nose-throat deposition by 16-, 6-, and 1.3-fold, respectively. In conclusion, a new patient interface design including a 3D rod array and rapid expansion dramatically improved transmission efficiency of a dry powder aerosol.
Collapse
|
13
|
Venkatesh S, Sakthivel M, Saranav H, Saravanan N, Rathnakumar M, Santhosh K. Performance investigation of the combined series and parallel arrangement cyclone separator using experimental and CFD approach. POWDER TECHNOL 2020. [DOI: 10.1016/j.powtec.2019.10.087] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/25/2022]
|
14
|
Bass K, Farkas D, Longest W. Optimizing Aerosolization Using Computational Fluid Dynamics in a Pediatric Air-Jet Dry Powder Inhaler. AAPS PharmSciTech 2019; 20:329. [PMID: 31676991 DOI: 10.1208/s12249-019-1535-4] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/24/2019] [Accepted: 09/12/2019] [Indexed: 11/30/2022] Open
Abstract
The objective of this study was to optimize the performance of a high-efficiency pediatric inhaler, referred to as the pediatric air-jet DPI, using computational fluid dynamics (CFD) simulations with supporting experimental analysis of aerosol formation. The pediatric air-jet DPI forms an internal flow pathway consisting of an inlet jet of high-speed air, capsule chamber containing a powder formulation, and outlet orifice. Instead of simulating full breakup of the powder bed to an aerosol in this complex flow system, which is computationally expensive, flow-field-based dispersion parameters were sought that correlated with experimentally determined aerosolization metrics. For the pediatric air-jet DPI configuration that was considered, mass median aerodynamic diameter (MMAD) directly correlated with input turbulent kinetic energy normalized by actuation pressure and flow kinetic energy. Emitted dose (ED) correlated best with input flow rate multiplied by the ratio of capillary diameters. Based on these dispersion parameters, an automated CFD process was used over multiple iterations of over 100 designs to identify optimal inlet and outlet capillary diameters, which affected system performance in complex and unexpected ways. Experimental verification of the optimized designs indicated an MMAD < 1.6 μm and an ED > 90% of loaded dose. While extrathoracic depositional loss will be determined in future studies, at an operating flow rate of 15 L/min, it is expected that pediatric mouth-throat or even nose-throat aerosol deposition fractions will be below 10% and potentially less than 5% representing a significant improvement in the delivery efficiency of dry powder pharmaceutical aerosols to children.
Collapse
|
15
|
Multi-objective optimization design of spiral demister with punched holes by combining response surface method and genetic algorithm. POWDER TECHNOL 2019. [DOI: 10.1016/j.powtec.2019.07.030] [Citation(s) in RCA: 12] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/22/2022]
|
16
|
Lavanya MN, Dutta S, Moses JA, Chinnaswamy A. Development of β‐carotene aerosol formulations using a modified spray dryer. J FOOD PROCESS ENG 2019. [DOI: 10.1111/jfpe.13233] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/05/2023]
Affiliation(s)
- Malur Narayanaswamy Lavanya
- Computational Modeling and Nanoscale Processing Unit, Indian Institute of Food Processing Technology (IIFPT)Ministry of Food Processing Industries, Govt. of India Tamil Nadu India
| | - Sayantani Dutta
- Computational Modeling and Nanoscale Processing Unit, Indian Institute of Food Processing Technology (IIFPT)Ministry of Food Processing Industries, Govt. of India Tamil Nadu India
| | - Jeyan Arthur Moses
- Computational Modeling and Nanoscale Processing Unit, Indian Institute of Food Processing Technology (IIFPT)Ministry of Food Processing Industries, Govt. of India Tamil Nadu India
| | - Anandharamakrishnan Chinnaswamy
- Computational Modeling and Nanoscale Processing Unit, Indian Institute of Food Processing Technology (IIFPT)Ministry of Food Processing Industries, Govt. of India Tamil Nadu India
| |
Collapse
|
17
|
Longest W, Farkas D, Bass K, Hindle M. Use of Computational Fluid Dynamics (CFD) Dispersion Parameters in the Development of a New DPI Actuated with Low Air Volumes. Pharm Res 2019; 36:110. [PMID: 31139939 DOI: 10.1007/s11095-019-2644-1] [Citation(s) in RCA: 20] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/29/2019] [Accepted: 05/10/2019] [Indexed: 11/26/2022]
Abstract
PURPOSE To determine the predictive power of computational fluid dynamics (CFD)-based dispersion parameters in the development of a new inline DPI that is actuated with low volumes of air. METHODS Four new versions of a dose aerosolization and containment (DAC)-unit DPI were created with varying inlet and outlet orifice sizes and analyzed with results from five previous designs. A concurrent in vitro and CFD analysis was conducted to predict the emitted dose (ED; as a % of loaded dose) and aerosol mass median aerodynamic diameter (MMAD) produced by each device when actuated with 10 ml air bursts. CFD simulations of device operation were used to predict flow field and particle-based dispersion parameters. RESULTS Comparisons of experimental and CFD results indicated that multiple flow field and particle-based dispersion parameters could be used to predict ED (minimum RMS Error = 4.9%) and MMAD (minimum RMS Error = 0.04 μm) to a high degree of accuracy. Based on experiments, the best overall device produced mean (standard deviation; SD) ED = 82.9(4.3)% and mean MMAD (SD) = 1.73(0.07)μm, which were in close agreement with the CFD predictions. CONCLUSIONS A unique relationship was identified in the DAC-unit DPI in which reducing turbulence also reduced the MMAD.
Collapse
Affiliation(s)
- Worth Longest
- Department of Mechanical and Nuclear Engineering, Virginia Commonwealth University, 401 West Main Street, P.O. Box 843015, Richmond, Virginia, 23284-3015, USA.
- Department of Pharmaceutics, Virginia Commonwealth University, Richmond, Virginia, USA.
| | - Dale Farkas
- Department of Mechanical and Nuclear Engineering, Virginia Commonwealth University, 401 West Main Street, P.O. Box 843015, Richmond, Virginia, 23284-3015, USA
| | - Karl Bass
- Department of Mechanical and Nuclear Engineering, Virginia Commonwealth University, 401 West Main Street, P.O. Box 843015, Richmond, Virginia, 23284-3015, USA
| | - Michael Hindle
- Department of Pharmaceutics, Virginia Commonwealth University, Richmond, Virginia, USA
| |
Collapse
|
18
|
Mucus transport and distribution by steady expiration in an idealized airway geometry. Med Eng Phys 2019; 66:26-39. [DOI: 10.1016/j.medengphy.2019.02.006] [Citation(s) in RCA: 16] [Impact Index Per Article: 3.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/15/2018] [Revised: 12/12/2018] [Accepted: 02/10/2019] [Indexed: 11/17/2022]
|
19
|
Longest W, Farkas D. Development of a New Inhaler for High-Efficiency Dispersion of Spray-Dried Powders Using Computational Fluid Dynamics (CFD) Modeling. AAPS JOURNAL 2019; 21:25. [PMID: 30734133 DOI: 10.1208/s12248-018-0281-y] [Citation(s) in RCA: 19] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 08/24/2018] [Accepted: 12/01/2018] [Indexed: 01/05/2023]
Abstract
Computational fluid dynamics (CFD) modeling offers a powerful tool for the development of drug delivery devices using a first principles approach but has been underutilized in the development of pharmaceutical inhalers. The objective of this study was to develop quantitative correlations for predicting the aerosolization behavior of a newly proposed dry powder inhaler (DPI). The dose aerosolization and containment (DAC) unit DPI utilizes inlet and outlet air orifices designed to maximize the dispersion of spray-dried powders, typically with low air volumes (~ 10 mL) and relatively low airflow rates (~ 3 L/min). Five DAC unit geometries with varying orifice outlet sizes, configurations, and protrusion distances were considered. Aerosolization experiments were performed using cascade impaction to determine mean device emitted dose (ED) and mass median aerodynamic diameter (MMAD). Concurrent CFD simulations were conducted to predict both flow field-based and particle-based dispersion parameters that captured different measures of turbulence. Strong quantitative correlations were established between multiple measures of turbulence and the experimentally observed aerosolization metrics of ED and MMAD. As expected, increasing turbulence produced increased ED with best case values reaching 85% of loaded dose. Surprisingly, decreasing turbulence produced an advantageous decrease in MMAD with values as low as approximately 1.6 μm, which is in contrast with previous studies. In conclusion, CFD provided valuable insights into the performance of the DAC unit DPI as a new device including a two-stage aerosolization process offering multiple avenues for future enhancements.
Collapse
Affiliation(s)
- Worth Longest
- Department of Mechanical and Nuclear Engineering, Virginia Commonwealth University, 401 West Main Street, P. O. Box 843015, Richmond, Virginia, 23284, USA. .,Department of Pharmaceutics, Virginia Commonwealth University, 410 North 12th Street, P.O. Box 980533, Richmond, Virginia, 23284, USA.
| | - Dale Farkas
- Department of Mechanical and Nuclear Engineering, Virginia Commonwealth University, 401 West Main Street, P. O. Box 843015, Richmond, Virginia, 23284, USA
| |
Collapse
|
20
|
Longest PW, Bass K, Dutta R, Rani V, Thomas ML, El-Achwah A, Hindle M. Use of computational fluid dynamics deposition modeling in respiratory drug delivery. Expert Opin Drug Deliv 2019; 16:7-26. [PMID: 30463458 PMCID: PMC6529297 DOI: 10.1080/17425247.2019.1551875] [Citation(s) in RCA: 43] [Impact Index Per Article: 8.6] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/30/2018] [Accepted: 11/20/2018] [Indexed: 10/27/2022]
Abstract
INTRODUCTION Respiratory drug delivery is a surprisingly complex process with a number of physical and biological challenges. Computational fluid dynamics (CFD) is a scientific simulation technique that is capable of providing spatially and temporally resolved predictions of many aspects related to respiratory drug delivery from initial aerosol formation through respiratory cellular drug absorption. AREAS COVERED This review article focuses on CFD-based deposition modeling applied to pharmaceutical aerosols. Areas covered include the development of new complete-airway CFD deposition models and the application of these models to develop a next-generation of respiratory drug delivery strategies. EXPERT OPINION Complete-airway deposition modeling is a valuable research tool that can improve our understanding of pharmaceutical aerosol delivery and is already supporting medical hypotheses, such as the expected under-treatment of the small airways in asthma. These complete-airway models are also being used to advance next-generation aerosol delivery strategies, like controlled condensational growth. We envision future applications of CFD deposition modeling to reduce the need for human subject testing in developing new devices and formulations, to help establish bioequivalence for the accelerated approval of generic inhalers, and to provide valuable new insights related to drug dissolution and clearance leading to microdosimetry maps of drug absorption.
Collapse
Affiliation(s)
- P. Worth Longest
- Department of Mechanical and Nuclear Engineering, Virginia Commonwealth University, Richmond, VA, USA
- Department of Pharmaceutics, Virginia Commonwealth University, Richmond, VA, USA
| | - Karl Bass
- Department of Mechanical and Nuclear Engineering, Virginia Commonwealth University, Richmond, VA, USA
| | - Rabijit Dutta
- Department of Mechanical and Nuclear Engineering, Virginia Commonwealth University, Richmond, VA, USA
| | - Vijaya Rani
- Department of Mechanical and Nuclear Engineering, Virginia Commonwealth University, Richmond, VA, USA
| | - Morgan L. Thomas
- Department of Mechanical and Nuclear Engineering, Virginia Commonwealth University, Richmond, VA, USA
| | - Ahmad El-Achwah
- Department of Mechanical and Nuclear Engineering, Virginia Commonwealth University, Richmond, VA, USA
| | - Michael Hindle
- Department of Pharmaceutics, Virginia Commonwealth University, Richmond, VA, USA
| |
Collapse
|
21
|
Bourke JM, Porter WR, Witmer LM. Convoluted nasal passages function as efficient heat exchangers in ankylosaurs (Dinosauria: Ornithischia: Thyreophora). PLoS One 2018; 13:e0207381. [PMID: 30566469 PMCID: PMC6300222 DOI: 10.1371/journal.pone.0207381] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/10/2018] [Accepted: 10/30/2018] [Indexed: 11/18/2022] Open
Abstract
Convoluted nasal passages are an enigmatic hallmark of Ankylosauria. Previous research suggested that these convoluted nasal passages functioned as heat exchangers analogous to the respiratory turbinates of mammals and birds. We tested this hypothesis by performing a computational fluid dynamic analysis on the nasal passages of two ankylosaurs: Panoplosaurus mirus and Euoplocephalus tutus. Our models predicted that Panoplosaurus and Euoplocephalus would have required 833 and 1568 thermal calories, respectively, to warm a single breath of air by 20°C. Heat recovery during exhalation resulted in energy savings of 65% for Panoplosaurus and 84% for Euoplocephalus. Our results fell well within the range of values for heat and water savings observed in extant terrestrial amniotes. We further tested alternate airway reconstructions that removed nasal passage convolutions or reduced nasal vestibule length. Our results revealed that the extensive elaboration observed in the nasal vestibules of ankylosaurs was a viable alternative to respiratory turbinates with regards to air conditioning. Of the two dinosaurs tested, Euoplocephalus repeatedly exhibited a more efficient nasal passage than Panoplosaurus. We suggest that the higher heat loads associated with the larger body mass of Euoplocephalus necessitated these more efficient nasal passages. Our findings further indicate that the evolution of complicated airways in dinosaurs may have been driven by the thermal requirements of maintaining cerebral thermal homeostasis.
Collapse
Affiliation(s)
- Jason M. Bourke
- Department of Biological Sciences, Ohio University, Athens, Ohio, United States of America
- * E-mail:
| | - Wm. Ruger Porter
- Department of Biomedical Sciences, Heritage College of Osteopathic Medicine, Ohio University, Athens, Ohio, United States of America
| | - Lawrence M. Witmer
- Department of Biomedical Sciences, Heritage College of Osteopathic Medicine, Ohio University, Athens, Ohio, United States of America
| |
Collapse
|
22
|
Das P, Nof E, Amirav I, Kassinos SC, Sznitman J. Targeting inhaled aerosol delivery to upper airways in children: Insight from computational fluid dynamics (CFD). PLoS One 2018; 13:e0207711. [PMID: 30458054 PMCID: PMC6245749 DOI: 10.1371/journal.pone.0207711] [Citation(s) in RCA: 31] [Impact Index Per Article: 5.2] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/14/2018] [Accepted: 11/03/2018] [Indexed: 11/28/2022] Open
Abstract
Despite the prevalence of inhalation therapy in the treatment of pediatric respiratory disorders, most prominently asthma, the fraction of inhaled drugs reaching the lungs for maximal efficacy remains adversely low. By and large drug delivery devices and their inhalation guidelines are typically derived from adult studies with child dosages adapted according to body weight. While it has long been recognized that physiological (e.g. airway sizes, breathing maneuvers) and physical transport (e.g. aerosol dynamics) characteristics are critical in governing deposition outcomes, such knowledge has yet to be extensively adapted to younger populations. Motivated by such shortcomings, the present work leverages in a first step in silico computational fluid dynamics (CFD) to explore opportunities for augmenting aerosol deposition in children based on respiratory physiological and physical transport determinants. Using an idealized, anatomically-faithful upper airway geometry, airflow and aerosol motion are simulated as a function of age, spanning a five year old to an adult. Breathing conditions mimic realistic age-specific inhalation maneuvers representative of Dry Powder Inhalers (DPI) and nebulizer inhalation. Our findings point to the existence of a single dimensionless curve governing deposition in the conductive airways via the dimensionless Stokes number (Stk). Most significantly, we uncover the existence of a distinct deposition peak irrespective of age. For the DPI simulations, this peak (∼ 80%) occurs at Stk ≈ 0.06 whereas for nebulizer simulations, the corresponding peak (∼ 45%) occurs in the range of Stk between 0.03-0.04. Such dimensionless findings hence translate to an optimal window of micron-sized aerosols that evolves with age and varies with inhalation device. The existence of such deposition optima advocates revisiting design guidelines for optimizing deposition outcomes in pediatric inhalation therapy.
Collapse
Affiliation(s)
- Prashant Das
- Department of Biomedical Engineering, Technion - Israel Institute of Technology, Haifa, Israel
| | - Eliram Nof
- Department of Biomedical Engineering, Technion - Israel Institute of Technology, Haifa, Israel
| | - Israel Amirav
- Department of Pediatrics, Faculty of Medicine & Dentistry, University of Alberta, Edmonton, Alberta, Canada
| | - Stavros C. Kassinos
- Computational Sciences Laboratory (UCY-CompSci), Department of Mechanical and Manufacturing Engineering, University of Cyprus, Kallipoleos Avenue 75, Nicosia 1678, Cyprus
| | - Josué Sznitman
- Department of Biomedical Engineering, Technion - Israel Institute of Technology, Haifa, Israel
| |
Collapse
|
23
|
Computational Fluid Dynamics Modeling of Respiratory Airflow in Tracheobronchial Airways of Infant, Child, and Adult. COMPUTATIONAL AND MATHEMATICAL METHODS IN MEDICINE 2018; 2018:9603451. [PMID: 30515236 PMCID: PMC6234445 DOI: 10.1155/2018/9603451] [Citation(s) in RCA: 14] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 09/11/2018] [Accepted: 10/16/2018] [Indexed: 12/19/2022]
Abstract
During human growth and development from infancy to adulthood, dramatic changes occur in the respiratory system. It is important to understand respiratory airflow in different age groups in age-specific treatment of respiratory disorders. This study numerically investigated the age-related effects on inspiratory and expiratory airflow dynamics in four-generation lung airway models under normal breathing conditions. Tracheobronchial airway models of infant (6 month old), child (5 years old), and adult (25 years old) from sixth to ninth generations were constructed for the study. Computational fluid dynamics (CFD) was used to solve the equations governing the airflow. Results of this study indicate that as age increases, airflow velocity, pressure, and wall shear stress decrease for both inspiration and expiration in this particular subregion of the respiratory tract. During inspiration, the splitting of velocity streamlines at bifurcations increases with age. The opposite situation merging happens during expiration, and it also increases with age. The level of splitting and merging of streamlines here reflects the influence of respiratory mechanics in the age groups. The computational models provide new information on characteristics and patterns of age-dependent respiratory airflow in the sixth to ninth generations of tracheobronchial airways and can be applied in other generations.
Collapse
|
24
|
An Analytical Review of the Numerical Methods used for Finite Element Modeling of Traumatic Brain Injury. Ann Biomed Eng 2018; 47:1855-1872. [DOI: 10.1007/s10439-018-02161-5] [Citation(s) in RCA: 51] [Impact Index Per Article: 8.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/13/2018] [Accepted: 10/22/2018] [Indexed: 01/24/2023]
|
25
|
Computational modelling and experimental validation of drug entrainment in a dry powder inhaler. Int J Pharm 2018; 553:37-46. [PMID: 30316002 DOI: 10.1016/j.ijpharm.2018.10.021] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/02/2018] [Revised: 10/02/2018] [Accepted: 10/07/2018] [Indexed: 11/22/2022]
Abstract
In a passive dry powder inhaler (DPI) a patient inhales to entrain drug powder. The goal of this study is to demonstrate experimentally that an Eulerian-Eulerian (EE CFD) computational fluid dynamics (CFD) method can accurately predict the entrainment of the dry powder formulation in DPIs. A CFD method that makes accurate predictions of the entrainment process can be applied in DPI design and optimization processes. Three different DPI entrainment geometries were tested. For each geometry, a transparent entrainment module was prepared. In each experiment, the chosen entrainment module was first filled with lactose powder and attached to an inhalation simulator (a computer controlled pump). The entrainment process was recorded with a high-speed camera. The resulting video footage was analysed and compared with CFD predictions. The observed distribution of powder in the entrainment compartment and the measured rate of drug entrainment were in good agreement with CFD predictions. Through a process of experimental validation, this study established the first demonstration that two-dimensional EE CFD methodology provides robust and accurate predictions of aerosol generation from DPI entrainment chambers. The findings support the wider application of EE CFD for the design optimization of DPI devices.
Collapse
|
26
|
Bass K, Longest PW. Recommendations for Simulating Microparticle Deposition at Conditions Similar to the Upper Airways with Two-Equation Turbulence Models. JOURNAL OF AEROSOL SCIENCE 2018; 119:31-50. [PMID: 30349146 PMCID: PMC6195318 DOI: 10.1016/j.jaerosci.2018.02.007] [Citation(s) in RCA: 22] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/22/2023]
Abstract
The development of a CFD model, from initial geometry to experimentally validated result with engineering insight, can be a time-consuming process that often requires several iterations of meshing and solver set-up. Applying a set of guidelines in the early stages can help to streamline the process and improve consistency between different models. The objective of this study was to determine both mesh and CFD solution parameters that enable the accurate simulation of microparticle deposition under flow conditions consistent with the upper respiratory airways including turbulent flow. A 90° bend geometry was used as a characteristic model that occurs throughout the airways and for which high-quality experimental aerosol deposition data is available in the transitional and turbulent flow regimes. Four meshes with varying degrees of near-wall resolution were compared, and key solver settings were applied to determine the parameters that minimize sensitivity to the near-wall (NW) mesh. The Low Reynolds number (LRN) k-ω model was used to resolve the turbulence field, which is a numerically efficient two-equation turbulence model, but has recently been considered overly simplistic. Some recent studies have used more complex turbulence models, such as Large Eddy Simulation (LES), to overcome the perceived weaknesses of two-equation models. Therefore, the secondary objective was to determine whether the more computationally efficient LRN k-ω model was capable of providing deposition results that were comparable to LES. Results show how NW mesh sensitivity is reduced through application of the Green-Gauss Node-based gradient discretization scheme and physically realistic near-wall corrections. Using the newly recommended meshing parameters and solution guidelines gives an excellent match to experimental data. Furthermore, deposition data from the LRN k-ω model compares favorably with LES results for the same characteristic geometry. In summary, this study provides a set of meshing and solution guidelines for simulating aerosol deposition in transitional and turbulent flows found in the upper respiratory airways using the numerically efficient LRN k-ω approach.
Collapse
Affiliation(s)
- Karl Bass
- Department of Mechanical Engineering, Virginia Commonwealth University, Richmond, VA
| | - P. Worth Longest
- Department of Mechanical Engineering, Virginia Commonwealth University, Richmond, VA
- Department of Pharmaceutics, Virginia Commonwealth University, Richmond, VA
- Dr. P. Worth Longest, PhD (Corresponding author), Virginia Commonwealth University, 401 West Main Street, P.O. Box 843015, Richmond, VA 23284-3015, Phone: (804)-827-7023, Fax: (804)-827-7030, , Mr. Karl Bass, Virginia Commonwealth University, 401 West Main Street, P.O. Box 843015, Richmond, VA 23284-3015, Phone: (804)-827-7023, Fax: (804)-827-7030,
| |
Collapse
|
27
|
Koullapis P, Hofemeier P, Sznitman J, Kassinos S. An efficient computational fluid-particle dynamics method to predict deposition in a simplified approximation of the deep lung. Eur J Pharm Sci 2018; 113:132-144. [DOI: 10.1016/j.ejps.2017.09.016] [Citation(s) in RCA: 25] [Impact Index Per Article: 4.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/15/2017] [Revised: 09/07/2017] [Accepted: 09/08/2017] [Indexed: 10/18/2022]
|
28
|
Deng Q, Ou C, Chen J, Xiang Y. Particle deposition in tracheobronchial airways of an infant, child and adult. THE SCIENCE OF THE TOTAL ENVIRONMENT 2018; 612:339-346. [PMID: 28854390 DOI: 10.1016/j.scitotenv.2017.08.240] [Citation(s) in RCA: 56] [Impact Index Per Article: 9.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 06/15/2017] [Revised: 08/22/2017] [Accepted: 08/22/2017] [Indexed: 05/04/2023]
Abstract
BACKGROUND Particle deposition in human airways is important for assessing both health effects of inhaled particles and therapeutic efficacy of inhaled drug aerosols, but is not well understood for infants and children. OBJECTIVE We investigate particle deposition in infants and children by using computational fluid dynamics (CFD), and compare this with particle deposition in adults. METHODS We chose three population age groups: 7-month infant, 4-year old child, and 20-year old adult. Both airway structures and breathing conditions are considered to vary as a human grows from infancy to adulthood. We investigated deposition of micron-size particles (1-10μm) in both the upper (G3-G6) and lower (G9-G12) tracheobronchial (TB) airways under sedentary conditions. RESULTS We found that particle deposition in both upper and lower airways is the highest in an infant, next in a child, and lowest in an adult. As age increases, particle deposition decreases in the upper airways but increases in the lower. For infants, inertial impaction is the dominant deposition mechanism, thus particles are deposited more in the upper airways than in the lower. However, particles are deposited more in the lower airways than in the upper in adults, as gravitational sedimentation is the dominant deposition mechanism. CONCLUSION Given the differences in the airway structure and particle deposition mechanisms, particle deposition in infants and children differs from that in adults, not only in the efficiency of deposition but also in the site. Our findings provide evidence that "children are not small adults".
Collapse
Affiliation(s)
- Qihong Deng
- School of Energy Science and Engineering, Central South University, Changsha, China; XiangYa School of Public Health, Central South University, Changsha, China.
| | - Cuiyun Ou
- School of Energy Science and Engineering, Central South University, Changsha, China
| | - Jiao Chen
- School of Energy Science and Engineering, Central South University, Changsha, China
| | - Yuguang Xiang
- School of Energy Science and Engineering, Central South University, Changsha, China
| |
Collapse
|
29
|
Numerical error analysis for three-dimensional CFD simulations in the two-room model containment THAI + : Grid convergence index, wall treatment error and scalability tests. NUCLEAR ENGINEERING AND DESIGN 2018. [DOI: 10.1016/j.nucengdes.2017.11.010] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/17/2022]
|
30
|
Numerical investigation of particle deposition in a triple bifurcation airway due to gravitational sedimentation and inertial impaction. POWDER TECHNOL 2018. [DOI: 10.1016/j.powtec.2017.09.050] [Citation(s) in RCA: 27] [Impact Index Per Article: 4.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/17/2022]
|
31
|
Ghaffari M, Tangen K, Alaraj A, Du X, Charbel FT, Linninger AA. Large-scale subject-specific cerebral arterial tree modeling using automated parametric mesh generation for blood flow simulation. Comput Biol Med 2017; 91:353-365. [PMID: 29126049 DOI: 10.1016/j.compbiomed.2017.10.028] [Citation(s) in RCA: 19] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/18/2017] [Revised: 10/23/2017] [Accepted: 10/23/2017] [Indexed: 10/18/2022]
Abstract
In this paper, we present a novel technique for automatic parametric mesh generation of subject-specific cerebral arterial trees. This technique generates high-quality and anatomically accurate computational meshes for fast blood flow simulations extending the scope of 3D vascular modeling to a large portion of cerebral arterial trees. For this purpose, a parametric meshing procedure was developed to automatically decompose the vascular skeleton, extract geometric features and generate hexahedral meshes using a body-fitted coordinate system that optimally follows the vascular network topology. To validate the anatomical accuracy of the reconstructed vasculature, we performed statistical analysis to quantify the alignment between parametric meshes and raw vascular images using receiver operating characteristic curve. Geometric accuracy evaluation showed an agreement with area under the curves value of 0.87 between the constructed mesh and raw MRA data sets. Parametric meshing yielded on-average, 36.6% and 21.7% orthogonal and equiangular skew quality improvement over the unstructured tetrahedral meshes. The parametric meshing and processing pipeline constitutes an automated technique to reconstruct and simulate blood flow throughout a large portion of the cerebral arterial tree down to the level of pial vessels. This study is the first step towards fast large-scale subject-specific hemodynamic analysis for clinical applications.
Collapse
Affiliation(s)
- Mahsa Ghaffari
- Department of Bioengineering, University of Illinois at Chicago, Chicago, IL, USA
| | - Kevin Tangen
- Department of Bioengineering, University of Illinois at Chicago, Chicago, IL, USA
| | - Ali Alaraj
- Department of Bioengineering, University of Illinois at Chicago, Chicago, IL, USA; Department of Neurosurgery, University of Illinois at Chicago, Chicago, IL, USA
| | - Xinjian Du
- Department of Neurosurgery, University of Illinois at Chicago, Chicago, IL, USA
| | - Fady T Charbel
- Department of Neurosurgery, University of Illinois at Chicago, Chicago, IL, USA
| | - Andreas A Linninger
- Department of Bioengineering, University of Illinois at Chicago, Chicago, IL, USA; Department of Neurosurgery, University of Illinois at Chicago, Chicago, IL, USA.
| |
Collapse
|
32
|
Guha A, Pradhan K. Secondary motion in three-dimensional branching networks. PHYSICS OF FLUIDS (WOODBURY, N.Y. : 1994) 2017; 29:063602. [PMID: 28713213 PMCID: PMC5489416 DOI: 10.1063/1.4984919] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 03/04/2017] [Accepted: 05/21/2017] [Indexed: 06/07/2023]
Abstract
A major aim of the present work is to understand and thoroughly document the generation, the three-dimensional distribution, and the evolution of the secondary motion as the fluid progresses downstream through a branched network. Six generations (G0-G5) of branches (involving 63 straight portions and 31 bifurcation modules) are computed in one go; such computational challenges are rarely taken in the literature. More than 30 × 106 computational elements are employed for high precision of computed results and fine quality of the flow visualization diagrams. The study of co-planar vis-à-vis non-planar space-filling configurations establishes a quantitative evaluation of the dependence of the fluid dynamics on the three-dimensional arrangement of the same individual branches. As compared to the secondary motion in a simple curved pipe, three distinctive features, viz., the change of shape and size of the flow-cross-section, the division of non-uniform primary flow in a bifurcation module, and repeated switchover from clockwise to anticlockwise curvature and vice versa in the flow path, make the present situation more complex. It is shown that the straight portions in the network, in general, attenuate the secondary motion, while the three-dimensionally complex bifurcation modules generate secondary motion and may alter the number, arrangement, and structure of vortices. A comprehensive picture of the evolution of quantitative flow visualizations of the secondary motion is achieved by constructing contours of secondary velocity [Formula: see text], streamwise vorticity [Formula: see text], and [Formula: see text] iso-surfaces. It is demonstrated, for example, that for in-plane configuration, the vortices on any plane appear in pair (i.e., for each clockwise rotating vortex, there is an otherwise identical anticlockwise vortex), whereas the vortices on a plane for the out-of-plane configuration may be dissimilar, and there may even be an odd number of vortices. We have formulated three new parameters (ES/P , [Formula: see text], and [Formula: see text]) for a quantitative description of the overall features of the secondary flow field. [Formula: see text] represents a non-uniformity index of the secondary flow in an individual branch, ES/P represents the mass-flow-averaged relative kinetic energy of the secondary motion in an individual branch, and [Formula: see text] provides a measure of the non-uniformity of the secondary flow between various branches of the same generation Gn. The repeated enhancement of the secondary kinetic energy in the bifurcation modules is responsible for the occurrence of significant values of ES/P even in generation G5. For both configurations, it is found that for any bifurcation module, the value of ES/P is greater in that daughter branch in which the mass-flow rate is greater. Even though the various contour plots of the complex secondary flow structure appear visually very different from one another, the values of [Formula: see text] are found to lie within a small range ([Formula: see text]) for the six-generation networks studied. It is shown that [Formula: see text] grows as the generation number Gn increases. It is established that the out-of-plane configuration, in general, creates more secondary kinetic energy (higher ES/P ), a similar level of non-uniformity in the secondary flow in an individual branch (similar [Formula: see text]), and a significantly lower level of non-uniformity in the distribution of secondary motion among various branches of the same generation (much lower [Formula: see text]), as compared to the in-plane arrangement of the same branches.
Collapse
Affiliation(s)
- Abhijit Guha
- Mechanical Engineering Department, Indian Institute of Technology Kharagpur, Kharagpur 721302, India
| | - Kaustav Pradhan
- Mechanical Engineering Department, Indian Institute of Technology Kharagpur, Kharagpur 721302, India
| |
Collapse
|
33
|
Longest PW, Tian G, Khajeh-Hosseini-Dalasm N, Hindle M. Validating Whole-Airway CFD Predictions of DPI Aerosol Deposition at Multiple Flow Rates. J Aerosol Med Pulm Drug Deliv 2016; 29:461-481. [PMID: 27082824 PMCID: PMC5220557 DOI: 10.1089/jamp.2015.1281] [Citation(s) in RCA: 43] [Impact Index Per Article: 5.4] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/27/2015] [Accepted: 02/05/2016] [Indexed: 11/13/2022] Open
Abstract
BACKGROUND The objective of this study was to compare aerosol deposition predictions of a new whole-airway CFD model with available in vivo data for a dry powder inhaler (DPI) considered across multiple inhalation waveforms, which affect both the particle size distribution (PSD) and particle deposition. METHODS The Novolizer DPI with a budesonide formulation was selected based on the availability of 2D gamma scintigraphy data in humans for three different well-defined inhalation waveforms. Initial in vitro cascade impaction experiments were conducted at multiple constant (square-wave) particle sizing flow rates to characterize PSDs. The whole-airway CFD modeling approach implemented the experimentally determined PSDs at the point of aerosol formation in the inhaler. Complete characteristic airway geometries for an adult were evaluated through the lobar bronchi, followed by stochastic individual pathway (SIP) approximations through the tracheobronchial region and new acinar moving wall models of the alveolar region. RESULTS It was determined that the PSD used for each inhalation waveform should be based on a constant particle sizing flow rate equal to the average of the inhalation waveform's peak inspiratory flow rate (PIFR) and mean flow rate [i.e., AVG(PIFR, Mean)]. Using this technique, agreement with the in vivo data was acceptable with <15% relative differences averaged across the three regions considered for all inhalation waveforms. Defining a peripheral to central deposition ratio (P/C) based on alveolar and tracheobronchial compartments, respectively, large flow-rate-dependent differences were observed, which were not evident in the original 2D in vivo data. CONCLUSIONS The agreement between the CFD predictions and in vivo data was dependent on accurate initial estimates of the PSD, emphasizing the need for a combination in vitro-in silico approach. Furthermore, use of the AVG(PIFR, Mean) value was identified as a potentially useful method for characterizing a DPI aerosol at a constant flow rate.
Collapse
Affiliation(s)
- P. Worth Longest
- Department of Mechanical and Nuclear Engineering, Virginia Commonwealth University, Richmond, Virginia
- Department of Pharmaceutics, Virginia Commonwealth University, Richmond, Virginia
| | - Geng Tian
- Department of Mechanical and Nuclear Engineering, Virginia Commonwealth University, Richmond, Virginia
| | | | - Michael Hindle
- Department of Pharmaceutics, Virginia Commonwealth University, Richmond, Virginia
| |
Collapse
|
34
|
Lotfi A, Simmons A, Barber T. Evaluation of Different Meshing Techniques for the Case of a Stented Artery. J Biomech Eng 2016; 138:4032502. [PMID: 26784359 DOI: 10.1115/1.4032502] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/18/2015] [Indexed: 02/01/2023]
Abstract
The formation and progression of in-stent restenosis (ISR) in bifurcated vessels may vary depending on the technique used for stenting. This study evaluates the effect of a variety of mesh styles on the accuracy and reliability of computational fluid dynamics (CFD) models in predicting these regions, using an idealized stented nonbifurcated model. The wall shear stress (WSS) and the near-stent recirculating vortices are used as determinants. The meshes comprise unstructured tetrahedral and polyhedral elements. The effects of local refinement, as well as higher-order elements such as prismatic inflation layers and internal hexahedral core, have also been examined. The uncertainty associated with individual mesh style was assessed through verification of calculations using the grid convergence index (GCI) method. The results obtained show that the only condition which allows the reliable comparison of uncertainty estimation between different meshing styles is that the monotonic convergence of grid solutions is in the asymptotic range. Comparisons show the superiority of a flow-adaptive polyhedral mesh over the commonly used adaptive and nonadaptive tetrahedral meshes in terms of resolving the near-stent flow features, GCI value, and prediction of WSS. More accurate estimation of hemodynamic factors was obtained using higher-order elements, such as hexahedral or prismatic grids. Incorporating these higher-order elements, however, was shown to introduce some degrees of numerical diffusion at the transitional area between the two meshes, not necessarily translating into high GCI value. Our data also confirmed the key role of local refinement in improving the performance and accuracy of nonadaptive mesh in predicting flow parameters in models of stented artery. The results of this study can provide a guideline for modeling biofluid domain in complex bifurcated arteries stented in regards to various stenting techniques.
Collapse
|
35
|
The effects of curvature and constriction on airflow and energy loss in pathological tracheas. Respir Physiol Neurobiol 2016; 234:69-78. [PMID: 27619197 DOI: 10.1016/j.resp.2016.09.002] [Citation(s) in RCA: 30] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/18/2016] [Revised: 08/01/2016] [Accepted: 09/01/2016] [Indexed: 11/23/2022]
Abstract
This paper considers factors that play a significant role in determining inspiratory pressure and energy losses in the human trachea. Previous characterisations of pathological geometry changes have focussed on relating airway constriction and subsequent pressure loss, however many pathologies that affect the trachea cause deviation, increased curvature, constriction or a combination of these. This study investigates the effects of these measures on tracheal flow mechanics, using the compressive goitre (a thyroid gland enlargement) as an example. Computational fluid dynamics simulations were performed in airways affected by goitres (with differing geometric consequences) and a normal geometry for comparison. Realistic airways, derived from medical images, were used because idealised geometries often oversimplify the complex anatomy of the larynx and its effects on the flow. Two mechanisms, distinct from stenosis, were found to strongly affect airflow energy dissipation in the pathological tracheas. The jet emanating from the glottis displayed different impingement and breakdown patterns in pathological geometries and increased loss was associated with curvature.
Collapse
|
36
|
Abstract
The effort required to inhale a breath of air is a critically important measure in assessing airway function. Although the contribution of the trachea to the total flow resistance of the airways is generally modest, pathological alterations in tracheal geometry can have a significant negative effect. This study investigates the mechanisms of flow energy loss in a healthy trachea and in four geometries affected by retrosternal goitre which can cause significant distortions of tracheal geometry including constriction and deviation with abnormal curvature. By separating out the component of energy loss related to the wall shear (frictional loss), striking differences are found between the patterns of energy dissipation in the normal and pathological tracheas. Furthermore the ratio of frictional to total loss is dramatically reduced in the pathological geometries.
Collapse
|
37
|
Bols J, Taelman L, De Santis G, Degroote J, Verhegghe B, Segers P, Vierendeels J. Unstructured hexahedral mesh generation of complex vascular trees using a multi-block grid-based approach. Comput Methods Biomech Biomed Engin 2015. [DOI: 10.1080/10255842.2015.1058925] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/23/2022]
|
38
|
|
39
|
|
40
|
Walenga RL, Tian G, Hindle M, Yelverton J, Dodson K, Longest PW. Variability in Nose-to-Lung Aerosol Delivery. JOURNAL OF AEROSOL SCIENCE 2014; 78:11-29. [PMID: 25308992 PMCID: PMC4187112 DOI: 10.1016/j.jaerosci.2014.08.003] [Citation(s) in RCA: 27] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/30/2023]
Abstract
Nasal delivery of lung targeted pharmaceutical aerosols is ideal for drugs that need to be administered during high flow nasal cannula (HFNC) gas delivery, but based on previous studies losses and variability through both the delivery system and nasal cavity are expected to be high. The objective of this study was to assess the variability in aerosol delivery through the nose to the lungs with a nasal cannula interface for conventional and excipient enhanced growth (EEG) delivery techniques. A database of nasal cavity computed tomography (CT) scans was collected and analyzed, from which four models were selected to represent a wide range of adult anatomies, quantified based on the nasal surface area-to-volume ratio (SA/V). Computational fluid dynamics (CFD) methods were validated with existing in vitro data and used to predict aerosol delivery through a streamlined nasal cannula and the four nasal models at a steady state flow rate of 30 L/min. Aerosols considered were solid particles for EEG delivery (initial 0.9 μm and 1.5 μm aerodynamic diameters) and conventional droplets (5 μm) for a control case. Use of the EEG approach was found to reduce depositional losses in the nasal cavity by an order of magnitude and substantially reduce variability. Specifically, for aerosol deposition efficiency in the four geometries, the 95% confidence intervals (CI) for 0.9 and 5 μm aerosols were 2.3-3.1 and 15.5-66.3%, respectively. Simulations showed that the use of EEG as opposed to conventional methods improved delivered dose of aerosols through the nasopharynx, expressed as penetration fraction (PF), by approximately a factor of four. Variability of PF, expressed by the coefficient of variation (CV), was reduced by a factor of four with EEG delivery compared with the control case. Penetration fraction correlated well with SA/V for larger aerosols, but smaller aerosols showed some dependence on nasopharyngeal exit hydraulic diameter. In conclusion, results indicated that the EEG technique not only improved lung aerosol delivery, but largely eliminated variability in both nasal depositional loss and lung PF in a newly developed set of nasal airway models.
Collapse
Affiliation(s)
- Ross L Walenga
- Department of Mechanical and Nuclear Engineering, Virginia Commonwealth University, Richmond, VA
| | - Geng Tian
- Department of Mechanical and Nuclear Engineering, Virginia Commonwealth University, Richmond, VA
| | - Michael Hindle
- Department of Pharmaceutics, Virginia Commonwealth University, Richmond, VA
| | - Joshua Yelverton
- Department of Otolaryngology – Head and Neck Surgery, Virginia Commonwealth University, Richmond, VA
| | - Kelley Dodson
- Department of Otolaryngology – Head and Neck Surgery, Virginia Commonwealth University, Richmond, VA
| | - P. Worth Longest
- Department of Mechanical and Nuclear Engineering, Virginia Commonwealth University, Richmond, VA
- Department of Pharmaceutics, Virginia Commonwealth University, Richmond, VA
| |
Collapse
|
41
|
Longest PW, Azimi M, Hindle M. Optimal delivery of aerosols to infants during mechanical ventilation. J Aerosol Med Pulm Drug Deliv 2014; 27:371-85. [PMID: 24299500 PMCID: PMC4227441 DOI: 10.1089/jamp.2013.1077] [Citation(s) in RCA: 15] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/09/2013] [Accepted: 10/17/2013] [Indexed: 02/06/2023] Open
Abstract
PURPOSE The objective of this study was to determine optimal aerosol delivery conditions for a full-term (3.6 kg) infant receiving invasive mechanical ventilation by evaluating the effects of aerosol particle size, a new wye connector, and timing of aerosol delivery. METHODS In vitro experiments used a vibrating mesh nebulizer and evaluated drug deposition fraction and emitted dose through ventilation circuits containing either a commercial (CM) or new streamlined (SL) wye connector and 3-mm endotracheal tube (ETT) for aerosols with mass median aerodynamic diameters of 880 nm, 1.78 μm, and 4.9 μm. The aerosol was released into the circuit either over the full inhalation cycle (T1 delivery) or over the first half of inhalation (T2 delivery). Validated computational fluid dynamics (CFD) simulations and whole-lung model predictions were used to assess lung deposition and exhaled dose during cyclic ventilation. RESULTS In vitro experiments at a steady-state tracheal flow rate of 5 L/min resulted in 80-90% transmission of the 880-nm and 1.78-μm aerosols from the ETT. Based on CFD simulations with cyclic ventilation, the SL wye design reduced depositional losses in the wye by a factor of approximately 2-4 and improved lung delivery efficiencies by a factor of approximately 2 compared with the CM device. Delivery of the aerosol over the first half of the inspiratory cycle (T2) reduced exhaled dose from the ventilation circuit by a factor of 4 compared with T1 delivery. Optimal lung deposition was achieved with the SL wye connector and T2 delivery, resulting in 45% and 60% lung deposition for optimal polydisperse (∼1.78 μm) and monodisperse (∼2.5 μm) particle sizes, respectively. CONCLUSIONS Optimization of selected factors and use of a new SL wye connector can substantially increase the lung delivery efficiency of medical aerosols to infants from current values of <1-10% to a range of 45-60%.
Collapse
Affiliation(s)
- P. Worth Longest
- Department of Mechanical and Nuclear Engineering, Virginia Commonwealth University, Richmond, VA
- Department of Pharmaceutics, Virginia Commonwealth University, Richmond, VA
| | - Mandana Azimi
- Department of Pharmaceutics, Virginia Commonwealth University, Richmond, VA
| | - Michael Hindle
- Department of Pharmaceutics, Virginia Commonwealth University, Richmond, VA
| |
Collapse
|
42
|
Walenga RL, Tian G, Longest PW. Development of characteristic upper tracheobronchial airway models for testing pharmaceutical aerosol delivery. J Biomech Eng 2014; 135:91010. [PMID: 23722698 DOI: 10.1115/1.4024630] [Citation(s) in RCA: 26] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Abstract] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/08/2013] [Accepted: 05/23/2013] [Indexed: 12/25/2022]
Abstract
Characteristic models of the upper conducting airways are needed to evaluate the performance of existing pharmaceutical inhalers and to develop new respiratory drug delivery strategies. Previous studies have focused on the development of characteristic mouth-throat (MT) geometries for orally inhaled products; however, characteristic upper tracheobronchial (TB) geometries are currently not available. In this study, a new characteristic model of the upper TB airways for an average adult male was developed based on an analysis of new and existing anatomical data. Validated computational fluid dynamics (CFD) simulations were used to evaluate the deposition of monodisperse and realistic polydisperse aerosols from multiple inhalers. Comparisons of deposition results between the new model and a simpler geometry were used to identify the effects of different anatomical features on aerosol deposition. The CFD simulations demonstrated a good match to regional pharmaceutical aerosol deposition from in vitro experiments in the same geometry. The deposition of both monodisperse and pharmaceutical aerosols was increased in the new TB geometry as a result of additional anatomical detail on a regional and highly localized basis. Tracheal features including an accurate coronal angle, asymmetry, and curvature produced a skewed laryngeal jet and significantly increased regional deposition. Branch curvature and realistic cross-sections increased deposition in the remainder of the TB model. A hexahedral mesh style was utilized to provide the best solution. In conclusion, a number of physiological features in the upper TB region were shown to influence deposition and should be included in a characteristic model of respiratory drug delivery.
Collapse
Affiliation(s)
- Ross L Walenga
- Department of Mechanical and Nuclear Engineering, Virginia Commonwealth University, Richmond, VA 23284, USA.
| | | | | |
Collapse
|
43
|
Oakes JM, Marsden AL, Grandmont C, Shadden SC, Darquenne C, Vignon-Clementel IE. Airflow and particle deposition simulations in health and emphysema: from in vivo to in silico animal experiments. Ann Biomed Eng 2014; 42:899-914. [PMID: 24318192 PMCID: PMC4092242 DOI: 10.1007/s10439-013-0954-8] [Citation(s) in RCA: 46] [Impact Index Per Article: 4.6] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/15/2013] [Accepted: 11/23/2013] [Indexed: 10/25/2022]
Abstract
Image-based in silico modeling tools provide detailed velocity and particle deposition data. However, care must be taken when prescribing boundary conditions to model lung physiology in health or disease, such as in emphysema. In this study, the respiratory resistance and compliance were obtained by solving an inverse problem; a 0D global model based on healthy and emphysematous rat experimental data. Multi-scale CFD simulations were performed by solving the 3D Navier-Stokes equations in an MRI-derived rat geometry coupled to a 0D model. Particles with 0.95 μm diameter were tracked and their distribution in the lung was assessed. Seven 3D-0D simulations were performed: healthy, homogeneous, and five heterogeneous emphysema cases. Compliance (C) was significantly higher (p = 0.04) in the emphysematous rats (C = 0.37 ± 0.14 cm(3)/cmH2O) compared to the healthy rats (C = 0.25 ± 0.04 cm(3)/cmH2O), while the resistance remained unchanged (p = 0.83). There were increases in airflow, particle deposition in the 3D model, and particle delivery to the diseased regions for the heterogeneous cases compared to the homogeneous cases. The results highlight the importance of multi-scale numerical simulations to study airflow and particle distribution in healthy and diseased lungs. The effect of particle size and gravity were studied. Once available, these in silico predictions may be compared to experimental deposition data.
Collapse
Affiliation(s)
- Jessica M Oakes
- Mechanical and Aerospace Engineering Department, University of California, San Diego, La Jolla, CA, 92093, USA
| | | | | | | | | | | |
Collapse
|
44
|
Tarjuelo-Gutierrez J, Rodriguez-Vila B, Pierce DM, Fastl TE, Verbrugghe P, Fourneau I, Maleux G, Herijgers P, Holzapfel GA, Gomez EJ. High-quality conforming hexahedral meshes of patient-specific abdominal aortic aneurysms including their intraluminal thrombi. Med Biol Eng Comput 2013; 52:159-68. [PMID: 24306943 DOI: 10.1007/s11517-013-1127-5] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/18/2013] [Accepted: 11/08/2013] [Indexed: 10/25/2022]
Abstract
In order to perform finite element (FE) analyses of patient-specific abdominal aortic aneurysms, geometries derived from medical images must be meshed with suitable elements. We propose a semi-automatic method for generating conforming hexahedral meshes directly from contours segmented from medical images. Magnetic resonance images are generated using a protocol developed to give the abdominal aorta high contrast against the surrounding soft tissue. These data allow us to distinguish between the different structures of interest. We build novel quadrilateral meshes for each surface of the sectioned geometry and generate conforming hexahedral meshes by combining the quadrilateral meshes. The three-layered morphology of both the arterial wall and thrombus is incorporated using parameters determined from experiments. We demonstrate the quality of our patient-specific meshes using the element Scaled Jacobian. The method efficiently generates high-quality elements suitable for FE analysis, even in the bifurcation region of the aorta into the iliac arteries. For example, hexahedral meshes of up to 125,000 elements are generated in less than 130 s, with 94.8 % of elements well suited for FE analysis. We provide novel input for simulations by independently meshing both the arterial wall and intraluminal thrombus of the aneurysm, and their respective layered morphologies.
Collapse
Affiliation(s)
- J Tarjuelo-Gutierrez
- Bioengineering and Telemedicine Centre, ETSI de Telecomunicación, Universidad Politécnica de Madrid, Madrid, Spain,
| | | | | | | | | | | | | | | | | | | |
Collapse
|
45
|
Tian G, Longest PW, Li X, Hindle M. Targeting aerosol deposition to and within the lung airways using excipient enhanced growth. J Aerosol Med Pulm Drug Deliv 2013; 26:248-65. [PMID: 23286828 PMCID: PMC3826577 DOI: 10.1089/jamp.2012.0997] [Citation(s) in RCA: 56] [Impact Index Per Article: 5.1] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/19/2012] [Accepted: 10/10/2012] [Indexed: 10/27/2022] Open
Abstract
BACKGROUND Previous studies have characterized the size increase of combination submicrometer particles composed of a drug and hygroscopic excipient when exposed to typical airway thermodynamic conditions. The objective of this study was to determine the deposition and size increase characteristics of excipient enhanced growth (EEG) aerosols throughout the tracheobronchial (TB) airways and to evaluate the potential for targeted delivery. METHODS Submicrometer particles composed of a poorly water-soluble drug (insulin) and hygroscopic excipient (sodium chloride) were considered at drug:excipient mass ratios of 50:50 and 25:75. A previously validated computational fluid dynamics model was used to predict aerosol size increase and deposition in characteristic geometries of the mouth-throat (MT), upper TB airways through the third bifurcation (B3), and remaining TB airways through B15. Additional validation experiments were also performed for albuterol sulfate:mannitol particles. Both growth of combination particles and deposition are reported throughout the conducting airways for characteristic slow and deep (SD) and quick and deep (QD) inhalations. RESULTS For all EEG cases considered, MT deposition was less than 1% of the drug dose, which is at least one order of magnitude lower than with state-of-the-art and conventional inhalers. Final aerosol sizes exiting the TB region and entering the alveolar airways were all greater than 3 μm. For SD inhalation, deposition fractions of 20% were achieved in the lower TB region of B8-B15, which is a factor of 20-30×higher than conventional delivery devices. With QD inhalation, maximum alveolar delivery of 90% was observed. CONCLUSIONS Increasing the dose delivered to the lower TB region by a factor of 20-30×or achieving 90% delivery to the alveolar airways was considered effective aerosol targeting compared with conventional devices. The trend of higher flow rates resulting in better alveolar delivery of aerosols is unique to EEG and may be used to design highly efficient dry powder inhalers.
Collapse
Affiliation(s)
- Geng Tian
- Department of Mechanical and Nuclear Engineering, Virginia Commonwealth University, Richmond, VA 23284
| | - P. Worth Longest
- Department of Mechanical and Nuclear Engineering, Virginia Commonwealth University, Richmond, VA 23284
- Department of Pharmaceutics, Virginia Commonwealth University, Richmond, VA 23284
| | - Xiang Li
- Department of Pharmaceutics, Virginia Commonwealth University, Richmond, VA 23284
| | - Michael Hindle
- Department of Pharmaceutics, Virginia Commonwealth University, Richmond, VA 23284
| |
Collapse
|
46
|
Hybrid mesh for nasal airflow studies. COMPUTATIONAL AND MATHEMATICAL METHODS IN MEDICINE 2013; 2013:727362. [PMID: 23983811 PMCID: PMC3747603 DOI: 10.1155/2013/727362] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 05/14/2013] [Revised: 06/24/2013] [Accepted: 06/25/2013] [Indexed: 11/24/2022]
Abstract
The accuracy of the numerical result is closely related to mesh density as well as its distribution. Mesh plays a very significant role in the outcome of numerical simulation. Many nasal airflow studies have employed unstructured mesh and more recently hybrid mesh scheme has been utilized considering the complexity of anatomical architecture. The objective of this study is to compare the results of hybrid mesh with unstructured mesh and study its effect on the flow parameters inside the nasal cavity. A three-dimensional nasal cavity model is reconstructed based on computed tomographic images of a healthy Malaysian adult nose. Navier-Stokes equation for steady airflow is solved numerically to examine inspiratory nasal flow. The pressure drop obtained using the unstructured computational grid is about 22.6 Pa for a flow rate of 20 L/min, whereas the hybrid mesh resulted in 17.8 Pa for the same flow rate. The maximum velocity obtained at the nasal valve using unstructured grid is 4.18 m/s and that with hybrid mesh is around 4.76 m/s. Hybrid mesh reported lower grid convergence index (GCI) than the unstructured mesh. Significant differences between unstructured mesh and hybrid mesh are determined highlighting the usefulness of hybrid mesh for nasal airflow studies.
Collapse
|
47
|
Jacob RE, Colby SM, Kabilan S, Einstein DR, Carson JP. In situ casting and imaging of the rat airway tree for accurate 3D reconstruction. Exp Lung Res 2013; 39:249-57. [PMID: 23786464 DOI: 10.3109/01902148.2013.801535] [Citation(s) in RCA: 9] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/13/2022]
Abstract
The use of anatomically accurate, animal-specific airway geometries is important for understanding and modeling the physiology of the respiratory system. One approach for acquiring detailed airway architecture is to create a bronchial cast of the conducting airways. However, typical casting procedures either do not faithfully preserve the in vivo branching angles or produce rigid casts that when removed for imaging are fragile and thus easily damaged. We address these problems by creating an in situ bronchial cast of the conducting airways in rats that can be subsequently imaged in situ using three-dimensional micro-CT imaging. We also demonstrate that deformations in airway branch angles resulting from the casting procedure are small, and that these angle deformations can be reversed through an interactive adjustment of the segmented cast geometry. Animal work was approved by the Institutional Animal Care and Use Committee of Pacific Northwest National Laboratory.
Collapse
Affiliation(s)
- Richard E Jacob
- Department of Systems Toxicology, Pacific Northwest National Laboratory, Richland, Washington 99352, USA.
| | | | | | | | | |
Collapse
|
48
|
|
49
|
Longest PW, Son YJ, Holbrook L, Hindle M. Aerodynamic factors responsible for the deaggregation of carrier-free drug powders to form micrometer and submicrometer aerosols. Pharm Res 2013; 30:1608-27. [PMID: 23471640 DOI: 10.1007/s11095-013-1001-z] [Citation(s) in RCA: 47] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/12/2012] [Accepted: 02/01/2013] [Indexed: 12/18/2022]
Abstract
PURPOSE To employ in vitro experiments combined with computational fluid dynamics (CFD) analysis to determine which aerodynamic factors were most responsible for deaggregating carrier-free powders to form micrometer and submicrometer aerosols from a capsule-based platform. METHODS Eight airflow passages were evaluated for deaggregation of the aerosol including a standard constricted tube, impaction surface, 2D mesh, inward radial jets, and newly proposed 3D grids and rod arrays. CFD simulations were implemented to evaluate existing and new aerodynamic factors for deaggregation and in vitro experiments were used to evaluate performance of each inhaler. RESULTS For the carrier-free formulation considered, turbulence was determined to be the primary deaggregation mechanism. A strong quantitative correlation was established between the mass median diameter (MMD) and newly proposed non-dimensional specific dissipation (NDSD) factor, which accounts for turbulent energy, inverse of the turbulent length scale, and exposure time. A 3D rod array design with unidirectional elements maximized NDSD and produced the best deaggregation with MMD<1 μm. CONCLUSIONS The new NDSD parameter can be used to develop highly effective dry powder inhalers like the 3D rod array that can efficiently produce submicrometer aerosols for next-generation respiratory drug delivery applications.
Collapse
Affiliation(s)
- P Worth Longest
- Department of Mechanical and Nuclear Engineering, Virginia Commonwealth University, Richmond, Virginia, USA.
| | | | | | | |
Collapse
|
50
|
Improving pharmaceutical aerosol delivery during noninvasive ventilation: effects of streamlined components. Ann Biomed Eng 2013; 41:1217-32. [PMID: 23423706 DOI: 10.1007/s10439-013-0759-9] [Citation(s) in RCA: 34] [Impact Index Per Article: 3.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/14/2012] [Accepted: 02/04/2013] [Indexed: 10/27/2022]
Abstract
Aerosol delivery efficiency during noninvasive ventilation (NIV) is known to be low (~10%) and is associated with poor outcomes of aerosol therapy. The objective of this study was to demonstrate the benefit of redesigning ventilation circuit components using a streamlining approach to improve aerosol delivery during nasal high flow therapy in adults with a conventional-sized aerosol from a mesh nebulizer. The ventilation circuit consisted of a humidifier, mesh nebulizer, mixing T-connector (with 90° angle), 10 mm tubing, and nasal cannula interface. In vitro experiments and computational fluid dynamics analyses were used to evaluate depositional losses in a system of existing components and a newly proposed streamlined T-connector and cannula at flow rates of 30 and 45 LPM. Streamlined designs reduced deposition in the T-connector by a factor of 4. In the nasal cannula, the streamlined designs reduced depositional losses by factors of 1.25-2.0. With the streamlined designs, the highest emitted dose achieved was >40% for a conventional-sized aerosol at 30 LPM. Streamlined geometries offer an effective method to significantly improve the delivery of aerosols through components of NIV systems. This increase in delivery efficiency is important for new inhaled medications with narrow therapeutic windows, increased costs, or long delivery times.
Collapse
|