Toward advanced subgrid models for Lattice-Boltzmann-based Large-eddy simulation: Theoretical formulations

This paper addresses the issue of developing advanced subgrid model for large-eddy simulations (LES) of turbulent flows based on Lattice Boltzmann methods (LBM). Most of already existing subgrid closures used in LES-LBM are straightforward extensions of the most crude model developed within the Navier–Stokes equations, namely the Smagorinsky eddy-viscosity model. In a first part, it is shown how to obtain an improved eddy-viscosity subgrid model for LBM. The original implementation of the Inertial-Range Consistent Smagorinsky model proposed by Dong and Sagaut for the D3Q19 scheme is used as an illustration. In a second step, an original extension of the Approximate Deconvolution Method proposed by Adams and Stolz for Navier–Stokes simulation is proposed. This new LBM-LES approach does not rely on the eddy-viscosity concept and is written directly within the LBM framework. It is shown that it can be implemented thanks to a trivial modification of the existing LBM solvers for Direct Numerical Simulation.     

Numerical studies of finite element variational multiscale methods for turbulent flow simulations

Different realizations of variational multiscale (VMS) methods within the framework of finite element methods are studied in turbulent channel flow simulations. One class of VMS methods uses bubble functions to model resolved small scales whereas the other class contains the definition of the resolved small scales by an explicit projection in its set of equations. All methods are employed with eddy viscosity models of Smagorinsky type. The simulations are performed on grids for which a Direct Numerical Simulation blows up in finite time.     

Direct numerical simulations of transverse and spanwise magnetic field effects on turbulent flow in a 2:1 aspect ratio rectangular duct

Magnetic fields are used extensively to direct liquid metal flows in material processing. Continuous casting of steel uses different configurations of magnetic fields to optimize turbulent flows in rectangular cross-sections to minimize defects in the solidified steel product. Realizing the importance of a magnetic field on turbulent flows in rectangular cross-sections, the present work is aimed at understanding the effect of a magnetic field on the turbulent metal flow at a nominal bulk Reynolds number of ∼5300 (based upon full duct height) (Re_τ = 170, based upon half duct height) and Hartmann numbers (based upon half duct height) of 0, 6.0 and 8.25 in a 2:1 aspect ratio rectangular duct. Direct numerical simulations in a non-MHD 2:1 aspect ratio duct followed by simulations with transverse and span-wise magnetic fields have been performed with 224 × 120 × 512 cells (∼13.7 million cells). The fractional step method with second order space and time discretization schemes has been used to solve the coupled Navier-Stokes-MHD equations. Instantaneous and time-averaged natures of the flow have been examined through distribution of velocities, various turbulence parameters and budget terms. Spanwise (horizontal) magnetic field reorganizes and suppresses secondary flows more strongly. Turbulence suppression and velocity flattening effects are stronger with transverse (vertical) magnetic field.     

A parallel, finite-volume algorithm for large-eddy simulation of turbulent flows

A parallel, finite-volume algorithm has been developed for large-eddy simulation (LES) of compressible turbulent flows. This algorithm includes piecewise linear least-square reconstruction, trilinear finite-element interpolation, Roe flux-difference splitting (FDS), and second-order MacCormack time marching. A systematic and consistent means of evaluating the surface and volume integrals of the control volume is described. Parallel implementation is done using the message-passing programming model. To validate the numerical method for turbulence simulation, LES of fully developed turbulent flow in a square duct is performed for a Reynolds number of 320 based on the average friction velocity and the hydraulic diameter of the duct. Direct numerical simulation (DNS) results are available for this test case, and the accuracy of this algorithm for turbulence simulations can be ascertained by comparing the LES solutions with the DNS results. For the first time, a finite volume method with Roe FDS was used for LES of turbulent flow in a square duct, and the effects of grid resolution, upwind numerical dissipation, and subgrid-scale dissipation on the accuracy of the LES are examined. Comparison with DNS results shows that the standard Roe FDS adversely affects the accuracy of the turbulence simulation. For accurate turbulence simulations, only 3–5% of the standard Roe FDS dissipation is needed.     

