Simulation of wing-body junction flows with hybrid RANS/LES methods

In this paper, flows past two wing-body junctions, the Rood at zero angle of attack and NASA TN D-712 at 12.5° angle of attack, are investigated with two Reynolds-Averaged Navier-Stokes (RANS) and large eddy simulation (LES) hybrid methods. One is detached eddy simulation (DES) and the other is delayed-DES, both are based on a weakly nonlinear two-equation k–ω model. While the RANS method can predict the mean flow behaviours reasonably accurately, its performance for the turbulent kinetic energy and shear stress, as compared with available experimental data, is not satisfactory. DES, through introducing a length scale in the dissipation terms of the turbulent kinetic energy equation, delivers flow separation, a vortex or the onset of vortex breakdown too early. DDES, with its delayed effect, shows a great improvement in flow structures and turbulence characteristics, and agrees well with measurements.     

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×10³ to 1.0×10⁵, 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–3 times vortex merging before it reattaches, and turbulence separated shear layer takes 5–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. The project supported by the National Natural Science Foundation of China and the Key Laboratory for Hydrodynamics of NDCST.     

Turbulent flow around a wall-mounted cube: A direct numerical simulation

An immersed-boundary method was employed to perform a direct numerical simulation (DNS) of flow around a wall-mounted cube in a fully developed turbulent channel for a Reynolds number Re = 5610, based on the bulk velocity and the channel height. Instantaneous results of the DNS of a plain channel flow were used as a fully developed inflow condition for the main channel. The results confirm the unsteadiness of the considered flow caused by the unstable interaction of a horseshoe vortex formed in front of the cube and on both its sides with an arch-type vortex behind the cube. The time-averaged data of the turbulence mean-square intensities, Reynolds shear stresses, kinetic energy and dissipation rate are presented. The negative turbulence production is predicted in the region in front of the cube where the main horseshoe vortex originates.     

A note on numerical simulation of vortical structures in shock diffraction

In numerical simulation of the Euler equations, the slipstream or shear layer that appears behind a diffracted shock wave may develop small discrete vortices using fine computational meshes. Similar phenomena were also observed in the simulation of a Mach reflection that is accompanied by a shear layer. However, these small vortices have never been observed in any shock-tube experiment, although the wave pattern and the shape of the main vortex agree very well with visualization results. Numerical solutions obtained with coarse grids may agree better with experimental photos than those with very fine grids because of the pollution of the small vortices. This note tries to investigate the effect of viscosity on the small vortices by comparing the solutions of the laminar Navier-Stokes equations and the turbulence model. It is found that the small vortices are still observed in the solution of the laminar Navier-Stokes equations, although they can be suppressed by using the turbulence model. Numerical and experimental factors that are responsible for the deviation of the laminar solutions from experimental results are discussed. The secondary vortex in shock diffraction is successfully simulated by solving the Navier-Stokes equations.Received: 28 March 2003, Accepted: 6 May 2003, Published online: 11 June 2003     

湍流一般机理及其应用

本文综述了湍流机理发展中最重要的一些文章,包括:Brown-Roshko的混合剪切层中大尺度涡的发展;Perry和Chong的湍流边界层中A形涡结构的湍流机理;以及笔者提出的关于3维湍流场中涡结构的涡量脉动对流和扩散并存的一般机理。在一般机理基础上建立了判别流场的准则,列举了在固壁边界附近的猝发现象并展望它的未来应用。      

Skin-friction drag reduction in turbulent channel flow based on streamwise shear control

