Numerical Solution of In-Viscid Burger Equation in the Application of Physical Phenomena: The Comparison between Three Numerical Methods

In this paper, upwind approach, Lax–Friedrichs, and Lax–Wendroﬀ schemes are applied for working solution of In-thick Burger equation in the application of physical phenomena and comparing their error norms. First, the given solution sphere is discretized by using an invariant discretization grid point. Next, by using Taylor series expansion, we gain discretized nonlinear diﬀerence scheme of given model problem. By rearranging this scheme, we gain three proposed schemes. To verify validity and applicability of proposed techniques, one model illustration with subordinated to three diﬀerent original conditions that satisfy entropy condition are considered, and solved it at each speciﬁc interior grid points of solution interval, by applying all of the techniques. The stability and convergent analysis of present three techniques are also worked by supporting both theoretical and numerical ﬁne statements. The accuracy of present techniques has been measured in the sense of average absolute error, root mean square error, and maximum absolute error norms. Comparisons of numerical gets crimes attained by these three methods are presented in table. Physical behaviors of numerical results are also presented in terms of graphs. As we can see from numerical results given in both tables and graphs, the approximate solution is good agreement with exact solutions. Therefore, the present systems approaches are relatively eﬀective and virtually well suited to approximate the solution of in-viscous Burger equation.


Introduction
Burgers equation is considered one of the most common nonlinear PDEs [1] and it arises in multitudinous physical problems including one-dimensional turbulence, sound swells and shock swells during a thick medium, swells in fluid-filled thick elastic tubes, and magneto hydrodynamic swells in a medium with finite electrical conductivity [2][3][4][5][6] and fluid dynamics [35]. Burger equations were first introduced by Bateman in 1915 [8][9][10][11][12] and he proposed the steady-state result of the matter [5]. It also arises in numerous physical problems, including shock overflows, business overflows [1], nonlinear surge propagation [1,13], swells of impact and sound swells in a thick medium [1], gas dynamics, hydrodynamics [14,15], heat conduction [14], isentropic inflow [16], solid-state drugs and optic filaments [15], boundary sub caste gets, shock surge conformation and mass transport [12]. Nonlinear kind of equation wich represents the Burgers equation, the warmth conduction equation, the nonlinear Schrödinger equation, the Navier-Stokes equation [5,14]. us, Burgers equation belongs to the class of Navier-Stokes equation and its fundamental PDEs, which find in different model, like fluid mechanics analyses [17]. is Navier-stoke equation is solved by using iteration free backward semi-Lagrangian scheme (see [36]). In this study, the spacial structure for system of equation associated with Ronaldo number is produced. In general, the simplification of this Navier-Stoke equation produces different types of nonlinear equation. Hence, one-dimensional Burger equations was introduced by the scientist Burgers in 1948, and he has attained this equation by simplify the Navier-Stokes equation with dropping its pressure term. So that, this equation, contains the simplest forms of nonlinear advection term and dispersions term [19]. e reason of these studies is, to capture some features of turbulence fluid in a channel caused by the commerce of the contrary goods of convection and diffusion [5,9,17]. In his study, he also illustrates the proposition of turbulences describe by using this two commerce. is condition arises in the proposition of shock swells, in turbulence problems, and non-stop stochastic processes [20], also see [21]. In addition, one-dimensional Burgers equation arises in oil force simulation, modeling of gas dynamics and business flux, turbulence [2,3,20,22], fluid mechanics, tube drugs, biology [15]. Most explicit type methods are best for overcoming nonlinearity but poor for stability [38]. It means that all implicit type methods require lots of computational costs due to nonlinearity but are good in the sense of stability and accuracy [38]. So that, studying result of Burger equation has been carried out for the last half-century and still, it is an active area of disquisition.
Accordingly, more effective and simpler numerical ways are demanded to break Burgers equation. In the literature review, several numerical ways have been applied to probe the numerical results of 1D Burgers equation like [8,11,24,25]. In addition, multitudinous researchers with different numerical approaches answer the coupled and 2D Burgers equations. Alam Khan et al. [12] Presented numerical results of time-fractional Burger equations. In this study, they present comparison between generalized conversion ways/techniques with homogony perturbation manner and take the results of model problem from coincident series terms by using GDTM. In GDTM, the equation and related original conditions are converted into a rush relation that eventually leads to the result of a system of algebraic equations as portions of a power series result. HPM requires a homotopy with an embedding parameter, which is accounted as a frail parameter.
Correspondingly, Yagmurluel.at. [15] presented numerical result of nonlinear coupled system thick burger's equation by using advanced order splitting system/techniques. Since they use this system/techniques within expression finite element scheme by unyoking the system of coupled-Burger equation into two subequations. Hence, the confluence of this system is based on grid point. Again, Khan et al. [13] presented numerical result of one-way surge propagations in Kerr-media.
is study is also use to probe, an operation of Burger equations for surge propagations in rail-life marvels. Azaria and Borhanifar [19] also presented the result of coupled and 2D Burgers equations by using the DTM. DTM is a semi-discretized numerical tactics of logical fashion that formalizes the Taylor series differently. Taylor series system is computationally time-consuming if we use it with large orders. In addition, it is high weakened round out and truncation crimes when its order increases. Ahmad et al. [26] presented numerical results of coupled Burgers' equations. Khater et al. [27] also presents the Chebyshev spectral ways for working the coupled Burgers equations.
is pseudo-spectral techniques watch must be taken with the round-off error issue when advanced derivatives or a large several points Mx and Nx is involved, see [5] and the edge of Chebyshev and wavelate wavelet fashion is worked [28]. In all Chebyshev family techniques, error are generated with the number of grid point are increases (see [35,37]). Mohamed [1] presented numerical result of one and two-dimensional unsteady Burgers' equation by using completely implicit finite difference schemes. Gowrisankar and Natesan [17] are also presented numerical result of singularly perturbed original-boundary Burgers equations by using an effective robust numerical approach. is fashion (techniques) also give A-uniformly convergent numerical results for singularly perturbed Burger's equation, due to actuality of small parameter [5,17].Özis et al. [21] presented Numerical result of Burgers' equation by quadratic B-spline finite elements. However, this system/technique has produced high truncation error, if the number of grid point is veritably large. Indeed this, the delicacy of those numerical ways needs to prove; because the treatments of the ways used to break the Burger equation aren't trivial distributions. Indeed though the delicacy of the forenamed numerical techniques is promising, they bear large memory and long computational time.
erefore, the treatments of these approaches' present severe difficulties that have to be addressed to ensure the accurateness of the numerical result.
To this end, the aims of this paper is, to apply three numerical techniques that are able of working solution of inviscid Burgers equation in the application of physical phenomena. e convergences of present numerical results have been shown (measured) in sense of average absolute error, pointwise absolute error and root mean square error norms, and so that original gets of the numerical result is captured exactly. e stability and confluence (convergency) of the present techniques are also delved by using Von-Neumann stability analysis ways.

