scholarly journals Incompressible Turbulent Flow Simulation Using theκ-ɛModel and Upwind Schemes

2007 ◽  
Vol 2007 ◽  
pp. 1-26 ◽  
Author(s):  
V. G. Ferreira ◽  
A. C. Brandi ◽  
F. A. Kurokawa ◽  
P. Seleghim Jr. ◽  
A. Castelo ◽  
...  

In the computation of turbulent flows via turbulence modeling, the treatment of the convective terms is a key issue. In the present work, we present a numerical technique for simulating two-dimensional incompressible turbulent flows. In particular, the performance of the high Reynoldsκ-ɛmodel and a new high-order upwind scheme (adaptative QUICKEST by Kaibara et al. (2005)) is assessed for 2D confined and free-surface incompressible turbulent flows. The model equations are solved with the fractional-step projection method in primitive variables. Solutions are obtained by using an adaptation of the front tracking GENSMAC (Tomé and McKee (1994)) methodology for calculating fluid flows at high Reynolds numbers. The calculations are performed by using the 2D version of theFreeflowsimulation system (Castello et al. (2000)). A specific way of implementing wall functions is also tested and assessed. The numerical procedure is tested by solving three fluid flow problems, namely, turbulent flow over a backward-facing step, turbulent boundary layer over a flat plate under zero-pressure gradients, and a turbulent free jet impinging onto a flat surface. The numerical method is then applied to solve the flow of a horizontal jet penetrating a quiescent fluid from an entry port beneath the free surface.

Author(s):  
Yan Jin

Abstract The turbulent flow in a compressor cascade is calculated by using a new simulation method, i.e., parameter extension simulation (PES). It is defined as the calculation of a turbulent flow with the help of a reference solution. A special large-eddy simulation (LES) method is developed to calculate the reference solution for PES. Then, the reference solution is extended to approximate the exact solution for the Navier-Stokes equations. The Richardson extrapolation is used to estimate the model error. The compressor cascade is made of NACA0065-009 airfoils. The Reynolds number 3.82 × 105 and the attack angles −2° to 7° are accounted for in the study. The effects of the end-walls, attack angle, and tripping bands on the flow are analyzed. The PES results are compared with the experimental data as well as the LES results using the Smagorinsky, k-equation and WALE subgrid models. The numerical results show that the PES requires a lower mesh resolution than the other LES methods. The details of the flow field including the laminar-turbulence transition can be directly captured from the PES results without introducing any additional model. These characteristics make the PES a potential method for simulating flows in turbomachinery with high Reynolds numbers.


1986 ◽  
Vol 108 (1) ◽  
pp. 71-75 ◽  
Author(s):  
O. K. Kwon ◽  
R. H. Pletcher

The viscous-inviscid interaction numerical procedure described in Part 1 is used to predict steady, two-dimensional turbulent flow over a rearward-facing step. The accuracy of predictions is observed to be quite sensitive to the specification of length scale in the turbulence modeling. The best results are observed when the length scale is specified algebraically downstream of the step using parameters characteristic of the step geometry. Predictions of mean flow quantities and reattachment length are shown to be in generally good agreement with measurements obtained over a range of channel expansion ratios.


2019 ◽  
Vol 24 (2) ◽  
pp. 387-410
Author(s):  
Md. Shahjada Tarafder ◽  
M. Al Mursaline

Abstract This paper deals with the numerical simulation of a turbulent flow around two-dimensional bodies by the finite volume method with non-orthogonal body-fitted grid. The governing equations are expressed in Cartesian velocity components and solution is carried out using the SIMPLE algorithm for collocated arrangement of scalar and vector variables. Turbulence is modeled by the k- ε turbulence model and wall functions are used to bridge the solution variables at the near wall cells and the corresponding quantities on the wall. A simplified pressure correction equation is derived and proper under-relaxation factors are used so that computational cost is reduced without adversely affecting the convergence rate. The numerical procedure is validated by comparing the computed pressure distribution on the surface of NACA 0012 and NACA 4412 hydrofoils for different angles of attack with experimental data. The grid dependency of the solution is studied by varying the number of cells of the C-type structured mesh. The computed lift coefficients of NACA 4412 hydrofoil at different angles of attack are also compared with experimental results to further substantiate the validity of the proposed methodology.


2005 ◽  
Vol 73 (3) ◽  
pp. 405-412 ◽  
Author(s):  
Hermann F. Fasel ◽  
Dominic A. von Terzi ◽  
Richard D. Sandberg

A flow simulation Methodology (FSM) is presented for computing the time-dependent behavior of complex compressible turbulent flows. The development of FSM was initiated in close collaboration with C. Speziale (then at Boston University). The objective of FSM is to provide the proper amount of turbulence modeling for the unresolved scales while directly computing the largest scales. The strategy is implemented by using state-of-the-art turbulence models (as developed for Reynolds averaged Navier-Stokes (RANS)) and scaling of the model terms with a “contribution function.” The contribution function is dependent on the local and instantaneous “physical” resolution in the computation. This physical resolution is determined during the actual simulation by comparing the size of the smallest relevant scales to the local grid size used in the computation. The contribution function is designed such that it provides no modeling if the computation is locally well resolved so that it approaches direct numerical simulations (DNS) in the fine-grid limit and such that it provides modeling of all scales in the coarse-grid limit and thus approaches a RANS calculation. In between these resolution limits, the contribution function adjusts the necessary modeling for the unresolved scales while the larger (resolved) scales are computed as in large eddy simulation (LES). However, FSM is distinctly different from LES in that it allows for a consistent transition between RANS, LES, and DNS within the same simulation depending on the local flow behavior and “physical” resolution. As a consequence, FSM should require considerably fewer grid points for a given calculation than would be necessary for a LES. This conjecture is substantiated by employing FSM to calculate the flow over a backward-facing step and a plane wake behind a bluff body, both at low Mach number, and supersonic axisymmetric wakes. These examples were chosen such that they expose, on the one hand, the inherent difficulties of simulating (physically) complex flows, and, on the other hand, demonstrate the potential of the FSM approach for simulations of turbulent compressible flows for complex geometries.


