Research Article | Open Access
Yueyue Pan, Lifei Wu, Xiaozhong Yang, "A New Class of Difference Methods with Intrinsic Parallelism for Burgers–Fisher Equation", Mathematical Problems in Engineering, vol. 2020, Article ID 9162563, 17 pages, 2020. https://doi.org/10.1155/2020/9162563
A New Class of Difference Methods with Intrinsic Parallelism for Burgers–Fisher Equation
This paper proposes a new class of difference methods with intrinsic parallelism for solving the Burgers–Fisher equation. A new class of parallel difference schemes of pure alternating segment explicit-implicit (PASE-I) and pure alternating segment implicit-explicit (PASI-E) are constructed by taking simple classical explicit and implicit schemes, combined with the alternating segment technique. The existence, uniqueness, linear absolute stability, and convergence for the solutions of PASE-I and PASI-E schemes are well illustrated. Both theoretical analysis and numerical experiments show that PASE-I and PASI-E schemes are linearly absolute stable, with 2-order time accuracy and 2-order spatial accuracy. Compared with the implicit scheme and the Crank–Nicolson (C-N) scheme, the computational efficiency of the PASE-I (PASI-E) scheme is greatly improved. The PASE-I and PASI-E schemes have obvious parallel computing properties, which show that the difference methods with intrinsic parallelism in this paper are feasible to solve the Burgers–Fisher equation.
The Burgers–Fisher equation is a typical model for describing diffusion propagation and convection conduction. It is an important partial differential equation in mathematical physics and widely applied not only in the study of gas dynamics and heat conduction, but also in explaining many physical phenomena such as elasticity [1–3]. In recent years, the numerical method for solving the Burgers–Fisher equation has been extensively concerned by researchers [4–6]. They study this equation both for conceptual understanding of physical flows and testing various numerical methods. So, the fast method of solving it has basic scientific significance and application value.
Up to now, researchers have given many numerical methods for solving the Burgers–Fisher equation. Ismail et al.  obtained the approximate solution of the Burgers–Fisher equation by using the Adomian decomposition method. Although the approximate solution had high accuracy, the method neither allows to choose different basis functions, nor could it adjust the convergence region and speed of the asymptotic solution. Babolian and Saeidian  applied the homotopy analysis method to give the analytical approximate solution of the Burgers–Fisher equation. This method is effective for solving such partial differential equations. However, auxiliary operators, auxiliary parameters, and auxiliary functions have been introduced in a wide range, and there is no strict mathematical theory to guide their selection to ensure the convergence of numerical solution series. Chandraker et al.  proposed a semi-implicit difference scheme and a modified C-N difference method to solve the Burgers–Fisher equation. For this method, numerical solution errors are much less, but the stability of the semi-implicit method is weak, and it is conditionally stable. Yadav and Jiwari  presented finite element analysis and approximation of the Burgers–Fisher equation, which proves that the scheme is convergent. However, the calculated amount of the method is very large, and its parallel calculation is not as intuitive as the finite difference method. And, the computational efficiencies, especially the computation time, of these methods in the above documents are low. For the high-dimensional Burgers–Fisher equation, even with high-performance computers, it is difficult to simulate a large computational domain.
With the rapid development of multicore and cluster technology, the parallel algorithm has become one of the mainstream technologies to improve the computing efficiency. It is well known that classical explicit difference schemes have ideal parallelism and are suitable for parallel computation, but they are conditionally stable. Especially in multidimensional problems, the time step of computing is severely restricted. The classical implicit difference scheme and the C-N difference scheme are absolutely stable, but they are not suitable for direct and effective application on parallel computers. For parabolic equations, Evans and Abdullah  proposed the idea of group explicit and designed an alternating group explicit (AGE) difference scheme, which not only ensures the stability of numerical calculations, but also has a good parallel property. Inspired by the construction of the AGE method, Zhang et al. [12, 13] proposed the idea of using the Saul’yev asymmetric scheme to construct a piecewise implicit scheme and established a class of alternating segment explicit-implicit (ASE-I) parallel difference methods and alternating segment C-N (ASC-N) parallel difference methods, where the stability and parallelism were both obtained. Han et al.  constructed pure alternating segment explicit-implicit (PASE-I) and pure alternating segment implicit-explicit (PASI-E) parallel difference schemes for the constant coefficient diffusion equation. The numerical solution had good stability and convergence. Academician Zhou  referred the explicit-implicit mixing scheme of the most general parabolic equation as the difference scheme with intrinsic parallelism. He studied the existence, uniqueness, convergence, and stability of differential decomposition and established the basic theory of the difference method with intrinsic parallelism for parabolic equations. Until now, the difference method with intrinsic parallelism has been extended to solve many evolution equations [16–18]. Qu and Wang  constructed an alternating segment explicit-implicit difference scheme for solving the KdV equation. The linear absolute stability of the scheme has also been demonstrated. Yuan et al.  proposed a class of parallel difference schemes with 2-order spatial accuracy and unconditional stability for the nonlinear parabolic system. Guo et al.  studied the difference method with intrinsic parallelism for the dispersive equation. The general alternating difference schemes with variable time steps are constructed and proved to be unconditionally stable. Namjoo et al.  presented the numerical solution of the generalized Burgers–Fisher equation on the basis of the nonstandard finite-difference (NSFD) scheme. The positivity, consistency, and boundedness of the scheme are discussed. The numerical results obtained by the NSFD scheme is compared with the exact solution and some available methods to verify the accuracy and efficiency of the NSFD scheme. For the nonlinear Leland equation, it is very difficult to obtain an analytical solution. Therefore, let the numerical solution of the C-N scheme approximately substitute the exact solution. Yan et al.  constructed PASE-I and PASI-E difference schemes, which combine the classical explicit scheme with the classical implicit scheme. The numerical experiments verify that the calculation accuracy of PASE-I and PASI-E schemes is better than that of the existing alternating segment C-N scheme and alternating segment explicit-implicit and implicit-explicit schemes. The main advantage of these difference schemes with intrinsic parallelism is that these methods can be directly applied to parallel computer systems with distributed memory and minimize the communication between processors. The algorithm only needs local message passing between adjacent processors, and the communication and computation involved are also local. It is easy to balance the load of computation and communication, so as to obtain good precision and extendibility of parallel computing.
For a long time, a large number of parallel difference schemes constructed are either conditionally stable or unconditionally stable, but the space has only 1-order accuracy . To obtain a parallel difference scheme with higher precision and more relaxed stability conditions, we study the parallel algorithm of the Burgers–Fisher equation from the point of view of the parallelization of traditional difference schemes. We hope to skip the difficulties of the numerical algebra and open up another way of parallelization. In this paper, a new class of difference method with intrinsic parallelism for the Burgers–Fisher equation is proposed by making piecewise processing of alternating whole explicit and implicit difference schemes. A new class of PASE-I and PASI-E parallel difference schemes are constructed by taking simple classical explicit and implicit schemes. Both numerical experiments and theoretical analysis show that PASE-I and PASI-E schemes have obvious parallel computational properties. The schemes show linear absolute stability and convergence with 2-order time accuracy and 2-order spatial accuracy. Numerical examples give that the computational efficiencies of PASE-I and PASI-E difference schemes are much higher than those of the implicit difference scheme. And, the parallel difference methods proposed in this paper are efficient and feasible for solving the Burgers–Fisher equation.
2. The Intrinsic Parallel Difference Schemes of Burgers–Fisher Equation
2.1. Burgers–Fisher Equation
The initial condition is
The boundary conditions are
Taking the determined , , and , respectively, as
We can get the analytic solution:
2.2. The Construction of PASE-I Parallel Difference Scheme
The solution region is divided into grids. Respectively, the space and time steps are and , where is a positive integer. Here, and . is the numerical solution, and is the analytic solution. To construct the PASE-I parallel difference scheme, the classical explicit difference scheme, the classical implicit difference scheme, and the C-N difference scheme of equation (1) will be given first.
At first, the classical explicit difference scheme is
Secondly, the classical implicit difference scheme is
Thirdly, the C-N difference scheme iswhere , , , , and . The discrete schemes of are in the explicit scheme and in the implicit scheme.
The classical explicit scheme has ideal parallelism, but it is conditionally stable. The classical implicit scheme is absolutely stable, but the inverse matrix of the tridiagonal matrix needs to be solved, so it is not convenient for us to get the results directly and quickly. So, combining the classical explicit difference scheme (6) with the classical implicit difference scheme (7), the specific approach is given as follows.
Let , where and are positive integers, is odd, and and . The points calculated at the same time layer are divided into segments, which are recorded in the order . Every segment of the odd time layer is arranged from left to right in the order of “classical explicit-classical implicit-classical explicit.” On the even time layer, the order of arrangement becomes “classical implicit-classical explicit-classical implicit.” The solution of each implicit segment relies on the calculation of the first or last point of the adjacent explicit segment to give its internal boundary value. See Figure 1 for the detail, where the classical explicit scheme is applied in the place and the classical implicit scheme is applied in the place.
For realizing the parallel computing of the PASE-I scheme, for , we consider the calculation of the explicit segment point . They are calculated with the classical explicit scheme (6), leading to the following explicit segment:
The implicit segment is
A complete calculation step of the PASE-I scheme is as follows: is the number of time layers, is the number of segments, and is the point on segment :(1)for (2) if mod (3) for (4) if mod (5) solve equation (9) to get ;(6) else(7) solve equation (10) to get ;(8) end(9) end(10) else(11) for (12) if mod (13) solve equation (9) to get ;(14) else(15) solve equation (10) to get ;(16) end(17) end(18) end(19)end for
The matrix form of the PASE-I scheme iswherein which , , and are the dimensional vectors. ; ; is the th-order zero matrix:
As can be seen from the piecewise pattern in Figure 1 and the expression of the matrix , the PASE-I parallel difference scheme turns a discrete problem of orders into some small independent problems to be solved. The computation at each time level is only to solve a low-order equations system at each implicit segment. Compared with the classical implicit scheme, there are fewer equation systems in the odd layer and fewer equation systems in the even layer. The computation is simple, and the parallel characteristic is obvious.
3. The Numerical Analysis of PASE-I Parallel Difference Method
3.1. Existence and Uniqueness of Solution to PASE-I Scheme
In order to discuss the existence and the uniqueness of the PASE-I scheme solution, we need to introduce the following two lemmas.
Lemma 1. (Kellogg lemma, see ). Let , and let be a nonnegative (or positive) definite matrix. Then, exists, and
Lemma 2. and in the PASE-I scheme (11) for solving the Burgers–Fisher equation are nonnegative definite matrices.
Proof. We only need to prove and are nonnegative definite matrices. That is to say, is a nonnegative definite matrix. The coefficient in formula (11) is assumed to be constant , calculate as follows:We consider the case of low-speed flow and might as well set . With , , , , and , is a diagonally dominant matrix, and the diagonal elements of are nonnegative real numbers. Therefore, is a nonnegative definite matrix. Similarly, is also a nonnegative definite matrix. Therefore, and are nonnegative definite matrices.
From the initial conditions and the boundary conditions of the Burgers–Fisher equation, we know the value of . Assuming that the value of the (2n)th time layer is known, the value of the (2n + 1)th time layer waits for calculating. From the PASE-I scheme (11), the matrix equation for calculating the value of the (2n + 1)th time layer isApparently, the right of equation (16) is known, and exists by Lemmas 1 and 2. Then, equation (16) has a unique solution.
In the same way, applying the PASE-I scheme to calculate the value of the (2n + 2)th time layer, the matrix equation isWe could also prove that the matrix equation (17) has a unique solution. Then, we could get the following.
Theorem 1. The solution of the PASE-I parallel difference scheme (11) for the Burgers–Fisher equation is existing and unique.
3.2. Linear Absolute Stability of PASE-I Scheme
Lemma 3. (see ). Let , and let be a nonnegative definite matrix. Then,
For any even number , there is
We can get the following inequality easily:where .
Therefore,where . Then, there is
Therefore, we have the following theorem.
Theorem 2. The PASE-I parallel difference scheme (11) for the Burgers–Fisher equation is linearly absolute stable.
3.3. Computational Accuracy of PASE-I Scheme
The classical explicit scheme is
The classical implicit scheme is
Let each point of the above two schemes be expanded as the Taylor series at the point . The truncation errors are, respectively, and . Then, we getwhere .
With the equation , we can get
In the expressions of and , the terms with the same absolute value but the opposite sign are contained, respectively. Piecewise explicit and piecewise implicit alternate not only at the same time level but also at different time levels. Some error terms will cancel each other. Therefore, the computational accuracy of the PASE-I scheme is .
Theorem 3. The computational accuracy of the PASE-I parallel difference scheme (11) for the Burgers–Fisher equation is .
3.4. Convergence of PASE-I Scheme
Assume that is the nominal solution of problem (30). An application of the quasilinearization process to the nonlinear problem (30) introduces a sequence of linear equations determined by the following recurrence relation:where is the iteration index. This is an application of the Newton–Raphson–Kantorovich approximation method in function space. We choose a reasonable initial guess satisfying the initial condition . For the sake of convenience, we let . Therefore, the above equation leads to the following initial boundary value problem:where
Further, we assume that the functions , , and are sufficiently smooth functions in the spatial direction with
By using quasilinearization process, we get the linear boundary value problem (32) for the function . For the solution of original nonlinear problem (30), we require thatwhereas numerically, we require thatwhere is the small prescribed value to terminate the computation. This is the requisite criterion for terminating the iteration, and the solution is used as the numerical solution of the nonlinear boundary value problem (30).
Let be the sequence of PASE-I scheme solutions at time level. To prove the convergence of the PASE-I scheme, we consider the following equation:where .
We assume to be the initial value. By using quasilinearization process and the PASE-I scheme, we obtain a sequence of linear equations determined by the following recurrence relation:
Thus, we have
Converting it into an integral function by using Green’s function, we havewhere Green’s function is defined by
The mean value theorem gives uswhere . Now, putting the value of in equation (40), we get the following estimate:
So, we obtain
Taking the maximum norm over the spatial domain and after some simplification, we have
Theorem 4. Let be the sequence of PASE-I scheme solutions at time level. Then, there exists a constant , independent of , such thati.e., the PASE-I scheme (11) for the Burgers–Fisher equation is convergent.
4. The PASI-E Parallel Difference Scheme for Burgers–Fisher Equation
Similar to the process of constructing the PASE-I parallel difference scheme, we can also construct the PASI-E parallel difference scheme of the Burgers–Fisher equation.
When solving the odd layer, we can calculate according to the rule of “classical implicit-classical explicit-classical implicit.” On the even time layer, the rule of calculation becomes “classical explicit-classical implicit-classical explicit.” We then get the following PASI-E parallel difference scheme for the Burgers–Fisher equation:with .
Imitating the numerical analysis process of the PASE-I parallel difference scheme, the following theorem is obtained.
Theorem 5. The pure alternating segment implicit-explicit (PASI-E) parallel difference scheme (48) for the Burgers–Fisher equation is uniquely solvable, linearly absolute stable, and convergent, besides that the order of convergence is .
5. Numerical Experiments
Numerical experiments were done in MATLAB R2014a, based on the Intel Core i5-4200 CPU@2.50 GHz.
Example 1. Let , , , and . We consider the Burgers–Fisher equation as follows [28, 29]:with the initial and boundary conditionsand the analytic solutionBy taking time layer , space layer , piecewise number , interior point number , and parameter , the analytic solution surface, implicit scheme solution surface, PASE-I scheme solution surface, and PASI-E scheme solution surface are given. Figures 2–5 show that these surfaces are of the same shape as a whole.
We take time layer , space layer , and parameter . At , the numerical solutions of the PASE-I scheme and the PASI-E scheme in this paper are compared with the analytic solution and the implicit scheme solution. The computing results are listed in Table 1. It is clearly shown that the numerical solutions of the PASE-I scheme and the PASI-E scheme are very close to the analytic solution. The value of u is between 0 and 1, which is appropriate and is consistent with the theoretical analysis.
Figure 6 is the comparison of the PASE-I scheme solution, PASI-E scheme solution, and analytic solution. Obviously, the solution curves of these two schemes are very close to the analytic solution. Therefore, the PASE-I and PASI-E parallel difference schemes are high-accuracy difference schemes for solving the Burgers–Fisher equation.
We treat the analytic solution as the control solution and the solution of the PASE-I (PASI-E) scheme as the perturbation solution. The definition of the node error (NE) is . The node error distributions of the implicit scheme solution, PASE-I scheme solution, and PASI-E scheme solution relative to the analytic solution are given below by taking time layer , space layer , and parameter .
Figures 7–9 are, respectively, the node error distributions of the implicit scheme solution, PASE-I scheme solution, and PASI-E scheme solution relative to the analytic solution. As can be seen from these figures, the computing errors of the PASE-I scheme and the PASI-E scheme are similar to that of the implicit scheme, and the maximum error is less than 6e − 04.
To verify the stability of the PASE-I scheme and the PASI-E scheme, the relative errors over time are given below. The definition of the sum of the relative error for every time level isIt can be seen from Figure 10 that when we take the parameter , the time layer , and the space layer , the SRET of the two scheme solutions remains stable and bounded with the increase in the time step. It also shows that the PASE-I and PASI-E methods of the Burgers–Fisher equation are computational stable.
Numerical experiments on the space and time convergence orders of the C-N scheme and the PASE-I scheme are given below. We define as the error, as the space-convergent order, and as the time-convergent order [30, 31]:At first, we verify the spatial accuracies of the C-N scheme and the PASE-I scheme by taking and . Table 2 shows that the spatial accuracies of the C-N scheme and the PASE-I scheme in this paper are approximately 2.
Secondly, we calculate the time accuracies of the C-N scheme and the PASE-I scheme by taking . We set the space step , that is, .
Table 3 shows that the time accuracies of the C-N scheme and the PASE-I scheme are approximately 2. Therefore, the PASE-I scheme constructed in this paper and the existing C-N scheme possess the same order of spatial accuracy and time accuracy. The numerical results agree well with theoretical analysis.
In order to better verify the accuracies of the PASE-I method and the PASI-E method, the distribution of errors in spatial lattices is also investigated. We define difference total energy asAs can be seen from Figure 11, when the time layer and the space layer , the DTE of the PASE-I and PASI-E scheme solutions is between 0 and 1.2e − 05. That is, the solutions of PASE-I and PASI-E schemes are very close to the analytic solution. Figure 11 shows that the PASE-I and PASI-E methods of the Burgers–Fisher equation have good accuracy.
The PASE-I scheme and the PASI-E scheme have superiority in the aspect of the computing time. But, this superiority can only be reflected when the amount of data is large. In case the number of grid points is less than a certain range, the effect of data communication on the cycle can reduce the computational efficiency, and the superiority of parallel computing is not obvious [13, 32]. When the amount of data is large, the impact of program loop execution is much greater than that of data communication, and the parallel scheme is more effective.
Figure 12 shows that the computing time of three schemes increases with the increase in the number of space grid points. When the number of space grid points is greater than 1000, the impact of program loop execution is much greater than that of data communication. The PASE-I scheme and the PASI-E scheme have obvious superiority with regard to the computing time.
To compare the computational efficiency of the implicit scheme, PASE-I scheme, and PASI-E scheme, the analysis of the scheme speedup ratio is done . The speedup ratios of the PASE-I and PASI-E schemes relative to the implicit scheme are, respectively, and . We fix the time layer as and take the number of space grids as . The computational results are shown in Table 4.
As we can see from Table 4, the computing time of three difference schemes is increasing with the increase in the number of space grids for the Burgers–Fisher equation. and are both greater than 7.6 and gradually increasing. It is shown that the PASE-I scheme and the PASI-E scheme have obvious advantages in saving the computing time. Therefore, the PASE-I scheme and the PASI-E scheme for the Burgers–Fisher equation are more efficient than the serial implicit scheme.
Example 2. Let , , and . Consider the following Burgers–Fisher equation:The initial and boundary conditions are determined by the analytic solution, and the analytic solution isFigures 13 and 14 are, respectively, the node error distributions of the C-N scheme solution and the PASE-I scheme solution with respect to the analytic solution with the time layer , the space layer , and the parameter . As can be seen from Figures 13 and 14, the maximum errors will not exceed 3e − 03, and the parameter has little effect on computational accuracy of the PASE-I scheme. It is shown that the PASE-I difference method with intrinsic parallelism in this paper is a high-precision difference method for solving the Burgers–Fisher equation.
The comparison and analysis of the computational efficiency of the two difference schemes are discussed as the spatial grid is continuously encrypted. For , Figure 15 shows the comparison of computing time between the PASE-I scheme and the C-N scheme. We can see that the PASE-I scheme requires the computing time as the same as the C-N scheme when . But, with the increase in the number of space grids, the computing time of the C-N scheme increases exponentially, and the computing time of the PASE-I scheme increases linearly. Compared with the C-N scheme, the advantage of the high efficiency of the PASE-I scheme is more and more obvious with the increase in the space grids.
The speedup , and efficiency ( is the computing time of C-N and is the computing time of parallel scheme) are defined [32, 33]. Using four cores for this numerical experiment, we fix the time layer as and take the number of space grid points as