An explicit algebraic Reynolds stress and heat flux model for incompressible turbulence: Part II Buoyant flow

This paper presents a derivation of an explicit algebraic model for two-dimensional (2-D) buoyant flows. It is an extension of the work reported in Part I (So et al. [27]). The tensor representation method of Jongen and Gatski [14] is extended to derive an explicit algebraic Reynolds stress model (EASM) for 2-D buoyant flow invoking the Boussinesq approximation. The projection methodology is further extended to treat the heat flux transport equation in the derivation of an explicit algebraic heat flux model (EAHFM) for buoyant flow. Again, the weak equilibrium assumption is invoked for the scaled Reynolds stress and scaled heat flux equation. An explicit algebraic model for buoyant flows is then formed with the EASM and EAHFM. From the derived EAHFM, an expression for the thermal diffusivity tensor in buoyant shear flows is deduced. Furthermore, a turbulent Prandtl number (Pr_T) for each of the three heat flux directions is determined. These directional Pr_T are found to be a function of the gradient Richardson number. Alternatively, a scalar Pr_T can be derived; its value is compared with the directional Pr_T. The EASM and EAHFM are used to calculate 2-D homogeneous buoyant shear flows and the results are compared with direct numerical simulation data and other model predictions, where good agreement is obtained. Dedicated to the memory of the late Professor Charles G. Speziale of Boston University     

An explicit algebraic Reynolds stress model in turbulence

A new algebraic Reynold stress model is constructed with recourse to the realizability constraints. Model coefficients are made a function of strain and vorticity invariants through calibration by reference to homogeneous shear flow data. The anisotropic production in near‐wall regions is accounted for substantially by modifying the model constants C_{ε(1, 2)} and adding a secondary source term in the ε equation. Hence, it reduces the kinetic energy and length scale magnitudes to improve prediction of adverse pressure gradient flows, involving flow separation and reattachment. To facilitate the evaluation of the turbulence model, some extensively used benchmark cases in the turbulence modelling are convoked. The comparisons demonstrate that the new model maintains qualitatively good agreement with the direct numerical simulation (DNS) and experimental data. Copyright © 2006 John Wiley & Sons, Ltd.     

Modifications for an explicit algebraic stress model

An extension of the explicit algebraic stress model, developed by Gatski and Speziale [Gatski TB, Speziale CG. On the explicit algebraic stress models for complex turbulent flows. Journal of Fluid Mechanics 1993; 254: 59–78] is proposed. The extension implicates some essential characteristics of second‐order closure models. The strain‐dependent coefficients are modified, resulting in an alleviation of the numerical instabilities involved in the model. A new near‐wall damping function f_μ in the eddy viscosity relation is introduced. To enhance dissipation in near‐wall regions, the model constant C_ϵ1 is modified and an extra positive source term is included in the dissipation equation. In addition, a realizable time scale is incorporated to remove the wall singularity. Computed results show that the modified Gatski–Speziale (MGS) model predictions are in better agreement with the direct numerical simulation (DNS) and experimental data than those of the original Gatski–Speziale (OGS) model. Copyright © 2001 John Wiley & Sons, Ltd.     

The correction of longitudinal heat flux in turbulence measurement

In this paper the influence of temperature on velocity signal in hot-wire measurement of turbulence is analysed. It is pointed out that when the temperature influence is small, the temperature influence on measured intensity of velocity fluctuations is second order small and negligible. However, the temperature influence on measuring longitudinal heat flux is of first order quantity, and must be corrected, or large error will occur. The method to correct the temperature influence on measuring ρ _θu and the procedure to decide experimentally temperature influence coefficient have been given.     

Fully explicit and self-consistent algebraic Reynolds stress model

A fully explicit, self-consistent algebraic expression (for Reynolds stress) which is the exact solution to the Reynolds stress transport equation in the weak-equilibrium limit for two-dimensional mean flows for all linear and some quasi-linear pressure-strain models, is derived. Current explicit algebraic Reynolds stress models derived by employing the weak-equilibrium assumption treat the production-to-dissipation (P/) ratio as a constant, resulting in an effective viscosity that can be singular away from the equilibrium limit. In this paper the set of simultaneous algebraic Reynolds stress equations in the weak-equilibrium limit are solved in the full nonlinear form and the eddy viscosity is found to be nonsingular. Preliminary tests indicate that the model performs adequately, even for three-dimensional mean-flow cases. Due to the explicit and nonsingular nature of the effective viscosity, this model should mitigate many of the difficulties encountered in computing complex turbulent flows with the algebraic Reynolds stress models.This research was supported by the National Aeronautics and Space Administration under NASA Contract No. NAS1-19480.     

