Next Article in Journal
A Novel Lightweight Real-Time Traffic Sign Detection Integration Framework Based on YOLOv4
Previous Article in Journal
An Information-Theoretic Analysis of the Cost of Decentralization for Learning and Inference under Privacy Constraints
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

General Propagation Lattice Boltzmann Model for the Boussinesq Equation

School of Mathematics and Information Science, North Minzu University, Yinchuan 750021, China
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Entropy 2022, 24(4), 486; https://doi.org/10.3390/e24040486
Submission received: 14 March 2022 / Revised: 27 March 2022 / Accepted: 28 March 2022 / Published: 30 March 2022

Abstract

:
A general propagation lattice Boltzmann model is used to solve Boussinesq equations. Different local equilibrium distribution functions are selected, and the macroscopic equation is recovered with second order accuracy by means of the Chapman–Enskog multi-scale analysis and the Taylor expansion technique. To verify the effectiveness of the present model, some Boussinesq equations with initial boundary value problems are simulated. It is shown that our model can remain stable and accurate, which is an effective algorithm worthy of promotion and application.

1. Introduction

Lattice Boltzmann method (LBM) is a popular numerical method developed in recent decades [1,2,3], which is based on kinetic theory. In the lattice Boltzmann method, in addition to the fluid being discretized into fluid particles, the physical region is also discretized into a series of lattices, and the time is discretized into a series of time steps. It has unique advantages, such as easy implementation, simple boundary condition processing, and easy parallel computing [4]. It can be used to simulate some fluid flow [5,6,7] and solve partial differential equations [8,9,10,11,12].
The Boussinesq equation [13,14] describes the motions of long waves in one-dimensional nonlinear lattices and in shallow water under gravity. The numerical solutions of the Boussinesq equation are extensively studied [15,16]. Ref. [17] proposes a lattice Boltzmann model with an optimization term for the generalized Boussinesq equation. Ref. [18] simulates Boussinesq equations with source terms.
A general propagation lattice Boltzmann (GPLB) scheme is more general than the standard lattice Bhatnagar–Gross–Krook (SLBGK) model [19]. Some studies on GPLB models are being promoted [9,12,19]. The stability and accuracy of simulating equations will be improved by GPLB models. In this paper, we develop a GPLB model to solve Boussinesq equations. The effectiveness and stability of our model are verified by comparing the simulation results and the exact solutions.
An outline of our paper is given by: in Section 2, we derive a GPLB model for the generalized nonlinear Boussinesq equation. In Section 3, numerical simulations are performed. Finally, conclusions will be summarized in Section 4.

2. GPLB Model for Boussinesq Equations

This paper studies Boussinesq equations with the following forms,
u t t = α u x x + p u x x + β u x x x x .
where the macro variable u = u ( x , t ) represents the free movement of the fluid surface. The depth of the flowing fluid and the characteristic velocity of the corresponding long wave determined α , β , respectively. Equation (1) is a well-known generalized Boussinesq equation.
For Equation (1), the evolution law of the particle distribution function can be the corresponding discrete velocity Boltzmann equation [20] with the Bhatnagar–Gross–Krook (BGK) collision operator [1],
f α t + ξ α f α x = 1 τ 0 [ f α f α e q ] ,
where f α ( x , t ) is a scalar function describing the particle distribution at position x and time t, ξ α , α = 0 , 1 , , n 1 is the set of discrete velocities. f α e q is the local equilibrium distribution function, and τ 0 is the single relaxation time.
In the LBM, the D d Q b model is often used to represent the dimensions and speed of the problem. Among them, d represents the dimension of the problem (1 represents one-dimension, 2 represents two-dimensions, 3 represents three-dimensions), b represents the number of lattice chains in the velocity model. In this paper, we use the D 1 Q 5 velocity model, in which the discrete velocities can be defined as,
ξ α = c e α = c 0 , 1 , 1 , 2 , 2 = 0 , c , c , 2 c , 2 c ,
where c = k Δ x / Δ t , representing the propagation speed of f α along the lattice chain. The discrete lattice time step and space step are Δ t and Δ x , respectively. We use k to change the propagation process, which is an important value in the current method.
Equation (2) can be decomposed into collision and propagation steps for each Δ t with applying the time-splitting method [19],
f α t = 1 τ 0 [ f α f α e q ]
f α t + e α f α x = 0 .
We can choose the appropriate numerical schemes method according to the features of Equations (4) and (5).

2.1. GPLB Model for Boussinesq Equations

