Three-Dimensional Numerical Simulations of Flows Past Smooth and Rough/Bare and Helically Straked Circular Cylinders Allowed to Undergo Two Degree-of-Freedom Motions

Author(s):  
Juan P. Pontaza ◽  
Raghu G. Menon ◽  
Hamn-Ching Chen

We simulate the flow past smooth and rough rigid circular cylinders that are either bare or outfitted with helical strakes. We consider operating conditions that correspond to high Reynolds numbers of 105 and 106, and allow for two degree-of-freedom motions such that the structure is allowed to respond to flow-induced cross-flow and in-line forces. The computations are performed using a parallelized Navier–Stokes in-house solver using overset grids. For smooth surface simulations at a Reynolds number of 105, we use a Smagorinsky large eddy simulation turbulence model and for the Reynolds number cases of 106 we make use of the unsteady Reynolds-averaged Navier–Stokes equations with a two-layer k-epsilon turbulence model. The rough surface modifications of the two-layer k-epsilon turbulence model due to Durbin et al. (2001, “Rough Wall Modification of Two-Layer k-Epsilon,” ASME J. Fluids Eng., 123, pp. 16–21) are implemented to account for surface roughness effects. In all our computations we aim to resolve the boundary layer directly by using adequate grid spacing in the near-wall region. The predicted global flow parameters under different surface conditions are in good agreement with experimental data, and significant vortex-induced vibration suppression is observed when using helically straked cylinders.

Author(s):  
Juan P. Pontaza ◽  
Raghu G. Menon ◽  
H. C. Chen

We simulate the flow past smooth and rough rigid circular cylinders that are either bare or outfitted with helical stakes. We consider operating conditions that correspond to high Reynolds numbers of 105 and 106, and allow for two degree-of-freedom motions when the structure is allowed to respond to vortex-induced cross flow and in-line forces. The computations are performed using a parallelized Navier-Stokes in-house solver using overset grids. For smooth surface simulations at a Reynolds number of 105, we use a Smagorinsky Large Eddy Simulation (LES) turbulence model and for the Reynolds number cases of 106 we make use of the unsteady Reynolds-Averaged Navier-Stokes (URANS) equations with a two-layer k-epsilon turbulence model. The rough surface modifications of the two-layer k-epsilon turbulence model due to Durbin et al. (ASME J. Fluids Eng., 2001) are implemented to account for surface roughness effects. In all our computations we aim to resolve the boundary layer directly by using adequate grid spacing in the near-wall region. The predicted global flow parameters under different surface conditions are in good agreement with experimental data and significant VIV suppression is observed when using helically straked cylinders.


Author(s):  
Juan P. Pontaza ◽  
Hamn-Ching Chen

In an effort to gain a better understanding of the VIV phenomena, we present three-dimensional numerical simulations of VIV of circular cylinders. We consider operating conditions that correspond to high Reynolds number flow, low structural damping, and allow for two-degree of freedom motion. The numerical implementation makes use of overset (Chimera) grids, in a multiple block environment where the workload associated with the blocks is distributed among multiple processors working in parallel. The three-dimensional grids around the cylinder are allowed to undergo arbitrary motions with respect to fixed background grids, eliminating the need for tedious grid regeneration at every time step.


2006 ◽  
Vol 129 (3) ◽  
pp. 158-164 ◽  
Author(s):  
Juan P. Pontaza ◽  
Hamn-Ching Chen

In an effort to gain a better understanding of vortex-induced vibrations (VIV), we present three-dimensional numerical simulations of VIV of circular cylinders. We consider operating conditions that correspond to a Reynolds number of 105, low structural mass and damping (m*=1.0, ζ*=0.005), a reduced velocity of U*=6.0, and allow for two degree-of-freedom (X and Y) motion. The numerical implementation makes use of overset (Chimera) grids, in a multiple block environment where the workload associated with the blocks is distributed among multiple processors working in parallel. The three-dimensional grid around the cylinder is allowed to undergo arbitrary motions with respect to fixed background grids, eliminating the need for grid regeneration as the structure moves on the fluid mesh.


2021 ◽  
Vol 36 (2) ◽  
pp. 101-115
Author(s):  
Roman S. Solomatin ◽  
Ilya V. Semenov

