Characterization of the Vertical Stiffness Gradient in Cadaveric Human and Excised Canine Larynges.

The elastic properties of the folds govern the characteristics of vocal fold vibrations. This study addresses existing gaps by investigating the Young's modulus along the anterior-posterior direction in excised canine and cadaveric human vocal folds. Micro-indentation testing was conducted on six excised canines and three cadaveric human larynges. Multiple points along the medial glottal wall were indented to determine force-displacement, stress-strain relationships, and Young's modulus as a function of Green's strain. A vertical stiffness gradient was consistently observed in both canine and human samples, with higher stiffness in the inferior aspect compared with the superior aspect. The stiffness increased toward both the anterior and posterior directions from the mid-coronal plane, with a more pronounced increase at the posterior edge. Human vocal folds generally exhibited lower stiffness at low strains but were comparable to canine vocal folds at higher strains. These findings suggest that the canine larynx model is a reasonable representation of the human laryngeal tissues' elastic property trends. This analysis of the vertical stiffness gradient in canine and human vocal folds provides valuable data for improving experimental and numerical models of phonation.

Comparison of Aerodynamic and Elastic Properties in Tissue and Synthetic Models of Vocal Fold Vibrations.

Three laryngeal models were used to investigate the aerodynamic and elastic properties of vocal fold vibration: cadaveric human, excised canine, and synthetic silicone vocal folds. The aim was to compare the characteristics of these models to enhance our understanding of phonatory mechanisms. Flow and medial glottal wall geometry were acquired via particle image velocimetry. Elastic properties were assessed from force-displacement tests. Relatively, the human larynges had higher fundamental frequency values, while canine and synthetic models exhibited greater flow rates. Canine models demonstrated the highest divergence angles and vertical stiffness gradients followed by the human model, both displaying flow separation vortices during closing. Synthetic models, whose advantage is their accessibility and repeatability, displayed the lowest glottal divergence angles and total circulation values compared to tissue models with no flow separation vortices. The elasticity tests revealed that tissue models showed significant hysteresis and vertical stiffness gradients, unlike the synthetic models. These results underscore the importance of model selection based on specific research needs and highlight the potential of canine and synthetic models for controlled experimental studies in phonation.

Acoustics and aerodynamic effects following glottal and infraglottal medialization in an excised larynx model.

OBJECTIVE: This study aimed to investigate the impact of the implant's vertical location during Type 1 Thyroplasty (T1T) on acoustics and glottal aerodynamics using excised canine larynx model, providing insights into the optimal technique for treating unilateral vocal fold paralysis (UVFP). METHODS: Measurements were conducted in six excised canine larynges using Silastic implants. Two implant locations, glottal and infraglottal, were tested for each larynx at low and high subglottal pressure levels. Acoustic and intraglottal flow velocity field measurements were taken to assess vocal efficiency (VE), cepstral peak prominence (CPP), and the development of intraglottal vortices. RESULTS: The results indicated that the implant's vertical location significantly influenced vocal efficiency (p = 0.045), with the infraglottal implant generally yielding higher VE values. The effect on CPP was not statistically significant (p = 0.234). Intraglottal velocity field measurements demonstrated larger glottal divergence angles and stronger vortices with the infraglottal implant. CONCLUSION: The findings suggest that medializing the paralyzed fold at the infraglottal level rather than the glottal level can lead to improved vocal efficiency. The observed larger divergence angles and stronger intraglottal vortices with infraglottal medialization may enhance voice outcomes in UVFP patients. These findings have important implications for optimizing T1T procedures and improving voice quality in individuals with UVFP. Further research is warranted to validate these results in clinical settings.

The Effect of an Increasing Subglottal Stenosis Constriction That Extends From the Vocal Folds to the Inferior Border of the Cricoid Cartilage.

