Viscous Flow Computations on a Smooth Cylinders: A Detailed Numerical Study With Validation

Author(s):  
Guilherme Vaz ◽  
Christophe Mabilat ◽  
Remmelt van der Wal ◽  
Paul Gallagher

The objective of this paper is to investigate several numerical and modelling features that the CFD community is currently using to compute the flow around a fixed smooth circular cylinder. Two high Reynolds numbers, 9 × 104 and 5 × 105, are chosen which are in the so called drag-crisis region. Using a viscous flow solver, these features are assessed in terms of quality by comparing the numerical results with experimental data. The study involves grid sensitivity, time step sensitivity, the use of different turbulence models, three-dimensional effects, and a RANS/DES (Reynolds Averaged Navier Stokes, Detached Eddy Simulation) comparison. The resulting drag forces and Strouhal numbers are compared with experimental data of different sources. Major flow features such as velocity and vorticity fields are presented. One of the main conclusions of the present study is that all models predict forces which are far from the experimental values, particularly for the higher Reynolds numbers in the drag-crisis region. Three-dimensional and unsteadiness effects are present, but are only fully captured by sophisticated turbulence models or by DES. DES seems to be the key to better solve the flow problem and obtain better agreement with experimental data. However, its considerable computational demands still do not allow to use it for engineering design purposes.

2019 ◽  
Author(s):  
Michael Thome ◽  
Jens Neugebauer ◽  
Ould el Moctar

Abstract The assessment of design loads acting on Liquefied Natural Gas (LNG) pump tower are widely based on Morison equation. However, the Morison equation lacks consideration of transverse flow, impact loads and the interaction between fluid and structure. Studies dealing with a direct simulation of LNG pump tower loads by means of Computational Fluid Dynamics (CFD), which can cover the aforementioned effects, are currently not available. A comparative numerical study on LNG pump tower loads is presented in this paper focusing on the following two questions: Are impact loads relevant for the structural design of LNG pump towers? In which way does the fluid-structure interaction influence the loads? Numerical simulations of the multiphase problem were conducted using field methods. Firstly, Unsteady Reynolds-Averaged Navier-Stokes (URANS) equations, extended by the Volume of Fluid (VoF) approach were used to simulate the flow inside a three-dimensional LNG tank in model scale without tower structure. The results were used to validate the numerical model against model tests. Motion periods and amplitudes were systematically varied. Velocities and accelerations along the positions of the main structural members of the pump tower were extracted and used as input data for load approximations with the Morison equation. Morison equation, URANS and Delayed Detached Eddy Simulation (DDES) computed tower loads were compared. Time histories as well as statistically processed data were used. Global loads acting on the full (with tower structure) and simplified structure (no tower structure, but using Morison equation) are in the same order of magnitude. However, their time evolution is different, especially at peaks, which is considered significant for the structural design.


2002 ◽  
Vol 124 (4) ◽  
pp. 924-932 ◽  
Author(s):  
Scott Morton ◽  
James Forsythe ◽  
Anthony Mitchell ◽  
David Hajek

An understanding of vortical structures and vortex breakdown is essential for the development of highly maneuverable vehicles and high angle of attack flight. This is primarily due to the physical limits these phenomena impose on aircraft and missiles at extreme flight conditions. Demands for more maneuverable air vehicles have pushed the limits of current CFD methods in the high Reynolds number regime. Simulation methods must be able to accurately describe the unsteady, vortical flowfields associated with fighter aircraft at Reynolds numbers more representative of full-scale vehicles. It is the goal of this paper to demonstrate the ability of detached-eddy Simulation (DES), a hybrid Reynolds-averaged Navier-Stokes (RANS)/large-eddy Simulation (LES) method, to accurately predict vortex breakdown at Reynolds numbers above 1×106. Detailed experiments performed at Onera are used to compare simulations utilizing both RANS and DES turbulence models.


Processes ◽  
2020 ◽  
Vol 8 (12) ◽  
pp. 1687
Author(s):  
Chao Yu ◽  
Xiangyao Xue ◽  
Kui Shi ◽  
Mingzhen Shao ◽  
Yang Liu

