scholarly journals NEW METHODS FOR STABILIZING RANS TURBULENCE MODELS WITH APPLICATION TO LARGE SCALE BREAKING WAVES

2020 ◽  
pp. 19
Author(s):  
Georgios Azorakos ◽  
Bjarke Eltard Larsen ◽  
David R. Fuhrman
Keyword(s):  
Turbulence Model ◽  
Large Scale ◽  
Eddy Viscosity ◽  
Turbulence Models ◽  
Breaking Waves ◽  
Wide Spread ◽  
Cfd Simulations ◽  
New Methods ◽  
Rans Turbulence Models ◽  
Scale Breaking

Recently, Larsen and Fuhrman (2018) have shown that seemingly all commonly used (both k-omega and k-epsilon variants) two-equation RANS turbulence closure models are unconditionally unstable in the potential flow beneath surface waves, helping to explain the wide-spread over-production of turbulent kinetic energy in CFD simulations, relative to measurements. They devised and tested a new formally stabilized formulation of the widely used k-omega turbulence model, making use of a modified eddy viscosity. In the present work, three new formally-stable k-omega turbulence model formulations are derived and tested in CFD simulations involving the flow and dynamics beneath large-scale plunging breaking waves.Recorded Presentation from the vICCE (YouTube Link): https://youtu.be/T2fFRgq3I8E

scholarly journals Assessment of Numerical Methods for Plunging Breaking Wave Predictions

2021 ◽  
Vol 9 (3) ◽  
pp. 264
Author(s):  
Shanti Bhushan ◽  
Oumnia El Fajri ◽  
Graham Hubbard ◽  
Bradley Chambers ◽  
Christopher Kees
Keyword(s):  
Solitary Wave ◽  
Wave Breaking ◽  
Large Scale ◽  
Turbulence Models ◽  
Navier Stokes ◽  
Wave Crest ◽  
Breaking Wave ◽  
Rans Turbulence Models ◽  
Run Up ◽  
Better Than

This study evaluates the capability of Navier–Stokes solvers in predicting forward and backward plunging breaking, including assessment of the effect of grid resolution, turbulence model, and VoF, CLSVoF interface models on predictions. For this purpose, 2D simulations are performed for four test cases: dam break, solitary wave run up on a slope, flow over a submerged bump, and solitary wave over a submerged rectangular obstacle. Plunging wave breaking involves high wave crest, plunger formation, and splash up, followed by second plunger, and chaotic water motions. Coarser grids reasonably predict the wave breaking features, but finer grids are required for accurate prediction of the splash up events. However, instabilities are triggered at the air–water interface (primarily for the air flow) on very fine grids, which induces surface peel-off or kinks and roll-up of the plunger tips. Reynolds averaged Navier–Stokes (RANS) turbulence models result in high eddy-viscosity in the air–water region which decays the fluid momentum and adversely affects the predictions. Both VoF and CLSVoF methods predict the large-scale plunging breaking characteristics well; however, they vary in the prediction of the finer details. The CLSVoF solver predicts the splash-up event and secondary plunger better than the VoF solver; however, the latter predicts the plunger shape better than the former for the solitary wave run-up on a slope case.

scholarly journals Assessment of CFD turbulence models for free surface flow simulation and 1-D modelling for water level calculations over a broad-crested weir floodway

Water SA ◽  
2019 ◽  
Vol 45 (3 July) ◽  
Author(s):  
Ahmed M Helmi
Keyword(s):  
Experimental Data ◽  
Free Surface ◽  
Water Level ◽  
Dimensional Analysis ◽  
Turbulence Model ◽  
Turbulence Models ◽  
Cfd Simulations ◽  
One Dimensional ◽  
The Mean ◽  
The One

