Prediction of Turbulent Source Flow Between Corotating Disks With an Anisotropic Two-Equation Turbulence Model

1988 ◽  
Vol 110 (2) ◽  
pp. 187-194 ◽  
Author(s):  
S. A. Shirazi ◽  
C. R. Truman

An anisotropic form of a low-Reynolds-number two-equation turbulence model has been implemented in a numerical solution for incompressible turbulent flow between corotating parallel disks. Transport equations for turbulent kinetic energy and dissipation rate were solved simultaneously with the governing equations for the mean-flow variables. Comparisons with earlier mixing-length predictions and with measurements are presented. Good agreement between the present predictions and the measurements of velocity components and turbulent kinetic energy was obtained. The low-Reynolds-number two-equation model was found to model adequately the near-wall region as well as the effects of rotation and streamline divergence, which required ad hoc assumptions in the mixing-length model.

1991 ◽  
Vol 230 ◽  
pp. 1-44 ◽  
Author(s):  
Lincoln P. Erm ◽  
Peter N. Joubert

An investigation was undertaken to improve our understanding of low-Reynolds-number turbulent boundary layers flowing over a smooth flat surface in nominally zero pressure gradients. In practice, such flows generally occur in close proximity to a tripping device and, though it was known that the flows are affected by the actual low value of the Reynolds number, it was realized that they may also be affected by the type of tripping device used and variations in free-stream velocity for a given device. Consequently, the experimental programme was devised to investigate systematically the effects of each of these three factors independently. Three different types of device were chosen: a wire, distributed grit and cylindrical pins. Mean-flow, broadband-turbulence and spectral measurements were taken, mostly for values of Rθ varying between about 715 and about 2810. It was found that the mean-flow and broadband-turbulence data showed variations with Rθ, as expected. Spectra were plotted using scaling given by Perry, Henbest & Chong (1986) and were compared with their models which were developed for high-Reynolds-number flows. For the turbulent wall region, spectra showed reasonably good agreement with their model. For the fully turbulent region, spectra did show some appreciable deviations from their model, owing to low-Reynolds-number effects. Mean-flow profiles, broadband-turbulence profiles and spectra were found to be affected very little by the type of device used for Rθ ≈ 1020 and above, indicating an absence of dependence on flow history for this Rθ range. These types of measurements were also compared at both Rθ ≈ 1020 and Rθ ≈ 2175 to see if they were dependent on how Rθ was formed (i.e. the combination of velocity and momentum thickness used to determine Rθ). There were noticeable differences for Rθ ≈ 1020, but these differences were only convincing for the pins, and there was a general overall improvement in agreement for Rθ ≈ 2175.


2019 ◽  
Vol 866 ◽  
pp. 897-928 ◽  
Author(s):  
P. Orlandi

Data available in the literature from direct numerical simulations of two-dimensional turbulent channels by Lee & Moser (J. Fluid Mech., vol. 774, 2015, pp. 395–415), Bernardini et al. (J. Fluid Mech., 742, 2014, pp. 171–191), Yamamoto & Tsuji (Phys. Rev. Fluids, vol. 3, 2018, 012062) and Orlandi et al. (J. Fluid Mech., 770, 2015, pp. 424–441) in a large range of Reynolds number have been used to find that $S^{\ast }$ the ratio between the eddy turnover time ($q^{2}/\unicode[STIX]{x1D716}$, with $q^{2}$ being twice the turbulent kinetic energy and $\unicode[STIX]{x1D716}$ the isotropic rate of dissipation) and the time scale of the mean deformation ($1/S$), scales very well with the Reynolds number in the wall region. The good scaling is due to the eddy turnover time, although the turbulent kinetic energy and the rate of isotropic dissipation show a Reynolds dependence near the wall; $S^{\ast }$, as well as $-\langle Q\rangle =\langle s_{ij}s_{ji}\rangle -\langle \unicode[STIX]{x1D714}_{i}\unicode[STIX]{x1D714}_{i}/2\rangle$ are linked to the flow structures, and also the latter quantity presents a good scaling near the wall. It has been found that the maximum of turbulent kinetic energy production $P_{k}$ occurs in the layer with $-\langle Q\rangle \approx 0$, that is, where the unstable sheet-like structures roll-up to become rods. The decomposition of $P_{k}$ in the contribution of elongational and compressive strain demonstrates that the two contributions present a good scaling. However, the good scaling holds when the wall and the outer structures are separated. The same statistics have been evaluated by direct simulations of turbulent flows in the presence of different types of corrugations on both walls. The flow physics in the layer near the plane of the crests is strongly linked to the shape of the surface and it has been demonstrated that the $u_{2}$ (normal to the wall) fluctuations are responsible for the modification of the flow structures, for the increase of the resistance and of the turbulent kinetic energy production.


2020 ◽  
Vol 24 (5 Part A) ◽  
pp. 2729-2741
Author(s):  
Zhenchuan Wang ◽  
Guoli Qi ◽  
Meijun Li

The turbulence model fails in supercritical fluid-flow and heat transfer simulation, owing to the drastic change of thermal properties. The inappropriate buoyancy effect model and the improper turbulent Prandtl number model are several of these factors lead to the original low-Reynolds number turbulence model unable to predict the wall temperature for vertically heated tubes under the deteriorate heat transfer conditions. This paper proposed a simplified improved method to modify the turbulence model, using the generalized gradient diffusion hypothesis approximation model for the production term of the turbulent kinetic energy due to the buoyancy effect, using a turbulence Prandtl number model for the turbulent thermal diffusivity instead of the constant number. A better agreement was accomplished by the improved turbulence model compared with the experimental data. The main reason for the over-predicted wall temperature by the original turbulence model is the misuse of the buoyancy effect model. In the improved model, the production term of the turbulent kinetic energy is much higher than the results calculated by the original turbulence model, especially in the boundary-layer. A more accurate model for the production term of the turbulent kinetic energy is the main direction of further modification for the low Reynolds number turbulence model.


2020 ◽  
Vol 78 (11) ◽  
pp. 674-695
Author(s):  
Yi Yang ◽  
Guanglin Qiang ◽  
Zhen Chen ◽  
Zhengqi Gu ◽  
Yong Zhang

Sign in / Sign up

Export Citation Format

Share Document