Acquired subglottal stenosis is an unpredicted complication that can occur in some patients who have undergone prolonged endotracheal intubation. It is a narrowing of the airway at the level of the cricoid cartilage that can restrict airflow and cause breathing difficulty. Stenosis is typically treated with endoscopic airway dilation, with some patients experiencing multiple recurrences. The study highlights the potential of computational fluid dynamics as a noninvasive method for monitoring subglottic stenosis, which can aid in early diagnosis and surgical planning. An anatomically accurate human laryngeal airway model was constructed from computerized tomography (CT) scans. The subglottis cross-sectional area was narrowed systematically using ≈10% decrements. A quadratic profile was used to interpolate the transformation of the airway geometry from its modified shape to the baseline geometry. The numerical results were validated by static pressure measurements conducted in a physical model. The results show that airway resistance follows a squared ratio that is inversely proportional to the size of the subglottal opening (R∝A-2). The study found that critical constriction occurs in the subglottal region at 70% stenosis (upper end of grade 2). Moreover, removing airway tissue below 40% stenosis during surgical intervention does not significantly decrease airway resistance.

The Effects of Negative Pressure Induced by Flow Separation Vortices on Vocal Fold Dynamics during Voice Production.

This study used a two-dimensional flow-structure-interaction computer model to investigate the effects of flow-separation-vortex-induced negative pressure on vocal fold vibration and flow dynamics during vocal fold vibration. The study found that negative pressure induced by flow separation vortices enhances vocal fold vibration by increasing aeroelastic energy transfer during vibration. The result showed that the intraglottal pressure was predominantly negative after flow separation before gradually recovering to zero at the glottis exit. When the negative pressure was removed, the vibration amplitude and flow rate were reduced by up to 20%, and the closing speed, flow skewness quotient, and maximum flow declination rate were reduced by up to 40%. The study provides insights into the complex interactions between flow dynamics, vocal fold vibration, and energy transfer during voice production.

Gels That Serve as Mucus Simulants: A Review.

Mucus is a critical part of the human body's immune system that traps and carries away various particulates such as anthropogenic pollutants, pollen, viruses, etc. Various synthetic hydrogels have been developed to mimic mucus, using different polymers as their backbones. Common to these simulants is a three-dimensional gel network that is physically crosslinked and is capable of loosely entrapping water within. Two of the challenges in mimicking mucus using synthetic hydrogels include the need to mimic the rheological properties of the mucus and its ability to capture particulates (its adhesion mechanism). In this paper, we review the existing mucus simulants and discuss their rheological, adhesive, and tribological properties. We show that most, but not all, simulants indeed mimic the rheological properties of the mucus; like mucus, most hydrogel mucus simulants reviewed here demonstrated a higher storage modulus than its loss modulus, and their values are in the range of that found in mucus. However, only one mimics the adhesive properties of the mucus (which are critical for the ability of mucus to capture particulates), Polyvinyl alcohol-Borax hydrogel.

Influence of aortic valve morphology on vortical structures and wall shear stress.

The aim of this paper is to assess the association between valve morphology and vortical structures quantitatively and to highlight the influence of valve morphology/orientation on aorta's susceptibility to shear stress, both proximal and distal. Four-dimensional phase-contrast magnetic resonance imaging (4D PCMRI) data of 6 subjects, 3 with tricuspid aortic valve (TAV) and 3 with functionally bicuspid aortic values (BAV) with right-left coronary leaflet fusion, were processed and analyzed for vorticity and wall shear stress trends. Computational fluid dynamics (CFD) has been used with moving TAV and BAV valve designs in patient-specific aortae to compare with in vivo shear stress data. Vorticity from 4D PCMRI data about the aortic centerline demonstrated that TAVs had a higher number of vortical flow structures than BAVs at peak systole. Coalescing of flow structures was shown to be possible in the arch region of all subjects. Wall shear stress (WSS) distribution from CFD results at the aortic root is predominantly symmetric for TAVs but highly asymmetric for BAVs with the region opposite the raphe (fusion location of underdeveloped leaflets) being subjected to higher WSS. Asymmetry in the size and number of leaflets in BAVs and TAVs significantly influence vortical structures and WSS in the proximal aorta for all valve types and distal aorta for certain valve orientations of BAV. Analysis of vortical structures using 4D PCMRI data (on the left side) and wall shear stress data using CFD (on the right side).

An Ex-vivo Model Examining Acoustics and Aerodynamic Effects Following Medialization With and Without Arytenoid Adduction.

