Influence of glottic aperture on the tracheal flow

The extra-thoracic mouth-throat area has a major influence on the aerosol delivery to the proximal or peripheral intra-thoracic airways. To characterize the particle deposition in this area, it is important to investigate first the flow structures in this crucial--in relation to the aerosol deposition--region. The glottis, which is delimited by the vocal cords and therefore has the narrowest passage, generates a laryngeal jet and a reverse flow downstream the glottis. It is generally assumed that the glottis has different shapes and cross-sectional areas at different moments during the respiratory cycle and also depends on the average inspiratory flow rate. Therefore, the influence of a circular glottal aperture, with a cross-sectional area of 90 mm2 and an elliptical and triangular shape, both with an area of 45 mm2, on the flow is investigated. However, the area of the circular aperture is twice as big as the area of the elliptical one, it has almost no influence on the flow structures. On the other hand, the triangular glottal aperture shifts the laryngeal jet in the direction of the posterior wall, and generates two pairs of counter rotating secondary vortices downstream the glottis, where the circular and elliptical only aperture generates one pair of vortices. The difference in pressure drop is more dominated by the cross-sectional area than by the shape of the glottis. This suggests the need for rendering geometry of future upper airway models even more realistic as the appropriate three-dimensional (3D) medical imaging techniques are becoming available.     

Tracheal stenosis: a flow dynamics study.

Patients referred for treatment of tracheal stenosis typically are asymptomatic until critical narrowing of the airway occurs, which then requires immediate intervention. To understand how tracheal stenosis affects local pressure drops and explore how a dramatic increase in pressure drop could possibly be detected at an early stage, a computational fluid dynamics (CFD) study was undertaken. We assessed flow patterns and pressure drops over tracheal stenoses artificially inserted into a realistic three-dimensional upper airway model derived from multislice computed tomography images obtained in healthy men. Solving the Navier-Stokes equations (with a Yang-shih k-epsilon turbulence model) for different degrees of tracheal constriction located approximately one tracheal diameter below the glottis, the simulated pressure drop over the stenosis (DeltaP) was seen to dramatically increase only when well over 70% of the tracheal lumen was obliterated. At 30 l/min, DeltaP increased from 7 Pa for a 50% stenosis to, respectively, 46 and 235 Pa for 80% and 90% stenosis. The pressure-flow relationship in the entire upper airway model (between mouth and end of trachea) in the flow range 0-60 l/min showed a power law relationship with best-fit flow exponent of 1.77 in the absence of stenosis. The exponent became 1.92 and 2.00 in the case of 60% and 85% constriction, respectively. The present simulations confirm that the overall pressure drop at rest is only affected in case of severe constriction, and the simulated flow dependence of pressure drop suggests a means of detecting stenosis at a precritical stage.     

Inertial deposition effects: a study of aerosol mechanics in the trachea using laser Doppler velocimetry and fluorescent dye

This study characterizes the axial velocity and axial turbulence intensity patterns noted in the tracheal portion of a cadaver-based throat model at two different steady flow rates (18.1 and 41.1 LPM.) This characterization was performed using Phase Doppler Interferometry (Laser Doppler Velocimetry). Deposition, as assessed qualitatively using fluorescent dye, is related to the position of the laryngeal jet within the trachea. The position of the jet is dependent on the downstream conditions of the model. It is proposed therefore that lung/airway conditions may have important effects on aerosol deposition within the throat. There is no correspondence noted between regions of high axial turbulence intensity and deposition.     

Three-dimensional Particle Tracking Velocimetry for Turbulence Applications: Case of a Jet Flow

3D-PTV is a quantitative flow measurement technique that aims to track the Lagrangian paths of a set of particles in three dimensions using stereoscopic recording of image sequences. The basic components, features, constraints and optimization tips of a 3D-PTV topology consisting of a high-speed camera with a four-view splitter are described and discussed in this article. The technique is applied to the intermediate flow field (5 <x/d <25) of a circular jet at Re ≈ 7,000. Lagrangian flow features and turbulence quantities in an Eulerian frame are estimated around ten diameters downstream of the jet origin and at various radial distances from the jet core. Lagrangian properties include trajectory, velocity and acceleration of selected particles as well as curvature of the flow path, which are obtained from the Frenet-Serret equation. Estimation of the 3D velocity and turbulence fields around the jet core axis at a cross-plane located at ten diameters downstream of the jet is compared with literature, and the power spectrum of the large-scale streamwise velocity motions is obtained at various radial distances from the jet core.     

