scholarly journals Numerical study of turbulent cavitating flows in thermal regime

2017 ◽  
Vol 27 (7) ◽  
pp. 1487-1503 ◽  
Author(s):  
Eric Goncalves ◽  
Dia Zeidan

Purpose The aim of this work is to quantify the relative importance of the turbulence modelling for cavitating flows in thermal regime. A comparison of various transport-equation turbulence models and a study of the influence of the turbulent Prandtl number appearing in the formulation of the turbulent heat flux are proposed. Numerical simulations are performed on a cavitating Venturi flow for which the running fluid is freon R-114 and results are compared with experimental data. Design/methodology/approach A compressible, two-phase, one-fluid Navier–Stokes solver has been developed to investigate the behaviour of cavitation models including thermodynamic effects. The code is composed by three conservation laws for mixture variables (mass, momentum and total energy) and a supplementary transport equation for the volume fraction of gas. The mass transfer between phases is closed assuming its proportionality to the mixture velocity divergence. Findings The influence of turbulence model as regard to the cooling effect due to the vaporization is weak. Only the k – ε Jones–Launder model under-estimates the temperature drop. The amplitude of the wall temperature drop near the Venturi throat increases with the augmentation of the turbulent Prandtl number. Originality/value The interaction between Reynolds-averaged Navier–Stokes turbulence closure and non-isothermal phase transition is rarely studied. It is the first time such a study on the turbulent Prandtl number effect is reported in cavitating flows.

2015 ◽  
Vol 776 ◽  
pp. 512-530 ◽  
Author(s):  
S. Leonardi ◽  
P. Orlandi ◽  
L. Djenidi ◽  
R. A. Antonia

Direct numerical simulations (DNS) are carried out to study the passive heat transport in a turbulent channel flow with either square bars or circular rods on one wall. Several values of the pitch (${\it\lambda}$) to height ($k$) ratio and two Reynolds numbers are considered. The roughness increases the heat transfer by inducing ejections at the leading edge of the roughness elements. The amounts of heat transfer and mixing depend on the separation between the roughness elements, an increase in heat transfer accompanying an increase in drag. The ratio of non-dimensional heat flux to the non-dimensional wall shear stress is higher for circular rods than square bars irrespectively of the pitch to height ratio. The turbulent heat flux varies within the cavities and is larger near the roughness elements. Both momentum and thermal eddy diffusivities increase relative to the smooth wall. For square cavities (${\it\lambda}/k=2$) the turbulent Prandtl number is smaller than for a smooth channel near the wall. As ${\it\lambda}/k$ increases, the turbulent Prandtl number increases up to a maximum of 2.5 at the crests plane of the square bars (${\it\lambda}/k=7.5$). With increasing distance from the wall, the differences with respect to the smooth wall vanish and at three roughness heights above the crests plane, the turbulent Prandtl number is essentially the same for smooth and rough walls.


Author(s):  
Hatem Kanfoudi ◽  
Ridha Zgolli

Cavitating flows inside a diesel injection nozzle hole were simulated using a mixture model. A two-dimensional (2D) numerical model is proposed in this paper to simulate steady cavitating flows. The Reynolds-averaged Navier–Stokes equations are solved for the liquid and vapor mixture, which is considered as a single fluid with variable density and expressed as a function of the vapor volume fraction. The closure of this variable is provided by the transport equation with a source term Transport-equation based methods (TEM). The processes of evaporation and condensation are governed by changes in pressure within the flow. The source term is implanted in the CFD code ANSYS CFX. The influence of numerical and physical parameters is presented in detail. The numerical simulations are in good agreement with the experimental data for steady flow.


2001 ◽  
Vol 124 (1) ◽  
pp. 73-93 ◽  
Author(s):  
R. P. Roy ◽  
S. Kang ◽  
J. A. Zarate ◽  
A. Laporta

Experiments and simulations were carried out in this investigation of turbulent subcooled boiling flow of Refrigerant-113 through a vertical annular channel whose inner wall only was heated. The measurements used, simultaneously, a two-component laser Doppler velocimeter for the liquid velocity field and a fast-response cold-wire for the temperature field, and a dual-sensor fiberoptic probe for the vapor fraction and vapor axial velocity. In the numerical simulation, the two-fluid model equations were solved by the solver ASTRID developed at Electricite´ de France. Wall laws for the liquid phase time-average axial velocity and temperature were developed from the experimental data, and the turbulent Prandtl number in the liquid was determined from the wall laws. The wall laws and turbulent Prandtl number were used in the simulations. The wall heat transfer model utilized the measured turbulent heat flux distribution in the liquid. Results from the simulations were compared with the measurements. Good agreement was found for some of the quantities while the agreement was only fair for others.


2019 ◽  
Vol 33 (20) ◽  
pp. 1950228
Author(s):  
Chunlai Tian ◽  
Tairan Chen ◽  
Tian Zou