Since there is no spatial derivative term, Equation (4) is discretized into the following form with the explicit Euler scheme,
f α ( x , t ) = ( 1 1 τ ) f α ( x , t ) + 1 τ f α e q ( x , t ) .
We can see that there is no difference between the collision process and that in the SLBGK models.
We treat Equation (5) as follows,
f α ( x , t + Δ t ) = m 0 f α ( x , t ) + m 1 f α ( x S i , t ) + m 1 f α ( x + S i , t ) , S i = Δ x · e i ,
where m 0 , m 1 and m 1 are free parameters satisfying,
m 0 + m 1 + m 1 = 1 , m 1 m 1 = Δ t · ξ α S i = Δ t c e α Δ x e α = k .
Thus, we can get from Equation (8),
m 0 = 1 n , m 1 = n + k 2 , m 1 = n k 2 .
We introduce one parameter n for Equation (9). From Equation (7), we improve the propagation process. Substituting Equation (9) into Equation (7), we can obtain:
f α ( x , t + Δ t ) = f α ( x , t ) + n 2 f α ( x + S i , t ) 2 f α ( x , t ) + f α ( x S i , t ) k 2 f α ( x + S i , t ) f α ( x S i , t ) .
We choose k 2 n 1 to ensure simulation stability [9].
So far, Equations (6) and (10) constitute the GPLB model for Equation (1). Next, we will recover Equation (1).

2.2. Recovery of Boussinesq Equations