Numerical 3D analysis of oscillatory flow in the time-varying laryngeal channel

Three-dimensional flow through an anatomically representative model of the human larynx has been numerically simulated. This model includes the vestibular folds, the vocal cords and the glottic and subglottic areas. Pseudo-time-varying glottic aperture and flow conditions have been considered during quiet breathing, with a peak volume flow rate of 0.75 l/s and a frequency of 0.25 Hz. Because of the severe constriction, jet-like configurations have been observed. Minor differences have been outlined between the inspiration and expiration profiles. Simulations demonstrated the presence of a backflow region which may extend to 60 mm from the glottis at peak inspiration and occupy 20% of the tracheal cross section. Because of its rolling, this backflow region appears in the sagittal plane close to the anterior wall, only one diameter from the laryngeal constriction and extends over about 40 mm. The evolution of the streamwise velocity contours and of the corresponding secondary vector plots at six critical stations, including the glottic section, has also been described. A double pair of counter-rotating vortices develops shortly downstream/upstream from the orifice respectively at inspiration/expiration and merges near the frontal plane about 25 mm from the glottis. The effect of the incoming flow has been evaluated by including the pharyngeal channel; no major difference has been observed in the computed flow patterns.     

Bronchial distribution of gas mixing in a model of the upper and central airways

Longitudinal mixing of He, O2, and sulfur hexafluoride boluses with air flowing through a three-generation tracheobronchial airway model was evaluated as the increase in volume variance of gas concentration distributions monitored at upstream and downstream sampling ports. Mixing was partitioned between tracheal and branched regions of the model at steady inspiratory and expiratory airflows from 0.044 to 0.884 l/s, both in the absence and in the presence of a removable larynx cast. During inspiration in the absence of the larynx, mixing increased as airflow increased, reaching a peak value at 0.2 l/s, and decreasing as airflow increased further. This mixing peak was higher in the branched region than in the tracheal region and was inversely related to the diffusion coefficient of the indicator gas-air mixture. During inspiration in the presence of the larynx, a mixing peak was observed in the tracheal region, but mixing peaks in the branched region were eliminated by turbulence propagated downstream from the larynx. During expiration, laryngeal turbulence was propagated far enough downstream (i.e., proximal to the trachea) that mixing peaks could be observed in both tracheal and branched regions whether or not the larynx cast was present.     

Computational simulations of airflow in an in vitro model of the pediatric upper airways

In order to understand mechanisms of gas and aerosol transport in the human respiratory system airflow in the upper airways of a pediatric subject (male aged 5) was calculated using Computational Fluid Dynamic techniques. An in vitro reconstruction of the subject's anatomy was produced from MRI images. Flow fields were solved for steady inhalation at 6.4 and 8 LPM. For validation of the numerical solution, airflow in an adult cadaver based trachea was solved using identical numerical methods. Comparisons were made between experimental results and computational data of the adult model to determine solution validity. It was found that numerical simulations can provide an accurate representation of axial velocities and turbulence intensity. Data on flow resistance, axial velocities, secondary velocity vectors, and turbulent kinetic energy are presented for the pediatric case. Turbulent kinetic energy and axial velocities were heavily dependant on flow rate, whereas turbulence intensity varied less over the flow rates studied. The laryngeal jet from an adult model was compared to the laryngeal jet in the pediatric model based on Tracheal Reynolds number. The pediatric case indicated that children show axial velocities in the laryngeal jet comparable to adults, who have much higher tracheal Reynolds numbers than children due to larger characteristic dimensions. The intensity of turbulence follows a similar trend, with higher turbulent kinetic energy levels in the pediatric model than would be expected from measurements in adults at similar tracheal Reynolds numbers. There was reasonable agreement between the location of flow structures between adults and children, suggesting that an unknown length scale correlation factor could exist that would produce acceptable predictions of pediatric velocimetry based off of adult data sets. A combined scale for turbulent intensity as well may not exist due to the complex nature of turbulence production and dissipation.     

Effects of oral airway geometry characteristics on the diffusional deposition of inhaled nanoparticles