Floodways, where a road embankment is permitted to be overtopped by flood water, are usually designed as broad-crested weirs. Determination of the water level above the floodway is crucial and related to road safety. Hydraulic performance of floodways can be assessed numerically using 1-D modelling or 3-D simulation using computational fluid dynamics (CFD) packages. Turbulence modelling is one of the key elements in CFD simulations. A wide variety of turbulence models are utilized in CFD packages; in order to identify the most relevant turbulence model for the case in question, 96 3-D CFD simulations were conducted using Flow-3D package, for 24 broad-crested weir configurations selected based on experimental data from a previous study. Four turbulence models (one-equation, k-ε, RNG k-ε, and k-ω) ere examined for each configuration. The volume of fluid (VOF) algorithm was adopted for free water surface determination. In addition, 24 1-D simulations using HEC-RAS-1-D were conducted for comparison with CFD results and experimental data. Validation of the simulated water free surface profiles versus the experimental measurements was carried out by the evaluation of the mean absolute error, the mean relative error percentage, and the root mean square error. It was concluded that the minimum error in simulating the full upstream to downstream free surface profile is achieved by using one-equation turbulence model with mixing length equal to 7% of the smallest domain dimension. Nevertheless, for the broad-crested weir upstream section, no significant difference in accuracy was found between all turbulence models and the one-dimensional analysis results, due to the low turbulence intensity at this part. For engineering design purposes, in which the water level is the main concern at the location of the flood way, the one-dimensional analysis has sufficient accuracy to determine the water level.

Temperature Corrected Turbulence Model for High Temperature Jet Flow

2004 ◽  
Vol 126 (5) ◽  
pp. 844-850 ◽  
Author(s):  
Khaled S. Abdol-Hamid ◽  
S. Paul Pao ◽  
Steven J. Massey ◽  
Alaa Elmiligui
Keyword(s):  
High Temperature ◽  
Turbulence Model ◽  
Room Temperature ◽  
Eddy Viscosity ◽  
Mixing Layer ◽  
Supersonic Jet ◽  
Turbulence Models ◽  
Jet Flows ◽  
Viscosity Term ◽  
Mixing Layer Flows

It is well known that the two-equation turbulence models under-predict mixing in the shear layer for high temperature jet flows. These turbulence models were developed and calibrated for room temperature, low Mach number, and plane mixing layer flows. In the present study, four existing modifications to the two-equation turbulence model are implemented in PAB3D and their effect is assessed for high temperature jet flows. In addition, a new temperature gradient correction to the eddy viscosity term is tested and calibrated. The new model was found to be in the best agreement with experimental data for subsonic and supersonic jet flows at both low and high temperatures.

Investigation of Annular Jet Flows in a Subsonic Air-Air Ejector With Multi-Ring Entraining Diffuser

2015 ◽  
Author(s):  
A. Namet-Allah ◽  
A. M. Birk
Keyword(s):  
Boundary Conditions ◽  
Turbulence Model ◽  
Nozzle Exit ◽  
Turbulence Models ◽  
Domain Size ◽  
Pressure Recovery ◽  
Cfd Simulations ◽  
The Core ◽  
Measured Flow ◽  
Air Ejector

The current paper presents a cold flow simulation study of a low Mach number air-air ejector with a four ring entraining diffuser that is used in a variety of applications including infrared (IR) suppression of exhaust from helicopters and fixed wing aircraft. The main objectives of this investigation were to identify key issues that must be addressed in successful CFD modelling of such devices, and recognize opportunities to improve the performance of these devices. Two-dimensional CFD simulations were carried out using commercial software, Ansys14. Studies of mesh and domain size sensitivity were made to ensure the CFD results were independent of both factors. A turbulence model independence study using k-ε, k-ω and RSM turbulence models was performed to figure out the appropriate turbulence model that produced the best agreement with the experimental data for several of ejector performance criteria. The measured flow properties in the annulus were used as input boundary conditions for the CFD simulations. However, in the comprehensive turbulence model study, the measured flow parameters at the nozzle exit were also applied as inlet boundary conditions for the CFD simulations. The measured flow velocity at the nozzle exit, at one centerline section inside the mixing tube and at the diffuser exit and the system pressure recovery were compared with the CFD predictions. The ejector pumping ratios, back pressure coefficient and diffuser gap velocities were also compared. It was found that the RANS-based CFD predictions were sensitive to the changes in the ejector domain size, mesh refinement and inlet boundary condition locations. With the annulus inlet boundary conditions, the tested turbulence models under predicted the size of the core separation downstream of the system, back pressure, pumping ratio and pressure recovery in the mixing tube and diffuser. However, the ability of the RNG turbulence model to predict the ejector performance parameters was better than that of the other turbulence models at all inlet flow conditions. Nevertheless, applying the inlet boundary conditions at the nozzle exit enhanced the capability of the RANS-based turbulence model particularly in predicting the ejector pumping ratios, pressure recovery and the size of the core separation. Finally, the acceptable agreement between the experimental data and the CFD predictions provides a valid tool to continue improving these devices using CFD techniques.

