Detached Eddy Simulation方法模拟不同类型翼型的失速特性

应用DES(DetachedEddySimulation)方法数值模拟了3种不同失速类型的翼型的升力特性。DES方法结合了RANS(ReynoldsaveragedNavierStokes)和LES(LargeEddySimulationapproaches)的优点。基于SpalartAllmaras湍流模型,在近壁面DES体现为RANS模型的特点而在远离物面处又具有LES的亚格子模型的特性。对此模型使用了LUSGS隐式格式求解。通过和实验结果对比,显示这种方法可以有效地预测翼型的失速特性。

Analysis of the Pump-turbine S Characteristics Using the Detached Eddy Simulation Method

Current research on pump-turbine units is focused on the unstable operation at off-design conditions, with the characteristic curves in generating mode being S-shaped. Unlike in the traditional water turbines, pump-turbine operation along the S-shaped curve can lead to difficulties during load rejection with unusual increases in the water pressure, which leads to machine vibrations. This paper describes both model tests and numerical simulations. A reduced scale model of a low specific speed pump-turbine was used for the performance tests, with comparisons to computational fluid dynamics（CFD） results. Predictions using the detached eddy simulation（DES） turbulence model, which is a combined Reynolds averaged Naviers-Stokes（RANS） and large eddy simulation（LES） model, are compared with the two-equation turbulence mode results. The external characteristics as well as the internal flow are for various guide vane openings to understand the unsteady flow along the so called S characteristics of a pump-turbine. Comparison of the experimental data with the CFD results for various conditions and times shows that DES model gives better agreement with experimental data than the two-equation turbulence model. For low flow conditions, the centrifugal forces and the large incident angle create large vortices between the guide vanes and the runner inlet in the runner passage, which is the main factor leading to the S-shaped characteristics. The turbulence model used here gives more accurate simulations of the internal flow characteristics of the pump-turbine and a more detailed force analysis which shows the mechanisms controlling of the S characteristics.

Comparison of Passive Control Methods on Cavity Aeroacoustic Using Delayed Detached Eddy Simulation

Delayed detached eddy simulation(DDES)is performed to investigate an open cavity at Ma0.85.Clean cavity and cavity with leading-edge saw tooth spoiler and flattop spoiler,are modeled.The results obtained from clean cavity prediction are compared with experimental sound pressure level(SPL)data from QinetiQ,UK.DDES results agree well with the experimental data.Furthermore,comparisons are made with the predicted SPL between the three configurations to find out the effect of different passive control methods.Both the spoilers can suppress the over-all SPL up to 8dB.The main focuses of this investigation are to exam the DDES method on cavity aeroacoustic analysis and test the noise suppression effect by saw tooth spoiler and flattop spoiler.

Detached eddy simulation on the structure of swirling jet flow field

The improved delayed detached eddy simulation method with shear stress transport model was used to analyze the evolution of vortex structure,velocity and pressure fields of swirling jet.The influence of nozzle pressure drop on vortex structure development and turbulence pulsation was investigated.The development of vortex structure could be divided into three stages:Kelvin-Helmholtz(K-H)instability,transition stage and swirling flow instability.Swirling flow could significantly enhance radial turbulence pulsation and increase diffusion angle.At the downstream of the jet flow,turbulence pulsation dissipation was the main reason for jet velocity attenuation.With the increase of pressure drop,the jet velocity,pulsation amplitude and the symmetry of velocity distribution increased correspondingly.Meanwhile the pressure pulsation along with the axis and vortex transport intensity also increased significantly.When the jet distance exceeded about 9 times the dimensionless jet distance,the impact distance of swirling jet could not be improved effectively by increasing the pressure drop.However,it could effectively increase the swirl intensity and jet diffusion angle.The swirling jet is more suitable for radial horizontal drilling with large hole size,coalbed methane horizontal well cavity completion and roadway drilling and pressure relief,etc.

Unsteady RANS and detached eddy simulation of the multiphase flow in a co-current spray drying

A detached eddy simulation(DES) and a k-ε-based Reynolds-averaged Navier–Stokes(RANS) calculation on the co-current spray drying chamber is presented. The DES used here is based on the Spalart–Allmaras(SA) turbulence model, whereas the standard k-ε(SKE) was considered here for comparison purposes. Predictions of the mean axial velocity, temperature and humidity profile have been evaluated and compared with experimental measurements. The effects of the turbulence model on the predictions of the mean axial velocity, temperature and the humidity profile are most noticeable in the(highly anisotropic) spraying region. The findings suggest that DES provide a more accurate prediction(with error less than 5%) of the flow field in a spray drying chamber compared with RANS-based k-ε models. The DES simulation also confirmed the presence of anisotropic turbulent flow in the spray dryer from the analysis of the velocity component fluctuations and turbulent structure as illustrated by the Q-criterion.

