scholarly journals Discussion on improved method of turbulence model for supercritical water flow and heat transfer

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.

1994 ◽  
Vol 116 (2) ◽  
pp. 369-380 ◽  
Author(s):  
P. Tekriwal

Standard and extended k–ε turbulence closure models have been employed for three-dimensional heat transfer calculations for radially outward flow in rectangular and square cooling passages rotating in orthogonal mode. The objective of this modeling effort is to validate the numerical model in an attempt to fill the gap between model predictions and the experimental data for heat transfer in rotating systems. While the trend of heat transfer predictions by the standard k–ε turbulence model is satisfactory, the differences between the data and the predictions are approximately 30 percent or so in the case of high rotation number flow. The extended k–ε turbulence model takes an approach where an extra “source” term based on a second time scale of the turbulent kinetic energy production rate is added to the equation for the dissipation rate of turbulent kinetic energy. This yields a more effective calculation of turbulent kinetic energy as compared to the standard k–ε turbulence model in the case of high rotation number and high density ratio flow. As a result, comparison with the experimental data available in the literature shows that an improvement of up to a significant 15 percent (with respect to data) in the heat transfer coefficient predictions is achieved over the standard k–ε model in the case of high rotation number flow. Comparisons between the results of the standard k–ε model and the extended formulation are made at different rotation numbers, different Reynolds numbers, and varying temperature ratio. The results of the extended k–ε turbulence model are either as good or better than those of the standard k–ε model in all these cases of parametric study. Thus, the extended k–ε turbulence model proves to be more general and reduces the discrepancy between the model predictions and the experimental data for heat transfer in rotating systems.


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.


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.


Author(s):  
Brian G. Wiedner ◽  
Cengiz Camci

The complex interaction between three-dimensional passage flow structure and endwall convective heat transfer in a square cross section, 90° turbulent duct flow has been experimentally investigated. Fine details of the momentum and heat transport process in a laboratory model that simulated a high Reynolds number three-dimensional passage flow are presented. The specific flow and heat transfer mechanisms are frequently encountered in the hot mainstream of axial flow turbines and internal coolant passages. Similar physical phenomena may also be observed in many other fluid machinery systems. The mean radius to duct width ratio was 2.3 and the Reynolds number based on inlet center line velocity, duct width, and ambient conditions was approximately 360,000. The complete Reynolds stress tensor was measured using a triple sensor hot wire. The turbulent normal and shear stresses, turbulent kinetic energy, and production of turbulent kinetic energy are presented. A steady state heat flux measurement technique and liquid crystal thermography were used to determine the character of the endwall heat transfer in the form of a high resolution heat transfer map. The flow field was dominated by strong counter rotating secondary flows characteristic of 90° turning ducts. The flow structure also included areas of strong streamwise accelerations and decelerations, high vorticity, local regions of significant total pressure loss, and a complex turbulent flow field structure. The development of the turbulent features of the 90° turning duct flow field and the influence of the turbulent flow field on the endwall convective heat transfer distribution are discussed. The multi-dimensional flow and high resolution heat transfer results are currently being incorporated in computational aerothermal models under development at Penn State University. The results are also available as a data base for future aerothermal model validation studies.


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 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.


2014 ◽  
Vol 18 (1) ◽  
pp. 1-16 ◽  
Author(s):  
Momir Sjeric ◽  
Darko Kozarac ◽  
Rudolf Tomic

The development of a two zone k-? turbulence model for the cycle-simulation software is presented. The in-cylinder turbulent flow field of internal combustion engines plays the most important role in the combustion process. Turbulence has a strong influence on the combustion process because the convective deformation of the flame front as well as the additional transfer of the momentum, heat and mass can occur. The development and use of numerical simulation models are prompted by the high experimental costs, lack of measurement equipment and increase in computer power. In the cycle-simulation codes, multi zone models are often used for rapid and robust evaluation of key engine parameters. The extension of the single zone turbulence model to the two zone model is presented and described. Turbulence analysis was focused only on the high pressure cycle according to the assumption of the homogeneous and isotropic turbulent flow field. Specific modifications of differential equation derivatives were made in both cases (single and two zone). Validation was performed on two engine geometries for different engine speeds and loads. Results of the cyclesimulation model for the turbulent kinetic energy and the combustion progress variable are compared with the results of 3D-CFD simulations. Very good agreement between the turbulent kinetic energy during the high pressure cycle and the combustion progress variable was obtained. The two zone k-? turbulence model showed a further progress in terms of prediction of the combustion process by using only the turbulent quantities of the unburned zone.


Author(s):  
Chung-Chu Chen ◽  
Tong-Miin Liou

Laser-Doppler velocimetry (LDV) measurements are presented of turbulent flow in a two-pass square-sectioned duct simulating the coolant passages employed in gas turbine blades under rotating and non-rotating conditions. For all cases studied, the Reynolds number characterized by duct hydraulic diameter (Dh) and bulk mean velocity (Ub) was fixed at 1 × 104. The rotating case had a range of rotation number (Ro = ΩDh/Ub) from 0 to 0.2. It is found that both the skewness of streamwise mean velocity and magnitude of secondary-flow velocity increase linearly, and the magnitude of turbulence intensity level increases non-linearly with increasing Ro. As Ro is increased, the curvature induced symmetric Dean vortices in the turn for Ro = 0 is gradually dominated by a single vortex most of which impinges directly on the outer part of leading wall. The high turbulent kinetic energy is closely related to the dominant vortex prevailing inside the 180-deg sharp turn. For the first time, the measured flow characteristics account for the reported spanwise heat transfer distributions in the rotating channels, especially the high heat transfer enhancement on the leading wall in the turn. For both rotating and non-rotating cases, the direction and strength of the secondary flow with respect to the wall are the most important fluid dynamic factors affecting local heat transfer distributions inside a 180-deg sharp turn. The role of the turbulent kinetic energy in affecting the overall enhancement of heat transfer is well addressed.


Author(s):  
Oguz Uzol ◽  
Cengiz Camci

A new concept for enhanced turbulent transport of heat in internal coolant passages of gas turbine blades is introduced. The new heat transfer augmentation component called “oscillator fin” is based on an unsteady flow system using the interaction of multiple unsteady jets and wakes generated downstream of a fluidic oscillator. Incompressible, unsteady and two dimensional solutions of Reynolds Averaged Navier-Stokes equations are obtained both for an oscillator fin and for an equivalent cylindrical pin fin and the results are compared. Preliminary results show that a significant increase in the turbulent kinetic energy level occur in the wake region of the oscillator fin with respect to the cylinder with similar level of aerodynamic penalty. The new concept does not require additional components or power to sustain its oscillations and its manufacturing is as easy as a conventional pin fin. The present study makes use of an unsteady numerical simulation of mass, momentum, turbulent kinetic energy and dissipation rate conservation equations for flow visualization downstream of the new oscillator fin and an equivalent cylinder. Relative enhancements of turbulent kinetic energy and comparisons of the total pressure field from transient simulations qualitatively suggest that the oscillator fin has excellent potential in enhancing local heat transfer in internal cooling passages without significant aerodynamic penalty.


Sign in / Sign up

Export Citation Format

Share Document