Scalar Conservation Laws

We present a theoretical aspect of conservation laws by using simplest scalar models with essential properties. We start by rewriting the general scalar conservation law as a quasilinear partial differential equation and solve it by method of characteristics. Here we come across with the notion of strong and weak solutions depending on the initial value of the problem. Taking into account a special initial data for the left and right side of a discontinuity point, we get the related Riemann problem. An illustration of this problem is provided by some examples. In the remaining part of the chapter, we extend this analysis to the gas dynamics given in the Euler system of equations in one dimension. The transformations of this system into the Lagrangian coordinates follow by applying a suitable change of coordinates which is one of the main issues of this section. We next introduce a first-order Godunov finite volume scheme for scalar conservation laws which leads us to write Godunov schemes in both Eulerian and Lagrangian coordinates in one dimension where, in particular, the Lagrangian scheme is reformulated as a finite volume method. Finally, we end up the chapter by providing a comparison of Eulerian and Lagrangian approaches.


Introduction
We present a general form of scalar conservation laws with further properties including some basic models and provide examples of computational methods for them.The equations described by in one dimension are known as scalar conservation laws where u ¼ u t; x ð Þ is the conserved quantity and f ¼ f u ð Þ is the associated flux function depending on t and x.Whenever an initial condition u 0; x ð Þ ¼ u 0 x ð Þ is attached to Eq. ( 1), the problem is called the Cauchy problem the solution of which is a content of this chapter.The outlook of chapter is as follows.We introduce basic concepts and provide particular examples of scalar conservation laws in the first part.The equation of gas dynamics in Eulerian coordinates in one dimension is the main issue of the second part.After providing further instruction for these equations, we provide a transformation of the Eulerian equations in the Lagrangian coordinates.In the final part, we give as an example of computational methods for conservation laws, the Godunov schemes for the Eulerian, and the Lagrangian coordinates, respectively.

Conservation laws: integral form and differential form
We start by investigating the relation of the equations in gas dynamics with conservation laws.We take into account the equation of conservation of mass in one dimension.The density and the velocity are assumed to be constant in the tube where x is the distance and ρ t; x ð Þis the density at the time t and at the point x.
The last equation is called integral form of conservation law.Integrating this expression in time from t 1 to t 2 , we get Using the fundamental theorem of calculus after reduction of Eq. ( 3), it follows that As a result, we get Here the end points of the integrations are arbitrary; that is, for any x 1 ; x 2 ½ and t 1 ; t 2 ½ , the integrant must be zero.It follows that the conservation of mass yields which is said to be the differential form of the conservation law.

A first-order quasilinear partial differential equations
A general solution to a quasilinear partial differential equation of the form where a, b, c are non-zero and smooth on a given domain D ∈ R 3 follows by the characteristic method where the characteristic curves are defined by By applying a parametrization of c, the relation ( 8) is transformed to a system of ordinary differential equation (ODE): In addition to these equations, if an initial condition u 0 ¼ u x 0 ð Þ is also given, then by the existence theorem of ODE, there is a unique characteristic curve passing from each point t 0 ; x 0 ; u 0 ð Þleading to an integral surface which is the solution to Eq. (7).
Observe that the scalar conservation law (1) is a particular example of Eq. ( 7) if we assign a t; x , and c t; x; u ð Þ¼0.The conserved quantity can be observed by integrating equation (1) over x 0 ; x 1 ½ .Indeed This means, the quantity u t; x ð Þ is conserved so that it depends on the difference of the flux functions between the points x 0 and x 1 :

