LES Study of Turbulent Boundary Layer Over a Smooth and a Rough 2D Hill Model

Large eddy simulation (LES) is carried out to investigate the turbulent boundary-layer flows over a hill-shaped model with a steep or relatively moderate slope at moderately high Reynolds numbers (Re = O(10³)) defined by the hill height and the velocity at the hill height. The study focuses on the effects of surface roughness and curvature. For Sub-grid Scale (SGS) modeling of LES, both the dynamic Smagorinsky model (DSM) and the dynamic mixed model (DMM) are applied. The behavior of the separated shear layer and the vortex motion are affected by the oncoming turbulence, such that the shear layer comes close to the ground surface, or the size of a separation region becomes small because of the earlier instability of the separated shear layer. Appropriate measures are required to generate the inflow turbulence. The methods of Lund et al. (J. Comput. Phys., 140:233–258, 1998) and Nozawa and Tamura (J. Wind Eng. Ind. Aerodyn., 90:1151–1162, 2002; The 4th European and African Conference on Wind Engineering, 1–6, 2005) are employed to simulate the smooth- and rough-wall turbulent boundary layers in order to generate time-sequential data of inflow turbulence. This paper discusses the unsteady phenomena of the wake flows over the smooth and rough 2D hill-shaped obstacles and aims to clarify the roughness effects on the flow patterns and the turbulence statistics distorted by the hill. Numerical validation is conducted by comparing the simulation results with wind tunnel experiment data for the same hill shape at almost the same Re. The applicability of DSM and DMM are discussed, focusing on the recirculation region behind a steep hill.     

A Hybrid RANS/LES Simulation of Turbulent Channel Flow

Hybrid models combining large eddy simulation (LES) with Reynolds-averaged Navier–Stokes (RANS) simulation are expected to be useful for wall modeling in the LES of high Reynolds number flows. Some hybrid simulations of turbulent channel flow have a common defect; the mean velocity profile has a mismatch between the RANS and LES regions due to a steep velocity gradient at the interface. This mismatch is reproduced and examined using a simple hybrid model; the Smagorinsky model is switched to a RANS model increasing the filter width. It is suggested that a rapid spatial variation in the eddy viscosity is responsible for an underestimate of the grid-scale shear stress and for the steep velocity gradient. To reduce the mean velocity mismatch a new scheme is proposed; additional filtering is introduced to define two kinds of velocity components at the interface between the two regions. The two components are used to remove inconsistency in the velocity equations due to a rapid variation in the filter width. Using the new scheme, simulations of channel flow are carried out with the simple hybrid model. It is shown that the grid-scale shear stress becomes large enough and most of the mean velocity mismatch is removed. Simulations for higher Reynolds numbers are carried out with the k–ε model and the one-equation subgrid-scale model. Although it is necessary to improve the turbulence models and the treatment of the buffer region, the new scheme is shown to be effective for reducing the mismatch and to be useful for developing better hybrid simulations. Received 5 April 2002 and accepted 8 January 2003 Published online 25 March 2003 Communicated by M.Y. Hussaini     

The State of the Art of Hybrid RANS/LES Modeling for the Simulation of Turbulent Flows

This review presents the state of the art of hybrid RANS/LES modeling for the simulation of turbulent flows. After recalling the modeling used in RANS and LES methodologies, we propose in a first step a theoretical formalism developed in the spectral space that allows to unify the RANS and LES methods from a physical standpoint. In a second step, we discuss the principle of the hybrid RANS/LES methods capable of representing a RANS-type behavior in the vicinity of a solid boundary and an LES-type behavior far away from the wall boundary. Then, we analyze the principal hybrid RANS/LES methods usually used to perform numerical simulation of turbulent flows encountered in engineering applications. In particular, we investigate the very large eddy simulation (VLES), the detached eddy simulation (DES), the partially integrated transport modeling (PITM) method, the partially averaged Navier-Stokes (PANS) method, and the scale adaptive simulation (SAS) from a physical point of view. Finally, we establish the connection between these methods and more precisely, the link between PITM and PANS as well as DES and PITM showing that these methods that have been built by different ways, practical or theoretical manners have common points of comparison. It is the opinion of the author to consider that the most appropriate method for a particular application will depend on the expectations of the engineer and the computational resources the user is prepared to expend on the problem.     

Hybrid RANS/LES computations of plane impinging jets with DES and PANS models