Study of jet precession, recirculation and vortex breakdown in turbulent swirling jets using LES

Large eddy simulations (LES) are used to investigate turbulent isothermal swirling flows with a strong emphasis on vortex breakdown, recirculation and instability behaviour. The Sydney swirl burner configuration is used for all simulated test cases from low to high swirl and Reynolds numbers. The governing equations for continuity and momentum are solved on a structured Cartesian grid, and a Smagorinsky eddy viscosity model with the localised dynamic procedure is used as the sub-grid scale turbulence model. The LES successfully predicts both the upstream first recirculation zone generated by the bluff body and the downstream vortex breakdown bubble. The frequency spectrum indicates the presence of low frequency oscillations and the existence of a central jet precession as observed in experiments. The LES calculations well captured the distinct precession frequencies. The results also highlight the precession mode of instability in the center jet and the oscillations of the central jet precession, which forms a precessing vortex core. The study further highlights the predictive capabilities of LES on unsteady oscillations of turbulent swirling flow fields and provides a good framework for complex instability investigations.     

Non-body-fitted Cartesian-mesh simulation of highly turbulent flows using multi-relaxation-time lattice Boltzmann method

This paper presents a lattice Boltzmann method (LBM) based study aimed at numerical simulation of highly turbulent and largely inclined flow around obstacles of curved geometry using non-body-fitted Cartesian meshes. The approach features (1) combining the interpolated bounce-back scheme with the LBM of multi-relaxation-time (MRT) type to enable the use of simple Cartesian mesh for the flow cases even with complex geometries; and (2) incorporating the Spalart–Allmaras (SA) turbulence model into LBM in order to represent the turbulent flow effect. The numerical experiments are performed corresponding to flows around an NACA0012 airfoil at Re=5×10⁵ and around a flat plate at Re=2×10⁴, respectively. The agreement between all simulation results obtained from this study and the data provided by other literature demonstrates the reliability of the enhanced LBM proposed in this paper for simulating, simply on Cartesian meshes, complex flows that may involve bodies of curved boundary, high Reynolds number, and large angle of attack.     

Depth-averaged Lattice Boltzmann and Finite Element methods for single-phase flows in fractures with obstacles

We use Lattice Boltzmann Method (LBM) MRT and Cumulant schemes to study the performance and accuracy of single-phase flow modeling for propped fractures. The simulations are run using both the two- and three-dimensional Stokes equations, and a 2.5D Stokes–Brinkman approximate model. The LBM results are validated against Finite Element Method (FEM) simulations and an analytical solution to the Stokes–Brinkman flow around an isolated circular obstacle. Both LBM and FEM 2.5D Stokes–Brinkman models are able to reproduce the analytical solution for an isolated circular obstacle. In the case of 2D Stokes and 2.5D Stokes–Brinkman models, the differences between the extrapolated fracture permeabilities obtained with LBM and FEM simulations for fractures with multiple obstacles are below 1%. The differences between the fracture permeabilities computed using 3D Stokes LBM and FEM simulations are below 8% . The differences between the 3D Stokes and 2.5 Stokes–Brinkman results are less than 7% for FEM study, and 8% for the LBM case. The velocity perturbations that are introduced around the obstacles are not fully captured by the parabolic velocity profile inherent to the 2.5D Stokes–Brinkman model.     

Free surface flow simulations on GPGPUs using the LBM

In this paper, we present the implementation of a volume-of-fluid-(VOF)-based algorithm for the simulation of free-surface flow problems on general purpose graphical processing units (GPGPUs). For the solution of the flow field and the additional advection equation for the VOF fill level, the lattice Boltzmann method on the basis of an MRT collision operator is used. A Smagorinsky LES model serves to capture the small-scale turbulent structures of the flow. We show that despite the additional non-local operations near the phase interface, we end up with an algorithm with good overall performance, which is suitable for the simulation of demanding real-world engineering applications.     

Flow simulation in Electro-Static-Precipitator (ESP) ducts with turning vanes

