Computational fluid dynamics analysis of cyclist aerodynamics: Performance of different turbulence-modelling and boundary-layer modelling approaches

This study aims at assessing the accuracy of computational fluid dynamics (CFD) for applications in sports aerodynamics, for example for drag predictions of swimmers, cyclists or skiers, by evaluating the applied numerical modelling techniques by means of detailed validation experiments. In this study, a wind-tunnel experiment on a scale model of a cyclist (scale 1:2) is presented. Apart from three-component forces and moments, also high-resolution surface pressure measurements on the scale model’s surface, i.e. at 115 locations, are performed to provide detailed information on the flow field. These data are used to compare the performance of different turbulence-modelling techniques, such as steady Reynolds-averaged Navier–Stokes (RANS), with several k–ε and k–ω turbulence models, and unsteady large-eddy simulation (LES), and also boundary-layer modelling techniques, namely wall functions and low-Reynolds number modelling (LRNM). The commercial CFD code Fluent 6.3 is used for the simulations. The RANS shear-stress transport (SST) k–ω model shows the best overall performance, followed by the more computationally expensive LES. Furthermore, LRNM is clearly preferred over wall functions to model the boundary layer. This study showed that there are more accurate alternatives for evaluating flow around bluff bodies with CFD than the standard k–ε model combined with wall functions, which is often used in CFD studies in sports.     

Turbulence model choice for the calculation of drag forces when using the CFD method

The aim of this work is to specify which model of turbulence is the most adapted in order to predict the drag forces that a swimmer encounters during his movement in the fluid environment. For this, a Computational Fluid Dynamics (CFD) analysis has been undertaken with a commercial CFD code (Fluent^®). The problem was modelled as 3D and in steady hydrodynamic state. The 3D geometry of the swimmer was created by means of a complete laser scanning of the swimmer’s body contour. Two turbulence models were tested, namely the standard k–ε model with a specific treatment of the fluid flow area near the swimmer’s body contour, and the standard k–ω model. The comparison of numerical results with experimental measurements of drag forces shows that the standard k–ω model accurately predicts the drag forces while the standard k–ε model underestimates their values. The standard k–ω model also enabled to capture the vortex structures developing at the swimmer’s back and buttocks in underwater swimming; the same vortices had been visualized by flow visualization experiments carried out at the INSEP (National Institute for Sport and Physical Education in Paris) with the French national swimming team.     

Effect of the velopharynx on intraluminal pressures in reconstructed pharynges derived from individuals with and without sleep apnea

The most collapsible part of the upper airway in the majority of individuals is the velopharynx which is the segment positioned behind the soft palate. As such it is an important morphological region for consideration in elucidating the pathogenesis of obstructive sleep apnea (OSA). This study compared steady flow properties during inspiration in the pharynges of nine male subjects with OSA and nine body-mass index (BMI)- and age-matched control male subjects without OSA. The k  –ω$\omega$SST turbulence model was used to simulate the flow field in subject-specific pharyngeal geometric models reconstructed from anatomical optical coherence tomography (aOCT) data. While analysis of the geometry of reconstructed pharynges revealed narrowing at velopharyngeal level in subjects with OSA, it was not possible to clearly distinguish them from subjects without OSA on the basis of pharyngeal size and shape alone. By contrast, flow simulations demonstrated that pressure fields within the narrowed airway segments were sensitive to small differences in geometry and could lead to significantly different intraluminal pressure characteristics between subjects. The ratio between velopharyngeal and total pharyngeal pressure drops emerged as a relevant flow-based criterion by which subjects with OSA could be differentiated from those without.     

Computational fluid dynamics simulation improves the design and characterization of a plug-flow-type scale-down reactor for microbial cultivation processes