The qualities of a DES (Detached Eddy Simulation) and a PANS (Partially-Averaged Navier–Stokes) hybrid RANS/LES model, both based on the k–ω RANS turbulence model of Wilcox (2008, “Formulation of the k–ω turbulence model revisited” AIAA J., 46: 2823–2838), are analysed for simulation of plane impinging jets at a high nozzle-plate distance (H/B = 10, Re = 13,500; H is nozzle-plate distance, B is slot width; Reynolds number based on slot width and maximum velocity at nozzle exit) and a low nozzle-plate distance (H/B = 4, Re = 20,000). The mean velocity field, fluctuating velocity components, Reynolds stresses and skin friction at the impingement plate are compared with experimental data and LES (Large Eddy Simulation) results. The k–ω DES model is a double substitution type, following Davidson and Peng (2003, “Hybrid LES–RANS modelling: a one-equation SGS model combined with a k–ω model for predicting recirculating flows” Int. J. Numer. Meth. Fluids, 43: 1003–1018). This means that the turbulent length scale is replaced by the grid size in the destruction term of the k-equation and in the eddy viscosity formula. The k–ω PANS model is derived following Girimaji (2006, “Partially-Averaged Navier–Stokes model for turbulence: a Reynolds-Averaged Navier–Stokes to Direct Numerical Simulation bridging method” J. Appl. Mech., 73: 413–421). The turbulent length scale in the PANS model is constructed from the total turbulent kinetic energy and the sub-filter dissipation rate. Both hybrid models change between RANS (Reynolds-Averaged Navier–Stokes) and LES based on the cube root of the cell volume. The hybrid techniques, in contrast to RANS, are able to reproduce the turbulent flow dynamics in the shear layers of the impacting jet. The change from RANS to LES is much slower however for the PANS model than for the DES model on fine enough grids. This delays the break-up process of the vortices generated in the shear layers with as a consequence that the DES model produces better results than the PANS model.     

DNS and Highly-Resolved LES of Heat and Mass Transfer in Two-Phase Counter-Current Condensing Flow

Flow, Turbulence and Combustion - A comprehensive study of direct-contact condensation heat transfer for turbulent, counter-current, liquid/vapour flow in a nearly horizontal channel at high...     

Compound Wall Treatment for RANS Computation of Complex Turbulent Flows and Heat Transfer

We present a generalised treatment of the wall boundary conditions for RANS computation of turbulent flows and heat transfer. The method blends the integration up to the wall (ItW) with the generalised wall functions (GWF) that include non-equilibrium effects. Wall boundary condition can thus be defined irrespective of whether the wall-nearest grid point lies within the viscous sublayer, in the buffer zone, or in the fully turbulent region. The computations with fine and coarse meshes of a steady and pulsating flow in a plane channel, in flow behind a backward-facing step and in a round impinging jet using the proposed compound wall treatment (CWT) are all in satisfactory agreement with the available experiments and DNS data. The method is recommended for computations of industrial flows in complex domains where it is difficult to generate a computational grid that will satisfy a priori either the ItW or WF prerequisites.     

DNS,LES and RANS of turbulent heat transfer in boundary layer with suddenly changing wall thermal conditions

The objectives of this study are to investigate a thermal field in a turbulent boundary layer with suddenly changing wall thermal conditions by means of direct numerical simulation (DNS), and to evaluate predictions of a turbulence model in such a thermal field, in which DNS of spatially developing boundary layers with heat transfer can be conducted using the generation of turbulent inflow data as a method. In this study, two types of wall thermal condition are investigated using DNS and predicted by large eddy simulation (LES) and Reynolds-averaged Navier–Stokes equation simulation (RANS). In the first case, the velocity boundary layer only develops in the entrance of simulation, and the flat plate is heated from the halfway point, i.e., the adiabatic wall condition is adopted in the entrance, and the entrance region of thermal field in turbulence is simulated. Then, the thermal boundary layer develops along a constant temperature wall followed by adiabatic wall. In the second case, velocity and thermal boundary layers simultaneously develop, and the wall thermal condition is changed from a constant temperature to an adiabatic wall in the downstream region. DNS results clearly show the statistics and structure of turbulent heat transfer in a constant temperature wall followed by an adiabatic wall. In the first case, the entrance region of thermal field in turbulence can be also observed. Thus, both the development and the entrance regions in thermal fields can be explored, and the effects upstream of the thermal field on the adiabatic region are investigated. On the other hand, evaluations of predictions by LES and RANS are conducted using DNS results. The predictions of both LES and RANS almost agree with the DNS results in both cases, but the predicted temperature variances near the wall by RANS give different results as compared with DNS. This is because the dissipation rate of temperature variance is difficult to predict by the present RANS, which is found by the evaluation using DNS results.     

Hybrid RANS/LES of flow and heat transfer in round impinging jets