Strong (classical) solutions
We consider the initial value problem where the initial data is assumed to be continuously differentiable, that is, Applying the chain rule to the relation (11), it follows that where we define characteristic curves of Eq. ( 12) to be the solution of Then a solution to the system (12) in a domain Ω ∈ R is said to be a strong (or classical) solution if it satisfies Eq. ( 11), and it is continuously differentiable on a domain Ω ∈ R: Let u be a strong solution and the initial data u 0 be differentiable.Observe that ( 12) is equivalent to a quasilinear form: Applying the method of characteristics to Eq. ( 13), the partial differential equation is transformed to a system of ordinary differential equations.We consider the characteristic curve passing through the point 0; x 0 ð Þ: Along this characteristic curve, is satisfied, that is, u is constant.Hence, the characteristic curves are straight lines satisfying Hence we can define smooth solutions by Þis increasing, then the slopes of the characteristics are decreasing.As a result, the characteristics do not intersect, and thus solution can be defined for all t which is greater than zero.On the other hand, if λ u 0 x ð Þ ð Þis decreasing, then the slopes of the characteristics will be increasing which implies that the characteristics intersect at some point.But at the intersection point, solution cannot take both values u 0 x 1 ð Þ and u 0 x 2 ð Þ.Therefore, we cannot define the strong solution for all t .0.

Linear advection equation
The basic example of the scalar conservation law is the linear advection equation.It can be obtained by setting a t; x; u ð Þ¼1, b t; x; u ð Þ¼λ, and c t; x; u ð Þ¼0 in Eq. (7).The flux function takes the form f u ð Þ ¼ λu where λ is a constant.Then the following quasilinear partial differential equation is a linear advection equation.Similar to Eqs. ( 11) and ( 12), an initial value problem for linear advection equation is described by Applying the method of characteristics, it follows that dt 1 ¼ dx λ ¼ du 0 or equivalently where c 1 and c 2 are constant and x À λt ¼ c 2 : As a conclusion, the solution is Here λ is the wave speed, and the characteristic lines x À λt ¼ c 2 are wavefronts which are constants.

Burgers' equation
Burgers' equation is the simplest nonlinear partial differential equation and is the one of the most common models used in the scalar conservation laws and fluid dynamics.The classical Burgers' equation is described by Advanced Computational Fluid Dynamics for Emerging Engineering… where ν∂ xx u is the viscosity term.Equation ( 21) can be considered as a combination of nonlinear wave motion and linear diffusion term so that it is balance between time evolution, nonlinearity, and diffusion.The term u∂ x u is a convection term that may have an effect to wave breaking, and the term ν∂ xx u is a diffusion term that may cause to efface the wave breaking and to flatten discontinuities, and thus we expect to achieve a smooth solution.We try to find a traveling wave solution of Eq. ( 21) of the form where g and λ are to be determined.Applying the chain rule, we get Plugging these terms in Eq. ( 21), we get Taking integration with respect to ξ gives Rewriting Eq. ( 25) by p .Supposing that g 1 , g 2 are real implies g 1 .g 2 .Using separation of variable and then integrating equation (26), we get As a result the explicit form of traveling wave solution of Eq. ( 21) becomes where λ ¼ 2 is the wave speed.We can observe that lim This means the solution g ξ ð Þ decreases monotonically with ξ from the value g 1 to g 2 .At ξ ¼ 0, u ¼ 2 ¼ λ, that is the wave form g ξ ð Þ travels from left to right with speed λ equal to the average value of its asymptotic values.The solution resembles to a shock form as it connects the asymptotic states g 1 and g 2 .Without the viscosity term, the solutions to Burgers equation allow shock forms which finally break.The diffusion term prevents incrementally deformation of the wave and its breaking by withstanding the nonlinearity.As a conclusion, there exists a balance between nonlinear advection term and the linear diffusion term.The wave form is notably affected by the diffusion coefficient ν.If ν is smaller, then the transition layer between two asymptotic values of solution is sharper.In the limit ν !0, the solutions converge to the step shock wave solutions to the inviscid Burgers' equation.
Scalar Conservation Laws DOI: http://dx.doi.org/10.5772/intechopen.83637 Remark.If the initial data is smooth and very small, then the u xx term is negligible compared to other terms before the beginning of wave breaking.As the wave breaking starts, the u xx term raises faster than u x term.After a while, the term u xx becomes comparable to the other terms so that it keeps the solution smooth, giving rise to avoid breakdown solutions.