To complete our general propagation lattice Boltzmann model for Equation (1), the multi-scale Chapman–Enskog [21] and Taylor expansions will be applied to obtain the specific expressions of the local equilibrium distribution function f α e q .
Firstly, applying the Taylor expansion to f α ( x + S i , t ) and f α ( x S i , t ) , and retaining the terms up to O ( Δ t 5 ) , we have:
f α ( x + S i , t ) = f α ( x , t ) + S i · x f α ( x , t ) + ( S i · x ) 2 2 f α ( x , t ) + ( S i · x ) 3 6 f α ( x , t ) + ( S i · x ) 4 24 f α ( x , t ) + O ( S i 5 ) = f α ( x , t ) + Δ t k ( ξ α · x ) f α ( x , t ) + Δ t 2 2 k 2 ( ξ α · x ) 2 f α ( x , t ) + Δ t 3 6 k 3 ( ξ α · x ) 3 f α ( x , t ) + Δ t 4 24 k 4 ( ξ α · x ) 4 f α ( x , t ) + O ( Δ t 5 ) ,
f α ( x S i , t ) = f α ( x , t ) Δ t k ( ξ α · x ) f α ( x , t ) + Δ t 2 2 k 2 ( ξ α · x ) 2 f α ( x , t ) Δ t 3 6 k 3 ( ξ α · x ) 3 f α ( x , t ) + Δ t 4 24 k 4 ( ξ α · x ) 4 f α ( x , t ) + O ( Δ t 5 ) .
Using Equations (10)–(12), we have:
f α ( x , t + Δ t ) = f α ( x , t ) Δ t ( ξ α · x ) f α ( x , t ) + Δ t 2 n 2 k 2 ( ξ α · x ) 2 f α ( x , t ) Δ t 3 6 k 3 ( ξ α · x ) 3 f α ( x , t ) + Δ t 4 n 24 k 4 ( ξ α · x ) 4 f α ( x , t ) + O ( Δ t 5 ) .
Taylor expansion of Equation (13) to O ( Δ t 5 ) ,
f α ( x , t + Δ t ) = f α ( x , t ) + Δ t t f α ( x , t ) + Δ t 2 2 t 2 f α ( x , t ) + Δ t 3 6 t 3 f α ( x , t ) + Δ t 4 24 t 4 f α ( x , t ) + O ( Δ t 5 ) ,
we have:
f α ( x , t ) + Δ t t f α ( x , t ) + Δ t 2 2 t 2 f α ( x , t ) + Δ t 3 6 t 3 f α ( x , t ) + + Δ t 4 24 t 4 f α ( x , t ) + O ( Δ t 5 ) = f α ( x , t ) Δ t ( ξ α · x ) f α ( x , t ) + Δ t 2 n 2 k 2 ( ξ α · x ) 2 f α ( x , t ) Δ t 3 6 k 3 ( ξ α · x ) 3 f α ( x , t ) + Δ t 4 n 24 k 4 ( ξ α · x ) 4 f α ( x , t ) + O ( Δ t 5 ) .
Multi-scale Chapman–Enskog expansion of the system in the following forms,
x = ε x 1 ,
t = ε t 1 + ε 2 t 2 + ε 3 t 2 + ε 4 t 4 ,
f α = n = 0 ϵ n f α ( n ) = f α 0 + ϵ f α ( 1 ) + ϵ 2 f α ( 2 ) + ϵ 3 f α ( 3 ) + .
Using Equations (6), (15)–(18), we obtain:
O ( ε 0 ) : f α 0 = ( 1 1 τ ) f α 0 + 1 τ f α e q , i . e . , f α 0 = f α e q ,
O ( ε 1 ) : f α 1 + Δ t t 1 f α 0 = ( 1 1 τ ) f α 1 Δ t ( ξ α · x 1 ) f α 0 ,
O ( ε 2 ) : f α 2 + Δ t t 2 f α 0 + Δ t t 1 f α 1 + Δ t 2 2 t 1 2 f α 0 = ( 1 1 τ ) f α 2 + n Δ t 2 2 k 2 ( ξ α · x 1 ) 2 f α 0 ,
O ( ε 3 ) : f α 3 + Δ t t 1 f α 2 + Δ t t 2 f α 1 + Δ t t 3 f α 0 + Δ t 2 2 t 1 2 f α 1 + Δ t 3 6 t 1 3 f α 0 = ( 1 1 τ ) f α 3 + Δ t ( ξ α · x 1 ) ( 1 1 τ ) f α 2 Δ t 3 6 k 3 ( ξ α · x 1 ) 3 f α 0 ,
O ( ε 4 ) : f α 4 + Δ t ( t 1 f α 3 + t 2 f α 2 + t 3 f α 1 + t 4 f α 0 ) + Δ t 2 2 ( t 1 2 f α 2 + t 2 2 f α 0 + 2 t 1 , t 2 2 f α 1 + 2 t 1 , t 3 2 f α 0 ) + Δ t 3 6 ( t 1 3 f α 1 + 3 t 2 t 1 2 f α 0 ) + Δ t 4 24 t 1 4 f α 0 = ( 1 1 τ ) f α 4 Δ t ( ξ α · x 1 ) ( 1 1 τ ) f α 3 + n Δ t 2 2 k 2 ( ξ α · x 1 ) 2 ( 1 1 τ ) f α 2 + n Δ t 4 24 k 4 ( ξ α · x 1 ) 4 f α 0 .
Simplifying Equations (19)–(23), we have:
f α 1 τ Δ t =   t 1 + ( ξ α · x 1 ) f α 0 ,
f α 2 τ Δ t =   Δ t 2 t 1 2 + t 2 n Δ t 2 k 2 ( ξ α · x 1 ) 2 f α 0 + t 1 + ( 1 1 τ ) ( ξ α · x 1 ) f α 1 ,
f α 3 τ Δ t =   Δ t 2 6 t 1 3 + Δ t t 1 , t 2 2 + t 3 + Δ t 2 6 k 2 ( ξ α · x 1 ) 3 f α 0 + Δ t 2 t 1 2 + t 2 n Δ t 2 k 2 ( 1 1 τ ) ( ξ α · x 1 ) 2 f α 1 + t 1 + ( 1 1 τ ) ( ξ α · x 1 ) f α 2 ,
f α 4 τ Δ t =   t 4 + Δ t 2 t 2 2 + Δ t t 1 , t 3 2 + Δ t 2 2 t 2 t 1 2 + Δ t 3 24 t 1 4 n Δ t 3 24 k 4 ( ξ α · x 1 ) 4 f α 0 + t 3 + Δ t t 1 , t 2 2 + Δ t 2 6 t 1 3 + Δ t 2 6 k 3 ( ξ α · x 1 ) 3 ( 1 1 τ ) f α 1 + t 2 + Δ t 2 t 1 2 n Δ t 2 k 2 ( ξ α · x 1 ) 2 ( 1 1 τ ) f α 2 + t 1 Δ t ( ξ α · x 1 ) ( 1 1 τ ) f α 3 .
From Equation (24), one can obtain:
f α 1 = τ Δ t t 1 + ( ξ α · x 1 ) f α 0 .
Substituting Equation (28) into (25), we have:
f α 2 = τ Δ t 2 ( 1 2 τ ) t 1 2 + 1 Δ t t 2 + ( 1 2 τ ) t 1 ( ξ α · x 1 ) τ 1 ( ξ α · x 1 ) 2 f α 0 .
Coupling Equations (26), (28) and (29), we obtain:
f α 3 = τ Δ t 3 τ 2 τ + 1 6 t 1 3 + 1 Δ t ( 1 2 τ ) t 1 , t 2 2 + 1 Δ t 2 t 3 + 3 τ 2 3 τ 2 1 t 1 2 + 1 Δ t ( 1 2 τ ) t 2 ξ α · x 1 + 3 τ 2 + n k 2 4 τ + 1 n 2 k 2 t 1 ξ α · x 1 2 + n k 2 1 τ + 1 n k 2 + 1 6 k 2 ξ α · x 1 3 f α 0 .
Lastly, substituting Equations (28)–(30) into Equation (27), we have:
f α 4 = τ Δ t 4 1 Δ t 3 t 4 + 1 2 Δ t 2 τ Δ t 2 t 2 2 + 1 Δ t 2 2 τ Δ t 2 t 1 , t 3 2 + 1 2 Δ t 3 τ Δ t + 3 τ 2 Δ t t 2 t 1 2 + τ 2 t 1 4 + 1 Δ t τ Δ t 2 τ + 1 Δ t t 3 + τ 3 t 1 , t 2 2 + τ 4 t 1 3 ξ α · x 1 + τ 5 ξ α · x 1 4 + τ 6 t 1 2 ξ α · x 1 3 + ( τ 7 t 1 2 + τ 8 t 2 ) ( ξ α · x 1 ) 2 f α 0 ,
where:
τ 1 = τ 1 + n 2 k 2 , τ 2 = τ 3 + 3 2 τ 2 Δ t 6 + 5 12 τ + 1 24 , τ 3 = 2 + 4 Δ t τ 2 3 Δ t + 3 τ + 1 , τ 4 = τ 3 Δ t + 2 τ 2 Δ t 7 Δ t 6 τ + Δ t 6 , τ 5 = n 2 k 2 n Δ t k 2 + Δ t τ 2 + 1 6 k 3 + n k 2 n 2 4 k 4 + 2 n Δ t k 2 2 Δ t Δ t 6 k 2 τ n 24 k 4 + 1 6 k 3 n 2 k 2 + n 2 4 k 4 n Δ t k 2 + Δ t 6 k 2 + Δ t , τ 6 = 3 Δ t τ 3 + 2 n k 2 + 7 Δ t n Δ t k 2 + 1 τ 2 + 5 n 2 k 2 2 6 k 3 + 3 n Δ t k 2 Δ t 1 + 1 6 k 3 n 2 k 2 3 Δ t n Δ t 2 k 2 , τ 7 = 9 2 3 n 2 k 2 9 Δ t 2 τ 2 + 3 Δ t 3 τ 3 + 3 n 2 k 2 + Δ t 2 3 2 τ + Δ t n 4 k 2 , τ 8 = 2 + 1 Δ t τ 2 + 3 + 1 Δ t + n k 2 Δ t τ + 1 n 2 k 2 Δ t .
f α and f α e q follow:
α f α = α f α e q = u t .
From Equation (19), we have:
α f α 0 = u t , α f α n = 0 , n > 0 .
In order to recover Equation (1) and satisfy the solvability, based on ensuring that the truncation error is minimized, it is determined that f α and f α e q should satisfy the following moment conditions:
α ξ α f α 0 = 0 ,
α ξ α 2 f α 0 = α u + p ( u ) τ 1 Δ t ,
α ξ α 3 f α 0 = 0 ,
α ξ α 4 f α 0 = β u τ 5 Δ t 3 .
Summing Equation (28) over α , and coupling with Equations (32) and (33), we obtain:
t 1 u t   = 0 .
Summing Equation (29) over α and coupling with Equations (33)–(35), we obtain:
t 2 u t   = 1 ε 2 α u x x + p ( u ) x x .
Summing Equation (30) over α and coupling with Equations (33)–(36), we obtain:
t 3 u t   = Δ t 2 ε 2 3 τ 2 + n k 2 4 τ + 1 n 2 k 2 t 1 α u x x + p ( u ) x x Δ t τ 1 .
Summing Equation (31) over α and using Equations (33)–(37), we obtain:
t 4 u t   = Δ t 3 ε 2 1 2 Δ t 2 τ Δ t 2 α u x x + p ( u ) x x + 1 ε 4 x 4 β u Δ t 3 ε 2 τ 8 t 2 α u x x + p ( u ) x x Δ t τ 1 .
Taking Equation (38) × ε + Equation (39) × ε 2 + Equation (40) × ε 3 + Equation (41) × ε 4 and assuming ε = Δ t , macroscopic equations, i.e., Equation (1), can be recovered with second order accuracy,
u t t α u x x p u x x β u x x x x = R .
Through error analysis, we have:
R = ε 2 τ 1 3 τ 2 + n k 2 4 τ + 1 n 2 k 2 α t 1 u x x + t 1 p ( u ) x x ε 3 1 2 τ α t 2 u x x + t 2 p ( u ) x x ε 3 τ 1 2 ε + 1 τ 2 + 3 ε + 1 + n k 2 ε τ + ε n 2 k 2 α t 2 u x x + t 2 p ( u ) x x = O ( ε 2 ) .

