10th joint ERCOFTAC (SIG-15) -IAHR-QNET/CFD Workshop on Refined Turbulence Modelling

 
 

CASE 10.1
Contra-rotating jets

Database: Experiments by S. Béharelle, J. Delville and J.-P. Bonnet
 
 


 

Fig. 1: Experimental setup and coordinate system
 

Fig. 2: Schematic streamwise (U) and azimuthal (W) profiles
 
 

Fig. 3: Wake/mixing layer interaction

















Flow description:

As shown in Fig. 1, the flow consists of two coaxial jets with opposite rotation.  The diameter of the inner jet is r0=120 mm and the outer jet is located between r and R=200 mm.  The flow is confined in a duct of the same diameter R as the outer jet.
As described in Fig. 2, the mean velocity U in x-direction has a peak at the centre, but the wake in r0 is surrounded by a region of constant U0 velocity (U0=26.5 m s-1).  In this region (approximately located between r1=70 mm and r2=170 mm), an azimuthal mixing layer develops, with =10 ms-1.

The region of interest is thus located between r1=70 mm and r2=170 mm, where an interaction between the wake of the inner nozzle and the azimuthal mixing layer takes place, as shown in Fig. 3.

Fig. 3 also describes the local axis system (x,y,z) used in the database, in which y=0 corresponds to r=r0 .

In this flow, a transition between a wake and a mixing layer occurs as the flow develops is x-direction.  In particular, the Reynolds stress tensor experiences a drastic change in its structure, from a dominant <u2>-component with two maxima to a dominant <w2>-component with a single maximum.  This behaviour is linked to the fact that the initial streamwise wake (U-velocity) is suddenly subjected to a transverse shear (W-velocity), which mainly produces turbulence on the <w2>-component.  Note that contrary to usual mixing layers, the direction in which the flow develops (x-direction) is normal to the mean velocity direction that creates the shear (W-velocity).  The flow thus tends to a spatially developing mixing layer very similar to a temporal mixing layer.
 

Calculation domain and definition of the flow:


Fig. 4: Definition of the calculation domain








This case is axisymmetric, but all the components are to be calculated (U, V, W, and, for Reynolds stress models, <u2>, <v2>, <w2>, <uv>, <uw>, <vw>).

The inlet is not located at x=0, but at x0/r0=10/120=0.08333.
The inlet profiles for the mean velocities and the Reynolds stresses are provided in the database (files *.0010).
 

The last experimental profiles are given at x/r0=1000/120=8.333, so the outlet must be placed well after that location (at least L/r0=10 is recommended).

The experiments are performed in air (kinematic viscosity=15x10-6).  The Reynolds number based on U and r0 is Re=212,000.
 

Database:

(All the quantities are made dimensionless using r0 =120 mm and U0=26.5 m s-1)

- Profiles:
The database contains profiles of the mean velocities and the Reynolds stresses at 15 locations, ranging from x=50 mm to x=1000 mm.
The files consist of two columns, for y and the variable under consideration.  They are named VAR.X, where VAR=U, W, uu, vv, ww, uv, uw or vw and X=0050, 0100, 0150, 0200, 0250, 0300, 0350, 0400, 0450, 0500, 0600, 0700, 0800, 0900 or 1000 (these numbers correspond to the location on the x-axis is mm).

Download the Mean Velocities (16k)
Download the Reynolds Stresses (33k)

- Inlet:
The files VAR.0010 are to be used as inlet boundary conditions.
V is very weak everywhere in the domain (V/U0<1%), so V=0 is to be imposed at the inlet.
The profile of the dissipation (file Eps.0010) has been evaluated from the turbulent kinetic energy profile using a constant mixing length.  This evaluation is very close (but smoother) to the one using the eddy-viscosity assumption with =0.09.

Download the Inlet profiles (4.0k)

- Budgets:
Some terms of the budgets of the Reynolds stresses are available.  The files consist of 4 columns for Y, convection, production and turbulent diffusion, respectively.  They are named budget.VAR.X, where VAR=k, uu, vv, ww, uv, uw, vw and X=0050, 0100, 0150, 0200, 0250, 0300, 0350, 0400, 0450, 0500, 0600, 0700, 0800 or 0900.

Download the Budgets (86k)
 
 

References:

Béharelle S., Delville J. and Bonnet J.-P.(516k), On the three dimensional evolution of a wake subjected to cross shear, Turbulence and Shear Flow 11, 8-11 sept 1997, Grenoble, France.

Béharelle S., Nayeri C., Delville J., Bonnet J.-P., Fiedler H.E., Influence of the transverse shear on the development of wake flows, Advances in Turbulence VI, S. Gravilakis et al. editors, 511-512.
 
 

Required results:

Participants are invited to submit results in the same form as the database, i.e.:

(Recall that all the quantities must be scaled by r0 =120 mm and U0=26.5 m s-1.)

- profiles in two-column files named VAR.X, where VAR=U, W, uu, vv, ww, uv, uw or vw and X=0050, 0100, 0150, 0200, 0250, 0300, 0350, 0400, 0450, 0500, 0600, 0700, 0800, 0900 or 1000.

- budgets of <uu>, <vv>, <ww>, <uv>, <uw>, <vw>  and k.  Irrelevant files are to be omitted (e.g., for a k-epsilon model, only k budgets are required).  The files are to be named budget.VAR.X, where VAR stands for uu, vv, ww, uv, uw, vw, k and X stands again for the location.

Each file starts with 8 comment lines beginning by a "#", consisting of:

#Case
#Name
#Affiliation
#Code/Numerical method
#Mesh type
#Turbulence model
#Any other useful information
#Column headings for the data

Examples:

File "U.0050":

#Case10.1
#Rémi Manceau
#Laboratoire d'études aérodynamiques, université de Poitiers
#Finite volumes
#Block structured
#Reynolds stress model (Speziale, Sarkar, Gatski, 1991)
#Without non-linear slow term, fine mesh
# y/r0     U/U0
 

File "budget.uu.0100"

#Case10.1
#Rémi Manceau
#Laboratoire d'études aérodynamiques, université de Poitiers
#Finite volumes
#Block structured
#Reynolds stress model (Speziale, Sarkar, Gatski, 1991)
#Without non-linear slow term, fine mesh
# Y/r0    convec   produc   t-diff    p-strain   p-diff    v-diff    dissip
 
 

File "budget.k.0900"

#Case10.1
#Rémi Manceau
#Laboratoire d'études aérodynamiques, université de Poitiers
#Finite volumes
#Block structured
#Reynolds stress model (Speziale, Sarkar, Gatski, 1991)
#Without non-linear slow term, fine mesh
#Y/r0    convec   produc   t-diff    p-strain    p-diff    v-diff    dissip
 

(N.B.: p-strain is of course supposed to be zero in the files budget.k.X.)
 

Moreover, participants are invited to submit a short description of their computations (mesh, numerical method, turbulence modelling, etc.)
Each participant will have 5 minutes to present these details (no results presentation).
 

To know how to submit the results, click here.
 
 


For additional information, please contact R. Manceau (remi.manceau@lea.univ-poitiers.fr)


Back to the Announcement Page

Back to the main page