旋转湍流中雷诺应力模型压力应变关联项的修正研究

利用张量的不变量理论,推导得出传统雷诺应力模型中压力应变关联项模型应用于旋转湍流模拟中的一些基本问题,即在纯旋转条件下,传统模型所描述的初始各向异性的湍流中雷诺应力张量演化规律是一个无衰减振荡过程,而快速畸变理论推导结果显示,其演化应是一个阻尼振荡衰减的过程。以衰减雷诺应力为目的,构造出包含旋转率张量高阶量的关联项。然后,结合变形率张量的高阶项,将修正模型扩展至椭圆形流线类型流动。最后,将修正模型应用于轴向旋转圆管内湍流流场的模拟,并将结果与实测结果进行了对比。     

Spatial artificial neural network model for subgrid-scale stress and heat flux of compressible turbulence

The subgrid-scale(SGS) stress and SGS heat flux are modeled by using an artificial neural network(ANN) for large eddy simulation(LES) of compressible turbulence. The input features of ANN model are based on the first-order and second-order derivatives of filtered velocity and temperature at different spatial locations. The proposed spatial artificial neural network(SANN)model gives much larger correlation coefficients and much smaller relative errors than the gradient model in an a priori analysis. In an a posteriori analysis, the SANN model performs better than the dynamic mixed model(DMM) in the prediction of spectra and statistical properties of velocity and temperature, and the instantaneous flow structures.     

Study on anisotropic buoyant turbulence model

Ｂｕｏｙａｎｔｆｌｏｗｉｓｏｎｅｏｆｔｈｅｆｕｎｄａｍｅｎｔａｌｆｌｏｗｓ．Ｔｈｅｄｉｆｆｅｒｅｎｃｅｏｆｄｅｎｓｉｔｙｂｅｔｗｅｅｎｄｉｓｃｈａｒｇｅｄｆｌｕｉｄａｎｄｔｈｅａｍｂｉｅｎｔｆｌｕｉｄｃａｎｃａｕｓｅｂｕｏｙａｎｔｔｕｒｂｕｌｅｎｔｆ...     

Macro-to-microchannel transition in two-phase flow: Part 2 - Flow boiling heat transfer and critical heat flux

This part of the paper presents the current experimental flow boiling heat transfer and CHF data acquired for R134a, R236fa and R245fa in single, horizontal channels of 1.03, 2.20 and 3.04 mm diameters over a range of experimental conditions. The aim of this study is to investigate the effects of channel confinement, heat flux, flow pattern, saturation temperature, subcooling and working fluid properties on the two-phase heat transfer and CHF. Experimentally, it was observed that the flow boiling heat transfer coefficients are a significant function of the type of two-phase flow pattern. Furthermore, the monotonically increasing heat transfer coefficients at higher vapor qualities, corresponding to annular flow, signifies convective boiling as the dominant heat transfer mechanism in these small scale channels. The decreasing heat transfer trend at low vapor qualities in the slug flow (coalescing bubble dominated regime) was indicative of thin film evaporation with intermittent dry patch formation and rewetting at these conditions. The coalescing bubble flow heat transfer data were well predicted by the three-zone model when setting the dryout thickness to the measured surface roughness, indicating for the first time a roughness effect on the flow boiling heat transfer coefficient in this regime. The CHF data acquired during the experimental campaign indicated the influence of saturation temperature, mass velocity, channel confinement and fluid properties on CHF but no influence of inlet subcooling for the conditions tested. When globally comparing the CHF values for R134a in the 0.51-3.04 mm diameter channels, a peak in CHF peak was observed lying in between the 0.79 (Co ≈ 0.99) and 1.03 (Co ≈ 0.78) mm channels. A new CHF correlation has been proposed involving the confinement number, Co that is able to predict CHF for R134a, R236fa and R245fa in single-circular channels, rectangular multichannels and split flow rectangular multichannels. In summary, the present flow boiling and CHF trends point to a macro-to-microscale transition as indicated by the results presented in Ong and Thome (2011) [1].     

Finite element implementation of two‐equation and algebraic stress turbulence models for steady incompressible flows

The main purpose of this paper is to describe a finite element formulation for solving the equations for k and ε of the classical k–ε turbulence model, or any other two‐equation model. The finite element discretization is based on the SUPG method together with a discontinuity capturing technique to deal with sharp internal and boundary layers. The iterative strategy consists of several nested loops, the outermost being the linearization of the Navier–Stokes equations. The basic k–ε model is used for the implementation of an algebraic stress model that is able to account for the effects of rotation. Some numerical examples are presented in order to show the performance of the proposed scheme for simulating directly steady flows, without the need of reaching the steady state through a transient evolution. Copyright © 1999 John Wiley & Sons, Ltd.     

