A Parallelized Coupled Navier-Stokes/Vortex-Panel Solver

2005 ◽  
Vol 127 (4) ◽  
pp. 475-487 ◽  
Author(s):  
Sven Schmitz ◽  
Jean-Jacques Chattot

A Navier-Stokes solver, CFX V5.6, is coupled with an in-house developed Vortex-Panel method for the numerical analysis of wind turbines. The Navier-Stokes zone is confined to the near-field around one wind turbine blade, the Vortex-Panel method models the entire vortex sheet of a two-bladed rotor and accounts for the far-field. This coupling methodology reduces both numerical diffusion and computational cost. The parallelized coupled solver (PCS) is applied to the NREL Phase VI rotor configuration under no-yaw conditions. Fully turbulent flow is assumed using the k - ϵ and k - ω turbulence models. Results obtained are very encouraging for fully attached flow. For separated and partially stalled flow, results are in good agreement with experimental data. Discrepancies observed between the turbulence models are attributed to different prediction of the onset of separation. This is revealed by two-dimensional (2D) results of the S809 airfoil.

Author(s):  
Sarun Benjanirat ◽  
Lakshmi N. Sankar ◽  
Guanpeng Xu

The performance of the NREL Phase VI horizontal axis wind turbine has been studied with a 3-D unsteady Navier-Stokes solver. This solver is third order accurate in space and second order accurate in time, and uses an implicit time marching scheme. Calculations were done for a range of wind conditions from 7 m/s to 25 m/s where the flow conditions ranged from attached flow to massively separated flow. A variety of turbulence models were studied: Baldwin-Lomax Model, Spalart-Allmaras one-equation model, and k-ε two equations model with and without wall corrections. It was found all the models predicted the normal forces and associated bending moments well, but most of them had difficulties in modeling the chord wise forces, power generation, and pitching moments. It was found that the k-ε model with near wall corrections did the best job of predicting most the quantities with acceptable levels of accuracy. Additional studies aimed at transition model development, and grid sensitivity studies in the tip region are deemed necessary to improve the correlation with experiments.


2006 ◽  
Vol 128 (6) ◽  
pp. 1172-1180 ◽  
Author(s):  
Stephen Mahon ◽  
Xin Zhang

The flow around an inverted double-element airfoil in ground effect was studied numerically, by solving the Reynolds averaged Navier-Stokes equations. The predictive capabilities of six turbulence models with regards to the surface pressures, wake flow field, and sectional forces were quantified. The realizable k−ε model was found to offer improved predictions of the surface pressures and wake flow field. A number of ride heights were investigated, covering various force regions. The surface pressures, sectional forces, and wake flow field were all modeled accurately and offered improvements over previous numerical investigations. The sectional forces indicated that the main element generated the majority of the downforce, whereas the flap generated the majority of the drag. The near field and far field wake development was investigated and suggestions concerning reduction of the wake thickness were offered. The main element wake was found to greatly contribute to the overall wake thickness with the contribution increasing as the ride height decreased.


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.


Author(s):  
J. Johansen ◽  
N. N. So̸rensen ◽  
J. A. Michelsen ◽  
S. Schreck

The Detached-Eddy Simulation model implemented in the computational fluid dynamics code, EllipSys3D, is applied on the flow around the NREL Phase-VI wind turbine blade. Results are presented for flow around a parked blade at fixed angle of attack and a blade pitching along the blade axis. Computed blade characteristics are compared with experimental data from the NREL/NASA Ames Phase-VI unsteady experiment. The Detached-Eddy Simulation model is a method for predicting turbulence in computational fluid dynamics computations, which combines a Reynolds Averaged Navier-Stokes method in the boundary layer with a Large Eddy Simulation in the free shear flow. The present study focuses on static and dynamic stall regions highly relevant for stall regulated wind turbines. Computations do predict force coefficients and pressure distributions fairly good and results using Detached-Eddy Simulation show considerably more three-dimensional flow structures compared to conventional two-equation Reynolds Averaged Navier-Stokes turbulence models, but no particular improvements are seen on the global blade characteristics.


Fluids ◽  
2021 ◽  
Vol 6 (8) ◽  
pp. 271
Author(s):  
Dustin Weaver ◽  
Sanja Mišković

This paper presents an analysis of linear viscous stress Favre averaged turbulence models for computational fluid dynamics (CFD) of fully turbulent round jets with a long straight tube geometry in the near field. Although similar work has been performed in the past with very relevant solutions, considerations were not given for the issues and limitations involved with coupling between an Eulerian and Lagrangian phase, such as in fully two-way coupled CFD-DEM. These issues include limitations on solution domain, mesh cell size, wall modelling, and momentum coupling between the two phases in relation to the particles size. Therefore, within these considerations, solutions are provided to the Navier–Stokes equations with various turbulence models using a three-dimensional wedge long straight tube geometry for fully developed turbulence flow. Simulations are performed with a Reynolds number of 13,000 and 51,000 using two different tube diameters. It is found that a modified k-ε turbulence model achieved the most agreeable results for both the velocity and turbulent flow fields between these two flow regimes, while a modified k-ω SST/BSL also provided suitable results.


Author(s):  
Oleg Schilling

