A new non-linear RANS model with enhanced near-wall treatment of turbulence anisotropy

2020 ◽  
Vol 82 ◽  
pp. 293-313
Author(s):  
H. Fadhila ◽  
H. Medina ◽  
S. Aleksandrova ◽  
S. Benjamin
Author(s):  
Benjamin H. Taylor ◽  
Tausif Jamal ◽  
D. Keith Walters

Abstract The presence of complex vortical structures, unsteady wakes, separated shear layers, and streamline curvature pose considerable challenges for traditional linear Eddy-Viscosity (LEV) models. Since Non-Linear Eddy Viscosity Models (NEV) models contain additional strain-rate and vorticity relationships, they can provide a better description for flows with Reynolds stress anisotropy and can be considered to be suitable alternatives to traditional EVMs in some cases. In this study, performance of a Non-Linear Explicit Algebraic Reynolds Stress Model (NEARSM) to accurately resolve flow over a surface mounted cube and a 3D axisymmetric hill is evaluated against existing experimental and numerical studies. Numerical simulations were performed using the SST k-ω RANS model, SST k-ω-NEARSM, SST-Multiscale LES model, and two variants of the Dynamic Hybrid RANS-LES (DHRL) model that include the SST k-ω and the SST k-ω-NEARSM as the RANS models. Results indicate that the SST k-ω RANS model fails to accurately predict the flowfield in the separated wake region and although the SST-NEARSM and SST-Multiscale LES models provide an improved description of the flow, they suffer from incorrect RANS-LES transition caused by Modeled Stress Depletion (MSD) and sensitivity to changes in grid resolution. The SST-DHRL and the SST-NEARSM-DHRL variants provide the best agreement to experimental and numerical data.


Author(s):  
Suad Jakirlic´ ◽  
Bjo¨rn Kniesner ◽  
Sanjin Sˇaric´ ◽  
Kemal Hanjalic´

A method of coupling a low-Reynolds-number k–ε RANS (Reynolds-Averaged Navier-Stokes) model with Large-Eddy Simulation (LES) in a two-layer Hybrid LES/RANS (HLR) scheme is proposed in the present work. The RANS model covers the near-wall region and the LES model the remainder of the flow domain. Two different subgrid-scale (SGS) models in LES were considered, the Smagorinsky model and the one-equation model for the residual kinetic energy (Yoshizawa and Horiuti, 1985), combined with two versions of the RANS ε equation, one governing the “isotropic” (ε˜; Chien, 1982) and the other the “homogeneous” dissipation rate (εh; Jakirlic and Hanjalic, 2002). Both fixed and self-adjusting interface locations were considered. The exchange of the variables across the interface was adjusted by smoothing the turbulence viscosity either by adjusting the RANS model parameters, such as Cμ (Temmerman et al., 2005), or by applying an additional forcing at the interface using a method of digital-filter-based generation of inflow data for spatially developing DNS and LES due to Klein et al. (2003). The feasibility of the method was illustrated against the available DNS, fine- and coarse grid LES, DES (Detached Eddy Simulation) and experiments in turbulent flow over a backward-facing step at a low (Yoshioka et al., 2001) and a high Re number (Vogel and Eaton, 1985), periodic flow over a series of 2-D hills (Fro¨hlich et al., 2005) and in a high-Re flow over a 2-D, wall-mounted hump (Greenblat et al, 2004). Prior to these computations, the method was validated in a fully-developed channel flow at a moderate Reynolds number Rem ≈ 24000 (Abe et al., 2004).


Author(s):  
S. Y. Jaw ◽  
R. R. Hwang

To improve the prediction of turbulent flows, a two-scale, non-linear Reynolds stress turbulence model is proposed in this study. It is known that for the near-wall low-Reynolds number turbulent flows, the Kolmogorov turbulence scale, based on the fluid kinematic viscosity and dissipation rate of turbulent kinetic energy (ν,ε), is the dominant turbulence scale, hence it is adopted to address the viscous effects and the rapid increase of dissipation rate in the near wall region. As a wall is approached, the turbulence scale transits smoothly from turbulent kinetic energy based (k, ε) scale to (ν,ε) scale. The damping functions of the low-Reynolds number models can thus be simplified and the near-wall turbulence characteristics, such as the ε distribution, are correctly reproduced. Furthermore, to improve the prediction of the anisotropic Reynolds stresses for complex flows, a nonlinear algebraic Reynolds stress model is incorporated. The same turbulence scales are adopted in the nonlinear algebraic Reynolds stress model. The developed two-scale non-linear Reynolds stress model is first calibrated with the DNS budgets of two-dimensional channel flows, and then applied to predict the separation flow behind a backward facing step. It is found that the proposed two-scale nonlinear Reynolds stress turbulence model is capable of providing satisfactory results without increasing much computation efforts or causing numerical stability problems.


Author(s):  
Alessandro Corsini ◽  
Franco Rispoli