Investigation of Inverse-Turbulent-Prandtl Number with Four RNG k-ε Turbulence Models on Compressor Discharge Pipe of Bioenergy Micro Gas Turbine

2016 ◽  
Vol 819 ◽  
pp. 392-400 ◽  
Author(s):  
Ahmad Indra Siswantara ◽  
Budiarso ◽  
Steven Darmawan
Keyword(s):  
Prandtl Number ◽  
Turbulence Model ◽  
Large Scale ◽  
Swirling Flow ◽  
Turbulent Viscosity ◽  
Turbulence Models ◽  
Small Scale ◽  
Turbulent Prandtl Number ◽  
Curved Pipe ◽  
Molecular Viscosity

Inverse-Turbulent Prandtl number (α) is an important parameter in RNG k-ε turbulence models since it affects the ratio of molecular viscosity and turbulent viscosity. In curved pipe, this highly affects the model prediction to a large range eddy-scale flow. According to Yakhot & Orzag, the α range from 1-1.3929 has not been investigated in detail in curved pipe flow (Yakhot & Orszag, 1986) and specific Re. This paper varied inverse-turbulent Prandtl number α to 1-1.3 in RNG k-ε turbulence model on cylindrical curved pipe in order to obtain the optimum value of α to predict unfully-developed flow in the curve with curve ratio R/D of 1.607. Analysis was conducted numericaly with inlet specified Re of 40900 which was generated from the experiment at α 1, 1.1, 1.2, 1.3. Wall surface roughness is not considered in this paper. With assumption that thermal diffusivity is always dominant to turbulent viscosity, higher Inverse-turbulent Prandtl number represent domination of turbulent viscosity to molecular viscosity of the flow and predict to have more interaction between large scale eddy to small scale eddy as well. The results show the use of α = 1.3 has increased the turbulent kinetic energy by 7% and the turbulent dissipation by 5% compared to general inverse-turbulent Prandtl number of 1. The value difference shows that the use of higher α on RNG turbulence model described more interaction between eddies in secondary and swirling flow at pipe curve at Re = 40900.

Comparison of CFD Predictions of Supercritical Carbon Dioxide Axial Flow Turbines Using a Number of Turbulence Models

2021 ◽  
Author(s):  
AbdElRahman AbdElDayem ◽  
Martin T. White ◽  
Abdulnaser I. Sayma
Keyword(s):  
Carbon Dioxide ◽  
Supercritical Carbon Dioxide ◽  
Eddy Viscosity ◽  
Axial Flow ◽  
Turbulence Models ◽  
Loss Assessment ◽  
Cfd Simulations ◽  
Loss Coefficients ◽  
Line Design ◽  
Supercritical Carbon

Abstract A detailed loss assessment of an axial turbine stage operating with a supercritical carbon dioxide (sCO2) based mixture, namely titanium tetrachloride (CO2-TiCl4 85-15%), is presented. To assess aerodynamic losses, computational fluid dynamics (CFD) simulations are conducted using a geometry generated using mean-line design equations which is part of the work delivered to the SCARABEUS project [1]. The CFD simulations are 3D steady state and employ a number of turbulence models to investigate various aerodynamic loss mechanisms. Two categories of turbulence models are used: Eddy Viscosity and Reynold’s Stress models (RSM). The Eddy Viscosity models are the k-ε, k-ε RNG, k-ω, k-ω SST and k-ω Generalized while the RSM models are BSL, LRR, w-RSM and k-ε EARSM. The comparison between different turbulence models showed minor deviations in mass-flow rate, power output and blade loading while significant deviations appear in the loss coefficients and the degree of reaction. It is noted that the k-ε model gives the highest loss coefficients and the lowest isentropic efficiencies while most of the RSM models indicate higher efficiencies and lower loss coefficients. At off-design conditions a sensitivity study revealed that the k-ε RNG model records the sharpest drop in the isentropic efficiency of 8.24% at low mass flowrate reaching 30% off-design. The efficiency sensitivity is found to be less for the other tested models getting 3.1% drop in efficiency for the LRR RSM model.

