scholarly journals Explicit filtering to obtain grid-spacing-independent and discretization-order-independent large-eddy simulation of compressible single-phase flow

2012 ◽  
Vol 697 ◽  
pp. 399-435 ◽  
Author(s):  
Senthilkumaran Radhakrishnan ◽  
Josette Bellan

AbstractIn large-eddy simulation (LES), it is often assumed that the filter width is equal to the grid spacing. Predictions from such LES are grid-spacing dependent since any subgrid-scale (SGS) model used in the LES equations is dependent on the resolved flow field which itself varies with grid spacing. Moreover, numerical errors affect the flow field, especially the smallest resolved scales. Thus, predictions using this approach are affected by both modelling and numerical choices. However, grid-spacing-independent LES predictions unaffected by numerical choices are necessary to validate LES models through comparison with a trusted template. First, such a template is created here through direct numerical simulation (DNS). Then, simulations are conducted using the conventional LES equations and also LES equations which are here reformulated so that the small-scale-producing nonlinear terms in these equations are explicitly filtered (EF) to remove scales smaller than a fixed filter width; this formulation is called EFLES. First, LES is conducted with four SGS models, then EFLES is performed with two of the SGS models used in LES; the results from all these simulations are compared to those from DNS and from the filtered DNS (FDNS). The conventional LES solution is both grid-spacing and spatial discretization-order dependent, thus showing that both of these numerical aspects affect the flow prediction. The solution from the EFLES equations is grid independent for a high-order spatial discretization on all meshes tested. However, low-order discretizations require a finer mesh to reach grid independence. With an eighth-order discretization, a filter-width to grid-spacing ratio of two is sufficient to reach grid independence, while a filter-width to grid-spacing ratio of four is needed to reach grid independence when a fourth- or a sixth-order discretization is employed. On a grid fine enough to be utilized in a DNS, the EFLES solution exhibits grid independence and does not converge to the DNS solution. The velocity-fluctuation spectra of EFLES follow those of FDNS independent of the grid spacing used, in concert with the original concept of LES. The reasons for the different predictions of conventional LES or EFLES according to the SGS model used, and the different characteristics of the EFLES predictions compared to those from conventional LES are analysed.

2021 ◽  
Vol 11 (15) ◽  
pp. 7167
Author(s):  
Liang Xu ◽  
Xu Zhao ◽  
Lei Xi ◽  
Yonghao Ma ◽  
Jianmin Gao ◽  
...  

Swirling impinging jet (SIJ) is considered as an effective means to achieve uniform cooling at high heat transfer rates, and the complex flow structure and its mechanism of enhancing heat transfer have attracted much attention in recent years. The large eddy simulation (LES) technique is employed to analyze the flow fields of swirling and non-swirling impinging jet emanating from a hole with four spiral and straight grooves, respectively, at a relatively high Reynolds number (Re) of 16,000 and a small jet spacing of H/D = 2 on a concave surface with uniform heat flux. Firstly, this work analyzes two different sub-grid stress models, and LES with the wall-adapting local eddy-viscosity model (WALEM) is established for accurately predicting flow and heat transfer performance of SIJ on a flat surface. The complex flow field structures, spectral characteristics, time-averaged flow characteristics and heat transfer on the target surface for the swirling and non-swirling impinging jets are compared in detail using the established method. The results show that small-scale recirculation vortices near the wall change the nearby flow into an unstable microwave state, resulting in small-scale fluctuation of the local Nusselt number (Nu) of the wall. There is a stable recirculation vortex at the stagnation point of the target surface, and the axial and radial fluctuating speeds are consistent with the fluctuating wall temperature. With the increase in the radial radius away from the stagnation point, the main frequency of the fluctuation of wall temperature coincides with the main frequency of the fluctuation of radial fluctuating velocity at x/D = 0.5. Compared with 0° straight hole, 45° spiral hole has a larger fluctuating speed because of speed deflection, resulting in a larger turbulence intensity and a stronger air transport capacity. The heat transfer intensity of the 45° spiral hole on the target surface is slightly improved within 5–10%.


Author(s):  
Yigang Luan ◽  
Lianfeng Yang ◽  
Bo Wan ◽  
Tao Sun

Gas turbine engines have been widely used in modern industry especially in the aviation, marine and energy fields. The efficiency of gas turbines directly affects the economy and emissions. It’s acknowledged that the higher turbine inlet temperatures contribute to the overall gas turbine engine efficiency. Since the components are subject to the heat load, the internal cooling technology of turbine blades is of vital importance to ensure the safe and normal operation. This paper is focused on exploring the flow and heat transfer mechanism in matrix cooling channels. In order to analyze the internal flow field characteristics of this cooling configuration at a Reynolds number of 30000 accurately, large eddy simulation method is carried out. Methods of vortex identification and field synergy are employed to study its flow field. Cross-sectional views of velocity in three subchannels at different positions have been presented. The results show that the airflow is strongly disturbed by the bending part. It’s concluded that due to the bending structure, the airflow becomes complex and disordered. When the airflow goes from the inlet to the turning, some small-sized and discontinuous vortices are formed. Behind the bending structure, the size of the vortices becomes big and the vortices fill the subchannels. Because of the structure of latticework, the airflow is affected by each other. Airflow in one subchannel can exert a shear force on another airflow in the opposite subchannel. It’s the force whose direction is the same as the vortex that enhances the longitudinal vortices. And the longitudinal vortices contribute to the energy exchange of the internal airflow and the heat transfer between airflow and walls. Besides, a comparison of the CFD results and the experimental data is made to prove that the numerical simulation methods are reasonable and acceptable.


