Detached-eddy Simulation for Time-dependent Turbulent Cavitating Flows

The Reynolds-averaged Navier-Stokes（RANS）,such as the original k-ω two-equation closures,have been very popular in providing good prediction for a wide variety of flows with presently available computational resource.But for cavitating flows,the above equations noticeably over-predict turbulent production and hence effective viscosity.In this paper,the detached eddy simulation（DES） method for time-dependent turbulent cavitating flows is investigated.To assess the state-of-the-art of computational capabilities,different turbulence models including the widely used RANS model and DES model are conducted.Firstly,in order to investigate the grid dependency in computations,different grid sizes are adopted in the computation.Furthermore,the credibility of DES model is supported by the unsteady cavitating flows over a 2D hydrofoil.The results show that the DES model can effectively reduce the eddy viscosities.From the experimental validations regarding the force analysis,frequency and the unsteady cavity visualizations,more favorable agreement with experimental visualizations and measurements are obtained by DES model.DES model is better able to capture unsteady phenomena including cavity length and the resulting hydrodynamic characteristics,reproduces the time-averaged velocity quantitatively around the hydrofoil,and yields more acceptable and unsteady dynamics features.The DES model has shown to be effective in improving the overall predictive capability of unsteady cavitating flows.

Detached-eddy simulation of flow around high-speed train on a bridge under cross winds

In order to describe an investigation of the flow around high-speed train on a bridge under cross winds using detached-eddy simulation(DES), a 1/8th scale model of a three-car high-speed train and a typical bridge model are employed, Numerical wind tunnel technology based on computational fluid dynamics(CFD) is used, and the CFD models are set as stationary models. The Reynolds number of the flow, based on the inflow velocity and the height of the vehicle, is 1.9×10~6. The computations are conducted under three cases, train on the windward track on the bridge(WWC), train on the leeward track on the bridge(LWC) and train on the flat ground(FGC). Commercial software FLUENT is used and the mesh sensitivity research is carried out by three different grids: coarse, medium and fine. Results show that compared with FGC case, the side force coefficients of the head cars for the WWC and LWC cases increases by 14% and 29%, respectively; the coefficients of middle cars for the WWC and LWC increase by 32% and 10%, respectively; and that of the tail car increases by 45% for the WWC whereas decreases by 2% for the LWC case. The most notable thing is that the side force and the rolling moment of the head car are greater for the LWC, while the side force and the rolling moment of the middle car and the tail car are greater for the WWC. Comparing the velocity profiles at different locations, the flow is significantly influenced by the bridge-train system when the air is close to it. For the three cases(WWC, LWC and FGC), the pressure on the windward side of train is mostly positive while that of the leeward side is negative. The discrepancy of train's aerodynamic force is due to the different surface area of positive pressure and negative pressure zone. Many vortices are born on the leeward edge of the roofs. Theses vortices develop downstream, detach and dissipate into the wake region. The eddies develop irregularly, leading to a noticeably turbulent flow at leeward side of train.

Detached-eddy simulation of supersonic turbulent flow over rearward facing step

Detailed formulation of the detached-eddy simulation （DES） based on shear stress transport model （SST） was deduced. Based on finite volume method, detached-eddy simulation of super- sonic turbulent flow over rearward facing step was conducted. Multi-block structured non-conform grid was designed. The obtained distribution of wall pressure coefficient along the downstream of the step was in good agreement with experiment. The structure of supersonic flow over rearward facing step is complicated and contains boundary layer-shock interaction, shock-shock interaction and sec- ondary flow separation. The predictions of the corner pressure and location of the reattachment shock are more exact than the Reynolds-averaged method. Results show that detached-eddy simulation can simulate the supersonic separated flows exactly and has advantages over Reynolds-averaged method. The partition of the detached-eddy simulation can be easily achieved by using multi-block non-conform grid, which reduces the requirement of the computational resources as well.

LARGE-EDDY AND DETACHED-EDDY SIMULATIONS OF THE SEPARATED FLOW AROUND A CIRCULAR CYLINDER

The separated turbulent flow around a circular cylinder is investigated using Large-Eddy Simulation （LES）, Detached-Eddy Simulation （DES, or hybrid RANS/LES methods）, and Unsteady Reynolds-Averaged Navier-Stokes （URANS）. The purpose of this study is to examine some typical simulation approaches for the prediction of complex separated turbulent flow and to clarify the capability of applying these approaches to a typical case of the separated turbulent flow around a circular cylinder. Several turbulence models, i.e. dynamic Sub-grid Scale （SGS） model in LES, the DES-based Spalart-Allmaras （S-A） and κ-ω Shear-Stress- Transport （SST） models in DES, and the S-A and SST models in URANS, are used in the calculations. Some typical results, e.g., the mean pressure and drag coefficients, velocity profiles, Strouhal number, and Reynolds stresses, are obtained and compared with previous computational and experimental data. Based on our extensive calculations, we assess the capability and performance of these simulation approaches coupled with the relevant turbulence models to predict the separated turbulent flow.

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.

Detached-eddy simulation of supersonic flow past a spike-tipped blunt nose

Space vehicle in atmosphere travels mostly at supersonic speed and generates a very strong bow shockwave around its blunt nose. Oblique shock and conical separated flow zone generated by a forward disk-tip spike significantly reduce the drag by reducing the high pressure area on the blunt nose. This study employs improved delayed detached eddy simulation to investigate the characteristic flow structures around a spike-tipped blunt nose at Mach number of 3 and Reynolds number（based on the blunt-body diameter） of 2.72x10;. The calculated time-averaged quantities agree well with experimental data. Characteristic frequencies in different flow regions are extracted using fast Fourier transform. It is found that two distinct instability modes exist: oscillation mode and pulsation mode. The former is related to the foreshock/turbulence interaction with nondimensional frequency at around 0.004. The latter corresponds to the interaction between turbulence and shock structures around the blunt nose, with a typical coherent structure shedding frequency at 0.092.

