scholarly journals Assessment Study of K-ɛ Turbulence Models and Near-Wall Modeling for Steady State Swirling Flow Analysis in Draft Tube UsingFluent

Author(s):  
Sergio Galván ◽  
Marcelo Reggio ◽  
Francois Guibault
Author(s):  
Sushant Dhiman ◽  
Hosein Foroutan ◽  
Savas Yavuzkurt

One of the major problems encountered in the operation of hydraulic turbines (such as Francis turbines) is a rotating vortex rope which forms in the draft tube under part load conditions. Overall goal of the present research is to study the formation of this oscillating vortex rope using CFD and understand the fundamental mechanisms governing this flow phenomenon. A systematic step by step CFD approach is chosen starting from the simplest to the most complicated flow. The current CFD study reported here therefore aims at studying flows in conical diffusers with and without swirl as a simplified draft tube flow. Two test cases are considered, one is flow with inlet swirl and the other without swirl in a conical diffuser. CFD simulations were carried out using five different turbulence models, namely standard, realizable and RNG k-ε (along with the enhanced wall treatment for near-wall region), SST k-ω and the Reynolds stress model (RSM). Wall pressure coefficient along the diffuser, streamwise and circumferential mean velocity, turbulent kinetic energy (TKE) and Reynolds stress profiles are compared with the experimental data as well as CFD results from literature. It is shown that the moderate levels of swirl cause improvement in the pressure recovery in the diffuser as much as 15%. Also, the standard k-ε and RSM models perform best in predicting turbulent swirling flow behavior. Profiles of the streamwise velocity obtained from these models are in relatively good agreement with the experimental data (with maximum deviation of 25%), while the predictions of the SST k-ω show as much as 60% difference. Also, there is only 8% difference between the level of Reynolds stress obtained from the standard k-ε model and those from the experimental data. Overall, however, all turbulence models need to be improved in order to fully capture the details of the swirling flow in a diffuser and certainly the flow in a draft tube of a hydroturbine where vortex rope breakdown and/or boundary layer separation occurs.


Author(s):  
Karsten Tawackolian ◽  
Martin Kriegel

AbstractThis study looks to find a suitable turbulence model for calculating pressure losses of ventilation components. In building ventilation, the most relevant Reynolds number range is between 3×104 and 6×105, depending on the duct dimensions and airflow rates. Pressure loss coefficients can increase considerably for some components at Reynolds numbers below 2×105. An initial survey of popular turbulence models was conducted for a selected test case of a bend with such a strong Reynolds number dependence. Most of the turbulence models failed in reproducing this dependence and predicted curve progressions that were too flat and only applicable for higher Reynolds numbers. Viscous effects near walls played an important role in the present simulations. In turbulence modelling, near-wall damping functions are used to account for this influence. A model that implements near-wall modelling is the lag elliptic blending k-ε model. This model gave reasonable predictions for pressure loss coefficients at lower Reynolds numbers. Another example is the low Reynolds number k-ε turbulence model of Wilcox (LRN). The modification uses damping functions and was initially developed for simulating profiles such as aircraft wings. It has not been widely used for internal flows such as air duct flows. Based on selected reference cases, the three closure coefficients of the LRN model were adapted in this work to simulate ventilation components. Improved predictions were obtained with new coefficients (LRNM model). This underlined that low Reynolds number effects are relevant in ventilation ductworks and give first insights for suitable turbulence models for this application. Both the lag elliptic blending model and the modified LRNM model predicted the pressure losses relatively well for the test case where the other tested models failed.


Entropy ◽  
2021 ◽  
Vol 23 (7) ◽  
pp. 901
Author(s):  
Davide Bertini ◽  
Lorenzo Mazzei ◽  
Antonio Andreini