Turbulence Input Parameters Correction Methodology in Water Lubricated Thrust Bearings

2018 ◽  
Author(s):  
Xin Deng ◽  
Harrison Gates ◽  
Brian Weaver ◽  
Houston Wood ◽  
Roger Fittro
Keyword(s):  
Shear Stress ◽  
Wall Shear Stress ◽  
Film Thickness ◽  
Turbulence Model ◽  
Thermal Deformation ◽  
Eddy Viscosity ◽  
Turbulence Models ◽  
Water Lubrication ◽  
Low Viscosity ◽  
Eddy Viscosity Model

Oil-lubricated bearings are widely used in high speed rotating machines such as those found in the aerospace and automotive industries. However, environmental issues and risk-averse operations are resulting in the removal of oil and the replacement of all sealed oil bearings with reliable water-lubricated bearings. Due to the different fluid properties between oil and water, the low viscosity of water increases Reynolds numbers drastically and therefore makes water-lubricated bearings prone to turbulence effects. This requires finer meshes when compared to oil-lubricated bearings as the low-viscosity fluid produces a very thin lubricant film. Analyzing water-lubricated bearings can also produce convergence and accuracy issues in traditional oil-based analysis codes. Thermal deformation largely affects oil-lubricated bearings, while having limited effects on water lubrication; mechanical deformation largely affects water lubrication, while its effects are typically lower than thermal deformation with oil. One common turbulence model used in these analysis tools is the eddy-viscosity model. Eddy-viscosity depends on the wall shear stress, therefore effective wall shear stress modeling is necessary in determining an appropriate turbulence model. Improving the accuracy and efficiency of modeling approaches for eddy-viscosity in turbulence models is of great importance. Therefore, the purpose of this study is to perform mesh refinement for water-lubricated bearings based on methodologies of eddy-viscosity modeling to improve their accuracy. According to Szeri [1], εm/v for the Boussinesq hypothesis is given by Reichardt’s formula. Fitting the velocity profile with experiments having a y+ in the range of 0–1,000 results in Ng-optimized Reichardt’s constants k = 0.4 and δ+ = 10.7. He clearly states that for y+ > 1000 theoretical predictions and experiments have a greater variance. Armentrout and others [2] developed an equation for δ+ as a function of the pivot Reynolds number, which they validated with CFD simulations. The definition of y+ can be used to approximate the first layer thickness calculated for a uniform mesh. Together with Armentrout’s equation, the number of required elements across the film thickness can be obtained. For typical turbulence models, the y+ must be within a certain range to be accurate. On the condition that the y+ is fixed to that of a standard oil bearing for which an oil bearing code was validated, the number of elements across the film thickness and coefficients used in the eddy-viscosity equation can be adjusted to allow for convergence with other fluids other than that which the traditional oil bearing code was designed for. In this study, the number of required elements across the film for improved prediction quality was calculated based on the proposed eddy-viscosity model mesh correction from the known literature. A comparison between water lubrication using the parameter correction and oil lubrication was also made. The results of this study could aid in improving future designs and models of water-lubricated bearings.

An improved turbulence model for predicting unsteady cavitating flows in centrifugal pump

2015 ◽  
Vol 25 (5) ◽  
pp. 1198-1213 ◽  
Author(s):  
Jian Wang ◽  
Yong Wang ◽  
Houlin Liu ◽  
Haoqin Huang ◽  
Linglin Jiang
Keyword(s):  
Centrifugal Pump ◽  
Turbulence Model ◽  
Eddy Viscosity ◽  
Turbulent Viscosity ◽  
Flow Characteristics ◽  
Turbulence Models ◽  
Content Type ◽  
Cavitation Inception ◽  
Cavitation Performance ◽  
Unsteady Rans