OBJECTIVES/HYPOTHESIS: Quantify differences in acoustics and intraglottal flow fields between Thyroplasty Type 1 (TT1) with and without arytenoid adduction (AA) using excised canine larynx model. STUDY DESIGN: Basic science experiments using excised larynges. METHODS: Surgical procedures were implemented in eight excised canine larynges. Acoustics and intraglottal flow measurements were taken at low and high subglottal pressures in each experimental setup. RESULTS: In all larynges, vocal efficiency (VE) and cepstrum peak prominence (CPP) were higher, and the mean phonatory flow rate was lower in TT1 with AA than without AA. The glottal asymmetry is reduced with AA and promotes the formation of stronger vortices in the glottal flow during the closing phase of the vibrating folds. CONCLUSIONS: Findings suggest a clear acoustic and aerodynamic benefit to the addition of AA when performing TT1. It shows significant improvement in CPP, translating to decreased breathiness and dysphonia and increased VE, leading to easier and more sustainable phonation. Stronger intraglottal vortices are known to be correlated with the loudness of voice produced by phonation. LEVEL OF EVIDENCE: N/A Laryngoscope, 133:621-627, 2023.

How Face Masks Affect Acoustic and Auditory Perceptual Characteristics of the Singing Voice.

Wearing a face mask has been accepted as one of the most effective ways for slowing the spread of COVID-19. Yet information regarding the degree to which masks affect acoustics and perception associated with voice performers is scarce. This study examines these effects with common face masks, namely a neck gaiter, disposable surgical mask, and N95 mask, as well as a novel material that could be used as a mask (acoustic foam). A recorded excerpt from the "Star-Spangled Banner" was played through a miniature speaker placed inside the mouth of a masked manikin. Experienced listeners were asked to rate perceptual qualities of these singing stimuli by blindly comparing them with the same recording captured without a mask. Acoustic analysis showed that face masks affected the sound by enhancing or suppressing different frequency bands compared to no mask. Acoustic energy around the singer's formant was reduced when using surgical and N95 masks, which matches observations that these masks are more detrimental to the perceptions of singing voice compared with neck gaiter or acoustic foam. It suggests that singers can benefit from masks designed for minimal impact on auditory perception of the singing voice while maintaining reasonable efficacy of filtering efficiency.

Effect of Blood Transfusion on Cerebral Hemodynamics and Vascular Topology Described by Computational Fluid Dynamics in Sickle Cell Disease Patients.

The main objective of this study was to demonstrate that computational fluid dynamics (CFD) modeling can be used to study the contribution of covert and overt vascular architecture to the risk for cerebrovascular disease in sickle cell disease (SCD) and to determine the mechanisms of response to therapy such as chronic red blood cell (cRBC) transfusions. We analyzed baseline (screening), pre-randomization and study exit magnetic resonance angiogram (MRA) images from 10 (5 each from the transfusion and observation arms) pediatric sickle SCD participants in the silent cerebral infarct transfusion (SIT) trial using CFD modeling. We reconstructed the intracranial portion of the internal carotid artery and branches and extracted the geometry using 3D Slicer. We cut specific portions of the large intracranial artery to include segments of the internal carotid, middle, anterior, and posterior cerebral arteries such that the vessel segment analyzed extended from the intracranial beginning of the internal carotid artery up to immediately after (~0.25 inches) the middle cerebral artery branching point. Cut models were imported into Ansys 2021R2/2022R1 and laminar and time-dependent flow simulation was performed. Change in time averaged mean velocity, wall shear stress, and vessel tortuosity were compared between the observation and cRBC arms. We did not observe a correlation between time averaged mean velocity (TAMV) and mean transcranial Doppler (TCD) velocity at study entry. There was also no difference in change in time average mean velocity, wall shear stress (WSS), and vessel tortuosity between the observation and cRBC transfusion arms. WSS and TAMV were abnormal for 2 (developed TIA) out of the 3 participants (one participant had silent cerebral infarctions) that developed neurovascular outcomes. CFD approaches allow for the evaluation of vascular topology and hemodynamics in SCD using MRA images. In this proof of principle study, we show that CFD could be a useful tool and we intend to carry out future studies with a larger sample to enable more robust conclusions.

Demonstration of mucus simulant clearance in a Bench-Model using acoustic Field-Integrated Intrapulmonary Percussive ventilation.

