# Multispeed Lattice Boltzmann Model with Space-Filling Lattice for Transcritical Shallow Water Flows

Authors: Y. Peng; J. P. Meng; J. M. Zhang
Journal: Mathematical Problems in Engineering
 (2017)

---

## Abstract

Inspired by the recent success of applying multispeed lattice Boltzmann models with a non-space-filling lattice for simulating transcritical shallow water flows, the capabilities of their space-filling counterpart are investigated in this work. Firstly, two lattice models with five integer discrete velocities are derived by using the method of matching hydrodynamics moments and then tested with two typical 1D problems including the dam-break flow over flat bed and the steady flow over bump. In simulations, the derived space-filling multispeed models, together with the stream-collision scheme, demonstrate better capability in simulating flows with finite Froude number. However, the performance is worse than the non-space-filling model solved by finite difference scheme. The stream-collision scheme with second-order accuracy may be the reason since a numerical scheme with second-order accuracy is prone to numerical oscillations at discontinuities, which is worthwhile for further study.


---

## Body

## 1. Introduction

The shallow water equations (SWEs) have been used to model free surface flows in rivers and coastal areas under the assumption of the hydrostatic pressure [1]. Numerically, the SWEs can be solved by using conventional numerical method such as finite difference methods [2], finite element methods [3], and finite volume methods [4]. Alternatively, the lattice Boltzmann method (LBM) [5–8] can be used to model SWEs at mesoscopic level [1, 9–18]. In particular, one-dimensional shallow water flows have been studied by Frandsen [19], Thang et al. [20], and Chopard et al. [21].A major limitation of LBM is its inability to model supercritical flows. For this reason, Chopard et al. [21] developed an asymmetric lattice Boltzmann model for one-dimensional flow flows which can simulate flows with Froude numbers larger than 1. La Rocca et al. [22] proposed a multispeed model with a non-space-filling lattice which is solved by finite difference scheme and successfully simulated supercritical flows. Here, by using the word “multispeed,” it means lattices that have more than one nonzero speed in one dimension; see, for example, La Rocca et al. [22] and Brownlee et al. [23]. For instance, the commonly used D1Q3, D2Q9, and D3Q19 lattices are not classified as multispeed lattice. Also, “non-space-filling”/“space-filling” means a lattice model that cannot/can fit into the standard steam-collision scheme.Inspired by the success of the non-space-filling multispeed model, the capability of their space-filling counterpart will be investigated on simulating flows with finite Froude number. Specifically, two space-filling models will be derived by matching hydrodynamic moments (see, e.g., [1, 22, 24]) and then tested by using two typical shallow water problems.


## 2. Multispeed Lattice Boltzmann Models for 1D Shallow Water Equation

The one-dimensional SWEs read(1)∂h∂t+∂hu∂x=0,∂hu∂t+∂hu2∂x=-g∂∂xh22+ν∂2hu∂x2-gh∂zb∂x+F,which describe the evolution of water depth h and depth-averaged velocity u. The force term F is used to model various interesting effects such as wind-induced surface stress and the bed gradient.SWEs can be modelled by using a mesoscopic evolution rule(2)fαx+eαΔt,t+Δt-fαx,t=-1τfα-fαeq+ΔtSαx+12eαΔt,twhich describes the fluid motion using a distribution function f. Important factors, such as wind-induced surface stress, are also modelled by a force term at the right hand side; that is, Sα=wαeaiFi/(cs)2 and Fi=Fix+(1/2)eαΔt,t. The weight factor is denoted by wα for a discrete velocity eα, and the sound speed cs can be calculated by using e/λ, where λ is a reference quantity varying with lattice.To successfully simulate shallow water flows, the key is to define an appropriate local equilibrium functionfαeq and the associated lattice. Here, the equilibrium distribution function in Zhou [1] is generalized to(3)fαeq=w0h1w0-Aghe2-uiui2cs2α=0,wαhBghe2+eαiuics2+12eαieαjuiujcs4-uiuics2α≠0,where the coefficients A and B can be adjusted to satisfy the conservation of mass and momentum, i.e.,(4)h=∑αfαeq,hu2+gh22=∑αeefαeqαα.By substituting (3) into (4), the following equations can be obtained:(5)w0A=∑αwαBα≠0,∑αeewαBααe2=12α≠0.The weights wα in (3) will depend on the choice of lattice. Here two sets of velocity lattices will be adopted, that is, discrete velocities (0, ±1, ±2) derived by Qian and Zhou [25] and (0, ±1, ±3) derived by Chikatamarla and Karlin [26]. They are named as D1Q5A and D1Q5B, respectively, and the relevant parameters are listed in Table 1. For convenience, the standard (0, ±1) (D1Q3) lattice also is listed and it will be used for comparison in the following simulations.Table 1
Parameters of three lattice Boltzmann models. The reference quantity and the weights for D1Q5A and D1Q5B are adopted from Qian and Zhou [25] and Chikatamarla and Karlin [26], respectively.
D1Q3
D1Q5A
D1Q5B
A
3/4
1/2
(
110
+
23
10
)
/
48
B
3/2
1/2
(
5
-
10
)
/
6
λ
3
1
(
5
-
10
)
/
3
w
0
2/3
1/2
(
16
-
4
10
)
/
45
w
±
1
1/6
1/6
(
24
+
3
10
)
/
80
w
±
2 or w±3
—
1/12
(
16
+
5
10
)
/
720