Inviscid Burgers' equation
Whenever ν ¼ 0, Eq. ( 21) is called the inviscid Burgers' equation.This equation can be obtained by substituting f u ð Þ ¼ u 2 =2 in the scalar conservation law (1), that is Observe that f u ð Þ is a nonlinear function of u; thus, the inviscid Burgers' equation is a nonlinear equation.Equation ( 29) is now equivalent to Eq. ( 17) with λ ¼ u.We know the solution of Eq. ( 17); so, plugging λ ¼ u into the relation (20) implies that the solution of Eq. ( 29) is Recall that the characteristic speed λ is constant for linear advection equation; that is, the characteristic curves become parallel for Eq. ( 17).In contrast, for the inviscid Burgers' equation ( 29), the characteristic speed λ ¼ u depends on u.As a result the characteristic lines are not parallel.If we apply the implicit function theorem to Eq. ( 29), the solution can be written as a function of t and x as u 0 is differentiable.More particularly, differentiating Eq. (30) with respect to t, we get and differentiating equation (30) with respect to x, we get Thus, substituting Eqs. ( 31) and ( 32) in (29), we can recover the inviscid Burgers' equation.Consequently, the relations (31) and (32) imply that the solutions of Eq. ( 1) and particularly of Eq. ( 29) depend on the initial value u 0 .It can be observed that whenever u 0 0 x ð Þ .0, then by Eq. (32), ∂ x u decreases in time because 1 þ u 0 0 t .0 for t .0. In other words, the profile of the wave flattens as time increases.On the other hand, whenever u 0 0 x ð Þ , 0, then ∂ x u increases in time as 1 þ u 0 t , 0: Hence u x in Eq. (32) tends to ∞ as 1 þ u 0 0 t approaches to zero.As a result, wave profile become sharp after some time.For further details on the Burgers' equations, we refer the reader to [12,13,22] and the references therein.

Shock waves
Let the constants u L and u R are given with a linear function, is a simple example of discontinuous solution of the conservation law (11).If u L 6 ¼ u R , the relation ( 33) is called a shock wave connecting u L to u R with shock speed λ.As an example, if we take into account the characteristics of the inviscid Burgers' equations which are of the form dx dt ¼ u t; x ð Þ, it follows that thus, the characteristics are straight lines.Depending on the behavior of these characteristics, we have two cases.If u L .u R , characteristics intersect, the solution will have an infinite slope, and the wave will break; as a result a shock is obtained.This is illustrated in Figure 1.On the other hand, if u R .u L , the characteristics do not intersect, and hence a region without characteristic will appear which is physically unacceptable.This is shown in Figure 2. We get rid of this by introducing the rarefaction waves.

Rarefaction waves
A rarefaction wave is a strong solution which is a union of characteristic lines.A rarefaction fan is a collection of rarefaction waves.These waves are constant on the characteristic line where u L and u R are the values of u at the edge of the rarefaction wave fan.If moreover f 0 is invertible, then the solution If, for instance, f is convex, then the rarefaction waves are increasing.If we consider again the inviscid Burgers' equation with the initial values, then the region without characteristics in Figure 2 will be covered by rarefaction solution which is described by  Scalar Conservation Laws DOI: http://dx.doi.org/10.5772/intechopen.83637 An illustration of rarefaction waves and rarefaction fan in Eq. ( 36) is given in Figure 3.
Remark.Whenever characteristics intersect, we may have multiple valued solution or no solution; but we have no more classical (strong) solution.To get rid of this situation, we introduce a more wide-ranging notion of solution, the weak solution, in the next part.By this arrangement, we may have non-differentiable and even discontinuous solutions.

Weak solution
Weak solutions occur whenever there is no smooth (classical) solution.These solutions may not be differentiable or even not continuous.Considering ϕ : R Â R þ !R as a smooth test function with a compact support and multiplying the scalar conservation law (1) by this test function ϕ, it follows after integration by parts that Putting the initial condition Observe that there are no more derivatives of u and f which may lead less smoothness.In other words, the smoothness requirement is reduced for finding a solution.Thus, the function u t; x ð Þ is said to be the weak solution of the initial value problem (11) if the relation (38) satisfied for all test function ϕ: Here it is significant to note that u needs not be smooth or continuous to satisfy Eq. (38).Consequently, by weak solutions, we extend the solutions so that discontinuous solutions may also be covered.However, in general weak solutions are not unique.We can also notice Advanced Computational Fluid Dynamics for Emerging Engineering… that strong solutions are also weak solutions and a weak solution which is continuous and piecewise differentiable is also strong solution.