A computational study is presented which investigates the predictive performance of two non-linear turbulence closures in simulating the physics pertinent to decelerating turbomachinery flows. The compared approaches are a cubic non-linear k-ε model and an algebraic Reynolds stress model. They have been considered as promising closures for improving the industry CFD state-of-the-art accounting for non-equilibrium effects. The authors adopt a parallel multi-grid algorithm, which is developed with a finite element formulation based on a highly accurate stabilized Petrov-Galerkin method. The finite element formulation is here applied on equal-order Q1-Q1 as well as mixed Q2-Q1 element pairs, and the accuracy of the latter approximation is assessed on near-wall flows simulation. The parallel solution algorithm for Reynolds Averaged Navier-Stokes modeling exploits an overlapping domain decomposition technique based on an “inexact explicit non linear Schwarz method”. The compressor flow considered for model benchmarking is highly challenging because of the transitional nature of the flow and the existence of significant leading- and trailing-edge separations. The potential of non-isotropic closures has been investigated. The algebraic stress model is shown to provide a better base-line for non-equilibrium effects simulation with respect to the cubic k-ε model. As it is shown for the studied compressor cascade, the cubic eddy-viscosity model exhibits some predictive weaknesses, among them an excessive turbulence attenuation that results in un-realistically delayed transition to turbulence.


Author(s):  
H. Sajjadi ◽  
M. Salmanzadeh ◽  
G. Ahmadi ◽  
S. Jafari

In this study the hybrid RANS/LES turbulence model within the framework of the Lattice Boltzmann method (LBM) was used to study turbulent indoor airflows. In this approach the near wall region was simulated by the RANS model, while the bulk of the domain was analyzed using the LES model with the LBM approach. In the near wall layer where RANS was used, the k-ε turbulence model was employed. For the k-ε turbulence model in conjunction with the LBM two population balance equations for k and ε were used. The present simulation results for the airflow showed good agreement with the experimental data and the earlier numerical results for the hybrid RANS/LES. The results showed that the hybrid model properly predicted the large scale turbulence fluctuation velocities in the bulk of the flow region. In addition, the computational time for the hybrid model is less than that of the LES method.


2001 ◽  
Author(s):  
A. Azzi ◽  
D. Lakehal

Abstract The paper exposes some recent new trends in modelling jets-in-crossflow with relevance to film-cooling of turbine blades. The aim is to compare two classes of turbulence models with respect to their predictive performance in reproducing flow physics. The study focuses on anisotropic eddyviscosity/diffusivity models and explicit algebraic stress models, up to cubic fragments of strain and vorticity tensors. The first class of models are DNS-based two-layer approaches transcending the conventional k–ε model by means of a non-isotropic representation of the turbulent transport coefficients; this is employed in connection with a near-wall one-equation model resolving the semi-viscous sublayer. The aspects of this new strategy are based on known DNS statistics of channel flows and boundary layers. The other class of models are quadratic and cubic explicit algebraic stress formulations rigorously derived from second-moment closures. The stress-strain relations are solved in the context of a two-layer strategy resolving the near-wall region by means of a non-linear one-equation model; the outer core flow is treated by use of the two-equation model. The models are tested for the film cooling of a flat plate, and are then extended to film cooling of a symmetrical turbine blade by a row of laterally injected jets. Comparison of the calculated and measured wall-temperature distributions shows that only the anisotropic eddy viscosity/diffusivity model can correctly predict the spanwise spreading of the temperature field and reduces the strength of the secondary vortices. The non-linear algebraic stress models were of a mixed quality in film cooling calculations.


2021 ◽  
Vol 5 ◽  
pp. 50-65
Author(s):  
Wei Sun ◽  
Liping Xu

In this paper, an investigation into the effect of explicit non-linear turbulence modelling on anisotropic turbulence flows is presented. Such anisotropic turbulence flows are typified in the corner separations in turbomachinery. The commonly used Reynolds-Averaged Navier-Stokes (RANS) turbulence closures, in which the Reynolds stress tensor is modelled by the Boussinesq (linear) constitutive relation with the mean strain-rate tensor, often struggle to predict corner separation with reasonable accuracy. The physical reason for this modelling deficiency is partially attributable to the Boussinesq hypothesis which does not count for the turbulence anisotropy, whilst in a corner separation, the flow is subject to three-dimensional (3D) shear and the effects due to turbulence anisotropy may not be ignored. In light of this, an explicit non-linear Reynolds stress-strain constitutive relation developed by Menter et al. is adopted as a modification of the Reynolds-stress anisotropy. Coupled with the Menter’s hybrid "k-ω" ⁄"k-ε" turbulence model, this non-linear constitutive relation gives significantly improved predictions for the corner separation flows within a compressor cascade, at both the design and off-design flow conditions. The mean vorticity field are studied to further investigate the physical reasons for these improvements, highlighting its potential for the widespread applications in the corner separation prediction.


Sign in / Sign up

Export Citation Format

Share Document