LATTICE BOLTZMANN SIMULATION OF TWO-DIMENSIONAL WALL BOUNDED TURBULENT FLOW

2010 ◽  
Vol 21 (05) ◽  
pp. 669-680 ◽  
Author(s):  
GÁBOR HÁZI ◽  
GÁBOR TÓTH

This paper reports on a numerical study of two-dimensional decaying turbulence in a square domain with no-slip walls. The generation of strong small-scale vortices near the no-slip walls have been observed in the lattice Boltzmann simulations just like in earlier pseudospectral calculations. Due to these vortices the enstrophy is not a monotone decaying function of time. Considering a number of simulations and taking their ensemble average, we have found that the decay of enstrophy and that of the kinetic energy can be described well by power-laws. The exponents of these laws depend on the Reynolds number in a similar manner than was observed before in pseudospectral simulations. Considering the ensemble averaged 1D Fourier energy spectra calculated along the walls, we could not find a simple power-law, which fits well to the simulation data. These spectra change in time and reveal an exponent close to -3 in the intermediate and an exponent -5/3 at low wavenumbers. On the other hand, the two-dimensional energy spectra, which remain almost steady in the intermediate decay stage, show clear power-law behavior with exponent larger than -3 depending on the initial Reynolds number.

2011 ◽  
Vol 667 ◽  
pp. 463-473 ◽  
Author(s):  
ANDREAS VALLGREN

High-resolution simulations of forced two-dimensional turbulence reveal that the inverse cascade range is sensitive to an infrared Reynolds number, Reα = kf/kα, where kf is the forcing wavenumber and kα is a frictional wavenumber based on linear friction. In the limit of high Reα, the classic k−5/3 scaling is lost and we obtain steeper energy spectra. The sensitivity is traced to the formation of vortices in the inverse energy cascade range. Thus, it is hypothesized that the dual limit Reα → ∞ and Reν = kd/kf → ∞, where kd is the small-scale dissipation wavenumber, will lead to a steeper energy spectrum than k−5/3 in the inverse energy cascade range. It is also found that the inverse energy cascade is maintained by non-local triad interactions.


2004 ◽  
Vol 15 (03) ◽  
pp. 435-445 ◽  
Author(s):  
SANTOSH ANSUMALI ◽  
SHYAM SUNDER CHIKATAMARLA ◽  
CHRISTOS EMMANOUIL FROUZAKIS ◽  
KONSTANTINOS BOULOUCHOS

Minimal Boltzmann kinetic models, such as lattice Boltzmann, are often used as an alternative to the discretization of the Navier–Stokes equations for hydrodynamic simulations. Recently, it was argued that modeling sub-grid scale phenomena at the kinetic level might provide an efficient tool for large scale simulations. Indeed, a particular variant of this approach, known as the entropic lattice Boltzmann method (ELBM), has shown that an efficient coarse-grained simulation of decaying turbulence is possible using these approaches. The present work investigates the efficiency of the entropic lattice Boltzmann in describing flows of engineering interest. In order to do so, we have chosen the flow past a square cylinder, which is a simple model of such flows. We will show that ELBM can quantitatively capture the variation of vortex shedding frequency as a function of Reynolds number in the low as well as the high Reynolds number regime, without any need for explicit sub-grid scale modeling. This extends the previous studies for this set-up, where experimental behavior ranging from Re ~O(10) to Re ≤1000 was predicted by a single simulation algorithm.1–5


2014 ◽  
Vol 15 (1) ◽  
pp. 265-284 ◽  
Author(s):  
Qiuxiang Li ◽  
Ning Hong ◽  
Baochang Shi ◽  
Zhenhua Chai

AbstractIn this paper, the power-law fluid flows in a two-dimensional square cavity are investigated in detail with multi-relaxation-time lattice Boltzmann method (MRT-LBM). The influence of the Reynolds number (Re) and the power-law index (n) on the vortex strength, vortex position and velocity distribution are extensively studied. In our numerical simulations, Re is varied from 100 to 10000, and n is ranged from 0.25 to 1.75, covering both cases of shear-thinning and shear-thickening. Compared with the Newtonian fluid, numerical results show that the flow structure and number of vortex of power-law fluid are not only dependent on the Reynolds number, but also related to power-law index.