The scale-up of bioprocesses remains one of the major obstacles in the biotechnology industry. Scale-down bioreactors have been identified as valuable tools to investigate the heterogeneities observed in large-scale tanks at the laboratory scale. Additionally, computational fluid dynamics (CFD) simulations can be used to gain information about fluid flow in tanks used for production. Here, we present the rational design and comprehensive characterization of a scale-down setup, in which a flexible and modular plug-flow reactor was connected to a stirred-tank bioreactor. With the help of CFD using the realizable k-ε model, the mixing time difference between a 20 and 4000 L bioreactor was evaluated and used as scale-down criterion. CFD simulations using a shear stress transport (SST) k-ω turbulence model were used to characterize the plug-flow reactor in more detail, and the model was verified using experiments. Additionally, the model was used to simulate conditions where experiments technically could not be performed due to sensor limitations. Nevertheless, verification is difficult in this case as well. This was the first time a scale-down setup was tested on high-cell-density Escherichia coli cultivations to produce industrially relevant antigen-binding fragments (Fab). Biomass yield was reduced by 11% and specific product yield was reduced by 20% during the scale-down cultivations. Additionally, the intracellular Fab fraction was increased by using the setup. The flexibility of the introduced scale-down setup in combination with CFD simulations makes it a valuable tool for investigating scale effects at the laboratory scale. More information about the large scale is still necessary to further refine the setup and to speed up bioprocess scale-up in the future.     

Large Eddy Simulation and Reynolds-Averaged Navier-Stokes modeling of flow in a realistic pharyngeal airway model: an investigation of obstructive sleep apnea

Computational fluid dynamics techniques employing primarily steady Reynolds-Averaged Navier-Stokes (RANS) methodology have been recently used to characterize the transitional/turbulent flow field in human airways. The use of RANS implies that flow phenomena are averaged over time, the flow dynamics not being captured. Further, RANS uses two-equation turbulence models that are not adequate for predicting anisotropic flows, flows with high streamline curvature, or flows where separation occurs. A more accurate approach for such flow situations that occur in the human airway is Large Eddy Simulation (LES). The paper considers flow modeling in a pharyngeal airway model reconstructed from cross-sectional magnetic resonance scans of a patient with obstructive sleep apnea. The airway model is characterized by a maximum narrowing at the site of retropalatal pharynx. Two flow-modeling strategies are employed: steady RANS and the LES approach. In the RANS modeling framework both k-epsilon and k-omega turbulence models are used. The paper discusses the differences between the airflow characteristics obtained from the RANS and LES calculations. The largest discrepancies were found in the axial velocity distributions downstream of the minimum cross-sectional area. This region is characterized by flow separation and large radial velocity gradients across the developed shear layers. The largest difference in static pressure distributions on the airway walls was found between the LES and the k-epsilon data at the site of maximum narrowing in the retropalatal pharynx.     

Computational fluid dynamics endpoints for assessment of adenotonsillectomy outcome in obese children with obstructive sleep apnea syndrome

Background

Improvements in obstructive sleep apnea syndrome (OSAS) severity may be associated with improved pharyngeal fluid mechanics following adenotonsillectomy (AT). The study objective is to use image-based computational fluid dynamics (CFD) to model changes in pharyngeal pressures after AT, in obese children with OSAS and adenotonsillar hypertrophy.

Methods

Three-dimensional models of the upper airway from nares to trachea, before and after AT, were derived from magnetic resonance images obtained during wakefulness, in a cohort of 10 obese children with OSAS. Velocity, pressure, and turbulence fields during peak tidal inspiratory flow were computed using commercial software. CFD endpoints were correlated with polysomnography endpoints before and after AT using Spearman?s rank correlation (r_s).

Results

Apnea hypopnea index (AHI) decreases after AT was strongly correlated with reduction in maximum pressure drop (dP_TAmax) in the region where tonsils and adenoid constrict the pharynx (r_s=0.78, P=0.011), and with decrease of the ratio of dP_TAmax to flow rate (r_s=0.82, P=0.006). Correlations of AHI decrease to anatomy, negative pressure in the overlap region (including nasal flow resistance), or pressure drop through the entire pharynx, were not significant. In a subgroup of subjects with more than 10% improvement in AHI, correlations between flow variables and AHI decrease were stronger than in all subjects.

Conclusions

The correlation between change in dP_TAmax and improved AHI suggests that dP_TAmax may be a useful index for internal airway loading due to anatomical narrowing, and may be better correlated with AHI than direct airway anatomic measurements.     

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.     

Influence of pool geometry on the biological efficiency of vertical slot fishways

The characteristics of the flow in vertical slot fishways depend mainly on the specific pool design, i.e., the geometry of the pool. This paper presents the results of a study on sixteen different designs of vertical slot fishways with two different slopes. The flow field is studied using both a physical-scale model and a numerical model which solves the 2D depth-averaged shallow water equations coupled with a k–? turbulence model. Several hydraulic parameters, such as the average velocity in the slot or the water depth, are considered in order to evaluate their biological efficiency. Based on the results, it appears that the length of the pool is the main geometric dimension that affects the flow in the fishway. The numerical and experimental agreement is very satisfactory.     