Abstract Turbulent mixing, ignition, and flame stabilization in the non-premixed supersonic hydrogen-air flow is numerically modelled in a near-wall region. Mixing algorithm based on the turbulence approach SARANS (Reynolds Averaged Navier–Stokes equations closed with Spalart–Allmaras turbulence model) with a diffusion model and a detailed kinetic model for hydrogen-air chemical reactions are employed. The system of governing equations that consists of basic conservation laws and the turbulence model equation is solved in a coupled manner with the LU–SGS–GMRES method. The model is applied to simulate the process of hydrogen injection into a M = 2.44 air flow with their subsequent mixing, ignition, and combustion in the Burrows– Kurkov chamber. The results are compared to available experimental and reference computational data. All calculations are carried out on the ‘MVS-10P’ JSCC RAS supercomputer cluster.


2021 ◽  
pp. 2150430
Author(s):  
Junjie Wu ◽  
Jiahua Li ◽  
Xiang Qiu ◽  
Xilin Xie ◽  
Yulu Liu

To address the closure problem of Reynolds-averaged Navier–Stokes in numerical simulations of turbulence, the method of solving Reynolds-averaged Navier–Stokes equations based on artificial neural network is introduced in this paper. We establish the nonlinear mapping relationship between the average flow field and the steady-state eddy viscosity field. The machine learning (ML) surrogate model for the shear stress transport turbulence model is constructed. The solution process of replacing the original turbulence model equations with the predicted field variables is realized by coupling the ML algorithm with the CFD solver. The classical backward facing step problem is selected in our study to verify the simulation accuracy of the surrogate model. The comparative analysis is carried out on the six backward facing step flows simulations at different Reynolds numbers. The results of simulations show that the testing flows with the Reynolds numbers closest training datasets Reynolds numbers can obtain the best simulation accuracy. Then for the Reynolds number that is lower than the training datasets, the simulation accuracy will decrease as the Reynolds number decreases. On the contrary, the simulation accuracy of the test flow will increase as the Reynolds number increases. These results indicate the feasibility of the ML surrogate model to simulate at higher Reynolds number. It shows the great potential of applying ML algorithms to Reynolds-averaged Navier–Stokes simulation (RANS) turbulence model and also provides a new idea for industrial simulations of turbulent flows.


2015 ◽  
Vol 137 (10) ◽  
Author(s):  
K. Supradeepan ◽  
Arnab Roy

Numerical simulations were performed for two-dimensional viscous incompressible flow past two stationary side-by-side rotating circular cylinders at Reynolds number (Re) 100 by varying center-to-center distance between the cylinders from 1.1 to 3.5 times the diameter (D) of a cylinder and different rotational speed ratio (α) = 0.5, 1.0, and 1.25. The incompressible Navier–Stokes equations were solved using consistent flux reconstruction (CFR) technique of Roy and Bandyopadhyay (2006, “A Finite Volume Method for Viscous Incompressible Flows Using a Consistent Flux Reconstruction Scheme,” Int. J. Numer. Methods Fluids, 52(3), pp. 297–319). Eight different flow regimes were observed within the investigated parametric space. An attempt has been made to characterize the different flow regimes using vorticity contours, λ2 criterion, and force coefficients. All these above stated methods confirm the existence of eight different regimes in the flow.


2021 ◽  
Vol 225 ◽  
pp. 108806
Author(s):  
Qunfeng Zou ◽  
Lin Ding ◽  
Rui Zou ◽  
Hao Kong ◽  
Haibo Wang ◽  
...  

1973 ◽  
Vol 59 (2) ◽  
pp. 391-396 ◽  
Author(s):  
N. C. Freeman ◽  
S. Kumar

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.


2013 ◽  
Vol 135 (7) ◽  
Author(s):  
A. Ghidoni ◽  
A. Colombo ◽  
S. Rebay ◽  
F. Bassi

In the last decade, discontinuous Galerkin (DG) methods have been the subject of extensive research efforts because of their excellent performance in the high-order accurate discretization of advection-diffusion problems on general unstructured grids, and are nowadays finding use in several different applications. In this paper, the potential offered by a high-order accurate DG space discretization method with implicit time integration for the solution of the Reynolds-averaged Navier–Stokes equations coupled with the k-ω turbulence model is investigated in the numerical simulation of the turbulent flow through the well-known T106A turbine cascade. The numerical results demonstrate that, by exploiting high order accurate DG schemes, it is possible to compute accurate simulations of this flow on very coarse grids, with both the high-Reynolds and low-Reynolds number versions of the k-ω turbulence model.


Sign in / Sign up

Export Citation Format

Share Document