Comparative assessment of SAS and DES turbulence modeling for massively separated flows

Numerical studies of the flow past a circular cylinder at Reynolds number 1.4 × 105 and NACA0021 airfoil at the angle of attack 60° have been carried out by scale-adaptive simulation （SAS） and detached eddy simu- lation （DES）, in comparison with the existing experimental data. The new version of the model developed by Egorov and Menter is assessed, and advantages and disadvantages of the SAS simulation are analyzed in detail to provide guidance for industrial application in the future. Moreover, the mechanism of the scale-adaptive characteristics in separated regions is discussed, which is obscure in previous analyses. It is con- cluded that： the mean flow properties satisfactorily agree with the experimental results for the SAS simulation, although the prediction of the second order turbulent statistics in the near wake region is just reasonable. The SAS model can produce a larger magnitude of the turbulent kinetic energy in the recir- culation bubble, and, consequently, a smaller recirculation region and a more rapid recovery of the mean velocity out- side the recirculation region than the DES approach with the same grid resolution. The vortex shedding is slightly less irregular with the SAS model than with the DES approach, probably due to the higher dissipation of the SAS simulation under the condition of the coarse mesh.

Applications of URANS on predicting unsteady turbulent separated flows

Accurate prediction of unsteady separated turbulent flows remains one of the toughest tasks and a practi cal challenge for turbulence modeling. In this paper, a 2D flow past a circular cylinder at Reynolds number 3,900 is numerically investigated by using the technique of unsteady RANS （URANS）. Some typical linear and nonlinear eddy viscosity turbulence models （LEVM and NLEVM） and a quadratic explicit algebraic stress model （EASM） are evaluated. Numerical results have shown that a high-performance cubic NLEVM, such as CLS, are superior to the others in simulating turbulent separated flows with unsteady vortex shedding.

NUMERICAL SOLUTION OF NAVIER-STOKES EQUATION OF UNSTEADY SEPARATED FLOWS DUE TO A SPOILER'S OSCILLATION

The finite difference method (FDM) is applied in the present paper to solve the unsteady NHS equations for incompressible fluids. ADI and SLOR methods are served for the vorticity equation and the Poisson equation for ψ respectively. The upwind scheme is used for the convective terms. The moving boundary conditions are specially treated, and the effects of outlet conditions on the flow field are abo examined. Numerical results obtained show that the spoiler's oscillation induces forming, growing and shedding of the vortices. The shedding frequency of vortices is equal to that of the spoiler's oscillation. The forced unsteady separated flows under the present investigation depend mainly on the reduced frequency. At low reduced frequency, the vortices shed from the spoiler interact weakly with each other, and move downstream at an almost uniform speed of 038 V∞. At high reduced frequency, the interaction between the adjacent vortices strengthens. They close up to and rotate around each other, and eventually, merge into one vortex.

INVESTIGATION OF PULSED-WIRE TECHNIQUE FOR WALL PARAMETERS MEASUREMENT IN SEPARATED FLOWS

A pulsed-wire wall probe measurement system was developed in this paper,which can be used for measuring wall Parameters in separated flow- The operating princi-ple was described and the way of probe calibration was given. Wall parameters of back-ward-facing and forward-facing step flow were measured,and the wall nows structure andcharacters were revealed.

Reconstruction of skin friction topology in complex separated flows

This paper describes a theoretical method for reconstruction of the skin friction topology in complex separated flows,which is developed based on the exact relation between skin friction and surface pressure through the boundary enstrophy flux(BEF).The key of this method is that a skin friction field is reconstructed from a surface pressure field as an inverse problem by applying a variational method.For applications,the approximate method is proposed,where the composite surface pressure field is given by a linear superposition of the base-flow surface pressure field and the surface pressure variation field and the base-flow BEF field is used as the first-order approximation.This approximate method is constructive in a mathematical sense since a complex skin friction field in separated flows can be reconstructed from some elemental skin friction structures(skin friction source/sink,vortex and their combinations)by a linear superposition of some simple surface pressure structures.The distinct topological features,such as critical points,separation lines and attachment lines,naturally occur as a result of such reconstruction.As examples,some elemental skin friction structures in separated flows are reconstructed in simulations,and the skin friction fields in shock-wave/boundary-layer interactions(SWBLIs)are reconstructed from pressure sensitive paint(PSP)images obtained in wind tunnel experiments.