Riemann problem
The Riemann problem is a Cauchy problem with a particular initial value which consists a conservation law together with piecewise constant data having a single discontinuity.We consider the Riemann problem for a convex flux described by The solution is a set of shock and rarefaction waves depending on the relation between u L and u R : There are two cases to investigate: Case 1: u L .u R ð ÞA shock is obtained because the left-hand side wave moves faster than the right-hand side one.Thus the solution is a shock wave satisfying the shock speed λ ¼ The solution given in Case 1 is also a solution for this case.In addition, we have rarefaction solutions of the form (36) illustrated by Figure 3.

Rankine-Hugoniot jump condition
A jump discontinuity along the characteristic line is controlled by the Rankine-Hugoniot jump condition.Integrating the scalar conservation law (1) in Suppose that there is a discontinuity at the point Þwhere u and u 0 are continuous on the , respectively.Suppose also that whenever , their limits exist.Next, Eq. ( 41) can be rewritten as By the fundamental theorem of calculus, the relations (41) and (42) yield Taking the limit whenever Scalar Conservation Laws DOI: http://dx.doi.org/10.5772/intechopen.83637 The relation ( 44) is said to be the Rankine-Hugoniot jump condition.Geometrical meaning of the Rankine-Hugoniot jump condition is that the shock speed is the slope of the secant line through the points u Þon the graph of f .

Entropy functions
Entropy and entropy flux are defined for attaining physically meaningful solutions.If u is the smooth solution of the conservation law (1), then the relation is satisfied for continuously differentiable functions G and F where the pair G; F ð Þis called as entropy pair so that G is entropy and F is entropy flux.If in addition u is smooth, then Eq. ( 45) becomes which looks like to the scalar conservation law (1).Indeed, if we multiply Eq. ( 1) by G 0 u ð Þ, it follows that It follows that Eqs. ( 46) and ( 47) are equivalent with Here the function u t; x ð Þ is said to be the entropy solution of Eq. ( 1) if

Entropy condition
Weak solutions to conservation laws may contain discontinuities as a result of a discontinuity in the initial data or of characteristics that cross each other or because of the jump conditions which are satisfied across the discontinuities.Although the Rankine-Hugoniot jump condition is satisfied, the uniqueness of the solution may always not be guaranteed.In order to eliminate the nonphysical solutions among the weak solutions, we need an additional condition, so-called entropy condition.It is described by the following: A discontinuity propagating with the characteristic speed λ given by the Rankine-Hugoniot jump condition satisfies the entropy condition if holds.
Example 1.1.The weak solutions to conservation laws need not be unique.If we write the inviscid Burgers' equation in quasilinear form and multiply by 2u, we obtain 2u∂ t u þ 2u 2 ∂ x u ¼ 0. In conservative form it becomes The inviscid Burgers' equation and Eq. ( 49) have exactly the same smooth solutions.But their weak solutions are different.A shock traveling speed for the inviscid Burgers' equation is λ ; however for Eq. ( 49), we have Applying the method of characteristics for t .0, it follows that Next if we integrate Eq. ( 51) with respect to t, we get the characteristic curves where c .0 and b are constants.Due to the discontinuity at the point x ¼ 0, there is no strong (classical) solution.The speed of propagation is λ ¼ u L þu R 2 ¼ 0:5: Moreover, the weak solution for t ≤ λ ¼ 0:5 becomes which satisfies both the jump condition and the entropy condition as u L ¼ 1 .u R ¼ 0. The characteristic curves can be observed in Figure 4.
Example 1.3.We now interchange the roles of u L and u R of the Example 1.2 so that u L , u R to get an initial value problem: By the method of characteristics, we obtain a solution Scalar Conservation Laws DOI: http://dx.doi.org/10.5772/intechopen.83637which is a classical (strong) solution on both sides of the characteristic line x t ¼ 1.Since it satisfies the Rankine-Hugoniot jump condition along the discontinuity curve, it is a weak solution.However, the entropy condition is not satisfied.It yields an empty region between the characteristic lines shown in Figure 4.In order to cover this empty state, we consider another solution described by which satisfies both jump and entropy conditions.Here we can observe the rarefaction fan arising on the interval 0 ≤ x t ≤ 1.An illustration of this solution is supplied in Figure 5.

