Mathematical Problems in Engineering
Volume 2017 (2017), Article ID 8917360, 5 pages
https://doi.org/10.1155/2017/8917360

Research Article

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

Y. Peng,¹ J. P. Meng,² and J. M. Zhang¹

¹State Key Laboratory of Hydraulics and Mountain River Engineering, Sichuan University, Chengdu 610065, China
²Scientific Computing Department, STFC Daresbury Laboratory, Warrington WA4 4AD, UK

Correspondence should be addressed to J. M. Zhang; zhangjianmin@scu.edu.cn

Received 21 February 2017; Revised 7 June 2017; Accepted 5 July 2017; Published 24 August 2017

Academic Editor: Ling Qian

Copyright © 2017 Y. Peng et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

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.

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

which describe the evolution of water depth and depth-averaged velocity . The force term 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

which describes the fluid motion using a distribution function . Important factors, such as wind-induced surface stress, are also modelled by a force term at the right hand side; that is, and . The weight factor is denoted by for a discrete velocity , and the sound speed can be calculated by using , where is a reference quantity varying with lattice.

To successfully simulate shallow water flows, the key is to define an appropriate local equilibrium function and the associated lattice. Here, the equilibrium distribution function in Zhou [1] is generalized to

where the coefficients and can be adjusted to satisfy the conservation of mass and momentum, i.e.,

By substituting (3) into (4), the following equations can be obtained:

The weights 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

3/4
1/2

3/2
1/2

1

2/3
1/2

1/6
1/6

or
—
1/12

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 ( m), described as follows:

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 is 1 m and 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.

(a)

(b)

Figure 1: Water depth and velocity profiles predicted by D1Q3, D1Q5A, and D1Q5B models at  s ( m and ).

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 Figure 2, 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.

(a)

(b)

Figure 2: Water depth and velocity profiles predicted by D1Q5B models at  s with  m and maximum Fr = 1.14 [(a) water depth, (b) water velocity].

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

The inlet boundary condition is specified by setting a constant discharge of  m²/s. The length of channel is 25 m. A fixed water depth will be used at the outlet boundary, which will be reduced gradually to find the maximum Fr that can be simulated. The space step is set to be 0.1 m to get the grid-independent solution according to Peng et al. [18] and the time step  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 with  m.