薄翼失速翼型前缘分离泡对失速特性的影响

DES方法结合了RANS(Reynolds-averaged Navier-Stokes)和LES(Large Eddy Simulation approaches)的优点。在近壁面它体现为RANS模型的特点而在远离物面处又起到LES的亚格子模型的特性。论文应用DES(Detached-EddySimulation)方法讨论了影响薄翼失速的分离泡对翼型的升力特性影响。

HYBRID SCHEME FOR COMPRESSIBLE TURBULENT FLOW AROUND CURVED SURFACE BODY

A hybrid central-upwind scheme is proposed. Two sub-schemes, the central difference scheme and the Roets flux difference splitting scheme, are hybridized by means of a binary sensor function. In order to examine the capability of the proposed hybrid scheme in computing compressible turbulent flow around a curved surface body, especially the flow involving shock wave, three typical eases are investigated by using detached-eddy simulation technique. Numerical results show good agreements with the experimental measurements. The present hybrid scheme can be applied to simulating the compressible flow around a curved surface body involving shock wave and turbulence.

Applications of scale-adaptive simulation technique based on one-equation turbulence model

A modified scale-adaptive simulation （SAS） technique based on the Spalart- Allmaras （SA） model is proposed. To clarify its capability in prediction of the complex turbulent flow, two typical cases are carried out, i.e., the subcritical flow past a circular cylinder and the transonic flow over a hemisphere cylinder. For comparison, the same cases are calculated by the detached-eddy simulation （DES）, the delayed-detached eddy simulation （DDES）, and the XY-SAS approaches. Some typical results including the mean pressure coefficient, velocity, and Reynolds stress profiles are obtained and compared with the experiments. Extensive calculations show that the proposed SAS technique can give better prediction of the massively separated flow and shock/turbulent-boundary-layer interaction than the DES and DDES methods. Furthermore, by the comparison of the XY-SAS model with the present SAS model, some improvements can be obtained.

Hybrid RANS-LES of Shaped Hole Film Cooling on an Adiabatic Flat Plate at Low Reynolds Number

Hybrid RANS-LES methods offer a means of reducing computational cost and setup time to simulate transitional flows. Several methods are evaluated in ANSYS CFX, including Scale-Adaptive Simulation (SAS), Shielded Detached Eddy Simulation (SDES), Stress-Blended Eddy Simulation (SBES), and Zonal Large Eddy Simulation (ZLES), along with a no-model laminar simulation. Each is used to simulate an adiabatic flat plate film cooling experiment of a shaped hole at low Reynolds number. Adiabatic effectiveness is calculated for Blowing Ratio (BR) = 1.5 and Density Ratio (DR) = 1.5. The ZLES method and laminar simulation most accurately match experimental lateral-average adiabatic effectiveness along the streamwise direction from the trailing edge of the hole to 35 hole diameters downstream of the hole (X/D = 0 to X/D = 35), with RMS deviations of 5.1% and 4.2%, and maximum deviations of 8% and 11%, respectively. The accuracy of these models is attributed to the resolution of turbulent structures in not only the mixing region but in the upstream boundary layer as well, where the other methods utilize RANS and do not switch to LES.

Constrained large-eddy simulation and detached eddy simulation of flow past a commercial aircraft at 14 degrees angle of attack

With the development of computational power and numerical algorithms,computational fluid dynamics(CFD) has become an important strategy for the design of aircraft,which significantly reduces the reliance on wind-tunnel and flight tests.In this paper,we conducted a numerical investigation on the flow past a full commercial aircraft at Mach number 0.2 and 14 degrees angle of attack by means of Reynolds-averaged Navier-Stokes(RANS),detached-eddy simulation(DES) and our newly developed constrained large-eddy simulation(CLES).The objective of this paper is to study the capability of these models in simulating turbulent flows.To our knowledge,this is the first large-eddy simulation method for full commercial aircraft simulation.The results show that the CLES can predict the mean statistical quantities well,qualitatively consistent with traditional methods,and can capture more small-scale structures near the surface of the aircraft with massive separations.Our study demonstrates that CLES is a promising alternative for simulating real engineering turbulent flows.

Vortex structures of dynamic pure yaw test using DDES approach and vortex identification method

Considered as the building blocks,vortex structures with variety of sizes and intensity are widely recognized in the viscous flow field around ship.In this paper,the computational fluid dynamics(CFD)solver,naoe-FOAM-SJTU,coupled with delayed detached-eddy simulation(DDES)is adopted to analyze the vortex structures around the benchmark model Yupeng Ship in dynamic pure yaw tests,which are captured by third generation of vortex identification method.The good agreement of the predicted force/moment by DDES method with the experimental data indicates that the present numerical schemes are reliable and robust.Three vortex identification methods,Q-criteria,Ω_(R) and Liutex,are used to capture the vortex structures around the hull.The large separated flow is able to be investigated by these three methods,in which more vortex structures are captured byΩ_(R) approach and Liutex method with scalar,vector and tensor form seems to be more suitable for analyzing the flow mechanism around the hull in dynamic pure yaw test.In general,each vortex structure corresponds to a dominant positive/negative axial Liutex and a bound vortex pair.The streamlines are spiral in the large separated flow,indicating that the flow in corresponding region is rotational.But the rotation of the flow is not directly related to the intensity of Liutex.