Complicated flow in tip flow field of a compressor tandem cascade using delayed detached eddy simulation

The complex flow characteristics in the tip region of a tandem cascade with tip clearance have been calculated and analyzed using Delayed Detached Eddy Simulation(DDES).The coherent mechanism of the vortex structures near the blade tip was discussed,and the unsteady behaviors and features in the tip flow field were analyzed.Additionally,the interaction between the tip leakage flow and the gap jet was revealed.The results show that,compared to the datum cascade,the blade tip load of the rear blade increases while that of the front blade decreases.Unsteady fluctuations of the tandem cascade are mainly caused by the interaction between the tip leakage flow and gap jet,and by the mixing of the vortex structures,but there is no essential change in the spectrum feature of the tip leakage flow.Finally,a detailed analysis of the development of vortices in the tip region is conducted by the topological structures of the flow field.Combined with the three-dimensional vortex structures,the schematic diagram of the vortex system of the datum single-row cascade and tandem cascade is summarized.

DETACHED EDDY SIMULATION OF HYDRAULIC CHARACTERISTICS ALONG THE SIDE-WALL AFTER A NEW ARRANGEMENT-SCHEME OF THE SUDDEN LATERAL ENLARGEMENT AND THE VERTICAL DROP

In order to study the cavitation damage in a side-wall when a sudden lateral enlargement and a vertical drop are imposed at the radial gate, a new arrangement-scheme is proposed, where the sudden lateral enlargement and the vertical drop can be imposed at the outlet of the gate chamber. The hydraulic characteristics along the side-wall are simulated by the detached eddy simulation and the Volume Of Fluid （VOF） method. The numerical results agree well with those of experiment. The experimental and numerical results show that the flow condition is smooth with only a weak water-wing appearing behind the lateral cavity, and the length of the lateral cavity becomes longer and is mainly affected by the size of the lateral enlargement and the zone of negative pressure after the water impacts the side-wall would disappear. The hydraulic characteristics of the new arrangement-scheme are beneficial to the prevention of the cavitation damage in the side-wall and the Detached Eddy Simulation （DES） with the VOF method can well predict the hydraulic characteristics after the new arrangement-scheme of the sudden lateral enlargement and the vertical drop.

Delayed Detached Eddy Simulation of Subcritical Flow past Generic Side Mirror

In the present study, the subcritical flow past a generic side mirror on a base plane is investigated at the Reynolds number of 5.2 × 10~5 using delayed detached eddy simulation(DDES) turbulence model. Asides from the capability of capturing main features of the large recirculation vortex in the wake of the side mirror and the front horseshoe vortex, the accuracy of DDES estimation of recirculation length is significantly increased by over20%, compared to the detached eddy simulation(DES) estimation using the same grid. And DDES prediction of pressure coefficient at the trailing edge of the mirror is in good agreement with the experiments, which is more accurate than both DES and large eddy simulation(LES) results. The results verify the capacity of DDES turbulence model to solve the turbulent flow around the side mirror. This is a key foundation for possible future study of full simulation of external flow field of vehicle.

Investigation of film cooling on the leading edge of turbine blade based on detached eddy simulation

In order to assess the influences of curved hole passage on cooling effectiveness and flow structure of turbine blade leading edge,the detached eddy simulation is applied to numerically investigate the AGTB turbine cascade under the condition of global blowing ratio M=0.7.The straight or curved cooling holes are located at either the pressure or suction side near the leading edge.The analysis and discussion focus on the local turbulence structure;influence of pressure gradient on the structure,and distribution of cooling effectiveness on the blade surface.The numerical results show that cooling hole with curved passage could bring positive impact on the increase of the local cooling effectiveness.On the suction side,the increased cooling effectiveness could be about 82% and about 77% on the pressure side,compared to the conventional straight hole.

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.

Study on the Fractal Characteristics of Cavitation Shedding over a Twisted Hydrofoil

Cavitation and cavitation erosion often occur and seriously threaten the safe and stable operation of hydraulic machinery.However,during the operation of hydraulic machinery,the cavitation flow field is often difficult to contact and measure,and the shedding and development characteristics of cavitation flow are unknown.This paper uses the Detached Eddy Simulation(DES)turbulence model and Zwart-Gerber-Belamri(ZGB)cavitation model to conduct numerical research on the cavitation flow of a twisted hydrofoil and verifies the effectiveness of numerical simulation by comparing it with experimental results.Then,based on the fractal dimension method,the number and fractal dimension of the main cavitation clusters in two different views are analyzed.The research results show that with the continuous flow of fluid,the cavitation at the tail of the hydrofoil periodically falls off,and the number of cavitation bubbles shows a strong randomness between 1 and 2.The fractal dimension method calculates that the frequency of cavitation bubbles falling off is 30 Hz.The fractal dimension of cavitation also fluctuates periodically.The fractal dimension of cavitation in the side view fluctuates between 1.30–1.38,and the fractal dimension of cavitation in the top view fluctuates between 1.31–1.42.The fractal dimension of the cavitation is closely related to the size and shape of the cavitation.By analyzing the relationship between the cavitation development process and the fractal dimension in the side view,we found the relationship between the development of cavitation in a cycle and the change of fractal dimension,which is of great significance to understanding the characteristics of cavitation flow through the image of the cavitation flow field.

