scholarly journals Large eddy simulation of the velocity-intermittency structure for flow over a field of symmetric dunes

2016 ◽  
Vol 805 ◽  
pp. 656-685 ◽  
Author(s):  
Christopher J. Keylock ◽  
Kyoungsik S. Chang ◽  
George S. Constantinescu

Owing to their frequent occurrence in the natural environment, there has been significant interest in refining our understanding of flow over dunes and other bedforms. Recent work in this area has focused, in particular, on their shear-layer characteristics and the manner by which flow structures are generated. However, field-based studies, are reliant on single-, or multi-point measurements, rather than delimiting flow structures from the velocity gradient tensor, as is possible in numerical work. Here, we extract pointwise time series from a well-resolved large eddy simulation as a means to connect these two approaches. The at-a-point analysis technique is termed the velocity-intermittency quadrant method and relates the fluctuating, longitudinal velocity, $u_{1}^{\prime }(t)$, to its fluctuating pointwise Hölder regularity, $\unicode[STIX]{x1D6FC}_{1}^{\prime }(t)$. Despite the difference in boundary conditions, our results agree very well with previous experiments that show the importance, in the region above the dunes, of a quadrant 3 ($u_{1}^{\prime }<0$, $\unicode[STIX]{x1D6FC}_{1}^{\prime }<0$) flow configuration. Our higher density of sampling beneath the shear layer and close to the bedforms relative to experimental work reveals a negative correlation between $u_{1}^{\prime }(t)$ and $\unicode[STIX]{x1D6FC}_{1}^{\prime }(t)$ in this region. This consists of two distinct layers, with quadrant 4 ($u_{1}^{\prime }>0$, $\unicode[STIX]{x1D6FC}_{1}^{\prime }<0$) dominant near the wall and quadrant 2 ($u_{1}^{\prime }<0$, $\unicode[STIX]{x1D6FC}_{1}^{\prime }>0$) dominant close to the lower part of the separated shear layer. These results are consistent with a near-wall advection of vorticity into a region downstream of a temporarily foreshortened reattachment region, and the entrainment of slow moving and quiescent fluid into a faster, more turbulent shear layer. A comparison of instantaneous vorticity fields to the velocity-intermittency analysis shows how the pointwise results reflect larger-scale organisation of the flow. We illustrate this using results from two instantaneous datasets. In the former, extreme velocity-intermittency events corresponding to a foreshortened recirculation region (and high pressures on the stoss slope of the dune immediately downstream) arise, and the development of intense flow structures occurs as a consequence. In the other case, development of a ‘skimming flow’ with relatively little exchange between the inner and outer regions results in exceedances because of the coherence associated with this high velocity, high turbulence outer region. Thus, our results shed further light on the characteristics of dune flow in the near-wall region and, importantly for field-based research, show that useful information on flow structure can be obtained from single-point single velocity component measurements.

2020 ◽  
Vol 24 (5 Part A) ◽  
pp. 2793-2808
Author(s):  
Hussain Al-Kayiem ◽  
Desmond Lim ◽  
Jundika Kurnia

Sharkskin-inspired riblets are widely adopted as a passive method for drag reduc?tion of flow over surfaces. In this research, large eddy simulation of turbulent flow over riblet-structured surface in a rectangular channel domain were performed at various Reynolds numbers, ranging from 4200-10000, to probe the resultant drag change, compared to smooth surface. The changes of mean streamwise velocity gradient in wall-normal direction at varied locations around riblet structures were also investigated to reduce mechanisms of streamlined riblet in reducing drag. The computational model is validated by comparing the simulation results against analytical and experimental data, for both smooth and riblet surfaces. Results in?dicating that the performance of the proposed streamlined riblet shows 7% drag reduction, as maximum, which is higher than the performance of L-shaped riblet with higher wetted surface area. The mean velocity profile analysis indicates that the streamlined riblet structures help to reduce longitudinal averaged velocity component rate in the normal to surface direction of near-wall region which leads to laminarization process as fluid-flows over riblet structures.


2014 ◽  
Vol 759 ◽  
pp. 579-611 ◽  
Author(s):  
Jae Hwan Jung ◽  
Hyun Sik Yoon