Intrapulmonary Percussive Ventilation (IPV) is a high-frequency airway clearance technique used to help in mucus transport for mechanically ventilated and unventilated patients. Despite the many years of usage, this technique does not provide clear evidence of its intended efficacy. This is mainly attributable to the lack of in vitro observations that show "mucokinesis" towards the direction of the mouth. In the current manuscript, we demonstrate and subsequently propose a mechanism that details the movement of a mucus simulant in the proximal (towards the mouthpiece) direction. Towards this end, a novel method utilizing a high-frequency acoustic field in addition to the conventional air pulsations brought forth by traditional IPV is proposed. Under these conditions, at certain parameter settings, it is shown that the simulant is broken down into much smaller parts and subsequently pushed in the upstream direction gradually over a period of half-hour.

Fluid-Structure Interaction Analysis of Aerodynamic and Elasticity Forces During Vocal Fold Vibration.

The effect of the intraglottal vortices on the glottal flow waveform was explored using flow-structure-interaction (FSI) modeling. These vortices form near the superior aspect of the vocal folds during the closing phase of the folds' vibration. The geometry of the vocal fold was based on the well-known M5 model. The model did not include a vocal tract to remove its inertance effect on the glottal flow. Material properties for the cover and body layers of the folds were set using curve fit to experimental data of tissue elasticity. A commercially available FSI solver was used to perform simulations at low and high values of subglottal input pressure. Validation of the FSI results showed a good agreement for the glottal flow and the vocal fold displacement data with measurements taken in the excised canine larynx model. The simulations result further support the hypothesis that intraglottal vortices can affect the glottal flow waveform, specifically its maximum flow declination rate (MFDR). It showed that MFDR occurs at the same phase when the highest intraglottal vortical strength and the negative pressure occur. It also showed that when MFDR occurs, the magnitude of the aerodynamic force acting on the glottal wall is greater than the elastic recoil force predicted in the tissue. These findings are significant because nearly all theoretical and computational models that study the vocal fold vibrations mechanism do not consider the intraglottal negative pressure caused by the vortices as an additional closing force acting on the folds.

Predicting critical closing pressure in children with obstructive sleep apnea using fluid-structure interaction.

Surgical treatment of obstructive sleep apnea (OSA) in children requires knowledge of upper airway dynamics, including the closing pressure (Pcrit), a measure of airway collapsibility. We applied a flow-structure interaction (FSI) computational model to estimate Pcrit in patient-specific upper airway models obtained from magnetic resonance imaging (MRI) scans. We sought to examine the agreement between measured and estimated Pcrit from FSI models in children with Down syndrome. We hypothesized that the estimated Pcrit would accurately reflect measured Pcrit during sleep and therefore reflect the severity of OSA as measured by the obstructive apnea-hypopnea index (AHI). All participants (n = 41) underwent polysomnography and sedated sleep MRI scans. We used Bland-Altman plots to examine the agreement between measured and estimated Pcrit. We determined associations between estimated Pcrit and OSA severity, as measured by AHI, using regression models. The agreement between passive and estimated Pcrit showed a fixed bias of -1.31 [confidence interval (CI) = -2.78, 0.15] and a nonsignificant proportional bias. A weaker agreement with active Pcrit was observed. A model including AHI, gender, an interaction term for AHI, and gender and neck circumference explained the largest variation (R2 = 0.61) in the relationship between AHI and estimated Pcrit (P < 0.0001). Overlap between the areas of the airway with the lowest stiffness, and areas of collapse on dynamic MRI, was 77.4 ± 30% for the nasopharyngeal region and 78.6 ± 33% for the retroglossal region. The agreement between measured and estimated Pcrit and the significant association with AHI supports the validity of Pcrit estimates from the FSI model.NEW & NOTEWORTHY We present a noninvasive method for estimating critical closing pressure (Pcrit) using fluid-structure interaction (FSI) simulations and magnetic resonance imaging (MRI) scans in patients with obstructive sleep apnea (OSA). We used patient-specific stiffness measures in our FSI model to account for any individual variability in the elasticity of soft tissues surrounding the upper airway. We validated this model by measuring the degree of agreement between measured and estimated Pcrit.