Purpose – The purpose of this paper is to study the unsteady caivitating flows in centrifugal pump, especially for improving the turbulence model to obtain highly resolution results-capable of predicting the cavitation inception, shedding off and collapse procedures. Design/methodology/approach – Both numerical simulations and experimental visualizations were performed in the present paper. An improved RCD turbulence models was proposed by considering three corrected methods: the rotating corrected method, the compressible corrected method and the turbulent viscosity corrected method. Unsteady RANS computations were conducted to compare with the experiments. Findings – The comparison of pump cavitation performance showed that the RCD turbulence model obtained better performance both in non-cavitation and cavitation conditions. The visualization of the cavitation evolution was recorded to validate the unsteady simulations. Good agreement was noticed between calculations and visualizations. It is indicated the RCD model can successfully capture the bubbles detachment and collapse at the rear of the cavity region, since it effectively reduces the eddy viscosity in the multiphase region of liquid and vapor. Furthermore, the eddy viscosity, the instantaneous pressure and density distribution were investigated. The effectiveness of the compressibility was found. Meanwhile, the influence of the rotating corrected method on prediction was explored. It is found that the RCD model solved more unsteady flow characteristics. Originality/value – The current work presented a turbulence model which was much more suitable for predicting the cavitating flow in centrifugal pump.

scholarly journals A Nonlineark-εTurbulence Model Applicable to High Pressure Gradient and Large Curvature Flow

2014 ◽  
Vol 2014 ◽  
pp. 1-9
Author(s):  
Xiyao Gu ◽  
Junlian Yin ◽  
Jintao Liu ◽  
Yulin Wu
Keyword(s):  
High Pressure ◽  
Pressure Gradient ◽  
Turbulence Model ◽  
Turbulence Models ◽  
Wall Region ◽  
Reynolds Stresses ◽  
Wall Functions ◽  
Large Curvature ◽  
Rans Turbulence Models ◽  
Calculation Results

Most of the RANS turbulence models solve the Reynolds stress by linear hypothesis with isotropic model. They can not capture all kinds of vortexes in the turbomachineries. In this paper, an improved nonlineark-εturbulence model is proposed, which is modified from the RNGk-εturbulence model and Wilcox'sk-ωturbulence model. The Reynolds stresses are solved by nonlinear methods. The nonlineark-εturbulence model can calculate the near wall region without the use of wall functions. The improved nonlineark-εturbulence model is used to simulate the flow field in a curved rectangular duct. The results based on the improved nonlineark-εturbulence model agree well with the experimental results. The calculation results prove that the nonlineark-εturbulence model is available for high pressure gradient flows and large curvature flows, and it can be used to capture complex vortexes in a turbomachinery.

Influence of Turbulence Models in Transient Simulations of Inducers Under Steady Operation

2019 ◽  
Author(s):  
Björn Gwiasda ◽  
Matthias Mohr ◽  
Dennis Herrmann-Verspagen ◽  
Martin Böhle
Keyword(s):  
Turbulence Model ◽  
Turbulence Models ◽  
Leading Edge ◽  
Geometrical Parameters ◽  
Transient Effects ◽  
Cfd Simulations ◽  
Transient Simulations ◽  
Models Comparison ◽  
Transient Data ◽  
Additional Influence

Abstract A significant uncertainty in full transient 3D-CFD simulations of inducers is the used turbulence model. To investigate the influence of different turbulence models (Spalart-Allmaras, SST, SSG) on the transient effects experiments and simulations are performed for two inducers. In this paper the results of simulations are represented that are performed under noncavitating conditions to exclude the additional influence of the cavitation model in order to gain an isolated understanding for the influence of turbulence models. Comparison of transient data from simulations as well as experiments determines the influence of the turbulence model. For the investigations two different inducers are designed with different leading edges. One inducer with a straight leading edge and one with a back swept leading edge. All other geometrical parameters are kept constant, such as the hub contour which is purely axial.

Sign in / Sign up

Export Citation Format

Share Document