Comparative predictions of turbulent non-isothermal flow of a viscoplastic fluid with yield stress

RANS simulation of turbulent non-isothermal flow in a pipe by transition of Newtonian fluid into a viscoplastic non-Newtonian fluid is carried out. Carrier phase turbulence was modeled using two-equation isotropic k‒ε ˜ - model and Reynolds stress transport model. Results of calculations of Newtonian and non-Newtonian power-law fluids were compared with data of DNS calculations of other authors. Comparisons with data from other works for isothermal Schwedoff-Bingham fluid were performed in this work using k‒ε ˜ - model. Satisfactory agreement was obtained with data from other studies for the axial averaged velocity and turbulent kinetic energy radial profiles (the difference is up to 15 %). Reynolds stress transport model showed significant anisotropy between streamwise and transverse velocity fluctuations (up to several times) and good agreement with DNS results of other authors. Averaged and pulsation profiles express the indicated transformation of the non-isothermal turbulent flow.

Computational investigation of NOx emission in biogas combustion

Landfills and anaerobic digesters in the waste treatment processes generate biogas. Biogas can be used as a fuel and excess biogas is typically burned in a flare to reduce the greenhouse effect. However, burning biogas produces several pollutants, including CO2, NOx, and SO2. To minimize these emissions, the amount of excess air used in the combustion process needs to be considered, which has a significant impact on NOx emissions. This study developed a Computational Fluid Dynamics (CFD) model to simulate a small-scale biogas combustion system and analyses the effect of excess air on heat output and NOx emissions during biogas combustion. The GRI-Mech reaction mechanism was used to simulate reactions, and the model was validated by comparing it to experimental data from the DLR-Stuttgart CH4/H2/N2 Jet Flame. To reduce computational costs, a Tabulation of Dynamic Adaptive Chemistry (TDAC) algorithm was used to dynamically adapt the reaction mechanism in real time. Turbulence in the DLR flame was simulated using Reynolds-Averaged Navier-Stokes (RANS). The CFD model used a co-flow of a natural draft to provide additional air, while the air was premixed with fuel. The CFD model was used to simulate various premixed equivalent ratios, and the resulting emissions and heat outputs were compared. The study found that the optimal premixed equivalent ratio for the studied system was between 0.85 and 1.1, as this range produced the highest temperature and lowest NOx emissions. This model facilitates emission analysis of gas-phase combustion systems.

Determination of optimal height of barriers to reduce the amount of pollution in the viaduct settings in an idealized urban canyon: a numerical study

In this work, we numerically investigate the process of atmospheric air pollution in idealized urban canyons along the road in the presence of a viaduct, taking into account different height of barriers. To solve this problem, the 3D Reynolds-averaged Navier-Stokes equations (RANS) were used. The closure of this system of equations was achieved by using various turbulent models. The verification of the mathematical model and the numerical algorithm was carried out using a test problem. The obtained results using various turbulent models were compared with experimental data and calculated results of other authors. The main problem considered in this work is characterized as follows: assessment of emissions of pollutants between buildings using barriers of various types in the presence of a viaduct. Computational results have shown that the barrier viaduct plays a large role in improving air quality in urban canyons. So, for example, a barrier erected on a viaduct with a height of 2 m reduces the concentration value to a cross-section x = 84 by more than 2 times in comparison with the case of a complete absence of protective barriers. A similar situation was observed with barriers erected above the earth's surface: located along the road, they also significantly reduce the value of the concentration of pollutants. Thus, the presence of barriers in both cases is necessary to prevent the dispersion and deposition of pollutants.

CFD simulation of gas pressure drop in porous packing for rotating packed beds (RPB) CO2 absorbers

Rotating packed bed (RPB) is a promising technology which can be used to intensify mass transfer in absorption processes. A better understanding of fluid dynamics is crucial to fill the gap in fundamental knowledge. Raising awareness on new technology and creating rules for process design and control are also very important. The experimental investigation of fluid in rotating beds is a very complex and difficult issue. What is more, the knowledge of the phase behavior in an RPB device is still insufficient. Therefore, an CFD (computational fluid dynamics) simulation is proposed as a tool for the study of gas phase flow inside porous packing. This study presents a three-dimensional numerical model for two fluid models: k-ε and RNG k-ε, for predicting dry pressure drop. The obtained simulation outcome was compared with the experimental results. The experimental dry pressure drop for porous packing was investigated for rotational speed in the range from 150 rpm to 1500 rpm and compared to the results from the CFD model. The comparison between the experimental and simulation results indicates very good consistency for the entire range of the rotational speed of interest. CFD modelling is recognised as an adequate tool leading to the better understanding of gas phase behaviour inside an RPB, filling an essential gap in our knowledge of the hydrodynamics of rotating packing, which allows to improve the design and performance of the process in RPB in terms of minimizing energy and material consumption.

