Development of Lattice Boltzmann Flux Solver for Simulation of Incompressible Flows

2014 ◽  
Vol 6 (4) ◽  
pp. 436-460 ◽  
Author(s):  
C. Shu ◽  
Y. Wang ◽  
C. J. Teo ◽  
J. Wu

AbstractA lattice Boltzmann flux solver (LBFS) is presented in this work for simulation of incompressible viscous and inviscid flows. The new solver is based on Chapman-Enskog expansion analysis, which is the bridge to link Navier-Stokes (N-S) equations and lattice Boltzmann equation (LBE). The macroscopic differential equations are discretized by the finite volume method, where the flux at the cell interface is evaluated by local reconstruction of lattice Boltzmann solution from macroscopic flow variables at cell centers. The new solver removes the drawbacks of conventional lattice Boltzmann method such as limitation to uniform mesh, tie-up of mesh spacing and time interval, limitation to viscous flows. LBFS is validated by its application to simulate the viscous decaying vortex flow, the driven cavity flow, the viscous flow past a circular cylinder, and the inviscid flow past a circular cylinder. The obtained numerical results compare very well with available data in the literature, which show that LBFS has the second order of accuracy in space, and can be well applied to viscous and inviscid flow problems with non-uniform mesh and curved boundary.

2015 ◽  
Vol 18 (3) ◽  
pp. 593-620 ◽  
Author(s):  
Yan Wang ◽  
Chang Shu ◽  
Chiang Juay Teo ◽  
Jie Wu ◽  
Liming Yang

AbstractA three-dimensional (3D) lattice Boltzmann flux solver (LBFS) is presented in this paper for the simulation of both isothermal and thermal flows. The present solver combines the advantages of conventional Navier-Stokes (N-S) solvers and lattice Boltzmann equation (LBE) solvers. It applies the finite volume method (FVM) to solve the N-S equations. Different from the conventional N-S solvers, its viscous and inviscid fluxes at the cell interface are evaluated simultaneously by local reconstruction of LBE solution. As compared to the conventional LBE solvers, which apply the lattice Boltzmann method (LBM) globally in the whole computational domain, it only applies LBM locally at each cell interface, and flow variables at cell centers are given from the solution of N-S equations. Since LBM is only applied locally in the 3D LBFS, the drawbacks of the conventional LBM, such as limitation to uniform mesh, tie-up of mesh spacing and time step, tedious implementation of boundary conditions, are completely removed. The accuracy, efficiency and stability of the proposed solver are examined in detail by simulating plane Poiseuille flow, lid-driven cavity flow and natural convection. Numerical results show that the LBFS has a second order of accuracy in space. The efficiency of the LBFS is lower than LBM on the same grids. However, the LBFS needs very less non-uniform grids to get grid-independence results and its efficiency can be greatly improved and even much higher than LBM. In addition, the LBFS is more stable and robust.


AIAA Journal ◽  
1985 ◽  
Vol 23 (4) ◽  
pp. 583-587 ◽  
Author(s):  
Ajay Kumar ◽  
M. D. Salas

1980 ◽  
Vol 31 (1) ◽  
pp. 42-55 ◽  
Author(s):  
H.N.V. Dutt ◽  
A.K. Sreekanth

SummaryA design procedure has been developed to generate aerofoil shapes for prescribed pressure distributions in an incompressible viscous attached flow. It is based on the method of singularities, originally proposed by Chen and later modified by Kennedy and Marsden, for inviscid flows. The classical approach of adding the displacement thickness of the boundary layer and wake to the aerofoil contour is used to account for viscous effects. Several numerical examples are worked out and are compared with the inviscid flow results. Significant changes in aerofoil contours due to viscous effects are observed and these are discussed.


2011 ◽  
Vol 105-107 ◽  
pp. 2307-2310
Author(s):  
Jian Ping Yu ◽  
Shu Rong Yu ◽  
Xing Wang Liu

Lattice Boltzmann methods (LBM) have become an alternative to conventional computational fluid dynamics (CFD) methods for various systems. In this paper, flow field of mean flow past a circular cylinder was simulated based on the lattice Boltzmann method. The streamline of air past the cylinder illuminated that the fluid adhere on the boundary and doesn’t separate from the surface of cylinder when Re number less than 5. When Re number equal 40, flow separated to form a pair of recirculating eddies can be observed. With the Re number increasing, the trailing vortex length is growth accordingly. When Re number come up to 80, the trailing vortex begin to shed regularly. This result is consistent with the experiment data. Drag coefficient that fluid act on the surface of cylinder was calculated. The calculated results were same as the experiment data. Simulation indicate that LBM can simulate the vortex taking place and shedding effectively.


1990 ◽  
Vol 220 ◽  
pp. 459-484 ◽  
Author(s):  
H. M. Badr ◽  
M. Coutanceau ◽  
S. C. R. Dennis ◽  
C. Ménard

