Transition and Turbulence Modeling for the Prediction of Cavitating Tip Vortices

2020 ◽  
Vol 143 (1) ◽  
Author(s):  
Rens Liebrand ◽  
Maarten Klapwijk ◽  
Thomas Lloyd ◽  
Guilherme Vaz

Abstract This study evaluates the influence of transition and turbulence modeling on the prediction of wetted and cavitating tip vortices for an elliptical wing, while investigating the numerical errors. Transition modeling increases the quality of numerical predictions since the assumption of a fully turbulent boundary layer, commonly found in literature, contributes to underprediction of the tip vortex cavity size. By applying the local correlation-based transition model (LCTM) and controlling the boundary layer thickness using different turbulent inflow conditions, the pressure in the vortex was found to reduce by 20% for an Angle of Attack (AoA) of 5 deg. The consequent increase in cavity size was found to be of a similar order of magnitude. At 9 deg AoA, transition always occurs just downstream of the leading edge, making this AoA more suitable to investigate the effect of different turbulence modeling approaches. Azimuthal and axial velocity fields are validated against stereographic-particle image velocimetry (S-PIV) measurements. The time-averaged velocity profiles predicted by delayed detached-eddy simulation (DDES) and improved delayed detached-eddy simulation (IDDES) are close to the experiments; however, no velocity fluctuations and vortex dynamics are observed around the vortex. A comparison of wetted and cavitating simulations shows that the cavity leads to a change in the balance between production and destruction of turbulence kinetic energy, which reduces the turbulent diffusion in and around the cavity compared to wetted flow conditions. Consequently, the vapor flow exhibits the characteristics of a potential flow. Whether this is physically plausible remains to be investigated.

Author(s):  
Purvic Patel ◽  
Yunchao Yang ◽  
Gecheng Zha

Abstract This paper utilizes the Improved Delayed Detached Eddy Simulation (IDDES) to investigate the non-synchronous vibration (NSV) mechanism of a 1.5 stage high-speed axial compressor. The NSV occurs at a part speed in the rig test. A low diffusion E-CUSP approximate Riemann solver with a third order Weighted Essentially Non-Oscillating (WENO) scheme for the inviscid flux and a second order central differencing scheme for the viscous flux are employed to solve the 3D time accurate Navier-Stokes equations. The fully conservative sliding boundary condition is used to preserve the wake-propagation. The aerodynamic instability in the tip region induces two alternating low pressure regions near the leading and the trailing edge on the suction side of the rotor blade. It is observed that the circumferential tip vortex motion in the rotor passage above 75 % span and its coupling forces cause NSV at the operating speed. This instability moves in the counter-rotating direction in the rotational frame. The NSV results using URANS simulation is also presented for comparison. The predicted frequency with the IDDES and URANS using rigid blades agrees well with the measured frequency in the rig test. In addition to the NSV, the IDDES solver also captures the dominant engine order frequencies. The tip flow structures show the vortex filament with one end on the suction side of the rotor blade and other side terminating on the casing or the pressure side of the rotor blade.


2021 ◽  
pp. 1-10
Author(s):  
Zifei Yin ◽  
Paul Durbin

Abstract The adaptive, l2-omega delayed detached eddy simulation model was selected to simulate the flow in the V103 linear compressor cascade. The Reynolds number based on axial chord length is 138,500. Varies inflow turbulent intensities from 0% to 10% were tested to evaluate the performance of the adaptive model. The adaptive model is capable of capturing the laminar boundary layer and the large scale perturbations inside it. The instability of large scale disturbances signals the switch to a hybrid simulation of turbulent boundary layer -- the transition front is thus predicted. In the case of separation-induced transition, the adaptive model, which uses eddy simulation in separated flow, can predict the separation bubble size accurately. Generally, the adaptive, delayed detached eddy simulation model can simulate the transitional separated flow in a linear compressor cascade, with a correct response to varying turbulent intensities.


2009 ◽  
Vol 131 (11) ◽  
Author(s):  
Daniel C. Lyons ◽  
Leonard J. Peltier ◽  
Frank J. Zajaczkowski ◽  
Eric G. Paterson

Separated flow past a hump in a turbulent boundary layer is studied numerically using detached-eddy simulation (DES), zonal detached-eddy simulation (ZDES), delayed detached-eddy simulation (DDES), and Reynolds-averaged Navier–Stokes (RANS) modeling. The geometry is smooth so the separation point is a function of the flow solution. Comparisons to experimental data show that RANS with the Spalart–Allmaras turbulence model predicts the mean-field statistics well. The ZDES and DDES methods perform better than the DES formulation and are comparable to RANS in most statistics. Analyses motivate that modeled-stress depletion near the separation point contributes to differences observed in the DES variants. The order of accuracy of the flow solver ACUSOLVE is also documented.


