Numerical Method for Non-Linear Conservation Laws: Inviscid Burgers Equation ()
1. Introduction
The Nonlinear conservation law is used to describe many physical phenomena mathematically, and numerical methods are required to solve these physical problems [1] [2] [3] [4]. The governing one-dimensional conservation law takes the form [5],
(1)
where f(u) is the flux function. The quantity u is the density of some fluid and
treated as mass of the fluid which is conserved in the conservation laws. In the paper, we will discuss integral form of conservation law, classical solution and weak solution for understanding the theoretical knowledge of the problem. To understand the shock wave and rarefaction wave, The Riemann problem also has been discussed. Finally, finite volume method is considered to approximate the numerical solution of the problem.
In finite difference method, the derivatives are approximated by finite differences and therefore, the differential equation does not hold near the discontinuities in the solution. The differential form of the conservation law may not be applied when the shock is generated. After shock formulation, shock actually travels with shock speed [6]. In this paper, finite volume method is considered, which is based on the integral form instead of differential equation. In this method, the domain is divided into grid cells and approximates the total integral of density, u over each grid cell average of density, u, which is the integral divided by the volume of the cell. On the other hand, in finite difference method, the pointwise approximation is needed at the grid point. The cell average is changed in each time step by the flux through the edges of the grid cell. Therefore, it is important to choose a good numerical flux function that approximates the correct fluxes reasonably well [5]. Therefore finite volume method is effective to describe the discontinuous solution. The main objective of this paper is to approximate the numerical solution of the non-linear conservation law by using the inviscid Burgers equation. In this paper, Godunov method is considered to approximate the numerical solution of the problem with continuous and discontinuous initial data. The Godunov method is a conservative numerical method, which is used for the solution of hyperbolic conservation laws. This method can be characterized by the solution (exact or approximate) of a Riemann problem within computational cells in order to obtain numerical fluxes. This method is first-order accurate in both space and time.
Integral Form of Conservation Laws
The integral form of the conservation law is obtained by integrating (1) within the interval
,
Suppose is a conserved quantity in
and defined as,
then,
Again, integrate the above equation over time
,
The above equation is the integral form of the conservation law which directly interpreted that the difference in the total amount of the state variable u in
between in a specific time is equal to the difference in the total flux through the boundary integrated over that time duration [6].
2. Classical Solution
To understand the classical solution, initially we will consider simple linear conservation law, known as linear advection equation [6] where
,
The solution of advection equation can be easily derived from the concept of method of characteristic. Here, the characteristic curve
, describe the straight line with slope a in x = x(t) curve and the solution is
Thus, linear conservation law transport the initial data in the x-direction and also it is exactly as smooth as the initial data.
Now we will consider a simple non-linear conservation laws, in particular, Burgers equation. The study of Burgers equation which is in the class of nonlinear conservation law has been continued to draw the attention of researchers as it represents various physical phenomenon such as theory of shock wave, rarefaction wave, dynamics of fluid or gas etc. [7] [8] [9] [10]. The classical Burgers equation is described by,
In this research, Burgers equation with zero viscosity term is considered which is known as inviscid Burgers equation. The inviscid Burgers equation is a basic case study as it has also the properties of nonlinear conservation law [11]. Therefore inviscid Burgers equation takes the form,
where flus function is,
. The characteristic curve of the equation is
The characteristics are straight lines with slope u, and the solution is
Thus, the solution is constant along the characteristic curves. In this case, discontinuity may form smooth initial conditions in finite time. Therefore, the study of solution of Burgers equation allows interesting phenomena such as shock waves and rarefaction waves.
3. Weak Solution
In classical sense, if u be a piecewise constant solution of the conservation law (1) and also discontinuous along the characteristic curve with any shock speed (s) then u is a weak solution of (1). The shock speed can be derived from the Rankine-Hugoniot jump condition [5] [12]. To get the shock speed, we need to integrate (1) in rectangle
,
Suppose shock splits the rectangle into two parts and
and
are two states immediately to the left and right side of the shock, then
If
, for s < 0 and then dividing by ∆t and taking
gives,
Which is known as the Rankine-Hugoniot Condition and s is the shock speed.
The shock speed for Inviscid Burgers equation is
, as
.
Therefore it is easy to calculate the shock speed. In general, weak solutions occur whenever there is no smooth or classical solution and these solutions may not be differential or even continuous. Therefore, in this case, Equation (1) is not valid and the integral form of conservation law does hold. In mathematically, the function
is a weak solution of conservative law (1) with initial condition if it satisfies the following, for all function of
[13],
However, the weak solution concept cannot guarantee the uniqueness of the solution, because of that addition condition need to be imposed, called entropy condition.
Entropy Condition
The entropy condition described by the following [12],
If a discontinuity propagate with the characteristic speed s, given by the Rankine-Hugoniot condition, satisfies the entropy condition if,
In particular, for the inviscid Burgers equation, if u satisfy the entropy condition then,
where
.
4. Riemann Problem
The Riemann problem is a particular initial value problem which contains a conservation law together with piecewise constant data having a single discontinuity. We consider the Riemann problems for the inviscid Burgers equation subject to special initial conditions, described by
The solution of the Riemann problem depends on
and
. If
then the characteristic curve is,
if
and
if
Then the solution becomes,
Figure 1 shows the characteristic curve for
. Therefore, for
, left wave travel faster than the right wave and as a result, a shock wave is generate. In this case, the solution satisfies the entropy condition also.
On the other hand, if
, the problem has two solutions. First one is,
But this solution does not follow the entropy condition, which is described above [12].
Another solution is weak solution. To get the weak solution suppose
, then the problem becomes,
Figure 1. Characteristic curves for shock wave
.
For
, that is,
On the other hand,
Therefore, weak solution becomes,
Thus, if
, the characteristic waves are spreading and produced rarefaction wave, head of the wave is faster than tail. Figure 2 shows the rarefaction wave.
The study of Riemann problem is very useful for understanding shock wave and rarefaction wave which appear as characteristics in the solution. In numerical methods, Riemann problem is also the basic tool for developing finite volume methods for the solution of conservation law equations due to the discreteness of the grid.
The Riemann problem has a similarity solution, a function of x/t alone; self-similar at different times. The solution
is constant along any ray x/t = constant through the origin.
5. Numerical Method
In this paper, finite volume method is considered for the numerical solution of the one dimensional conservation law. Initially, we have discussed the discretization the computational domain in both space and time of the problem (1).
Figure 2. Characteristic curves for rarefaction wave
.
For simplicity, we have considered the uniform discretization for the special and time domain. Suppose the special domain is
and then, the discretized point is
for
where
. Thus, the i-th grid cell or control volume is
.
The discretization in time with step size ∆t is
.
Now, suppose the value
is the cell average over the i-th interval at time tn. Then,
Similarly,
is the average value at the interface
between tn and tn+1,
Now, integrating the conservation law (1) over
, we get
Therefore,
Dividing both side by ∆x, then
This equation represents the conservation law directly, that is, the rate of change of the cell average is computed by the flux differences through the boundary of the cell. Therefore, this equation is the exact form of conservation law.
In finite volume method, different way of approximation of flux function will give the different methods. Next, the Godunov method is presented by approximating flux function.
5.1. Godunov Method
For Godunov method, we need to approximate the numerical flux
at each interface
. As the cell average
is constant in each cell
at each time level, it may define at each cell interface
a Riemann problem,
According to this method, initially, need to reconstruct a piecewise polynomial function
, from the cell average
. In the simplest case, this is a piecewise constant function that takes
in the i-th grid cell, i.e.
.
Next, after time ∆t, obtain
, with initial data. In this way, need to average this function over every grid cell to calculate new cell averages
. The full process is then repeated in the next time step.
Thus,
is constant over the time interval
. As the Riemann problem has similarity solution. So, Riemann problem centered at
has a similarity solution along
= constant. The value along
gives the value
. Denote this value by
. Therefore, the numerical flux
become
(say)
Therefore, Godunov method is obtained by solving the Riemann problem between states
and
in order to determine the flux
.
Thus, the Godunov flux for conservation law is [5],
Consequently, there are five possible form of Riemann solution and it contains a single shock or rarefaction wave. The solution
is either
if the shock or rarefaction wave moving entirely to the left and on the other hand,
is
if the shock or rarefaction wave moving entirely to the right. Moreover, if the rarefaction wave spreads partly to the left and also partly to the right, at this point, the solution
is neither
nor
and the solution will be
, where
is the unique value for which
where
.
Where s is the shock speed and
.
The formula can also be written as [5] [14],
Additionally, it guarantees that the solution satisfies the entropy condition. However, the necessary condition for consistency of the numerical flux function is that the wave speed should be bounded [14], i.e.
, where L is any constant.
Moreover, the CFL condition for the given problem is
.
5.2. Examples by Godunov Method
Here, the initial value problem for inviscid Burgers equation using Godunov method will be presented by some examples. Both discontinuous and continuous initial value problems are considered.
5.2.1. Discontinuous Initial Value Problems
Initially considering inviscid Burgers equation with the following discontinuous initial value problem,
(2)
The initial condition and the numerical solution using Godunov method is presented in Figure 3 and Figure 4. The numerical solution approximated very well. In this problem, computational time is 0.040 section and also
and
. Figure 4 also shows the wave at different time, where n is represent the number of iteration. The numerical solution does not show any instability. Here it is clear that shock wave propagates with the initial condition.
Again, considering the following discontinuous initial value problem,
(3)
Figure 4. Solution of Burgers equation shows shock wave.
The initial condition and the numerical solution shows the rarefaction wave using Godunov method is presented in Figure 5 and Figure 6, where,
and
. Figure 6 also shows the wave at different time, where n is represent the number of iteration. The computational time is 0.40 second.
5.2.2. Continuous Initial Value Problems
Next, considering inviscid Burgers equation with following continuous initial condition is as smooth as possible. Figures 7-10 show the initial conditions and the numerical solutions by Godunov method of (4) and (5) respectively, where, x – axis is only differ from other.
(4)
(5)
In Figure 8,
and
. In Figure 10,
and
. The computational time is 0.5 second. The speed of propagation of inviscid Burgers equation is u itself. In this case, the initial value moves with variable speed in x. As the sine wave is positive at
, and negative at
, the solution is therefore moves first forward and the backward direction, and develop a discontinuity at
.
Figure 6. Solution of Burgers equation shows Rarefactions wave.
Figure 8. Solution of (4) shows shock wave.
In Godunov method, the solution computed on a very fine mesh. Figure 8 shows initially sine wave compress and develop a shock, and then it expands. In Figure 10, the solution is called N-wave also. Therefore, Godunov method shows a good approximation of the solution.
6. Conclusion
In this research, inviscid Burgers equation is considered as an example of nonlinear conservation law. Here we presented a theoretical feature of nonlinear conservation law by using the common model, inviscid Burgers equation. Some important features such as integral form of conservation law, classical solution, and weak solution with entropy condition also have been presented. The Riemann problem is discussed to understand the shock wave and rarefaction wave. Finite volume method is considered as numerical method to approximate the solution, as it is effective to describe the discontinuous solution. Taking into account a special initial value problem such as continuous and discontinuous initial value problems are solved by Godunov method. In conclusion, Godunov method shows a good approximation of the problem.