Large-eddy simulation of flow over a cylinder with from to : a skin-friction perspective

2017 ◽  
Vol 820 ◽  
pp. 121-158 ◽  
Author(s):  
W. Cheng ◽  
D. I. Pullin ◽  
R. Samtaney ◽  
W. Zhang ◽  
W. Gao
Keyword(s):  
Boundary Layer ◽  
Friction Coefficient ◽  
Large Eddy Simulation ◽  
Skin Friction ◽  
Flow Separation ◽  
Skin Friction Coefficient ◽  
Cylinder Surface ◽  
Mean Flow ◽  
Eddy Simulation ◽  
Large Eddy

We present wall-resolved large-eddy simulations (LES) of flow over a smooth-wall circular cylinder up to$Re_{D}=8.5\times 10^{5}$, where$Re_{D}$is Reynolds number based on the cylinder diameter$D$and the free-stream speed$U_{\infty }$. The stretched-vortex subgrid-scale (SGS) model is used in the entire simulation domain. For the sub-critical regime, six cases are implemented with$3.9\times 10^{3}\leqslant Re_{D}\leqslant 10^{5}$. Results are compared with experimental data for both the wall-pressure-coefficient distribution on the cylinder surface, which dominates the drag coefficient, and the skin-friction coefficient, which clearly correlates with the separation behaviour. In the super-critical regime, LES for three values of$Re_{D}$are carried out at different resolutions. The drag-crisis phenomenon is well captured. For lower resolution, numerical discretization fluctuations are sufficient to stimulate transition, while for higher resolution, an applied boundary-layer perturbation is found to be necessary to stimulate transition. Large-eddy simulation results at$Re_{D}=8.5\times 10^{5}$, with a mesh of$8192\times 1024\times 256$, agree well with the classic experimental measurements of Achenbach (J. Fluid Mech., vol. 34, 1968, pp. 625–639) especially for the skin-friction coefficient, where a spike is produced by the laminar–turbulent transition on the top of a prior separation bubble. We document the properties of the attached-flow boundary layer on the cylinder surface as these vary with$Re_{D}$. Within the separated portion of the flow, mean-flow separation–reattachment bubbles are observed at some values of$Re_{D}$, with separation characteristics that are consistent with experimental observations. Time sequences of instantaneous surface portraits of vector skin-friction trajectory fields indicate that the unsteady counterpart of a mean-flow separation–reattachment bubble corresponds to the formation of local flow-reattachment cells, visible as coherent bundles of diverging surface streamlines.

scholarly journals Large-eddy simulation of the zero-pressure-gradient turbulent boundary layer up to Reθ = O(1012)

2011 ◽  
Vol 686 ◽  
pp. 507-533 ◽  
Author(s):  
M. Inoue ◽  
D. I. Pullin
Keyword(s):  
Boundary Layer ◽  
Friction Coefficient ◽  
Turbulent Boundary Layer ◽  
Large Eddy Simulation ◽  
Pressure Gradient ◽  
Skin Friction Coefficient ◽  
Smooth Wall ◽  
Eddy Simulation ◽  
Large Eddy ◽  
Near Wall

AbstractA near-wall subgrid-scale (SGS) model is used to perform large-eddy simulation (LES) of the developing, smooth-wall, zero-pressure-gradient flat-plate turbulent boundary layer. In this model, the stretched-vortex, SGS closure is utilized in conjunction with a tailored, near-wall model designed to incorporate anisotropic vorticity scales in the presence of the wall. Large-eddy simulations of the turbulent boundary layer are reported at Reynolds numbers ${\mathit{Re}}_{\theta } $ based on the free-stream velocity and the momentum thickness in the range ${\mathit{Re}}_{\theta } = 1{0}^{3} \text{{\ndash}} 1{0}^{12} $. Results include the inverse square-root skin-friction coefficient, $ \sqrt{2/ {C}_{f} } $, velocity profiles, the shape factor $H$, the von Kármán ‘constant’ and the Coles wake factor as functions of ${\mathit{Re}}_{\theta } $. Comparisons with some direct numerical simulation (DNS) and experiment are made including turbulent intensity data from atmospheric-layer measurements at ${\mathit{Re}}_{\theta } = O(1{0}^{6} )$. At extremely large ${\mathit{Re}}_{\theta } $, the empirical Coles–Fernholz relation for skin-friction coefficient provides a reasonable representation of the LES predictions. While the present LES methodology cannot probe the structure of the near-wall region, the present results show turbulence intensities that scale on the wall-friction velocity and on the Clauser length scale over almost all of the outer boundary layer. It is argued that LES is suggestive of the asymptotic, infinite Reynolds number limit for the smooth-wall turbulent boundary layer and different ways in which this limit can be approached are discussed. The maximum ${\mathit{Re}}_{\theta } $ of the present simulations appears to be limited by machine precision and it is speculated, but not demonstrated, that even larger ${\mathit{Re}}_{\theta } $ could be achieved with quad- or higher-precision arithmetic.

