Effect of different degrees of adenoid hypertrophy on pediatric upper airway aerodynamics: a computational fluid dynamics study

Biomechanics and Modeling in Mechanobiology - To improve the diagnostic accuracy of adenoid hypertrophy (AH) in children and prevent further complications in time, it is important to study and...     

Model reduction methodology for computational simulations of endovascular repair

Endovascular aneurysm repair (EVAR) is a current alternative treatment for thoracic and abdominal aortic aneurysms, but is still sometimes compromised by possible complications such as device migration or endoleaks. In order to assist clinicians in preventing these complications, finite element analysis (FEA) is a promising tool. However, the strong material and geometrical nonlinearities added to the complex multiple contacts result in costly finite-element models. To reduce this computational cost, we establish here an alternative and systematic methodology to simplify the computational simulations of stent-grafts (SG) based on FEA. The model reduction methodology relies on equivalent shell models with appropriate geometrical and mechanical parameters. It simplifies significantly the contact interactions but still shows very good agreement with a complete reference finite-element model. Finally, the computational time for EVAR simulations is reduced of a factor 6–10. An application is shown for the deployment of a SG during thoracic endovascular repair, showing that the developed methodology is both effective and accurate to determine the final position of the deployed SG inside the aneurysm.     

Assessment of human left ventricle flow using statistical shape modelling and computational fluid dynamics

Blood flow patterns in the human left ventricle (LV) have shown relation to cardiac health. However, most studies in the literature are limited to a few patients and results are hard to generalize. This study aims to provide a new framework to generate more generalized insights into LV blood flow as a function of changes in anatomy and wall motion. In this framework, we studied the four-dimensional blood flow in LV via computational fluid dynamics (CFD) in conjunction with a statistical shape model (SSM), built from segmented LV shapes of 150 subjects. We validated results in an in-vitro dynamic phantom via time-resolved optical particle image velocimetry (PIV) measurements. This combination of CFD and the SSM may be useful for systematically assessing blood flow patterns in the LV as a function of varying anatomy and has the potential to provide valuable data for diagnosis of LV functionality.     

Airflow behavior changes in upper airway caused by different head and neck positions: Comparison by computational fluid dynamics

The feasibility of computational fluid dynamics (CFD) to evaluate airflow characteristics in different head and neck positions has not been established. This study compared the changes in volume and airflow behavior of the upper airway by CFD simulation to predict the influence of anatomical and physiological airway changes due to different head–neck positions on mechanical ventilation. One awake volunteer with no risk of difficult airway underwent computed tomography in neutral position, extension position (both head and neck extended), and sniffing position (head extended and neck flexed). Three-dimensional airway models of the upper airway were reconstructed. The total volume (V) and narrowest area (A_min) of the airway models were measured. CFD simulation with an Spalart–Allmaras model was performed to characterize airflow behavior in neutral, extension, and sniffing positions of closed-mouth and open-mouth ventilation. The comparison result for V was neutral <extension≈sniffing, and for A_min was neutral<extension<sniffing. A_min in sniffing position was nearly 3.0 times that in neutral position and 1.7 times that in extension position. The pressure drop and velocity increasing were more obvious in neutral than sniffing or extension position at the same airflow rate. In sniffing position, pressure differences decreased and velocity remained almost constant. Recirculation airflow was generated near the subglottic region in neutral and extension positions. Sniffing position improves airway patency by increasing airway volume and decreasing airway resistance, suggesting that sniffing position may be the optimal choice for mask ventilation.     

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

Computational fluid dynamic (CFD) analysis was used to model the effect of airway geometry on internal pressure in the upper airway of three children with obstructive sleep apnea syndrome (OSAS), and three controls. Model geometry was reconstructed from magnetic resonance images obtained during quiet tidal breathing, meshed with an unstructured grid, and solved at normative peak resting flow. The unsteady Reynolds-averaged Navier-Stokes equations were solved with steady flow boundary conditions in inspiration and expiration, using a two-equation low-Reynolds number turbulence model. Model results were validated using an in-vitro scale model, unsteady flow simulation, and reported nasal resistance measurements in children. Pharynx pressure drop strongly correlated to airway area restriction. Inspiratory pressure drop was primarily proportional to the square of flow, consistent with pressure losses due to convective acceleration caused by area restriction. On inspiration, in OSAS pressure drop occurred primarily between the choanae and the region where the adenoids overlap the tonsils (overlap region) due to airway narrowing, rather than in the nasal passages; in controls the majority of pressure drop was in the nasal passages. On expiration, in OSAS the majority of pressure drop occurred between the oropharynx (posterior to the tongue) and overlap region, and local minimum pressure in the overlap region was near atmospheric due to pressure recovery in the anterior nasopharynx. The results suggest that pharyngeal airway shape in children with OSAS significantly affects internal pressure distribution compared to nasal resistance. The model may also help explain regional dynamic airway narrowing during expiration.     

Hemodynamics in stented vertebral artery ostial stenosis based on computational fluid dynamics simulations

Hemodynamic factors may affect the potential occurrence of in-stent restenosis (ISR) after intervention procedure of vertebral artery ostial stenosis (VAOS). The purpose of the present study is to investigate the influence of stent protrusion length in implantation strategy on the local hemodynamics of the VAOS. CTA images of a 58-year-old female patient with posterior circulation transient ischemic attack were used to perform a 3D reconstruction of the vertebral artery. Five models of the vertebral artery before and after the stent implantation were established. Model 1 was without stent implantation, Model 2–5 was with stent protruding into the subclavian artery for 0, 1, 2, 3 mm, respectively. Computational fluid dynamics simulations based on finite element analysis were employed to mimic the blood flow in arteries and to assess hemodynamic conditions, particularly the blood flow velocity and wall shear stress (WSS). The WSS and the blood flow velocity at the vertebral artery ostium were reduced by 85.33 and 35.36% respectively after stents implantation. The phenomenon of helical flow disappeared. Hemodynamics comparison showed that stent struts that protruded 1 mm into the subclavian artery induced the least decrease in blood speed and WSS. The results suggest that stent implantation can improve the hemodynamics of VAOS, while stent struts that had protruded 1 mm into the subclavian artery would result in less thrombogenesis and neointimal hyperplasia and most likely decrease the risk of ISR.     

Effect of cut-off distance used in molecular dynamics simulations on fluid properties

The effect of cut-off distance used in molecular dynamics (MD) simulations on fluid properties was studied systematically in both canonical (NVT) and isothermal–isobaric (NPT) ensembles. Results show that the cut-off distance in the NVT ensemble plays little role in determining the equilibrium structure of fluid if the ensemble has a high density. However, pressures calculated in the same NVT ensembles strongly depend on the cut-off distance used. In the NPT ensemble, cut-off distance plays a key role in determining fluid equilibrium structure, density and self-diffusion coefficient. The characteristic of the radial distribution function of fluid in NPT ensembles depending on the cut-off distance used in MD simulations means that the WCA theory (a perturbation theory developed by Weeks, Chandler and Andersen) is not suitable for NPT ensembles because the assumption (the effect of the attractive force in determining the liquid structure is negligible) used in the WCA theory is not valid. The dependence of fluid properties on the cut-off distance also indicates that using the WCA potential (the repulsive part of the intermolecular potential proposed in the WCA theory) to calculate fluid transport in heterogeneous systems could lead to significant errors or incorrect results.     

Two-dimensional intraventricular flow pattern visualization using the image-based computational fluid dynamics