Author(s):  
Heinz Herwig ◽  
Yan Jin ◽  
Marc-Florian Uth ◽  
Andrey V. Kuznetsov

Some fundamental issues with respect to turbulent flows through a porous matrix are addressed by analyzing DNS results (DNS: direct numerical simulation, i.e. no turbulence modeling). In a porous matrix with pore sizes of micro or even nano dimensions turbulent flow may occur when the (local) Reynolds number is sufficiently large. An open question, however, is whether turbulence structures are restricted in size by the pore size dimensions or not. This is an important aspect that immediately affects the way turbulence has to be modelled. In order to find out which influence the solid matrix has on the turbulence a generic matrix built from a large number of bars with square cross sections is investigated. Two different DNS approaches are used, a finite volume one and a Lattice-Boltzmann approach. From both DNS calculations detailed flow field information about the influence of the solid matrix on the turbulence structure are obtained. Finally the extension of Darcy’s friction law by the Forchheimer term is investigated with respect to the question whether this extended law may be used in the fully turbulent flow regime.


Author(s):  
Xu Zhang ◽  
Dan Stanescu ◽  
Jonathan W. Naughton

This paper describes a turbulent flow simulation method, which is based on combination of spectral element and large eddy simulation (LES) technique. The robust, high-order discontinuous Galerkin (DG) spectral element method for large-eddy simulation of compressible flows allows for arbitrary order of accuracy and has excellent stability properties. A local spectral discretization in terms of Legendre polynomials is used on each element of the (possibly unstructured) mesh, which allows for high-accurate simulations of turbulent flows. Discontinuities across the interfaces of the elements are resolved using a Riemann solver. An isoparametric representation of the geometry is implemented, with boundaries of the domain discretized to the same order of accuracy as the solution, and explicit low-storage Runge-Kutta methods are used for time integration. Large eddy simulation has proven to be a valuable technique for the calculation of turbulent flows. An element based filtering technique is used in conjunction with the standard Smagorinsky eddy viscosity model to estimate the effect of sub-grid scales stresses in this paper. The recently developed nonlinear model [1] will also be added in the future. The final aim of this project is to use the LES methodology in swirling jet flow simulation. As a first step towards these simulations, simulations of compressible turbulent mixing layer and back-facing step are also performed to evaluate the robust method. Initial results based on both DNS and large eddy simulations are presented in this paper. Future work will be to validate the code.


2006 ◽  
Vol 129 (1) ◽  
pp. 120-128 ◽  
Author(s):  
Zhaogao Luan ◽  
M. M. Khonsari

Turbulent flow inside the seal chamber of a pump operating at high Reynolds number is investigated. The K−ε turbulence model posed in cylindrical coordinates was applied for this purpose. Simulations are performed using the fractional approach method. The results of the computer code are verified by using the FLUENT and by comparing to published results for turbulent Taylor Couette flow. Numerical results of four cases including two rotational speeds with four flush rates are reported. Significant difference between the laminar and the turbulence flow in the seal chamber is predicted. The behavior of the turbulent flows with very high Reynolds number was also investigated. The physical and practical implications of the results are discussed.


A numerical technique is presented for the analysis of turbulent flow associated with combustion. The technique uses Chorin’s random vortex method (r.v.m .), an algorithm capable of tracing the action of elementary turbulent eddies and their cumulative effects without imposing any restriction upon their motion. In the past, the r.v.m . has been used with success to treat non-reacting turbulent flows, revealing in particular the mechanics of large-scale flow patterns, the so-called coherent structures. Introduced here is a flame propagation algorithm , also developed by Chorin, in conjunction with volume sources modelling the mechanical effects of the exothermic process of combustion. As an illustration of its use, the technique is applied to flow in a combustion tunnel w here the flame is stabilized by a back-facing step. Solutions for both non-reacting and reacting flow fields are obtained. Although these solutions are restricted by a set of far-reaching idealizations, they nonetheless mimic quite satisfactorily the essential features of turbulent combustion in a lean propane—air mixture that were observed in the laboratory by means of high speed schlieren photography.


2019 ◽  
Vol 21 (4) ◽  
pp. 632-648 ◽  
Author(s):  
Vesselin Krassimirov Krastev ◽  
Alessandro d’Adamo ◽  
Fabio Berni ◽  
Stefano Fontanesi

A zonal hybridization of the RNG [Formula: see text]-[Formula: see text] URANS model is proposed for the simulation of turbulent flows in internal combustion engines. The hybrid formulation is able to act as URANS, DES or LES in different zones of the computational domain, which are explicitly set by the user. The resulting model has been implemented in a commercial computational fluid dynamics code and the LES branch of the modified RNG [Formula: see text]-[Formula: see text] closure has been initially calibrated on a standard homogeneous turbulence box case. Subsequently, the full zonal formulation has been tested on a fixed intake valve geometry, including comparisons with third-party experimental data. The core of the work is represented by a multi-cycle analysis of the TCC-III experimental engine configuration, which has been compared with the experiments and with prior full-LES computational studies. The applicability of the hybrid turbulence model to internal combustion engine flows is demonstrated, and PIV-like flow statistics quantitatively validate the model performance. This study shows a pioneering application of zonal hybrid models in engine-relevant simulation campaigns, emphasizing the relevance of hybrid models for turbulent engine flows.


Sign in / Sign up

Export Citation Format

Share Document