Fluid flow and convective heat transfer predictions are presented of round impinging jets for several combinations of nozzle-plate distances H/D = 2, 6 and 13.5 (where D is the nozzle diameter) and Reynolds numbers Re = 5000, 23,000 and 70,000 with the newest version of the k-ω model of Wilcox (2008) and three hybrid RANS/LES models. In the RANS mode of the hybrid RANS/LES models, the k-ω model is recovered. Three formulations are considered to activate the LES mode. The first model is similar to the hybrid models of Davidson and Peng (2003) and Kok et al. (2004). The turbulent length scale is replaced by the grid size in the destruction term of the k-equation and in the definition of the RANS eddy viscosity. As grid size, a maximum measure of the hexahedral grid cell is used. The second model has the same k-equation, but the eddy viscosity is the minimum of the k-ω eddy viscosity and the Smagorinsky eddy viscosity, following a proposal by Batten et al. (2004). The Smagorinsky eddy viscosity is formed with the cube root of the cell volume. The third model has, again, the same k-equation, but has an eddy viscosity which is an intermediate between the eddy viscosities of the first and second models. This is reached by using the cube root of the cell volume in the eddy viscosity formula of the first model.The simulation results are compared with experimental data for the high Reynolds number cases Re = 23,000 and Re = 70,000 and LES data for the low-Reynolds number case Re = 5000. The Reynolds numbers are defined with the nozzle diameter and the bulk velocity at nozzle outlet. At low nozzle-plate distance (the impingement plate is in the core of the jet), turbulent kinetic energy is overpredicted by RANS in the stagnation flow region. This leads to overprediction of the heat transfer rate along the impingement plate in the impact zone. At high nozzle-plate distance (the impingement plate is in the mixed-out region of the jet), the turbulence mixing is underpredicted by RANS in the shear layer of the jet which gives a too high length of the jet core. This also results in overprediction of the heat transfer rate in the impingement zone caused by too big temperature gradients at impingement.All hybrid RANS/LES models are able to correct the heat transfer overprediction of the RANS model. For good predictions at low nozzle-plate distance, it is necessary to sufficiently resolve the formation and development of the near-wall vortices in the jet impingement region. At high nozzle-plate distance, the essence is to capture the evolution and breakup of the flow unsteadiness in the shear layer of the jet, so that accurate mean and fluctuating velocity profiles are obtained in the impingement region. Although the models have a quite different theoretical justification and generate a quite different eddy viscosity in some flow regions, their overall results are very comparable. The reason is that in zones that are crucial for the results, the models behave similarly.     

Non-Equilibrium Mixing of Turbulence Scales Using a Continuous Hybrid RANS/LES Approach: Application to the Shearless Mixing Layer

A new subgrid-scale (SGS) model based on partially integrated transport method (PITM) is applied to the case of a turbulent spectral non-equilibrium flow created by the mixing of two turbulence fields of differing scales: the shearless mixing layer. The method can be viewed as a continuous hybrid RANS/LES approach. In this model the SGS length scale is no longer given by the size of the discretization step, but is dynamically estimated using an additional transport equation for the dissipation rate. The results are compared to those corresponding to the classical model of Smagorinsky and to the experimental data of Veeravalli and Warhaft. A method for creating an anisotropic analytical pseudo-random field for inflow conditions is also proposed. This approach based on subgrid-scale transport modelling combined with anisotropic inlet conditions gives better results for the prediction of the shearless mixing layer.     

Dynamic Evaluation of Mesh Resolution and Its Application in Hybrid LES/RANS Methods

In this work, we investigate a resolution evaluation criterion based on the ratio between turbulent length-scales and grid spacing within the context of dynamic resolution evaluation in hybrid LES/RANS simulations. A modified version of the commonly used length-scale criterion is adopted. The modified length-scale criterion is evaluated for a plane channel flow and compared to the criterion based on two-point correlations. Simulation results show qualitative agreement between the two criteria and physical predictions from both resolution indicators. These observations are confirmed by simulations of flows over periodic hills. It is further demonstrated that the length-scale based criterion is relatively less sensitive on variation of model parameters compared to criteria based on resolved percentage of turbulent quantities. The improved resolution criterion is applied in a dual-mesh hybrid LES/RANS solver. Numerical simulations with the hybrid solver suggest that the interactions between the length-scale resolution indicator and the solution are moderate, and that favorable comparisons with benchmark results are obtained. In summary, we demonstrate that the modified length-scale based resolution indicator performs satisfactorily in both pure LES and hybrid simulations. Therefore, it is selected as a promising candidate to provide reliable predictions of resolution adequacy for individual cells in hybrid LES/RANS simulations.     

Nonsimilar Boundary Layer Analysis of Free Convection Heat Transfer Over a Vertical Cylinder in Bidisperse Porous Media

This work presents a boundary layer analysis for the free convection heat transfer from a vertical cylinder in bidisperse porous media with constant wall temperature. A boundary layer analysis and the two-velocity two-temperature formulation are used to derive the nonsimilar governing equations. The transformed governing equations are solved by the cubic spline collocation method to yield computationally efficient numerical solutions. The effects of inter-phase heat transfer parameter, modified thermal conductivity ratio, and permeability ratio on the heat transfer and flow characteristics are studied. Results show that an increase in the modified thermal conductivity ratio and the permeability ratio can effectively enhance the free convection heat transfer of the vertical cylinder in a bidisperse porous medium. Moreover, the thermal nonequilibrium effects are strong for low values of the inter-phase heat transfer parameter.     