It is known that stretching and intensification of a hairpin vortex by mean shear play an important role to create a hairpin vortex packet, which generates the large Reynolds shear stress associated with skin-friction drag in wall-bounded turbulent flows. In order to suppress the mean shear at the wall for high efficient drag reduction (DR), in the present study, we explore an active flow control concept using streamwise shear control (SSC) at the wall. The longitudinal control surface is periodically spanwise-arranged with no-control surface while varying the structural spacing, and an amplitude parameter for imposing the strength of the actuating streamwise velocity at the wall is introduced to further enhance the skin-friction DR. Significant DR is observed with an increase in the two parameters with an accompanying reduction of the Reynolds stresses and vorticity fluctuations, although a further increase in the parameters amplifies the turbulence activity in the near-wall region. In order to study the direct relationship between turbulent vortical structures and DR under the SSC, temporal evolution with initial eddies extracted by conditional averages for Reynolds-stress-maximizing Q2 events are examined. It is shown that the generation of new vortices is dramatically inhibited with an increase in the parameters throughout the flow, causing fewer vortices to be generated under the control. However, when the structural spacing is sufficiently large, the generation of new vortex is not suppressed over the no-control surface in the near-wall region, resulting in an increase of the second- and fourth-quadrant Reynolds shear stresses. Although strong actuating velocity intensifies the near-wall turbulence, the increase in the turbulence activity is attributed to the generation of counter-clockwise near-wall vortices by the increased vortex transport.     

Vortex reynolds number in turbulent boundary layers

The effects of vortex Reynolds number on the statistics of turbulence in a turbulent boundary layer have been investigated. Vortex Reynolds number is defined as the ratio of circulation around the vortex structure to the fluid viscosity. The vortex structure of the outer region was modeled and a full numerical simulation was then conducted using a high-order spectral method. A unit domain of the outer region of a turbulent boundary layer was assumed to be composed of essentially three elements: a wall, a Blasius mean shear, and an elliptic vortex inclined at 45° to the flow direction. The laminar base-flow Reynolds number is roughly in the same range as that of a turbulent boundary layer based on eddy viscosity, and the vortex-core diameter based on the boundary-layer thickness is nearly the same as the maximum mixing length in a turbulent boundary layer. The computational box size, namely, 500, 150, and 250 wall units in the streamwise, surface-normal, and spanwise directions, respectively, is approximately the same as the measured quasi-periodic spacings of the near-wall turbulence-producing events in a turbulent boundary layer. The effects of vortex Reynolds number and the signs of the circulation on the moments of turbulence were examined. The signs mimic the ejection and sweep types of organized motions of a turbulent boundary layer. A vortex Reynolds number of 200 describes the turbulence moments in the outer layer reasonably well.     

Linear Dynamics of Perturbations in Flows with Constant Horizontal Shear

The dynamics of perturbations in shallow water and incompressible stratified fluid flows with constant horizontal shear are described using the nonmodal analysis. It is shown that the shear flow perturbations can be divided into two classes on the basis of the potential vorticity: rapidly oscillating wave perturbations with zero potential vorticity and slow vortex perturbations with nonzero potential vorticity. In the cases of weak and strong shear the main features of the dynamics of wave and vortex perturbations are studied analytically (using the WKBG method) and numerically. It is shown that for large times the wave perturbation energy increases linearly, i.e., the shear flow is algebraically unstable due to the growth of rapid wave perturbations. This instability can be of importance in processes of turbulence development and surface and internal wave generation.     

Higher‐order and adaptive discontinuous Galerkin methods with shock‐capturing applied to transonic turbulent delta wing flow

Discontinuous Galerkin (DG) methods allow high‐order flow solutions on unstructured or locally refined meshes by increasing the polynomial degree and using curved instead of straight‐sided elements. DG discretizations with higher polynomial degrees must, however, be stabilized in the vicinity of discontinuities of flow solutions such as shocks. In this article, we device a consistent shock‐capturing method for the Reynolds‐averaged Navier–Stokes and k‐ ω turbulence model equations based on an artificial viscosity term that depends on element residual terms. Furthermore, the DG method is combined with a residual‐based adaptation algorithm that targets at resolving all flow features. The higher‐order and adaptive DG method is applied to a fully turbulent transonic flow around the second Vortex Flow Experiment (VFE‐2) configuration with a good resolution of the vortex system.Copyright © 2013 John Wiley & Sons, Ltd.     