Flow simulation in inlet ducts along with several turning vanes used in electrostatic precipitator (ESP) are analysed to understand the flow pattern at its exit locations. The geometry of inlet duct has been extracted from the Plant Design Manufacturing System (PDMS) and refined with turning vanes placed at several locations. The domain of duct geometry around turning vanes are decomposed with several volumes and filled with hexahedral elements with the help of stat-of-art mesh generator - Hypermesh. The resulting computational grid has been used in TASCflow solver to predict its flow pattern in the duct. Simulation for the specified conditions predicts uneven flow distribution in the ESP inlet duct. Due to large flow recirculation and turbulent losses in the duct, non-uniform averaged mass flow rates are noticed at duct exit locations. Simulation results suggest that the improvement of flow distribution in the duct through optimization can be tried by placing more turning/splitter vanes in the inlet duct. In order to ensure that the results obtained from TASCflow are meaningful and in right direction, in the absence of measurement data, simulation was benchmarked with other industry standard commercial flow solvers. The observations made from these popular solvers confirm the findings obtained using the TASCflow solver. The analysis with multiple solvers indicates that Fluent provides quick results, while better visualization can be made using CFX solver. The Star-CD solver, which captures the turbulent losses accurately takes more time for convergence provides alternatives.     

Numerical prediction of sound generated from flows with a low Mach number

Numerical computations of sound generated from flows with a low Mach number are presented based on Lighthill’s acoustic analogy with an assumption that sound does not alter the flow field from which it is generated. The source fluctuations of the flow field are computed by a large-eddy simulation (LES) with Dynamic Smagorinsky Model (DSM) and they are fed to the following acoustical computation as input data. An explicit/implicit finite element method with second order accuracy both in time and space is used for flow field discretization. The method is applied to the prediction of sound in three different classes of problems: far-field sound generated from flow around a bluff body, sound resulting from blade-stator interaction of turbomachinery and sound due to a turbulent boundary layer on an aerofoil. The computed frequency spectra of the sound show a fairly good agreement with the measured spectra for all the cases.     

Application of Compact Schemes to Large Eddy Simulation of Turbulent Jets

We present 3-D large eddy simulation (LES) results for a turbulent Mach 0.9 isothermal round jet at a Reynolds number of 100,000 (based on jet nozzle exit conditions and nozzle diameter). Our LES code is part of a Computational Aeroacoustics (CAA) methodology that couples surface integral acoustics techniques such as Kirchhoff's method and the Ffowcs Williams– Hawkings method with LES for the far field noise estimation of turbulent jets. The LES code employs high-order accurate compact differencing together with implicit spatial filtering and state-of-the-art non-reflecting boundary conditions. A localized dynamic Smagorinsky subgrid-scale (SGS) model is used for representing the effects of the unresolved scales on the resolved scales. A computational grid consisting of 12 million points was used in the present simulation. Mean flow results obtained in our simulation are found to be in very good agreement with the available experimental data of jets at similar flow conditions. Furthermore, the near field data provided by the LES is coupled with the Ffowcs Williams–Hawkings method to compute the far field noise. Far field aeroacoustics results are also presented and comparisons are made with experimental measurements of jets at similar flow conditions. The aeroacoustics results are encouraging and suggest further investigation of the effects of inflow conditions on the jet acoustic field.     

A parallel lattice Boltzmann method for large eddy simulation on multiple GPUs

To improve the simulation efficiency of turbulent fluid flows at high Reynolds numbers with large eddy dynamics, a CUDA-based simulation solution of lattice Boltzmann method for large eddy simulation (LES) using multiple graphics processing units (GPUs) is proposed. Our solution adopts the “collision after propagation” lattice evolution way and puts the misaligned propagation phase at global memory read process. The latest GPU platform allows a single CPU thread to control up to four GPUs that run in parallel. In order to make use of multiple GPUs, the whole working set is evenly partitioned into sub-domains. We implement Smagorinsky model and Vreman model respectively to verify our multi-GPU solution. These two LES models have different relaxation time calculation behavior and lead to different CUDA implementation characteristics. The implementation based on Smagorinsky model achieves 190 times speedup over the sequential implementation on CPU, while the implementation based on Vreman model archives more than 90 times speedup. The experimental results show that the parallel performance of our multi-GPU solution scales very well on multiple GPUs. Therefore large-scale (up to 10,240 $\times $ 10,240 lattices) LES–LBM simulation becomes possible at a low cost, even using double-precision floating point calculation.     

