Direct simulations of low-Reynolds-number turbulent flow in a rotating channel

1993 ◽  
Vol 256 ◽  
pp. 163-197 ◽  
Author(s):  
Reidar Kristoffersen ◽  
Helge I. Andersson
Keyword(s):  
Reynolds Number ◽  
Turbulent Flow ◽  
Stokes Equations ◽  
Navier Stokes ◽  
Complete Suppression ◽  
Reynolds Stresses ◽  
Mean Velocity ◽  
Computational Mesh ◽  
Developed Pressure ◽  
Rotating Channel

Direct numerical simulations of fully developed pressure-driven turbulent flow in a rotating channel have been performed. The unsteady Navier–Stokes equations were written for flow in a constantly rotating frame of reference and solved numerically by means of a finite-difference technique on a 128 × 128 × 128 computational mesh. The Reynolds number, based on the bulk mean velocity Um and the channel half-width h, was about 2900, while the rotation number Ro = 2|Ω|h/Um varied from 0 to 0.5. Without system rotation, results of the simulation were in good agreement with the accurate reference simulation of Kim, Moin & Moser (1987) and available experimental data. The simulated flow fields subject to rotation revealed fascinating effects exerted by the Coriolis force on channel flow turbulence. With weak rotation (Ro = 0.01) the turbulence statistics across the channel varied only slightly compared with the nonrotating case, and opposite effects were observed near the pressure and suction sides of the channel. With increasing rotation the augmentation and damping of the turbulence along the pressure and suction sides, respectively, became more significant, resulting in highly asymmetric profiles of mean velocity and turbulent Reynolds stresses. In accordance with the experimental observations of Johnston, Halleen & Lezius (1972), the mean velocity profile exhibited an appreciable region with slope 2Ω. At Ro = 0.50 the Reynolds stresses vanished in the vicinity of the stabilized side, and the nearly complete suppression of the turbulent agitation was confirmed by marker particle trackings and two-point velocity correlations. Rotational-induced Taylor-Görtler-like counter-rotating streamwise vortices have been identified, and the simulations suggest that the vortices are shifted slightly towards the pressure side with increasing rotation rates, and the number of vortex pairs therefore tend to increase with Ro.

Parameter Extension Simulation of Turbulent Flows in a Compressor Cascade With a High Reynolds Number

2020 ◽  
Author(s):  
Yan Jin
Keyword(s):  
Reynolds Number ◽  
Turbulent Flow ◽  
Turbulent Flows ◽  
Stokes Equations ◽  
Simulation Method ◽  
Potential Method ◽  
Navier Stokes ◽  
Reference Solution ◽  
Compressor Cascade ◽  
High Reynolds

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.

scholarly journals Model-based scaling of the streamwise energy density in high-Reynolds-number turbulent channels

2013 ◽  
Vol 734 ◽  
pp. 275-316 ◽  
Author(s):  
Rashad Moarref ◽  
Ati S. Sharma ◽  
Joel A. Tropp ◽  
Beverley J. McKeon
Keyword(s):  
Reynolds Number ◽  
Stokes Equations ◽  
Reynolds Numbers ◽  
Low Rank ◽  
Navier Stokes ◽  
Mean Velocity ◽  
Navier Stokes Equations ◽  
Self Similar ◽  
The Mean ◽  
High Reynolds

AbstractWe study the Reynolds-number scaling and the geometric self-similarity of a gain-based, low-rank approximation to turbulent channel flows, determined by the resolvent formulation of McKeon & Sharma (J. Fluid Mech., vol. 658, 2010, pp. 336–382), in order to obtain a description of the streamwise turbulence intensity from direct consideration of the Navier–Stokes equations. Under this formulation, the velocity field is decomposed into propagating waves (with single streamwise and spanwise wavelengths and wave speed) whose wall-normal shapes are determined from the principal singular function of the corresponding resolvent operator. Using the accepted scalings of the mean velocity in wall-bounded turbulent flows, we establish that the resolvent operator admits three classes of wave parameters that induce universal behaviour with Reynolds number in the low-rank model, and which are consistent with scalings proposed throughout the wall turbulence literature. In addition, it is shown that a necessary condition for geometrically self-similar resolvent modes is the presence of a logarithmic turbulent mean velocity. Under the practical assumption that the mean velocity consists of a logarithmic region, we identify the scalings that constitute hierarchies of self-similar modes that are parameterized by the critical wall-normal location where the speed of the mode equals the local turbulent mean velocity. For the rank-1 model subject to broadband forcing, the integrated streamwise energy density takes a universal form which is consistent with the dominant near-wall turbulent motions. When the shape of the forcing is optimized to enforce matching with results from direct numerical simulations at low turbulent Reynolds numbers, further similarity appears. Representation of these weight functions using similarity laws enables prediction of the Reynolds number and wall-normal variations of the streamwise energy intensity at high Reynolds numbers (${Re}_{\tau } \approx 1{0}^{3} {\unicode{x2013}} 1{0}^{10} $). Results from this low-rank model of the Navier–Stokes equations compare favourably with experimental results in the literature.

Similarity in non-rotating and rotating turbulent pipe flows