2.3. Equilibrium Distribution Functions

In this paper, the D 1 Q 5 velocity model is used, i.e., b = 5 . The discrete velocity set is ξ α = 0 , c , c , 2 c , 2 c . Coupling with Equations (34)–(37), we can derive the equilibrium distribution functions,
  f 0 0 = u t 5 4 α u + p ( u ) c 2 Δ t τ 1 + 1 4 β u c 4 Δ t τ 5 ,
f 1 0 = 2 3 α u + p ( u ) c 2 Δ t τ 1 1 6 β u c 4 Δ t 3 τ 5 ,
f 2 0 = 2 3 α u + p ( u ) c 2 Δ t τ 1 1 6 β u c 4 Δ t 3 τ 5 ,
f 3 0 = 1 24 α u + p ( u ) c 2 Δ t τ 1 β u c 4 Δ t 3 τ 5 ,
    f 4 0 = 1 24 α u + p ( u ) c 2 Δ t τ 1 β u c 4 Δ t 3 τ 5 ,
where f 1 0 = f 2 0 , f 3 0 = f 4 0 .

3. Numerical Simulations

Different values of k and n determine the relationship between the GPLB and SLBGK models:
(I)
k = n = 1 , the SLBGK scheme;
(II)
n = k 2 , the LW scheme;
(III)
k 2 < n < k , here we choose n = k + k 2 / 2 ;
(IV)
n = k , the FP scheme;
(V)
k < n < 1 , here, we choose n = k + 0.1 .
Among them, the LW scheme (II) performs better than other schemes.
To illustrate the GPLB model constructed by the combination of Equations (6) and (10), the numerical simulations of Equation (1) are developed. With taking the LW scheme, the non-equilibrium extrapolation scheme is used to treat the boundary condition [22]. In each of the test examples, the exact solution determines the initial and boundary condition. We use E 2 (root-mean-square error), E (maximum absolute error) and GRE (global relative error) to measure accuracy of our model by comparing the LBM solutions with the exact solutions [23].
      E 2 = 1 M i = 1 M u x i , t u x i , t 2 ,
      E = max i = 1 , 2 , , M u x i , t u x i , t ,
      G R E = i = 1 M u x i , t u x i , t i = 1 M u x i , t .