Computational Fluid Dynamics is a fundamental tool to simulate the flow field and the multi-physics nature of the phenomena involved in gas turbine combustors, supporting their design since the very preliminary phases. Standard steady state RANS turbulence models provide a reasonable prediction, despite some well-known limitations in reproducing the turbulent mixing in highly unsteady flows. Their affordable cost is ideal in the preliminary design steps, whereas, in the detailed phase of the design process, turbulence scale-resolving methods (such as LES or similar approaches) can be preferred to significantly improve the accuracy. Despite that, in dealing with multi-physics and multi-scale problems, as for Conjugate Heat Transfer (CHT) in presence of radiation, transient approaches are not always affordable and appropriate numerical treatments are necessary to properly account for the huge range of characteristics scales in space and time that occur when turbulence is resolved and heat conduction is simulated contextually. The present work describes an innovative methodology to perform CHT simulations accounting for multi-physics and multi-scale problems. Such methodology, named U-THERM3D, is applied for the metal temperature prediction of an annular aeroengine lean burn combustor. The theoretical formulations of the tool are described, together with its numerical implementation in the commercial CFD code ANSYS Fluent. The proposed approach is based on a time de-synchronization of the involved time dependent physics permitting to significantly speed up the calculation with respect to fully coupled strategy, preserving at the same time the effect of unsteady heat transfer on the final time averaged predicted metal temperature. The results of some preliminary assessment tests of its consistency and accuracy are reported before showing its exploitation on the real combustor. The results are compared against steady-state calculations and experimental data obtained by full annular tests at real scale conditions. The work confirms the importance of high-fidelity CFD approaches for the aerothermal prediction of liner metal temperature.


Author(s):  
Jeffrey D. Ferguson ◽  
Dibbon K. Walters ◽  
James H. Leylek

For the first time in the open literature, code validation quality data and a well-tested, highly reliable computational methodology are employed to isolate the true performance of seven turbulence treatments in discrete jet film cooling. The present research examines both computational and high quality experimental data for two length-to-diameter ratios of a row of streamwise injected, cylindrical film holes. These two cases are used to document the performance of the following turbulence treatments: 1) standard k-ε model with generalized wall functions; 2) standard k-ε model with non-equilibrium wall functions: 3) Renormalization Group k-ε (RNG) model with generalized wall functions; 4) RNG model with non-equilibrium wall functions: 51 standard k-ε model with two-layer turbulence wall treatment; 6) Reynolds Stress Model (RSM) with generalized wall functions; and 7) RSM with non-equilibrium wall functions. Overall, the standard k-ε turbulence model with the two-layer near-wall treatment, which resolves the viscous sublayer, produces results that are more consistent with experimental data.


2017 ◽  
Vol 782 ◽  
pp. 012006
Author(s):  
A Gavrilov ◽  
A Dekterev ◽  
A Minakov ◽  
D Platonov ◽  
A Sentyabov

1999 ◽  
Vol 2 (03) ◽  
pp. 271-280 ◽  
Author(s):  
Ekrem Kasap ◽  
Kun Huang ◽  
Than Shwe ◽  
Dan Georgi

