Numerical predictions of channel ows with uid injection through a porous wall are performed by solving the time-dependent Navier-Stokes equations using a Reynolds-stress turbulent model. In uence of the turbulence injected uid is investigated. Numerical results with experimental data indicate that the ows evolve signi cantly vs the distance from the front wall such that different regimes of ow development can be observed. In the rst regime the velocity eld is developed in accordance with the laminar theory. The second regime is characterized by the development of turbulence, which occurs at different locations in the channel because of the presence of impermeable and permeable walls, and by the transition process of the mean axial velocity when a critical turbulence threshold is attained. Computed results are compared with existing experimental data including axial mean velocity pro les and full turbulent stresses. As a result for the simulations, the Reynolds-stress model predicts the mean velocity pro les, the transition process, and the turbulent stresses, in good agreement with experimental data.
Nomenclature= total speci c energy, m 2 /s 2 , J/kg H = total speci c enthalpy, h C u i u i =2, m 2 /s 2 h = speci c enthalpy, m 2 /s 2 J i j = tensor of diffusion for the Reynolds stress ¿ i j k = speci c turbulent kinetic energy, ¿ ii =2, m 2 /s 2 L = channel length, m M = Mach number m = injection mass ux, kg/(m 2 ¢ s) n i = normal to the wall P i j = production rate of ¿ i j caused by mean shear P rt = turbulent Prandtl number p = static pressure, Pa q i = total heat ux vector, W/m 2 R s = injection Reynolds number, ½ s u s ±=¹ s R t = turbulent Reynolds number k 2 =º² R u = universal gas constant S i j = strain-rate tensor s = speci c entropy, J/(kg ¢ K)velocity vector, m/s u m = bulk velocity, m/s u s = injection velocity, m/s u ¿ = friction velocity, m/s x i = Cartesian coordinate, m x C i = dimensionless distance from walls, x i u ¿ =º ® = coef cient for planar or axisymmetric geometryratio of speci c heats ± = channel height, m ± i j = Kronecker tensor ² = dissipation rate, m 2 /s 3 ² i j k = permutation tensor · = thermal conductivity, W/(m ¢ K) ¹ = dynamic viscosity, kg/(m ¢ s) º = kinematic viscosity, m/s 2 ½ = density, kg/m 3 6 i j = total stress tensor ¾ i j = viscous stress tensor ¾ s = surface-generatedpseudoturbulence,. ] u 00 2 u 00turbulent stress tensor, ] u 00 i u 00 j 8 i j = pressure-strain uctuations, p 0 S 00 i j ! i = vorticity tensor, ² i jk @u k =@ x j , (1/s) Subscripts m = bulk mean quantity s = condition at injection surface w = wall Superscripts N = Reynolds averaged of variable Q = Favre averaged of variable 0 = Reynolds turbulent uctuating value of variable 00 = Favre turbulent uctuating value of variable