Multigrid Computation of Incompressible Flows Using Two-Equation Turbulence Models: Part I—Numerical Method

1997 ◽  
Vol 119 (4) ◽  
pp. 893-899 ◽  
Author(s):  
X. Zheng ◽  
C. Liao ◽  
C. Liu ◽  
C. H. Sung ◽  
T. T. Huang

A highly efficient numerical approach based on multigrid and preconditioning methods is developed for modeling 3-D incompressible turbulent flows. The incompressible Reynolds-averaged Navier-Stokes equations are written in pseudo-compressibility from, then a preconditioning method is used to reduce the wave speed disparity. The k-ω and k-ε turbulence models are used to estimate the effects of turbulence. The model equations are solved together with the N-S equations in a strongly-coupled way, and all the acceleration techniques originally developed for N-S equations are also used for the turbulence model equations. A point-implicit technique is developed to improve the efficiency of the solution of the turbulence model equations.

1997 ◽  
Vol 119 (4) ◽  
pp. 900-905 ◽  
Author(s):  
X. Zheng ◽  
C. Liao ◽  
C. Liu ◽  
C. H. Sung ◽  
T. T. Huang

In this paper, computational results are presented for three-dimensional high-Reynolds number turbulent flows over a simplified submarine model. The simulation is based on the solution of Reynolds-Averaged Navier-Stokes equations and two-equation turbulence models by using a preconditioned time-stepping approach. A multiblock method, in which the block loop is placed in the inner cycle of a multi-grid algorithm, is used to obtain versatility and efficiency. It was found that the calculated body drag, lift, side force coefficients and moments at various angles of attack or angles of drift are in excellent agreement with experimental data. Fast convergence has been achieved for all the cases with large angles of attack and with modest drift angles.


2014 ◽  
Vol 757 ◽  
Author(s):  
Stephen B. Pope

AbstractFor inhomogeneous turbulent flows at high Reynolds number, it is shown that the redistribution term in Reynolds-stress turbulence models can be determined from the velocity–acceleration correlation. It is further shown that the drift coefficient in the generalized Langevin model (which is used in probability density function (PDF) methods) can be determined from the Reynolds stresses and the velocity–acceleration correlation. These observations are valuable, since the second moments of velocity and acceleration can be measured in experiments, in direct numerical simulations and in well-resolved large-eddy simulations (LES), and hence these turbulence-model quantities can be determined. The redistribution is closely related to the pressure–rate-of-strain, and the unknown in the PDF equation is closely related to the conditional mean pressure gradient (conditional on velocity). In contrast to the velocity–acceleration moments, these pressure statistics are much more difficult to obtain, and our knowledge of them is quite limited. It is also shown that the generalized Langevin model can be re-expressed to provide a direct connection between the drift term and the fluid acceleration. All of these results are first obtained using the constant-property Navier–Stokes equations, but it is then shown that the results are simply extended to variable-density flows.


2009 ◽  
Vol 62 (4) ◽  
Author(s):  
Giancarlo Alfonsi

The approach of Reynolds-averaged Navier–Stokes equations (RANS) for the modeling of turbulent flows is reviewed. The subject is mainly considered in the limit of incompressible flows with constant properties. After the introduction of the concept of Reynolds decomposition and averaging, different classes of RANS turbulence models are presented, and, in particular, zero-equation models, one-equation models (besides a half-equation model), two-equation models (with reference to the tensor representation used for a model, both linear and nonlinear models are considered), stress-equation models (with reference to the pressure-strain correlation, both linear and nonlinear models are considered) and algebraic-stress models. For each of the abovementioned class of models, the most widely-used modeling techniques and closures are reported. The unsteady RANS approach is also discussed and a section is devoted to hybrid RANS/large methods.


2020 ◽  
Vol 9 (1) ◽  
pp. 1402-1419 ◽  
Author(s):  
Nejmeddine Chorfi ◽  
Mohamed Abdelwahed ◽  
Luigi C. Berselli

Abstract In this paper we propose some new non-uniformly-elliptic/damping regularizations of the Navier-Stokes equations, with particular emphasis on the behavior of the vorticity. We consider regularized systems which are inspired by the Baldwin-Lomax and by the selective Smagorinsky model based on vorticity angles, and which can be interpreted as Large Scale methods for turbulent flows. We consider damping terms which are active at the level of the vorticity. We prove the main a priori estimates and compactness results which are needed to show existence of weak and/or strong solutions, both in velocity/pressure and velocity/vorticity formulation for various systems. We start with variants of the known ones, going later on to analyze the new proposed models.


2011 ◽  
Vol 9 (3) ◽  
pp. 740-755 ◽  
Author(s):  
Dinesh A. Shetty ◽  
Jie Shen ◽  
Abhilash J. Chandy ◽  
Steven H. Frankel

AbstractThe rotational incremental pressure-correction (RIPC) scheme, described in Timmermans et al. [Int. J. Numer. Methods. Fluids., 22 (1996)] and Shen et al. [Math. Comput., 73 (2003)] for non-rotational Navier-Stokes equations, is extended to rotating incompressible flows. The method is implemented in the context of a pseudo Fourier-spectral code and applied to several rotating laminar and turbulent flows. The performance of the scheme and the computational results are compared to the so-called diagonalization method (DM) developed by Morinishi et al. [Int. J. Heat. Fluid. Flow., 22 (2001)]. The RIPC predictions are in excellent agreement with the DM predictions, while being simpler to implement and computationally more efficient. The RIPC scheme is not in anyway limited to implementation in a pseudo-spectral code or periodic boundary conditions, and can be used in complex geometries and with other suitable boundary conditions.


2008 ◽  
Vol 18 (08) ◽  
pp. 1443-1479 ◽  
Author(s):  
CHRISTOPHE BERTHON ◽  
FRÉDÉRIC COQUEL

The present work is devoted to an extension of the Navier–Stokes equations where the fluid is governed by two independent pressure laws. Several turbulence models typically enter this framework. The striking novelty over the usual Navier–Stokes equations stems from the impossibility to recast equivalently the system of interest in full conservation form. Opposing to systems of conservation laws, where the end states of the viscous shock are completely characterized by jump relations, the lack of conservation implies the absence of jump relations. We analyze the traveling wave behaviors according to the ratio of viscosities, and we show that the traveling wave solutions of our system tend to the traveling wave solutions of a fully conservative system. This result is used to exhibit asymptotic expansions of the end states. Such an asymptotic behavior achieves a deep physical interpretation when illustrated in the case of compressible turbulent flows.


1994 ◽  
Vol 116 (2) ◽  
pp. 363-369 ◽  
Author(s):  
Y. Mao ◽  
M. Buffat ◽  
D. Jeandel

This paper presents numerical simulations of turbulent flows during the intake and the compression strokes of a model engine. The Favre average Navier-Stokes equations are solved with a k-ε turbulence model. The numerical procedure uses a time dependent semi-implicit scheme and a finite element method with a moving mesh (Buffat, 1991, Mao, 1990). Results of 2-D axisymmetrical calculations with and without inlet swirl are presented and compared to experimental data (Lance et al., 1991). The influence of different turbulence models and the numerical precision of the simulations are also discussed.


Sign in / Sign up

Export Citation Format

Share Document