Computational Fluid Dynamics Study for Improvement of Prediction of Various Thermally Stratified Turbulent Boundary Layers

2017 ◽  
Vol 139 (5) ◽  
Author(s):  
Hirofumi Hattori ◽  
Tomoya Houra ◽  
Amane Kono ◽  
Shota Yoshikawa

The objectives of this study are to reconstruct a turbulence model of both the large Eddy simulation (LES) and the Reynolds-averaged Navier–Stokes simulation (RANS) which can predict wind synopsis in various thermally stratified turbulent boundary layers over any obstacles. Hence, the direct numerical simulation (DNS) of various thermally stratified turbulent boundary layers with/without forward-step, two-dimensional block, or two-dimensional hill is carried out in order to obtain detailed turbulent statistics for the construction of a database for the evaluation of a turbulence model. Also, DNS clearly reveals the characteristics of various thermally stratified turbulent boundary layers with/without forward-step, two-dimensional block, or two-dimensional hill. The turbulence models employed in LES and RANS are evaluated using the DNS database we obtained. In the LES, an evaluated turbulence model gives proper predictions, but the quantitative agreement of Reynolds shear stress with DNS results is difficult to predict. On the other hand, the nonlinear eddy diffusivity turbulence models for Reynolds stress and turbulent heat flux are also evaluated using DNS results of various thermally stratified turbulent boundary layers over a forward-step in which the turbulence models are evaluated using an a priori method. Although the evaluated models do not make it easy to properly predict the Reynolds shear stresses in all cases, the turbulent heat fluxes can be qualitatively predicted by the nonlinear eddy diffusivity for a heat turbulence model. Therefore, the turbulence models of LES and RANS should be improved in order to adequately predict various thermally stratified turbulent boundary layers over an obstacle.

Author(s):  
H. Hattori ◽  
Y. Nagano

Direct numerical simulations (DNS) of boundary layers with various thermal stratifications are carried out to investigate the turbulent structures of these flows. The present DNSs quantitatively provide the characteristics of thermally stratified turbulent boundary layers. In particular, the counter gradient diffusion phenomenon is found in a strong, stable stratified boundary layer. On the other hand, in order to adequately predict turbulent boundary layers with various thermal stratifications, an appropriate turbulence model should be employed in the calculation. Thus, using a database obtained by DNS, the strict assessment of turbulent heat transfer model is made so as to construct a reliable advanced turbulence model. The results of in-depth turbulent model evaluation are indicated, in which we have explored the prediction potential of the proposed nonlinear eddy diffusivity models for momentum and heat in both stable and unstable stratified boundary layers.


1968 ◽  
Vol 19 (1) ◽  
pp. 1-19 ◽  
Author(s):  
H. McDonald

SummaryRecently two authors, Nash and Goldberg, have suggested, intuitively, that the rate at which the shear stress distribution in an incompressible, two-dimensional, turbulent boundary layer would return to its equilibrium value is directly proportional to the extent of the departure from the equilibrium state. Examination of the behaviour of the integral properties of the boundary layer supports this hypothesis. In the present paper a relationship similar to the suggestion of Nash and Goldberg is derived from the local balance of the kinetic energy of the turbulence. Coupling this simple derived relationship to the boundary layer momentum and moment-of-momentum integral equations results in quite accurate predictions of the behaviour of non-equilibrium turbulent boundary layers in arbitrary adverse (given) pressure distributions.


Author(s):  
Ralph J. Volino ◽  
Terrence W. Simon

A technique called “octant analysis” was used to examine the eddy structure of turbulent and transitional heated boundary layers on flat and curved surfaces. The intent was to identify important physical processes that play a role in boundary layer transition on flat and concave surfaces. Octant processing involves the partitioning of flow signals into octants based on the instantaneous signs of the fluctuating temperature, t′; streamwise velocity, u′; and cross-stream velocity, v′. Each octant is associated with a particular eddy motion. For example, u′<0, v′>0, t′>0 is associated with an ejection or “burst” of warm fluid away from a heated wall. Within each octant, the contribution to various quantities of interest (such as the turbulent shear stress, −u′v′, or the turbulent heat flux, v′t′) can be computed. By comparing and contrasting the relative contributions from each octant, the importance of particular types of motion can be determined. If the data within each octant is further segregated based on the magnitudes of the fluctuating components so that minor events are eliminated, the relative importance of particular types of motion to the events that are important can also be discussed. In fully-developed, turbulent boundary layers along flat plates, trends previously reported in the literature were confirmed. A fundamental difference was observed in the octant distribution between the transitional and fully-turbulent boundary layers, however, showing incomplete mixing and a lesser importance of small scales in the transitional boundary layer. Such observations were true on both flat and concave walls. The differences are attributed to incomplete development of the turbulent kinetic energy cascade in transitional flows. The findings have potential application to modelling, suggesting the utility of incorporating multiple length scales in transition models.


