Large-eddy simulation and parameterization of buoyant plume dynamics in stratified flow

2016 ◽  
Vol 794 ◽  
pp. 798-833 ◽  
Author(s):  
Di Yang ◽  
Bicheng Chen ◽  
Scott A. Socolofsky ◽  
Marcelo Chamecki ◽  
Charles Meneveau

Characteristics of laboratory-scale bubble-driven buoyant plumes in a stably stratified quiescent fluid are studied using large-eddy simulation (LES). As a bubble plume entrains stratified ambient water, its net buoyancy decreases due to the increasing density difference between the entrained and ambient fluids. A large fraction of the entrained fluid eventually detrains and falls along an annular outer plume from a height of maximum rise (peel height) to a neutral buoyancy level (trap height), during which less buoyant scalars (e.g. small droplets) are trapped and dispersed horizontally, forming quasi-horizontal intrusion layers. The inner/outer double-plume structure and the peel/intrusion process are found to be more distinct for cases with small bubble rise velocity, while weak and unstable when the slip velocity is large. LES results are averaged to generate distributions of mean velocity and turbulent fluxes. These distributions provide data for assessing the performance of previously developed closures used in one-dimensional integral plume models. In particular, the various LES cases considered in this study yield consistent behaviour for the entrainment coefficients for various plume cases. Furthermore, a new continuous peeling model is derived based on the insights obtained from LES results. Comparing to previous peeling models, the new model behaves in a more self-consistent manner, and it is expected to provide more reliable performance when applied in integral plume models.

2010 ◽  
Vol 652 ◽  
pp. 75-103 ◽  
Author(s):  
B. J. DEVENISH ◽  
G. G. ROONEY ◽  
D. J. THOMSON

We consider large-eddy simulation (LES) of buoyant plumes in uniform and stably stratified environments. We show that in the former case the results agree well with the simple plume model of Morton, Taylor & Turner (Proc. R. Soc. Lond. A, vol. 234, 1956, p. 1). In particular, we calculate an entrainment constant which is consistent with laboratory and field measurements and find no significant difference between the radial spreading rates of vertical velocity and buoyancy. In a stably stratified environment, the LES plume shows better agreement with Morton et al. (1956) below the level at which the buoyancy first vanishes than above this level. Above the level of neutral buoyancy, the LES plume is characterized by an ascending core of negative buoyancy surrounded by a descending annulus of positive buoyancy. We compare the LES data with the model of Bloomfield & Kerr (J. Fluid Mech., vol. 424, 2000, p. 197), which explicitly accounts for these coherent motions. The model exhibits many qualitative aspects of the LES plume and quantitative agreement can be improved by adjusting the downward volume flux relative to the upward volume flux in a manner consistent with the LES plume. This simple adjustment, along with revised values of the entrainment constants, represents the combined effects of an overturning region at the top of the plume (where a fluid element reverses direction), ‘plume-top’ entrainment (whereby the plume entrains ambient fluid above the plume) as well as lateral entrainment and detrainment processes (both external and internal) occurring above the top of the model plume.


Author(s):  
N Kharoua ◽  
L Khezzar

Large eddy simulation of turbulent flow around smooth and rough hemispherical domes was conducted. The roughness of the rough dome was generated by a special approach using quadrilateral solid blocks placed alternately on the dome surface. It was shown that this approach is capable of generating the roughness effect with a relative success. The subgrid-scale model based on the transport of the subgrid turbulent kinetic energy was used to account for the small scales effect not resolved by large eddy simulation. The turbulent flow was simulated at a subcritical Reynolds number based on the approach free stream velocity, air properties, and dome diameter of 1.4 × 105. Profiles of mean pressure coefficient, mean velocity, and its root mean square were predicted with good accuracy. The comparison between the two domes showed different flow behavior around them. A flattened horseshoe vortex was observed to develop around the rough dome at larger distance compared with the smooth dome. The separation phenomenon occurs before the apex of the rough dome while for the smooth dome it is shifted forward. The turbulence-affected region in the wake was larger for the rough dome.


2021 ◽  
Vol 932 ◽  
Author(s):  
Changping Yu ◽  
Zelong Yuan ◽  
Han Qi ◽  
Jianchun Wang ◽  
Xinliang Li ◽  
...  

