AN IMPROVED LATTICE KINETIC SCHEME FOR INCOMPRESSIBLE VISCOUS FLUID FLOWS

2013 ◽  
Vol 25 (01) ◽  
pp. 1340017 ◽  
Author(s):  
KOSUKE SUZUKI ◽  
TAKAJI INAMURO

The lattice Boltzmann method (LBM) is an explicit numerical scheme for the incompressible Navier–Stokes equations (INSE) without integrating the Poisson equation for the pressure. In spite of its merit, the LBM has some drawbacks in accuracy. First, we review drawbacks for three numerical methods based on the LBM. The three methods are the LBM with the Bhatnagar–Gross–Krook model (LBGK), the lattice kinetic scheme (LKS) and the link-wise artificial compressibility method (LWACM). Second, in order to remedy the drawbacks, we propose an improved LKS. The present method incorporates (i) the scheme used in the LWACM for determining the kinematic viscosity, (ii) an iterative calculation of the pressure and (iii) a semi-implicit algorithm, while preserving the simplicity of the algorithm of the original LKS. Finally, in simulations of test problems, we find that the improved LKS eliminates the drawbacks and gives more accurate and stable results than LBGK, LKS and LWACM.

2017 ◽  
Vol 20 (01) ◽  
pp. 1650064 ◽  
Author(s):  
Luigi C. Berselli ◽  
Stefano Spirito

We prove that suitable weak solutions of 3D Navier–Stokes equations in bounded domains can be constructed by a particular type of artificial compressibility approximation.


Author(s):  
Joris C. G. Verschaeve

By means of the continuity equation of the incompressible Navier–Stokes equations, additional physical arguments for the derivation of a formulation of the no-slip boundary condition for the lattice Boltzmann method for straight walls at rest are obtained. This leads to a boundary condition that is second-order accurate with respect to the grid spacing and conserves mass. In addition, the boundary condition is stable for relaxation frequencies close to two.


2010 ◽  
Vol 132 (10) ◽  
Author(s):  
Minoru Watari

Relationship between accuracy and number of velocity particles in velocity slip phenomena was investigated by numerical simulations and theoretical considerations. Two types of 2D models were used: the octagon family and the D2Q9 model. Models have to possess the following four prerequisites to accurately simulate the velocity slip phenomena: (a) equivalency to the Navier–Stokes equations in the N-S flow area, (b) conservation of momentum flow Pxy in the whole area, (c) appropriate relaxation process in the Knudsen layer, and (d) capability to properly express the mass and momentum flows on the wall. Both the octagon family and the D2Q9 model satisfy conditions (a) and (b). However, models with fewer velocity particles do not sufficiently satisfy conditions (c) and (d). The D2Q9 model fails to represent a relaxation process in the Knudsen layer and shows a considerable fluctuation in the velocity slip due to the model’s angle to the wall. To perform an accurate velocity slip simulation, models with sufficient velocity particles, such as the triple octagon model with moving particles of 24 directions, are desirable.


Water ◽  
2019 ◽  
Vol 11 (1) ◽  
pp. 82
Author(s):  
Getnet Kebede Demeke ◽  
Dereje Hailu Asfaw ◽  
Yilma Seleshi Shiferaw

Hydraulic structures are often complex and in many cases their designs require attention so that the flow behavior around hydraulic structures and their influence on the environment can be predicted accurately. Currently, more efficient computational fluid dynamics (CFD) codes can solve the Navier–Stokes equations in three-dimensions and free surface computation in a significantly improved manner. CFD has evolved into a powerful tool in simulating fluid flows. In addition, CFD with its advantages of lower cost and greater flexibility can reasonably predict the mean characteristics of flows such as velocity distributions, pressure distributions, and water surface profiles of complex problems in hydraulic engineering. In Ethiopia, Tendaho Dam Spillway was constructed recently, and one flood passed over the spillway. Although the flood was below the designed capacity, there was an overflow due to superelevation at the bend. Therefore, design of complex hydraulic structures using the state-of- art of 3D hydrodynamic modelling enhances the safety of the structures. 3D hydrodynamic modelling was used to verify the safety of the spillway using designed data and the result showed that the constructed hydraulic section is not safe unless it is modified.


1969 ◽  
Vol 37 (4) ◽  
pp. 727-750 ◽  
Author(s):  
Gareth P. Williams

A method of numerically integrating the Navier-Stokes equations for certain three-dimensional incompressible flows is described. The technique is presented through application to the particular problem of describing thermal convection in a rotating annulus. The equations, in cylindrical polar co-ordinate form, are integrated with respect to time by a marching process, together with the solving of a Poisson equation for the pressure. A suitable form of the finite difference equations gives a computationally-stable long-term integration with reasonably faithful representation of the spatial and temporal characteristics of the flow.Trigonometric interpolation techniques provide accurate (discretely exact) solutions to the Poisson equation. By using an auxiliary algorithm for rapid evaluation of trigonometric transforms, the proportion of computation needed to solve the Poisson equation can be reduced to less than 25% of the total time needed to’ advance one time step. Computing on a UNIVAC 1108 machine, the flow can be advanced one time-step in 2 sec for a 14 × 14 × 14 grid upward to 96 sec for a 60 × 34 × 34 grid.As an example of the method, some features of a solution for steady wave flow in annulus convection are presented. The resemblance of this flow to the classical Eady wave is noted.


2019 ◽  
Vol 17 (06) ◽  
pp. 977-1003
Author(s):  
H. B. de Oliveira

The purpose of this work is to study the generalized Navier–Stokes equations with nonlinear viscosity that, in addition, can be fully anisotropic. Existence of very weak solutions is proved for the associated initial and boundary-value problem, supplemented with no-slip boundary conditions. We show that our existence result is optimal in some directions provided there is some compensation in the remaining directions. A particular simplification of the problem studied here, reduces to the Navier–Stokes equations with (linear) anisotropic viscosity used to model either the turbulence or the Ekman layer in atmospheric and oceanic fluid flows.


1976 ◽  
Vol 98 (2) ◽  
pp. 214-222 ◽  
Author(s):  
J. E. Zweig ◽  
H. J. Sneck

The general hydrodynamic behavior at small clearance Reynolds numbers of two fluids of different density and viscosity occupying the finite annular space between a rotating and stationary disk is explored using a simplified version of the Navier-Stokes equations which retains only the centrifugal force portion of the inertia terms. A criterion for selecting the annular flow fields that are compatible with physical reservoirs is established and then used to determine the conditions under which two-fluid flows in the annulus might be expected for specific fluid combinations.


Sign in / Sign up

Export Citation Format

Share Document