In our model, u ( x i , t ) represents the LBM solution, u ( x i , t ) represents the exact solution.The number of lattices is M.
Example 1.
When α = 0 , β = 1 , and p ( u ) = u 2 , and Equation (1) is as follows,
u t t u 2 x x + u x x x x = 0 .
Ref. [24] gives the exact solution u ( x , t ) of Equation (51)
u ( x , t ) = 1 2 b 2 4 s 2 + 6 s 2 tanh 2 [ s ( x b t ) ] ,
where b and s are arbitrary constants.
We set x [ 10 , 10 ] , b = 0.2 , s = 0.3 , k = 0.9 , n = k 2 , τ = 4 , Δ x = 0.025 , and Δ t = 0.00025 . Figure 1 shows the evolution process of numerical solution and exact solution with time. It can be seen from Figure 1 that the soliton propagates along the negative direction of the x-axis. The space-time evolution graph of the LBM solution and exact solution are listed in Figure 2. Table 1 lists the error of the LBM solutions at different times. The LBM solutions agree with the exact solutions well.
Example 2.
When α = 1 , β = 1 , and p ( u ) = 3 u 2 , and Equation (1) is as follows,
u t t u x x 3 u 2 x x u x x x x = 0 .
Ref. [25] gives the exact solution u ( x , t ) of Equation (53).
u ( x , t ) = 2 a r 2 exp r x + r 1 + r 2 t 1 + a exp r x + r 1 + r 2 t 2 .
We set x [ 80 , 40 ] , a = 0.2 , r = 0.2 , k = 0.5 , n = k 2 , τ = 3.9 , Δ x = 0.25 , and Δ t = 0.00025 . Figure 3 shows the evolution process of the numerical solution and exact solution with time. It can be seen from Figure 3 that the solution propagates along the negative direction of the x-axis. The space-time evolution graph of the LBM solution and exact solution is shown in Figure 4. Table 2 lists the error of the LBM solutions at different times. The LBM solutions agree with the exact solutions well.
Example 3.
We set α = 1 , β = 1 , and p ( u ) = u 2 , Equation (1) becomes the good nonlinear Boussinesq equation, which is of the following form,
u t t = u x x + u 2 x x u x x x x .
The exact solution u ( x , t ) of Equation (55) is as follows,
u ( x , t ) = A s e c h 2 A 6 x + 1 2 A 3 t + x 0 .
In the computational domain of simulation, x [ 20 , 20 ] . A = 0.5 , x 0 = 0 , τ = 1.8 , Δ x = 0.1 , and Δ t = 0.0001 are setted. Figure 5 shows the evolution process of the numerical solution and exact solution with time. It can be seen from Figure 5 that the solution propagates along the negative direction of the x-axis. The space-time evolution graph of the LBM solution and exact solution are shown in Figure 6. Table 3 lists the error of the LBM solutions at different times. The LBM solutions agree with the exact solutions well.
Example 4.
We set α = 1 , β = 1 , and p ( u ) = u 2 , Equation (1) becomes the bad Boussinesq equation, which is of the following form,
u t t = u x x + u 2 x x + u x x x x .
The exact solution u ( x , t ) of Equation (57) is as follows,
u ( x , t ) = A s e c h 2 A 6 x + 1 + 2 A 3 t + x 0 .
In the computational domain of the simulation, x [ 30 , 30 ] . A = 0.08 , x 0 = 0 , τ = 16 , Δ x = 0.065 , and Δ t = 0.0001 are setted. Figure 7 shows the evolution process of numerical solution and exact solution with time. It can be seen from Figure 7 that the solution propagates along the negative direction of the x-axis. The space-time evolution graph of the LBM solution and exact solution are shown in Figure 8. Table 4 lists the error of the LBM solutions at different times. The LBM solutions agree with the exact solutions well.