The image-based computational fluid dynamics (IB-CFD) technique, as the combination of medical images and the CFD method, is utilized in this research to analyze the left ventricle (LV) hemodynamics. The research primarily aims to propose a semi-automated technique utilizing some freely available and commercial software packages in order to simulate the LV hemodynamics using the IB-CFD technique. In this research, moreover, two different physiological time-resolved 2D models of a patient-specific LV with two different types of aortic and mitral valves, including the orifice-type valves and integrated with rigid leaflets, are adopted to visualize the process of developing intraventricular vortex formation and propagation. The blood flow pattern over the whole cardiac cycle of two models is also compared to investigate the effect of utilizing different valve types in the process of the intraventricular vortex formation. Numerical findings indicate that the model with integrated valves can predict more complex intraventricular flow that can match better the physiological flow pattern in comparison to the orifice-type model.     

Effect of inlet velocity profiles on patient-specific computational fluid dynamics simulations of the carotid bifurcation

Patient-specific computational fluid dynamics (CFD) is a powerful tool for researching the role of blood flow in disease processes. Modern clinical imaging technology such as MRI and CT can provide high resolution information about vessel geometry, but in many situations, patient-specific inlet velocity information is not available. In these situations, a simplified velocity profile must be selected. We studied how idealized inlet velocity profiles (blunt, parabolic, and Womersley flow) affect patient-specific CFD results when compared to simulations employing a "reference standard" of the patient's own measured velocity profile in the carotid bifurcation. To place the magnitude of these effects in context, we also investigated the effect of geometry and the use of subject-specific flow waveform on the CFD results. We quantified these differences by examining the pointwise percent error of the mean wall shear stress (WSS) and the oscillatory shear index (OSI) and by computing the intra-class correlation coefficient (ICC) between axial profiles of the mean WSS and OSI in the internal carotid artery bulb. The parabolic inlet velocity profile produced the most similar mean WSS and OSI to simulations employing the real patient-specific inlet velocity profile. However, anatomic variation in vessel geometry and the use of a nonpatient-specific flow waveform both affected the WSS and OSI results more than did the choice of inlet velocity profile. Although careful selection of boundary conditions is essential for all CFD analysis, accurate patient-specific geometry reconstruction and measurement of vessel flow rate waveform are more important than the choice of velocity profile. A parabolic velocity profile provided results most similar to the patient-specific velocity profile.     

Investigation of inhalation and exhalation flow pattern in a realistic human upper airway model by PIV experiments and CFD simulations

Biomechanics and Modeling in Mechanobiology - In this study, flow field characteristics in the trachea region in a realistic human upper airway model were firstly measured by particle image...     

Effects of cricothyroid muscle contraction on upper airway flow dynamics in dogs.

We studied the effects of cricothyroid muscle (CT) contraction on upper airway flow dynamics in eight prone open-mouth anesthetized dogs. Animals were mechanically ventilated via a tracheostomy while a constant airflow (Vuaw) passed through the isolated upper airway. Nasal airflow (Vn) was monitored using a nasal mask and pneumotachograph. Bilateral CT contraction was induced by electrical stimulation of the external branches of the superior laryngeal nerves. During CT contraction with Vuaw of 100-443 ml/s in the inspiratory direction, total upper airway resistance (Ruaw) fell by 49.1 +/- 5.4% (SE) while supraglottic resistance fell by 63.6 +/- 3.6%; simultaneously Vn fell by 55.3 +/- 3.8% and Vuaw increased by 7.2 +/- 1.7%. Similar results were obtained when Vuaw was in the expiratory direction. In three dogs in which the attachments of the CT to either the thyroid or cricoid cartilage were severed, superior laryngeal nerve stimulation had no systematic effect on Ruaw. Because visual assessment during CT contraction consistently revealed dilation of the piriform recesses, we suggest that CT contraction is associated with pharyngeal dilation, which in open-mouth dogs (with overlapping soft palate and epiglottis) redistributes flow to the oral route with a net reduction in Ruaw. Thus the CT may have a respiratory role as a pharyngeal dilator.     