The deposition of ultrafine aerosols in the respiratory tract presents a significant health risk due to the increased cellular-level response that these particles may invoke. However, the effects of geometric simplifications on local and regional nanoparticle depositions remain unknown for the oral airway and throughout the respiratory tract. The objective of this study is to assess the effects of geometric simplifications on diffusional transport and deposition characteristics of inhaled ultrafine aerosols in models of the extrathoracic oral airway. A realistic model of the oral airway with the nasopharynx (NP) included has been constructed based on computed tomography scans of a healthy adult in conjunction with measurements reported in the literature. Three other geometries with descending degrees of physical realism were then constructed with successive geometric simplifications of the realistic model. A validated low Reynolds number k-omega turbulence model was employed to simulate laminar, transitional, and fully turbulent flow regimes for the transport of 1-200 nm particles. Results of this study indicate that the geometric simplifications considered did not significantly affect the total deposition efficiency or maximum local deposition enhancement of nanoparticles. However, particle transport dynamics and the underlying flow characteristics such as separation, turbulence intensity, and secondary motions did show an observable sensitivity to the geometric complexity. The orientation of the upper trachea was shown to be a major factor determining local deposition downstream of the glottis and should be retained in future models of the respiratory tract. In contrast, retaining the NP produced negligible variations in airway dynamics and could be excluded for predominantly oral breathing conditions. Results of this study corroborate the use of existing diffusion correlations based on a circular oral airway model. In comparison to previous studies, an improved correlation for the deposition of nanoparticles was developed based on a wider range of particle sizes and flow rates, which captures the dependence of the Sherwood number on both Reynolds and Schmidt numbers.     

Velocity field studies at surgically imposed arterial stenoses on the abdominal aorta in pigs

In order to describe velocity profiles and the size of deterministic and non-deterministic velocity disturbances at arterial stenoses, symmetrical and asymmetrical stenoses with intended area reductions of 50% (‘moderate’) and 85% (‘severe’) were applied on the abdominal aorta in six pigs. Blood velocities were registered by hot-film anemometry in 21 measuring points distributed across the vessel cross-sectional area in one pre-stenotic and three post-stenotic positions. Signal analysis included ensemble averaging, the high-pass filtering technique, and three-dimensional visualization. None of the stenoses affected the pre-stenotic velocity field. Downstream moderate stenoses flow separation and vortex formation were present. Moderate asymmetric stenoses induced turbulence in the post-stenotic velocity field. Immediately downstream of severe stenoses a prominent post-stenotic jet was present. Farther downstream, a multitude of coherent vortices and turbulence dominated the flow field. The transverse distribution of turbulence intensity parallelled with the peak systolic velocity profile, whereas transverse profiles of the relative turbulence intensity (turbulence intensity/mean velocity) revealed peak values in flow field locations with high velocity gradients. Velocity parameters for symmetric and asymmetric severe stenoses were highly comparable. However, the exact degree of stenosis was significantly higher for symmetrical (85%) than for asymmetrical (76%) stenoses. Therefore, recalling that stenosis severity strongly influences the development of velocity disturbances, this indicates that asymmetry of a stenosis is a predictor for blood velocity disturbances.     

Microscale Gaseous Slip Flow in the Insect Trachea and Tracheoles

An analytical investigation into compressible gas flow with slight rarefactions through the insect trachea and tracheoles during the closed spiracle phase is undertaken, and a complete set of asymptotic analytical solutions is presented. We first obtain estimates of the Reynolds and Mach numbers at the channel terminal ends where the tracheoles directly deliver respiratory gases to the cells, by comparing the magnitude of the different forces in the compressible gas flow. The 2D Navier–Stokes equations with a slip boundary condition are used to investigate compressibility and rarefied effects in the trachea and tracheoles. Expressions for the velocity components, pressure gradients and net flow inside the trachea are then presented. Numerical simulations of the tracheal compressible flow are performed to validate the analytical results from this study. This work extends previous work of Arkilic et al. (J Microelectromech Syst 6(2):167–178, 1997) on compressible flows through a microchannel. Novel devices for microfluidic compressible flow transport may be invented from results obtained in this study.     

Anterior adhesive areas and adjacent secretions in the parasitic flatworms Decacotyle lymmae and D. tetrakordyle (Monogenea: Monocotylidae) from the gills of stingrays