Hybrid URANS/LES Method of Flow Fields in Axial-flow Compressor Rotor Rotor

Accurate and efficient prediction of the aerodynamic performance and flow details of axial-flow com-pressors is of great engineering application value for the aerodynamic design and flow control of axial-flow compres-sors.In this work,a delayed detached eddy simulation method is developed and applied to numerically simulate the tur-bulent channel flow and the aerodynamic performance of NASA Rotor 35.Several acceleration techniques including parallel implementation are also used to speed up the iteration convergence.The mean velocity distribution and Reyn-olds stress distribution in the boundary layer of turbulent channel flow and the aerodynamic performance curve of NASA Rotor 35 are predicted.The good agreement between the present delayed detached eddy simulation results and the available direct numerical simulation results or experimental data confirms the effectiveness of the developed meth-od in the accurate and efficient prediction of complex flow in turbomachinery.

翼型失速特性的RANS方法与DES方法数值模拟对比分析

RAN S(R eyno lds-averaged N av ier-Stokes)加湍流模型是当前计算飞机粘性流场的最常用方法,数值实践说明要计算大分离流动,需要更高级的方法例如LES(Large Eddy S im u lation)或DN S(D irect N S S im u lation)。然而实际雷诺数下,LES和DN S对网格的要求太高,以至目前还难以应用。DES(D etached-Eddy S im u lation)方法结合了RAN S和LES的优点,通过对Spalart-A llm aras湍流模型中长度尺度的修正,在近壁面它体现为RAN S模型的特点,而在远离物面处又保持LES的亚格子模型的特性。论文对比了采用RAN S和DES方法数值模拟翼型失速特性的能力,并与实验结果进行了对比。结果表明,对大分离流动的数值模拟,DES方法体现出更强的能力。

Numerical Prediction of Marine Propeller Noise in Non-Uniform Inflow

A numerical study on the acoustic radiation of a propeller interacting with non-uniform inflow has been conducted. Real geometry of a marine propeller DTMB 4118 is used in the calculation, and sliding mesh technique is adopted to deal with the rotational motion of the propeller. The performance of the DES （Detached Eddy Simulation） approach at capturing the unsteady forces and moments on the propeller is compared with experiment. Far-field sound radiation is predicted by the formation 1A developed by Farassat, an integral solution of FW-H （Ffowcs Williams-Hawkings） equation in time domain. The sound pressure and directivity patterns of the propeller operating in two specific velocity distributions are discussed.

Location of anemometer along Lanzhou-Xinjiang railway

Using structured mesh to discretize the calculation region, the wind velocity and pressure distribution in front of the wind barrier under different embankment heights are investigated based on the Detached Eddy Simulation(DES) with standard SpalartAllmaras(SA) model. The Reynolds number is 4.0×105 in this calculation. The region is three-dimensional. Since the wind barrier and trains are almost invariable cross-sections, only 25 m along the track is modeled. The height of embankment ranges from 1 m to 5 m and the wind barrier is 3 m high. The results show that the wind speed changes obviously before the wind barrier on the horizontal plane, which is 4.5 m high above the track. The speed of wind reduces gradually while approaching the wind barrier. It reaches the minimum value at a distance about 5 m before the wind barrier, and increases dramatically afterwards. The speed of wind at this location is linear with the speed of far field. The train aerodynamic coefficients decrease sharply with the increment of the embankment height. And they take up the monotonicity. Meanwhile, when the height increases from 3 m to 5 m, they just change slightly. It is concluded that the optimum anemometer location is nearly 5 m in front of the wind barrier.

Detached-eddy simulation of wing-tip vortex in the near field of NACA 0015 airfoil

In the present study, the formation of the wing-tip vortex from a rectangular NACA0015 wing with a square tip at the Reynolds number of 1.8× 105 and the angles of attack （AOA） α = 8° and 10° were simulated with an incompressible detached eddy simulation （DES） method and the Reynolds averaged Navier-Stokes （RANS） equations with the SA model respectively. Numerical results were compared with experimental results to validate the capability of the employed methods in resolving tip vortex flows. The results show that DES model could capture the complicated three-dimensional structures in the vortex, and the streamwise vorticity and the cross-flow velocity agree with the experiment results quite well, but RANS-SA model with the same grid as that of DES failed to capture the correct structures and under-predicted the streamwise vorticity in the vortex by 40%. The present study suggests that under the same calculation cost, DES but not RANS-SA could be used to effectively predict the flow characteristics in tip vortex.