This paper compares the performances of three Computational Fluid Dynamics (CFD) turbulence models, Reynolds Average Navier-Stokes (RANS), Detached Eddy Simulation (DES), and Large Eddy Simulation (LES), for simulating the flow field of a wheel loader engine compartment. The distributions of pressure fields, velocity fields, and vortex structures in a hybrid-grided engine compartment model are analyzed. The result reveals that the LES and DES can capture the detachment and breakage of the trailing edge more abundantly and meticulously than RANS. Additionally, by comparing the relevant calculation time, the feasibility of the DES model is proved to simulate the three-dimensional unsteady flow of engine compartment efficiently and accurately. This paper aims to provide a guiding idea for simulating the transient flow field in the engine compartment, which could serve as a theoretical basis for optimizing and improving the layout of the components of the engine compartment.


2018 ◽  
Vol 42 (3) ◽  
pp. 243-251 ◽  
Author(s):  
Edison H Caicedo ◽  
Muhammad S Virk

This article describes a multiphase computational fluid dynamics–based numerical study of the aeroacoustics response of symmetric and asymmetric wind turbine blade profiles in both normal and icing conditions. Three different turbulence models (Reynolds-averaged Navier–Stokes, detached eddy simulation, and large eddy simulation) have been used to make a comparison of numerical results with the experimental data, where a good agreement is found between numerical and experimental results. Detached eddy simulation turbulence model is found suitable for this study. Later, an extended computational fluid dynamics–based aeroacoustics parametric study is carried out for both normal (clean) and iced airfoils, where the results indicate a significant change in sound levels for iced profiles as compared to clean.


2005 ◽  
Vol 127 (2) ◽  
pp. 214-222 ◽  
Author(s):  
F. Bertagnolio ◽  
N. N. Sørensen ◽  
F. Rasmussen

The objective of this paper is an improved understanding of the physics of the aeroelastic motion of wind turbine blades in order to improve the numerical models used for their design. Two- and three-dimensional Navier–Stokes calculations of the flow around a wind turbine airfoil using the k−ω SST and Detached Eddy Simulation (DES) turbulence models, as well as an engineering semiempirical dynamic stall model, are conducted. The computational results are compared to the experimental results that are available for both the static airfoil and the pitching airfoil. It is shown that the Navier–Stokes simulations can reproduce the main characteristic features of the flow. The DES model seems to be able to reproduce most of the details of the unsteady aerodynamics. Aerodynamic work computations indicate that a plunging motion of the airfoil can become unstable.


Author(s):  
Yannis Kallinderis ◽  
Hyung Taek Ahn

Numerical prediction of vortex-induced vibrations requires employment of the unsteady Navier-Stokes equations. Current Navier-Stokes solvers are quite expensive for three-dimensional flow-structure applications. Acceptance of Computational Fluid Dynamics as a design tool for the offshore industry requires improvements to current CFD methods in order to address the following important issues: (i) stability and computation cost of the numerical simulation process, (ii) restriction on the size of the allowable time-step due to the coupling of the flow and structure solution processes, (iii) excessive number of computational elements for 3-D applications, and (iv) accuracy and computational cost of turbulence models used for high Reynolds number flow. The above four problems are addressed via a new numerical method which employs strong coupling between the flow and the structure solutions. Special coupling is also employed between the Reynolds-averaged Navier-Stokes equations and the Spalart-Allmaras turbulence model. An element-type independent spatial discretization scheme is also presented which can handle general hybrid meshes consisting of hexahedra, prisms, pyramids, and tetrahedral.


1997 ◽  
Vol 119 (1) ◽  
pp. 77-84 ◽  
Author(s):  
R. S. Abhari ◽  
M. Giles

An unsteady, compressible, two-dimensional, thin shear layer Navier–Stokes solver is modified to predict the motion-dependent unsteady flow around oscillating airfoils in a cascade. A quasi-three-dimensional formulations is used to account for the stream-wise variation of streamtube height. The code uses Ni’s Lax–Wendroff algorithm in the outer region, an implicit ADI method in the inner region, conservative coupling at the interface, and the Baldwin–Lomax turbulence model. The computational mesh consists of an O-grid around each blade plus an unstructured outer grid of quadrilateral or triangular cells. The unstructured computational grid was adapted to the flow to better resolve shocks and wakes. Motion of each airfoil was simulated at each time step by stretching and compressing the mesh within the O-grid. This imposed motion consists of harmonic solid body translation in two directions and rotation, combined with the correct interblade phase angles. The validity of the code is illustrated by comparing its predictions to a number of test cases, including an axially oscillating flat plate in laminar flow, the Aeroelasticity of Turbomachines Symposium Fourth Standard Configuration (a transonic turbine cascade), and the Seventh Standard Configuration (a transonic compressor cascade). The overall comparison between the predictions and the test data is reasonably good. A numerical study on a generic transonic compressor rotor was performed in which the impact of varying the amplitude of the airfoil oscillation on the normalized predicted magnitude and phase of the unsteady pressure around the airfoil was studied. It was observed that for this transonic compressor, the nondimensional aerodynamic damping was influenced by the amplitude of the oscillation.


