Interaction between a simplified soft palate and compressible viscous flow

Fluid-structure interaction in a simplified 2D model of the upper airways is simulated to study flow-induced oscillation of the soft palate in the pharynx. The goal of our research has been a better understanding of the mechanisms of the Obstructive Sleep Apnea Syndrome and snoring by taking into account compressible viscous flow. The inspiratory airflow is described by the 2D compressible Navier-Stokes equations, and the soft palate is modeled as a flexible plate by the linearized Euler-Bernoulli thin beam theory. Fluid-structure interaction is handled by the arbitrary Lagrangian-Eulerian formulation. The fluid flow is computed by utilizing 4th order accurate summation by parts difference operators and the 4th order accurate classical Runge-Kutta method which lead to very accurate simulation results. The motion of the cantilevered plate is solved numerically by employing the Newmark time integration method. The numerical schemes for the structure are verified by comparing the computed frequencies of plate oscillation with the associated second mode eigenfrequency in vacuum. Vortex dynamics is assessed for the coupled fluid-structure system when both airways are open and when one airway is closed. The effect of mass ratio, rigidity and damping coefficient of the plate on the oscillatory behavior is investigated. An acoustic analysis is carried out to characterize the acoustic wave propagation induced by the plate oscillation. It is observed that the acoustic wave corresponding to the quarter wave mode along the length of the duct is the dominant frequency. However, the frequency of the plate oscillation is recognizable in the acoustic pressure when reducing the amplitude of the quarter wave mode.