Influence of ribs on train aerodynamic performances

The influence of ribs on the train aerodynamic performance was computed using detached eddy simulation(DES), and the transient iteration was solved by the dual-time step lower-upper symmetric Gauss-Seidel(LU-SGS) method. The results show that the ribs installed on the roof have a great effect on the train aerodynamic performance. Compared with trains without ribs, the lift force coefficient of the train with convex ribs changes from negative to positive, while the side force coefficient increases by 110%and 88%, respectively. Due to the combined effect of the lift force and side force, the overturning moment of the train with convex ribs and cutting ribs increases by 140% and 106%, respectively. There is larger negative pressure on the roof of the train without ribs than that with ribs. The ribs on the train would disturb the flow structure and contribute to the air separation, so the separation starts from the roof, while there is no air separation on the roof of the train without ribs. The ribs can also slow down the flow speed above the roof and make the air easily sucked back to the train surface. The vortices at the leeward side of the train without ribs are small and messy compared with those of the train with convex or cutting ribs.

CFD Simulation of Local and Global Mixing Time in an Agitated Tank

The Issue of mixing efficiency in agitated tanks has drawn serious concern in many industrial processes. The turbulence model is very critical to predicting mixing process in agitated tanks. On the basis of computational fluid dynamics（CFD） software package Fluent 6.2, the mixing characteristics in a tank agitated by dual six-blade-Rushton-turbines（6-DT） are predicted using the detached eddy simulation（DES） method. A sliding mesh（SM） approach is adopted to solve the rotation of the impeller. The simulated flow patterns and liquid velocities in the agitated tank are verified by experimental data in the literature. The simulation results indicate that the DES method can obtain more flow details than Reynolds-averaged Navier-Stokes（RANS） model. Local and global mixing time in the agitated tank is predicted by solving a tracer concentration scalar transport equation. The simulated results show that feeding points have great influence on mixing process and mixing time. Mixing efficiency is the highest for the feeding point at location of midway of the two impellers. Two methods are used to determine global mixing time and get close result. Dimensionless global mixing time remains unchanged with increasing of impeller speed. Parallel, merging and diverging flow pattern form in the agitated tank, respectively, by changing the impeller spacing and clearance of lower impeller from the bottom of the tank. The global mixing time is the shortest for the merging flow, followed by diverging flow, and the longest for parallel flow. The research presents helpful references for design, optimization and scale-up of agitated tanks with multi-impeller.

Detached-eddy simulation of turbulent coherent structures around groynes in a trapezoidal open channel

The hydrodynamics in a straight open channel with a multiple-embayment groyne field was investigated using the detached-eddy simulation(DES).A series of short groynes were included on a 1:3 side slope of the channel.This work focuses on the turbulent coherent structures around groynes on an uneven bottom.Flows around groyne fields are characterized by massive separation and highly unsteady vortices.DES can capture a wide spectrum of eddies at a lower computational cost than the large eddy simulation(LES)or direct numerical simulation(DNS).In the present work,a zonal DES model(ZDES)was used to simulate the flow around groynes.The ZDES model is a modified version of the DES designed to overcome the model-stress depletion(MSD)of the RANS/LES hybrid model.The vortex system consists of the horseshoe vortex(HV)formed at the base of the obstructions,the necklace vortex(NV)that wrapped the groyne tips near the free surface,and the shedding vortex(SV)underneath the free surface.The effects of the incident flow and local topography on the vortex evolution were investigated by analyzing the mean flow structures and the instantaneous turbulent flow fields.Some important vortices cannot be captured because of the averaging process,while some flow structures cannot be observed in the instantaneous flow.The mean flow is only a reflection of the averaging process when complex vortices are present.

Effect of car-body lower-center rolling on aerodynamic performance of a high-speed train

The interaction between the car-body vibration and aerodynamic performance of the train becomes more prominent motivated by the vehicle’s light-weighting design.To address this topic,this study firstly analyzes the posture characteristics of the car-body based on the previous full-scale test results.And then the aerodynamic performance under different vibration cases(different car-body roll angles)is studied with an improved delayed detached eddy simulation(IDDES).The results revealed that car-body rolling had a significant impact on the aerodynamic behavior of bogies,which significantly increased the lateral force and yaw moment of a bogie and further may have aggravated the operational instability of the train.The unbalanced distribution of the longitudinal pressure on both sides of the bogie caused by the car-body rolling motion was the primary cause for the bogie yaw moment increase.The tail vortex of the train was also affected by the car-body rolling,resulting in vertical jitter.