Author(s):  
A. A. Boretti

The paper presents a computer code for steady and unsteady, three-dimensional, compressible, turbulent, viscous flow simulations. The mathematical model is based on the Favre-averaged Navier-Stokes conservation equations, closed by a statistical model of turbulence. Turbulence effects are represented by using a low Reynolds number K-ω model. The numerical model makes use of a finite difference approximation in generalized coordinates for space discretization. The solution of time-dependent, three-dimensional, non-homogeneous, partial differential equations is obtained by solving, in a prescribed, symmetric pattern, three time-dependent, one-dimensional, homogeneous partial differential equations, representing convection and diffusion along each generalized coordinate direction, and one ordinary differential equation, representing generation and destruction. An explicit, multi-step, dissipative, Runge-Kutta scheme is finally adopted for time discretization. The code is applied to simulate the flow through a linear cascade of turbine rotor blades, where detailed experimental data are available. Blade aerodynamic and heat transfer are computed, at variable Reynolds and Mach numbers and turbulence levels, and compared with experimental data. While the aerodynamic prediction is relatively unaffected by the properties of both mathematical and numerical models, the heat transfer prediction proves to be extremely sensitive to models details. Low Reynolds number K-ω turbulence models theoretically reproduce laminar, turbulent and transitional boundary layers. However, their practical use in a Navier-Stokes code does not allow to entirely capture the effects of turbulence intensity and Mach and Reynolds numbers on blade heat transfer.


Mathematics ◽  
2022 ◽  
Vol 10 (2) ◽  
pp. 272
Author(s):  
Chenyu Wu ◽  
Haoran Li ◽  
Yufei Zhang ◽  
Haixin Chen

The accuracy of an airfoil stall prediction heavily depends on the computation of the separated shear layer. Capturing the strong non-equilibrium turbulence in the shear layer is crucial for the accuracy of a stall prediction. In this paper, different Reynolds-averaged Navier–Stokes turbulence models are adopted and compared for airfoil stall prediction. The results show that the separated shear layer fixed k−v2¯−ω (abbreviated as SPF k−v2¯−ω) turbulence model captures the non-equilibrium turbulence in the separated shear layer well and gives satisfactory predictions of both thin-airfoil stall and trailing-edge stall. At small Reynolds numbers (Re~105), the relative error between the predicted CL,max of NACA64A010 by the SPF k−v2¯−ω model and the experimental data is less than 3.5%. At high Reynolds numbers (Re~106), the CL,max of NACA64A010 and NACA64A006 predicted by the SPF k−v2¯−ω model also has an error of less than 5.5% relative to the experimental data. The stall of the NACA0012 airfoil, which features trailing-edge stall, is also computed by the SPF k−v2¯−ω model. The SPF k−v2¯−ω model is also applied to a NACA0012 airfoil, which features trailing-edge stall and an error of CL relative to the experiment at CL>1.0 is smaller than 3.5%. The SPF k−v2¯−ω model shows higher accuracy than other turbulence models.


Author(s):  
Mehran Masoumifar ◽  
Suyash Verma ◽  
Arman Hemmati

Abstract This study evaluates how Reynolds-Averaged-Navier-Stokes (RANS) models perform in simulating the characteristics of mean three-dimensional perturbed flows in pipes with targeted wall-shapes. Capturing such flow features using turbulence models is still challenging at high Reynolds numbers. The principal objective of this investigation is to evaluate which of the well-established RANS models can best predict the flow response and recovery characteristics in perturbed pipes at moderate and high Reynolds numbers (10000-158000). First, the flow profiles at various axial locations are compared between simulations and experiments. This is followed by assessing the well-known mean pipeflow scaling relations. The good agreement between our computationally predicted data using Standard k-epsilon model and those of experiments indicated that this model can accurately capture the pipeflow characteristics in response to introduced perturbation with smooth sinusoidal axial variations.


Sign in / Sign up

Export Citation Format

Share Document