Investigation of the interaction between NS-DBD plasma-induced vortexes and separated flow over a swept wing

The effect of nanosecond pulsed dielectric barrier discharge(NS-DBD) plasma flow separation control is closely related to the actuation frequency,because it involves the interaction between plasma-induced vortexes and separated flow.In order to study the mechanism of NS-DBD plasma flow separation control over a swept wing,especially the influence of the actuation frequency,at first,experimental studies of the actuation frequencies at 100 Hz are conducted to validate the numerical simulation method.Then,numerical studies of different actuation frequencies which are 50 Hz,100 Hz,160 Hz,200 Hz,500 Hz,and 1000 Hz,respectively are conducted.The interaction between the plasma-induced vortexes and the separated flow is analyzed.Results show that there is a range of the actuation frequency which includes the frequency(160 Hz) calculated by the average aerodynamic chord length to make the control effect good,but when the actuation frequencies are too low(50 Hz) or too high(1000 Hz),the control effect will get worse.The former is because plasmainduced vortexes disappear in a period within an actuation cycle;the latter is because plasma-induced vortexes cannot develop completely,resulting in a weak vortex intensity.

Analysis and experimental study on resistance-increasing behavior of composite high efficiency autonomous inflow control device

Bottom water coning is the main reason to reduce the recovery of horizontal bottom water reservoir. By water coning, we mean the oil-water interface changes from a horizontal state to a mound-shaped cone and breaks through to the wellbore. Autonomous inflow control device(AICD) is an important instrument maintain normal production after bottom water coning, however, the resistance increasing ability of the swirl type AICD is insufficient at present, which seriously affects the water control effect. Aiming this problem, this paper designs a multi-stage resistance-increasing and composite type AICD. The separation mechanism of oil-water two phases in this structure, the resistance form of oil-water single phase and the resistance-increasing principle of water phase are analyzed. Establishing the dual-phase multi-stage separation and resistance-increasing model, and verified by measuring the throttling pressure drop and oil-water volume fraction of the AICD, it is found that the composite type AICD has the effect of ICD and AICD at the same time, which can balance the production rate of each well section at the initial stage of production, delay the occurrence of bottom water coning. In the middle and later stages of production, water-blocking can be effectively increased to achieve water control and stable production.After structural sensitivity analysis, the influence law of various structural parameters on the water control performance of composite AICD was obtained. The simulation calculation results show that,compared with the existing swirl type AICD, composite AICD has higher sensitivity to moisture content,the water phase throttling pressure drop is increased by 4.5 times on average. The composite AICD is suitable for the entire stage of horizontal well production.

A Lagrangian criterion of unsteady flow separation for two-dimensional periodic flows

The present paper proposes a Lagrangian criterion of unsteady flow separation for two-dimensional periodic flows based on the principle of weighted averaging zero skin-friction given by Haller （HALLER, G. Exact theory of unsteady separation for two-dimensional flows. Journal of Fluid Mechanics, 512, 257-311 （2004））. By analyzing the distribution of the finite-time Lyapunov exponent （FTLE） along the no-slip wall, it can be found that the periodic separation takes place at the point of the zero FTLE. This new criterion is verified with an analytical solution of the separation bubble and a numerical simulation of lid-driven cavity flows.

TURBULENT SEPARATED REATTACHED FLOW IN A TWO-DIMENSIONAL CURVED-WALL DIFFUSER

A turbulent separation-rcattachment flow in a two-dimensional asymmetrical curved-wall diffuser is studied by a two-dimensional laser doppler velocimeter.The turbulent boundary layer separates on the lower curved wall under strong pressure gradient and then reattaches on a parallel channel.At the inlet of the diffuser,Reynolds number based on the diffuser height is 1.2×10~5 and the velocity is 25.2m/s.The re- sults of experiments are presented and analyzed in new defined streamline-aligned coordinates.The experiment shows that after Transitory Detachment Reynolds shear stress is negative in the near-wall backflow region. Their characteristics are approximately the same as in simple turbulent shear layers near the maximum Reynolds shear stress.A scale is formed using the maximum Reynolds shear stresses.It is found that a Reynolds shear stress similarity exists from separation to reattachment and the Schofield-Perry velocity law ex- ists in the forward shear flow.Both profiles are used in the experimental work that leads to the design of a new eddy-viscosity model.The length scale is taken from that developed by Schofield and Perry.The composite velocity scale is formed by the maximum Reynolds shear stress and the Schofield Perry velocity scale as well as the edge velocity of the boundary layer.The results of these experiments are presented in this paper

