Starting Flow and Heat Transfer Downstream of a Backward-Facing Step

1991 ◽  
Vol 113 (3) ◽  
pp. 583-589 ◽  
Author(s):  
F. K. Tsou ◽  
Shih-Jiun Chen ◽  
Win Aung

Experiments are performed to study the starting process of heat transfer downstream of a backward-facing step. A Ludwieg tube wind tunnel is employed to produce the incompressible flow, which accelerates from a zero velocity to a steady state value with an accelerating period of 7 ms and a steady-state period of 12 ms. Hot-wire anemometry and heat flux gages are used to measure the flow and heat transfer history, respectively. The onset of transition in the free shear layer shows that the disturbance originates from the top corner of the step, then propagating to the free stream. The velocity and turbulence profiles in the free shear layer reach steady-state values after the leading edge disturbance traverses to the measurement locations. In regions upstream and far downstream of the step, heat flux history data suggest the transformation of the flow from laminar to transitional and finally to turbulent flow. Hot-wire anemometry measurements indicate high-frequency turbulence with a short characteristic time. In the recirculating region, however, a longer characteristic time is observed because of the existence of large-scale eddies. The dimensionless reattachment length (xr/H) is shown to increase with time from the bottom corner (xr/H = 0) in the laminar regime to a maximum value of 13.6 in the transitional regime, and decreases to a constant values of 7.6 in the turbulent regime. The steady-state flow field and heat transfer compare favorably with existing data obtained using steady-state techniques.

Author(s):  
Alexander V. Mirzamoghadam ◽  
Zhenhua Xiao

Flow and heat transfer in the row-1 upstream rotor-stator disc cavity of a large 3600-rpm industrial gas turbine was investigated using an integrated approach. A 2D axisymmetric transient thermal analysis using aero engine-based correlations was performed to predict the steady state metal temperatures and hot running seal clearances at ISO rated power condition. The cooling mass flow and the flow pattern assumption for the thermal model were obtained from the steady state 2D axisymmetric CFD study. The CFD model with wall heat transfer was validated using cavity steady state air temperatures and static pressures measured at inlet to the labyrinth seal and four cavity radial positions in an engine test which included the mean annulus static pressure at hub radius. The predicted wall temperature distribution from the matched thermal model was used in the CFD model by incorporating wall temperature curve-fit polynomial functions. Results indicate that although the high rim seal effectiveness prevents ingestion from entering the cavity, the disc pumping flow draws air from within the cavity to satisfy entrainment leading to an inflow along the stator. The supplied cooling flow exceeds the minimum sealing flow predicted from both the rotational Reynolds number-based correlation and the annulus Reynolds number correlation. However, the minimum disc pumping flow was found to be based on a modified entrainment expression with a turbulent flow parameter of 0.08. The predicted coefficient of discharge (Cd) of the industrial labyrinth seal from CFD was confirmed by modifying the carry-over effect of a correlation reported recently in the literature. Moreover, the relative effects of seal windage and heat transfer were obtained and it was found that contrary to what was expected, the universal windage correlation was more applicable than the aero engine-based labyrinth seal windage correlation. The CFD predicted disc heat flux profile showed reasonably good agreement with the free disc calculated heat flux. The irregular cavity shape and high rotational Reynolds number (in the order of 7×107) leads to entrance effects that produce a thicker turbulent boundary layer profile compared to that predicted by the 1/7 power velocity profile assumption.


2000 ◽  
Vol 124 (1) ◽  
pp. 125-132 ◽  
Author(s):  
A. V. Mirzamoghadam ◽  
Z. Xiao

Flow and heat transfer in the row-1 upstream rotor-stator disk cavity of a large 3600-rpm industrial gas turbine was investigated using an integrated approach. A two dimensional axisymmetric transient thermal analysis using aeroengine-based correlations was performed to predict the steady-state metal temperatures and hot running seal clearances at ISO rated power condition. The cooling mass flow and the flow pattern assumption for the thermal model were obtained from the steady-state two dimensional axisymmetric CFD study. The CFD model with wall heat transfer was validated using cavity steady-state air temperatures and static pressures measured at inlet to the labyrinth seal and four cavity radial positions in an engine test which included the mean annulus static pressure at hub radius. The predicted wall temperature distribution from the matched thermal model was used in the CFD model by incorporating wall temperature curve-fit polynomial functions. Results indicate that although the high rim seal effectiveness prevents ingestion from entering the cavity, the disk pumping flow draws air from within the cavity to satisfy entrainment leading to an inflow along the stator. The supplied cooling flow exceeds the minimum sealing flow predicted from both the rotational Reynolds-number-based correlation and the annulus Reynolds number correlation. However, the minimum disk pumping flow was found to be based on a modified entrainment expression with a turbulent flow parameter of 0.08. The predicted coefficient of discharge (Cd) of the industrial labyrinth seal from CFD was confirmed by modifying the carryover effect of a correlation reported recently in the literature. Moreover, the relative effects of seal windage and heat transfer were obtained and it was found that contrary to what was expected, the universal windage correlation was more applicable than the aeroengine-based labyrinth seal windage correlation. The CFD predicted disk heat flux profile showed reasonably good agreement with the free disk calculated heat flux. The irregular cavity shape and high rotational Reynolds number (in the order of 7×107) leads to entrance effects that produce a thicker turbulent boundary layer profile compared to that predicted by the 1/7 power velocity profile assumption.