Investigation of Advanced Turbulence Models for the Flow in a Generic Wing-Body Junction

The predictive performance of several turbulence models, among them formulations based on non-linear stress-strain relationships and on stress-transport equations, is examined in a collaborative university-industry study directed towards a generic wing-body junction. The geometry consists of a variation of the symmetric NACA 0020 aerofoil mounted on a flat plate, with the oncoming stream aligned with the aerofoil's symmetry plane. The dominant feature of this flow is a pronounced horseshoe vortex evolving in the junction region following separation ahead of the aerofoil's leading edge. This case is one of 6 forming a broad programme of turbulence-model validation by UMIST, Loughborough University, BAE Systems, Aircraft Research Association, Rolls-Royce plc and DERA. Key aspects of this collaboration were a high level of interaction between the partners, the use of common grids and boundary conditions, and numerical verifications aimed at maximizing confidence in the validity of the computational solutions. In total, 12 turbulence models were studied by four partners. Model performance is judged by comparing solutions with experimental data for pressure fields on the plane wall and around the aerofoil; for velocity, turbulence energy, shear stress and streamwise normal stress in the upstream symmetry plane; and for velocity, turbulence energy and shear stress in cross-flow planes downstream of the aerofoil leading edge. The emphasis of the study is on the structure of the horseshoe vortex and its effects on the forward flow. The main finding of the study is that, for this particular 3D flow, second-moment closure offers predictive advantages over the other models examined, especially in terms of the far-field structure of the horse-shoe vortex, although no model achieves close agreement with the experimental data in respect of both mean flow and turbulence quantities. This revised version was published online in July 2006 with corrections to the Cover Date.     

High-order strand grid methods for shock turbulence interaction

ABSTRACT

In this work, we examine the flux correction method for three-dimensional transonic turbulent flows on strand grids. Building upon previous work, we treat flux derivatives along strands with high-order summation-by-parts operators and penalty-based boundary conditions. A finite-volume like limiting strategy is implemented in the flux correction algorithm in order to sharply capture shocks. To achieve turbulence closure in the Reynolds-Averaged Navier–Stokes equations, a robust version of the Spalart–Allmaras turbulence model is employed that accommodates negative values of the turbulence working variable. Validation studies are considered which demonstrate the flux correction method achieves a high degree of accuracy for turbulent shock interaction flows.     

Vortices in time-periodic shear flow

Vortices emerging in geophysical turbulence may experience deformations due to the non-uniform ambient flow induced by neighbouring vortices. At first approximation this ambient flow is modeled by a linear shear flow. It is well known from previous studies that the vortex may be (partially) destructed through removal of weak vorticity at the vortex edge—a process referred to as ‘stripping’. While most previous studies considered a stationary external shear flow, we have examined the behaviour of the vortex embedded in a linear shear flow whose strength changes harmonically in time. Aspects of the vortex dynamics and the (chaotic) transport of tracers have been studied by both laboratory experiments and numerical simulations based on a simple kinematical model.     

A study of the energetic turbulence structures during stall delay

This study revealed the three-dimensional instantaneous topologies of the large-scale turbulence structures in the separated flow on the suction surface of wind turbine’s blade during stall delay. These structures are the major contributors to the first two POD (proper orthogonal decomposition) modes. The two kinds of instantaneous flow structures as major contributors to the first POD mode are: (1) extended regions of downwash flow with an upstream upward flow beside it and a compact vortex pair closer to the blade’s leading edge; (2) a large-scale clockwise vortex with strong induced flows. The two kinds of flow structures contributing significantly to the second POD mode are: (1) large counter-rotating vortices inducing strong upward velocities and a series of small vortices; (2) strong downwash flow coming from the leading-edge shear layer with a large and strong vortex on the left side and small vortices upstream. The statistical impacts of these large-scale and energetic structures on the turbulence have also been studied. It was observed that when these turbulence structures were removed from the flow, the peak values of some statistics were significantly reduced.     