Kinetic energy flux (KEF) is an important physical quantity that characterizes cascades of kinetic energy in turbulent flows. In large-eddy simulation (LES), it is crucial for the subgrid-scale (SGS) model to accurately predict the KEF in turbulence. In this paper, we propose a new eddy-viscosity SGS model constrained by the properly modelled KEF for LES of compressible wall-bounded turbulence. The new methodology has the advantages of both accurate prediction of the KEF and strong numerical stability in LES. We can obtain an approximate KEF by the tensor-diffusivity model, which has a high correlation with the real value. Then, using the artificial neural network method, the local ratios between the real KEF and the approximate KEF are accurately modelled. Consequently, the SGS model can be improved by the product of that ratio and the approximate KEF. In LES of compressible turbulent channel flow, the new model can accurately predict mean velocity profile, turbulence intensities, Reynolds stress, temperature–velocity correlation, etc. Additionally, for the case of a compressible flat-plate boundary layer, the new model can accurately predict some key quantities, including the onset of transitions and transition peaks, the skin-friction coefficient, the mean velocity in the turbulence region, etc., and it can also predict the energy backscatters in turbulence. Furthermore, the proposed model also shows more advantages for coarser grids.


Author(s):  
Mehmet Kırtas¸ ◽  
Nayan Patel ◽  
Vaidyanathan Sankaran ◽  
Suresh Menon

Large-eddy simulation (LES) of a lean-direct injection (LDI) combustor is reported in this paper. The full combustor and all the six swirl vanes are resolved and both cold and reacting flow simulations are performed. Cold flow predictions with LES indicate the presence of a broad central recirculation zone due to vortex breakdown phenomenon near the dump plane and two corner recirculation zones at the top and bottom corner of the combustor. These predicted features compare well with the experimental non-reacting data. Reacting case simulated a liquid Jet-A fuel spray using a Lagrangian approach. A three-step kinetics model that included CO and NO is used for the chemistry. Comparison of mean velocity field predicted in the reacting LES with experiments shows reasonable agreement. Comparison with the non-reacting case shows that the centerline recirculation bubble is shorter but more intense in the reacting case.


1998 ◽  
Vol 120 (2) ◽  
pp. 243-252 ◽  
Author(s):  
S. A. Jordan ◽  
S. A. Ragab

The formation and the downstream transport of the Strouhal vortices in the near wake of a circular cylinder are investigated using the large-eddy simulation (LES) method. The governing equations are formulated in curvilinear coordinates to accommodate a nonorthogonal grid with formal development of a dynamic model to account for the subgrid turbulent scales. Results were produced with and without use of the model. The focus of the investigation is at a subcritical Reynolds number of 5600. Using the dynamic model, the LES results compared best to the published experimental data in terms of both the global and local wake characteristics such as the drag and base pressure coefficients, shedding and detection frequencies, peak vorticity, and the downstream mean velocity-defect and Reynolds stresses. The results further showed streamwise filaments that connect subsequent Strouhal vortices. Qualitatively, the time-averaged Reynolds stresses of the formation region revealed similar symmetric characteristics over the range 525 ≤ Re ≤ 140,000.


2009 ◽  
Vol 643 ◽  
pp. 233-266 ◽  
Author(s):  
BISHAKHDATTA GAYEN ◽  
SUTANU SARKAR ◽  
JOHN R. TAYLOR

A numerical study based on large eddy simulation is performed to investigate a bottom boundary layer under an oscillating tidal current. The focus is on the boundary layer response to an external stratification. The thermal field shows a mixed layer that is separated from the external stratified fluid by a thermocline. The mixed layer grows slowly in time with an oscillatory modulation by the tidal flow. Stratification strongly affects the mean velocity profiles, boundary layer thickness and turbulence levels in the outer region although the effect on the near-bottom unstratified fluid is relatively mild. The turbulence is asymmetric between the accelerating and decelerating stages. The asymmetry is more pronounced with increasing stratification. There is an overshoot of the mean velocity in the outer layer; this jet is linked to the phase asymmetry of the Reynolds shear stress gradient by using the simulation data to examine the mean momentum equation. Depending on the height above the bottom, there is a lag of the maximum turbulent kinetic energy, dissipation and production with respect to the peak external velocity and the value of the lag is found to be influenced by the stratification. Flow instabilities and turbulence in the bottom boundary layer excite internal gravity waves that propagate away into the ambient. Unlike the steady case, the phase lines of the internal waves change direction during the tidal cycle and also from near to far field. The frequency spectrum of the propagating wave field is analysed and found to span a narrow band of frequencies clustered around 45°.


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.


Sign in / Sign up

Export Citation Format

Share Document