Numerical Modeling of Wind Turbine Wakes

2002 ◽  
Vol 124 (2) ◽  
pp. 393-399 ◽  
Author(s):  
Jens No̸rkær So̸rensen ◽  
Wen Zhong Shen

An aerodynamical model for studying three-dimensional flow fields about wind turbine rotors is presented. The developed algorithm combines a three-dimensional Navier-Stokes solver with a so-called actuator line technique in which the loading is distributed along lines representing the blade forces. The loading is determined iteratively using a blade-element approach and tabulated airfoil data. Computations are carried out for a 500 kW Nordtank wind turbine equipped with three LM19.1 blades. The computed power production is found to be in good agreement with measurements. The computations give detailed information about basic features of wind turbine wakes, including distributions of interference factors and vortex structures. The model serves in particular to analyze and verify the validity of the basic assumptions employed in the simple engineering models.

2000 ◽  
Vol 122 (4) ◽  
pp. 653-660 ◽  
Author(s):  
M. Greiner ◽  
R. J. Faulkner ◽  
V. T. Van ◽  
H. M. Tufo ◽  
P. F. Fischer

Navier-Stokes simulations of three-dimensional flow and augmented convection in a channel with symmetric, transverse grooves on two opposite walls were performed for 180⩽Re⩽1600 using the spectral element technique. A series of flow transitions was observed as the Reynolds number was increased, from steady two-dimensional flow, to traveling two and three-dimensional wave structures, and finally to three-dimensional mixing. Three-dimensional simulations exhibited good agreement with local and spatially averaged Nusselt number and friction factor measurements over the range 800⩽Re⩽1600. [S0022-1481(00)00904-X]


2012 ◽  
Vol 2012 ◽  
pp. 1-11 ◽  
Author(s):  
Radoslav Bozinoski ◽  
Roger L. Davis

This paper describes a detached-eddy simulation (DES) for the flow over a wall-mounted hump. The Reynolds number based on the hump chord isRec=9.36×105with an in-let Mach number of 0.1. Solutions of the three-dimensional Reynolds-averaged Navier-Stokes (RANS) procedure are obtained using the Wilcoxk−ωequations. The DES results are obtained using the model presented by Bush and Mani and are compared with RANS solutions and experimental data from NASA's 2004 Computational Fluid Dynamics Validation on Synthetic Jets and Turbulent Separation Control Workshop. The DES procedure exhibited a three-dimensional flow structure in the wake, with a 13.65% shorter mean separation region compared to RANS and a mean reattachment length that is in good agreement with experimental measurements. DES predictions of the pressure coefficient in the separation region also exhibit good agreement with experiment and are more accurate than RANS predictions.


Author(s):  
B. W. Swanson

A modification of Moore’s (12) method has been developed for solving the three-dimensional Navier Stokes equations to calculate the flow in a log spiral impeller. A complete mathematical development of this method is presented. The parabolic finite-difference marching code used to make the calculations is an extensive revision of a CATHY3 code obtained from D. B. Spalding. Calculations are made for Vr, Vθ and Vz on a non-orthogonal grid that is ideally suited for impellers with back-swept blades. An inviscid solution is in good agreement with the elliptic solution and validates the computational procedure. The method is applied to calculate the viscous three-dimensional flow in a log spiral impeller.


Author(s):  
Qiangqiang Huang ◽  
Xinqian Zheng ◽  
Aolin Wang

Air often flows into compressors with inlet prewhirl, because it will obtain a circumferential component of velocity via inlet distortion or swirl generators such as inlet guide vanes. A lot of research has shown that inlet prewhirl does influence the characteristics of components, but the change of the matching relation between the components caused by inlet prewhirl is still unclear. This paper investigates the influence of inlet prewhirl on the matching of the impeller and the diffuser and proposes a flow control method to cure mismatching. The approach combines steady three-dimensional Reynolds-averaged Navier-Stokes (RANS) simulations with theoretical analysis and modeling. The result shows that a compressor whose impeller and diffuser match well at zero prewhirl will go to mismatching at non-zero prewhirl. The diffuser throat gets too large to match the impeller at positive prewhirl and gets too small for matching at negative prewhirl. The choking mass flow of the impeller is more sensitive to inlet prewhirl than that of the diffuser, which is the main reason for the mismatching. To cure the mismatching via adjusting the diffuser vanes stagger angle, a one-dimensional method based on incidence matching has been proposed to yield a control schedule for adjusting the diffuser. The optimal stagger angle predicted by analytical method has good agreement with that predicted by computational fluid dynamics (CFD). The compressor is able to operate efficiently in a much broader flow range with the control schedule. The flow range, where the efficiency is above 80%, of the datum compressor and the compressor only employing inlet prewhirl and no control are just 25.3% and 31.8%, respectively. For the compressor following the control schedule, the flow range is improved up to 46.5%. This paper also provides the perspective of components matching to think about inlet distortion.