A numerical study on the transition from laminar to turbulent of two-dimensional fuel jet flames developed in a co-flowing air stream was made by adopting the flame surface model of infinite chemical reaction rate and unit Lewis number. The time dependent compressible Navier–Stokes equation was solved numerically with the equation for coupling function by using a finite difference method. The temperature-dependence of viscosity and diffusion coefficient were taken into account so as to study effects of increases of these coefficients on the transition. The numerical calculation was done for the case when methane is injected into a co-flowing air stream with variable injection Reynolds number up to 2500. When the Reynolds number was smaller than 1000 the flame, as well as the flow, remained laminar in the calculated domain. As the Reynolds number was increased above this value, a transition point appeared along the flame, downstream of which the flame and flow began to fluctuate. Two kinds of fluctuations were observed, a small scale fluctuation near the jet axis and a large scale fluctuation outside the flame surface, both of the same origin, due to the Kelvin–Helmholtz instability. The radial distributions of density and transport coefficients were found to play dominant roles in this instability, and hence in the transition mechanism. The decreased density in the flame accelerated the instability, while the increase in viscosity had a stabilizing effect. However, the most important effect was the increase in diffusion coefficient. The increase shifted the flame surface, where the large density decrease occurs, outside the shear layer of the jet and produced a thick viscous layer surrounding the jet which effectively suppressed the instability.


Author(s):  
Hiroshi Suzuki ◽  
Shinpei Maeda ◽  
Yoshiyuki Komoda

Two-dimensional numerical computations have been performed in order to investigate the development characteristics of flow and thermal field in a flow between parallel plates swept by a visco-elastic fluid. In the present study, the effect of the cavity number in the domain and of Reynolds number was focused on when the geometric parameters were set constant. From the results, it is found that the flow penetration into the cavities effectively causes the heat transfer augmentation in the cavities in any cavity region compared with that of water case. It is also found that the development of thermal field in cases of the present visco-elastic fluid is quicker compared with that of water cases. The present heat transfer augmentation technique using Barus effect of a visco-elastic fluid is effective in the range of low Reynolds number.


2013 ◽  
Vol 5 (3) ◽  
pp. 435-445
Author(s):  
M. S. I. Mallik ◽  
M. A. Uddin ◽  
M. A. Rahman

Direct numerical simulation (DNS) in two-dimensional homogeneous isotropic turbulence is performed by using the Spectral method at a Reynolds number Re = 1000 on a uniformly distributed grid points. The Reynolds number is low enough that the computational grid is capable of resolving all the possible turbulent scales. The statistical properties in the computed flow field show a good agreement with the qualitative behavior of decaying turbulence. The behavior of the flow structures in the computed flow field also follow the classical idea of the fluid flow in turbulence. Keywords: Direct numerical simulation, Isotropic turbulence, Spectral method. © 2013 JSR Publications. ISSN: 2070-0237 (Print); 2070-0245 (Online). All rights reserved. doi:http://dx.doi.org/10.3329/jsr.v5i3.12665 J. Sci. Res. 5 (3), 435-445 (2013)  


2020 ◽  
Vol 555 ◽  
pp. 124402
Author(s):  
Yuxian Xia ◽  
Xiang Qiu ◽  
Jianping Lou ◽  
Yuehong Qian

2018 ◽  
Vol 7 (4.35) ◽  
pp. 617
Author(s):  
P. Mathupriya ◽  
L. Chan ◽  
H. Hasini ◽  
A. Ooi

The numerical study of the flow over a two-dimensional cylinder which is symmetrically confined in a plane channel is presented to study the characteristics of vortex shedding. The numerical model has been established using direct numerical simulation (DNS) based on the open source computational fluid dynamics (CFD) code named OpenFOAM. In the present study, the flow fields have been computed at blockage ratio, β of 0.5 and at Reynolds number, Re of 200 and 300. Two-dimensional simulations investigated on the effects of Reynolds number based on the vortex formation and shedding frequency. It was observed that the presence of two distinct shedding frequencies appear at higher Reynolds number due to the confinement effects where there is strong interactions between boundary layer, shear layer and the wake of the cylinder. The range of simulations conducted here has shown to produce results consistent with that available in the open literature. Therefore, OpenFOAM is found to be able to accurately capture the complex physics of the flow.


Sign in / Sign up

Export Citation Format

Share Document