Numerical Simulation and Experimental Investigation of 3 D Separated Flow Field around a Blunt Body

To study the effects of factors such as electromagnetic field and electron temperature on the ion extraction characteristics (extraction time, collision loss ratio) in atomic vapor laser isotope separation (AVLIS), 2 D electron equilibrium model was used to analyze and simulate the ion extraction process. The collisions between particles mainly considered charge exchanges between isotopic ions and atoms, which were treated by using cross section. The results show that the electric field and electron temperature have significant effects on the ion extraction characteristics: the stronger the electric field is, the higher the electron temperature is; the shorter the extraction time is, the less the collision loss ratio is, and moderate constant magnetic field has few effects on them. Key words atomic vapor laser isotope separation (AVLIS); ion extraction; charge exchange cross section; collision loss ratio

ON SEPARATED SHEAR LAYER OF BLUNT CIRCULAR CYLINDER

Separated shear layer of blunt circular cylinder has been experimentally investigated for the Reynolds numbers (based on the diameter) ranging from 2.8 x 10(3) to 1.0 x 10(5), with emphasis on evolution of separated shear layer, its structure and distribution of Reynolds shear stress and turbulence kinetic energy. The results demonstrate that laminar separated shear layer experiences 2 similar to 3 times vortex merging before it reattaches, and turbulence separated shear layer takes 5 similar to 6 times vortex merging. In addition, relationship between dimensionless initial frequencies of K-H instability and Reynolds numbers is identified, and reasons for the decay of turbulence kinetic energy and Reynolds shear stress in reattachment region are discussed.

SEPARATED FLOW AROUND BENT-NOSE BICONIC IN HYPERSONIC FLOW

This article investigates the separated flow around the bent-nose biconic, in which the second ordered TVD(Total Variation Diminishing) method is applied. The governing cqua-tion is a thin layer Navier-Stokes equation. The surface pressure distribution at the leeward side is a little higher in the vicinity of the symmetrical plane. This seems to be caused by the development of the vortex at the leeward side due to the separated flow. The location of a separation line and the size are in very good agreement with the experiment. Also the secondary separation is captured very clearly.

Investigation of Separated Flow around a Curved Air Intake

This paper prcsents an investigation of two-dimensional separated flow in and around a submerged curved intake by experiment and computation. The eniptic Navier-Stokes equations are employed which are discreted in body-fitted coordinate system by the SIMPLE method. The results show that the flow separation exists near and behind the inlet and the duct flow is distorted. Comparison between the computational and experimental data is fairly satisfactory.

FLUID DYNAMICS IN A GAS-STIRRED LADLE——A Separated Flow Model with Stochastical Trajectories

A separated flow model with stochastical trajectories has been developed to describe the fluid flow in a bubble stirred ladle.The bubble dispersion,turbulent characteristics and gas-liquid interactions can be predicted by this mathematical model.The bubble flow as a dispersed phase is treated in a Lagrangian frame of reference and the analysis of the turbulent flow for liquid phase is conducted in a Eulerian field.The interactions between bubbles and liquid phases are considered as a bubble source term in the control equation for a continuous phase. The Monte Carlo sampling method is used to determine the bubble trajectories.The homoge- neous flow model is also taken into consideration so that it can be compared with the sepa- rated flow model.Numerical predictions using a water model of a ladle show that the pre- dicted results of the separated flow model agree satisfactorily with the experimental results, but the prediction of the homogeneous flow model are not in good agreement with the experi- mental results.

NATURE OF THE SURFACE HEAT TRANSFER FLUCTUATION IN A HYPERSONIC SEPARATED TURBULENT FLOW

