Prediction of separation flows around a 6:1 prolate spheroid using RANS/LES hybrid approaches

This paper presents hybrid Reynolds-averaged Navier-Stokes （RANS） and large-eddy-simulation （LES） methods for the separated flows at high angles of attack around a 6：1 prolate spheroid. The RANS/LES hybrid meth- ods studied in this work include the detached eddy simulation （DES） based on Spalart-Allmaras （S-A）, Menter＇s k-ω shear-stress-transport （SST） and k-o9 with weakly nonlinear eddy viscosity formulation （Wilcox-Durbin＋, WD＋） models and the zonalANS/LES methods based on the SST and WD＋ models. The switch from RANS near the wall to LES in the core flow region is smooth through the implementation of a flow-dependent blending function for the zonal hybrid method. All the hybrid methods are designed to have a RANS mode for the attached flows and have a LES behavior for the separated flows. The main objective of this paper is to apply the hybrid methods for the high Reynolds number separated flows around prolate spheroid at high-incidences. A fourth-order central scheme with fourth-order artificial viscosity is applied for spatial differencing. The fully implicit lower-upper symmetric-Gauss-Seidel with pseudo time sub-iteration is taken as the temporal differentiation. Comparisons with available measurements are carried out for pressure distribution, skin friction, and profiles of velocity, etc. Reasonable agreement with the experiments, accounting for the effect on grids and fundamental turbulence models, is obtained for the separation flows.

基于RANS/LES混合方法的近岸水域波浪增阻影响的数值模拟

分析船舶在近岸水域航行时所受波浪增阻大小的影响因素,以寻求降低波浪增阻的方法。基于近岸水域中船舶阻力理论,采用雷诺平均(RANS)和大涡模拟(LES)相结合的方法进行数值求解,并运用两相流方法和造波消波技术建立数值水池模型。用Fluent求解器对船舶在不同工况下的波浪增阻进行数值模拟。试验结果表明:用该方法对近岸水域航行的船舶进行波浪增阻的数值模拟可有效提高计算机运算效率,在粗糙网格层面上得到的计算结果精度良好;船舶在近岸水域中航行时的波浪增阻随波高的增加而增加,并且与船体湿润面积呈正相关;船舶在近岸水域以中低速航行时所受波浪增阻随弗劳德数的增加而增加。该研究为近岸水域通航船舶减小航行阻力提供了可行的方法与思路。

Stress-Blended Eddy Simulation of the Turbulent Flow Around Circular Cylinder with Dynamic Hybrid RANS/LES Model

In order to verify the effectiveness and superiority of the dynamic hybrid RANS/LES(DHRL)model,the flow around a cylinder with sinusoidal fluctuating velocity at the inlet was used as the test case.The latest computational fluid dynamics(CFD)model can flexibly choose any existing large-eddy simulation(LES)method combined with RANS method to calculate the flow field.In addition,the DLES model and DDES model are selected as typical representatives of the turbulence model to compare the capture ability of the flow field mechanism.The internal flow field including the y+value,velocity distribution,turbulent kinetic energy and vortex structures is comprehensively analyzed.Finally,the results show that the new model has enough sensitivity to capture the information of the flow field and has more consistent velocity distribution with the experimental value,which shows its potential in practical engineering applications to some extent.

Numerical study on turbulent mixed convection in a vertical plane channel using hybrid RANS/LES and LES models

Two Delayed-Detached Eddy Simulation(DDES) models, and a Large-Eddy Simulation(LES) model are used to investigate the turbulent flows and mixed convection between a hot plate and a cold plate via the software FLUENT. The two DDES models include Production-limited DDES(PL-DDES) and Improved DDES(IDDES) models.The Wall-Adapting Local Eddy-Viscosity(WALE) model is the used LES model. The numerical computations are performed at Reynolds number Reb= 4494 and different Richardson numbers Ri = 0.025, 0.048, 0.1. The comparing data is from the Direct Numerical Simulation(DNS) at Reb= 4494 and Ri = 0.048. The comparison reveals that the two DDES models have better performance in predicting time-averaged parameters than the WALE model in the aiding flow. The best predicted time-averaged results are obtained by the PL-DDES model in the opposing flow. Furthermore, the results of different Ri obtained by the PL-DDES model agree well with the DNS data.

