In predicting the fully developed turbulent mixed flow of forced and natural convection, the one-equation model of turbulence for low Reynolds number flows is adopted. Finite-element solutions are carried out in the upward turbulent flow between vertical parallel plates at different wall temperatures. In this condition, the aiding and opposing flows arise simultaneously on the heated and cooled sides, respectively, and a fully developed condition is established. The solutions with and without the effect of buoyancy force are predicted and the comparisons are discussed. In these two situations, the mean velocity and temperature profiles, the locations of the maximum velocity, eddy diffusivities of momentum, Nusselt number, friction factor, and turbulence kinetic energy are presented. The substantial effects of the buoyancy force are also confirmed.eddy diffusivity for heat 8 = global source vector b = distance between parallel plates bj = element of global source vector b + = dimensionless distance between parallel plates = bu*/v c p = specific heat at constant pressure f c -friction factor = 2 t c \/p (u >? f H = friction factor = 2 i; H \/p <« >Ĝ + = buoyancy parameter = gf$(T H -T c )v/u* 3 Gr -Grashof number based on temperature difference = *gp(T H -T c )b 3 /v 2 g = gravitational acceleration K = turbulence kinetic energy K + = dimensionless turbulence kinetic energy = K/u* 2 k = von Karman constant / = element length / 1Z) = mixing length in dissipation ID -dimensionless mixing length in dissipation = /i jD w*/v £ lv = mixing length in diffusivity £i~v = dimensionless mixing length in diffusivity = £ lv u*/v N f = shape function Nu c = Nusselt number = 4|#|<5/«r> c -rjM Nu H = Nusselt number = 4\q\(b -d)/(T H -P = turbulence production P u = viscous production P 0 = buoyancy production Pr = Prandtl number p = real pressure p s = static pressure p + = dimensionless pressure = /pw* 2 q = heat flux -Reynolds number = 2£/v = local turbulent Reynolds number = dimensionless local turbulent Reynolds number = temperature = dimensionless temperature normalized by the cooled-wall parameter = (T -T c )pc p u*/\q\ = dimensionless temperature normalized by the heated-wall parameter = (T H -T)pc p Uff/\q\ = time-averaged velocity in the x direction = dimensionless velocity = u/u* = friction velocity = ^/\r c \/p = velocity in y direction = coordinate in the flow direction = dimensionless coordinate in the flow direction = xu*/v = distance from the cooled wall = dimensionless distance from the cooled = correction factor defined by Eq. (40) = thermal expansion coefficient = distance from the cooled wall to the maximum velocity location = dimensionless distance = du*/v = turbulent dissipation = dimensionless temperature = (r -T c )j(T H -T c ) = thermal conductivity = kinematic viscosity = eddy diffusivity for momentum = dimensionless eddy diffusivity for momentum = V T /V = density = turbulent Prandtl number = shear stress Subscripts c H m o = at the cooled wall = at the heated wall = momentum = forced convection < )...