This paper presents the results of an experimental study of the unsteady nature of a hypersonic sepa- rated turbulent flow.The nominal test conditions were a freestream Mach number of 7.8 and a unit Reynolds number of 3.5x10^7/m.The separated flow was generated using finite span forward facing steps.An array of flush mounted high spatial resolution and fast response platinum film resistance thermometers was used to make mul- ti-channel measurements of the fluctuating surface heat trtansfer within the separated flow.Conditional sampling ana- lysis of the signals shows that the root of separation shock wave consists of a series of compression wave extending over a streamwise length about one half of the incoming boundary layer thickness.The compression waves con- verge into a single leading shock beyond the boundary layer.The shock structure is unsteady and undergoes large-scale motion in the streamwise direction.The length scale of the motion is about 22 percent of the upstream influence length of the separation shock wave.There exists a wide band of frequency of oscillations of the shock system.Most of the frequencies are in the range of 1-3 kHz.The heat transfer fluctuates intermittently between the undisturbed level and the disturbed level within the range of motion of the separation shock wave.This inter mittent phenomenon is considered as the consequence of the large-scale shock system oscillations.Downstream of the range of shock wave motion there is a separated region where the flow experiences continuous compression and no intermittency phenomenon is observed.

Use of Targeted Orographic Smoothing in Very High Resolution Simulations of a Downslope Windstorm and Rotor in a Sub-tropical Highland Location

Nested simulations of a downslope windstorm over Cangshan mountain,Yunnan,China,have been used to demonstrate a method of topographic smoothing that preserves a relatively large amount of terrain detail compared to typical smoothing procedures required for models with terrain-following grids to run stably.The simulations were carried out using the Met Office Unified Model(MetUM)to investigate downslope winds.The smoothing method seamlessly blends two terrain datasets to which uniform smoothing has been applied—one with a minimum of smoothing,the other smoothed more heavily to remove gradients that would cause model instabilities.The latter dataset dominates the blend where the steepest slopes exist,but this is localised and recedes outside these areas.As a result,increased detail is starkly apparent in depictions of flow simulated using the blend,compared to one using the default approach.This includes qualitative flow details that were absent in the latter,such as narrow shooting flows emerging from roughly 1-2 km wide leeside channels.Flow separation is more common due to steeper lee slopes.The use of targeted smoothing also results in increased lee side temporal variability at a given point during the windstorm,including over flat areas.Low-/high-pass filtering of the wind perturbation field reveals that relative spatial variability above 30 km in scale(reflecting the background flow)is similar whether or not targeting is used.Beneath this scale,when smoothing is targeted,relative flow variability decreases at the larger scales,and increases at lower scales.This seems linked to fast smaller scale flows disturbing more coherent flows(notably an along-valley current over Erhai Lake).Spatial variability of winds in the model is unsurprisingly weaker at key times than is observed across a local network sampling mesoscale variation,but results are compromised due to relatively few observation locations sampling the windstorm.Only when targeted smoothing is applied does the model capture the downslope windstorm's extension over the city of Dali at the mountain's foot,and the peak mean absolute wind.

Effect of leading-edge tubercles on the flow over low-aspect-ratio wings at low Reynolds number

Two-dimensional time-resolved particle image velocimetry(TR-PIV)and stereographic particle image velocimetry(SPIV)techniques were used to investigate the effect of leading-edge tubercles on the flow over low-aspect-ratio wing models.The angle of attack is fixed at 10°,and the Reynolds number based on chord length is 5.8×10^(3).It is shown that the leading-edge tubercles can effectively mitigate flow separation in the model and also reduce the contribution of wake vortex to the fluctuating energy of flow.Counter-rotating vortex pairs(CVPs)initiated from the peak of leading-edge tubercles can promote nearby momentum exchange,enhance mixing of the flow and increase the energy contained in the boundary layer,which results in resisting the larger adverse pressure gradient.Therefore,it is concluded that CVPs play an important role in mitigating the flow separation for wings with leading-edge tubercles.

Conjugate Heat Transfer Predictions of Gas Turbine Hot Walls Jets Cooling: Influence of Short Hole Grid Resolutions Using Computational Fluid Dynamics*