Evaluating the biomechanical characteristics of cuffed-tracheostomy tubes using finite element analysis.

The objective of this study was to perform finite element analysis (FEA) of cuff inflation within an anatomically accurate model of an adult trachea in four different cuffed-tracheostomy tube designs. The leakage quantified by the distance between the cuff and trachea was largest for the Tracoe cuff and smallest for the Portex cuff. The smooth muscle stresses were greatest for the Portex and least for the Distal cuff, respectively. The proposed FEA model offers a promising approach to virtually evaluate the sealing efficacy of cuffed-tracheostomy tubes and the tracheal wall stresses induced by cuff inflation, prior to application.

Influence of Material Model and Aortic Root Motion in Finite Element Analysis of Two Exemplary Cases of Proximal Aortic Dissection.

The risk of type-A dissection is increased in subjects with connective tissue disorders and dilatation of the proximal aorta. The location and extents of vessel wall tears in these patients could be potentially missed during prospective imaging studies. The objective of this study is to estimate the distribution of systolic wall stress in two exemplary cases of proximal dissection using finite element analysis (FEA) and evaluate the sensitivity of the distribution to the choice of anisotropic material model and root motion. FEA was performed for predissection aortas, without prior knowledge of the origin and extents of vessel wall tear. The stress distribution was evaluated along the wall tear in the postdissection aortas. The stress distribution was compared for the Fung and Holzapfel models with and without root motion. For the subject with spiral dissection, peak stress coincided with the origin of the tear in the sinotubular junction. For the case with root dissection, maximum stress was obtained at the distal end of the tear. The FEA predicted tear pressure was 20% higher for the subject with root dissection as compared to the case with spiral dissection. The predicted tear pressure was higher (9-11%) for root motions up to 10 mm. The Holzapfel model predicted a tear pressure that was lower (8-15%) than the Fung model. The FEA results showed that both material response and root motion could potentially influence the predicted dissection pressure of the proximal aorta at least for conditions tested in this study.

Quantification of the Intraglottal Pressure Induced by Flow Separation Vortices Using Large Eddy Simulation.

The greatest rate of change in the glottal flow rate during phonation is a rapid decrease that occurs during the latter part of the glottal closing. Previous works showed that intraglottal flow separation vortices form in a divergent glottis, produce negative gauge pressures (below atmospheric) during closing. It is hypothesized here that flow separation vortices contribute to the rapid closing mechanism of the true vocal folds during phonation. Four idealized static models (M5) of the human larynx were investigated using large eddy simulation: 2 models featured parallel folds that did not enable flow separation in the glottis and 2 models involved a divergent glottis. The influence of the ventricular gap (narrow/wide) is evaluated. An unsteady pressure inlet representing a voicing cycle was applied to the sub-glottal region to mimic the time-varying glottal flow. Intraglottal vortex structures formed downstream of the separation point in a divergent glottis. Their existence caused a higher closing force that was applied onto the vocal folds. A narrow ventricular gap strengthens this effect. Strength of the intraglottal vortices increased with the maximum flow declination rate. Therefore, a more divergent shape of the glottis during glottal closing will be one of the main contributors to the voice quality.

Effects of False Vocal Folds on Intraglottal Velocity Fields.

Previous models have theorized that, during phonation, skewing of the glottal waveform (which is correlated with acoustic intensity) occurred because of inertance of the vocal tract. Later, we reported that skewing of the flow rate waveform can occur without the presence of a vocal tract in an excised canine larynx. We hypothesized that in the absence of a vocal tract, the skewing formed when dynamic pressures acted on the glottal wall during the closing phase; such pressures were greatly affected by formation of intraglottal vortices. In this study, we aim to identify how changes in false vocal folds constriction can affect the acoustics and intraglottal flow dynamics. The intraglottal flow measurements were made using particle image velocimetry in an excised canine larynx where a vocal tract model was placed above the larynx and the constriction between the false vocal folds was varied. Our results show that for similar values of subglottal pressures, the skewing of the glottal waveform, strength of the intraglottal vortices, and acoustic energy increased as the constriction between the false vocal folds was increased. These preliminary findings suggest that acoustic intensity during phonation can be increased by the addition of a vocal tract with false fold constriction.