The gas dynamic equations in one dimension
The equation of fluid dynamics can be represented in Eulerian and Lagrangian forms.Eulerian coordinates are related to the coordinates of a fixed observer.On the other hand, Lagrangian coordinates are in usual related to the local flow velocity.That is, due to the velocity taking different values in different parts of the fluid, the change of coordinates is different from one point to another one.

Eulerian coordinates
The equations of gas dynamics in Eulerian coordinates can be written in the following conservative forms: where we ignored the heat conduction.If we denote then Eq. ( 57) can be written by where ρ is density, p is pressure, u is velocity, and e is the specific internal energy.

Hyperbolicity of the Euler system
If we do not neglect the heat conduction, then the U and F terms in Eq. ( 59) become where Þρ , and for perfect gases δ ¼ c p =c v is the ratio of specific heats.Rewriting Eq. ( 59) in quasilinear form, we get where where a is the sound speed given by a ¼ ffiffiffi ffi δp ρ q .
Moreover the corresponding eigenvectors are which are real, and the eigenvectors are linearly independent implying that the Euler equations for perfect gases are hyperbolic.

Rankine-Hugoniot conditions for the Euler system
Using the results in the previous part, the Rankine-Hugoniot jump conditions for the Euler system will be of the form where the indices 1 and 2 refer to the left and right of the shock, respectively, and s denotes the wave speed.

Riemann problem for the Euler system
The Riemann problem for the one-dimensional Euler equation ( 57) is represented by The reader is addressed to the references [18,24] for further details.

Lagrangian coordinates
We aim to transform the equations of gas dynamics (57) given in the Eulerian coordinates into the Lagrangian coordinates for one-dimensional case.We start denoting by u ¼ u t; x ð Þ the velocity field of the fluid flow and consider the differential system We set the following change of coordinates from Euler coordinates to Lagrange coordinates for space and time as t Þare the Lagrangian coordinates associated with the velocity field u.We set which gives It follows by some algebraic manipulations that the gas dynamic equations become In order to derive a more convenient form of the system (69), we derive firstly the equation of conservation of mass: where Assuming that ρ .0, we introduce the specific volume τ ¼ 1=ρ, and by using Eq. ( 68) we get which yields Hence the second and third equations of Eq. ( 69) become Advanced Computational Fluid Dynamics for Emerging Engineering… Moreover, we define a mass variable m by Finally, using Eqs.( 69) and (73), the Euler system (57) can be written in Lagrangian coordinates with the mass variable in the form where . 0, we obtain a scalar conservation law of the form which is strictly hyperbolic.This can be verified by checking the Jacobian matrix of the flux calculated with respect to the variables τ; so that they are all distinct, and thus the system is strictly hyperbolic.In fact there are different versions of the gas dynamics in Lagrangian coordinates.In this part we followed the approaches stated in [9,10,12].For further details we cite these works with references therein.

Rankine-Hugoniot conditions for the Lagrangian system
Similarly as in the Euler system, the Rankine-Hugoniot jump conditions for the Lagrangian system (79) are of the form where σ denotes the speed of propagation of the discontinuity with respect to the mass variable.
Remark.The Eulerian and Lagrangian Rankine-Hugoniot relations are equivalent.Moreover, Eulerian entropy relations are equivalent to all Lagrangian entropy relations (see [9] for further detail).
Example 2.1.For simplicity of notation, we take t; x ð Þ as the Lagrangian coordinates.Then the system of equations is a one-dimensional isentropic gas dynamics in Lagrangian coordinates which is also known as p-system.It is the simplest nontrivial example of a nonlinear system of conservation laws.Here τ is the specific volume, u is the velocity, and the pressure p ¼ p τ ð Þ is given as a function of τ by The system (79) is equivalent to where τ .0 and τ; u ð Þ∈ R 2 : If we assume that p 0 τ ð Þ , 0, it follows that the Jacobian matrix of f . In other words, the system (81) is strictly hyperbolic.On the other hand, for the case p 0 τ ð Þ .0, it becomes elliptic.Moreover, one can verify that the solutions of the p-system (79) and the Euler system (57) are equivalent.