Abstract. The monogeneans Decacotyle lymmae and D. tetrakordyle (Monocotylidae: Decacotylinae), from gills of the dasyatid stingrays Taeniura lymma and Pastinachus sephen , respectively, have a single aperture for adhesive secretion on each side of the anterior ventrolateral region. Rod-shaped bodies (S1) and electron-dense spherical secretion (S2) exit through specialised ducts opening adjacent to one another within these apertures. The S1 bodies are 230 ± 11 nm wide and ≥4 μm long in D. lymmae and 240 ± 9 nm wide and ≥3.3 μm long in D. tetrakordyle . The S2 bodies have a diameter of 88 ± 7 nm in D. lymmae and 65 ± 6 nm in D. tetrakordyle . The apertures are unusual in being extremely small (internal diameter, 3–5 μm). Each aperture has a slit-like surface opening as small as 160 nm wide, surrounded by muscle fibres indicating that they may be opened and closed. The aperture is also surrounded and underlain by muscle fibres that may aid in secretion from, or even eversion of, the tissue within the aperture. Sensilla/cilia are also found within the apertures. Additional secretions from anteromedian and anterolateral glands (body glands), each containing granular secretions, occur in profusion and exit anteriorly and posteriorly to the position of the apertures, through duct openings in the general body tegument. These granular secretions do not appear to be associated with anterior adhesion. Both species show similarities in aperture, underlying tissue, sense organ, and secretion detail, in accordance with findings from other monogenean genera, and which supports the importance of such data for phylogenetic studies.     

Respiratory flow in a realistic tracheostenosis model

The possible mechanism of wheeze generation in tracheostenosis was identified by measuring inspiratory and expiratory flow in a "morphological and distensible" realistic tracheostenosis model. The shape of the model was based on CT (Computed Tomography) images of a patient that had tracheostenosis. A trachea consists of tracheal cartilage rings and smooth muscle. Spatial variation of wall distensibility was achieved in the model by varying the wall thickness based on the elastic modulus measured in pig airways. The spatial variation influenced the flow in the airway and the turbulence production rate decreased faster at smooth muscles. Using the model, we investigated the mechanism of wheeze generation by focusing on the turbulence intensity. The turbulence intensity in expiratory flow was about twice that in inspiratory flow, and larger vortices existed in post-stenosis in expiratory flow, and thus might contribute to wheeze generation.     

Simulation of swallowing dysfunction and mechanical ventilation after a Montgomery T-tube insertion

The Montgomery T-tube is used as a combined tracheal stent and airway after laryngotracheoplasty, to keep the lumen open and prevent mucosal laceration from scarring. It is valuable in the management of upper and mid-tracheal lesions, while invaluable in long and multisegmental stenting lesions. Numerical simulations based on real-patient-tracheal geometry, experimental tissue characterization, and previous numerical estimation of the physiological swallowing force are performed to estimate the consequences of Montgomery T-tube implantation on swallowing and assisted ventilation: structural analysis of swallowing is performed to evaluate patient swallowing capacity, and computational fluid dynamics simulation is carried out to analyze related mechanical ventilation. With an inserted Montgomery T-tube, vertical displacement (Z-axis) reaches 8.01 mm, whereas in the Y-axis, it reaches 6.63 mm. The maximal principal stress obtained during swallowing was 1.6 MPa surrounding the hole and in the upper contact with the tracheal wall. Fluid flow simulation of the mechanical ventilation revealed positive pressure for both inhalation and exhalation, being higher for inspiration. The muscular deflections, considerable during normal breathing, are nonphysiological, and this aspect results in a constant overload of the tracheal muscle. During swallowing, the trachea ascends producing a nonhomogeneous elongation. This movement can be compromised when prosthesis is inserted, which explains the high incidence of glottis close inefficiency. Fluid simulations showed that nonphysiological pressure is established inside the trachea due to mechanical ventilation. This may lead to an overload of the tracheal muscle, explaining several related problems as muscle thinning or decrease in contractile function.     

Numerical investigation of the three-dimensional flow in a human lung model

The flow field at inspiration and expiration in the upper human airways consisting of the trachea down to the sixth generation of the bronchial tree is numerically simulated. The three-dimensional steady flow at a hydraulic diameter-based Reynolds number Re(D)=1250 is computed via a lattice-Boltzmann method (LBM). The simulation is validated by the experimental data based on particle-image velocimetry (PIV) measurements. The good agreement between numerical and experimental results is evidenced by comparing velocity contours and distributions in a defined reference plane. The results show the LBM to be an accurate tool to numerically predict flow structures in the human lung. Using an automatic Cartesian grid generator, the overall process time from meshing to a steady-state solution is <12h. Moreover, the numerical simulation allows a closer analysis of the secondary flow structures than in the experimental investigation. The three-dimensional streamline patterns reveal some insight on the air exchange mechanism at inspiration and expiration. At inspiration, the slower near-wall tracheal flow enters through the right principal bronchus into the right upper lobar bronchus. The bulk mass flux in the trachea is nearly evenly distributed over the left upper, center and lower lobar bronchi and the right center and lower bronchi. At expiration, the air from the right upper lobar bronchus enters the right center of the trachea and displaces the airflow from the lower and center right bronchi such that the tracheal positions of the streamlines at inspiration and expiration are switched. The flow in the left bronchi does not show this kind of switching. The findings emphasize the impact of the asymmetry of the lung geometry on the respiratory air exchange mechanism.     

