Research Article | Open Access
Jie Zhao, Hong Li, Zhichao Fang, Xue Bai, "Numerical Solution of Burgers’ Equation Based on Mixed Finite Volume Element Methods", Discrete Dynamics in Nature and Society, vol. 2020, Article ID 6321209, 13 pages, 2020. https://doi.org/10.1155/2020/6321209
Numerical Solution of Burgers’ Equation Based on Mixed Finite Volume Element Methods
In this article, mixed finite volume element (MFVE) methods are proposed for solving the numerical solution of Burgers’ equation. By introducing a transfer operator, semidiscrete and fully discrete MFVE schemes are constructed. The existence, uniqueness, and stability analyses for semidiscrete and fully discrete MFVE schemes are given in detail. The optimal a priori error estimates for the unknown and auxiliary variables in the norm are derived by using the stability results. Finally, numerical results are given to verify the feasibility and effectiveness.
In this article, we consider the following one-dimensional Burgers’ equation:with initial and boundary conditionswhere , with , α is a positive constant, is the viscosity coefficient, and is the given initial function.
Burgers’ equation is a famous nonlinear evolution equation which was first derived by Bateman  in 1915. Burgers  utilized this equation to model turbulence behavior. Since its appearance, this equation has been widely concerned by researchers because of its various practical applications, such as gas dynamics, shock theory, traffic flows, viscous flow, and turbulence. The exact solutions can be expressed as a Fourier series expansion by introducing a Hopf–Cole transformation [3, 4]. Benton and Platzman  gave the exact solutions of Burgers’ equation in one-dimensional spatial regions for some different initial functions. In the past several decades, many numerical techniques had been constructed and tested to solve Burgers’ equation [6–13], such as finite element methods, finite difference methods, least-squares finite element methods, and spectral methods.
In recent years, the mixed finite element (MFE) methods and finite volume (FV) methods have been used to solve Burgers’ equation by many researchers. Luo and Liu  proposed an MFE method to solve one-dimensional Burgers’ equation by introducing a flux function as an auxiliary variable and gave the existence, uniqueness, and error analyses for the discrete solutions. Chen and Jiang  constructed a characteristic MFE scheme to solve one-dimensional Burgers’ equation and obtained the optimal error estimates for the velocity and flux (gradient) in the norm. Pany et al.  applied an -Galerkin MFE method to approximate the velocity and flux of one-dimensional Burgers’ equation and gave a priori error estimates and numerical experiments. Shi et al.  provided a low-order least-squares nonconforming characteristic MFE scheme to solve two-dimensional Burgers’ equation and obtained the optimal order error estimates in the norm. Hu et al.  constructed a Crank–Nicolson time discretization MFE scheme to treat two-dimensional Burgers’ equation by using the pair and gave the optimal error analysis and numerical experiments. Nascimento et al.  applied a Fourier pseudospectral method and an FV method to solve one-dimensional Burgers’ equation and gave the numerical comparison. Guo et al.  proposed a fifth-order FV weighted compact scheme to solve one-dimensional Burgers’ equation and gave numerical experiments. Sheng and Zhang  proposed an FV method to solve two-dimensional Burgers’ equation and obtained the optimal error estimate in the norm.
The aim of this article is to develop mixed finite volume element (MFVE) methods to solve one-dimensional Burger’s equation by combining the MFE methods [22–25] with the finite volume element (FVE) methods [26–30]. The MFVE methods, also called mixed covolume methods, were first proposed by Russell  to solve the elliptic equation. Now, the methods have been applied to solve second-order elliptic equations [32–34], integrodifferential equations , parabolic equations [36, 37], time-fractional partial differential equations , and so on. In this article, we introduce a flux function as an auxiliary variable, rewrite (1) as the first-order system, and construct the semidiscrete and nonlinear backward Euler fully discrete MFVE schemes by using a transfer operator. We give the theoretical analysis for semidiscrete and fully discrete schemes in detail, including existence, uniqueness, and stability. In particular, in the analysis of the fully discrete scheme, we apply the Brouwer fixed-point theorem to prove the existence and use the Sobolev embedding theorem and the inverse inequality to prove the stability. Making use of the stability results of discrete solutions, we obtain the optimal a priori error estimates for the velocity and flux in the norm. Furthermore, we give some numerical results to verify the feasibility and effectiveness of the MFVE scheme.
The rest of this article is organized as follows: In Section 2, we use a flux function as an auxiliary variable and give the mixed variational formulation and the semidiscrete MFVE scheme. The existence, uniqueness, stability, and convergence analyses for semidiscrete and fully discrete schemes are given in Section 3 and 4, respectively. In Section 5, a numerical example is given to verify the theoretical results. In this article, the standard definitions and notations of the Sobolev spaces as in  are used. Furthermore, we use the symbol C to represent a generic constant which is independent of the space and time mesh parameters h and .
2. Semidiscrete MFVE Scheme
The mixed variational formulation of (3) is to find such that
Now, we construct the primal mesh for the interval with the nodes , where N is some positive integer. Then, the primal mesh is denoted by , and the diameter of the primal mesh is defined by , where . We assume that the mesh satisfies the quasiuniform condition for some constant .
Next, the corresponding dual mesh is constructed by the nodes , where . Denote , , and , then the dual mesh is defined by .
We choose the mixed finite element space as the trial function space, where
Let and be the basis of the spaces and , respectively, where is the piecewise linear polynomial defined in  and is the characteristic function of the set .
The system (3) is integrated as follows:
Now, we define the transfer operator (see ) as follows:
Thus, we have . The range of is used as the test function space. Then, we can rewrite (6) as
Note thatBy a simple calculation, it is easy to have that , , . Then, we get the semidiscrete MFVE scheme to find such thatwhere satisfies
3. Theoretical Analysis for the Semidiscrete MFVE Scheme
3.1. Some Lemmas
For theoretical analysis, we first give some properties of the transfer operator and two projection operators.
Lemma 1 (see ). The transfer operator is bounded, that is,
Lemma 2 (see ). The transfer operator satisfies the following symmetry relation:
Lemma 3 (see ). The transfer operator satisfies the following positivity:
Lemma 4 (see ). Let I be an identity operator, then the transfer operator satisfies the following properties:
Lemma 5 (see ). There exists a constant such that
Now, the elliptic projection operator is introduced below, which satisfies
At the same time, the orthogonal projection operator is introduced to satisfy
Lemma 6. There exists a constant such that, for ,
3.2. Existence, Uniqueness, and Stability Analyses
Theorem 1. There exists a unique discrete solution for the semidiscrete MFVE scheme (10).
Proof. Obviously, there exists a unique solution for the scheme (11). Let and be the basis functions of the spaces and , respectively, then and can be expressed as follows:Choosing and in (10), we rewrite the semidiscrete scheme (10) in the matrix form to find such thatwhereIt is easy to know that matrices A and D are symmetrically positive definite. Then, the system (21) can be rewritten asAccording to the theory of differential equation, we can see that the system (23) has a unique solution, which shows that there exists a unique discrete solution for the semidiscrete MFVE scheme (10).
Theorem 2. Let be the discrete solution of the semidiscrete scheme (10), then there exists a constant such that
Proof. Choosing and in (10), we haveNoting that , we rewrite (25) asIntegrating (26) from 0 to t, we getApplying Lemma 1 and Lemma 3, we obtainCalculating the derivative of (10) with respect to t, we getSetting in (10) and in (29), we haveApplying the Sobolev embedding theorem to estimate the right-hand side of (30), we obtainMaking use of (31) and Lemma 3 in (30), we haveIntegrating (32) from 0 to t, we obtainApply the Gronwall lemma and (28) in (33), we obtainNext, choosing in (10), we getNoting that , we haveSubstituting (34) into the above inequality, we haveThus, applying the Sobolev embedding theorem, we obtainNow, we estimate and . We choose in (11) to obtainChoosing in (11), we haveMaking use of (39), we easily getThus, we obtain the proof of Theorem 2.
3.3. A Priori Error Estimates
Theorem 3. Let and be the solutions of the systems (4) and (10), respectively. Assume that the initial solution satisfies (11), then there exists a constant such thatwhere represents the function of and .
Proof. Taking and in (43), we haveApplying Theorem 2, we have . We apply the Cauchy–Schwarz inequality and the Young inequality in (46) to obtainIntegrating (47) from 0 to t, and applying Lemma 3, we obtainApplying the Gronwall lemma in (48) yieldsNext, calculating the derivative of (43) with respect to t, we obtainChoosing in (50) and in (43), we haveApplying the stability results in Theorem 2, we havewhere represents the function of and . Applying Lemmas 3–5 in (51), and making use of (52), we obtainIntegrating (53) from 0 to t, we getApplying the Gronwall lemma in (54), we obtainWe estimate by choosing in (44) and noting that :Applying Lemma 6 in (56), we haveFinally, we apply Lemma 6 and the triangle inequality to obtain the desired conclusion.
4. Fully Discrete MFVE Scheme and Its Theoretical Analysis
4.1. Fully Discrete MFVE Scheme
Let be an equidistant partition of the time interval , where M is a positive integer. Denote to represent the step length and . And denote and for a function ψ.
Let and be the fully discrete solutions of u and p at , respectively. We can obtain the following nonlinear backward Euler MFVE scheme to find :where the initial value satisfies the following equations:
Remark 1. The fully discrete MFVE scheme (58) is implicit in time. In the actual calculation of , we need to make predictions first by using the linear backward Euler MFVE scheme defined as follows:Similar to the proof process of Theorem 4 in Reference , we can also obtain that there exists a unique discrete solution for the linear MFVE scheme (60).
4.2. Existence, Uniqueness, and Stability Analyses
Lemma 7. Let be a finite-dimensional inner product space with a norm . Let be a continuous operator. Assume that there exists such that for with . Then, there exists such that and .
Theorem 4. Suppose that has been given, then there exists a unique fully discrete solution for the nonlinear backward Euler MFVE scheme (58).
Proof. Choosing in (58), we haveThe operator is defined as follows:It is easy to know that the operator is continuous. Setting in (62), we haveSelecting to satisfy , we have . Making use of Lemma 7, we know that there exists such that . Thus, we choose to satisfy (61). Furthermore, selecting in (58), we getThus, there exists a solution which satisfies (64). Then, it is obviously known that satisfies the scheme (58), which proves the existence of the fully discrete solutions.
Next, we give the uniqueness for the fully discrete scheme (58). Let be another solution of the scheme (58) with the initial value , then we haveLet , then we haveThe following proof is based on the mathematical induction. First, ; next assume and then choose in (66) to obtainApplying the stability results in Theorem 5 (because the existence of the discrete solutions has been proved), we haveMultiplying (68) by , and making use of Lemma 3, we obtainSelecting to satisfy , we have and , which indicate that , and then is known from (58). Thus, we complete the proof of Theorem 4.
Proof. Similarly to the proof process of Theorem 2, we can obtain the estimate of the initial solution . We choose and in (58) to obtainTaking note of and in (72), we obtainMultiplying (73) by , summing from 1 to n, and applying Lemma 1 and Lemma 3, we haveNext, we make use of (58) and (59) to obtainChoosing in (58) and in (75), we get