2014 ◽  
Vol 694 ◽  
pp. 187-192
Author(s):  
Jin Xiang Wu ◽  
Jian Sun ◽  
Xiang Gou ◽  
Lian Sheng Liu

The three-dimensional coupled explicit Reynolds Averaged Navier–Stokes (RANS) equations and the two equation shear-stress transport k-w (SST k-w) model has been employed to numerically simulate the cold flow field in a special-shaped cavity-based supersonic combustor. In a cross-section shaped rectangular, hypersonic inlet with airflow at Mach 2.0 chamber, shock structures and flow characteristics of a herringbone-shaped boss and a herringbone-shaped cavity models were discussed, respectively. The results indicate: Firstly, according to the similarities of bevel-cutting shock characteristics between the boss case and the cavity case, the boss structure can serve as an ideal alternative model for shear-layer. Secondly, the eddies within cavity are composed of herringbone-spanwise vortexes, columnar vortices in the front and main-spanwise vortexes in the rear, featuring tilting, twisting and stretching. Thirdly, the simulated bottom-flow of cavity is in good agreement with experimental result, while the reverse flow-entrainment resulting from herringbone geometry and pressure gradient. However, the herringbone-shaped cavity has a better performance in fuel-mixing.


2000 ◽  
Vol 31 (1) ◽  
pp. 57-72 ◽  
Author(s):  
N. R. B. Olsen ◽  
D. K. Lysne

A three-dimensional numerical model was used to model water circulation and spatial variation of temperature in Lake Sperillen in Norway. A winter situation was simulated, with thermal stratification and ice cover. The numerical model solved the Navier-Stokes equations on a 3D unstructured non-orthogonal grid with hexahedral cells. The SIMPLE method was used for the pressure coupling and the k-ε model was used to model turbulence, with a modification for density stratification due to the vertical temperature profile. The results were compared with field measurements of the temperature in the lake, indicating the location of the water current. Reasonably good agreement was found.


2019 ◽  
Vol 9 (22) ◽  
pp. 4919 ◽  
Author(s):  
Wei Zhong ◽  
Tong Guang Wang ◽  
Wei Jun Zhu ◽  
Wen Zhong Shen

The Actuator Disc/Navier-Stokes (AD/NS) method has played a significant role in wind farm simulations. It is based on the assumption that the flow is azimuthally uniform in the rotor plane, and thus, requires a tip loss correction to take into account the effect of a finite number of blades. All existing tip loss corrections were originally proposed for the Blade-Element Momentum Theory (BEMT), and their implementations have to be changed when transplanted into the AD/NS method. The special focus of the present study is to investigate the performance of tip loss corrections combined in the AD/NS method. The study is conducted by using an axisymmetric AD/NS solver to simulate the flow past the experimental NREL Phase Ⅵ wind turbine and the virtual NREL 5MW wind turbine. Three different implementations of the widely used Glauert tip loss function F are discussed and evaluated. In addition, a newly developed tip loss correction is applied and compared with the above implementations. For both the small and large rotors under investigation, the three different implementations show a certain degree of difference to each other, although the relative difference in blade loads is generally no more than 4%. Their performance is roughly consistent with the standard Glauert correction employed in the BEMT, but they all tend to make the blade tip loads over-predicted. As an alternative method, the new tip loss correction shows superior performance in various flow conditions. A further investigation into the flow around and behind the rotors indicates that tip loss correction has a significant influence on the velocity development in the wake.


2012 ◽  
Vol 23 (04) ◽  
pp. 1250030 ◽  
Author(s):  
FAYÇAL HAMMAMI ◽  
NADER BEN-CHEIKH ◽  
ANTONIO CAMPO ◽  
BRAHIM BEN-BEYA ◽  
TAIEB LILI

In this work, a numerical study devoted to the two-dimensional and three-dimensional flow of a viscous, incompressible fluid inside a lid-driven cavity is undertaking. All transport equations are solved using the finite volume formulation on a staggered grid system and multi-grid acceleration. Quantitative aspects of two and three-dimensional flows in a lid-driven cavity for Reynolds number Re = 1000 show good agreement with benchmark results. An analysis of the flow evolution demonstrates that, with increments in Re beyond a certain critical value Rec, the steady flow becomes unstable and bifurcates into unsteady flow. It is observed that the transition from steadiness to unsteadiness follows the classical Hopf bifurcation. The time-dependent velocity distribution is studied in detail and the critical Reynolds number is localized for both 2D and 3D cases. Benchmark solutions for 2D and 3D lid-driven cavity flows are performed for Re = 1500 and 6000.


Sign in / Sign up

Export Citation Format

Share Document