Short hole investigations relevant to gas turbine (GT) hot walls cooling heat transfer techniques, were carried out using computational fluid dynamics (CFD) combined with conjugate heat transfer (CHT) code. The CFD software are commercial ones: ICEM for grid modelling and ANSYS Fluent for the numerical calculation, where symmetrical application prevails. The CFD CHT predictions were undertaken for Nimonic-75 metal walls with square (152.4 mm) arrays of 10 holes, whereby the lumped heat capacitance method was applied in order to determine the surface average heat transfer coefficient (HTC), h (W/m2 K) and the dimensionless Nusselt number, Nu. The major parameters considered for the short hole geometries are the pitch to diameter, X/D and length to diameter, L/D ratios and both were varied with range of D values, but X of 15.24 mm and L of 6.35 mm kept constant. Also applied, are variable mass flux, G (kg/s∙m2) and were used in predicting the flow aerodynamics in the short holes. The predictions were for classic thermal entry length into a round hole, as vena contracta, flow separation and reattachment dominates the holes, hence the development of thermal profile through the depth of the GT hot walls. Additionally, the acceleration of the flow along the wall surfaces as it approaches the holes, was a significant part of the overall heat transfer. This was shown to be independent of the hole length, even though the L/D parameter is a critical component to enhanced heat transfer. The CFD CHT predictions showed that validation of the HTC h, Nu and pressure loss, ∆P are in better agreement with measured data and within reasonable acceptance. The ∆P agreement signifies that the aerodynamics were predicted correctly, which is also the reason why the HTC expressed per wall hole approach surface area and Nu were better predicted. This illustrates how effective and efficient the wall internal heat transfer cooling is for gas turbine hot wall heat transfer using airflow jets cooling.

Analysis of the Oscillatory Flows of Multiple Shock Waves in a Constant Area Duct

The oscillatory response of multiple shock waves to downstream perturbations in a supersonic flow is studied numerically in a rectangular duct.Multiple shock waves are formed inside the duct at a shock Mach number of 1.75.The duct has an exit height of H,and the effect of duct resonance on multiple shock oscillations is investigated by attaching exit ducts of lengths 0H,50H,and 150H.The downstream disturbance frequency varied from 10 Hz to 200 Hz to explore the oscillation characteristics of the multiple shock waves.The oscillatory response of shock waves under self-excited and forced oscillation conditions are analyzed in terms of wall static pressure,shock train leading-edge location,shock train length,and the size of the separation bubble.The extent of the initial shock location increases with an increase in exit duct length for the self-excited oscillation condition.The analysis of the shock train leading edge and the spectral analysis of wall static pressure variations are conducted.The variation in the shock train length is analyzed using the pressure ratio method for self-excited as well as forced oscillations.The RMS amplitude of the normalized shock train length(ζ_(ST))increases with an increase in the exit duct length for the self-excited oscillation condition.When the downstream perturbation frequency is increased,ζ_(ST)is decreased for exit duct configurations.For all exit duct designs and downstream forcing frequencies,the size of the separation bubble grows and shrinks during the shock oscillations,demonstrating the dependence on duct resonance and forced oscillations.

Flow control of micro-ramps on supersonic forward-facing step flow

The effects of the micro-ramps on supersonic turbulent flow over a forward-facing step（FFS） was experimentally investigated in a supersonic low-noise wind tunnel at Mach number 3 using nano-tracer planar laser scattering（NPLS）and particle image velocimetry（PIV） techniques. High spatiotemporal resolution images and velocity fields of supersonic flow over the testing model were captured. The fine structures and their spatial evolutionary characteristics without and with the micro-ramps were revealed and compared. The large-scale structures generated by the micro-ramps can survive the downstream FFS flowfield. The micro-ramps control on the flow separation and the separation shock unsteadiness was investigated by PIV results. With the micro-ramps, the reduction in the range of the reversal flow zone in streamwise direction is 50% and the turbulence intensity is also reduced. Moreover, the reduction in the average separated region and in separation shock unsteadiness are 47% and 26%, respectively. The results indicate that the micro-ramps are effective in reducing the flow separation and the separation shock unsteadiness.