An experimental investigation of a large-scale structure of a two-dimensional backward-facing step by using advanced multi-point LDV

Measurements of spatio–temporal velocity fields at the separated shear layer and reattachment region of a two-dimensional backward-facing step flow are carried out simultaneously using a multi-point LDV. The objective of this paper is to clarify experimentally the structure of a large-scale structure of this flow field using a space and time correlation and conditional average. From the results of the correlation of the velocity fluctuation, the moving path of the vortex shedding from the separated shear layer to the reattachment region exhibits two patterns which it moves to near the wall region or the middle of the step height at the reattachment region. Especially, it moves to near the wall region when it grows larger in the separated shear layer. Moreover, the turbulence concerned with reattachment phenomenon transports from the reattachment region to a separated shear layer by recirculation flow. According to these transports of turbulence, a model for large-scale fluctuation is proposed as a self-excitation motion.     

Kinematics and dynamics of small-scale vorticity and strain-rate structures in the transition from isotropic to shear turbulence

We applied a technique that defines and extracts “structures” from a DNS dataset of a turbulence variable in a way that allows concurrent quantitative and visual analysis. Local topological and statistical measures of enstrophy and strain-rate structures were compared with global statistics to determine the role of mean shear in the dynamical interactions between fluctuating vorticity and strain-rate during transition from isotropic to shear-dominated turbulence. We find that mean shear adjusts the alignment of fluctuating vorticity, fluctuating strain-rate in principal axes, and mean strain-rate in a way that (1) enhances both global and local alignments between vorticity and the second eigenvector of fluctuating strain-rate, (2) two-dimensionalizes fluctuating strain-rate, and (3) aligns the compressional components of fluctuating and mean strain-rate. Shear causes amalgamation of enstrophy and strain-rate structures, and suppresses the existence of strain-rate structures in low-vorticity regions between enstrophy structures. A primary effect of shear is to enhance “passive” strain-rate fluctuations, strain-rate kinematically induced by local vorticity concentrations with negligible enstrophy production, relative to “active,” or vorticity-generating strain-rate fluctuations. Enstrophy structures separate into “active” and “passive” based on the level of the second eigenvalue of fluctuating strain-rate. We embedded the structure-extraction algorithm into an interactive visualization-based analysis system from which the time evolution of a shear-induced hairpin enstrophy structure was visually and quantitatively analyzed. The structure originated in the initial isotropic state as a vortex sheet, evolved into a vortex tube during a transitional period, and developed into a well-defined horseshoe vortex in the shear-dominated asymptotic state.     

Vortex dynamics of a cylinder wake in proximity to a wall

Large-eddy simulations (LES) are used to investigate the modifications of wake dynamics and turbulence characteristics behind a circular cylinder placed near a wall for varying gap-to-diameter (G/D) ratios (where G signifies the gap between the wall and the cylinder, and D the cylinder diameter). The three-dimensional (3-D), time-dependent, incompressible Navier–Stokes equations with a dynamic subgrid-scale model are solved using a symmetry-preserving finite-difference scheme of second-order spatial and temporal accuracy. The immersed boundary (IB) method is employed to impose the no-slip boundary condition on the cylinder surface. Flow visualizations along with turbulence statistics are presented to gain insight into the flow structures that are due to interaction between the shear layers and the approaching boundary layer. Apart from the vortex shedding mechanism, the paper illustrates the physics involving the shear layer transition, stretching, breakdown and turbulence generation, either qualitatively or quantitatively, in the presence of a wall for a Reynolds number of Re=1440 (based on D and the inlet free-stream velocity U_∞).     

Verification of three-dimensional laser Doppler velocimeter measurements