Velocity and turbulence distribution around lotic macrophytes

Flow velocity and turbulence patterns were measured in and around a common lotic macrophyte, Ranunculus penicillatus subspecies pseudofluitans (stream water-crowfoot), using a two-dimensional electromagnetic current meter (EMCM). Due to the high shooting density of this species, there was a sharp velocity gradient at the plant boundary, with velocities dropping to a constantly low value after no more than 5 cm into the plant, thus forcing most of the flow over and around the macrophyte. There was a dead-water zone immediately downstream of the plant, beyond which the current moved in from the sides to allow flow under the trailing shoots. High turbulence intensities were recorded for both downstream and cross-stream velocity components at the lateral margins and downstream of the plant. Meanwhile, pulses of water upstream of the plant produced turbulence in the downstream component, but not in the cross-stream component.     

Role of the epithelial layer in the generation of superoxide anion by the guinea-pig isolated trachea.

The lucigenin-dependent chemiluminescence generation by guinea-pig isolated tracheal two rings preparations was studied. Tracheal preparations stimulated with phorbol myristate acetate (PMA) or opsonized zymosan generated chemiluminescence. The total amount of chemiluminescence generated in 120 min was 754+/-63 mV x min for PMA and 4832+/-396 mV x min for zymosan. Generation of chemiluminescence was decreased by more than 50% when the tissues were co-incubated with superoxide dismutase (100 U/ml). Also, addition of direct donors of nitric oxide diminished chemiluminescence generation by zymosan-activated tracheal rings significantly by about 50%. However, the presence of the precursor or of inhibitors of nitric oxide synthase did not influence zymosan-induced chemiluminescence. Removal of the epithelial layer from tracheal rings caused an approximately 90% decrease in chemiluminescence response. However, isolated epithelial cell suspensions did not generate chemiluminescence. Histologic examination showed that the number of eosinophils in the tracheal tissue was reduced from 56+/-7 to 18+/-8 per mm basal membrane when the epithelial layer was removed. These results indicated that (1) superoxide anion formation can take place in the guinea-pig trachea, (2) eosinophils in the epithelial and submucosal layers of guinea-pig trachea are likely candidates for superoxide generation although other cell types can also be involved, and (3) besides relaxing airway smooth muscle, nitric oxide donors may also affect superoxide in the airways.     

Spectrum analysis of turbulence in the canine ascending aorta measured with a hot-film anemometer

We measured turbulence velocity in the canine ascending aorta using a hot-film anemometer. Blood flow velocity was measured at various points across the ascending aorta approximately 1.5-2 times the diameter downstream from the aortic valve. The turbulence spectrum was calculated and its characteristics were examined in connection with the mean Reynolds number and/or measuring positions. In the higher wave number range the values of the turbulence spectra were higher at larger mean Reynolds number. In the higher wave number range, the values of the turbulence spectra were higher at points closer to the centerline of the aorta, when the mean Reynolds number was relatively large. The patterns of the turbulence spectra at various points outside the boundary layer on the aortic wall were similar.     

Particle deposition in a CT-scanned human lung airway

The particle deposition in a computerized tomography (CT)-scanned human lung was numerically investigated. The five-generation airway is extracted from the trachea to segmental bronchi of a 60-year-old Chinese male patient. Computations were carried out in the flow rate range of 210–630 ml/s (Reynolds number range of 1000–3000) and particle size of 2–10 μm (Stokes number range of 0.0007–0.049). To count the effect of laryngeal jet on trachea inlet, the trachea was extended and modified to simulate the larynx, consequently the inlet velocity profile is biased towards the rear wall. The laryngeal jet-induced turbulence was simulated using low Reynolds number (LRN) κ–ω turbulent model. Particle deposition patterns, deposition efficiency and deposition factor were studied in detail. The turbulent flow has significant effect on the particle deposition, and the present deposition factor is compared well with the available data.