Hybrid RANS/LES Simulations of Two and Three Dimensional Supersonic Cavity Flows

Supersonic unsteady flows of two and three dimensional cavities, which have the same basic measures, length, depth and aft wall angle, are investigated numerically by using hybrid RANS/LES （Reynolds-Averaged Navier-Stokes/Large-Eddy Simulation） method. In particular, the mass exchange and oscillation characteristics of the cavities are analyzed and compared. It is shown that the two and three dimensional cavities have almost the same residence time except that the three dimensional one has a little larger mass decay at beginning, which may attribute to the influence of the streamwise vortices. The two dimensional cavity has three dominant frequencies while the three dimensional one has only one oscillation dominant frequency, that is, the three dimensional effects simplify the oscillation modes in the cavity. The distributions of oscillation energy are approximately universal in the transverse and spanwise directions. However, the oscillation energy in the streamwise direction shows a hybrid monotone/periodic distribution, which may be caused by the streamwise-propagating pressure waves.

Hybrid LES/RANS Simulations of Compressible Flow in a Linear Cascade of Flat Blade Profiles

The paper reports on 3D numerical simulations of unsteady compressible airflow in a blade cascade consisting of flat profiles using a hybrid LES/RANS approach including a transition model.As a first step towards simulation of blade flutter in turbomachinery,various incidence angle offsets of the middle blade were modeled.All simulations were run for the flow regime characterized by outlet isentropic Mach number M_(is)=0.5and zero incidence.The results of the LES/RANS simulations(pressure and Mach number distributions)were compared to a baseline RANS model,and to experimental data measured in a high-speed wind tunnel.The numerical results show that both methods overpredict flow separation taking place at the leading edge.In this regard,the hybrid LES/RANS method does not provide superior results compared to the traditional RANS simulations.Nevertheless,the LES/RANS results also capture vortex shedding from the blunt trailing edge.The frequency of the trailing edge vortex shedding in CFD simulations matches perfectly the spectral peak recorded during wind tunnel measurements.

Interaction mechanisms of shock waves with the boundary layer and wakes in a highly-loaded NGV using hybrid RANS/LES

Accurate predictions of Shock Waves and Boundary Layer Interaction(SWBLI)and strong Shock Waves and Wake Vortices Interaction(SWWVI)in a highly-loaded turbine propose challenges to the currently widely used Reynolds-Averaged Navier-Stokes(RANS)model.In this work,the SWBLI and the SWWVI in a highly-loaded Nozzle Guide Vane(NGV)are studied using a hybrid RANS/LES strategy.The Turbulence Kinetic Energy(TKE)budget and the Proper Orthogonal Decomposition(POD)method are used to analyze flow mechanisms.Results show that this hybrid RANS/LES method can obtain detailed flow structures for flow mechanisms analysis.Strong shock waves induce boundary layer separation,while the presence of a separation bubble can in turn lead to a Mach reflection phenomenon.The shock waves cause trailing-edge vortices to break clearly,and the wakes,in turn,can change the shocks intensity and direction.Furthermore,the Entropy Generation Rate(EGR)is used to analyze the irreversible loss.It turns out that the SWWVI can reduce the flow field loss.There are several weak shock waves in the NGV flow field,which can increase the irreversible loss.This work offers flow mechanisms analysis and presents the EGR distribution in SWBLI and SWWVI areas in a transonic turbine blade.

A hybrid RANS/LES model for simulating time-dependent cloud cavitating flow around a NACA66 hydrofoil