Impact of Vertical Stiffness Gradient on the Maximum Divergence Angle.

INTRODUCTION: During vocal fold vibration, the medial surface of both folds forms a convergent shape during opening and a divergent shape during closing. A greater maximum divergence angle is associated with greater closing forces which will increase the closing speed of the glottis. An increased closing speed results in a greater acoustic intensity and greater vocal efficiency. Indentation testing showed that as the strain increases, the inferior aspect of the folds becomes stiffer than the superior aspect, resulting in the vertical stiffness gradient (VSG). We hypothesize that a reduction of the vertical stiffness gradient will reduce the maximum divergence angle. METHODS: Four excised canine larynges were tested. Stress-strain curves of the superior and inferior aspects of the fold in the mid membranous plane of the baseline larynges were taken using the indentation method. Calcium hydroxylapatite (CaHA) crystals were then injected into the superior aspect of the fold. The stress-strain tests were repeated. Particle imaging velocimetry (PIV) of the intraglottal velocity fields was performed in three larynges at different subglottal pressures in the mid coronal plane for the baseline and CaHA-injected larynges. RESULTS: CaHA injection reduced the inferior-superior stiffness gradient in all larynges. The maximal divergence angle was markedly reduced. In some cases, there was not a divergent angle. DISCUSSION: Marked reduction of the vertical stiffness gradient significantly reduces the maximum divergence angle. Clinical implications will be discussed. LEVEL OF EVIDENCE: NA Laryngoscope, 131:E1934-E1940, 2021.

Medial Surface Dynamics as a Function of Subglottal Pressure in a Canine Larynx Model.

During vocal fold vibration, there may be a mucosal wave in the superior-inferior (vertical) direction, resulting in a convergent shape during opening and a divergent shape during closing. Most of our understanding of the converging/diverging shape of the glottis has come from studies in a hemilarynx model. Previous work has shown that vibratory patterns in the full excised larynx are different than the hemilarynx. This study characterized the dynamics of the medial glottal wall geometry during vibrations in the full excised canine larynx model. Using particle image velocimetry, the intraglottal geometry was measured at the midmembranous coronal plane in an excised canine larynx model. Measurements of the glottal area were taken simultaneously using high-speed imaging. The results show that skewing of the glottal area waveform occurs without the presence of a vocal tract and that the phase-lag of the superior edge relative to the inferior edge is smaller than reported and depends on the subglottal pressure. In addition, it shows that the glottal divergence angle during closing is proportional to the magnitude of the acoustic intensity and the intraglottal negative pressure. This preliminary data suggests that more studies are needed to determine the important mechanisms determining the relationship between intraglottal flow, intraglottal geometry, and acoustics.

Hemodynamics and tissue biomechanics of the thoracic aorta with a trileaflet aortic valve at different phases of valve opening.

In a normal cardiac cycle, the trileaflet aortic valve opening is progressive, which correlates with the phasic blood flow. Therefore, we aimed to determine the impact of including an anatomically accurate reconstructed trileaflet aortic valve within a fluid-structure interaction (FSI) simulation model and determine the cyclical hemodynamic forces imposed on the thoracic aortic walls from aortic valve opening to closure. A pediatric patient with a normal trileaflet valve was recruited. Using the Cardiac Magnetic Resonance Data (CMR), a 3D model of the aortic valve and thoracic aorta was reconstructed. FSI simulations were employed to assess the tissue stress during a cardiac cycle as the result of changes in the valve opening. The blood flow was simulated as a mixture of blood plasma and red blood cells to account for non-Newtonian effects. The computation was validated with phase-contrast CMR. Windkessel boundary conditions were employed to ensure physiological pressures during the cardiac cycle. The leaflets' dynamic motion during the cardiac cycle was defined with an analytic grid velocity function. At the beginning of the valve opening a thin jet is developing. From mid-open towards full opening the stress level increases where the jet impinges the convex wall. At peak systole two counter-rotating Dean-like vortex cores manifest in the ascending aorta, which correlates with increased integrated mean stress levels. An accurate trileaflet aortic valve is needed for capturing of both primary and secondary flow features that impact the forces on the thoracic aorta wall. Omitting the aortic valve underestimates the biomechanical response.