Large Eddy Simulation of the Flow Past Pitching NACA0012 Airfoil at 1E5 Reynolds Number

Author(s):  
Venkata Ravishankar Kasibhotla ◽  
Danesh Tafti

The paper is concerned with the prediction and analysis of dynamic stall of flow past pitching NACA-0012 airfoil at 105 Reynolds number based on the chord length of the airfoil and at reduced frequency of 0.188 in a three dimensional flow field. The turbulence in the flow field is resolved using large eddy simulations with dynamic Smagorinsky model at the sub grid scale. The lift hysteresis plots indicate closer match to experimental results, although discrepancies exist during the downstroke. The development of dynamic stall vortex, vortex shedding and reattachment as predicted by the present study are discussed in detail. This study has shown that the downstroke phase of the pitching motion is strongly three dimensional and is highly complex, whereas the flow is practically two dimensional during the upstroke.

Author(s):  
Venkata Ravishankar Kasibhotla ◽  
Danesh Tafti

The paper is concerned with the prediction and analysis of dynamic stall of flow past a pitching NACA0012 airfoil at 1 million Reynolds number based on the chord length of the airfoil and at reduced frequency of 0.25 in a three dimensional flow field. The turbulence in the flow field is resolved using large eddy simulations with the dynamic Smagorinsky model at the sub grid scale. The development of dynamic stall vortex, shedding and reattachment as predicted by the present study are discussed in detail. This study has shown that the downstroke phase of the pitching motion is strongly three dimensional and is highly complex, whereas the flow is practically two dimensional during the upstroke. The lift coefficient agrees well with the measurements during the upstroke. However, there are differences during the downstroke. The computed lift coefficient undergoes a sharp drop during the start of the downstroke as the convected leading edge vortex moves away from the airfoil surface. This is followed by a recovery of the lift coefficient with the formation of a secondary trailing edge vortex. While these dynamics are clearly reflected in the predicted lift coefficient, the experimental evolution of lift during the downstroke maintains a fairly smooth and monotonic decrease in the lift coefficient with no lift recovery. The simulations also show that the reattachment process of the stalled airfoil is completed before the start of the upstroke in the subsequent cycle due to the high reduced frequency of the pitching cycle.


2021 ◽  
Vol 158 (A1) ◽  
Author(s):  
S Kim ◽  
P A Wilson ◽  
Z Chen

The effect of the spanwise discretisation on numerical calculations of the turbulent flow around a circular cylinder is systematically assessed at a subcritical Reynolds number of 10000 in the frame of three-dimensional large-eddy simulation. The eddy-viscosity k-equation subgrid scale model is implemented to evaluate unsteady turbulent flow field. Large-eddy simulation is known to be a reliable method to resolve such a challenging flow field, however, the high computational efforts restrict to low Reynolds number flow or two-dimensional calculations. Therefore, minimum spatial density in the spanwise direction or cylinder axis direction needs to be carefully evaluated in order to reduce high computational resources. In the present study, the influence of the spanwise resolutions to satisfactorily represent three- dimensional complex flow features is discussed in detail and minimum spatial density for high Reynolds flow is suggested.


2012 ◽  
Vol 232 ◽  
pp. 471-476 ◽  
Author(s):  
Rui Zhao ◽  
Chao Yan

The flow past a circular cylinder at a subcritical Reynolds number 3900 was simulated by the method of detached-eddy simulation (DES). The objective of this present work is not to investigate the physical phenomena of the flow but to study modeling as well as numerical aspects which influence the quality of DES solutions in detail. Firstly, four typical spanwise lengths are chosen and the results are systematically compared. The trend of DES results along the span increment is different from previous large-eddy simulation (LES) investigation. A wider spanwise length does not necessary improve the results. Then, the influence of mesh resolution is studied and found that both too coarse and over refined grids will deteriorate the performance of DES. Finally, different orders of numerical schemes are applied in the inviscid fluxes and the viscous terms. The discrepancies among different schemes are found tiny. However, the instantaneous flow structures produced by 5th order WENO with 4th order central differencing scheme are more abundant than the others. That is, for the time-averaged quantities, the second-order accurate schemes are effective enough, whereas the higher-order accurate methods are needed to resolve the transient characteristics of the flow.


2015 ◽  
Vol 2015 ◽  
pp. 1-7
Author(s):  
Yan Xu ◽  
Zunce Wang ◽  
Lin Ke ◽  
Sen Li ◽  
Jinglong Zhang

Reynolds Stress Model and Large Eddy Simulation are used to respectively perform numerical simulation for the flow field of a hydrocyclone. The three-dimensional hexahedral computational grids were generated. Turbulence intensity, vorticity, and the velocity distribution of different cross sections were gained. The velocity simulation results were compared with the LDV test results, and the results indicated that Large Eddy Simulation was more close to LDV experimental data. Large Eddy Simulation was a relatively appropriate method for simulation of flow field within a hydrocyclone.


2014 ◽  
Vol 1008-1009 ◽  
pp. 906-909
Author(s):  
Xiao Bing Wang

Unsteady flow of a single six-blade ruston turbine in stirred tank is numerically simulated by using the large eddy simulation. Then the effect of the turbine installation position on mixing flow field is studied. The result shows that with a relatively low paddle installation position, mixing effect at the bottom of tank is obvious, while which go against the materials at the top layer mixing. When the paddle is installed at the top of the stirred tank, liquid splash and a concave downward liquid surface are easily caused. Finally the cavitation phenomenon is generated. When the paddle is installed from 1/3H to 2/3H, there are a uniform flow field distribution and higher average velocity flow. Large size vortex structures at the top and the bottom of the paddle are obvious which is beneficial to mix the materials.


Author(s):  
Mojtaba Honarmand ◽  
Mohammad Hassan Djavareshkian ◽  
Behzad Forouzi Feshalami ◽  
Esmaeil Esmaeilifar

In this research, viscous, unsteady and turbulent fluid flow is simulated numerically around a pitching NACA0012 airfoil in the dynamic stall area. The Navier-Stokes equations are discretized based on the finite volume method and are solved by the PIMPLE algorithm in the open source software, namely OpenFOAM. The SST k - ω model is used as the turbulence model for Low Reynolds Number flows in the order of 105. A homogenous dynamic mesh is used to reduce cell skewness of mesh to prevent non-physical oscillations in aerodynamic forces unlike previous studies. In this paper, the effects of Reynolds number, reduced frequency, oscillation amplitude and airfoil thickness on aerodynamic force coefficients and dynamic stall delay are investigated. These parameters have a significant impact on the maximum lift, drag, the ratio of aerodynamic forces and the location of dynamic stall. The most important parameters that affect the maximum lift to drag coefficient ratio and cause dynamic stall delaying are airfoil thickness and reduced frequency, respectively.


Sign in / Sign up

Export Citation Format

Share Document