Results for case Asymmetric plane diffuser

Code: Code_Saturne

Version: 1.4.0

Authors: Juan Uribe

Method and Numerical Options

The k-eps model is treated with scalable wall function (IDEUCH=2).

Centred scheme is used for the velocities (BLENCV = 1) and upwind for the turbulent variables (BLENCV=0).

The time scheme was unsteady IDTVAR = 0.

Models

k-eps

Mesh

The mesh is refined near the walls to have y+=1 and solve the viscous sublayer.

In v1.4.0 the cgns mesh was not read properly by the code because of the lack of colour definition for the boundaries. The mesh was imported into StarCD and then exported with colours as a .ccm file.

Description of the results files

The files contain measures at the following planes to match the experiments by Buice and Eaton.

Section code	x/H	U / Ub	uu / Ub2	uv / Ub2	vv / Ub2
-6	-5.87	x	x	x	x
03	2.59	x	x	x	x
06	5.98	x	x	x	x
13	12.75			x	x
14	13.56	x	x
16	16.14			x	x
17	16.93	x	x
19	19.53			x	x
20	20.32	x	x
23	22.91			x	x
24	23.71	x	x
26	26.3			x	x
27	27.09	x	x
29	29.69			x	x
30	30.48	x	x
33	33.07			x	x
34	33.87	x	x
40	39.85	x	x	x	x
47	46.62	x	x	x	x
53	53.39	x	x	x	x
60	60.17	x	x	x	x
67	66.94	x	x	x	x
74	73.71	x	x	x	x

• Pressure coefficient Cp = (p-px=-1.7) / (0.5 rho Ub2) along the inclined bottom wall (g_cp%.y00) and straight top wall (g_cp%.y10). Here, px=-1.7 means static pressure at x/H = -1.7.

• Skin friction coefficient cf = tauw / (0.5 rho Ub2) along the inclined bottom wall (g_cf%.y00) and the straight top wall (g_cf%.y10).

• U-component of the mean velocity scaled by Ub.

• Reynolds stresses: uu, uv, and vv scaled by Ub2

Boundary conditions

The inlet boundary conditions are taken from a precursor 1D channel flow simulation with the corresponding turbulence model.

The wall are treated as no-slip but wall function is used for k-eps

The sides are set to symmetry.

Reference Publications

Additional information (links, pictures, etc.)

• profiles at x/h =-6:
  

• profiles at x/h =3:
  

• profiles at x/h = 14:
  

• profiles at x/h =20:
  

• profiles at x/h =27:
  

• profiles at x/h =34:
  

• profiles at x/h =47:
  

• profiles at x/h = 60:
  

---