Numerical and Experimental Analysis of Drug Inhalation in Realistic Human Upper Airway Model.

The demand for a more efficient and targeted method for intranasal drug delivery has led to sophisticated device design, delivery methods, and aerosol properties. Due to the complex nasal geometry and measurement limitations, numerical modeling is an appropriate approach to simulate the airflow, aerosol dispersion, and deposition for the initial assessment of novel methodologies for better drug delivery. In this study, a CT-based, 3D-printed model of a realistic nasal airway was reconstructed, and airflow pressure, velocity, turbulent kinetic energy (TKE), and aerosol deposition patterns were simultaneously investigated. Different inhalation flowrates (5, 10, 15, 30, and 45 L/min) and aerosol sizes (1, 1.5, 2.5, 3, 6, 15, and 30 µm) were simulated using laminar and SST viscous models, with the results compared and verified by experimental data. The results revealed that from the vestibule to the nasopharynx, the pressure drop was negligible for flow rates of 5, 10, and 15 L/min, while for flow rates of 30 and 40 L/min, a considerable pressure drop was observed by approximately 14 and 10%, respectively. However, from the nasopharynx and trachea, this reduction was approximately 70%. The aerosol deposition fraction alongside the nasal cavities and upper airway showed a significant difference in pattern, dependent on particle size. More than 90% of the initiated particles were deposited in the anterior region, while just under 20% of the injected ultrafine particles were deposited in this area. The turbulent and laminar models showed slightly different values for the deposition fraction and efficiency of drug delivery for ultrafine particles (about 5%); however, the deposition pattern for ultrafine particles was very different.

Investigation of airflow at different activity conditions in a realistic model of human upper respiratory tract.

In the present study, the turbulent flows inside a realistic model of the upper respiratory tract were investigated numerically and experimentally. The airway model included the geometrical details of the oral cavity to the end of the trachea that was based on a series of CT-scan images. The topological data of the respiratory tract were used for generating the computational model as well as the 3D-printed model that was used in the experimental pressure drop measurement. Different airflow rates of 30, 45, and 60 L/min, which correspond to the light, semi-light, and heavy activity breathing conditions, were investigated numerically using turbulence and transition models, as well as experimentally. Simulation results for airflow properties, including velocity vectors, pressure drops, streamlines, eddy viscosity, and turbulent kinetic energy contours in the oral-trachea airway model, were presented. The simulated pressure drop was compared with the experimental data, and reasonable agreement was found. The obtained results showed that the maximum pressure drop occurs in the narrowest part of the larynx region. A comparison between the numerical results and experimental data showed that the transition (γ-Reθ) SST model predicts higher pressure losses, especially at higher breathing rates. Formations of the secondary flows in the oropharynx and trachea regions were also observed. In addition, the simulation results showed that in the trachea region, the secondary flow structures dissipated faster for the flow rate of 60 L/min compared to the lower breathing rates of 30 and 45 L/min.

Laminar-to-turbulence and relaminarization zones detection by simulation of low Reynolds number turbulent blood flow in large stenosed arteries.

Laminar, turbulent, transitional, or combine areas of all three types of viscous flow can occur downstream of a stenosis depending upon the Reynolds number and constriction shape parameter. Neither laminar flow solver nor turbulent models for instance the k-ω (k-omega), k-ε (k-epsilon), RANS or LES are opportune for this type of flow. In the present study attention has been focused vigorously on the effect of the constriction in the flow field with a unique way. It means that the laminar solver was employed from entry up to the beginning of the turbulent shear flow. The turbulent model (k-ω SST Transitional Flows) was utilized from starting of turbulence to relaminarization zone while the laminar model was applied again with onset of the relaminarization district. Stenotic flows, with 50 and 75% cross-sectional area, were simulated at Reynolds numbers range from 500 to 2000 employing FLUENT (v6.3.17). The flow was considered to be steady, axisymmetric, and incompressible. Achieving results were reported as axial velocity, disturbance velocity, wall shear stress and the outcomes were compared with previously experimental and CFD computations. The analogy of axial velocity profiles shows that they are in acceptable compliance with the empirical data. As well as disturbance velocity and wall shear stresses anticipated by this new approach, part by part simulation, are reasonably valid with the acceptable experimental studies.