## 3. Case  1: 1D Dam-Break Flow

To test the ability of the derived space-filling models, firstly a simple but challenging 1D dam-break flow over flat bed is used. Initially, the water is stationary but there is a discontinuity in the middle of computational region (x=1000 m), described as follows:(6)hx=h1,0≤x≤1000mh2,1000m<x≤2000m.The discontinuity will then propagate stream-wise and change its shape. The sharp gradient at the wave front is difficult to be captured smoothly.In simulations, 2000 cells are used and thereforeΔx is 1 m and Δt is 0.1 s. The initial upstream water depth is set to be 5 m for all tests. Firstly, in order to verify the two multispeed lattice models, the water depth of 3 m is specified in the downstream, where the maximum Froude (Fr) number is 0.26. In Figure 1, the predictions of water depth and velocity of two models and the D1Q3 model are compared with the corresponding analytical solutions [27]. It is shown that all three models are satisfactory, which is as expected since the Fr is relatively low.Figure 1
Water depth and velocity profiles predicted by D1Q3, D1Q5A, and D1Q5B models att=80 s (h1=5 m and h2=3).
(a)
(b)To examine the capability of simulating flow with higher Froude number, a series of tests will be carried out by adjusting the downstream depth. It is found, as shown in Figure2, that the D1Q5B model can simulate flows with Fr = 1.14 if the relaxation time is set to be 1.99, although there are errors near the discontinuity perhaps due to the numerical oscillations. The D1Q5A is able to simulate flows with Fr = 1.07, while the D1Q3 model fails to simulate Froude number beyond unity. The difference between two multispeed models shows that the capability can also be impacted by the lattice structure.Figure 2
Water depth and velocity profiles predicted by D1Q5B models att=80 s with h2=0.61 m and maximum Fr = 1.14 [(a) water depth, (b) water velocity].
(a)
(b)


## 4. Case  2: 1D Steady Flow over Bump

A steady flow over a bump bed will be utilized to further test the multispeed models. For this flow problem, the subcritical flow over the same bump has been studied by Goutal and Maurel [28], Zhou [1], and Peng et al. [18].For testing, the bed bottom is defined as(7)Zbx=0.2-0.05x-62,8m<x<12m0,else.The inlet boundary condition is specified by setting a constant discharge of Q=4.42 m2/s. The length of channel is 25 m. A fixed water depth h0 will be used at the outlet boundary, which will be reduced gradually to find the maximum Fr that can be simulated. The space step Δx is set to be 0.1 m to get the grid-independent solution according to Peng et al. [18] and the time step Δt=0.002 s. A relaxation time of 1.99 is employed for all simulations except for the cases shown in Figure 3. Besides, the bed friction is ignored here.Figure 3
Water surface profiles withh0=2 m.First, the case with a downstream water depth of 2 m is tested to verify the multispeed models. As shown in Figure3, all models agree well with the analytical solution, which again confirms capability of multispeed model in the relative low Froude number.Then the downstream water depth is adjusted to test the capability of higher Froude number, that is, to find the minimal downstream water depth. The corresponding results are presented in Figure4. It is found that the maximum Fr for D1Q3 is 0.69, which is consistent with the previous findings by Zhou [1]. The D1Q5B model can simulate flows with a Froude number of 1.09, while the D1Q5A model fails to simulate Froude number beyond 0.91, which also shows the importance of lattice structure.Figure 4
Water surface and Froude number profiles with minimum water depth for three models (h0=1.93, 1.77, and 1.71 and Fr=0.69, 0.91, and 1.09 for D1Q3, D1Q5A, and D1Q5B models, resp.).


## 5. Concluding Remarks

In this study, the space-filling multispeed lattice Boltzmann model has been investigated to simulate the transcritical shallow water flows. By matching hydrodynamics moments, two lattice models with five integer lattice velocities have been derived, which are tested together with the standard model with three discrete velocities for a 1D dam-break and a 1D steady flow over bump.It is encouraging to find that the space-filling multispeed model can improve the capability of simulating flows with higher Fr, particularly compared with the standard D1Q3 model. The different performances between two multispeed models show that capability can be impacted by the lattice structure.In comparison with the non-space-filling model proposed by La Rocca et al. [22], the space-filling versions here demonstrate limited capability in terms of finite Froude number although similar numbers of discrete velocities are used in these two studies. This is potentially due to the fact that the stream-collision scheme is a second-order scheme. It is well known that a numerical scheme of second-order accuracy is prone to numerical oscillations at discontinuities. The impact of the stream-collision scheme will be studied further in the future.



---

*Source: 8917360-2017-08-24.xml*