Human phonation involves the flow-induced vibrations of the vocal folds (VFs) that result from the interaction with airflow through the larynx. Most voice dysfunctions correspond with the fluid–structure interaction (FSI) features as well as the local changes in perfusion within the VF tissue. This study aims to develop a multiphysics computational framework to simulate the interstitial fluid flow dynamics in vibrating VFs using a biphasic description of the tissue and FSI methodology. The integration of FSI and a permeable VF model presents a novel approach to capture phonation physics' complexity and investigate VF tissue's porous nature. The glottal airflow is modeled by the unsteady, incompressible Navier–Stokes equations, and the Brinkman equation is employed to simulate the flow through the saturated porous medium of the VFs. The computational model provides a prediction of tissue deformation metrics and pulsatile glottal flow, in addition to the interstitial fluid velocity and flow circulation within the porous structure. Furthermore, the model is used to characterize the effects of variation in subglottal lung pressure and VF permeability coefficient by conducting parametric studies. Subsequent investigations to quantify the relationships between these input variables, flow perfusion, pore pressure, and vibration amplitude are presented. A linear relationship is found between the vibration amplitude, pore pressure, and filtration flow with subglottal pressure, whereas a nonlinear dependence between the filtration velocity and VF permeability coefficient is detected. The outcomes highlight the importance of poroelasticity in phonation models.
more »
« less
Gender in human phonation: Fluid–structure interaction and vocal fold morphology
This paper aims to examine the effects of variations in the vocal fold (VF) morphological features associated with gender on glottal aerodynamics and tissue deformation. Nine three-dimensional geometries of the VFs in the larynx are created with various VF lengths, thicknesses, and depths to perform a parametric analysis according to gender-related geometrical parameters. The computational model is incorporated in a fluid–structure interaction methodology by adopting the transient Navier–Stokes equations to model airflow through the larynx and considering a linear elasticity model for VF dynamics. The model predictions, such as aerodynamic data through the larynx, glottal airflow, and VF deformations, are analyzed. The comparison of the simulation results for the nine cases supports the hypothesis that gender differences in laryngeal dimensions remarkably influence the glottal airflow and deformation of the VFs. Decreasing VF thickness and increasing its length corresponds to a noticeable increase in maximum tissue displacement, while variations in depth affect the flow rate significantly in the small and large larynges. Conversely, we observed that the pressure drop at the glottis is nearly independent of the VF length. A comparison of the glottal area with published imaging data illustrated a direct correlation between the glottal configuration and the morphology of the VFs.
more »
« less
- Award ID(s):
- 2138225
- PAR ID:
- 10472770
- Publisher / Repository:
- AIP Publishing
- Date Published:
- Journal Name:
- Physics of Fluids
- Volume:
- 35
- Issue:
- 4
- ISSN:
- 1070-6631
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Jetté, Marie (Ed.)Current research approaches employ traditional tissue engineering strategies to promote vocal fold (VF) tissue regeneration, whereas recent novel advances seek to use principles of developmental biology to guide tissue generation by mimicking native developmental cues, causing tissue or allogenic/autologous progenitor cells to undergo the regeneration process. To address the paucity of data to direct VF differentiation and subsequent new tissue formation, we characterize structure-proliferation relationships and tissue elastic moduli over embryonic development using a murine model. Growth, cell proliferation, and tissue biomechanics were taken at E13.5, E15.5, E16.5, E18.5, P0, and adult time points. Quadratic growth patterns were found in larynx length, maximum transverse diameter, outer dorsoventral diameter, and VF thickness; internal VF length was found to mature linearly. Cell proliferation measured with EdU in the coronal and transverse planes of the VFs was found to decrease with increasing age. Exploiting atomic force microscopy, we measured significant differences in tissue stiffness across all time points except between E13.5 and E15.5. Taken together, our results indicate that as the VF mature and develop quadratically, there is a concomitant tissue stiffness increase. Greater gains in biomechanical stiffness at later prenatal stages, correlated with reduced cell proliferation, suggest that extracellular matrix deposition may be responsible for VF thickening and increased biomechanical function, and that the onset of biomechanical loading (breathing) may also contribute to increased stiffness. These data provide a profile of VF biomechanical and growth properties that can guide the development of biomechanically-relevant scaffolds and progenitor cell differentiation for VF tissue regeneration.more » « less
-
Abstract Extracellular matrix (ECM) is a complex structure composed of bioactive molecules representative of the local tissue microenvironment. Decellularized ECM biomaterials harness these biomolecules for regenerative medicine applications. One potential therapeutic application is the use of vocal fold (VF) specific ECM to restore the VFs after injury. ECM scaffolds are derived through a process of decellularization, which aims to remove unwanted immunogenic biomolecules (e.g. DNA) while preserving the composition of the ECM. The effectiveness of the decellularization is typically assessed at the end by quantifying ECM attributes such as final dsDNA content. However, batch-to-batch variability in ECM manufacturing remains a significant challenge for the standardization, cost-effectiveness, and scale-up process. The limited number of tools available for in-process control heavily restricts the uncovering of the correlations between decellularization process parameters and ECM attributes. In this study, we developed a technique applicable to both the classical batch method and semi-continuous decellularization systems to trace the decellularization of two laryngeal tissues in real-time. We hypothesize that monitoring the bioreactor’s effluent absorbance at 260 nm as a function of time will provide a representative DNA release profile from the tissue and thus allow for process optimization. The DNA release profiles were obtained for laryngeal tissues and were successfully used to optimize the derivation of VF lamina propria-ECM (auVF-ECM) hydrogels. This hydrogel had comparable rheological properties to commonly used biomaterials to treat VF injuries. Also, the auVF-ECM hydrogel promoted the down-regulation of CCR7 by THP-1 macrophages upon lipopolysaccharide stimulation in vitro suggesting some anti-inflammatory properties. The results show that absorbance profiles are a good representation of DNA removal during the decellularization process thus providing an important tool to optimize future protocols.more » « less
-
Growing evidence suggests that intracranial pressure (ICP) plays an important role in the pathophysiology of glaucoma, especially in normal-tension glaucoma (NTG) patients. Controversial results exist about ICP’s relationship to visual field (VF) changes. With the aim to assess the relationship between ICP and VF zones in NTG patients, 80 NTG patients (age 59.5 (11.6) years) with early-stage glaucoma were included in this prospective study. Intraocular pressure (IOP) (Goldmann), visual perimetry (Humphrey) and non-invasive ICP (via a two-depth Transcranial Doppler, Vittamed UAB, Lithuania) were evaluated. Translaminar pressure difference (TPD) was calculated according to the formula TPD = IOP − ICP. The VFs of each patient were divided into five zones: nasal, temporal, peripheral, central, and paracentral. The average pattern deviation (PD) scores were calculated in each zone. The level of significance p < 0.05 was considered significant. NTG patients had a mean ICP of 8.5 (2.4) mmHg. Higher TPD was related with lower mean deviation (MD) (p = 0.01) and higher pattern standard deviation (PSD) (p = 0.01). ICP was significantly associated with the lowest averaged PD scores in the nasal VF zone (p < 0.001). There were no significant correlations between ICP and other VF zones with the most negative mean PD value. (p > 0.05). Further studies are needed to analyze the involvement of ICP in NTG management.more » « less
-
Abstract A primary cilium, made of nine microtubule doublets enclosed in a cilium membrane, is a mechanosensing organelle that bends under an external mechanical load and sends an intracellular signal through transmembrane proteins activated by cilium bending. The nine microtubule doublets are the main load‐bearing structural component, while the transmembrane proteins on the cilium membrane are the main sensing component. No distinction was made between these two components in all existing models, where the stress calculated from the structural component (nine microtubule doublets) was used to explain the sensing location, which may be totally misleading. For the first time, we developed a microstructure‐based primary cilium model by considering these two components separately. First, we refined the analytical solution of bending an orthotropic cylindrical shell for individual microtubule, and obtained excellent agreement between finite element simulations and the theoretical predictions of a microtubule bending as a validation of the structural component in the model. Second, by integrating the cilium membrane with nine microtubule doublets and simulating the tip‐anchored optical tweezer experiment on our computational model, we found that the microtubule doublets may twist significantly as the whole cilium bends. Third, besides being cilium‐length‐dependent, we found the mechanical properties of the cilium are also highly deformation‐dependent. More important, we found that the cilium membrane near the base is not under pure in‐plane tension or compression as previously thought, but has significant local bending stress. This challenges the traditional model of cilium mechanosensing, indicating that transmembrane proteins may be activated more by membrane curvature than membrane stretching. Finally, we incorporated imaging data of primary cilia into our microstructure‐based cilium model, and found that comparing to the ideal model with uniform microtubule length, the imaging‐informed model shows the nine microtubule doublets interact more evenly with the cilium membrane, and their contact locations can cause even higher bending curvature in the cilium membrane than near the base.more » « less
An official website of the United States government