2019 ◽  
Vol 213 ◽  
pp. 02076
Author(s):  
Jan Sip ◽  
Frantisek Lizal ◽  
Jakub Elcner ◽  
Jan Pokorny ◽  
Miroslav Jicha

The velocity field in the area behind the automotive vent was measured by hot-wire anenemometry in detail and intensity of turbulence was calculated. Numerical simulation of the same flow field was performed using Computational fluid dynamics in commecial software STAR-CCM+. Several turbulence models were tested and compared with Large Eddy Simulation. The influence of turbulence model on the results of air flow from the vent was investigated. The comparison of simulations and experimental results showed that most precise prediction of flow field was provided by Spalart-Allmaras model. Large eddy simulation did not provide results in quality that would compensate for the increased computing cost.


2019 ◽  
Vol 885 ◽  
Author(s):  
A. E. Tejada-Martínez ◽  
A. Hafsi ◽  
C. Akan ◽  
M. Juha ◽  
F. Veron


2015 ◽  
Vol 2015 ◽  
pp. 1-7
Author(s):  
Yan Xu ◽  
Zunce Wang ◽  
Lin Ke ◽  
Sen Li ◽  
Jinglong Zhang

Reynolds Stress Model and Large Eddy Simulation are used to respectively perform numerical simulation for the flow field of a hydrocyclone. The three-dimensional hexahedral computational grids were generated. Turbulence intensity, vorticity, and the velocity distribution of different cross sections were gained. The velocity simulation results were compared with the LDV test results, and the results indicated that Large Eddy Simulation was more close to LDV experimental data. Large Eddy Simulation was a relatively appropriate method for simulation of flow field within a hydrocyclone.


2015 ◽  
Vol 40 (7) ◽  
pp. 3098-3109 ◽  
Author(s):  
M.A. Abdel-Raheem ◽  
S.S. Ibrahim ◽  
W. Malalasekera ◽  
A.R. Masri

2011 ◽  
Vol 679 ◽  
pp. 288-314 ◽  
Author(s):  
W. ANDERSON ◽  
C. MENEVEAU

Many flows especially in geophysics involve turbulent boundary layers forming over rough surfaces with multiscale height distribution. Such surfaces pose special challenges for large-eddy simulation (LES) when the filter scale is such that only part of the roughness elements of the surface can be resolved. Here we consider LES of flows over rough surfaces with power-law height spectra Eh(k) ~ kβs (−3 ≤ βs < −1), as often encountered in natural terrains. The surface is decomposed into resolved and subgrid-scale height contributions. The effects of the unresolved small-scale height fluctuations are modelled using a local equilibrium wall model (log-law or Monin–Obukhov similarity), but the required hydrodynamic roughness length must be specified. It is expressed as the product of the subgrid-scale root-mean-square of the height distribution and an unknown dimensionless quantity, α, the roughness parameter. Instead of specifying this parameter in an ad hoc empirical fashion, a dynamic methodology is proposed based on test-filtering the surface forces and requiring that the total drag force be independent of filter scale or resolution. This dynamic surface roughness (DSR) model is inspired by the Germano identity traditionally used to determine model parameters for closing subgrid-scale stresses in the bulk of a turbulent flow. A series of LES of fully developed flow over rough surfaces are performed, with surfaces built using random-phase Fourier modes with prescribed power-law spectra. Results show that the DSR model yields well-defined, rapidly converging, values of α. Effects of spatial resolution and spectral slopes are investigated. The accuracy of the DSR model is tested by showing that predicted mean velocity profiles are approximately independent of resolution for the dynamically computed values of α, whereas resolution-dependent results are obtained when using other, incorrect, α values. Also, strong dependence of α on βs is found, where α ranges from α ~ 0.1 for βs = −1.2 to α ~ 10−5 for βs = −3.


2015 ◽  
Vol 93 (10) ◽  
pp. 1124-1130 ◽  
Author(s):  
T. Wang ◽  
P. Li ◽  
J.S. Bai ◽  
G. Tao ◽  
B. Wang ◽  
...  

The subgrid-scale (SGS) terms of turbulence transport are modelled by the stretched-vortex SGS stress model, and a large-eddy simulation code multi-viscous fluid and turbulence (MVFT) is developed to investigate the MVFT problems. Then one AWE shock tube experiment of interface instability is simulated numerically by MVFT code, which reproduces the development process of the interface. The obtained numerical images of interface evolution and wave structures in flow field are consistent with the experimental results. The evolution of perturbed interface and propagation of shock waves in flow field and their interactions are analyzed in detail. The statistics features of turbulence mixing in the form of finer quantities, such as the turbulent kinetic energy, enstrophy, density variance, and turbulent mass flux are investigated, which also proves that the SGS model has a key role in large-eddy simulation. The turbulent kinetic energy and enstrophy decay with time as a power law.


Sign in / Sign up

Export Citation Format

Share Document