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
This content will become publicly available on December 19, 2024
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
- NSF-PAR ID:
- 10482086
- Publisher / Repository:
- AIP Publishing
- 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
-
more » « less
Abstract Modern approaches to modelling cardiac perfusion now commonly describe the myocardium using the framework of poroelasticity. Cardiac tissue can be described as a saturated porous medium composed of the pore fluid (blood) and the skeleton (myocytes and collagen scaffold). In previous studies fluid–structure interaction in the heart has been treated in a variety of ways, but in most cases, the myocardium is assumed to be a hyperelastic fibre‐reinforced material. Conversely, models that treat the myocardium as a poroelastic material typically neglect interactions between the myocardium and intracardiac blood flow. This work presents a poroelastic immersed finite element framework to model left ventricular dynamics in a three‐phase poroelastic system composed of the pore blood fluid, the skeleton, and the chamber fluid. We benchmark our approach by examining a pair of prototypical poroelastic formations using a simple cubic geometry considered in the prior work by Chapelle et al (2010). This cubic model also enables us to compare the differences between system behaviour when using isotropic and anisotropic material models for the skeleton. With this framework, we also simulate the poroelastic dynamics of a three‐dimensional left ventricle, in which the myocardium is described by the Holzapfel–Ogden law. Results obtained using the poroelastic model are compared to those of a corresponding hyperelastic model studied previously. We find that the poroelastic LV behaves differently from the hyperelastic LV model. For example, accounting for perfusion results in a smaller diastolic chamber volume, agreeing well with the well‐known wall‐stiffening effect under perfusion reported previously. Meanwhile differences in systolic function, such as fibre strain in the basal and middle ventricle, are found to be comparatively minor.
-
Determination of abdominal aortic aneurysm (AAA) rupture risk involves the accurate prediction of mechanical stresses acting on the arterial tissue, as well as the wall strength which has a correlation with oxygen supply within the aneurysmal wall. Our laboratory has previously reported the significance of an intraluminal thrombus (ILT) presence and morphology on localized oxygen deprivation by assuming a uniform consistency of ILT. The aim of this work is to investigate the effects of ILT structural composition on oxygen flow by adopting a multilayered porous framework and comparing a two-layer ILT model with one-layer models. Three-dimensional idealized and patient-specific AAA geometries are generated. Numerical simulations of coupled fluid flow and oxygen transport between blood, arterial wall, and ILT are performed, and spatial variations of oxygen concentrations within the AAA are obtained. A parametric study is conducted, and ILT permeability and oxygen diffusivity parameters are individually varied within a physiological range. A gradient of permeability is also defined to represent the heterogenous structure of ILT. Results for oxygen measures as well as filtration velocities are obtained, and it is found that the presence of any ILT reduces and redistributes the concentrations in the aortic wall markedly. Moreover, it is found that the integration of a porous ILT significantly affects the oxygen transport in AAA and the concentrations are linked to ILT’s permeability values. Regardless of the ILT stratification, maximum variation in wall oxygen concentrations is higher in models with lower permeability, while the concentrations are not sensitive to the value of the diffusion coefficient. Based on the observations, we infer that average one-layer parameters for ILT material characteristics can be used to reasonably estimate the wall oxygen concentrations in aneurysm models.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
