scholarly journals Simulation of turbulent flow and heat transfer over a backward facing step with ribs turbulators

2011 ◽  
Vol 15 (1) ◽  
pp. 245-255 ◽  
Author(s):  
Khudheyer Mushatet

Simulation is presented for a backward facing step flow and heat transfer inside a channel with ribs turbulators. The problem was investigated for Reynolds numbers up to 32000. The effect of a step height, the number of ribs and the rib thickness on the flow and thermal field were investigated. The computed results are presented as streamlines counters, velocity vectors and graphs of Nusselt number and turbulent kinetic energy variation. A control volume method employing a staggered grid techniques was imposed to discretize the governing continuity, full Navier Stockes and energy equations. A computer program using a SIMPLE algorithm was developed to handle the considered problem. The effect of turbulence was modeled by using a k-? model with its wall function formulas. The obtained results show that the strength and size of the re-circulation zones behind the step are increased with the increase of contraction ratio(i.e. with the increase of a step height). The size of recirculation regions and the reattachment length after the ribs are decreased with increasing of the contraction ratio. Also the results show that the Reynolds number and contraction ratio have a significant effect on the variation of turbulent kinetic energy and Nusselt number.

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.


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.


Author(s):  
Sogol Pirbastami ◽  
Samir Moujaes

A Computational Fluid Dynamics (CFD) study of heat enhancement in helically grooved tubes was carried out by using a 3-dimensional simulation with the STARCCM+ simulation package software. The k-ε model selected for turbulent flow simulation and the governing equations were solved by using the finite volume method. Geometric models of the current study include 3 rectangular grooved tubes with different groove width (w) and depth (e) which varies from 0.2 mm to 0.6 mm for the same tube length of 2.0m and diameter of 7.1 mm. The simulations were performed in the Reynolds number (Re) range of 4000–10000 with a uniform wall heat flux of 3150 w/m2 applied as a boundary condition on the surface of each tube. The purpose of this research is to investigate the effect of different groove dimensions on the thermal performance and pressure drop of water inside the grooved tubes and clarify the structural nature of the flow in regards to flow swirl and turbulent kinetic energy distributions. It was found that the highest performance belongs to the groove with these dimensions (w = 0.2 mm and e = 0.2 mm) which was considered for further study. Then, for these same groove dimensions four pitch size to tube diameter (p/D) ratios ranging from 1 to 18 were simulated for the same 2.0 m length tube. The results for Nusselt number (Nu) and friction factor (f) showed that by increasing the (p/D) ratio both the Nu numbers and the friction factors (f) values decrease. With a smaller pitch length (p) the turbulence intensity generated by the internal groove was also found to increase. The physical behavior of the turbulent flow and heat transfer characteristics were observed by contour plots which showed an increasing swirl flow and turbulent kinetic energy as p/D decreases. With an increase of the Nu number for smaller p/D ratio, a penalty of a higher pressure drop was obtained. The results were validated with a previous experimental work and the average error between the experimental and CFD Nu numbers and f were 13% and 8% respectively. A higher level of turbulent kinetic energy is observed near the grooves, as compared to the smooth areas of the pipe surface away from the grooves, which are expected to lead to higher levels of heat transfer. The effect of pitch length (p) on the flow pattern were plotted by streamlines along the tubes, by decreasing the pitch size (p/D ratio) an increase in the swirl is noticed as evidenced by the plots of the path lines. Finally, empirical correlations for Nusselt number and friction factor were provided as a function of p/D and Re number. This study indicates that the incorporation of the internal groove, of particular dimensions, can lead to an improvement of performance in heat exchanger devices. A limited variation of the groove dimensions was conducted and it was found that the values of Nu and f do not improve with an increase of (w) nor with that of (e) from 0.2–0.6 mm.


2020 ◽  
Vol 184 ◽  
pp. 01027
Author(s):  
B Ch Nookaraju

Computational investigation of steady, two-dimensional heat transfer attributes for forced convective chaotic discharge in a vertical channel of cluster of heated rectangular sections is performed. The discharge is deemed to be periodic fully developed so that the issue is determined for two extending zone and explanation is developed to more number of sections. This structure reproduces the driven convective cooling of a cluster of engraved circuit panels confronted in computerize belongings. Two mathematical statements for k- ℇ model is used for modeling for the turbulence and the finite volume methodology is used. Computations are performed for Reynolds numbers ranging from 6000-12000, Prandtl number of 0.7 and various geometric parameters characterizing the problem. As Reynolds number steps up the Nusselt Number increases. Re-circulations undermine the local Nusselt number when matched with comparing variation from a identical plate. The velocity contours, temperature distributions, variation of turbulent kinetic energy and kinetic energy dissipation rates in a vertical channel is found. With the blocks in the cluster, pressure fall is higher in resemblance to plane duct.


Author(s):  
Domenico Borello ◽  
Alessandro Salvagni