Summary The formation-rate-analysis (FRASM) technique is introduced. The technique is based on the calculated formation rate by correcting the piston rate with fluid compressibility. A geometric factor is used to account for irregular flow geometry caused by probe drawdown. The technique focuses on the flow from formation, is applicable to both drawdown and buildup data simultaneously, does not require long buildup periods, and can be implemented with a multilinear regression, from which near-wellbore permeability, p * and formation fluid compressibility are readily determined. The field data applications indicate that FRA is much less amenable to data quality because it utilizes the entire data set. Introduction A wireline formation test (WFT) is initiated when a probe from the tool is set against the formation. A measured volume of fluid is then withdrawn from the formation through the probe. The test continues with a buildup period until pressure in the tool reaches formation pressure. WFTs provide formation fluid samples and produce high-precision vertical pressure profiles, which, in turn, can be used to identify formation fluid types and locate fluid contacts. Wireline formation testing is much faster compared with the regular pressure transient testing. Total drawdown time for a formation test is just a few seconds and buildup times vary from less than a second (for permeability of hundreds of millidarcy) to half a minute (for permeability of less than 0.1 md), depending on system volume, drawdown rate, and formation permeability. Because WFT tested volume can be small (a few cubic centimeters), the details of reservoir heterogeneity on a fine scale are given with better spatial resolution than is possible with conventional pressure transient tests. Furthermore, WFTs may be preferable to laboratory core permeability measurements since WFTs are conducted at in-situ reservoir stress and temperature. Various conventional analysis techniques are used in the industry. Spherical-flow analysis utilizes early-time buildup data and usually gives permeability that is within an order of magnitude of the true permeability. For p* determination, cylindrical-flow analysis is preferred because it focuses on late-time buildup data. However, both the cylindrical- and spherical-flow analyses have their drawbacks. Early-time data in spherical-flow analysis results in erroneous p* estimation. Late-time data are obtained after long testing times, especially in low-permeability formations; however, long testing periods are not desirable because of potential tool "sticking" problems. Even after extended testing times, the cylindrical-flow period may not occur or may not be detectable on WFTs. When it does occur, permeability estimates derived from the cylindrical-flow period may be incorrect and their validity is difficult to judge. New concepts and analysis techniques, combined with 3-D numerical studies, have recently been reported in the literature.1–7 Three-dimensional numerical simulation studies1–6 have contributed to the diagnosis of WFT-related problems and the improved analysis of WFT data. The experimental studies7 showed that the geometric factor concept is valid for unsteady state probe pressure tests. This study presents the FRA technique8 that can be applied to the entire WFT where a plot for both drawdown and buildup periods renders straight lines with identical slopes. Numerical simulation studies were used to generate data to test both the conventional and the FRA techniques. The numerical simulation data are ideally suited for such studies because the correct answer is known (e.g., the input data). The new technique and the conventional analysis techniques are also applied to the field data and the results are compared. We first review the theory of conventional analysis techniques, then present the FRA technique for combined drawdown and buildup data. A discussion of the numerical results and the field data applications are followed by the conclusions. Analysis Techniques It has been industry practice to use three conventional techniques, i.e., pseudo-steady-state drawdown (PSSDD), spherical and cylindrical-flow analyses, to calculate permeability and p* Conventional Techniques Pseudo-Steady-State Drawdown (PSSDD). When drawdown data are analyzed, it is assumed that late in the drawdown period the pressure drop stabilizes and the system approaches to a pseudo-steady state when the formation flow rate is equal to the drawdown rate. PSSDD permeability is calculated from Darcy's equation with the stabilized (maximum) pressure drop and the flowrate resulting from the piston withdrawal:9–11 $$k {d}=1754.5\left({q\mu \over r {i}\Delta p {{\rm max}}}\right),\eqno ({\rm 1})$$where kd=PSSDD permeability, md. The other parameters are given in Nomenclature.


Author(s):  
Alexander Führing ◽  
Subha Kumpaty ◽  
Chris Stack

In external and internal fluid flow analysis using numerical methods, most attention is paid to the properties of the flow assuming absolute rigidity of the solid bodies involved. However, this is often not the case for water flow or other fluids with high density. The pressure forces cause the geometry to deform which in turn changes the flow properties around it. Thus, a one-way and two-way Fluid-Structure Interaction (FSI) coupling is proposed and compared to a CFD analysis of a windsurfing fin in order to quantify the differences in performance data as well as the properties of the flow. This leads to information about the necessity of the use of FSI in comparison to regular CFD analysis and gives indication of the value of the enhanced results of the deformable analysis applied to water flow around an elastically deformable hydrofoil under different angles of attack. The performance data and flow property evaluation is done in ANSYS Fluent using the k-ω SST and k-ε model with a y+ of 1 and 35 respectively in order to be able to compare the behavior of both turbulence models. It is found that the overall lift coefficient in general is lower and that the flow is less turbulent because of softer transition due to the deformed geometry reducing drag forces. It is also found that the deformation of the tip of the hydrofoil leads to vertical lift forces. For the FSI analysis, one-way and two-way coupling were incorporated leading to the ability to compare results. It has been found that one-way coupling is sufficient as long as there is no stall present at any time.


Sign in / Sign up

Export Citation Format

Share Document