Wall-Modeled Large-Eddy Simulations of Flows With Curvature and Mild Separation

2008 ◽  
Vol 130 (10) ◽  
Author(s):  
Senthilkumaran Radhakrishnan ◽  
Ugo Piomelli ◽  
Anthony Keating
Keyword(s):  
Friction Coefficient ◽  
Large Eddy Simulation ◽  
Skin Friction ◽  
Outer Region ◽  
Skin Friction Coefficient ◽  
Stochastic Forcing ◽  
Mean Velocity ◽  
Eddy Simulation ◽  
Buffer Region ◽  
Large Eddy

The performance of wall-modeled large-eddy simulation (WMLES) based on hybrid models, in which the inner region is modeled by Reynolds-averaged Navier–Stokes (RANS) equation and the outer region is resolved by large-eddy simulation (LES), can make the application of LES attainable at high Reynolds numbers. In previous work by various authors, it was found that in most cases a buffer region exists between the RANS and LES zones, in which the velocity gradient is too high; this leads to an inaccurate prediction of the skin-friction coefficient. Artificially perturbing the RANS∕LES interface has been demonstrated to be effective in removing the buffer region. In this work, WMLES has been performed with stochastic forcing at the interface, following the previous work by our group on two nonequilibrium complex flows. From the two flows studied, we conclude that the application of stochastic forcing results in improvements in the prediction of the skin-friction coefficient in the equilibrium regions of these flows, a better agreement with the experiments of the Reynolds stresses in the adverse pressure gradient and the recovery region, and a good agreement of the mean velocity field with experiments in the separation region. Some limitations of this method, especially in terms of CPU requirements, will be discussed.

High-resolution atmospheric boundary layer simulations using WRF-LES

2020 ◽  
Author(s):  
Gokhan Kirkil
Keyword(s):  
Boundary Layer ◽  
Large Eddy Simulation ◽  
Cross Sections ◽  
Wrf Model ◽  
Field Measurements ◽  
Mean Flow ◽  
Spatial And Temporal Scales ◽  
Eddy Simulation ◽  
Wide Range ◽  
Large Eddy

<p>WRF model provides a potentially powerful framework for coupled simulations of flow covering a wide range of<br>spatial and temporal scales via a successive grid nesting capability. Nesting can be repeated down to turbulence<br>solving large eddy simulation (LES) scales, providing a means for significant improvements of simulation of<br>turbulent atmospheric boundary layers. We will present the recent progress on our WRF-LES simulations of<br>the Perdigao Experiment performed over mountainous terrain. We performed multi-scale simulations using<br>WRF’s different Planetary Boundary Layer (PBL) parameterizations as well as Large Eddy Simulation (LES)<br>and compared the results with the detailed field measurements. WRF-LES model improved the mean flow field<br>as well as second-order flow statistics. Mean fluctuations and turbulent kinetic energy fields from WRF-LES<br>solution are investigated in several cross-sections around the hill which shows good agreement with measurements.</p>

PIV Study of Shallow Open Channel Flow Over d- and k-Type Transverse Ribs

2007 ◽  
Vol 129 (8) ◽  
pp. 1058-1072 ◽  
Author(s):  
M. F. Tachie ◽  
K. K. Adane
Keyword(s):  
Boundary Layer ◽  
Friction Coefficient ◽  
Skin Friction ◽  
Cross Section ◽  
Open Channel ◽  
Skin Friction Coefficient ◽  
Mixing Length ◽  
Mean Flow ◽  
Turbulent Statistics ◽  
The Mean