Statement of the Problem.
Consider that the following inviscid Burger equation given by: Subjected to initial and boundary conditions respectively where f(x), g(t)) and q(t) are assumed to be sufficiently smooth functions for actuality and existence of result [30,31]. Result U may represent a temperature for heat transfer or a species attention for mass transfer at position x and time t with the advection haste U [14]. Now we define constant of grid points by drawing horizontal and vertical lines of distance 'h' and 'k' singly in 'x' and 't' direction respectively. ese lines are called refinement of mesh line and the point at which they interacting is known as mesh point. Mesh point that lies at end of sphere is called boundary point. e points that lie inside region are called interiors points. erefore, the aim is to compute numerical results of U at each these interior mesh points by applying present three methods and comparing those numerical results obtained at corresponding interior points. Hence to do this, first we discretized the result sphere (interval) as where x j+1 � x j + h and t n+1 � t n + k, j � 0(1)Mx, n � 0(1)Nx. Both Mx and Nt are maximum numbers of grid points in x and t direction, respectively. In addition, the present paper is organized as follows. Section 2 is a 2 International Journal of Mathematics and Mathematical Sciences description of numerical styles/techniques. Section 3 is Stability and confluence (convergence) analysis. Section 4 is Numerical results. Section 5 is Discussions of numerical experiments and Section 6 is the conclusion.