The effect of turbulence modelling on the assessment of platelet activation

Pathological platelet activation by abnormal shear stresses is regarded as a main clinical complication in recipients of cardiovascular mechanical devices. In order to improve their performance computational fluid dynamics (CFD) are used to evaluate flow fields and related shear stresses. CFD models are coupled with mathematical models that describe the relation between fluid dynamics variables, and in particular shear stresses, and the platelet activation state (PAS). These models typically use a Lagrangian approach to compute the shear stresses along possible platelet trajectories. However, in the case of turbulent flow, the choice of the proper turbulence closure is still debated for both concerning its effect on shear stress calculation and Lagrangian statistics. In this study different numerical simulations of the flow through a mechanical heart valve were performed and then compared in terms of Eulerian and Lagrangian quantities: a direct numerical simulation (DNS), a large eddy simulation (LES), two Reynolds-averaged Navier-Stokes (RANS) simulations (SST k-ω and RSM) and a "laminar" (no turbulence modelling) simulation. Results exhibit a large variability in the PAS assessment depending on the turbulence model adopted. "Laminar" and RSM estimates of platelet activation are about 60% below DNS, while LES is 16% less. Surprisingly, PAS estimated from the SST k- ω velocity field is only 8% less than from DNS data. This appears more artificial than physical as can be inferred after comparing frequency distributions of PAS and of the different Lagrangian variables of the mechano-biological model of platelet activation. Our study indicates how much turbulence closures may affect platelet activation estimates, in comparison to an accurate DNS, when assessing blood damage in blood contacting devices.

VITOM-3D-assisted retroauricular neck surgery (RANS-3D): preliminary experience at Candiolo Cancer Institute

Objective

The recent introduction of 3D exoscopic surgery has engendered interesting technical improvements in head and neck surgery. The main goal of this study was to describe the application of 3D exoscopic technology on a wide range of pathologies of the neck, benign and malignant, through a minimally invasive retroauricular approach.

Methods

In the period January-December, 2019, 40 consecutive patients underwent neck surgery with a retroauricular approach, enhanced by using a 3D exoscope at the Head and Neck Oncological Unit of Candiolo Cancer Institute.

Results

Data regarding time to drain removal, length of hospitalisation, degree of pain experienced, need for opioid drugs during hospitalisation and after discharge, and intra-operative and post-operative complications were collected. All patients were followed for a minimum of 90 days with possible complications evaluated at each post-operative visit. Post-operative outcomes were evaluated at 3 months after surgery.

Conclusions

The current study indicates that VITOM-3D-assisted retroauricular neck surgery (RANS-3D) may be an interesting approach for neck surgery. The hybrid execution of neck dissection under direct and exoscopic vision represents a valid alternative to video-assisted endoscopic- and robot-assisted techniques.

Large eddy simulation as a fast and accurate engineering approach for the simulation of rotary blood pumps

An accurate representation of the flow field in blood pumps is important for the design and optimization of blood pumps. The primary turbulence modeling methods applied to blood pumps have been the Reynolds-averaged Navier-Stokes (RANS) or URANS (unsteady RANS) method. Large eddy simulation (LES) method has been introduced to simulate blood pumps. Nonetheless, LES has not been widely used to assist in the design and optimization of blood pumps to date due to its formidable computational cost. The purpose of this study is to explore the potential of the LES technique as a fast and accurate engineering approach for the simulation of rotary blood pumps. The performance of "Light LES" (using the same time and spatial resolutions as the URANS) and LES in two rotary blood pumps was evaluated by comparing the results with the URANS and extensive experimental results. This study showed that the results of both "Light LES" and LES are superior to URANS, in terms of both performance curves and key flow features. URANS could not predict the flow separation and recirculation in diffusers for both pumps. In contrast, LES is superior to URANS in capturing these flows, performing well for both design and off-design conditions. The differences between the "Light LES" and LES results were relatively small. This study shows that with less computational cost than URANS, "Light LES" can be considered as a cost-effective engineering approach to assist in the design and optimization of rotary blood pumps.

Boundary layer transition modeling on leading edge inflatable kite airfoils