The need to study flow and heat transfer in turbine blade cooling design calls to develop appropriate modelling approaches able to return accurate predictions at a reduced computational costs. Here we propose and scrutinize a quadratic version of the well-known k-ε-ζ-f RANS turbulence models, aiming at sensitizing the model to the effect of rotation in configurations mimicking the flow in turbine internal cooling. Starting from the evidence that rotation modified turbulent flow through a turbulence suppression (enhancement) on the stabilized (destabilized) surface, we modified the Cμ coefficient present in the formulation of turbulent viscosity introducing a dependence on the strain and vorticity tensors, the latter explicitly including solid body rotation. The proposed model was tested on plane channel and square-sectioned duct flows, and then used for simulating a rib-duct rotating channel. Results are assessed against DNS literature data and properly developed LES computations, by examining flow variables, heat transfer and turbulence budgets. We demonstrate that, as for the channel flows, the proposed quadratic model is able to accurately reproduce velocity, temperature and turbulent variables at various angular velocity regimes. In the duct flow the flow is subjected to the mutual influence vorticity induced by rotation and turbulence anisotropy developing close the walls. In particular, the non-linear rotation-sensitized model is able to reproduce the near-wall turbulent kinetic energy distribution close to the suction side, returning a zero value in the mid-span and a small peak close to the wall on the suction side. Turbulent kinetic energy and temperature budgets analysis demonstrates the capabilities of the model in describing all the terms in the equations. Also if some tuning of the model is required, these analysis showed very encouraging results. In fact if the basic mechanisms of turbulence and heat transfer are properly predicted, then it can be expected that the model can be successfully applied to a set of different cases. For such reason, the model was applied to the analysis of flow and heat transfer in a rotating ribduct with reasonably results.


2012 ◽  
Vol 468-471 ◽  
pp. 2014-2018
Author(s):  
Khudheyer S. Mushatet

In this paper, A numerical study has been done to predict the heat transfer of the staggered ribbed backward facing step flow with inclined impinging jet flow. The impinging jet flow was inclined towards the main cross flow and the angle of inclination is ranged from 30ْ to 90ْ . The ribs were in staggered arrangement and aligned after the slot jet in normal direction to the main cross flow. The effect of angle of inclination and contraction ratio on thermal field was studied for jet and channel Reynolds number of 20000 and 16000 respectively. The aim of the present study is to verify how adding staggered ribs with inclined jet flow to the problem of backward facing can enhance the rate of heat transfer. The governing partial differential equations of continuity, Navier-Stockes and energy was discretised on non-uniform staggered grid by using finite volume method. The discretised algebraic equations were solved by using a built home computer program based on simple algorithm. The conducted results show that the considered flow geometry increased noticeably the rate of heat transfer for the studied contraction ratios and angle of inclination. It was observed that the rate of heat transfer is decreased as angle of inclination increases.


2020 ◽  
Vol 98 (3) ◽  
pp. 297-302 ◽  
Author(s):  
Ali Asghar Sedighi ◽  
Zeynab Deldoost ◽  
Bahram Mahjoob Karambasti

The flow and heat transfer of Al2O3–water nanofluid in a channel partially filled with porous media is investigated numerically. The turbulence effect in the porous media is taken under consideration in this article. A simple case is simulated first to evaluate the accuracy of the results in comparison with the available data. The turbulent kinetic energy profile is investigated at a flow cross section. The results show that the maximum turbulent kinetic energy occurs in the clear fluid region in the vicinity of the porous media region. The turbulent kinetic energy is a decreasing function of the porosity of the porous medium. The effect of porosity on the variation of turbulent kinetic energy decreases with the increase in the porosity of the porous medium. The turbulent kinetic energy in clear fluid and porous media regions decreases with the increase in nanofluid concentration from 0.01 to 0.03, and it increases with the increase in nanofluid concentration from 0.03 to 0.05. The temperature of the nanofluid increases with the increase in the nanofluid concentration and decrease in the porosity of porous media. It is shown that for this case, with the increase in nanofluid concentration and porosity of porous media, the skin friction coefficient increases and the Nusselt number decreases.


2021 ◽  
pp. 334-334
Author(s):  
Ali Hajimohammadi ◽  
Mehran Zargarabadi ◽  
Javad Mohammadpour

A computational study is carried out of the three-dimensional flow field and heat transfer under a turbulent intermittent circular jet impingement on a concave surface. The control-volume procedure with the SIMPLE algorithm is employed to solve the unsteady RANS (use full form) equations. The RNG k-? model is implemented to simulate turbulence due to its success in predicting similar flows. The numerical results are validated by comparing them with the experimental data. The effects of jet Reynolds number and oscillation frequency on the flow and heat transfer are evaluated. The profiles of instantaneous and time-averaged Nusselt numbers exhibit different trends in axial (x) and circumferential (s) directions. It is found that increasing frequency from 50 to 200 Hz results in considerable time-averaged Nusselt number enhancement in both axial and curvature directions. The intermittent jet at a frequency of 200 Hz enhances the total average Nusselt number by 51.4%, 40%, and 33.7% compared to the steady jet values at jet Reynolds numbers of 10000, 23000, and 40000, respectively. In addition, a correlation for the average Nusselt number is proposed depending on the Reynolds number and the Strouhal number.


Sign in / Sign up

Export Citation Format

Share Document