Description of Numerical Methods
In this paper, upwind, Lax-Friedrichs, and Lax-Wendroff styles are used to break (solve) nonlinear in-viscid parabolic partial differential equation given in (1) and comparing the result with each other. ese nonlinear parabolic partial differential equations have the first-order outgrowth in time and first-order outgrowth in the spatial variable are called in thick (in-viscous) Burger equation. Now assuming that U(x, t) has non-stop advanced order partial outgrowth on the region (0, 1)x (0, T. For the sake of simplicity, consider that U(x j , t n ) � U jn , z p U/zx p � z p x U j n and z p U/zt p � z p t U j n for p ≥ 1 is p th order partial derivatives. By using Taylor series expansion, we have

Upwind Methods.
Considering that, nonlinear Burger equation in (1) and using the Taylor expansion in (4) we obtain finite difference discretized scheme of Burger equation in (1) given as: where α � k/h and τ 1 � 1/2(hz 2 x + kz 2 t )U j n is their principal local truncation error respectively in spatial and time discretization of given partial differential equation. Now we refer to (5) as the upwind nonconservative scheme. Although this method is consistent and is adequate for smooth solutions see [16]. But it will not converge in general to a discontinuous weak solution as the grid is refined. Hence, to more understand how this scheme is consistent and convergent simultaneously, you can refer [16]. Now, a time in this work, let us consider that We obtain the general upwind nonconservative scheme of the form: where ψ some function of p + q + 1 argument is called numerical flux function. Methods that conform to this scheme are called conservative methods. If we use the standard finite difference discretizations, we obtain the conservation method called upwind conservative, namely, For p � 0, q � 1, (8) becomes (7). Hence, for ψ(U) � U 2 /2, the upwind conservative scheme for the Burger equation becomes us, general local truncation error for upwind scheme is And it is first-order accuracy which means O(h + k).

Lax-Friedrichs Scheme.
Considering that, nonlinear Burger equation in (1) and using the Taylor series expansion in (4) with (6), we obtain new finite difference discretized scheme of Burger equation given as where β � k/2h, τ 2 � k/2z 2 t U j n and τ 3 � h 2 /6z 3 x U j n are local truncation error term in temporally and spatial direction respectively. In (11), replacing U j,n into the average value of U j− 1,n and U j+1,n and using the definition of ψ(U) given in (6), the Lax-Friedrichs scheme is us the general local truncation error for upwind scheme is And, it is first-order accuracy in time variable and second-order accuracy in the spatial variable which means O(h 2 + k).

Lax-Wendroff Scheme.
Lax-Wendroff's method is an explicit method that uses for the improvement of accuracy in time, see [16]. is method is same times called an explicit of time integration techniques. Now considering in-viscous Burger equation in (1) and using the Taylor series expansion in (4) with refinement of weak solution term (a flux function) given in (6), we obtain new finite difference discretized scheme of Burger equation given as International Journal of Mathematics and Mathematical Sciences where c � k 2 /2h 2 and H j±1/2 is the Jacobean matrix Since for given in-viscous Hence, the required Lax-Wendroff's for given equation becomes where is local truncation error of the proposed scheme. Now, truncating this local truncation error Lax-Wendroff's is as follows: erefore, this method is the third-order accuracy in the time variable and the second-order accuracy in the spatial variable, and this order is given by O(h 2 + k 3 ).

Stability Analysis and Convergence of the Presented Method
Consistency is one of the most important criteria to measure the quantity of numerical computational methods. If the truncation error of the discretized equation is closed to zero, when time and special step sizes are very small, this shows that the discretization equation is consistent. In the numerical discretization scheme, if truncation error contains the expression Δt p /Δx q (i.e., in our case k p /h q , k � Δt and h � Δx where p, q are constantly greater than zero), consistency is satisfied when Δt p is a higher infinitesimal of Δx q . is kind of discrete scheme is the Du-Fort-Frankel scheme [32] and is regarded to be conditional consistency. When truncation error contains Δx q /Δt p , the Lax-Friedrich scheme [33] is conditional consistency under the circumstance that Δx q is the higher-order infinitesimal of Δt p . Hull and White in [34] observe that explicit schemes do not require upper and lower boundary conditions to achieve convergence. is property makes explicit schemes desirable from a computational point of view, see [22]. Again the consistency of Lax-Wondroff and Lax-Friedrich scheme is briefly worked in [16].
Definition: we say the discretization scheme in (9), (12), and (17) are consistent if the limit of truncation error in (10), (11), and (16) is going to zero for which the grid size h, k ⟶ 0. Hence, the considerations of the scheme for the given problem that worked in this paper with their local truncation error are consistence. It means obliviously we can see that for all proposed schemes lim|τ j,n | ⟶ 0 as both mesh-size h, k ⟶ 0.

Stability of Scheme.
e stability of the proposed numerical method is investigated by using Von-Neumann stability analysis. To avoid nonlinear in the proposed scheme, we assumed that V j n � U j,n U j,n . en, without losing generality, we obtain the linear system of a difference equation. Now, we can know to inquire about the eigenvalues of the Mx by Mx linear system of a difference equation. To obtain this eigenvalue, as the working of 4 International Journal of Mathematics and Mathematical Sciences sterility analysis in references [35][36][37][38], we first assume that a trial solution is as follows: Further, to investigate the stability of the proposed method by Von Neumann techniques, we assume that the trial function ϕ(x) defined by the following: where λ a ∈ C, i � �� � − 1 √ , K a � aπ and a � 1(1)Mx , and K is a Fourier number or amplification factor. Now, substituting (18) and (19) into in, and (9) and (12), (17), we obtain the criteria of stability, respectively, as follows, for V j n � U j,n U j,n : λ n+1 a e ijhK a � λ n a e ijhK a − α 2 λ n a e ijhK a − λ n a e ihK a (j− 1) Hence, from (20), we obtain: |λ a | ≪ 1. is shows that the proposed up-wind scheme in (9) is stable with first-order accuracy. Again from (12), we obtain the following: λ n+1 a e ijhK a � 1 2 λ n a e ihK a (j− 1) + λ n a e ihK a (j+1) − β 2 λ n a e ihK a (j+1) − λ n a e ihK a (j− 1) From (21), we have − 1 ≤ cos(hK a ), sin(hK a ) ≤ 1. Hence, |λ a | � |cos(hK a ) − βisin(hK a )| ≤ |cos(hK a )| + β|sin(hK a )| < 1. erefore, the correspondence Lax-Friedrichs scheme is stable with first-order accuracy in the time variable and second-order accuracy in the spatial variable. At last From (17), we obtain the following: λ n+1 a e ijhK a � λ n a e ihK a (j+1) − β 2 λ n a e ihK a (j+1) − λ n a e ihK a (j− 1) λ n a e ijhK a + λ n a e ihK a (j+1) λ n a e ijhK a (j+1) − λ n a e ijhK a λ n a e ijhK a + λ n a e ihK a (j− 1) λ n a e ijhK a − λ n a e ijhK a (j− 1) Since for any value of amplification factor K a , − 1 ≤ cos(hK a ), sin(hK a ) ≤ 1. Hence from (22) us, the correspondence lax-Wendorff scheme is stable with third-order accuracy in the time variable and secondorder accuracy in the spatial variable. Generally, in all cases, the proposed schemes are convergent. erefore, in this paper, the average absolute error (AAE), root mean square error (L 2 ), and maximum absolute error (L ∞ ) norms are used to measure the accuracy of the proposed method, and they are calculated as follows: where U 2j 2n and U j n are numerical solutions of given inviscose Burger equation obtained at the grid point (x j , t n ) where U 2j 2n is a numerical solution obtained at the selected point by doubling the maximum number of grid-point in solution interval.

Results of Numerical Experiments
To test the validity of proposed methods, we have considered one model problem subjected to three different initial conditions associated with different desperation parameters given as follows. (1), which is subjected to discontinuous (shock) initial condition of the following form:

Example 1. Consider the Riemann problem or classical inviscosity Burger equation given in
Example 2. Consider the Riemann problem or classical inviscosity Burger equation given in (1), which is subjected to the piecewise continuous initial condition of the following form: Example 3. Consider the Riemann problem or classical inviscosity Burger equation given in (1), which is subjected to the Gaussian initial condition of the following form:

Discussion
In this paper, the upwind method, Lax-Friedrichs scheme, and Lax-Wendroff scheme are applied to solving onedimensional first-order in-voices Burger equation in the application of physical phenomena. To demonstrate competence of these three methods, one model example subjected to three different initial conditions is solved by taking different values for step size h, and time step k. e computation of numerical results obtained by present methods has been presented in terms of average absolute error, root mean square error, and pointwise maximum absolute error that calculated the result obtained from between two successive iterated numerical solutions, and the results are summarized in tables and graphs. In the present numerical computation, the result presented in Table 1 shows that average absolute error (AAE), roots mean square error (L 2 ), and point wise maximum absolute error norm (L ∞ ) decrease as the number of mesh points Mx and Nt increases in both directions. However, the accuracy of the Lax-Wendroff method rapidly increases and its superior to the accuracy of both the up-wind method and the Lax-Friedrichs scheme. us, an accurate solution is obtained after the breaking solution interval (either time or spatial direction). It means that a solution after the breaking time is to allow refinement of discontinuities of solution U. Such discontinuity is called a shock. For a better understanding about discontinuities (shock) of solution for nonlinear partial differential equation, you can see reference in [16]. Figure 1 shows the physical background of discontinuity solution within the traditional form of a conservation law (the rate of change of the total amount of U in any sections of 0 ≤ x ≤ 1 must be balanced by the net inflow across 0 to 1). To obtain the satisfaction of this condition, we need some shock conditions relating to the jumps of U across the discontinuity, see [16]. Again the result presented in Table 2 also shows that all average absolute error (AAE), roots mean square error norm (L 2 ), and maximum absolute error norm (L ∞ ) decrease rapidly as step length and time step decreases (one of them is constant wise verse) when the model problem is solved subjected to the piecewise continuous initial condition. In addition, in this case, the accuracy of the present method is rapidly increased. However, when we compare the method, the accuracy of Lax-Wendroff is superior to the accuracy of both up-wind method and Lax-Friedrichs scheme. Further, as shown in Figure 2, the proposed method approximates the solution very well for different values of step length h and time step. e result presented in Table 3 also shows that all average absolute error (AAE), roots mean square error norm (L 2 ), and maximum absolute error norm (L ∞ ) decrease rapidly as step length and time step decrease (one of them is constant wise verse) when the model problem is solved with being subjected to Gaussian initial condition. In addition, in this case the accuracy of present method is rapidly increased. 6 International Journal of Mathematics and Mathematical Sciences But when we compare the accuracy of present three methods, the accuracy of all three schemes are uniformly increasing as step length and time step decreases (one of them is constant wise verse), and no superior between the accuracy of this three methods. Hence, also, in this case, Figure 3 shows that the proposed method approximates the      solution very well for Mx � 100 & the time step is Δt � 0.001 . Generally, the solution satisfies the entropy condition and specifying the behavior solution obtained by present methods and the effects of mesh sizes in the solution domain are acceptable. Hence, in all cases the entropy of the conservation laws of system is nondecreasing. us, to further accomplish the accuracy of all three present methods, both the theoretical and numerical error bounds have been established. en, the established error bound for all three method is acceptable. Because the results in the Tables and graphs are further confirmed, computational rate of convergence and this established error bound for all three cases of the applied initial condition. erefore, applications of inviscous-Burger equations are further simulated by these three methods.

Conclusion
Approaches, the upwind, Lax-Friedrichs and Lax-Wendroff schemes are applied to solve in-viscous-Barger equation numerically presented in this study. Generally, the comparison of the accuracy of the numerical solution obtained by these three presents' methods with each other reveals that the Lax-Wendroff method is more convenient, reliable, and effective than both the upwind method, Lax-Friedrichs scheme. As it can be seen from comparisons of the numerical result obtained by all those three numerical techniques in terms of error norm listed in tables, the accuracy of the numerical result obtained by the Lax-Wendroff scheme is faster increase either both Mx and Nt increase or one of them is fixed (either Mx or Nt is fixed) and it is superior to the accuracy of both upwind method and Lax-Friedrichs scheme. Because either both Mx and Nt increase or one of them is fixed (either Mx or Nt is fixed), the error norms of numerical result obtained by Lax-Wendroff scheme rapidly decrease. is shows that the accuracy of the numerical result increases and the accuracy of this numerical result is improved as we compare the accuracy of the numerical result obtained by both upwind method and the Lax-Friedrichs scheme. In summary, the upwind, Lax-Friedrichs, and Lax-Wendroff schemes are capable of solving in-viscous Barger equations. Based on the findings, Lax-Wendroff scheme gives well approximate solutions for the given model problem and gives better accuracy of the numerical solution with a fixed step size h, and decreasing time step k. Step length k:

Abbreviations
Time step Mx: Maximum number grid point in spatial direction.

Data Availability
No data were used to support the findings of this study.

Conflicts of Interest
e author declares that there are no conflicts of interest.

Authors' Contributions
e author planned to apply the present three schemes, initiated developing the Research paper, suggested the experiment, conducted the experiments, and analyzed the empirical results by using MATLAB. e author also reviewed the results and approved the final version of the manuscript.