Unsteady cavitating flow with high Reynolds number and significant instability commonly exists in fluid machinery and engineering system. The high-resolution approaches, such as direct numerical simulation and large eddy simulation, are not practical for engineering issues due to the significant cost in the computational resource. The objective of this paper is to provide the approach with Detached-Eddy Simulation (DES) model based on the Reynolds-averaged Navier–Stokes (RANS) equations for predicting unsteady cavitating flows. The credibility of the approach is validated by a set of numerical examples of its application: the unsteady cavitating flows around the two-dimensional (2D) Clark-Y hydrofoil and the three-dimensional (3D) blunt body. It is found that the calculated cavity shapes, cavity lengths and unsteady characteristics by DES model agree well with the experimental measurements and observations. Further analysis indicates that the turbulent eddy viscosity around the cavity and wake region is well predicted by the DES model, which results in the development of large-scale vortexes, and further cavitation instability. The DES model, which exhibits a significantly unsteady 3D behavior, is a more comprehensive turbulence model for unsteady cavitating flows.


1999 ◽  
Vol 121 (2) ◽  
pp. 322-325 ◽  
Author(s):  
C. U. Buice ◽  
J. K. Eaton

The recovering boundary layer downstream of a separation bubble is known to have a highly perturbed turbulence structure which creates difficulty for turbulence models. The present experiment addressed the effect of this perturbed structure on turbulent heat transport. The turbulent diffusion of heat downstream of a heated wire was measured in a perturbed channel flow and compared to that in a simple, fully developed channel flow. The turbulent diffusivity of heat was found to be more than 20 times larger in the perturbed flow. The turbulent Prandtl number increased to 1.7, showing that the turbulent eddy viscosity was affected even more strongly than the eddy thermal diffusivity. This result corroborates previous work showing that boundary layer disturbances generally have a stronger effect on the eddy viscosity, rendering prescribed turbulent Prandtl number models ineffective in perturbed flows.


Author(s):  
Dawid Taler

Purpose The purpose of this paper is to develop new semi-empirical heat transfer correlations for turbulent flow of liquid metals in the tubes, and then to compare these correlations with the experimental data. The Prandtl and Reynolds numbers can vary in the ranges: 0.0001 ≤ Pr ≤ 0.1 and 3000 ≤ Re ≤ 106. Design/methodology/approach The energy conservation equation averaged by Reynolds was integrated using the universal velocity profile determined experimentally by Reichardt for the turbulent tube flow and four different models for the turbulent Prandtl number. Turbulent heat transfer in the circular tube was analyzed for a constant heat flux at the inner surface. Some constants in different models for the turbulent Prandtl number were adjusted to obtain good agreement between calculated and experimentally obtained Nusselt numbers. Subsequently, new correlations for the Nusselt number as a function of a Peclet number was proposed for different models of the turbulent Prandtl number. Findings The inclusion of turbulent Prandtl number greater than one and the experimentally determined velocity profile of the fluid in the tube while solving the energy conservation equation improved the compatibility of calculated Nusselt numbers, with Nusselt numbers determined experimentally. The correlations proposed in the paper have a sound theoretical basis and give Nusselt number values that are in good agreement with the experimental data. Research limitations/implications Heat transfer correlations proposed in this paper were derived assuming a constant heat flux at the inner surface of the tube. However, they can also be used for a constant wall temperature, as for the turbulent flow (Re > 3,000), the relative difference between the Nusselt number for uniform wall heat flux and uniform wall temperature is very low. Originality/value Unified, systematic approach to derive correlations for the Nusselt number for liquid metals was proposed in the paper. The Nusselt number was obtained from the solution of the energy conservation equation using the universal velocity profile and eddy diffusivity determined experimentally, and various models for the turbulent Prandtl number. Four different relationships for the Nusselt number proposed in the paper were compared with the experimental data.


1981 ◽  
Vol 103 (1) ◽  
pp. 153-158 ◽  
Author(s):  
H. W. Coleman ◽  
R. J. Moffat ◽  
W. M. Kays

Heat transfer behavior of a fully rough turbulent boundary layer subjected to favorable pressure gradients was investigated experimentally using a porous test surface composed of densely packed spheres of uniform size. Stanton numbers and profiles of mean temperature, turbulent Prandtl number, and turbulent heat flux are reported. Three equilibrium acceleration cases (one with blowing) and one non-equilibrium acceleration case were studied. For each acceleration case of this study, Stanton number increased over zero pressure gradient values at the same position or enthalpy thickness. Turbulent Prandtl number was found to be approximately constant at 0.7–0.8 across the layer, and profiles of the non-dimensional turbulent heat flux showed close agreement with those previously reported for both smooth and rough wall zero pressure gradient layers.


Volume 4 ◽  
2004 ◽  
Author(s):  
X. Cheng ◽  
A. Batta ◽  
H. Y. Chen ◽  
N. I. Tak

The present paper gives a brief literature review on turbulent heat transfer in heavy liquid metals (HLM), especially liquid lead-bismuth eutectic (LBE). Some models available in the open literature on heat transfer and turbulent Prandtl number are assessed. In addition, CFD analysis is carried out for circular tube geometries. The effect of turbulence models, mesh structure and turbulent Prandtl number on the numerical results is studied. Application of ε-type turbulence models with scalable wall function shows less dependence of the numerical results on mesh structure than the ω-type turbulence models with automatic wall treatment. The turbulent Prandtl number affects strongly the heat transfer performance. Comparison between the CFD results, heat transfer correlations and heat transfer test data reveals a decrease in turbulent Prandtl number by increasing Reynolds number. Based on the results achieved, recommendations are made on correlations of heat transfer and turbulent Prandtl number for LBE flows.


Sign in / Sign up

Export Citation Format

Share Document