AbstractWe consider a twisted cylinder that was designed by rotating the elliptic cross-section along the spanwise direction, resulting in a passive control. The flow over the twisted cylinder is investigated at a subcritical Reynolds number (Re) of 3000 using large eddy simulation based on the finite volume method. For comparison, the flow past smooth and wavy cylinders is also calculated. The twisted cylinder achieves reductions of approximately 13 and 5 % in mean drag compared with smooth and wavy cylinders, respectively. In particular, the root mean square (r.m.s.) value of the lift fluctuation of the twisted cylinder shows a substantial decrease of approximately 96 % compared with the smooth cylinder. The shear layer of the twisted cylinder covering the recirculation region is more elongated than those of the smooth and wavy cylinders, and vortex shedding from the twisted cylinder is considerably suppressed. Consequently, the elongation of the shear layer from the body and the near disappearance of vortex shedding in the near wake with weak vortical strength contributes directly to the reduction of drag and lift oscillation. Various fundamental mechanisms that affect the flow phenomena, three-dimensional separation, pressure coefficient, vortex formation length and turbulent kinetic energy are examined systematically to demonstrate the effect of the twisted cylinder surface. In addition, for the twisted cylinder at $\mathit{Re}=3000$, the effect of the cross-sectional aspect ratio is investigated from 1.25 to 2.25 to find an optimal value that can reduce the drag and lift forces. Moreover, the effect of the Reynolds number on the aerodynamic characteristics is investigated in the range of $3\times 10^{3}\leqslant \mathit{Re}\leqslant 1\times 10^{4}$. We find that as Re increases, the mean drag and the r.m.s. lift coefficient of the twisted cylinder increase, and the vortex formation length decreases.


Author(s):  
Sachin S. Badarayani ◽  
Kyle D. Squires

Large Eddy Simulation (LES) of high-Reynolds-number wall-bounded turbulent flows is prohibitively expensive if the energy-containing eddies in the near-wall region are resolved. This motivates the use of wall-layer models in which an approximate solution of the near wall dynamics is bridged to an LES of the outer flow. The main interest of the present work are wall-modeling strategies based on Detached Eddy Simulation (DES). In these approaches, the near-wall solution is closed using a Reynolds-averaged Navier Stokes model with a subgrid closure applied to the outer flow. As is well known, the original DES formulation applied directly as a wall model results in a shift in the velocity profile, corresponding to an under-estimation of the skin friction. A new formulation is proposed in this contribution in which the wall-parallel components of the modeled stress are reduced in order to lower the influence of the model and increase the resolved stress. The effectiveness of the new model is evaluated via comparison against DES predictions using the original and recently-proposed versions of the method. The effect of grid resolution and model parameters are also assessed using computations of turbulent channel flow at a Reynolds number based on friction velocity and channel halfwidth of 5000. The predictions show that the anisotropic form of the model stress yields an improved prediction of the mean velocity profile in better agreement with the logarithmic law and with larger resolved stress in the near-wall region.


2011 ◽  
Vol 133 (11) ◽  
Author(s):  
Dalibor Cavar ◽  
Knud Erik Meyer

A large eddy simulation (LES) study of turbulent non-equilibrium boundary layer flow over 2D Bump, at comparatively low Reynolds number Reh=U∞h/ν=1950, was conducted. A well-known LES issue of obtaining and sustaining turbulent flow inside the computational domain at such low Re, is addressed by conducting a precursor calculation of the spatially developing boundary layer flow. Those results were subsequently used as turbulent inflow database for the main non-equilibrium boundary layer flow computation. The Sagaut (Rech. Aero., pp. 51–63, 1996) sub grid scale (SGS) turbulence model, based on a local estimate of the subgrid scale turbulent kinetic energy ksgs and implicit damping of turbulent SGS viscosity νt(sgs) in the near-wall region, was selected as a suitable basis for the present LES computations due to the fact that block structured MPI parallelized CFD code used in the current computations did not provide a direct possibility for wall-damping of, e.g., the Smagorinsky constant in the near-wall region. The grid utilized in the main calculation consisted of approximately 9.4 × 106 grid points and the boundary layer flow results obtained, regarding both mean flow profiles and turbulence quantities, showed a good agreement with the available laser Doppler anemometry (LDA) measurements. Analysis of the flow was directly able to identify and confirm the existence of internal layers at positions related to the vicinity of the upstream and downstream discontinuities in the surface curvature and also partially confirm a close interdependency between generation and evolution of internal layers and the abrupt changes in the skin friction, previously reported in the literature.


2009 ◽  
Vol 631 ◽  
pp. 281-309 ◽  
Author(s):  
D. CHUNG ◽  
D. I. PULLIN