Cloud cavitating flow is highly turbulent and dominated by coherent large-scale anisotropic vortical structures. For the numer- ical investigation of such a class of flow, large eddy simulation （LES） is a reliable method but it is computationally extremely costly in engineering applications. An efficient approach to reduce the computational cost is to combine Reynolds-averaged Navier-Stokes （RANS） equations with LES used only in the parts of interest, such as massively separated flow regions. A new hybrid RANS/LES model, the modified filter-based method （FBM）, is proposed in the present study which can perform RANS or LES depending on the numerical resolution. Compared to the original FBM, the new method has three modifications： the state-of-the-art shear stress transport （SST） model replaces the k-c model as a baseline RANS model. A shielding function is introduced to obviate the switch from RANS to LES occurring inside the boundary layer. An appropriate threshold controlling the switch from RANS to LES is added to achieve an optimal predictive accuracy. The new model is assessed for its predictive capability of highly unsteady cavitating flows in a typical case of cloud cavitation around a NACA66 hydrofoil. The new mod- el results are compared with data obtained from the Smagorinsky LES and SST model based on the same homogeneous Zwart cavitation model. It is found that the modified FBM method has significant advantages over SST model in all aspects of pre- dicted instantaneous and mean flow field, and its predictive accuracy is comparable to the Smagorinsky LES model even using a much coarser grid in the simulations.

RANS and hybrid LES/RANS simulations of flow over a square cylinder

Unsteady RANS(URANS),hybrid LES/RANS and IDDES simulations were conducted to numerically investigate the velocity field around,and pressures distribution and forces over a square cylinder immersed in a uniform,steady oncoming flow with Reynolds number Re=21,400.The vortex shedding responses in terms of Strouhal number,the pressure distribution,the velocity profile and the velocity fluctuations obtained by numerical simulations are compared with experimental data.Compared with 2D URANS simulation,3D simulations using hybrid LES/RANS and IDDES models provide more accurate prediction on the responses in the wake,including mean streamwise velocity profile and rms velocity fluctuations.This also results in more accurate prediction of time-averaged surface pressure coefficient on the rear surface obtained by 3D hybrid LES/RANS and IDDES simulations than by URANS simulation.When a hybrid LES/RANS model or IDDES model is used,a more accurate prediction for either pressure coefficient or velocity profile(especially in the far wake region)is not guaranteed by increasing the mesh resolution along the spanwise direction of the square cylinder.

Application of the hybrid RANS/LES method on the hydraulic dynamic performance of centrifugal pumps

The nonlinear hybrid RANS/LES method has been developed considering the simulation efficiency and accuracy to simulate the hydraulic dynamic characteristics of the centrifugal pump. It has been proved that the nonlinear hybrid RANS/LES method could effectively capture the unsteady flow structure and pressure pulsation in 2?D centrifugal pumps. To further investigate the hydraulic dynamic performance of a redesigned 3-D-gap drainage centrifugal pump, the nonlinear hybrid RANS/LES method is employed. Both numerical simulation and experimental results show that the hydraulic performance and pressure pulsation characteristics of the 3-D-gap drainage impeller centrifugal pump are significantly enhanced.

An hybrid RANS/LES model for simulation of complex turbulent flow

A nonqinear eddy viscosity model （NLEVM） and a scalable hybrid Reynolds averaged Navier-Stokes/large eddy simula- tion （RANS/LES） strategy are developed to improve the capability of the eddy viscosity model （EVM） to simulate complex flows featuring separations and unsteady motions. To study the performance of the NLEVM, numerical simulations around S809 airfoil are carried out and the results show that the NLEVM performs much better when a large separation occurs. Calculated results of the flow around a triangular cylinder show that the NLEVM can improve the precision of the flow fields to some extents, but the error is still considerable, and the small turbulence structures can not be clearly captured as the EVM. Whereas the scalable hybrid RANS/LES model based on the NLEVM is fairy effective on resolving the turbulent structures and can give more satisfactory predictions of the flow fields.

Improved delayed detached-eddy simulation of massive separation around triple cylinders