2017 ◽  
Vol 34 (2) ◽  
pp. 231-241
Author(s):  
M. Zhang ◽  
Y. K. Wang ◽  
S. Fu

AbstractThe formation and evolution of wingtip vortex system generated from three wing configurations are simulated with the improved delayed detached eddy simulation (IDDES) method. Numerical results show that each layout produces an interacting wingtip vortex system. These three corresponding vortical interactions are, respectively, the interaction between wingtip vortex and its counter-rotating vortex, winglet-tip vortex, and winglet four-vortex system. The fluid entrainment of ambient fluid and vortical impulse transport resulted from inductive effect have been founded generally existing in its formation and evolution. These two dominated mechanisms account for induced drag generation. On one hand, the winglet with toed-out angle is considered capable of changing the flow field around the winglet, and decomposing the winglet-tip vortex into four small vortices. Due to quite few fluid entrainment effects, this typical four-vortex system that cannot merge and only dissipate in the near wake scarcely contributes to the induced drag. On the other hand, a potential drag reduction method is also indicated that a lower induced drag can be obtained when the merger of wingtip and winglet-tip vortex is controlled and eliminated. This investigation will offer a novel perspective to guide the design of wingtip device and method of crusing resistance reduction for aircrafts.


Author(s):  
Daniel C. Lyons ◽  
Leonard J. Peltier ◽  
Frank J. Zajaczkowski ◽  
Eric G. Paterson

Turbulent flow past the Glauert-Goldschmied body, a flow-control hump in a turbulent boundary layer, is studied numerically using detached-eddy simulation (DES), zonal detached-eddy simulation (ZDES), delayed detached-eddy simulation (DDES), and Reynolds-Averaged Navier-Stokes (RANS) modeling. The geometry is smooth so the downstream separation point is not set by facets of the geometry but is a function of the pressure gradient, a challenging condition for turbulence models. Comparisons to experimental data show that RANS with the Spalart-Allmaras turbulence model predicts the mean-field statistics well. The ZDES and DDES methods perform better than the DES formulation and are comparable to RANS in most statistics. An analysis of model behavior indicates that modeled stress depletion in the detached shear layer shortly after separation leads to loss of accuracy in the DES variants.


2017 ◽  
Vol 34 (2) ◽  
pp. 123-134 ◽  
Author(s):  
L. Zhang ◽  
J. Li ◽  
Y. F. Mou ◽  
H. Zhang ◽  
W. B. Shi ◽  
...  

AbstractAccurate prediction of the flow around multi-element airfoil is a prerequisite for improving aerodynamic performance, but its complex flow features impose high demands on turbulence modeling. In this work, delayed detached-eddy-simulation (DDES) and zonal detached-eddy-simulation (ZDES) was applied to simulate the flow past a three-element airfoil. To investigate the effects of numerical dissipation of spatial schemes, the third-order MUSCL and the fifth-order interpolation based on modified Roe scheme were implemented. From the comparisons between the calculations and the available experimental result, third-order MUSCL-Roe can provide satisfactory mean velocity profiles, but the excessive dissipation suppresses the velocity fluctuations level and eliminates the small-scale structures; DDES cannot reproduce the separation near the trailing edge of the flap which lead to the discrepancy in mean pressure coefficients, while ZDES result has better tally with the experiment.


Author(s):  
Yong-chen Pan ◽  
Jian-wei Yao ◽  
Rui Xu ◽  
Tao Liu ◽  
Jun Zheng ◽  
...  

Two flow cases for high speed train models with different lengths have been numerically computed by performing the improved delayed detached-eddy simulation. Based on the Omega method and turbulence production (TP) distribution, the relations between the shear flow and vortices in the near turbulent wake of a high speed train have been comparatively analyzed. First, in the wake region immediately close to the tail, the boundary layer separation plays significant roles. And the mechanism makes shear deformation prominent in the region with the formed vortices. Moreover, the shear layers are pertinent to the boundary-layer thicknesses and help to the TP distribution. However, the shear-dominated region is very limited due to high dissipation. One the other hand, a vast majority of the vortices captured with the parameter Omega increasing in the downstream region away from the tail. And the TP distributions are stable at different streamwise positions, though obviously decreased. They are greatly attributed to the mean strain rate in the horizontal plane. Meanwhile, the vortical vorticity is thought to be the dominant component inside the vortex cores, although the shear becomes weaker. And thus the turbulence itself can be spatially sustained due to the relatively stable vortex structure. Moreover, the weak shear is believed to depend upon the interaction between the vortices and the ground.


Sign in / Sign up

Export Citation Format

Share Document