A particle image velocimetry was used to study shallow open channel turbulent flow over d-type and k-type transverse ribs of square, circular, and semi-circular cross sections. The ratio of boundary layer thickness to depth of flow varied from 50% to 90%. The mean velocities and turbulent quantities were evaluated at the top plane of the ribs to characterize interaction between the cavities and overlying boundary layer. It was found that the overlying boundary layer interacts more strongly with k-type cavities than observed for d-type cavities. The profiles of the mean velocities and turbulent statistics were then spatially averaged over a pitch, and these profiles were used to study the effects of rib type and cross section on the flow field. The mean velocity gradients were found to be non-negligible across the boundary layer, and the implications of this observation for momentum transport, eddy viscosity, and mixing length distributions are discussed. The results show that the skin friction coefficient, Reynolds stresses and mixing length distributions are independent of rib cross section for d-type. For the k-type ribs, significant variations in skin friction coefficient values, mean flow, and turbulence fields are observed between square ribs and circular/semi-circular ribs.

Large Eddy Simulation of Flow and Heat Transfer in the 180‐Deg Bend Region of a Stationary Gas Turbine Blade Ribbed Internal Cooling Duct

2005 ◽  
Vol 128 (4) ◽  
pp. 763-771 ◽  
Author(s):  
Evan A. Sewall ◽  
Danesh K. Tafti
Keyword(s):  
Heat Transfer ◽  
Large Eddy Simulation ◽  
Friction Factor ◽  
Flow Separation ◽  
Mean Flow ◽  
Dean Vortices ◽  
Heat Transfer Augmentation ◽  
Eddy Simulation ◽  
High Heat Transfer ◽  
Large Eddy

Large eddy simulation of the 180 deg bend in a stationary ribbed duct is presented. The domain studied includes three ribs upstream of the bend region and three ribs downstream of the bend with an outflow extension added to the end, using a total of 8.4 million cells. Two cases are compared to each other: one includes a rib in the bend and the other does not. The friction factor, mean flow, turbulence, and heat transfer are compared in the two cases to help explain the benefits and disadvantages of the wide number of flow effects seen in the bend, including flow separation at the tip of the dividing wall, counter-rotating Dean vortices, high heat transfer at areas of flow impingement, and flow separation at the upstream and downstream corners of the bend. Mean flow results show a region of separated flow at the tip of the dividing region in the case with no rib in the bend, but no separation region is observed in the case with a rib. A pair of counter-rotating Dean vortices in the middle of the bend is observed in both cases. Turbulent kinetic energy profiles show a 30% increase in the midplane of the bend when the rib is added. High gradients of heat transfer augmentation are observed on the back wall and downstream outside wall, where mean flow impingement occurs. This heat transfer is increased with the presence of a rib. Including a rib in the bend increases the friction factor in the bend by 80%, and it increases the heat transfer augmentation by approximately 20%, resulting in a trade-off between pressure drop and heat transfer.

Turbulence Modeling Using Z-F, RSM, LES and WMLES for Flow Analysis in Z-Shape Ducts

2020 ◽  
Author(s):  
Mohammed Karbon ◽  
Ahmad K. Sleiti
Keyword(s):  
Experimental Data ◽  
Steady State ◽  
Large Eddy Simulation ◽  
Flow Separation ◽  
Temporal Variations ◽  
Mean Flow ◽  
Mean Velocity ◽  
Eddy Simulation ◽  
Large Eddy ◽  
Rsm Model