An integrated physical and biological model for anaerobic lagoons

A computational fluid dynamics (CFD) model that integrates physical and biological processes for anaerobic lagoons is presented. In the model development, turbulence is represented using a transition k-ω model, heat conduction and solar radiation are included in the thermal model, biological oxygen demand (BOD) reduction is characterized by first-order kinetics, and methane yield rate is expressed as a linear function of temperature. A test of the model applicability is conducted in a covered lagoon digester operated under tropical climate conditions. The commercial CFD software, ANSYS-Fluent, is employed to solve the integrated model. The simulation procedures include solving fluid flow and heat transfer, predicting local resident time based on the converged flow fields, and calculating the BOD reduction and methane production. The simulated results show that monthly methane production varies insignificantly, but the time to achieve a 99% BOD reduction in January is much longer than that in July.     

Aerodynamic study of different cyclist positions: CFD analysis and full-scale wind-tunnel tests

Three different cyclist positions were evaluated with Computational Fluid Dynamics (CFD) and wind-tunnel experiments were used to provide reliable data to evaluate the accuracy of the CFD simulations. Specific features of this study are: (1) both steady Reynolds-averaged Navier–Stokes (RANS) and unsteady flow modelling, with more advanced turbulence modelling techniques (Large-Eddy Simulation – LES), were evaluated; (2) the boundary layer on the cyclist’s surface was resolved entirely with low-Reynolds number modelling, instead of modelling it with wall functions; (3) apart from drag measurements, also surface pressure measurements on the cyclist’s body were performed in the wind-tunnel experiment, which provided the basis for a more detailed evaluation of the predicted flow field by CFD. The results show that the simulated and measured drag areas differed about 11% (RANS) and 7% (LES), which is considered to be a close agreement in CFD studies. A fair agreement with wind-tunnel data was obtained for the predicted surface pressures, especially with LES. Despite the higher accuracy of LES, its much higher computational cost could make RANS more attractive for practical use in some situations. CFD is found to be a valuable tool to evaluate the drag of different cyclist positions and to investigate the influence of small adjustments in the cyclist’s position. A strong advantage of CFD is that detailed flow field information is obtained, which cannot easily be obtained from wind-tunnel tests. This detailed information allows more insight in the causes of the drag force and provides better guidance for position improvements.     

Can sites prone to flow induced vascular complications in a-v fistulas be assessed using computational fluid dynamics?

Arterio-venous fistulas (shunts between arteries and veins) are the preferred vascular access for hemodialysis. Despite their superior patency, compared with synthetic tubes and grafts, functional problems and inadequate flow rates are the common complications. Local flow conditions, in particular low and oscillating wall shear stresses (WSS), are central to vascular problems and a robust framework for analyzing flow conditions in vascular structures could provide an understanding of the mechanisms leading to vascular complications, such as stenoses, aneurisms, and thromboses.We hypothesize that a validated computational fluid dynamics (CFD) framework can be used to identify critical fistula configurations with elevated risk of complications. Therefore, the aim of the present study was to develop a CFD framework for analyzing fluid flow in complex vascular structures, such as arterio-venous fistulas validated by comparisons of in vitro volume flows with CFD results and flow fields from ultrasound scans with CFD simulations.Volume flows measured in vitro and CFD data differed quantitatively. However, good relative correlations exist between the data using logarithmic scales. Qualitatively, visual comparisons between ultrasound and CFD images showed good agreement between the two methods. In addition, WSS levels and the oscillatory shear index (OSI) were calculated and visualized on the model surface. The method was successfully validated and the method is deemed suitable for more thorough investigations into the field of vascular complications in a-v fistulas.     

Contribution of upper airway geometry to convective mixing