First, the case with a downstream water depth of 2 m is tested to verify the multispeed models. As shown in Figure 3, 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 Figure 4. 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 (, 1.77, and 1.71 and , 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.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

Y. Peng would like to acknowledge the financial support of National Natural Science Foundation of China (Grant nos. 51409183, 51579166, and 51611130203). J. P. Meng is grateful to the support of the Open Fund from the State Key Laboratory of Hydraulics and Mountain River Engineering in Sichuan University, China (Grant no. SKHL1410). Both Y. Peng and J. P. Meng would like to acknowledge the support of International Exchange Programme of Royal Society, UK (IE151119) and NSFC, China (51611130203).

References

J. G. Zhou, Lattice Boltzmann Methods for Shallow Water Flows, Springer, Berlin, Germany, 2004.

V. Casulli and R. T. Cheng, “Semi‐implicit finite difference methods for three‐dimensional shallow water flow,” International Journal for Numerical Methods in Fluids, vol. 15, no. 6, pp. 629–648, 1992.

T. W. H. Sheu and C. C. Fang, “High resolution finite-element analysis of shallow water equations in two dimensions,” Computer Methods in Applied Mechanics and Engineering, vol. 190, no. 20-21, pp. 2581–2601, 2001.

A. H. Daoud, K. A. Rakha, and A. G. Abul-azm, “A two-dimensional finite volume hydrodynamic model for coastal areas: model development and validation,” Ocean Engineering, vol. 35, no. 1, pp. 150–164, 2008.

F. J. Higuera and S. Succi, “Simulating the flow around a circular cylinder with a lattice boltzmann equation,” Europhysics letters, vol. 8, no. 6, pp. 517–521, 1989.

R. Benzi, S. Succi, and M. Vergassola, “The lattice Boltzmann equation: theory and applications,” Physics Reports, vol. 222, no. 3, pp. 145–197, 1992.

S. Chen and G. D. Doolen, “Lattice Boltzmann method for fluid flows,” Annual Review of Fluid Mechanics, vol. 30, pp. 329–364, 1998.

C. K. Aidun and J. R. Clausen, “Lattice-boltzmann method for complex flows,” Annual Review of Fluid Mechanics, vol. 42, pp. 439–472, 2010.

R. Salmon, “The lattice Boltzmann method as a basis for ocean circulation modeling,” Journal of Marine Research, vol. 57, no. 3, pp. 503–535, 1999.

P. J. Dellar, “Nonhydrodynamic modes and a priori construction of shallow water lattice Boltzmann equations,” Physical Review E—Statistical, Nonlinear, and Soft Matter Physics, vol. 65, no. 3, Article ID 036309, pp. 1–12, 2002.

H. Liu, J. G. Zhou, and R. Burrows, “Numerical modeling of turbulent compound channel flow using the lattice Boltzmann method,” International Journal for Numerical Methods in Fluids, vol. 59, no. 7, pp. 753–765, 2009.

H. Liu, J. G. Zhou, and R. Burrows, “Multi-block lattice Boltzmann simulations of subcritical flow in open channel junctions,” Computers and Fluids, vol. 38, no. 6, pp. 1108–1117, 2009.

C. Biscarini, S. D. Francesco, and M. Mencattini, “Application of the lattice Boltzmann method for large-scale hydraulic problems,” International Journal of Numerical Methods for Heat and Fluid Flow, vol. 21, no. 5, pp. 584–601, 2011.

K. R. Tubbs and F. T. Tsai, “GPU accelerated lattice Boltzmann model for shallow water flow and mass transport,” International Journal for Numerical Methods in Engineering, vol. 86, no. 3, pp. 316–334, 2011.

C. Biscarini, S. D. Francesco, F. Nardi, and P. Manciola, “Detailed simulation of complex hydraulic problems with macroscopic and mesoscopic mathematical methods,” Mathematical Problems in Engineering, vol. 2013, Article ID 928309, 14 pages, 2013.

Y. Peng, J. G. Zhou, and R. Burrows, “Modelling solute transport in shallow water with the lattice Boltzmann method,” Computers and Fluids, vol. 50, pp. 181–188, 2011.

Y. Peng, J. G. Zhou, J. M. Zhang, and H. Liu, “Lattice Boltzmann modeling of shallow water flows over discontinuous beds,” International Journal for Numerical Methods in Fluids, vol. 75, no. 8, pp. 608–619, 2014.

Y. Peng, J. M. Zhang, and J. G. Zhou, “Lattice boltzmann model using two relaxation times for shallow-water equations,” Journal of Hydraulic Engineering, vol. 142, no. 2, Article ID 06015017, 2015.

J. B. Frandsen, “A simple LBE wave runup model,” Progress in Computational Fluid Dynamics, vol. 8, no. 1-4, pp. 222–232, 2008.

P. V. Thang, B. Chopard, L. Lefévre, D. A. Ondo, and E. Mendes, “Study of the 1D lattice Boltzmann shallow water equation and its coupling to build a canal network,” Journal of Computational Physics, vol. 229, no. 19, pp. 7373–7400, 2010.

B. Chopard, P. V. Thang, and L. Lefèvre, “Asymmetric lattice Boltzmann model for shallow water flows,” Computers and Fluids, vol. 88, pp. 225–231, 2013.

M. La Rocca, A. Montessori, P. Prestininzi, and S. Succi, “A multispeed Discrete Boltzmann Model for transcritical 2D Shallow water flows,” Journal of Computational Physics, vol. 284, pp. 117–132, 2015.

R. A. Brownlee, A. N. Gorban, and J. Levesley, “Stable multispeed lattice Boltzmann methods,” Tech. Rep. MA-06-026, University of Leicester, Leicester, UK, 2006, https://arxiv.org/abs/cond-mat/0611616.

H. W. Zheng, C. Shu, Y. T. Chew, and J. Qiu, “A platform for developing new lattice Boltzmann models,” International Journal of Modern Physics C. Computational Physics and Physical Computation, vol. 16, no. 1, pp. 61–84, 2005.

Y.-H. Qian and Y. Zhou, “Complete Galilean-invariant lattice BGK models for the Navier-Stokes equation,” Europhysics Letters, vol. 42, no. 4, pp. 359–364, 1998.

S. S. Chikatamarla and I. V. Karlin, “Entropy and Galilean invariance of lattice Boltzmann theories,” Physical Review Letters, vol. 97, no. 19, Article ID 190601, 2006.

O. Delestre, C. Lucas, P.-A. Ksinant et al., “SWASHES: a compilation of shallow water analytic solutions for hydraulic and environmental studies,” International Journal for Numerical Methods in Fluids, vol. 72, no. 3, pp. 269–300, 2013.

N. Goutal and F. Maurel, Proceedings of the 2nd Workshop on Dam-break Wave Simulation, Tech. Rep. HE-43/97/016/B, Département Laboratoire National d’Hydraulique, Groupe Hydraulique Fluviale, Electricité de France, France, 1997.