Godunov schemes
The Godunov scheme deals with solving the Riemann problem forward in time for each grid cell and then taking the mean value over these cells.The Riemann problem is solved per mesh point at each time step iteratively.If there are no strong shock discontinuities, this process may cost much and will not be effective.To get rid of such a situation, we establish approximate Riemann solvers that are easier to implement and also low cost to use.Eulerian and Lagrangian Godunov schemes are current Godunov scheme in literature.Both have advantages and disadvantages depending on the structure of the problem.A brief comparison of the method for these two approaches is presented in the last part of the chapter.In this work we will not go further in numerical examples and details of these methods; instead, we aim to present a general form of Godunov schemes for gas dynamics in Eulerian and Lagrangian coordinate.Before introducing these, we present a first-order Godunov scheme for scalar conservation laws.

First-order Godunov scheme
Consider the scalar conservation law (1).Godunov scheme is a numerical scheme which takes advantage of analytical solutions of the Riemann problem for the conservation law (1).The numerical flux functions are evaluated at the spatial steps x jÀ1=2 and x jþ1=2 by handling the solutions of the Riemann problem.On each grid cell where λ is the wave propagation speed.Hence the numerical flux for Godunov's method can be generalized by For numerical illustration of Godunov schemes, we cite the articles [14, 20, 27].

Godunov method in Eulerian coordinates
We consider Eq. ( 59) with (60).The eigenvalues of Then the Riemann problem at the point x iþ1=2 between the states U i and U iþ1 which is solved by the Godunov scheme can be written by

Godunov method in Lagrangian coordinates
Consider the initial condition for a quantity v given by the mean value The eigenvalues satisfy σ 1 , σ 2 ¼ 0 , σ 3 : Setting u iþ1=2 and p iþ1=2 as the values of u and p at the contact discontinuity between V n i and V n iþ1 , it follows that Then Godunov scheme for the Lagrangian coordinates takes the form 18 Advanced Computational Fluid Dynamics for Emerging Engineering… where If we now consider the moving coordinates, Godunov scheme can also be derived equivalently by the following.Setting x iþ1=2 ¼ ξ iþ1=2 with the approximation of u ¼ dx=dt, it follows that the Eulerian coordinate x iþ1=2 of the interface ξ iþ1=2 at t n is upgraded with respect to Next we deduce by a simple induction process.Hence the Lagrangian Godunov schemes become with Notice that the Lagrangian Godunov schemes can be reformulated as a finite volume method.Equation (100) can be written in conservative form: If we integrate these equations on ξ Here we omit the dependency of f , φ and x on t.Moreover, if we suppose that φ is constant in each cell ξ iÀ1=2 ; ξ iÀ1=2 À Á , it follows by an explicit one-step method that is Moreover, if ρ; u; e ð Þare constant in each cell with v ¼ u, we get the Godunov scheme: Scalar Conservation Laws DOI: http://dx.doi.org/10.5772/intechopen.83637 provided u n iþ1=2 ; p n iþ1=2 are determined by the solution of the Riemann problem, which is the desired result.

Comparison of Eulerian and Lagrangian schemes
In the literature there are two types of Godunov schemes: the Eulerian and Lagrangian.To compare one with the other, both have advantages and disadvantages.These are briefly listed in the following:

Eulerian approach
It is more nature; that is the properties of a flow field are described as functions of the coordinates which are in the natural physical space and time.The flow is determined by examining the behavior of the functions.Eulerian coordinates correspond to the coordinates of a fixed observer.This approach is ease of implementation and computation.The computational grids derived from the geometry constraints are generated in advance.The computational cells are fixed in space, and the fluid particles move across the cell interfaces.Since the Eulerian schemes consider the implementation at the nodes of a fixed grid, this may lead to spurious oscillations for the problems like diffusion-dominated transport equations.By adding artificial diffusion, one can get rid of these oscillations; however the nature of the problem may differ from the original one.Besides, refining the grids may also lead to remove numerical oscillations, but this process may augment the computation cost.Besides, while refining the grids, it may cause restriction of the size of time step which is limited by CFL condition.This restriction does not occur in Lagrangian case.