2020 ◽  
Vol 2020 ◽  
pp. 1-13
Author(s):  
Zhanwei Liu ◽  
Xinyu Li ◽  
Tenglong Cong ◽  
Rui Zhang ◽  
Lingyun Zheng ◽  
...  

The prediction of flow and heat transfer characteristics of liquid sodium with CFD technology is of significant importance for the design and safety analysis of sodium-cooled fast reactor. The accuracies and uncertainties of the CFD models should be evaluated to improve the confidence of the numerical results. In this work, the uncertainties from the turbulent model, boundary conditions, and physical properties for the flow and heat transfer of liquid sodium were evaluated against the experimental data. The results of uncertainty quantization show that the maximum uncertainties of the Nusselt number and friction coefficient occurred in the transition zone from the inlet to the fully developed region in the circular tube, while they occurred near the reattachment point in the backward-facing step. Furthermore, in backward-facing step flow, the maximum uncertainty of temperature migrated from the heating wall to the geometric center of the channel, while the maximum uncertainty of velocity occurred near the vortex zone. The results of sensitivity analysis illustrate that the Nusselt number was negatively correlated with the thermal conductivity and turbulent Prandtl number, while the friction coefficient was positively correlated with the density and Von Karman constant. This work can be a reference to evaluate the accuracy of the standard k-ε model in predicting the flow and heat transfer characteristics of liquid sodium.


Computation ◽  
2021 ◽  
Vol 9 (6) ◽  
pp. 65
Author(s):  
Aditya Dewanto Hartono ◽  
Kyuro Sasaki ◽  
Yuichi Sugai ◽  
Ronald Nguele

The present work highlights the capacity of disparate lattice Boltzmann strategies in simulating natural convection and heat transfer phenomena during the unsteady period of the flow. Within the framework of Bhatnagar-Gross-Krook collision operator, diverse lattice Boltzmann schemes emerged from two different embodiments of discrete Boltzmann expression and three distinct forcing models. Subsequently, computational performance of disparate lattice Boltzmann strategies was tested upon two different thermo-hydrodynamics configurations, namely the natural convection in a differentially-heated cavity and the Rayleigh-Bènard convection. For the purposes of exhibition and validation, the steady-state conditions of both physical systems were compared with the established numerical results from the classical computational techniques. Excellent agreements were observed for both thermo-hydrodynamics cases. Numerical results of both physical systems demonstrate the existence of considerable discrepancy in the computational characteristics of different lattice Boltzmann strategies during the unsteady period of the simulation. The corresponding disparity diminished gradually as the simulation proceeded towards a steady-state condition, where the computational profiles became almost equivalent. Variation in the discrete lattice Boltzmann expressions was identified as the primary factor that engenders the prevailed heterogeneity in the computational behaviour. Meanwhile, the contribution of distinct forcing models to the emergence of such diversity was found to be inconsequential. The findings of the present study contribute to the ventures to alleviate contemporary issues regarding proper selection of lattice Boltzmann schemes in modelling fluid flow and heat transfer phenomena.


Author(s):  
Ratan Kumar Chanda ◽  
Mohammad Sanjeed Hasan ◽  
Md. Mahmud Alam ◽  
Rabindra Nath Mondal

2021 ◽  
pp. 875608792110258
Author(s):  
Azhar Ali ◽  
Dil Nawaz Khan Marwat ◽  
Aamir Ali

Flows and heat transfer over stretching/shrinking and porous surfaces are studied in this paper. Unusual and generalized similarity transformations are used for simplifying governing equations. Current model includes all previous cases of stretched/shrunk flows with thermal effects discussed so far. Moreover, we present three different cases of thermal behavior (i) prescribed surface temperature (ii) Variable/uniform convective heat transfer at plat surface and (iii) prescribed variable/uniform heat flux. Stretching/shrinking velocity Uw(x), porosity [Formula: see text], heat transfer [Formula: see text], heat flux [Formula: see text] and convective heat transfer at surface are axial coordinate dependent. Boundary layer equations and boundary conditions are transformed into nonlinear ODEs by introducing unusual and generalized similarity transformations for the variables. These simplified equations are solved numerically. Final ODEs represent suction/injection, stretching/shrinking, temperature, heat flux, convection effects and specific heat. This current problem encompasses all previous models as special cases which come under the scope of above statement (title). The results of classical models are scoped out as a special case by assigning proper values to the parameters. Numerical result shows that the dual solutions can be found for different possible values of the shrinking parameter. A stability analysis is accomplished and apprehended in order to establish a criterion for determining linearly stable and physically compatible solutions. The significant features and diversity of the modeled equations are scrutinized by recovering the previous problems of fluid flow and heat transfer from a uniformly heated sheet of variable (uniform) thickness with variable (uniform) stretching/shrinking and injection/suction velocities.


Sign in / Sign up

Export Citation Format

Share Document