Dynamic coefficient evaluation for an algebraic subgrid stress model

Static model coefficients for an algebraic subgrid stress (SGS) model are determined using a dynamic approach, based on results from simulations of isotropic decaying turbulence. The study was motivated by the discrepancies in energy transfer predictions using the previously documented coefficients (Bhushan and Warsi, Int. J. Numer. Meth. Fluids 2005; 49 : 489–519). The discrepancies are identified to be due to inconsistent filter functions used in the analytic estimates and the simulations. The study emphasizes that SGS model development should use filter functions compatible with those inherent in CFD application solvers. The dynamic approach predicts consistent model and transfer coefficients for different grid resolutions and is judged to be a reliable basis for model coefficient adjustments. The predicted Leonard's stress coefficient and associated energy transfer coefficients agree very well with the analytic estimates using a Gaussian/cutoff combination filter. This suggests that the modeling of Leonard's stress term using a truncated Taylor series expansion is robust and may not benefit significantly from dynamic modeling. Validation simulations were performed for turbulent channel flow at Re_τ = 180 and 590. The dynamic approach was found to be reliable only for the lower log‐layer of the Re_τ = 590 case, where the scale invariance condition was satisfied. Nonetheless, in this narrow range, the model and transfer coefficients compare well with the isotropic case. The static coefficient algebraic model with new adjusted coefficients shows improved predictions compared with the previous coefficients, for both isotropic decaying turbulence and channel flow cases. Copyright © 2013 John Wiley & Sons, Ltd.     

Cattaneo-Christov heat flux model for third-grade fluid flow towards exponentially stretching sheet

The Cattaneo-Christov heat flux in the two-dimensional (2D) flow of a third-grade fluid towards an exponentially stretching sheet is investigated. The energy equation is considered through thermal relaxation. Similarity transformations are accounted to obtain the ordinary differential systems. The converted non-dimensional equations are solved for the series solutions. The convergence analysis of the computed solutions is reported. The graphical results of the velocity and temperature profiles are plotted and elaborated in detail. The results show that the thermal relaxation enhances the temperature gradient while reduces the temperature profile.     

Numerical analysis of developing turbulent flow in a 180° bend tube by an algebraic Reynolds stress model

A numerical analysis has been performed for three‐dimensional developing turbulent flow in a 180° bend tube with straight inlet and outlet section used by an algebraic Reynolds stress model. To our knowledge, numerical investigations, which show the detailed comparison between calculated results and experimental data including distributions of Reynolds stresses, are few and far between. From this point of view, an algebraic Reynolds stress model in conjunction with boundary‐fitted co‐ordinate system is applied to a 180° bend tube in order to predict the anisotropic turbulent structure precisely. Calculated results are compared with the experimental data including distributions of Reynolds stresses. As a result of this analysis, it has been found that the calculated results show a comparatively good agreement with the experimental data of the time‐averaged velocity and the secondary vectors in both the bent tube and straight outlet sections. For example, the location of the maximum streamwise velocity, which appears near the top or bottom wall in the bent tube, is predicted correctly by the present method. As for the comparison of Reynolds stresses, the present method has been found to simulate many characteristic features of streamwise normal stress and shear stresses in the bent tube qualitatively and has a tendency to under‐predict its value quantitatively. Judging from the comparison between the calculated and the experimental results, the algebraic Reynolds stress model is applicable to the developing turbulent flow in a bent tube that is known as a flow with a strong convective effect. Copyright © 2005 John Wiley & Sons, Ltd.     

Numerical analysis of turbulent structure in compound meandering open channel by algebraic Reynolds stress model

The turbulent flow in a compound meandering channel with a rectangular cross section is one of the most complicated turbulent flows, because the flow behaviour is influenced by several kinds of forces, including centrifugal forces, pressure‐driven forces and shear stresses generated by momentum transfer between the main channel and the flood plain. Numerical analysis has been performed for the fully developed turbulent flow in a compound meandering open‐channel flow using an algebraic Reynolds stress model. The boundary‐fitted coordinate system is introduced as a method for coordinate transformation in order to set the boundary conditions along the complicated shape of the meandering open channel. The turbulence model consists of transport equations for turbulent energy and dissipation, in conjunction with an algebraic stress model based on the Reynolds stress transport equations. With reference to the pressure–strain term, we have made use of a modified pressure–strain term. The boundary condition of the fluctuating vertical velocity is set to zero not only for the free surface, but also for computational grid points next to the free surface, because experimental results have shown that the fluctuating vertical velocity approaches zero near the free surface. In order to examine the validity of the present numerical method and the turbulent model, the calculated results are compared with experimental data measured by laser Doppler anemometer. In addition, the compound meandering open channel is clarified somewhat based on the calculated results. As a result of the analysis, the present algebraic Reynolds stress model is shown to be able to reasonably predict the turbulent flow in a compound meandering open channel. Copyright © 2005 John Wiley & Sons, Ltd.     

Explicit algebraic Reynolds stress model for shock-dominated flows

Shock waves drastically alter the nature of Reynolds stresses in a turbulent flow, and conventional turbulence models cannot reproduce this effect. In the present study, we employ explicit algebraic Reynolds stress model (EARSM) to predict the Reynolds stress anisotropy generated by a shockwave. The model by Wallin and Johansson (2000) is used as the baseline model. It is found to over-predict the post-shock Reynolds stresses in canonical shock turbulence interaction. The budget of the transport equation of Reynolds stresses computed using linear interaction analysis shows that the unsteady shock distortion mechanism and the pressure–velocity correlations are important. We propose improvement to the baseline model using linear interaction analysis results and redistribute the turbulent kinetic energy between the principle Reynolds stresses. The new model matches DNS data for the amplification of Reynolds stresses across the shock and their post-shock evolution, for a range of Mach numbers. It is applied to oblique shock/boundary-layer interaction at Mach 5. Significant improvements are observed in predicting surface pressure and skin friction coefficient, with respect to experimental measurements.     

Numerical analysis of turbulent flow separation in a rectangular duct with a sharp 180‐degree turn by algebraic Reynolds stress model

Turbulent flow in a rectangular duct with a sharp 180‐degree turn is difficult to predict numerically because the flow behavior is influenced by several types of forces, including centrifugal force, pressure‐driven force, and shear stress generated by anisotropic turbulence. In particular, this type of flow is characterized by a large‐scale separated flow, and it is difficult to predict the reattachment point of a separated flow. Numerical analysis has been performed for a turbulent flow in a rectangular duct with a sharp 180‐degree turn using the algebraic Reynolds stress model. A boundary‐fitted coordinate system is introduced as a method for coordinate transformation to set the boundary conditions next to complicated shapes. The calculated results are compared with the experimental data, as measured by a laser‐Doppler anemometer, in order to examine the validity of the proposed numerical method and turbulent model. In addition, the possibility of improving the wall function method in the separated flow region is examined by replacing the log‐law velocity profile for a smooth wall with that for a rough wall. The analysis results indicated that the proposed algebraic Reynolds stress model can be used to reasonably predict the turbulent flow in a rectangular duct with a sharp 180‐degree turn. In particular, the calculated reattachment point of a separated flow, which is difficult to predict in a turbulent flow, agrees well with the experimental results. In addition, the calculation results suggest that the wall function method using the log‐law velocity profile for a rough wall over a separated flow region has some potential for improving the prediction accuracy. Copyright © 2007 John Wiley & Sons, Ltd.     

Large eddy simulation of turbulent channel flow using an algebraic model

In this paper an algebraic model from the constitutive equations of the subgrid stresses has been developed. This model has an additional term in comparison with the mixed model, which represents the backscatter of energy explicitly. The proposed model thus provides independent modelling of the different energy transfer mechanisms, thereby capturing the effect of subgrid scales more accurately. The model is also found to depict the flow anisotropy better than the linear and mixed models. The energy transfer capability of the model is analysed for the isotropic decay and the forced isotropic turbulence. The turbulent plane channel flow simulation is performed over three Reynolds numbers, Re_τ=180, 395 and 590, and the results are compared with that of the dynamic model, Smagorinsky model, and the DNS data. Both the algebraic and dynamic models are in good agreement with the DNS data for the mean flow quantities. However, the algebraic model is found to be more accurate for the turbulence intensities and the higher‐order statistics. The capability of the algebraic model to represent backscatter is also demonstrated. Copyright © 2005 John Wiley & Sons, Ltd.     

Explicit algebraic stress modelling of homogeneous and inhomogeneous flows

Capability of the explicit algebraic stress models to predict homogeneous and inhomogeneous shear flows are examined. The importance of the explicit solution of the production to dissipation ratio is first highlighted by examining the algebraic stress models performance at purely irrotational strain conditions. Turbulent recirculating flows within sudden expanding pipes are further simulated with explicit algebraic stress model and anisotropic eddy viscosity model. Both models predict better stress–strain interactions, showing reasonable shear layer developments. The anisotropic stress field are also accurately predicted by the models, though the anisotropic eddy viscosity model of Craft et al. returns marginally better results. Copyright © 2005 John Wiley & Sons, Ltd.