4. Conclusions

In this paper, we have developed a general propagation lattice Boltzmann model for the generalized nonlinear Boussinesq equation. The macroscopic equation is recovered correctly from our model with the second-order accuracy through the Chapman–Enskog analysis. By applying the D1Q5 velocity model, Boussinesq equations are simulated. The numerical results agree well with the exact solutions with selecting the appropriate parameters that affect the propagation process. The results show that our model can remain stable and accurate, which is an effective algorithm worthy of promotion and application.

Author Contributions

W.Y. contributed to the conceptualization and the writing. C.L. helped reviewing the article. All authors have read and agreed to the published version of the manuscript.

Funding

This research received funding from the National Natural Science Foundation of China (Grant No.11761005).

Data Availability Statement

Not applicable.

Acknowledgments

We express our sincere thanks to the Editors, Referees.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Mohamad, A.A. Lattice Boltzmann Method; Electronic Industry Press: Beijing, China, 2015. [Google Scholar]
  2. Chen, S.Y.; Doolen, G.D. Lattice Boltzmann method for fluid flows. Annu. Rev. Fluid Mech. 1998, 30, 329–364. [Google Scholar] [CrossRef] [Green Version]
  3. Qian, Y.H.; Succi, S.; Orszag, S.A. Recent advances in lattice Boltzmann computing. Annu. Rev. Comput. Phys. 1995, 3, 195–242. [Google Scholar]
  4. Krüger, T.; Kusumaatmaja, H.; Kuzmin, A.; Shardt, O.; Silva, G.; Viggen, E.M. The Lattice Boltzmann Method: Principles and Practice; Springer: Cham, Switzerland, 2017. [Google Scholar]
  5. Inamuro, T.; Ogata, T.; Tajima, S.; Konishi, N. A lattice Boltzmann method for incompressible two-phase flow with large density differences. J. Comput. Phys. 2004, 198, 628–644. [Google Scholar] [CrossRef]
  6. Luo, L.S.; Girimaji, S.S. Theory of the lattice Boltzmann method: Two-fluid model for binary mixtures. Phys. Rev. E 2003, 67, 036302. [Google Scholar] [CrossRef] [PubMed]
  7. Guo, Z.L.; Zhao, T.S. Lattice Boltzmann model for incompressible flows through porous media. Phys. Rev. E 2002, 66, 036304. [Google Scholar] [CrossRef] [PubMed]
  8. Shi, B.C.; Guo, Z.L. Lattice Boltzmann model for nonlinear convection-diffusion equations. Phys. Rev. E 2009, 79, 016701. [Google Scholar] [CrossRef]
  9. Guo, X.Y.; Shi, B.C.; Chai, Z.H. General propagation lattice Boltzmann model for nonlinear advection-diffusion equations. Phys. Rev. E 2018, 97, 043310. [Google Scholar] [CrossRef]
  10. Chai, Z.H.; Shi, B.C.; Guo, Z.L. A multiple-relaxation-time lattice Boltzmann model for general nonlinear anisotropic convection-diffusion equations. J. Sci. Comput. 2016, 69, 355. [Google Scholar] [CrossRef]
  11. Chai, Z.H.; He, N.Z.; Guo, Z.L.; Shi, B.C. Lattice Boltzmann model for high-order nonlinear partial differential equations. Phys. Rev. E 2018, 97, 013304. [Google Scholar] [CrossRef] [Green Version]
  12. Lan, Z.Z.; Hu, W.Q.; Guo, B.L. General propagation lattice Boltzmann model for a variable-coefficient compound KdV-Burgers equation. Appl. Math. Model. 2019, 73, 695–714. [Google Scholar] [CrossRef]
  13. Hirota, R. Exact envelope-soliton solutions of a nonlinear wave. J. Math. Phys. 1973, 14, 805–809. [Google Scholar] [CrossRef]
  14. Hirota, R. Exact N-soliton solutions of the wave equation of long waves in shallow-water and in nonlinear lattices. J. Math. Phys. 1973, 14, 810–814. [Google Scholar] [CrossRef]
  15. Hajji, M.A.; Al-Khaled, K. Analytic studies and numerical simulations of the generalized Boussinesq equation. Appl. Math. Comput. 2007, 191, 320–333. [Google Scholar] [CrossRef]
  16. Hu, W.P.; Deng, Z.C. Multi-symplectic method for generalized Boussinesq equation. Appl. Math. Mech. 2008, 29, 927–932. [Google Scholar] [CrossRef]
  17. Liu, F.; Shi, W.P.; Wu, F.F. A lattice Boltzmann model for the generalized Boussinesq equation. Appl. Math. Comput. 2016, 274, 331–342. [Google Scholar]
  18. He, Y.B.; Dong, X.L.; Lin, X.Y. Numerical analysis and simulation of solutions to a class of Boussinesq systems with source terms. Appl. Math. Mech. 2018, 39, 961–978. [Google Scholar]
  19. Guo, Z.L.; Zheng, C.G.; Zhao, T.S. A lattice BGK scheme with general propagation. J. Sci. Comput. 2002, 16, 569. [Google Scholar] [CrossRef]
  20. Hu, W.Q.; Li, Z.H. Investigation on different discrete velocity quadrature rules in gas-kinetic unified algorithm solving Boltzmann model equation. Comput. Math. Appl. 2018, 75, 4179. [Google Scholar] [CrossRef]
  21. Chapman, S.; Cowling, T.G. The Mathematical Theory of Non-Uniform Gases, 3rd ed.; Cambridge University Press: Cambridge, UK, 1970. [Google Scholar]
  22. Guo, Z.L.; Zheng, C.G.; Shi, B.C. Non-equilibrium extrapolation method for velocity and pressure boundary conditions in the lattice Boltzmann method. Chin. Phys. 2002, 11, 366. [Google Scholar]
  23. Burden, R.L.; Faires, J.D. Numerical Analysis, 7th ed.; International Thomson Publishing: Belmont, CA, USA, 2001. [Google Scholar]
  24. Fu, Z.T.; Liu, S.K.; Liu, S.D. The JEFE method and periodic solutions of two kinds of nonlinear wave equations. Commun. Nonlinear Sci. Numer. Simul. 2003, 8, 67–75. [Google Scholar] [CrossRef]
  25. Wazwaz, A.M. Construction of soliton solutions and periodic solutions of the Boussinesq equation by the modified decomposition method. Chaos Solitons Fractals 2001, 12, 1549–1556. [Google Scholar] [CrossRef]
