Abstract
This work presents numerical prediction for the turbulent flow field confined in a circular duct past a segment of gradually varying cross section. Both expanding and contracting sections are investigated. Equations of boundary-layer type are used and the linear k–ε model, in its high Reynolds form, is applied. A new correlation for treating the grid point closest to the wall is proposed. A marching-forward method is employed for sweeping the computational domain. Computations are first performed for developing and fully developed constant-area ducts in order to assess the reliability of the code. Results are then presented for contractions and diffusers, where comparisons with experimental data for air and water are carried out. Turbulence damping in contractions and its enhancement in diffusers are calculated correctly. Further, for contractions with angles of up to 21°, the use of a parabolic solver shows good agreement with experimental values for the mean and statistical quantities. For diffusers, adverse pressure gradient along the flow limits the quality of the predictions as the angle and length of diffuser increase past 5° and 10 duct radii, respectively.
The authors are thankful to CNPq, Brazil, for financial support during the course of this work.