Large eddy simulation of free-shear flows using an algebraic model

Large eddy simulation of the free-shear flows has been performed using an algebraic subgrid stress model. The test cases considered are the temporal planar jet and mixing layer. The results are compared with that of the Smagorinsky, dynamic Smagorinsky, mixed models and the DNS data. For the mean flow quantities none of the models perform exceptionally better over the others. However, both the algebraic and mixed models predict better statistics for the turbulent quantities than the linear models. The results also demonstrate the capability of the algebraic model in depicting backscatter of energy explicitly, which is consistent in both the cases considered.     

Globally conservative high-order filters for large-eddy simulation and computational aero-acoustics

In computational aero-acoustics, large-eddy simulations (LES) or direct numerical simulations (DNS) are often employed for flow computations in the source region. As part of the numerical implementation or required modeling, explicit spatial filters are frequently employed. For instance, in LES spatial filters are employed in the formulation of various subgrid-scale (SGS) models such as the dynamic model or the variational multi-scale (VMS) Smagorinsky model; both in LES or DNS, spatial high-pass filters are often used to remove undesired grid-to-grid oscillations. Though these type of spatial filters adhere to local accuracy requirements, in practice, they often destroy global conservation properties in the presence of non-periodic boundaries conditions. This leads to the incorrect prediction of the flow properties near hard boundaries, such as walls. In the current work, we present globally conservative high-order accurate filters, which combine traditional filters at the internal points with one-sided conservative filters near the wall boundary. We test these filters to remove grid-to-grid oscillations both in a channel-flow case and in 2D cavity flow. We find that the use of a non-conservative filter leads to erroneous predictions of the skin friction in channel flows up to 30%. In the cavity-flow simulations, the use of non-conservative filters to remove grid-to-grid oscillations leads to important shifts in the Strouhal number of the dominant mode, and a change of the flow pattern inside the cavity. In all cases, the use of conservative high-order filter formulations to remove grid-to-grid oscillations lead to very satisfactory results. Finally, in our channel-flow test case, we also illustrate the importance of using conservative filters for the formulation of the VMS Smagorinsky model.     

Highly interactive computational steering for coupled 3D flow problems utilizing multiple GPUs

Most computational fluid dynamics (CFD) simulations require massive computational power which is usually provided by traditional High Performance Computing (HPC) environments. Although interactivity of the simulation process is highly appreciated by scientists and engineers, due to limitations of typical HPC environments, present CFD simulations are usually executed non interactively. A recent trend is to harness the parallel computational power of graphics processing units (GPUs) for general purpose applications. As an alternative to traditional massively parallel computing, GPU computing has also gained popularity in the CFD community, especially for its application to the lattice Boltzmann method (LBM). For instance, Tölke and others presented very efficient implementations of the LBM for 2D as well as 3D space (Toelke J, in Comput Visual Sci. (2008); Toelke J and Krafczk M, in Int J Comput Fluid Dyn 22(7): 443–456 (2008)). In this work we motivate the use of GPU computing to facilitate interactive CFD simulations. In our approach, the simulation is executed on multiple GPUs instead of traditional HPC environments, which allows the integration of the complete simulation process into a single desktop application. To demonstrate the feasibility of our approach, we show a fully bidirectional fluid-structure-interaction for self induced membrane oscillations in a turbulent flow. The efficiency of the approach allows a 3D simulation close to realtime.     

On the interaction between dynamic model dissipation and numerical dissipation due to streamline upwind/Petrov–Galerkin stabilization