Abstract Turbulent flow in Z-shape duct configuration is investigated and analyzed using Reynolds Stress Model (RSM), Large Eddy Simulation (LES), ζ-f Model, and Wall-Modeled Large Eddy Simulation (WMLES). The results are validated and compared to experimental data. Both RSM and ζ-f models are based on steady-state RANS solutions, while LES and WMLES models account for temporal variations transient behavior of the flow turbulence. The focus was on regions where RSM has over or under predicted the flow and regions where there are flow separations and high turbulence. LES simulation results have shown under-prediction and over-prediction in the flow separation and re-attachment regions. It is found that the turbulent kinetic energy production in ζ equation is much easier to reproduce accurately than other models. Both mean velocity gradient and local turbulent stress terms are also much easier to resolve properly. The current research has found that ζ-f model not only takes less time to complete the simulation but also the mean flow velocity profile results are in better agreement with experimental data than RSM model despite both are coupled steady-state RANS. ζ-f model numerically resolved both the flow separation and re-attachment regions better than RSM model. WMLES model is employed to investigate the SGS model impact on the small eddies dissipated from the large eddies. Such WMLES model produces much better results than the LES model, however the SGS viscosity damps the energy of the flow.

Large-eddy simulation of flow over a rotating cylinder: the lift crisis at

2018 ◽  
Vol 855 ◽  
pp. 371-407 ◽  
Author(s):  
W. Cheng ◽  
D. I. Pullin ◽  
R. Samtaney
Keyword(s):  
Boundary Layer ◽  
Large Eddy Simulation ◽  
Large Scale ◽  
Free Stream ◽  
Lift Coefficient ◽  
Small Scale ◽  
Cylinder Surface ◽  
Eddy Simulation ◽  
Large Eddy ◽  
Attached Flow

We present wall-resolved large-eddy simulation (LES) of flow with free-stream velocity $\boldsymbol{U}_{\infty }$ over a cylinder of diameter $D$ rotating at constant angular velocity $\unicode[STIX]{x1D6FA}$ , with the focus on the lift crisis, which takes place at relatively high Reynolds number $Re_{D}=U_{\infty }D/\unicode[STIX]{x1D708}$ , where $\unicode[STIX]{x1D708}$ is the kinematic viscosity of the fluid. Two sets of LES are performed within the ( $Re_{D}$ , $\unicode[STIX]{x1D6FC}$ )-plane with $\unicode[STIX]{x1D6FC}=\unicode[STIX]{x1D6FA}D/(2U_{\infty })$ the dimensionless cylinder rotation speed. One set, at $Re_{D}=5000$ , is used as a reference flow and does not exhibit a lift crisis. Our main LES varies $\unicode[STIX]{x1D6FC}$ in $0\leqslant \unicode[STIX]{x1D6FC}\leqslant 2.0$ at fixed $Re_{D}=6\times 10^{4}$ . For $\unicode[STIX]{x1D6FC}$ in the range $\unicode[STIX]{x1D6FC}=0.48{-}0.6$ we find a lift crisis. This range is in agreement with experiment although the LES shows a deeper local minimum in the lift coefficient than the measured value. Diagnostics that include instantaneous surface portraits of the surface skin-friction vector field $\boldsymbol{C}_{\boldsymbol{f}}$ , spanwise-averaged flow-streamline plots, and a statistical analysis of local, near-surface flow reversal show that, on the leeward-bottom cylinder surface, the flow experiences large-scale reorganization as $\unicode[STIX]{x1D6FC}$ increases through the lift crisis. At $\unicode[STIX]{x1D6FC}=0.48$ the primary-flow features comprise a shear layer separating from that side of the cylinder that moves with the free stream and a pattern of oscillatory but largely attached flow zones surrounded by scattered patches of local flow separation/reattachment on the lee and underside of the cylinder surface. Large-scale, unsteady vortex shedding is observed. At $\unicode[STIX]{x1D6FC}=0.6$ the flow has transitioned to a more ordered state where the small-scale separation/reattachment cells concentrate into a relatively narrow zone with largely attached flow elsewhere. This induces a low-pressure region which produces a sudden decrease in lift and hence the lift crisis. Through this process, the boundary layer does not show classical turbulence behaviour. As $\unicode[STIX]{x1D6FC}$ is further increased at constant $Re_{D}$ , the localized separation zone dissipates with corresponding attached flow on most of the cylinder surface. The lift coefficient then resumes its increasing trend. A logarithmic region is found within the boundary layer at $\unicode[STIX]{x1D6FC}=1.0$ .

Sign in / Sign up

Export Citation Format

Share Document