An experimental comparison method is proposed for the verification of mean flow and turbulence measurements obtained with a three-dimensional laser Doppler velocimeter system. Such measurements can include large errors caused by problems unique to three-dimensional systems. Direct comparisons of laser and hot cross-wire measurements obtained in two-dimensional flows, as is the common practice, will not bear out all the errors associated with three-dimensional laser systems. It is proposed here that the errors may be adequately quantified by making the direct comparisons in a weakly three-dimensional turbulent shear flow. The weak three-dimensional flow ensures high accuracy of the cross-wire data while still generating sufficiently strong secondary mean flow and Reynolds shear stresses so that all the laser measurements may be fully verified. This type of shear flow is easily generated by introducing a weak streamwise vortex into a nominally two-dimensional turbulent shear layer.     

不定原点喷嘴射流中流体物理的实验研究(英文)

为了能够更好地了解不定源喷嘴(indeterminate origin nozzle)射流中的物理过程,本文应用平面激光诱导荧光技术对一个大尺度的水射流进行了实验研究。流场显示的实验结果表明不定源喷嘴在射流的剪切层引入了蘑菇形反向旋转的涡对。这些涡的矢量方向与射流方向相同或相反,被称为流向涡(streamwise vortex)。由于射流中存在开尔文-亥姆霍兹不稳定,每当一个横向涡(spanwise vortex,即涡的矢量方向与射流方向垂直)从喷嘴脱流时会产生瞬时的低压,该瞬时低压促使向内发展的流向涡对在喷嘴的凹槽处生成。这些涡对在向下游流动的过程中会重组并在喷嘴的尖峰面生成向外发展的涡对。这些流向涡极大地影响了射流的发展。流向涡与横向涡的相互作用促使射流更早地发展成为湍流。由于流向涡同时也在射流中引入了径向的剪切流动,因此导致了更多的湍流生成从而增强了射流与周围流体的混合。     

Transition from vortex to wall driven turbulence production in the Taylor–Couette system with a rotating inner cylinder

Axisymmetrically stable turbulent Taylor vortices between two concentric cylinders are studied with respect to the transition from vortex to wall driven turbulent production. The outer cylinder is stationary and the inner cylinder rotates. A low Reynolds number turbulence model using the k‐ω formulation, facilitates an analysis of the velocity gradients in the Taylor–Couette flow. For a fixed inner radius, three radius ratios 0.734, 0.941 and 0.985 are employed to identify the Reynolds number range at which this transition occurs. At relatively low Reynolds numbers, turbulent production is shown to be dominated by the outflowing boundary of the Taylor vortex. As the Reynolds number increases, shear driven turbulence (due to the rotating cylinder) becomes the dominating factor. For relatively small gaps turbulent flow is shown to occur at Taylor numbers lower than previously reported. Copyright © 2002 John Wiley & Sons, Ltd.     

不定原点喷嘴射流中流体物理的实验研究

为了能够更好地了解不定源喷嘴（indeterminate origin nozzle）射流中的物理过程，本文应用平面激光诱导荧光技术对一个大尺度的水射流进行了实验研究。流场显示的实验结果表明不定源喷嘴在射流的剪切层引入了蘑菇形反向旋转的涡对。这些涡的矢量方向与射流方向相同或相反，被称为流向涡（streamwise vortex）。由于射流中存在开尔文一亥姆霍兹不稳定，每当一个横向涡（spanwisevortex，即涡的矢量方向与射流方向垂直）从喷嘴脱流时会产生瞬时的低压，该瞬时低压促使向内发展的流向涡对在喷嘴的凹槽处生成。这些涡对在向下游流动的过程中会重组并在喷嘴的尖峰面生成向外发展的涡对。这些流向涡极大地影响了射流的发展。流向涡与横向涡的相互作用促使射流更早地发展成为湍流。由于流向涡同时也在射流中引入了径向的剪切流动，因此导致了更多的湍流生成从而增强了射流与周围流体的混合。