We present a computational fluid dynamic analysis of boundary layer transition on leading edge inflatable kite airfoils used for airborne wind energy generation. Because of the operation in pumping cycles, the airfoil is generally subject to a wide range of Reynolds numbers. The analysis is based on the combination of the shear stress transport turbulence model with the γ - R ˜ e θ t transition model, which can handle the laminar boundary layer and its transition to turbulence. The implementation of both models in OpenFOAM is described. We show a validation of the method for a sailwing (ie, a wing with a membrane) airfoil and an application to a leading edge inflatable kite airfoil. For the sailwing airfoil, the results computed with transition model agree well with the existing low Reynolds number experiment over the whole range of angles of attack. For the leading edge inflatable kite airfoil, the transition modeling has both favorable and unfavorable effects on the aerodynamics. On the one hand, the aerodynamics suffer from the laminar separation. But, on the other hand, the laminar boundary layer thickens slower than the turbulent counterpart, which, in combination with transition, delays the separation. The results also indicate that the aerodynamics of the kite airfoil could be improved by delaying the boundary layer transition during the traction phase and tripping the transition in the retraction phase.

Numerical study of crude oil batch mixing in a long channel

The main objective of this work is to predict the mixing of two different miscible oils in a very long channel. The background to this problem relates to the mixing of heavy and light oil in a pipeline. As a first step, a 2D channel with an aspect ratio of 250 is considered. The batch-mixing of two miscible crude oils with different viscosities and densities is modeled using an unsteady laminar model and unsteady RANS model available in the commercial CFD solver ANSYS-Fluent. For a comparison, a LES model was used for a 3D version of the 2D channel. The distinguishing feature of this work is the Lagrangian coordinate system utilized to set no-slip wall boundary conditions. The global CFD model has been validated against classical analytical solutions. Excellent agreement has been achieved. Simulations were carried out for a Reynolds number of 6300 (calculated using light oil properties) and a Schmidt number of 10 4 . The results show that, in contrast to the unsteady RANS model, the LES and unsteady laminar models produce comparable mixing dynamics for two oils in the channel. Analysis of simulations also shows that, for a channel length of 100 m and a height of 0.4 m, the complete mixing of two oils across the channel has not been achieved. We showed that the mixing zone consists of the three different mixing sub-zones, which have been identified using the averaged mass fraction of the heavy oil along the flow direction. The first sub-zone corresponds to the main front propagation area with a length of several heights of the channel. The second and third sub-zones are characterized by so-called shear-flow-driven mixing due to the Kelvin-Helmholtz vortices occurring between oils in the axial direction. It was observed that the third sub-zone has a steeper mass fraction gradient of the heavy oil in the axial direction in comparison with the second sub-zone, which corresponds to the flow-averaged mass fraction of 0.5 for the heavy oil.

Evaluating hydraulic and disinfection efficiencies of a full-scale ozone contactor using a RANS-based modeling framework

The capability of predicting hydraulic and disinfection efficiencies of ozone disinfection contactors is essential for evaluating existing contactors and improving future designs. Previous attempts based on ideal and non-ideal models for the hydraulics and simplified mechanisms for chemical reaction modeling have resulted in low accuracy and are restricted to contactors with simple geometries. This manuscript develops a modeling framework for the ozonation process by combining computational fluid dynamics (CFD) with a kinetics-based reaction modeling for the first time. This computational framework has been applied to the full-scale ozone contactor operated by the City of Tampa Water Department. Flow fields, residence time distribution, ozone concentration distribution, and concentration-contact time (CT) distribution within the contactor have been predicted via the computational framework. The predictions of ozone and bromate concentrations at sample points agree well with physical experimental data measured in the contactor. The predicted CT values at the contactor outlet demonstrate that the disinfection performance of the ozone contactor operated by the City of Tampa Water Department is sufficient to meet regulation requirements. The impact of seasonal flow rate change on disinfection performance is found to be significant and deserves attention during the management and operation of a water treatment plant.

Validation of OpenFoam for heavy gas dispersion applications

In the present paper heavy gas dispersion calculations were performed with OpenFoam. For a wind tunnel test case, numerical data was validated with experiments. For a full scale numerical experiment, a code to code comparison was performed with numerical results obtained from Fluent. The validation was performed in a gravity driven environment (slope), where the heavy gas induced the turbulence. For the code to code comparison, a hypothetical heavy gas release into a strongly turbulent atmospheric boundary layer including terrain effects was selected. The investigations were performed for SF6 and CO2 as heavy gases applying the standard k-ɛ turbulence model. A strong interaction of the heavy gas with the turbulence is present which results in a strong damping of the turbulence and therefore reduced heavy gas mixing. Especially this interaction, based on the buoyancy effects, was studied in order to ensure that the turbulence-buoyancy coupling is the main driver for the reduced mixing and not the global behaviour of the turbulence modelling. For both test cases, comparisons were performed between OpenFoam and Fluent solutions which were mainly in good agreement with each other. Beside steady state solutions, the time accuracy was investigated. In the low turbulence environment (wind tunnel test) which for both codes (laminar solutions) was in good agreement, also with the experimental data. The turbulent solutions of OpenFoam were in much better agreement with the experimental results than the Fluent solutions. Within the strong turbulence environment, both codes showed an excellent comparability.