Pulsatile flow computational simulations of mitral regurgitation.

The noninvasive quantification of mitral regurgitation remains an important clinical goal. Recently, the flow convergence method was developed to estimate the regurgitant flow rate. This study used three-dimensional pulsatile flow computational simulations to evaluate the accuracy of the flow convergence method in the presence of complicating factors such as ventricular confinement, noncircular orifice shape, and the presence of aortic outflow. Results showed that in the absence of aortic outflow and ventricular confinement, there was a plateau zone where the calculated flow rate by the hemispheric formula approximated the true flow rate, independent of the orifice shape. In the presence of aortic outflow and in chambers of physiologic dimensions, there was no longer a clear zone where the hemispheric formula was valid. The hemi-elliptic modification of the flow convergence method worked in all cases, independent of the degree of ventricular confinement or the presence of aortic outflow. Therefore, application of the hemi-elliptic formula should be considered in future clinical studies.     

MRI-based quantification of outflow boundary conditions for computational fluid dynamics of stenosed human carotid arteries

Accurate assessment of wall shear stress (WSS) is vital for studies on the pathogenesis of atherosclerosis. WSS distributions can be obtained by computational fluid dynamics (CFD) using patient-specific geometries and flow measurements. If patient-specific flow measurements are unavailable, in- and outflow have to be estimated, for instance by using Murray’s Law. It is currently unknown to what extent this law holds for carotid bifurcations, especially in cases where stenoses are involved. We performed flow measurements in the carotid bifurcation using phase-contrast MRI in patients with varying degrees of stenosis. An empirical relation between outflow and degree of area stenosis was determined and the outflow measurements were compared to estimations based on Murray’s Law. Furthermore, the influence of outflow conditions on the WSS distribution was studied.For bifurcations with an area stenosis smaller than 65%, the outflow ratio of the internal carotid artery (ICA) to the common carotid artery (CCA) was 0.62±0.12 while the outflow ratio of the external carotid artery (ECA) was 0.35±0.13. If the area stenosis was larger than 65%, the flow to the ICA decreased linearly to zero at 100% area stenosis. The empirical relation fitted the flow data well (R²=0.69), whereas Murray’s Law overestimated the flow to the ICA substantially for larger stenosis, resulting in an overestimation of the WSS. If patient-specific flow measurements of the carotid bifurcation are unavailable, estimation of the outflow ratio by the presented empirical relation will result in a good approximation of calculated WSS using CFD.     

On the choice of outlet boundary conditions for patient-specific analysis of aortic flow using computational fluid dynamics

Boundary conditions (BCs) are an essential part in computational fluid dynamics (CFD) simulations of blood flow in large arteries. Although several studies have investigated the influence of BCs on predicted flow patterns and hemodynamic wall parameters in various arterial models, there is a lack of comprehensive assessment of outlet BCs for patient-specific analysis of aortic flow. In this study, five different sets of outlet BCs were tested and compared using a subject-specific model of a normal aorta. Phase-contrast magnetic resonance imaging (PC-MRI) was performed on the same subject and velocity profiles extracted from the in vivo measurements were used as the inlet boundary condition. Computational results obtained with different outlet BCs were assessed in terms of their agreement with the PC-MRI velocity data and key hemodynamic parameters, such as pressure and flow waveforms and wall shear stress related indices. Our results showed that the best overall performance was achieved by using a well-tuned three-element Windkessel model at all model outlets, which not only gave a good agreement with in vivo flow data, but also produced physiological pressure waveforms and values. On the other hand, opening outlet BCs with zero pressure at multiple outlets failed to reproduce any physiologically relevant flow and pressure features.     

Sawfishes stealth revealed using computational fluid dynamics