We investigated the axial dispersive effect of the upper airway structure (comprising mouth cavity, oropharynx, and trachea) on a traversing aerosol bolus. This was done by means of aerosol bolus experiments on a hollow cast of a realistic upper airway model (UAM) and three-dimensional computational fluid dynamics (CFD) simulations in the same UAM geometry. The experiments showed that 50-ml boluses injected into the UAM dispersed to boluses with a half-width ranging from 80 to 90 ml at the UAM exit, across both flow rates (250, 500 ml/s) and both flow directions (inspiration, expiration). These experimental results imply that the net half-width induced by the UAM typically was 69 ml. Comparison of experimental bolus traces with a one-dimensional Gaussian-derived analytical solution resulted in an axial dispersion coefficient of 200-250 cm(2)/s, depending on whether the bolus peak and its half-width or the bolus tail needed to be fully accounted for. CFD simulations agreed well with experimental results for inspiratory boluses and were compatible with an axial dispersion of 200 cm(2)/s. However, for expiratory boluses the CFD simulations showed a very tight bolus peak followed by an elongated tail, in sharp contrast to the expiratory bolus experiments. This indicates that CFD methods that are widely used to predict the fate of aerosols in the human upper airway, where flow is transitional, need to be critically assessed, possibly via aerosol bolus simulations. We conclude that, with all its geometric complexity, the upper airway introduces a relatively mild dispersion on a traversing aerosol bolus for normal breathing flow rates in inspiratory and expiratory flow directions.     

Equation for model intracranial aneurysm with consideration of small dissipation term

An electric circuit model of a ‘model intracranial aneurysm’ has been developed from previous experimental work showing a non linear pressure-volume relationship. This leads to a revised formulation of the following non linear differential equation which simulates the response of the model aneurysm to a pulsatile driving function

$\text{V}\text{?}\text{+?}\text{υ}\text{?}\text{+ζmυ+ξnυ}{}^2\text{=F cos(ωt + B)}$

, where V = velocity of flow within the aneurysm.Solution of this equation using the principle of harmonic balance shows that discontinuous jumps in amplitude of velocity of flow leading to turbulence or increased turbulence can occur with an increase in pulse rate or pulse pressure. It is postulated that turbulence is dangerous due to both the large increase in velocity of flow and also due to the possible stasis effect on large molecules and platelets from eddy formation.     

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.     

Numerical simulation of aneurysmal haemodynamics with calibrated porous-medium models of flow-diverting stents

Modelling flow-diverting (FD) stents as porous media (PM) markedly improves the efficiency of computational fluid dynamics (CFD) simulations in the study of intracranial aneurysm treatment. Nonetheless, the parameters of PM models adopted for simulations up until now were rarely calibrated to match the represented FD structure. We therefore sought to evaluate the PM parameters for a representative variety of commercially available stents, so characterising the flow-diversion behaviours of different FD devices on the market.We generated fully-resolved geometries for treatments using PED, Silk+, FRED, and dual PED stents. We then correspondingly derived the calibrated PM parameters—permeability (k) and inertial resistance factor (C₂)—for each stent design from CFD simulations, to ensure the calibrated PM model has identical flow resistance to the FD stent it represents. With each of the calibrated PM models respectively deployed in two aneurysms, we studied the flow-diversion effects of these stent configurations.This work for the first time reported several sets of parameters for PM models, which is vital to address the current knowledge gap and rectify the errors in PM model simulations, thereby setting right the modelling protocol for future studies using PM models. The flow resistance parameters were strongly affected by porosity and effective thickness of the commercial stents, and thus accounted for in the PM models. Flow simulations using the PM stent models revealed differences in aneurysmal mass flowrate (MFR) and energy loss (EL) between various stent designs.This study improves the practicability of FD simulation by using calibrated PM models, providing an individualised method with improved simulation efficiency and accuracy.     

Low Reynolds number turbulence modeling of blood flow in arterial stenoses.

Moderate and severe arterial stenoses can produce highly disturbed flow regions with transitional and or turbulent flow characteristics. Neither laminar flow modeling nor standard two-equation models such as the kappa-epsilon turbulence ones are suitable for this kind of blood flow. In order to analyze the transitional or turbulent flow distal to an arterial stenosis, authors of this study have used the Wilcox low-Re turbulence model. Flow simulations were carried out on stenoses with 50, 75 and 86% reductions in cross-sectional area over a range of physiologically relevant Reynolds numbers. The results obtained with this low-Re turbulence model were compared with experimental measurements and with the results obtained by the standard kappa-epsilon model in terms of velocity profile, vortex length, wall shear stress, wall static pressure, and turbulence intensity. The comparisons show that results predicted by the low-Re model are in good agreement with the experimental measurements. This model accurately predicts the critical Reynolds number at which blood flow becomes transitional or turbulent distal an arterial stenosis. Most interestingly, over the Re range of laminar flow, the vortex length calculated with the low-Re model also closely matches the vortex length predicted by laminar flow modeling. In conclusion, the study strongly suggests that the proposed model is suitable for blood flow studies in certain areas of the arterial tree where both laminar and transitional/turbulent flows coexist.     

