An official website of the United States government
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
Biotransport in human phonation: Porous vocal fold tissue and fluid–structure interaction
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
- Award ID(s):
- 2138225
- PAR ID:
- 10589105
- Publisher / Repository:
- American Institute of Physics
- Date Published:
- Journal Name:
- Physics of Fluids
- Volume:
- 35
- Issue:
- 12
- ISSN:
- 1070-6631
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Abstract This paper proposes a deep-learning based generalized empirical flow model (EFM) that can provide a fast and accurate prediction of the glottal flow during normal phonation. The approach is based on the assumption that the vibration of the vocal folds can be represented by a universal kinematics equation (UKE), which is used to generate a glottal shape library. For each shape in the library, the ground truth values of the flow rate and pressure distribution are obtained from the high-fidelity Navier-Stokes (N-S) solution. A fully-connected deep neural network (DNN) is then trained to build the empirical mapping between the shapes and the flow rate and pressure distributions. The obtained DNN based EFM is coupled with a finite-element method (FEM) based solid dynamics solver for flow-structure-interaction (FSI) simulation of phonation. The EFM is evaluated by comparing to the N-S solutions in both static glottal shapes and FSI simulations. The results demonstrate a good prediction performance in accuracy and efficiency.more » « less
-
Microfluidic devices offer well-defined physical environments that are suitable for effective cell seeding and in vitro three-dimensional (3D) cell culture experiments. These platforms have been employed to model in vivo conditions for studying mechanical forces, cell–extracellular matrix (ECM) interactions, and to elucidate transport mechanisms in 3D tissue-like structures, such as tumor and lymph node organoids. Studies have shown that fluid flow behavior in microfluidic slides (µ-slides) directly influences shear stress, which has emerged as a key factor affecting cell proliferation and differentiation. This study investigates fluid flow in the porous channel of a µ-slide using computational fluid dynamics (CFD) techniques to analyze the impact of perfusion flow rate and porous properties on resulting shear stresses. The model of the µ-slide filled with a permeable biomaterial is considered. Porous media fluid flow in the channel is characterized by adding a momentum loss term to the standard Navier–Stokes equations, with a physiological range of permeability values. Numerical simulations are conducted to obtain data and contour plots of the filtration velocity and flow-induced shear stress distributions within the device channel. The filtration flow is subsequently measured by performing protein perfusions into the slide embedded with native human-derived ECM, while the flow rate is controlled using a syringe pump. The relationships between inlet flow rate and shear stress, as well as filtration flow and ECM permeability, are analyzed. The findings provide insights into the impact of shear stress, informing the optimization of perfusion conditions for studying tissues and cells under fluid flow.more » « less
-
Reactive transport modeling of subsurface environments plays an important role in addressing critical problems of geochemical processes, such as dissolution and precipitation of minerals. Current transport models for porous media span various scales, ranging from pore-scale to continuum-scale. In this study, we established an upscaling method connecting pore-scale and continuum-scale models by employing a deep learning methodology of Convolutional Neural Networks (CNNs). We applied Darcy-Brinkmann-Stokes (DBS) method to simulate the fluid flow and reactive transport in pore-scale models, which would act as constituents of a continuum-scale model. The datasets of spatial pore distribution of subvolume samples were used as the input for the deep learning model, and the continuum (Darcy)-scale parameters such as permeability, effective surface area, and effective diffusion coefficient were figured out as outputs (i.e., labels). By applying the trained models of the subvolumes in the entire sample volume, we generated the initial field of porosity, permeability, effective diffusion coefficient, and effective surface area for continuum-scale simulation of a mineral dissolution problem. We took an acid dissolution case as an example to utilize the outcomes of trained deep learning models as input data in the continuum-scale simulation. This work presents a comprehensive upscaling workflow, as bridging the findings of microscale simulations to the continuum-scale simulations of a reactive transport problem.more » « less
-
In this paper, we report a novel symmetry-breaking phenomenon that occurs in turbulent convective flow in periodic porous media in the intermediate porosity flow regime for values of porosities in between 0.8 and 0.9. Large eddy simulation is used to numerically simulate the momentum and thermal transport inside the porous medium at the microscale level. The phenomenon is observed to occur for periodically repeating porous media consisting of an in-line arrangement of circular cylinder solid obstacles, such as typically found in heat exchangers. The transition from symmetric to asymmetric flow occurs in between the pore scale Reynolds numbers of 37 (laminar) and 100 (turbulent), and asymmetric flow patterns are reported for Reynolds numbers up to 1,000. A Hopf bifurcation resulting in unsteady oscillatory laminar flow marks the origin of a secondary flow instability arising from the interaction of the shear layers around the solid obstacle. In turbulent flow, stochastic phase difference in the vortex wake oscillations caused by the secondary flow instability results in asymmetrical velocity and temperature distributions in the pore space. Consequently, high and low velocity flow channels are formed in the pore space that leads to the asymmetrical velocity and pressure distributions. At the macroscale level, symmetry-breaking results in residual transverse drag force components acting on the solid obstacle surfaces. The vortex wake oscillations caused by the secondary flow instability promote attached flow on the solid obstacle surface, which improves the surface averaged heat flux from the solid obstacle surface.more » « less