Abstract A numerical implementation of a large number of Reynolds-averaged Navier–Stokes (RANS) models based on two-, three-, four-equation, and Reynolds stress turbulence models (using either the turbulent kinetic energy dissipation rate or the turbulent lengthscale) in an Eulerian, finite-difference shock-capturing code is described. The code uses third-order weighted essentially nonoscillatory (WENO) reconstruction of the advective fluxes, and second- or fourth-order central difference derivatives for the computation of spatial gradients. A third-order TVD Runge–Kutta time-evolution scheme is used to evolve the fields in time. Improved closures for the turbulence production terms, compressibility corrections, mixture transport coefficients, and a consistent initialization methodology for the turbulent fields are briefly summarized. The code framework allows for systematic comparisons of detailed predictions from a variety of turbulence models of increasing complexity. Applications of the code with selected K–ε based models are illustrated for each of the three instabilities. Simulations of Rayleigh–Taylor unstable flows for Atwood numbers 0.1–0.9 are shown to be consistent with previous implicit LES (ILES) results and with the expectation of increased asymmetry in the mixing layer characteristics with increasing stratification. Simulations of reshocked Richtmyer–Meshkov turbulent mixing corresponding to experiments with light-to-heavy transition in air/sulfur hexafluoride and incident shock Mach number Mas = 1.50, and heavy-to-light transition in sulfur hexafluoride/air with Mas = 1.45 are shown to be in generally good agreement with both pre- and post-reshock mixing layer widths. Finally, simulations of the seven Brown–Roshko Kelvin–Helmholtz experiments with various velocity and density ratios using nitrogen, helium, and air are shown to give mixing layer predictions in good agreement with data. The results indicate that the numerical algorithms and turbulence models are suitable for simulating these classes of inhomogeneous turbulent flows.


Author(s):  
J Larsson

A full Navier-Stokes solver is used to calculate external heat transfer in two linear two-dimensional turbine cascades, one subsonic and one transonic. Heat transfer results obtained with two low-Reynolds k-∊ models (Chien and Launder-Sharma) and two k-ω models (Wilcox standard and transition) are compared with measurements. Good agreement is found in some regions, but the suction side transition and the leading edge are not predicted correctly. Problems with turbulence levels that are too high in the leading edge region are investigated. The numerical quality of the results is investigated and a few general guidelines about the numerics are given. Grid and scheme independence is also demonstrated.


Fluids ◽  
2019 ◽  
Vol 4 (3) ◽  
pp. 148 ◽  
Author(s):  
Chunhui Zhang ◽  
Charles Patrick Bounds ◽  
Lee Foster ◽  
Mesbah Uddin

In today’s road vehicle design processes, Computational Fluid Dynamics (CFD) has emerged as one of the major investigative tools for aerodynamics analyses. The age-old CFD methodology based on the Reynolds Averaged Navier–Stokes (RANS) approach is still considered as the most popular turbulence modeling approach in automotive industries due to its acceptable accuracy and affordable computational cost for predicting flows involving complex geometries. This popular use of RANS still persists in spite of the well-known fact that, for automotive flows, RANS turbulence models often fail to characterize the associated flow-field properly. It is even true that more often, the RANS approach fails to predict correct integral aerodynamic quantities like lift, drag, or moment coefficients, and as such, they are used to assess the relative magnitude and direction of a trend. Moreover, even for such purposes, notable disagreements generally exist between results predicted by different RANS models. Thanks to fast advances in computer technology, increasing popularity has been seen in the use of the hybrid Detached Eddy Simulation (DES), which blends the RANS approach with Large Eddy Simulation (LES). The DES methodology demonstrated a high potential of being more accurate and informative than the RANS approaches. Whilst evaluations of RANS and DES models on various applications are abundant in the literature, such evaluations on full-car models are relatively fewer. In this study, four RANS models that are widely used in engineering applications, i.e., the realizable k - ε two-layer, Abe–Kondoh–Nagano (AKN) k - ε low-Reynolds, SST k - ω , and V2F are evaluated on a full-scale passenger vehicle with two different front-end configurations. In addition, both cases are run with two DES models to assess the differences between the flow predictions obtained using RANS and DES.


2019 ◽  
Vol 64 (4) ◽  
pp. 1-13
Author(s):  
Jannik Petermann ◽  
Yong Su Jung ◽  
James Baeder ◽  
Jürgen Rauleder

Time-accurate numerical predictions of the interactional aerodynamics between NASA's generic ROBIN fuselage and its four-bladed rotor were performed using the recently developed Reynolds-averaged Navier–Stokes solver HAMSTR. Two stencil-based reconstruction schemes (MUSCL, WENO), a second-order temporal accuracy, and the Spalart–Allmaras turbulence model were used. Three-dimensional volume meshes were created in a robust manner from two-dimensional unstructured surface grids using Hamiltonian paths and strands on nearbody domains. Grid connectivity was established between nearbody and background domains in an overset fashion. Two previously researched operational conditions were reproduced, i. e., a near-hover case and a medium-speed forward flight case at an advance ratio of 0.151. The results were compared with various experimental and numerical references and were found to be in good agreement with both. The comparison included the analysis of the rotor wake structure, tip-vortex trajectories, steady and dynamic fuselage pressure distributions in longitudinal and lateral directions, and rotor inflow predictions.


Author(s):  
N. N. So̸rensen ◽  
J. A. Michelsen ◽  
S. Schreck

The application of an incompressible Reynolds Averaged Navier-Stokes solver to cases from the NREL/NASA Ames wind tunnel test is described. Six cases of the NREL PHASE-VI rotor in the upwind configuration under zero yaw and zero degrees tip pitch are computed. Favorable comparison of the computed results with measurements in the form of shaft torque, root moments, spanwise force distributions, and pressure distributions are shown. The good agreement documents the feasibility of 3D CFD computations in connection with prediction of the performance of new rotors. Additionally it is shown how CFD computations can be used to determine the three dimensional effects in rotor flows.


Sign in / Sign up

Export Citation Format

Share Document