Pulsatile spiral blood flow through arterial stenosis

Pulsatile spiral blood flow in a modelled three-dimensional arterial stenosis, with a 75% cross-sectional area reduction, is investigated by using numerical fluid dynamics. Two-equation k-ω model is used for the simulation of the transitional flow with Reynolds numbers 500 and 1000. It is found that the spiral component increases the static pressure in the vessel during the deceleration phase of the flow pulse. In addition, the spiral component reduces the turbulence intensity and wall shear stress found in the post-stenosis region of the vessel in the early stages of the flow pulse. Hence, the findings agree with the results of Stonebridge et al. (2004). In addition, the results of the effects of a spiral component on time-varying flow are presented and discussed along with the relevant pathological issues.     

Direct numerical simulation of transitional flow in a patient-specific intracranial aneurysm

In experiments turbulence has previously been shown to occur in intracranial aneurysms. The effects of turbulence induced oscillatory wall stresses could be of great importance in understanding aneurysm rupture. To investigate the effects of turbulence on blood flow in an intracranial aneurysm, we performed a high resolution computational fluid dynamics (CFD) simulation in a patient specific middle cerebral artery (MCA) aneurysm using a realistic, pulsatile inflow velocity. The flow showed transition to turbulence just after peak systole, before relaminarization occurred during diastole. The turbulent structures greatly affected both the frequency of change of wall shear stress (WSS) direction and WSS magnitude, which reached a maximum value of 41.5Pa. The recorded frequencies were predominantly in the range of 1-500Hz. The current study confirms, through properly resolved CFD simulations that turbulence can occur in intracranial aneurysms.     

Numerical modeling of pulsatile turbulent flow in stenotic vessels

Pulsatile turbulent flow in stenotic vessels has been numerically modeled using the Reynolds-averaged Navier-Stokes equation approach. The commercially available computational fluid dynamics code (CFD), FLUENT, has been used for these studies. Two different experiments were modeled involving pulsatile flow through axisymmetric stenoses. Four different turbulence models were employed to study their influence on the results. It was found that the low Reynolds number k-omega turbulence model was in much better agreement with previous experimental measurements than both the low and high Reynolds number versions of the RNG (renormalization-group theory) k-epsilon turbulence model and the standard k-epsilon model, with regard to predicting the mean flow distal to the stenosis including aspects of the vortex shedding process and the turbulent flow field. All models predicted a wall shear stress peak at the throat of the stenosis with minimum values observed distal to the stenosis where flow separation occurred.     

Nonlinear Compliance Modulates Dynamic Bronchoconstriction in a Multiscale Airway Model

The role of breathing and deep inspirations (DI) in modulating airway hyperresponsiveness remains poorly understood. In particular, DIs are potent bronchodilators of constricted airways in nonasthmatic subjects but not in asthmatic subjects. Additionally, length fluctuations (mimicking DIs) have been shown to reduce mean contractile force when applied to airway smooth muscle (ASM) cells and tissue strips. However, these observations are not recapitulated on application of transmural pressure (P_TM) oscillations (that mimic tidal breathing and DIs) in isolated intact airways. To shed light on this paradox, we have developed a biomechanical model of the intact airway, accounting for strain-stiffening due to collagen recruitment (a large component of the extracellular matrix (ECM)), and dynamic actomyosin-driven force generation by ASM cells. In agreement with intact airway studies, our model shows that P_TM fluctuations at particular mean transmural pressures can lead to only limited bronchodilation. However, our model predicts that moving the airway to a more compliant point on the static pressure-radius relationship (which may involve reducing mean P_TM), before applying pressure fluctuations, can generate greater bronchodilation. This difference arises from competition between passive strain-stiffening of ECM and force generation by ASM yielding a highly nonlinear relationship between effective airway stiffness and P_TM, which is modified by the presence of contractile agonist. Effectively, the airway at its most compliant may allow for greater strain to be transmitted to subcellular contractile machinery. The model predictions lead us to hypothesize that the maximum possible bronchodilation of an airway depends on its static compliance at the P_TM about which the fluctuations are applied. We suggest the design of additional experimental protocols to test this hypothesis.