Figure 1. LBM solution (a) and exact solution (b) for the propagation of the soliton from t = 0 to t = 2 for Example 1.
Figure 1. LBM solution (a) and exact solution (b) for the propagation of the soliton from t = 0 to t = 2 for Example 1.
Entropy 24 00486 g001
Figure 2. LBM solutions versus exact solutions at (a) t = 0.5 ; (b) t = 1 ; (c) t = 1.5 ; (d) t = 2 for Example 1.
Figure 2. LBM solutions versus exact solutions at (a) t = 0.5 ; (b) t = 1 ; (c) t = 1.5 ; (d) t = 2 for Example 1.
Entropy 24 00486 g002
Figure 3. LBM solution (a) and exact solution (b) for the propagation of the solution from t = 0 to t = 40 for Example 2.
Figure 3. LBM solution (a) and exact solution (b) for the propagation of the solution from t = 0 to t = 40 for Example 2.
Entropy 24 00486 g003
Figure 4. LBM solutions versus exact solutions at (a) t = 10 ; (b) t = 20 ; (c) t = 30 ; (d) t = 40 for Example 2.
Figure 4. LBM solutions versus exact solutions at (a) t = 10 ; (b) t = 20 ; (c) t = 30 ; (d) t = 40 for Example 2.
Entropy 24 00486 g004
Figure 5. LBM solution (a) and exact solution (b) for the propagation of the solution from t = 0 to t = 2 for Example 3.
Figure 5. LBM solution (a) and exact solution (b) for the propagation of the solution from t = 0 to t = 2 for Example 3.
Entropy 24 00486 g005
Figure 6. LBM solutions versus exact solutions at (a) t = 0.5; (b) t = 1; (c) t = 1.5; (d) t = 2 for Example 3.
Figure 6. LBM solutions versus exact solutions at (a) t = 0.5; (b) t = 1; (c) t = 1.5; (d) t = 2 for Example 3.
Entropy 24 00486 g006
Figure 7. LBM solution (a) and exact solution (b) for the propagation of the solution from t = 0 to t = 4 for Example 4.
Figure 7. LBM solution (a) and exact solution (b) for the propagation of the solution from t = 0 to t = 4 for Example 4.
Entropy 24 00486 g007
Figure 8. LBM solutions versus exact solutions at (a) t = 1 ; (b) t = 2 ; (c) t = 3 ; (d) t = 4 for Example 4.
Figure 8. LBM solutions versus exact solutions at (a) t = 1 ; (b) t = 2 ; (c) t = 3 ; (d) t = 4 for Example 4.
Entropy 24 00486 g008
Table 1. Comparison of LBM solutions and exact solutions for Example 1 at different times.
Table 1. Comparison of LBM solutions and exact solutions for Example 1 at different times.
t = 0.5 t = 1 t = 1.5 t = 2
E 2 1.7818 × 10 4 7.1242 × 10 4 1.6018 × 10 3 2.8446 × 10 3
E 5.0682 × 10 4 2.0260 × 10 3 4.5535 × 10 3 8.0828 × 10 3
GRE 7.7264 × 10 4 3.0896 × 10 3 6.9474 × 10 3 1.2341 × 10 2
Table 2. Comparison of LBM solutions and exact solutions for Example 2 at different times.
Table 2. Comparison of LBM solutions and exact solutions for Example 2 at different times.
t = 10 t = 20 t = 30 t = 40
E 2 3.0848 × 10 6 6.1918 × 10 6 9.2945 × 10 6 1.2367 × 10 5
E 2.2720 × 10 4 4.3032 e × 10 4 6.2525 × 10 4 8.2072 × 10 4
GRE 1.2330 × 10 2 2.5070 × 10 2 3.7928 × 10 2 5.1775 × 10 2
Table 3. Comparison of LBM solutions and exact solutions for Example 3 at different times.
Table 3. Comparison of LBM solutions and exact solutions for Example 3 at different times.
t = 0.5 t = 1 t = 1.5 t = 2
E 2 1.3459 × 10 5 4.9778 × 10 5 1.0154 × 10 4 1.6501 × 10 4
E 2.1978 × 10 3 2.8161 × 10 3 5.3003 × 10 3 7.5997 × 10 3
GRE 2.6839 × 10 3 8.2844 × 10 3 1.6677 × 10 2 2.7141 × 10 2
Table 4. Comparison of LBM solutions and exact solutions for Example 4 at different times.
Table 4. Comparison of LBM solutions and exact solutions for Example 4 at different times.
t = 1 t = 2 t = 3 t = 4
E 2 8.4609 × 10 7 3.2913 × 10 6 7.2237 × 10 6 1.2485 × 10 5
E 4.8385 × 10 4 1.8602 × 10 4 4.0126 × 10 4 6.9719 × 10 4
GRE 9.3621 × 10 4 3.6597 × 10 3 8.0689 × 10 3 1.4017 × 10 2
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Yang, W.; Li, C. General Propagation Lattice Boltzmann Model for the Boussinesq Equation. Entropy 2022, 24, 486. https://doi.org/10.3390/e24040486

AMA Style

Yang W, Li C. General Propagation Lattice Boltzmann Model for the Boussinesq Equation. Entropy. 2022; 24(4):486. https://doi.org/10.3390/e24040486

Chicago/Turabian Style

Yang, Wei, and Chunguang Li. 2022. "General Propagation Lattice Boltzmann Model for the Boussinesq Equation" Entropy 24, no. 4: 486. https://doi.org/10.3390/e24040486

APA Style

Yang, W., & Li, C. (2022). General Propagation Lattice Boltzmann Model for the Boussinesq Equation. Entropy, 24(4), 486. https://doi.org/10.3390/e24040486

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop