An elliptic blending differential flux model for natural,mixed and forced convection

Several modifications are introduced to the Elliptic Blending Differential Flux Model proposed by Shin et al. (2008) to account for the influence of wall blockage on the turbulent heat flux. These modifications are introduced in order to reproduce, in association with the most recent version of the EB-RSM, the full range of regimes, from forced to natural convection, without any case-specific modification. The interest of the new model is demonstrated using analytical arguments, a priori tests and computations in channel flows in the different convection regimes, as well as in a differentially heated cavity.     

宽速域下神经网络对雷诺应力各向异性张量的预测

基于Pope修正的有效黏度假设,张量基神经网络(tensor based neural network,TBNN)构建了从雷诺平均方程湍流模型(RANS)的平均应变率张量和平均旋转率张量到高精度数值解的雷诺应力各向异性张量的映射.将高精度数值解用于TBNN的训练,从而使TBNN根据RANS求解的湍动能、湍流耗散率和速度...     

Rotation invariant constitutive relation for Reynolds stress structure parameter

A new Reynolds stress constitutive formula is constructed using the firstorder statistics of turbulent fluctuations instead of the mean strain rate. It includes zero empirical coefficients. The formula is validated with the direct numerical simulation(DNS) data of turbulent channel flow at Reτ =180. The Reynolds stresses given by the proposed formula agree very well with the DNS results. The good agreement persists even after the multi-angle rotation of the coordinate system, indicating the rotation invariance of the formula. The autocorrelation of the fluctuating velocity rather than the mean strain rate is close to the essence of the Reynolds stress.     

Application of Reynolds stress transport turbulence closure models to flows affected by Lorentz and buoyancy forces

The effect of magnetic field strength and orientation on two types of electromagnetically influenced turbulent flows was studied numerically under the Reynolds averaged Navier–Stokes (RANS) framework. Previous work (Wilson et al., 2014) used an electromagnetically extended linear eddy-viscosity model, whilst the current paper focuses on the performance of a more advanced Reynolds stress transport type model both with and without electromagnetic modifications proposed by Kenjereš et al. (2004). First, a fully-developed 2D channel flow is considered with a magnetic field imposed in either the wall-normal or streamwise direction. Both forms of the RSM gave good agreement with the DNS data for the wall-normal magnetic field across the range of Hartmann numbers with the additional electromagnetic terms providing a small, but noticeable, difference. For the streamwise magnetic field, where electromagnetic influence is only through the turbulence, the electromagnetically extended RSM performed well at moderate Hartmann numbers but returned laminar flow at the highest Hartmann number considered, contrary to the DNS. The RSM results were, however, significantly better than the previous eddy-viscosity model predictions. The second case is that of unsteady 3D Rayleigh–Bénard convection with a magnetic field imposed in either a horizontal or vertical direction. Results revealed that a significant reorganization of the flow structures is predicted to occur. For a vertically oriented magnetic field, the plume structures increase in number and become thinner and elongated along the magnetic field lines, leading to an increase in thermal mixing within the core in agreement with Hanjalić and Kenjereš (2000). With a horizontal magnetic field, the structures become two-dimensional and a striking realignment of the roll cells’ axes with the magnetic field lines occurs. The results demonstrate the capability of the Reynolds stress transport approach in modelling MHD flows that are relevant to industry and offer potential for those wishing to control levels of turbulence, heat transfer or concentration without recourse to mechanical means.     

Low‐Reynolds number $k-\tilde {\varepsilon}$ modelling of flows with transpiration

An improved low‐Reynolds number model was adopted to predict the dynamic and thermal fields in flows with transpiration. The performance of the adopted model was first contrasted with the direct numerical simulation (DNS) data of channel flow with uniform wall injection and suction. The validity of the present model applied to flows with a high level of transpiration was further examined. To explore the model's performance in complex environments, the model was applied to simulate a transpired developing channel flow. By contrasting the predictions with DNS data and measurements, the results indicated that the present model reproduced correctly the deceleration and acceleration of the flow caused by the injection and suction from the permeable part of the wall. The turbulence structure of transpired flows was also well captured and the superior performance of the adopted model was reflected by the predicted correct level of ε with the maximum being located at both the injection and the suction walls. The predicted thermal field by the present model also compared favourably with the DNS data and measurements. Copyright © 2000 John Wiley & Sons, Ltd.     

Finite element simulation of turbulent Couette–Poiseuille flows using a low Reynolds number k–ϵ model

Developing Couette–Poiseuille flows at Re=5000 are studied using a low Reynolds number k–ϵ two‐equation model and a finite element formulation. Mesh‐independent solutions are obtained using a standard Galerkin formulation and a Galerkin/least‐squares stabilized method. The predictions for the velocity and turbulent kinetic energy are compared with available experimental results and to the DNS data. Second moment closure's solutions are also compared with those of the k–ϵ model. The deficiency of eddy viscosity models to predict dissymmetric low Reynolds number channel flows has been demonstrated. Copyright © 1999 John Wiley & Sons, Ltd.     

Evaluation of Extended Weak-Equilibrium Conditions for Fully Developed Rotating Channel Flow

The weak-equilibrium condition, which is the basis in the development of algebraic Reynolds stress models is assessed in the prediction of fully developed turbulent channel flow under the influence of system rotation. The budget of the various terms of the exact transport equation for the anisotropy tensor is evaluated by using a DNS database. Two diffusive transport constraints are evaluated by using the DNS data. The results show that neither of them can hold for the near-wall region. An asymptotic analysis of the near-wall behavior is performed and an alternative form of the diffusive transport constraint is proposed. The analysis shows that the proposed alternative diffusive transport constraint has the potential to improve the predictive ability of the resultant ARSM.     

A priori evaluations and least-squares optimizations of turbulence models for fully developed rotating turbulent channel flow

The present study involves a priori tests of pressure-strain and dissipation rate tensor models using data from direct numerical simulations (DNS) of fully developed turbulent channel flow with and without spanwise system rotation. Three different pressure-strain rate models are tested ranging from a simple quasi-linear model to a realizable fourth order model. The evaluations demonstrate the difficulties of developing RANS-models that accurately describe the flow for a wide range of rotation numbers. Furthermore, least-squares based tensor representations of the exact pressure-strain and dissipation rate tensors are derived pointwise in space. The relation obtained for the rapid pressure-strain rate is exact for general 2D mean flows. Hence, the corresponding distribution of the optimized coefficients show the ideal behaviour. The corresponding representations for the slow pressure-strain and dissipation rate tensors are incomplete but still optimal in a least-squares sense. On basis of the least-squares analysis it is argued that the part of the representation that is tensorially linear in the Reynolds stress anisotropy is the most important for these parts.     

LES-based filter-matrix lattice Boltzmann model for simulating fully developed turbulent channel flow

In this paper, a three-dimensional filter-matrix lattice Boltzmann (FMLB) model based on large eddy simulation (LES) was verified for simulating wall-bounded turbulent flows. The Vreman subgrid-scale model was employed in the present FMLB–LES framework, which had been proved to be capable of predicting turbulent near-wall region accurately. The fully developed turbulent channel flows were performed at a friction Reynolds number Re_τ of 180. The turbulence statistics computed from the present FMLB–LES simulations, including mean stream velocity profile, Reynolds stress profile and root-mean-square velocity fluctuations greed well with the LES results of multiple-relaxation-time (MRT) LB model, and some discrepancies in comparison with those direct numerical simulation (DNS) data of Kim et al. was also observed due to the relatively low grid resolution. Moreover, to investigate the influence of grid resolution on the present LES simulation, a DNS simulation on a finer gird was also implemented by present FMLB–D3Q19 model. Comparisons of detailed computed various turbulence statistics with available benchmark data of DNS showed quite well agreement.     

Predicting Noninertial Effects with Linear and Nonlinear Eddy-Viscosity, and Algebraic Stress Models

Three types of turbulence models which account for rotational effects in noninertial frames of reference are evaluated for the case of incompressible, fully developed rotating turbulent channel flow. The different types of models are a Coroiolis-modified eddy-viscosity model, a realizable nonlinear eddy-viscosity model, and an algebraic stress model which accounts for dissipation rate anisotropies. A direct numerical simulation of a rotating channel flow is used for the validation of the turbulence models. This simulation differs from previous studies in that significantly higher rotation numbers are investigated. Flows at these higher rotation numbers are characterized by a relaminarization on the cyclonic or suction side of the channel, and a linear velocity profile on the anticyclonic or pressure side of the channel. The predictive performance of the three types of models are examined in detail, and formulation deficiencies are identified which cause poor predictive performance for some of the models. Criteria are identified which allow for accurate prediction of such flows by algebraic stress models and their corresponding Reynolds stress formulations.     

Direct Numerical Simulation of a Recirculating, Swirling Flow

A direct numerical simulation (DNS) of a recirculating, swirling flow is performed at a Reynolds number of 5000. Detailed one and two point statistics are presented in this paper. Flow visualization and frequency analysis are used to identify a precessing vortex core and to characterize its position, extent and influence on the flow field. The results are compared with laser Doppler velocimetry (LDV) measurements as well as large eddy simulation (LES) data reported in the literature. The present work constitutes a first step in setting up a DNS data base for complex flows.     

Numerical investigation of turbulent channel flow controlled by spatially oscillating spanwise Lorentz force

A formulation of the skin-friction drag related to the Reynolds shear stress in a turbulent channel flow is derived. A direct numerical simulation (DNS) of the turbulent control is performed by imposing the spatially oscillating spanwise Lorentz force. Under the action of the Lorentz force with several proper control parameters, only the periodically well-organized streamwise vortices are finally observed in the near-wall region. The Reynolds shear stress decreases dramatically, especially in the near-wall area, resulting in a drag reduction.     

Turbulent flow around a wall-mounted cube: A direct numerical simulation

An immersed-boundary method was employed to perform a direct numerical simulation (DNS) of flow around a wall-mounted cube in a fully developed turbulent channel for a Reynolds number Re = 5610, based on the bulk velocity and the channel height. Instantaneous results of the DNS of a plain channel flow were used as a fully developed inflow condition for the main channel. The results confirm the unsteadiness of the considered flow caused by the unstable interaction of a horseshoe vortex formed in front of the cube and on both its sides with an arch-type vortex behind the cube. The time-averaged data of the turbulence mean-square intensities, Reynolds shear stresses, kinetic energy and dissipation rate are presented. The negative turbulence production is predicted in the region in front of the cube where the main horseshoe vortex originates.     

Passive Scalar in a Turbulent Channel Flow with Wall Velocity Disturbances

Direct Numerical Simulations (DNS) of a passive scalar in a turbulent channel flow with a normal velocity disturbance on the lower wall are presented for high and low Reynolds numbers. The aim is to reproduce the complex physics of turbulent rough flows without dealing with the geometric complexity. In addition, isothermal walls that cannot be easily assigned in an experiment, are considered. The paper explains the increase of heat transfer through the changes of the velocity and thermal structures. As in real rough flows, the transpiration produces an isotropization of the turbulence near the wall.     

Examination of wall damping for the k-ε turbulence model using direct simulations of turbulent channel flow

Handler, Hendricks and Leighton have recently reported results for the direct numerical simulation (DNS) of a turbulent channel flow at moderate Reynolds number. These data are used to evaluate the terms in the exact and modelled transport equations for the turbulence kinetic energy k and the isotropic dissipation function ε. Both modelled transport equations show significant imbalances in the high-shear region near the channel walls. The model for the eddy viscosity is found to yield distributions for the production terms which do not agree well with the distributions calculated from the DNS data. The source of the imbalance is attributed to the wall-damping function required in eddy viscosity models for turbulent flows near walls. Several models for the damping function are examined, and it is found that the models do not vary across the channel as does the damping when evaluated from the DNS data. The Lam-Bremhorst model and the standard van Driest model are found to give reasonable agreement with the DNS data. Modification of the van Driest model to include an effective origin yields very good agreement between the modelled production and the production calculated from the DNS data, and the imbalance in the modelled transport equations is significantly reduced.     

Implicit LES of free and wall‐bounded turbulent flows based on the discontinuous Galerkin/symmetric interior penalty method

This paper presents the second validation step of a compressible discontinuous Galerkin solver with symmetric interior penalty (DGM/SIP) for the direct numerical simulation (DNS) and the large eddy simulation (LES) of complex flows. The method has already been successfully validated for DNS of an academic flow and has been applied to flows around complex geometries (e.g. airfoils and turbomachinery blades). During these studies, the advantages of the dissipation properties of the method have been highlighted, showing a natural tendency to dissipate only the under‐resolved scales (i.e the smallest scales present on the mesh), leaving the larger scales unaffected. This phenomenon is further enhanced as the polynomial order is increased. Indeed, the order increases the dissipation at the largest wave numbers, while its range of impact is reduced. These properties are spectrally compatible with a subgrid‐scale model, and hence DGM may be well suited to be used for an implicit LES (ILES) approach. A validation of this DGM/ILES approach is here investigated on canonical flows, allowing to study the impact of the discretisation on the turbulence for under‐resolved computations. The first test case is the LES of decaying homogeneous isotropic turbulence (HIT) at very high Reynolds number. This benchmark allows to assess the spectral behaviour of the method for implicit LES. The results are in agreement with theory and are even slightly more accurate than other numerical results from literature, obtained using a pseudo‐spectral (PS) method with a state‐of‐the‐art subgrid‐scale model. The second benchmark is the LES of the channel flow. Three Reynolds numbers are considered: Re_τ=395, 590 and 950. The results are compared with DNS of Moser et al. and Hoyas et al., also using PS methods. Both averaged velocity and fluctuations are globally in good agreement with the reference, showing the ability of the method to predict equilibrium wall‐bounded flow turbulence. To show that the method is able to perform accurate DNS, a DNS of HIT at Re_λ=64 and a DNS of the channel flow at Re_τ=180 are also performed. The effects of the grid refinement are investigated on the channel flow at Re_τ=395, highlighting the improvement of the results when refining the mesh in the spanwise direction. Finally, the modification of the ILES parameters, that is the Riemann solver and of the SIP coefficient, is studied on both cases, showing a significant influence on the choice of the Riemann solver. Copyright © 2015 John Wiley & Sons, Ltd.     

Anisotropic Turbulence and Coherent Structures in Eccentric Annular Channels

Eccentric annular pipe flows represent an ideal model for investigating inhomogeneous turbulent shear flows, where conditions of turbulence production and transport vary significantly within the cross-section. Moreover recent works have proven that in geometries characterized by the presence of a narrow gap, large-scale coherent structures are present. The eccentric annular channel represents, in the opinion of the present authors, the prototype of these geometries. The aim of the present work is to verify the capability of a numerical methodology to fully reproduce the main features of the flow field in this geometry, to verify and characterize the presence of large-scale coherent structures, to examine their behavior at different Reynolds numbers and eccentricities and to analyze the anisotropy associated to these structures. The numerical approach is based upon LES, boundary fitted coordinates and a fractional step algorithm. A dynamic Sub Grid Scale (SGS) model suited for this numerical environment has been implemented and tested. An additional interest of this work is therefore in the approach employed itself, considering it as a step into the development of an effective LES methodology for flows in complex channel geometries. Agreement with previous experimental and DNS results has been found good overall for the streamwise velocity, shear stress and the rms of the velocity components. The instantaneous flow field presented large-scale coherent structures in the streamwise direction at low Reynolds numbers, while these are absent or less dominant at higher Reynolds and low eccentricity. After Reynolds averaging is performed over a long integration time the existence of secondary flows in the cross session is proven. Their shape is found to be constant over the Reynolds range surveyed, and dependent on the geometric parameters. The effect of secondary flows on anisotropy is studied over an extensive Reynolds range through invariant analysis. Additional insight on the mechanics of turbulence in this geometry is obtained.     

Large eddy simulation of gas-particle turbulent channel flow with momentum exchange between the phases

This paper presents results of a large eddy simulation (LES) combined with Lagrangian particle tracking and a point-force approximation for the feedback effect of particles on the downward turbulent gaseous flow in a vertical channel. The LES predictions are compared with the results obtained by direct numerical simulation (DNS) of a finer computational mesh. A parametric study is conducted for particles with two response times in simulations with and without streamwise gravitational settling and elastic, binary interparticle collisions. It is shown that the classical and the dynamic Smagorinsky turbulence models adequately predict the particle-induced changes in the mean streamwise velocity and the Reynolds stresses of the carrier phase for the range of parameters studied. However, the largest discrepancies between the LES and DNS results are found in the cases of particle-laden flows. Conditional sampling of the instantaneous resolved flow fields indicates that the mechanisms by which particles directly oppose the production of momentum and vorticity of the organized fluid motions are also observed in the LES results. However, the geometric features of the near-wall quasistreamwise vortices are overestimated by the use of both turbulence models compared to the DNS predictions.     

Evidences of persisting thermal structures in Couette flows

DNS of passive thermal turbulent Couette flow at several friction Reynolds numbers (180, 250, and 500), and the Prandtl number of air are presented. The time averaged thermal flow shows the existence of long and wide thermal structures never described before in Couette flows. These thermal structures, named CTFS (Couette Thermal Flow Superstructures), are defined as coherent regions of hot and cold temperature fluctuations. They are intrinsically linked to the velocity structures present in Couette flows. Two different 2D symmetries can be recognized, which get stronger with the Reynolds number. These structures do not affect the mean flow or mean quantities as the Nusselt number. However, turbulent intensities and thermal fluxes depend on the width of the structures, mainly far from the walls. Since the width of the structures is related to the channel width, the statistics of thermal Couette flow are to some point box-dependent.     

An eddy viscosity model with near‐wall modifications

An extended version of the isotropic k–ε model is proposed that accounts for the distinct effects of low‐Reynolds number (LRN) and wall proximity. It incorporates a near‐wall correction term to amplify the level of dissipation in nonequilibrium flow regions, thus reducing the kinetic energy and length scale magnitudes to improve prediction of adverse pressure gradient flows, involving flow separation and reattachment. The eddy viscosity formulation maintains the positivity of normal Reynolds stresses and the Schwarz' inequality for turbulent shear stresses. The model coefficients/functions preserve the anisotropic characteristics of turbulence. The model is validated against a few flow cases, yielding predictions in good agreement with the direct numerical simulation (DNS) and experimental data. Comparisons indicate that the present model is a significant improvement over the standard eddy viscosity formulation. Copyright © 2005 John Wiley & Sons, Ltd.