1999 ◽  
Vol 379 ◽  
pp. 1-22 ◽  
Author(s):  
MARTIN OBERLACK
Keyword(s):  
Reynolds Number ◽  
Scaling Laws ◽  
Stokes Equations ◽  
Scaling Law ◽  
Navier Stokes ◽  
Invariant Solutions ◽  
Mean Velocity ◽  
Navier Stokes Equations ◽  
Pipe Flows ◽  
The Mean

The Lie group approach developed by Oberlack (1997) is used to derive new scaling laws for high-Reynolds-number turbulent pipe flows. The scaling laws, or, in the methodology of Lie groups, the invariant solutions, are based on the mean and fluctuation momentum equations. For their derivation no assumptions other than similarity of the Navier–Stokes equations have been introduced where the Reynolds decomposition into the mean and fluctuation quantities has been implemented. The set of solutions for the axial mean velocity includes a logarithmic scaling law, which is distinct from the usual law of the wall, and an algebraic scaling law. Furthermore, an algebraic scaling law for the azimuthal mean velocity is obtained. In all scaling laws the origin of the independent coordinate is located on the pipe axis, which is in contrast to the usual wall-based scaling laws. The present scaling laws show good agreement with both experimental and DNS data. As observed in experiments, it is shown that the axial mean velocity normalized with the mean bulk velocity um has a fixed point where the mean velocity equals the bulk velocity independent of the Reynolds number. An approximate location for the fixed point on the pipe radius is also given. All invariant solutions are consistent with all higher-order correlation equations. A large-Reynolds-number asymptotic expansion of the Navier–Stokes equations on the curved wall has been utilized to show that the near-wall scaling laws for at surfaces also apply to the near-wall regions of the turbulent pipe flow.

Steady Flow Structures and Pressure Drops in Wavy-Walled Tubes

1987 ◽  
Vol 109 (3) ◽  
pp. 255-261 ◽  
Author(s):  
M. E. Ralph
Keyword(s):  
Reynolds Number ◽  
Pressure Gradient ◽  
Stokes Equations ◽  
Reynolds Numbers ◽  
Creeping Flow ◽  
Navier Stokes ◽  
Mean Velocity ◽  
Pressure Drops ◽  
Mean Pressure ◽  
Reynolds Number Flows

Solutions of the Navier-Stokes equations for steady axisymmetric flows in tubes with sinusoidal walls were obtained numerically, for Reynolds numbers (based on the tube radius and mean velocity at a constriction) up to 500, and for varying depth and wavelength of the wall perturbations. Results for the highest Reynolds numbers showed features suggestive of the boundary layer theory of Smith [23]. In the other Reynolds number limit, it has been found that creeping flow solutions can exhibit flow reversal if the perturbation depth is large enough. Experimentally measured pressure drops for a particular tube geometry were in agreement with computed predictions up to a Reynolds number of about 300, where transitional effects began to disturb the experiments. The dimensionless mean pressure gradient was found to decrease with increasing Reynolds number, although the rate of decrease was less rapid than in a straight-walled tube. Numerical results showed that the mean pressure gradient decreases as both the perturbation wavelength and depth increase, with the higher Reynolds number flows tending to be more influenced by the wavelength and the lower Reynolds number flows more affected by the depth.

scholarly journals Supersonic turbulent flow simulation using a scalable parallel modal discontinuous Galerkin numerical method

2019 ◽  
Vol 9 (1) ◽  
Author(s):  
Tomas Houba ◽  
Arnob Dasgupta ◽  
Shivasubramanian Gopalakrishnan ◽  
Ryan Gosse ◽  
Subrata Roy
Keyword(s):  
Reynolds Number ◽  
Turbulent Flow ◽  
Discontinuous Galerkin ◽  
Stokes Equations ◽  
Parallel Computer ◽  
Navier Stokes ◽  
Finite Difference Schemes ◽  
Performance Study ◽  
Navier Stokes Equations ◽  
Flow Simulations

Abstract The scalability and efficiency of numerical methods on parallel computer architectures is of prime importance as we march towards exascale computing. Classical methods like finite difference schemes and finite volume methods have inherent roadblocks in their mathematical construction to achieve good scalability. These methods are popularly used to solve the Navier-Stokes equations for fluid flow simulations. The discontinuous Galerkin family of methods for solving continuum partial differential equations has shown promise in realizing parallel efficiency and scalability when approaching petascale computations. In this paper an explicit modal discontinuous Galerkin (DG) method utilizing Implicit Large Eddy Simulation (ILES) is proposed for unsteady turbulent flow simulations involving the three-dimensional Navier-Stokes equations. A study of the method was performed for the Taylor-Green vortex case at a Reynolds number ranging from 100 to 1600. The polynomial order P = 2 (third order accurate) was found to closely match the Direct Navier-Stokes (DNS) results for all Reynolds numbers tested outside of Re = 1600, which had a normalized RMS error of 3.43 × 10−4 in the dissipation rate for a 603 element mesh. The scalability and performance study of the method was then conducted for a Reynolds number of 1600 for polynomials orders from P = 2 to P = 6. The highest order polynomial that was tested (P = 6) was found to have the most efficient scalability using both the MPI and OpenMP implementations.

