Subfilter-scale enrichment of planetary boundary layer large eddy simulation using discrete Fourier–Gabor modes

2017 ◽  
Vol 819 ◽  
pp. 494-539 ◽  
Author(s):  
Aditya S. Ghate ◽  
Sanjiva K. Lele

A new multiscale simulation methodology is introduced to facilitate computationally efficient simulations of high Reynolds number turbulence seen in wall-bounded flows. The scale splitting methodology uses traditional large eddy simulation (LES) with a wall model to simulate the larger scales which are subsequently enriched using a space–time compatible kinematic simulation. Computational feasibility and robustness of the methodology are investigated using two idealized problems that emulate turbulence within the planetary boundary layer (PBL), and a finite Reynolds number channel flow problem which serves to validate the methodology against direct numerical simulation. The space–time correlations and spectra generated using enriched LES show excellent agreement with LES conducted at high resolution for all three problems; thereby demonstrating the potential of this approach for high resolution PBL simulations with a drastic reduction in the computational costs when compared to the conventional approach.

2019 ◽  
Vol 12 (6) ◽  
pp. 2523-2538 ◽  
Author(s):  
Sadiq Huq ◽  
Frederik De Roo ◽  
Siegfried Raasch ◽  
Matthias Mauder

Abstract. Large-eddy simulation (LES) has become a well-established tool in the atmospheric boundary layer research community to study turbulence. It allows three-dimensional realizations of the turbulent fields, which large-scale models and most experimental studies cannot yield. To resolve the largest eddies in the mixed layer, a moderate grid resolution in the range of 10 to 100 m is often sufficient, and these simulations can be run on a computing cluster with a few hundred processors or even on a workstation for simple configurations. The desired resolution is usually limited by the computational resources. However, to compare with tower measurements of turbulence and exchange fluxes in the surface layer, a much higher resolution is required. In spite of the growth in computational power, a high-resolution LES of the surface layer is often not feasible: to fully resolve the energy-containing eddies near the surface, a grid spacing of O(1 m) is required. One way to tackle this problem is to employ a vertical grid nesting technique, in which the surface is simulated at the necessary fine grid resolution, and it is coupled with a standard, coarse, LES that resolves the turbulence in the whole boundary layer. We modified the LES model PALM (Parallelized Large-eddy simulation Model) and implemented a two-way nesting technique, with coupling in both directions between the coarse and the fine grid. The coupling algorithm has to ensure correct boundary conditions for the fine grid. Our nesting algorithm is realized by modifying the standard third-order Runge–Kutta time stepping to allow communication of data between the two grids. The two grids are concurrently advanced in time while ensuring that the sum of resolved and sub-grid-scale kinetic energy is conserved. We design a validation test and show that the temporally averaged profiles from the fine grid agree well compared to the reference simulation with high resolution in the entire domain. The overall performance and scalability of the nesting algorithm is found to be satisfactory. Our nesting results in more than 80 % savings in computational power for 5 times higher resolution in each direction in the surface layer.


1997 ◽  
Vol 336 ◽  
pp. 151-182 ◽  
Author(s):  
BRANKO KOSOVIĆ

It has been recognized that the subgrid-scale (SGS) parameterization represents a critical component of a successful large-eddy simulation (LES). Commonly used linear SGS models produce erroneous mean velocity profiles in LES of high-Reynolds-number boundary layer flows. Although recently proposed approaches to solving this problem have resulted in significant improvements, questions about the true nature of the SGS problem in shear-driven high-Reynolds-number flows remain open.We argue that the SGS models must capture inertial transfer effects including backscatter of energy as well as its redistribution among the normal SGS stress components. These effects are the consequence of nonlinear interactions and anisotropy. In our modelling procedure we adopt a phenomenological approach whereby the SGS stresses are related to the resolved velocity gradients. We show that since the SGS stress tensor is not frame indifferent a more general nonlinear model can be applied to the SGS parameterization. We develop a nonlinear SGS model capable of reproducing the effects of SGS anisotropy characteristic for shear-driven boundary layers. The results obtained using the nonlinear model for the LES of a neutral shear-driven atmospheric boundary layer show a significant improvement in prediction of the non-dimensional shear and low-order statistics compared to the linear Smagorinsky-type models. These results also demonstrate a profound effect of the SGS model on the flow structures.


Author(s):  
Susumu Teramoto ◽  
Takuya Ouchi ◽  
Hiroki Sanada ◽  
Koji Okamoto

Fully resolved large eddy simulation (LES) is applied to two simple geometry flowfields with well-defined boundary conditions. The LES results are compared with simulations based on a Reynolds-averaged Navier-Stokes (RANS) model with turbulence, and pros and cons of using high-resolution LES for turbomachinery flows are discussed. One flow is a linear compressor cascade flow composed of the tip section of GE rotor B at Rec = 4 × 105 with a clearance, and the other is a Mach 1.76 supersonic turbulent boundary layer at Reδ = 5000 that laminerizes through a 12-degree expansion corner. The grids are prepared fine enough to resolve the turbulent boundary layer through a grid sensitivity study. The liner cascade result shows that all the turbulent shear layers and boundary layers including those in the small tip clearance are well resolved with 800 million grid points. The Reynolds stress derived from the LES results are compared directly with those predicted from the Spalart-Allmaras one-equation RANS turbulence model. The two results agreed qualitatively well except for the shear layer surrounding the tip leakage vortex, demonstrating that the RANS model performs well at least for flowfields near the design condition. From the simulation of the turbulent boundary layer experiencing sudden expansion, noticeable decreases of both Reynolds stress and local friction coefficient were observed, showing that the turbulent boundary layer has relaminarized through the sudden expansion. The boundary layer downstream of the expansion exhibits a nonequilibrium condition and was different from the laminar boundary layer.


1996 ◽  
Vol 80 (1-2) ◽  
pp. 167-202 ◽  
Author(s):  
Peter P. Sullivan ◽  
James C. McWilliams ◽  
Chin-Hoh Moeng

Sign in / Sign up

Export Citation Format

Share Document