2019 ◽  
Vol 880 ◽  
pp. 239-283 ◽  
Author(s):  
Christoph Wenzel ◽  
Tobias Gibis ◽  
Markus Kloker ◽  
Ulrich Rist

A direct numerical simulation study of self-similar compressible flat-plate turbulent boundary layers (TBLs) with pressure gradients (PGs) has been performed for inflow Mach numbers of 0.5 and 2.0. All cases are computed with smooth PGs for both favourable and adverse PG distributions (FPG, APG) and thus are akin to experiments using a reflected-wave set-up. The equilibrium character allows for a systematic comparison between sub- and supersonic cases, enabling the isolation of pure PG effects from Mach-number effects and thus an investigation of the validity of common compressibility transformations for compressible PG TBLs. It turned out that the kinematic Rotta–Clauser parameter $\unicode[STIX]{x1D6FD}_{K}$ calculated using the incompressible form of the boundary-layer displacement thickness as length scale is the appropriate similarity parameter to compare both sub- and supersonic cases. Whereas the subsonic APG cases show trends known from incompressible flow, the interpretation of the supersonic PG cases is intricate. Both sub- and supersonic regions exist in the boundary layer, which counteract in their spatial evolution. The boundary-layer thickness $\unicode[STIX]{x1D6FF}_{99}$ and the skin-friction coefficient $c_{f}$, for instance, are therefore in a comparable range for all compressible APG cases. The evaluation of local non-dimensionalized total and turbulent shear stresses shows an almost identical behaviour for both sub- and supersonic cases characterized by similar $\unicode[STIX]{x1D6FD}_{K}$, which indicates the (approximate) validity of Morkovin’s scaling/hypothesis also for compressible PG TBLs. Likewise, the local non-dimensionalized distributions of the mean-flow pressure and the pressure fluctuations are virtually invariant to the local Mach number for same $\unicode[STIX]{x1D6FD}_{K}$-cases. In the inner layer, the van Driest transformation collapses compressible mean-flow data of the streamwise velocity component well into their nearly incompressible counterparts with the same $\unicode[STIX]{x1D6FD}_{K}$. However, noticeable differences can be observed in the wake region of the velocity profiles, depending on the strength of the PG. For both sub- and supersonic cases the recovery factor was found to be significantly decreased by APGs and increased by FPGs, but also to remain virtually constant in regions of approximated equilibrium.


2013 ◽  
Vol 726 ◽  
pp. 196-225 ◽  
Author(s):  
Jin Lee ◽  
Seo Yoon Jung ◽  
Hyung Jin Sung ◽  
Tamer A. Zaki

AbstractDirect numerical simulations (DNS) of turbulent boundary layers over isothermally heated walls were performed, and the effect of viscosity stratification on the turbulence statistics and skin friction were investigated. An empirical relation for temperature-dependent viscosity for water was adopted. Based on the free-stream temperature (30°C), two wall temperatures (70°C and 99°C) were selected. In the heated flows, the turbulence energy diminishes in the buffer layer, but increases near the wall. The reduction in turbulence kinetic energy in the buffer layer is accompanied by smaller levels of Reynolds shear stresses and, hence, weaker turbulence production. The enhanced turbulence energy near the wall is attributed to enhanced transfer of energy via additional diffusion-like terms due to the viscosity stratification. Despite the lower fluid viscosity near the wall, dissipation is also increased owing to the augmented near-wall fine-scale motion. Wall heating results in reduction in the skin-friction coefficient by up to 26 %. An evaluation of the different contributions to the skin friction demonstrates that drag reduction is primarily due to the changes in the Reynolds shear stresses across the boundary layer. Quadrant and octant analyses showed that ejections (Q2) and sweeps (Q4) are significantly reduced, a result further supported by an examination of outer vortical structures from linear stochastic estimation of the ejection events and spanwise vortices.


Sign in / Sign up

Export Citation Format

Share Document