scholarly journals Lattice Boltzmann method for 2D flows in curvilinear coordinates

2012 ◽  
Vol 14 (3) ◽  
pp. 772-783 ◽  
Author(s):  
Ljubomir Budinski

In order to improve efficiency and accuracy, while maintaining an ease of modeling flows with the lattice Boltzmann approach in domains having complex geometry, a method for modeling equations of 2D flow in curvilinear coordinates has been developed. Both the transformed shallow water equations and the transformed 2D Navier-Stokes equations in the horizontal plane were synchronized with the equilibrium distribution function and the force term in the rectangular lattice. Since the solution of these equations takes place in the classical rectangular lattice environment, boundary conditions are modeled in the standard form of already existing simple methods (bounce-back), not requiring any additional functions. Owing to this and to the fact that the proposed method ensures a more accurate fitting of equations, even to domains of interest having complex geometry, the accuracy of solution is significantly increased, while the simplicity of the standard lattice Boltzmann approach is maintained. For the shallow water equations transformed in curvilinear coordinates, the proposed procedure is verified in three different hydraulic problems, all characterized by complex geometry.

2011 ◽  
Vol 675 ◽  
pp. 1-26 ◽  
Author(s):  
R. ROTUNNO ◽  
J. B. KLEMP ◽  
G. H. BRYAN ◽  
D. J. MURAKI

Nearly all analytical models of lock-exchange flow are based on the shallow-water approximation. Since the latter approximation fails at the leading edges of the mutually intruding fluids of lock-exchange flow, solutions to the shallow-water equations can be obtained only through the specification of front conditions. In the present paper, analytic solutions to the shallow-water equations for non-Boussinesq lock-exchange flow are given for front conditions deriving from free-boundary arguments. Analytic solutions are also derived for other proposed front conditions – conditions which appear to the shallow-water system as forced boundary conditions. Both solutions to the shallow-water equations are compared with the numerical solutions of the Navier–Stokes equations and a mixture of successes and failures is recorded. The apparent success of some aspects of the forced solutions of the shallow-water equations, together with the fact that in a real fluid the density interface is a free boundary, shows the need for an improved theory of lock-exchange flow taking into account non-hydrostatic effects for density interfaces intersecting rigid boundaries.


Author(s):  
S. A. Hosseini ◽  
N. Darabiha ◽  
D. Thévenin

With growing interest in the simulation of compressible flows using the lattice Boltzmann (LB) method, a number of different approaches have been developed. These methods can be classified as pertaining to one of two major categories: (i) solvers relying on high-order stencils recovering the Navier–Stokes–Fourier equations, and (ii) approaches relying on classical first-neighbour stencils for the compressible Navier–Stokes equations coupled to an additional (LB-based or classical) solver for the energy balance equation. In most cases, the latter relies on a thermal Hermite expansion of the continuous equilibrium distribution function (EDF) to allow for compressibility. Even though recovering the correct equation of state at the Euler level, it has been observed that deviations of local flow temperature from the reference can result in instabilities and/or over-dissipation. The aim of the present study is to evaluate the stability domain of different EDFs, different collision models, with and without the correction terms for the third-order moments. The study is first based on a linear von Neumann analysis. The correction term for the space- and time-discretized equations is derived via a Chapman–Enskog analysis and further corroborated through spectral dispersion–dissipation curves. Finally, a number of numerical simulations are performed to illustrate the proposed theoretical study. This article is part of the theme issue ‘Fluid dynamics, soft matter and complex systems: recent results and new methods’.


Author(s):  
Hiroshi Dan ◽  
◽  
Hiroshi Kanayama

In this paper, viscous shallow-water equations are derived from three-dimensional Navier-Stokes equations under the hydrostatic assumption. The viscous shallow-water equations are approximated by the finite element method based on our numerical scheme developed in 1978. This approach is used to simulate a tsunami in Hakata Bay. Results show a reasonable estimate of the tsunami arrival time.


2010 ◽  
Vol 20 (07) ◽  
pp. 1049-1087 ◽  
Author(s):  
BORIS HASPOT

In this paper, we consider the compressible Navier–Stokes equation with density-dependent viscosity coefficients and a term of capillarity introduced formally by van der Waals in Ref. 51. This model includes at the same time the barotropic Navier–Stokes equations with variable viscosity coefficients, shallow-water system and the model introduced by Rohde in Ref. 46. We first study the well-posedness of the model in critical regularity spaces with respect to the scaling of the associated equations. In a functional setting as close as possible to the physical energy spaces, we prove global existence of solutions close to a stable equilibrium, and local in time existence of solutions with general initial data. Uniqueness is also obtained.


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.


2017 ◽  
Vol 145 (3) ◽  
pp. 751-772 ◽  
Author(s):  
Michael D. Toy ◽  
Ramachandran D. Nair

An energy and potential enstrophy conserving finite-difference scheme for the shallow-water equations is derived in generalized curvilinear coordinates. This is an extension of a scheme formulated by Arakawa and Lamb for orthogonal coordinate systems. The starting point for the present scheme is the shallow-water equations cast in generalized curvilinear coordinates, and tensor analysis is used to derive the invariant conservation properties. Preliminary tests on a flat plane with doubly periodic boundary conditions are presented. The scheme is shown to possess similar order-of-convergence error characteristics using a nonorthogonal coordinate compared to Cartesian coordinates for a nonlinear test of flow over an isolated mountain. A linear normal mode analysis shows that the discrete form of the Coriolis term provides stationary geostrophically balanced modes for the nonorthogonal coordinate and no unphysical computational modes are introduced. The scheme uses centered differences and averages, which are formally second-order accurate. An empirical test with a steady geostrophically balanced flow shows that the convergence rate of the truncation errors of the discrete operators is second order. The next step will be to adapt the scheme for use on the cubed sphere, which will involve modification at the lateral boundaries of the cube faces.


Sign in / Sign up

Export Citation Format

Share Document