DNS Analysis of Wall Heat Transfer and Combustion Regimes in a Turbulent Non-premixed Wall-jet Flame

Understanding the heat-release effects on the wall heat transfer in turbulent reacting flows, i.e. heat transfer with or without significant density variation, is essential for a wide variety of industrial flows, especially combustion problems. The present study focuses on the wall heat transfer and the near-wall reaction characteristics. The heat-release effects on the wall heat transfer and skin-friction coefficients are investigated using three-dimensional direct numerical simulations of a turbulent reacting wall-jet flow with and without heat release. Reductions in the skin-friction coefficient are observed in the exothermic case, compared to the isothermal one, and the underlying mechanism is explained. The absolute wall heat flux also increases, while the corresponding Nusselt number decreases with increasing heat release. Furthermore, the wall effects on the near-wall average burning rate are assessed. It is found that the isothermal cold wall results in an appreciable decrease of the burning rate in the exothermic cases. We observed indications that the wall increases the chances for the development of the premixed mode and its occurrence is very fast in the wall-normal direction.     

Some Features of a Turbulent Separated Flow and Heat Transfer Behind a Step and a Rib. 2. Heat Transfer in a Separated Flow

Results of an experimental study of heat transfer in a separated flow behind a step and a rib are presented. The influence of the obstacle height (H = 6–30 mm) on heat and mass transfer and the structure of the thermal boundary layer is studied. The features of heat transfer in recirculation and relaxation zones of the separated flow are analyzed, and the effect of separation on intensification and suppression of turbulent heat transfer is determined.     

Numerical Study of Natural Convection Heat Transfer in an Inclined Porous Cavity with Time-Periodic Boundary Conditions

Kalabin et al. (Numer. Heat Transfer A 47, 621-631, 2005) studied the unsteady natural convection for the sinusoidal oscillating wall temperature on one side wall and constant average temperature on the opposing side wall. The present article is on the unsteady natural convective heat transfer in an inclined porous cavity with similar temperature boundary conditions as those of Kalabin et al. The inclined angle of the cavity is varied from 0° to 80°. The flow field is modeled with the Brinkman-extended Darcy model. The combined effects of inclination angle of the enclosure and oscillation frequency of wall temperature are studied for Ra* = 10³, Da = 10⁻³, , and Pr=1. Some results are also obtained with the Darcy–Brinkman–Forchheimer model and Darcy’s law and are compared with the present Brinkman-extended Darcy model. The maximal heat transfer rate is attained at the oscillating frequency f = 46.7π and the inclined angle .     

Stability of the Boundary Layer with Internal Heat Release and Gas Injection through a Porous Wall

Fluid Dynamics - The problem of stability of a subsonic boundary layer is solved under the conditions of heat supply inside the boundary layer with injection of a homogeneous gas through a porous...     

Experimental and Numerical Study of the Scalar Turbulent/Non-Turbulent Interface Layer in a Jet Flow

Based on two large-eddy simulations (LES) of a non-reacting turbulent round jet with a nozzle based Reynolds number of 8,610 with the same configuration as the one that has recently been investigated experimentally (Gampert et al., 2012; J Fluid Mech, 2012; J Fluid Mech 724:337, 2013), we examine the scalar turbulent/non-turbulent (T/NT) interface layer in the mixture fraction field of the jet flow between ten and thirty nozzle diameters downstream. To this end, the LES—one with a coarse grid and one with a fine grid—are in a first step validated against the experimental data using the axial decay of the mean velocity and the mean mixture fraction as well as based on radial self-similar profiles of mean and root mean square values of these two quantities. Then, probability density functions (pdf) of the mixture fraction at various axial and radial positions are compared and the quality of the LES is discussed. In general, the LES results are consistent with the experimental data. However, in the flow region where the imprint of the T/NT interface layer is dominant in the mixture fraction pdf, discrepancies are observed. In a next step, statistics of the T/NT interface layer are studied, where a satisfactory agreement for the pdf of the location of the interface layer from the higher resolved LES with the experimental data is observed, while the one with the coarse grid exhibits considerable deviations. Finally, the mixture fraction profile across the interface is investigated where the same trend as for the pdf of the location is present. In particular, it is found that the sharp interface that is present in experimental studies (Gampert et al., J Fluid Mech, 2013; Westerweel et al., J Fluid Mech 631:199, 2009) is less distinct in the LES results and rather diffused in radial direction outside of the T/NT interface layer.