We report large-eddy simulation (LES) of turbulent channel flow. This LES neither resolves nor partially resolves the near-wall region. Instead, we develop a special near-wall subgrid-scale (SGS) model based on wall-parallel filtering and wall-normal averaging of the streamwise momentum equation, with an assumption of local inner scaling used to reduce the unsteady term. This gives an ordinary differential equation (ODE) for the wall shear stress at every wall location that is coupled with the LES. An extended form of the stretched-vortex SGS model, which incorporates the production of near-wall Reynolds shear stress due to the winding of streamwise momentum by near-wall attached SGS vortices, then provides a log relation for the streamwise velocity at the top boundary of the near-wall averaged domain. This allows calculation of an instantaneous slip velocity that is then used as a ‘virtual-wall’ boundary condition for the LES. A Kármán-like constant is calculated dynamically as part of the LES. With this closure we perform LES of turbulent channel flow for Reynolds numbers Reτ based on the friction velocity uτ and the channel half-width δ in the range 2 × 103 to 2 × 107. Results, including SGS-extended longitudinal spectra, compare favourably with the direct numerical simulation (DNS) data of Hoyas & Jiménez (2006) at Reτ = 2003 and maintain an O(1) grid dependence on Reτ.


Author(s):  
Chuang Jin ◽  
Giovanni Coco ◽  
Rafael O. Tinoco ◽  
Pallav Ranjan ◽  
Jorge San Juan ◽  
...  

Fluids ◽  
2018 ◽  
Vol 3 (3) ◽  
pp. 65 ◽  
Author(s):  
Arne Heinrich ◽  
Guido Kuenne ◽  
Sebastian Ganter ◽  
Christian Hasse ◽  
Johannes Janicka

Combustion will play a major part in fulfilling the world’s energy demand in the next 20 years. Therefore, it is necessary to understand the fundamentals of the flame–wall interaction (FWI), which takes place in internal combustion engines or gas turbines. The FWI can increase heat losses, increase pollutant formations and lowers efficiencies. In this work, a Large Eddy Simulation combined with a tabulated chemistry approach is used to investigate the transient near wall behavior of a turbulent premixed stoichiometric methane flame. This sidewall quenching configuration is based on an experimental burner with non-homogeneous turbulence and an actively cooled wall. The burner was used in a previous study for validation purposes. The transient behavior of the movement of the flame tip is analyzed by categorizing it into three different scenarios: an upstream, a downstream and a jump-like upstream movement. The distributions of the wall heat flux, the quenching distance or the detachment of the maximum heat flux and the quenching point are strongly dependent on this movement. The highest heat fluxes appear mostly at the jump-like movement because the flame behaves locally like a head-on quenching flame.


Author(s):  
Lara Schembri Puglisevich ◽  
Gary Page

Unsteady Large Eddy Simulation (LES) is carried out for the flow around a bluff body equipped with an underbody rear diffuser in close proximity to the ground, representing an automotive diffuser. The goal is to demonstrate the ability of LES to model underbody vortical flow features at experimental Reynolds numbers (1.01 × 106 based on model height and incoming velocity). The scope of the time-dependent simulations is not to improve on Reynolds-Averaged Navier Stokes (RANS), but to give further insight into vortex formation and progression, allowing better understanding of the flow, hence allowing more control. Vortical flow structures in the diffuser region, along the sides and top surface of the bluff body are successfully modelled. Differences between instantaneous and time-averaged flow structures are presented and explained. Comparisons to pressure measurements from wind tunnel experiments on an identical bluff body model shows a good level of agreement.


2015 ◽  
Vol 137 (8) ◽  
Author(s):  
Yuchuan Wang ◽  
Lei Tan ◽  
Binbin Wang ◽  
Shuliang Cao ◽  
Baoshan Zhu

Large eddy simulation (LES) approach was used to investigate jumps of primary frequency of shear layer flow over a cavity. Comparisons between computational results and experimental data show that LES is an appropriate approach to accurately capturing the critical values of velocity and cavity length of a frequency jump, as well as characteristics of the separated shear layer. The drive force of the self-sustained oscillation of impinging shear layer is fluid injection and reinjection. Flow patterns in the shear layer and cavity before and after the frequency jump demonstrate that the frequency jump is associated with vortex–corner interaction. Before frequency jump, a mature vortex structure is observed in shear layer. The vortex is clipped by impinging corner at approximately half of its size, which induces strong vortex–corner interaction. After frequency jump, successive vortices almost escape from impinging corner without the generation of a mature vortex, thereby indicating weaker vortex–corner interaction. Two wave peaks are observed in the shear layer after the frequency jump because of: (1) vortex–corner interaction and (2) centrifugal instability in cavity. Pressure fluctuations inside the cavity are well regulated with respect to time. Peak values of correlation coefficients close to zero time lags indicate the existence of standing waves inside the cavity. Transitions from a linear to a nonlinear process occurs at the same position (i.e., x/H = 0.7) for both velocity and cavity length variations. Slopes of linear region are solely the function of cavity length, thereby showing increased steepness with increased cavity length.


Sign in / Sign up

Export Citation Format

Share Document