Lagrangian approach
It is based on the notion of mass coordinate denoted by m ξ ð Þ.An important feature of the mass coordinate is that two segments have the same length if the mass contained in these segments is the same.This leads to face with a disadvantage; that is, at each iteration time step, the problem has to be converted from the natural coordinate system to the mass coordinate system.Once the solution at the next step is known, it has to be remapped into the natural coordinate system.As a result, this process raises the cost of the computation.Lagrangian coordinates are associated to the local flow velocity.In other words, as the velocity has different values in different parts of the fluid, then the change of coordinates is different from one point to another one in Lagrangian coordinates.Thus Lagrangian coordinates are equivalent to the Eulerian coordinates at another time.Lagrangian description states the motions and properties of the given fluid particles as they travel to different locations.Hence the computational grid points are precisely fluid particles.Since the particle paths in steady flow coincide with the streamlines, no fluid particles will cross the streamlines.Hence, there is no convective flux across cell boundaries, and the numerical diffusion is minimized.As a result, Godunov method in a Lagrangian grid is easier to handle.Moreover, in the case of higher schemes, the Advanced Computational Fluid Dynamics for Emerging Engineering… subsonic character of the flow makes the transformation much easier than in Eulerian schemes.Lagrangian schemes consider the implementation in a grid that moves with the flow which is an advantage for the problems like the transport equations since the advective and diffusion terms can separately be examined.
Apart from the two main approaches, there is another method which is a combination of both, so-called Eulerian-Lagrangian methods.It combines the advantages and eliminates disadvantages of both approaches to get a more efficient method.For further details we address the reader to the reference in the next part.

Notes
We have tried to present only the theoretical aspects of scalar conservation laws with some basic models and provide some examples of computational methods for the scalar models.There are plenty of contributors to the subject; however, we just cite some important of these and the references therein.Scalar conservation laws are thoroughly studied in particular in [12]; for a more general introduction including systems, see [13,15,18,19,22] and the references therein.There are some important works related to the concept of entropy provided by [7,15,16].A more precise study of the shock and rarefaction waves can be found in [23].A simple analysis for inviscid Burgers' equation is done by [21].The readers who are deeply interested in systems of conservation laws and the Riemann problem should see [8,13,15,22,24].A well-ordered work of the propagation and the interaction of nonlinear waves are provided by [26].We refer the reader to the papers [1,17] for the theory of hyperbolic conservation laws on spacetime geometries and finite volume analysis with different aspects.A widely introductory material for finite difference and finite volume schemes to scalar conservation laws can be found in [18].In this chapter we have studied the one-dimensional gas dynamics on the Eulerian and Lagrangian coordinates.For the detail on the Lagrangian conservation laws, we refer [10]; moreover for both Eulerian and Lagrangian conservation laws, we cite [11].The proof of the equivalency of the Euler and Lagrangian equations for weak solutions is given in [25].There are several numerical works for Lagrangian approach; some of the basic works on Lagrangian schemes are given in [2][3][4][5][6].We refer the reader to the book [7] for a detailed analysis of the mathematical standpoint of compressible flows.Moreover Godunov-type schemes are precisely analyzed in [14,27]; whereas, Lagrangian Godunov schemes can be found in [2,12,20].As a last word, we must cite [9] as a recent and more general book consisting of scalar and system approaches of both Eulerian and Lagrangian conservation laws with theoretical and numerical parts which can be a basic source for the curious readers.

Figure 1 .
Figure 1.For the initial value u L .u R , characteristics, and shock wave.

Figure 2 .
Figure 2.For the initial value u R .u L , characteristics and rarefaction waves.

Figure 4 .
Figure 4.For initial value u L .u R , the characteristic solutions.
Then if we integrate the density on x 1 ; x 2 and thus these two equations have different weak solutions.Example 1.2.We first consider the initial value problem for u L .u R given by we have a piecewise constant function.The Riemann problem for (1) for the left and right sides of I i are described by Scalar Conservation Laws DOI: http://dx.doi.org/10.5772/intechopen.83637for each u n j .Next, if assigning u * as the intermediate value over the grid I i in the Riemann solution, it implies that