Computational fluid dynamics modeling of the upper airway of children with obstructive sleep apnea syndrome in steady flow

Numerical investigation on the flow characteristics and aerodynamic force of the upper airway of patient with obstructive sleep apnea using computational fluid dynamics

Developing a mathematical model to predict the abnormal flow characteristics that are produced by obstructive sleep apnea is an important step in learning the pathophysiology of the obstructive sleep apnea (OSA) disease. The present study provides detailed calculations of flow in the pharyngeal airway of a patient with obstructive sleep apnea. To achieve this goal, a computational fluid dynamics model was constructed using raw data from three-dimensional computed tomogram (CT) images of an OSA patient. To reproduce the important transition from laminar to turbulent flow in the pharyngeal airway, the low Reynolds number k-epsilon model was adopted and successfully validated using previous open literature. The results show that the flow in the pharyngeal airway of patients with OSA comprises a turbulent jet formed by area restriction at the velopharynx. This turbulent jet causes higher shear and pressure forces in the vicinity of the velopharynx. From the results, It may be deduced that the most collapsible area in the pharyngeal airway of OSA patients is the velopharynx where minimum intraluminal pressure and maximum aerodynamic force lie.

Development of a CFD urethral model to study flow-generated vortices under different conditions of prostatic obstruction

A novel, non-invasive method to diagnose bladder outlet obstruction involves the recording of noise with a contact microphone pressed against the perineum (between anus and scrotum). This noise results from flow-generated vortices caused by prostatic obstruction. We developed a computational fluid dynamic (CFD) urethral model including urethral geometry to study the relation between generated noise and the degree of obstruction. This model comprised a bladder, bladder neck, prostate and urethra. Calculations were carried out at four bladder pressures, five degrees of obstruction and three obstruction shapes. For each of the sixty simulations, the velocity and pressure distributions along the urethra were calculated including wall shear stresses to localize flow transition from disturbed to normal. Negative pressures at the obstruction outlet induced recirculation of flow. The location of transition was independent of the applied bladder pressure, but it depended primarily on the degree and secondarily on the shape of the obstruction. Based on the presented results, we hypothesize that the location of the maximum amplitude of perineal noise mainly depends on the degree and shape of the prostatic obstruction. Our future aim is to test our hypothesis in male patients and to extend the presented model to 3D with a viscoelastic urethral wall to calculate the fluid-wall interaction.