scholarly journals Application of Different Turbulence Models for Duct Flow Simulation With Reduced and Full Navier-Stokes Equations

1995 ◽  
Author(s):  
Charles Hirsch ◽  
Andrei E. Khodak
Keyword(s):  
Reynolds Number ◽  
Reynolds Stress ◽  
Stokes Equations ◽  
Flow Simulation ◽  
Turbulence Models ◽  
Navier Stokes ◽  
Duct Flow ◽  
Mean Velocity ◽  
Navier Stokes Equations ◽  
Wall Functions

An S-shaped diffusing duct flow is analysed with various turbulence models: standard high-Reynolds-number k-ε model with wall functions, a low-Reynolds-number k-ε model, and an explicit non-linear algebraic Reynolds stress closure model (ASM). In addition, computations were obtained with parabolic and partially-parabolic approximations along with solutions of the full Navier-Stokes equations. Results of the numerical simulation are compared with LDA measurements of the turbulent flow as reported by Whitelaw and Yu (1993). Detailed comparisons of mean velocity and Reynolds stress distributions are presented. It is shown that the partially-parabolic approach gives a significant improvement over parabolic approximation in predicting mean velocities and pressure distributions. However, the turbulence models used are still not able to reproduce all the observed flow features, although the ASM results appear to be slightly but consistently closer to the experimental data.

Magnetic levitation of liquid metals

1982 ◽  
Vol 117 ◽  
pp. 27-43 ◽  
Author(s):  
A. J. Mestel
Keyword(s):  
Reynolds Number ◽  
Stokes Equations ◽  
Magnetic Levitation ◽  
Liquid Metals ◽  
High Surface ◽  
Spherical Shape ◽  
Magnetic Reynolds Number ◽  
Surface Shape ◽  
Navier Stokes ◽  
Mean Velocity

The process of levitation melting of metals is examined analytically and numerically for the case of axisymmetric toroidal high-frequency currents. The governing equations for the mean-velocity field and associated free-surface shape are derived under the assumption of low magnetic Reynolds number and the neglect of thermal effects. The form of the solution for high Reynolds number is discussed in general, and particularized to the case of high surface tension, in which limit a perturbation analysis about a spherical shape is presented. Finite-difference techniques are used to solve the Navier–Stokes equations in the sphere, and the surface perturbation is calculated. The asymptotic behaviour of the potential vorticity is illustrated by the numerical experiments.

A note on the solution of the Navier-Stokes equations for a spherically symmetric expansion into a very low pressure

1973 ◽  
Vol 59 (2) ◽  
pp. 391-396 ◽  
Author(s):  
N. C. Freeman ◽  
S. Kumar
Keyword(s):  
Shock Wave ◽  
Reynolds Number ◽  
Stokes Equation ◽  
Stokes Equations ◽  
Low Pressure ◽  
Navier Stokes ◽  
Spherically Symmetric ◽  
Navier Stokes Equation ◽  
Navier Stokes Equations ◽  
Type Flow

It is shown that, for a spherically symmetric expansion of a gas into a low pressure, the shock wave with area change region discussed earlier (Freeman & Kumar 1972) can be further divided into two parts. For the Navier–Stokes equation, these are a region in which the asymptotic zero-pressure behaviour predicted by Ladyzhenskii is achieved followed further downstream by a transition to subsonic-type flow. The distance of this final region downstream is of order (pressure)−2/3 × (Reynolds number)−1/3.

scholarly journals Finite-amplitude steady waves in plane viscous shear flows

1985 ◽  
Vol 160 ◽  
pp. 281-295 ◽  
Author(s):  
F. A. Milinazzo ◽  
P. G. Saffman
Keyword(s):  
Boundary Layer ◽  
Reynolds Number ◽  
Stokes Equations ◽  
Reynolds Numbers ◽  
Finite Amplitude ◽  
Linear Instability ◽  
Navier Stokes ◽  
Unidirectional Flow ◽  
Viscous Shear ◽  
Finite Amplitude Waves

Computations of two-dimensional solutions of the Navier–Stokes equations are carried out for finite-amplitude waves on steady unidirectional flow. Several cases are considered. The numerical method employs pseudospectral techniques in the streamwise direction and finite differences on a stretched grid in the transverse direction, with matching to asymptotic solutions when unbounded. Earlier results for Poiseuille flow in a channel are re-obtained, except that attention is drawn to the dependence of the minimum Reynolds number on the physical constraint of constant flux or constant pressure gradient. Attempts to calculate waves in Couette flow by continuation in the velocity of a channel wall fail. The asymptotic suction boundary layer is shown to possess finite-amplitude waves at Reynolds numbers orders of magnitude less than the critical Reynolds number for linear instability. Waves in the Blasius boundary layer and unsteady Rayleigh profile are calculated by employing the artifice of adding a body force to cancel the spatial or temporal growth. The results are verified by comparison with perturbation analysis in the vicinity of the linear-instability critical Reynolds numbers.

Sign in / Sign up

Export Citation Format

Share Document