Numerical study of the impact of glottis properties on the airflow field in the human trachea using V-LES

Effects of vocal fold lesions on particle deposition in a mouth-throat model: A numerical study

OBJECTIVE

Dysphonia is very common worldwide and aerosol drug inhalation is an important treatment for patients with dysphonia. This study aimed to explore the effects of vocal fold (VF) lesions on the particle deposition pattern using computational modeling.

METHODS

A realistic mouth-throat (MT) model of a healthy adult was constructed based on computed tomography images. Small and large vocal fold lesions were incorporated in the original model. A steady inhalation flowrate of 15 and 30 liter per minute (LPM) was used as the velocity inlet and monodisperse particles with diameters of 5 to 10 µm were simulated.

RESULTS

Particles of larger size are more likely to be deposited in MT models, most of them distributed in oral cavity, oropharynx and supraglottis. The ideal sizes at 30 LPM ranged over 7-10 µm for healthy VFs and 6-8 µm for VF lesions. The best sizes at 15 LPM ranged over 6-8 µm for healthy VFs and 8-9 µm for VF lesions.

CONCLUSION

Based on this study, VF lesions influence the deposition pattern in the glottis obviously. The ideal sizes differ at the flow rates of 15 and 30 LPM.

Assessing airflow unsteadiness in the human respiratory tract under different expiration conditions

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.

Effects of vocal fold adduction on the particle deposition in the glottis: A numerical analysis and in vitro assessment

BACKGROUND

The efficacy of inhalation therapy depends on the drug deposition in the human respiratory tract. This study investigates the effects of vocal fold adduction on the particle deposition in the glottis.

METHODS

A realistic mouth-throat (MT) geometry was built based on CT images of a healthy adult (MT-A). Mild (MT-B) and great (MT-C) vocal fold (VF) adduction were incorporated in the original model. Monodisperse particles range in size from 3 to 12 μm were simulated at inspiration flow rates of 15, 30 and 45 L per minute (LPM). The regional deposition of drug aerosols was performed in 3D-printed models and quantified using high-performance liquid chromatography.

RESULTS

Both the numerical analysis and in vitro experiments show that most particles are deposited in the mouth, pharynx and supraglottis, while few are deposited in the glottis and subglottis. For most cases in MT-A, the particle quantity in glottis is lower than 0.02 N/mm2 at 15 and 30 LPM while they increase dramatically at 45 LPM. It peaked at 0.347 N/mm2 for 5-μm particles at 45 LPM in MT-B and 2.324 N/mm2 for 6-μm particles at 30 LPM in MT-C. The lowest drug mass faction in the glottis in vitro were found at 15 LPM for MT-A and MT-C, and at 30 LPM for MT-B, whereas it peaked at 45 LPM for all MT models, 0.71% in MT-A, 1.16% in MT-B, and 2.53% in MT-C, respectively.

CONCLUSION

Based on the results of this study, larger particles are more likely to be deposited in the oral cavity, oropharynx, and supraglottis than in the glottis. However, particle deposition in the glottis generally increases with VF adduction and greater inspiratory flow rates.

Computational analysis of human upper airway aerodynamics

There is a considerable interest in understanding transient human upper airway aerodynamics, especially in view of assessing the effects of various ventilation therapies. Experimental analyses in a patient-specific manner pose challenges as the upper airway consists of a narrow confined region with complex anatomy. Pressure measurements are feasible, but, for example, PIV experiments require special measures to accommodate for the light refraction by the model. Computational fluid dynamics can bridge the gap between limited experimental data and detailed flow features. This work aims to validate the use of combined lattice Boltzmann method and a large eddy scale model for simulating respiration, and to identify clinical features of the flow and show the clinical potential of the method. Airflow was computationally analyzed during a realistic, transient, breathing profile in an upper airway geometry ranging from nose to trachea, and the resulting pressure calculations were compared against in vitro experiments. Simulations were conducted on meshes containing about 1 billion cells to ensure accuracy and to capture intrinsic flow features. Airway pressures obtained from simulations and in vitro experiments are in good agreement both during inhalation and exhalation. High velocity pharyngeal and laryngeal jets and recirculation in the region of the olfactory cleft are observed. Graphical Abstract The Lattice-Boltzmann Method combined with Large Eddy Simulations was used to compute the aerodynamics in a human upper airway geometry. The left side of this graphical abstract shows the velocity and vorticity (middle figure in bottom row, and right figure of the right bottom figure) profiles at peak exhalation. The simulations were validated against experiments on a 3D-print of the geometry (shown in the top figures on the right hand side). The pressure drop (right bottom corner) shows a good agreement between experiments and simulations.

Effects of face shield on an emitter during a cough process: A large-eddy simulation study

Face shield is a common personal protection equipment for pandemic. In the present work, three-dimensional computational fluid dynamic (CFD) method is used to simulate a cough jet from an emitter who wears a face shield. A realistic manikin model with a simplified mouth cavity is employed. A large eddy simulation with a dynamic structure subgrid scale model is applied to model the turbulence. An Eulerian-Lagrangian approach is adopted to model the two-phase flows, with which the droplets are represented by a cloud of particles. The droplet breakup, evaporation, dispersion, drag force, and wall impingement are considered in this model. An inlet velocity profile that is based on a variable mouth opening area is considered. Special attentions have been put the vortex structure and droplet re-distribution induced by the face shield. It is found that the multiple vortices are formed when the cough jet impinges on the face shield. Some droplets move backward and others move downward after the impinging. It is also found that a small modification of the face shield significantly modifies the flow field and droplet distribution. We conclude that face shield significantly reduces the risk factor in the front of the emitter, meanwhile the risk factor in the back of the emitter increases. When the receiver standing in front of the emitter is shorter than the emitter, the risk is still very high. More attentions should be paid on the design of the face field, clothes cleaning and floor cleaning of the emitters with face shields. Based on the predicted droplet trajectory, a conceptual model for droplet flux is proposed for the scenario with the face shield.