The massively separated flow past triple cylin- ders （TriC） in tandem arrangement is simulated using the improved delayed detached-eddy simulation （IDDES） method based on the shear stress transport （SST） model, coupled with the high order adaptive dissipation scheme. The spacing between adjacent cylinders is sub-critical （1.435D）. IDDES prediction of two cylinders （TC） with the same spacing is compared to experimental data for validation, and the numerical results agree well with the available measurements, except for the asymmetry in the gap region. The flow past TriC is investigated using the same method. Generally, the mean flow quantities past TriC, such as the velocity, pressure, and vorticity, are similar to the corresponding components of TC. However, the pressure fluctuations on the TriC surface are uniformly larger than those on TC. Meanwhile, the instantaneous flows past TriC are much more complex. The periodical blockage in the first gap region is found in the TriC case and leads to the up-and-down movement of shear layer in the second gap region.

Hybrid Reynolds-averaged Navier-Stokes/large-eddy simulation of jet mixing in a supersonic crossflow

A sonic under-expanded transverse jet injection into a Ma 1.6 supersonic crossflow is investigated numerically using our hybrid RANS/LES (Reynolds-averaged Navier-Stokes/large eddy simulation) method. First, a calculation is carried out to validate the code, where both the instantaneous and statistical results show good agreement with the existing experimental data. Then the jet-mixing characteristics are analyzed. It is observed that the large-scale vortex on the windward portion of the jet boundary is formed mainly by the intermittent impingement of the incoming high-speed fluid on the relatively low-speed region of the upstream jet boundary, where the interaction between the upstream separated region and the jet supplies a favorable pressure condition for the sustaining acceleration of the high-speed fluid during the vortex forming, associated with which the incoming fluid is entrained into the jet boundary and large-scale mixing occurs. Meanwhile, the secondary recirculation zone between the upstream separated region and the jet is observed to develop evidently during the vortex forming, inducing the entrainment of jet fluid into the upstream separated region. Moreover, effects of the incoming boundary layer on the jet mixing are addressed.

Towards an entropy-based detached-eddy simulation

A concept of entropy increment ratio （ s- ） is introduced for compressible turbulence simulation through a series of direct nu- merical simulations （DNS）. s- represents the dissipation rate per unit mechanical energy with the benefit of independence of freestream Mach numbers. Based on this feature, we construct the shielding function f, to describe the boundary layer region and propose an entropy-based detached-eddy simulation method （SDES）. This approach follows the spirit of delayed de- tached-eddy simulation （DDES） proposed by Spalart et al. in 2005, but it exhibits much better behavior after their performanc- es are compared in the following flows, namely, pure attached flow with thick boundary layer （a supersonic fiat-plate flow with high Reynolds number）, fully separated flow （the supersonic base flow）, and separated-reattached flow （the supersonic cavity-ramp flow）. The Reynolds-averaged Navier-Stokes （RANS） resolved region is reliably preserved and the modeled stress depletion （MSD） phenomenon which is inherent in DES and DDES is partly alleviated. Moreover, this new hybrid strategy is simple and general, making it applicable to other models related to the boundary layer predictions.

A rescaling method for correcting log-layer mismatch in detached eddy simulation

In order to correct the unphysical log-layer mismatch commonly encountered in detached eddy simulation(DES) of flows with attached boundary layers,a function M,ML,which has a multi-layer structure with scaling laws in each layer and a plateau related to the Kármán constant,is defined.The height of this plateau is found to be crucial for obtaining the correct log-layer.A target scaling function is designed which equals M,ML in the near-wall region,but with the height of plateau determined analytically from the classical log-law.This scaling function is used as a target function according to which the resolved turbulent fluctuations are renormalized,in order to recover the height of plateau prescribed by the log-law.The renormalization procedure guarantees the height of M,ML required by log-law,resulting in correct log layer slope.The method is also shown to maintain similar turbulent properties in the large eddy simulation(LES) region of DES method.Hence it predicts the turbulent intensity correctly.The results demonstrate the relationship between constant M,ML and log-law profile of mean velocity,and relate the Kármán constant to turbulent fluctuations,implying a complete description of turbulent structural ensemble dynamics.The proposed method can be extended to more general flows with log layers since it uses only the log-law with Kármán constant as the input,while the intercept of log layer depends on the solution of Spalart-Allmaras(SA) model in the near-wall field,where Reynolds-averaged Navier-Stokes(RANS) solutions are accurate.