The unsteady flow past a circular cylinder which starts translating and rotating impulsively from rest in a viscous fluid is investigated both theoretically and experimentally in the Reynolds number range 103 [les ] R [les ] 104 and for rotational to translational surface speed ratios between 0.5 and 3. The theoretical study is based on numerical solutions of the two-dimensional unsteady Navier–Stokes equations while the experimental investigation is based on visualization of the flow using very fine suspended particles. The object of the study is to examine the effect of increase of rotation on the flow structure. There is excellent agreement between the numerical and experimental results for all speed ratios considered, except in the case of the highest rotation rate. Here three-dimensional effects become more pronounced in the experiments and the laminar flow breaks down, while the calculated flow starts to approach a steady state. For lower rotation rates a periodic structure of vortex evolution and shedding develops in the calculations which is repeated exactly as time advances. Another feature of the calculations is the discrepancy in the lift and drag forces at high Reynolds numbers resulting from solving the boundary-layer limit of the equations of motion rather than the full Navier–Stokes equations. Typical results are given for selected values of the Reynolds number and rotation rate.


2010 ◽  
Vol 132 (12) ◽  
Author(s):  
Sunil Lakshmipathy ◽  
Sharath S. Girimaji

The objective of this study is to evaluate the capability of the partially averaged Navier–Stokes (PANS) method in a moderately high Reynolds number (ReD 1.4×105) turbulent flow past a circular cylinder. PANS is a bridging closure model purported for use at any level of resolution ranging from Reynolds-averaged Navier–Stokes to direct numerical simulations. The closure model is sensitive to the length-scale cut-off via the ratios of unresolved-to-total kinetic energy (fk) and unresolved-to-total dissipation (fε). Several simulations are performed to study the effect of the cut-off length-scale on computed closure model results. The results from various resolutions are compared against experimental data, large eddy simulation, and detached eddy simulation solutions. The quantities examined include coefficient of drag (Cd), Strouhal number (St), and coefficient of pressure distribution (Cp) along with the mean flow statistics and flow structures. Based on the computed results for flow past circular cylinder presented in this paper and analytical attributes of the closure model, it is reasonable to conclude that the PANS bridging method is a theoretically sound and computationally viable variable resolution approach for practical flow computations.


1995 ◽  
Vol 303 ◽  
pp. 215-232 ◽  
Author(s):  
H. M. Badr ◽  
S. C. R. Dennis ◽  
S. Kocabiyik ◽  
P. Nguyen

The transient flow field caused by an infinitely long circular cylinder placed in an unbounded viscous fluid oscillating in a direction normal to the cylinder axis, which is at rest, is considered. The flow is assumed to be started suddenly from rest and to remain symmetrical about the direction of motion. The method of solution is based on an accurate procedure for integrating the unsteady Navier–Stokes equations numerically. The numerical method has been carried out for large values of time for both moderate and high Reynolds numbers. The effects of the Reynolds number and of the Strouhal number on the laminar symmetric wake evolution are studied and compared with previous numerical and experimental results. The time variation of the drag coefficients is also presented and compared with an inviscid flow solution for the same problem. The comparison between viscous and inviscid flow results shows a better agreement for higher values of Reynolds and a Strouhal numbers. The mean flow for large times is calculated and is found to be in good agreement with previous predictions based on boundary-layer theory.


2012 ◽  
Vol 9 (2) ◽  
pp. 135-152 ◽  
Author(s):  
Sreekar Gomatam ◽  
S Vengadesan ◽  
S K Bhattacharyya

Three dimensional (3D) flow past an Autonomous Underwater Vehicle (AUV) is simulated using a Computational Fluid Dynamics (CFD) approach at a Reynolds (Re) number of 2.09x106. A non-linear k-? (NLKE) turbulence model is used for solving the Reynolds Averaged Navier-Stokes (RANS) equations. The effect of control surfaces over the flow, the flow interaction between the hull and the appendages at various Angles of Attack (AoA) and the effect of the symmetry plane is studied. Flow structure, variation of flow variables and force distribution for various AoA are presented and discussed in detail.DOI: http://dx.doi.org/10.3329/jname.v9i2.12567 Journal of Naval Architecture and Marine Engineering 9(2012) 135-152


1960 ◽  
Vol 9 (4) ◽  
pp. 613-620 ◽  
Author(s):  
R. J. Swigart

The inviscid flow behind a cylindrical blast wave and its analogy with hypersonic flow past blunt-nosed cylinders is considered. Sakurai (1953, 1954) obtained a solution for the flow field behind a propagating blast wave by expanding the flow variables in power series of 1/M2, where M is the blast wave Mach number, and determining the coefficients of the first two terms in the series. Here the work is extended to include third-order terms. Third-order theory is shown to improve the prediction of shock wave shapes and surface pressure distribution on hemisphere-cylinder configurations at M∞ = 7·7 and 17·18.


Sign in / Sign up

Export Citation Format

Share Document