Detailed computational fluid dynamics simulations for the rostrum of three species of sawfish (Pristidae) revealed that negligible turbulent flow is generated from all rostra during lateral swipe prey manipulation and swimming. These results suggest that sawfishes are effective stealth hunters that may not be detected by their teleost prey's lateral line sensory system during pursuits. Moreover, during lateral swipes, the rostra were found to induce little velocity into the surrounding fluid. Consistent with previous data of sawfish feeding behaviour, these data indicate that the rostrum is therefore unlikely to be used to stir up the bottom to uncover benthic prey. Whilst swimming with the rostrum inclined at a small angle to the horizontal, the coefficient of drag of the rostrum is relatively low and the coefficient of lift is zero.     

Developmental changes in upper airway dynamics.

Normal children have a less collapsible upper airway in response to subatmospheric pressure administration (P(NEG)) during sleep than normal adults do, and this upper airway response appears to be modulated by the central ventilatory drive. Children have a greater ventilatory drive than adults. We, therefore, hypothesized that children have increased neuromotor activation of their pharyngeal airway during sleep compared with adults. As infants have few obstructive apneas during sleep, we hypothesized that infants would have an upper airway that was resistant to collapse. We, therefore, compared the upper airway pressure-flow (V) relationship during sleep between normal infants, prepubertal children, and adults. We evaluated the upper airway response to 1). intermittent, acute P(NEG) (infants, children, and adults), and 2). hypercapnia (children and adults). We found that adults had a more collapsible upper airway during sleep than either infants or children. The children exhibited a vigorous response to both P(NEG) and hypercapnia during sleep (P < 0.01), whereas adults had no significant change. Infants had an airway that was resistant to collapse and showed a very rapid response to P(NEG). We conclude that the upper airway is resistant to collapse during sleep in infants and children. Normal children have preservation of upper airway responses to P(NEG) and hypercapnia during sleep, whereas responses are diminished in adults. Infants appear to have a different pattern of upper airway activation than older children. We speculate that the pharyngeal airway responses present in normal children are a compensatory response for a relatively narrow upper airway.     

Flow in the well: computational fluid dynamics is essential in flow chamber construction

A perfusion system was developed to generate well defined flow conditions within a well of a standard multidish. Human vein endothelial cells were cultured under flow conditions and cell response was analyzed by microscopy. Endothelial cells became elongated and spindle shaped. As demonstrated by computational fluid dynamics (CFD), cells were cultured under well defined but time varying shear stress conditions. A damper system was introduced which reduced pulsatile flow when using volumetric pumps. The flow and the wall shear stress distribution were analyzed by CFD for the steady and unsteady flow field. Usage of the volumetric pump caused variations of the wall shear stresses despite the controlled fluid environment and introduction of a damper system. Therefore the use of CFD analysis and experimental validation is critical in developing flow chambers and studying cell response to shear stress. The system presented gives an effortless flow chamber setup within a 6-well standard multidish.     

Analysis of a swimmer's hand and arm in steady flow conditions using computational fluid dynamics

Propulsive forces generated by swimmers' hands and arms have, to date, been determined strictly through experimental testing. As an alternative to these complex and costly experiments, the present research has applied the numerical technique of computational fluid dynamics (CFD) to calculate the steady flow around a swimmer's hand and arm at various angles of attack. Force coefficients computed for the hand and arm compared well with steady-state coefficients determined experimentally. The simulations showed significant boundary layer separation from the arm and hand, suggesting that Bernoulli's equation should not be used to mathematically describe the lift generated by a swimmer. Additionally, "2D" lift was shown to be inaccurate for the arm at all angles of attack and for the hand near angles of attack of 90 degrees. Such simulations serve to validate the chosen CFD techniques, and are an important first step towards the use of CFD methods for determining swimming hydrodynamic forces in more complex unsteady flow conditions.     

Design and analysis of flow velocity distribution inside a raceway pond using computational fluid dynamics

Bioprocess and Biosystems Engineering - Open raceway ponds are widely adopted for cultivating microalgae on a large scale. Working depth of the raceway pond is the major component to be analysed...