Here we investigate the roles of physical and numerical subgrid-scale modeling. The subgrid-scales are represented by a physical large-eddy simulation model, namely the popular dynamic Smagorinsky model (or simply dynamic model), as well as by a numerical model in the form of the well-known streamline upwind/Petrov–Galerkin stabilization for finite element discretizations of advection–diffusion systems. The latter is not a physical model, as its purpose is to provide sufficient algorithmic dissipation for a stable, consistent, and convergent numerical method. We study the interaction between the physical and numerical models by analyzing energy dissipation associated to the two. Based on this study, a modification to the dynamic model is proposed as a way to discount the numerical method’s algorithmic dissipation from the total subgrid-scale dissipation. The modified dynamic model is shown to be successful in simulations of turbulent channel flow.     

Improved coupling of time integration and hydrodynamic interaction in particle suspensions using the lattice Boltzmann and discrete element methods

This paper introduces improvements to the simulation of particle suspensions using the lattice Boltzmann method (LBM) and the discrete element method (DEM). First, the benefit of using a two-relaxation-time (TRT) collision operator, instead of the popular Bhatnagar–Gross–Krook (BGK) collision operator, is demonstrated. Second, a modified solid weighting function for the partially saturated method (PSM) for fluid–solid interaction is defined and tested. Results are presented for a range of flow configurations, including sphere packs, duct flows, and settling spheres, with good accuracy and convergence observed. Past research has shown that the drag, and consequently permeability, predictions of the LBM exhibit viscosity-dependence when used with certain boundary conditions such as bounce-back or interpolated bounce-back, and this is most pronounced when the BGK collision operator is employed. The improvements presented here result in a range of computational viscosities, and therefore relaxation parameters, within which drag and permeability predictions remain invariant. This allows for greater flexibility in using the relaxation parameter to adjust the LBM timestep, which can subsequently improve synchronisation with the time integration of the DEM. This has significant implications for the simulation of large-scale suspension phenomena, where the limits of computational hardware persistently constrain the resolution of the LBM lattice.     

Direct numerical simulation of turbulent channel flows using a stabilized finite element method

Direct numerical simulations (DNS) of incompressible turbulent channel flows at Re_τ = 180 and 395 (i.e., Reynolds number, based on the friction velocity and channel half-width) were performed using a stabilized finite element method (FEM). These simulations have been motivated by the fact that the use of stabilized finite element methods for DNS and LES is fairly recent and thus the question of how accurately these methods capture the wide range of scales in a turbulent flow remains open. To help address this question, we present converged results of turbulent channel flows under statistical equilibrium in terms of mean velocity, mean shear stresses, root mean square velocity fluctuations, autocorrelation coefficients, one-dimensional energy spectra and balances of the transport equation for turbulent kinetic energy. These results are consistent with previously published DNS results based on a pseudo-spectral method, thereby demonstrating the accuracy of the stabilized FEM for turbulence simulations.     

LES of turbulent square jet flow using an MRT lattice Boltzmann model

In this paper we consider the application of multiple-relaxation-time (MRT) lattice Boltzmann equation (LBE) for large-eddy simulation (LES) of turbulent flows. The implementation is discussed in the context of 19-velocity (D3Q19) MRT-LBE model in conjunction with the Smagorinsky subgrid closure model. The MRT-LBE-LES is then tested in the turbulent square jet flow case. We compare MRT-LBE-LES results with (a) single-relaxation-time (SRT) or BGK LBE results and (b) experimental data. Computed results include the distribution of centerline mean streamwise velocity, jet spread, and spanwise profiles of mean streamwise velocity in the near-field region. The phenomenon of axis switching is investigated. The advantages of MRT over SRT are demonstrated. Reasonable agreement between our numerical results and experimental data demonstrate that the MRT-LBE is a potentially viable tool for LES of turbulence.     

直方槽道中重力对颗粒相二次流的影响

采用Eulerian／Lagrangian方法模拟直方槽道中气粒两相流动过程。气相采用大涡模拟方法,直接求解大尺度涡运动,小尺度涡采用标准的Smagorinsky亚格子模式模拟,壁面采用幂次率应力模型代替无滑移边界条件。颗粒相采用轨道模型求解。大涡模拟预报的气相平均速度与DNS结果相吻合。结果表明,在直方槽道流向截面,气相存在二次流现象。受气相二次流的作用,颗粒相也存在类似于气相的二次流现象,并考察了重力对颗粒相二次流的影响。