Galilean invariance and the conservative difference schemes for scalar laws

Transcription

1 RESEARCH Open Access Galilean invariance and the conservative difference schemes for scalar laws Zheng Ran Correspondence: edu.cn Shanghai Institute of Applied Mathematics and Mechanics, Shanghai University, Shanghai 0007, P. R. China Abstract Galilean invariance for general conservative finite difference schemes is presented in this article. Two theorems have been obtained for first- and second-order conservative schemes, which demonstrate the necessity conditions for Galilean preservation in the general conservative schemes. Some concrete application has also been presented. Keywords: difference scheme, symmetry, shock capturing method 1. Introduction For gas dynamics, the non-invariance relative to Galilean transformation of a difference scheme which approximates the equations results in non-physical fluctuations, that has been marked in the 1960s of the past century [1]. In 1970, Yanenko and Shokin [] developed a method of differential approximations for the study of the group properties of difference schemes for hyperbolic systems of equations. They used the first differential approximation to perform a group analysis. A more recent series of articles was devoted to the Lie point symmetries of differential difference equations on [3]. In a series of more recent articles, the author of this article has used Lie symmetry analysis method to investigate some noteworthy properties of several difference schemes for nonlinear equations in shock capturing [4,5]. It is well known that as for Navier-Stokes equations, the intrinsic symmetries, except for the scaling symmetries, are just macroscopic consequences of the basic symmetries of Newton s equations governing microscopic molecular motion (in classical approximation). Any physical difference scheme should inherit the elementary symmetries (at least for Galilean symmetry) from the Navier-Stokes equations. This means that Galilean invariance has been an important issue in computational fluid dynamics (CFD). Furthermore, we stress that Galilean invariance is a basic requirement that is demanded for any physical difference scheme. The main purpose of this article is to make differential equations discrete while preserving their Galilean symmetries. Two important questions on numerical analysis, especially important for shock capturing methods, are discussed from the point view of group theory below. (1) Galilean preservation in first- second-order conservative schemes; () Galilean symmetry preservation and Harten s entropy enforcement condition [6]. The structure of this article is as follows. First, the general remarks on scalar conservation law and its numerical approximation are very briefly discussed in Section, 011 Ran; licensee Springer. This is an Open Access article distributed under the terms of the Creative Commons Attribution License ( which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

2 Page of 16 while Section 3 is devoted to the theory of symmetries of differential equations. The following sections are devoted to a complete development of Lie symmetry analysis method proposed here and its application to some special cases of interest. The final section contains concluding remarks.. Scalar conservation laws and its numerical approximation In this article, we consider numerical approximations to weak solutions of the initial value problem (IVP) for hyperbolic systems of conservation laws [6,7] u t + f (u) x =0,u(x,0)=u 0 (x), < x < +. (:1) where u(x, t) is a column vector of m unknowns, and f(u), the flux, is a scalar valued function. Equation.1 can be written as u t + a(u)u x =0,a(u) = df du, (:) which asserts that u is constant along the characteristic curves x = x(t), where dx dt = a (u). (:3) The constancy of u along the characteristic combined with (.3) implies that the characteristics are straight lines. Their slope, however, depends upon the solution and therefore they may intersect, and where they do, no continuous solution can exist. To get existence in the large, i.e., for all time, we admit weak solutions which satisfy an integral version of (.1) 0 [ wt u + w x f (u) ] dxdt + w (x,0) u 0 (x) dx =0. (:4) for every smooth test function w(x, t) of compact support. If u is piecewise continuous weak solution, then it follows from (.4) that across the line of discontinuity the Rankine-Hugoniot relation f (u R ) f (u L ) = s (u R u L ). (:5) holds, where s is the speed of propagation of the discontinuity, and u L and u R are the states on the left and on the right of the discontinuity, respectively. The class of all weak solutions is too wide in the sense that there is no uniqueness for the IVP, and an additional principle is needed for determining a physically relevant solution. Usually this principle identifies the physically relevant solution as a limit of solutions with some dissipation, namely u t + f (u) x = ε[β (u) u x ] x. (:6) Oleinik [8] has shown that discontinuities of such admissible solutions can be characterized by the following condition: f (u) f (u L ) u u L s f (u) f (u R) u u R. (:7) for all u between u L and u R ; this is called the entropy condition, or Condition E. Oleinik has shown that weak solutions satisfying Condition E are uniquely determined

3 Page 3 of 16 by their initial data. We shall discuss numerical approximations to weak solutions of (.1) which are obtained by (K+1) -point explicit schemes in conservation form u n+1 j = u n j λ f n 1 f n 1, (:8) j+ j where f n j+ 1 = f ( ) u n j K+1,..., un j+k. (:9) where u n j = u ( j x, n t ), and f is a numerical flux function. We require the numerical flux function to be consistent with the flux f(u) in the following sense: f (u,..., u) = f (u). (:10) We note that f is a continuous function of each of its arguments. Let f r = f r, r = K + 1,..., K (:11) f K =0, (:1) f k+1 =0. (:13) Equation.8 can be written as follows: u n+1 j = u n j λ f n 1 f ( ) n 1 G u n j K,..., un j+k. (:14) j+ j It follows from (.14) that ( ) ( ( ) ( )) G u n j,..., un j = u n j λ f u n j f u n j = u n j. (:15) Suppose that G is a smooth function of its all arguments, then G r = G r, (:16) G rs = G r s. (:17) At last, one can derive the conservation form scheme approximation solutions of the viscous modified equation [9,10]

4 Page 4 of 16 where u t + f (u) x = 1 t x [β (u, λ) u x]. (:18) β (u, λ) = 1 λ K r= K ( ) f r G r. (:19) We claim that, except in a trivial case, b(u, l) 0andb(u, l) 0; this shows that the scheme in conservative form is of first-order accuracy [9-11]. 3. Mathematical preliminaries on Lie group analysis All the problems to be addressed here can be described by a general system of nonlinear differential equations of the nth order ν ( x, u (n) ) =0, (3:1) where v = 1,...,l and x =(x 1,...,x p ) Î X are independent variables, u =(u 1,...,u q ) Î U are dependent variables, and Δ v (x, u (n) )=(Δ 1 (x, u (n) ),..., Δ l (x, u (n) )) is a smoothing function that depends on x, u and derivatives of u up to order n with respect to x 1,...,x p.if we define a jet space X U (n) as a space whose coordinates are independent variables, dependent variables and derivatives of dependent variables up to order n then Δ is a smoothing mapping : X U (n) R l. (3:) Before studying the symmetries of difference schemes, let us briefly review the theory of symmetries for differential equations. For all details, proofs, and further information, we refer to the many excellent books on the subject, e.g., [1-14]. Here, we follow the style of [1], but the Lie symmetry description is made concise by emphasizing the significant points and results. In order to provide the reader with a relatively quick and painless introduction to Lie symmetry theory, some important concepts must be introduced. ThemaintoolusedinLiegrouptheoryandworkingwithtransformationgroupsis infinitesimal transformation. In order to present this, we need first to develop the concept of a vector field on a manifold. We begin with a discussion of tangent vectors. Suppose C is a smooth curve on a manifold M, parameterized by φ : I M, (3:3) where I is a subinterval of R. In local co-ordinates x = x 1,...,x p, C is given by p smoothing functions φ (ε) = ( φ 1 (ε),..., φ p (ε) ), (3:4) of the real variable ε. At each point x = j(ε) of C the curve has a tangent vector, namely the derivative. φ = dφ ( dφ 1 ) dε = dφp,...,. (3:5) dε dε

5 Page 5 of 16 In order to distinguish between tangent vectors and local coordinate expressions for a point on the manifold, we adopt the notation V = dϕ dε = dϕ1 dε dϕp x1 dε x p (3:6) for the vector tangential to C at x = j(ε) The collection of all tangent vectors to all possible curves passing through a given point x in M is called the tangent space to M at x, and is denoted by TM. AvectorfieldV on M assigns a tangent vector V Î TM to each point x Î M, withv varying smoothly from point to point. In local coordinates, a vector field has the form V = ξ 1 (x) x ξ p (x) x p. (3:7) where each ζ i (x) is a smoothing function of x. If V is a vector field, we denote the parameterized maximal integral curve passing through x in M by Ψ(ε, x) and call Ψ the flow generated by V. Thus for each x in M, and ε in some interval I x containing 0, Ψ(ε, x) is a point on the integral curve passing through x in M. The flow of a vector field has the basic properties: (δ, (ε, x)) = (δ + ε, x), (3:8) for all δ, ε Î R such that both sides of equation are defined, and (0, x) = x, (3:9) d (ε, x) = V (3:10) dε for all ε where defined. We see that the flow generated by a vector field is the same as a local group action of the Lie group on the manifold M, often called a one parameter group of transformations. The vector field V is called the infinitesimal generator of the action since by Taylor s theorem, in local coordinates (ε, x) = x + εξ (x) + O ( ε ), (3:11) where ζ =(ζ 1,..., ζ p )arethecoefficientsofv. The orbits of the one-parameter group action are the maximal integral curves of the vector field V. Definition 1: A symmetry group of Equation 3.1 is a one-parameter group of transformations G, actingonx U, such that if u = f(x) is an arbitrary solution of (3.1) and g ε Î G then g ε f(x) is also a solution of (3.1). The infinitesimal generator of a symmetry group is called an infinitesimal symmetry. Infinitesimal generators are used to formulate the conditions for a group G to make it a symmetry group. Working with infinitesimal generators is simple. First, we define a prolongation of a vector field. The symmetry group of a system of differential equations is the largest local group of transformations acting on the independent and dependent variables of the system such that it can transform one system solution to another. The main goal of Lie symmetry theory is to determine a useful, systematic, computational method that explicitly determines the symmetry group of any given

6 Page 6 of 16 system of differential equations. The search for the symmetry algebra L of a system of differential equations is best formulated in terms of vector fields acting on the space X U of independent and dependent variables. The vector field tells us how the variables x, u transform. We also need to know how the derivatives, that is u x, u xx,..., transform. This is given by the prolongation of the vector field V. Combining these, we have [[1], p. 110, Theorem.36]. Theorem 1 Let V = p q ξ i (x, u) x i + η a (x, u) u a i=1 a=1 be a vector defined on an open subset M X U. Thenth prolongation of the original vector filed is the vector field: pr (n) V = V + q ηa J ( x, u (n) ) a J a=1 J defined on the corresponding jet space M (n) X U (n). The second summation here is over all (unordered) multi-indices J =(j 1, j,...,j k ), with 1 j k p, 1 k n,. The coefficient functions φa J of pr(n) V are given by the following formula: ηa J ( x, u (n) ) ( ) p p = D J η a ξ i u a i + ξ i u a J,i i=1 i=1 where u a i = a x i,and u a J,i = a J x i,andd J are the total derivative of h with respect to x j. In the following analysis, we only deal with one-dimensional scalar differential equations that are assumed to be differentiable up to the necessary order. Consider the special case, where p =,q = 1 in the prolongation formula, so that we are looking at a partial differential equation involving the function u = f(x, t). A general vector field on X U R R then takes the form [[1], p. 114] V = ξ (x, t, u) x + τ (x, t, u) t + η (x, t, u) (3:1) The first prolongation of V is the vector field: where pr (1) V = V + [η x ] + [η t ] (3:13) x t and [η x ] = η x +(η u ξ x )u x τ x u t ξ u u x τ uu x u t [η t ] = η t +(η u τ t )u t ξ t u x τ u u t ξ u u x u t

7 Page 7 of 16 The subscripts on h, ζ, τ denote partial derivatives. Similarly, where pr () V = pr (1) V + [η xx ] xx + [η xt ] xt + [η tt ] tt (3:14) [η xx ] = η xx +(η xu ξ xx )u x τ xx u t +(η uu ξ xu )u x τ xuu x u t ξ uu u 3 x τ uu u x u t +(η u ξ x )u xx τ x u xt 3ξ u u xx u x τ u u xx u t τ u u xt u t [η xt ] = η xt +(η xu τ tx )u t +(η tu ξ tx )u x τ xu u t +(η uu ξ xu τ ut )u x u t ξ tu u x τ uuu x u t ξ uu u t u x τ xu tt +(η u ξ x τ t )u xt ξ t u xx τ u u t u xt ξ u u x u xt τ u u x u tt ξ u u t u xx [η tt ] = η tt +(η tu τ tt )u t ξ tt u x +(η uu τ tu )u t ξ tu u x u t τ uu u 3 x τ uu u t u x +(η u τ t )u tt ξ t u xt 3τ u u tt u t ξ u u tt u x ξ u u xt u t From here on analysis of difference equations only concerns modified equations, which have third prolongation of the vector field. From work in CFD, we know that the right-hand side of the modified equation is written entirely in terms of x derivatives. So, investigation can be limited to the terms of the spatial derivatives in the following analysis. The coefficients of the various monomials in the third-order partial derivatives of u are given in the following: where, pr (3) V = pr () V + [η xxx ] xxx + [η xxt ] xxt + [η xtt ] xtt + [η ttt ] ttt (3:15) [η xxx ] = η xxx +(3η xxu ξ xxx )u x τ xxx u t +3(η xuu ξ xxu )u x 3τ xxu u x u t +(η uuu 3ξ xuu )(u x ) 3 +3(η xu ξ xx )u xx 3τ xx u xt 3τ xuu (u x ) u t +3(η uu 3ξ xu )u x u xx 3τ xu u t u xx 6τ xu u xt u x 3τ x u xxt +(η u 3ξ x )u xxx ξ xxx (u x ) 4 6ξ uu (u x ) u xx 3τ uu (u x ) u xt τ uuu (u x ) 3 u t 3ξ u (u xx ) 3τ u u xxt u x 3τ u u xt u xx 3τ uu u xx u x u t 4ξ u u xxx u x τ u u xxx u t Suppose we are given an nth order system of differential equations, or, equivalently, a subvariety of the jet space M (n) X U (n). A symmetry group of this system is a local transformation G acting on M X U. which transforms solutions of the system to other solutions. We can reduce the important infinitesimals condition for a group G to be a symmetry group of a given system of differential equations. The following theorem [[1], p. 104, Theorem.31] provides the infinitesimal conditions for a group G to be a symmetry group. Theorem Suppose ν ( x, u (n) ) =0,ν = 1,,..., l is a system of differential equations of maximal rank defined over M X U. IfG is a local group of transformations acting on M, and

8 Page 8 of 16 pr (n) V [ ν ( x, u (n) )] =0,ν = 1,,..., l whenever ν ( x, u (n) ) =0, for every infinitesimal generator V of G, then G is a symmetry group of the system. In the following sections, this theorem is used to deduce explicitly different infinitesimal conditions for specific problems. It must be remembered, however, that, in all cases, though only the scalar differential problem is being discussed, Δ v is still used to denote different differential equations. 4. Galilean group and its prolongation It is well known that as for Navier-Stokes equations, the intrinsic symmetries, except for the scaling symmetries, are just macroscopic consequences of the basic symmetries of Newton s equations governing microscopic molecular motion (in classical approximation). Any physical difference scheme should inherit the elementary symmetries (at least for Galilean symmetry) from the Navier-Stokes equations. This means that Galilean invariance has been an important issue in CFD. Furthermore, we stress that Galilean invariance is a basic requirement that is demanded for any physical difference scheme. We have the Galilean transformation x = x + tε t = t (4:1) u = u + ε Thus, the vector of the Galilean transformation is V = t x + u (4:) According to Theorem 1, we have pr (1) V = V ρ x u x p x (4:3) ρ t t p t pr () V = pr (1) V ρ xx u xx p xx (4:4) ρ xt xt p xt 5. Galilean invariance of first-order conservative form scheme The main prototype equation here is the modified equation. Equation.18 can be recast into 1 u t + uu x 1 tβ (u, λ) u xx 1 tβ uu x u x =0. (5:1) Based on the prolongation formula presented in Section 4, the Galilean invariance condition reads 1 =0. (5:)

9 Page 9 of 16 pr () V 1 =0. (5:3) Before beginning the group analysis, some detailed but mechanical calculations must be performed: d 1 = u 1 = u x 1 tβ uu xx 1 tβ uuu x u x. (5:4) d = ut 1 =1. (5:5) With these formulas, it is clear from Equation 5.3 that the invariance condition reduces into β u u xx + β uu u x u x =0. (5:6) Hence, we have u x u x = β u β uu u xx. (5:7) we can then write the model equation as u t + uu x 1 [ t β β ] uβ u u xx =0. (5:8) β uu with [ 1 t β β ] uβ u = ν 1. (5:9) β uu This manipulation yields the Burgers equation as following u t + uu x = ν 1 u xx. (5:10) where v 1 = constant. Based on the analysis of Equation 5.9, one have β = β 0 exp (αu) + ν 1 t. (5:11) where b 0, a are some parameters. Here, it is useful to list some well-known first-order conservative schemes to show their unified character.

10 Page 10 of Lax-Friedrichs scheme f n 1 = 1 [ ( ) f u n j+1 + f j+ ( ) u n j 1 ( ) ] u n j+1 λ un j. (5:1) G 1 = 1 G 1 = 1 [ 1 λ f ], (5:13) [ 1+λ f ]. (5:14) β (u, λ) = 1 ( ) f λ = 1 λ u. (5:15) point monotonicity scheme (Godunov, 1959) K u n+1 j = C l u n j+r, (5:16) r= K G r = C r (5:17) β (u, λ) = 1 λ K r= K ( ) f r G r = 1 K λ r C r u. r= K (5:18) 5.3. General 3-point conservation scheme f ( ) 1 u j, u j+1 = f ( ) ( ) 1 u j + f uj+1 λ Q λā 1 1 u, (5:19) j+ j+ where ā j+ 1 = f ( u ) ( ) j+1 f uj 1 u j+, when 1 u 0, j+ ā 1 = a ( ) u j, when 1 u =0, j+ j+ Here Q(x) is some function, which is often referred to as the coefficient of numerical viscosity. Harten s lemma. Let Q(x) in (5.19) satisfy the inequalities x Q (x) 1for0 x μ 1;

11 Page 11 of 16 then finite-difference scheme is TVNI under the CFL-like restriction λmax ā j 1 μ. j+ The coefficient of numerical viscosity could be expressed in terms of the b as follows β (u, λ) = 1 Therefore, one can have [ Q (u) u ]. (5:0) Q (u) = u +β 0 exp (αu) + 4ν 1 t. (5:1) If we choose β 0 << 1 (5:) Then Equation 5.1 is consistence with the results of Harten s. In summary, we obtain Theorem 3 If we let the coefficients in (.18) satisfy the equality β = β 0 exp (αu) + ν 1 t. where b 0, a is a dimensionless constant, then the first-order conservative finite different scheme satisfies Galilean invariant condition. 6. Galilean invariance of second-order conservative scheme The same manipulation could be conducted for the case of the second-order conservative scheme. The main prototype equation here is [15] 1 u t + uu x = 1 6 t3 x [γ (u, λ) u xx + δ (u, λ) u x u x ]. (6:1) where γ (u, λ) = 1 λ 3 δ (u, λ) = 1 λ 3 K r= K K r,s= K According to Theorem, one have ( ) f 3 r 3 G r +, (6:) 1 4 (r + s) [ ( ) rs (r s) ] f ( ) f G rs +3. (6:3) u t + uu x 1 6 t3 [ γ u xxx + (γ u +δ) u x u xx + δ u (u x ) 3]. (6:4) Before beginning the group analysis, some detailed but mechanical calculations must be performed: d 1 = u = u x 1 6 t3 [ γ u u xxx + (γ uu +δ u ) u x u xx + δ uu (u x ) 3]. (6:5)

12 Page 1 of 16 d = ut =1. (6:6) pr (3) V = t x + u u x ut. (6:7) The corresponding Galilean invariant condition reads: =0. (6:8) pr (3) V =0. (6:9) The substitution leads u x 1 6 t3 [ γ u u xxx + (γ uu +δ u ) u x u xx + δ uu (u x ) 3] u x =0 (6:10) Hence, we have γ u u xxx + (γ uu +δ u ) u x u xx + δ uu (u x ) 3 =0. (6:11) It is clear that u x u xx = δ uu (u x ) 3 γ u u xxx. (6:1) γ uu +δ u γ uu +δ u After some manipulation, one could obtain the model equation as follows u t + uu x 1 ( 6 x3 γ γ ) u +δ γ u.u xxx γ uu +δ u 1 ( 6 x3 δ u γ ) u +δ δ uu (u x ) 3. γ uu +δ u (6:13) In order to obtain the non-oscillation solution of shock, we could let the term of (u x ) 3 to be zero, then we have δ u γ u +δ γ uu +δ u δ uu =0. (6:14) This equation can be rewritten as below δ uu δ u = γ uu +δ u γ u +δ (6:15) d [ ] d [ ln δu = ln (γu +δ) ]. (6:16) du du δ u = γ u +δ. (6:17) or γ u = δ u +δ. (6:18)

13 Page 13 of 16 The coefficient of the term of u xxx could be rewritten as If we set γ γ u +δ γ uu +δ u γ u = γ δ u δ uu γ u = γ δ u δ uu (δ u δ). (6:19) δ = au n + b (6:0) where a, b, m are some parameters. It is easy to show that γ = an n 1 un a n 1 un+1 where g 0 = constant. In summary, we obtain Theorem 4 If we let the coefficients in (6.1) satisfy the equality δ (u, λ) = au n + b, γ (u, λ) = an n 1 un a n 1 un+1 b n 1 u + γ 0. (6:1) b n 1 u + γ 0. where a, b, n, g 0 is a dimensionless constant, then the second-order conservative finite different scheme satisfies Galilean invariant condition. Here,wecouldgivethedetailsofthecorresponding analysis by using the Lax- Wendroff scheme. It is well known that the Lax-Wendroff difference approximation to (.1) is defined by u n+1 j t u n j + f (un j+1 ) f (un j 1 ) x = 1 t a j+ 1 (u n j+1 un j ) a j 1 (u n j u n j 1 ) x (6:) where u n j = u ( n t, j x ) a (u) = f (u) a j+ 1 u n j+ 1 = a(u n j+ 1 ) = 1 (un j+1 + un j ) It is well known that the Lax-Wendroff scheme was designed to have the following desirable computational features [16-18]: conservation of form; to have a three-point

14 Page 14 of 16 scheme; second-order accuracy on smooth solutions. Numerical spikes and downstream oscillations are generated in the vicinity of the shock. The numerical flux for Lax-Wendroff scheme can be written as f n j+ 1 = 1 f ( ) ( ) u n j+1 + f u n j t ( x a 1 u n j+1 un j j+ ). (6:3) Using the general method presented in Section, one can obtain G 1 = 1 λ f + 1 ( ) f λ (6:4) G 1 = 1 λ f + 1 ( ) f λ (6:5) G 11 = λ ( 1 λ f ) f, (6:6) G 00 = λ f f, (6:7) G 1, 1 = λ ( 1+λ f ) f. (6:8) and for the other case of the index r, s G rs =0. (6:9) γ LW (u, λ) = 1 λ ( δ LW (u, λ) = 1 λ ( 1 λ ( f 1 3λ ( f ) ) ) ) f, (6:30) f, (6:31) γ LW (u, λ) = 1 λ u + u3, (6:3) δ LW (u, λ) = 1 λ +3u. (6:33) Using the results of Theorem 4, we have δ G (u, λ) = au n + b = δ LW. (6:34) This leads the following corresponding relation for Lax-Wendroff scheme n =, (6:35)

15 Page 15 of 16 a =3,b = 1 λ. (6:36) γ G (u, λ) = an n 1 un a n 1 un+1 =6u 6u 3 + λ u + γ 0. b n 1 u + γ 0 (6:37) The comparison of Equations 6.3 and 6.37 help us to draw some conclusion as follows: based on this analysis, we have known that the well-known Lax-Wendroff scheme can recover the Galilean symmetry approximately. 7. Conclusions It is known that the numerical solutions calculated by finite difference schemes are always associated with numerical dissipation and dispersion. Such errors can lead to undesirable numerical effects, especially for shock capturing. A full understanding of the nature of this odd numerical phenomenon is still lacking. Regardless of definition, spurious oscillations and overshoots are the most common symptoms of numerical stability. In one natural interpretation, these numerical phenomena are due to nonlinear stability, which links some symmetry breaking. This article uses a Lie symmetry analysis method to investigate Galilean invariance properties of several difference schemes for nonlinear equations. Two theorems have been obtained, which have demonstrated that the properties of Galilean invariance from a modification equation, can serve as the positive constrain condition for general conservative finite difference schemes. It should be pointed out that the conclusions presented in this article have preliminary character and demand further study. Acknowledgements The author appreciates the numerous corrections and very useful comments from the referees on the manuscript. The study was supported by the National Natural Science Foundation of China (Grant Nos , ), and also supported by the State Key Laboratory for Turbulence and Complex System. Competing interests The authors declare that they have no competing interests. Received: 4 November 010 Accepted: 7 November 011 Published: 7 November 011 References 1. Harlow, FH: The particle-in-cell computing method for fluid dynamics. Methods Comput Phys. 3, 319 (1964). Yanenko, NN, Shokin, YI: On the group classification of difference schemes for systems of equations in gas dynamics. Proceedings of the second International Conference on Numerical Methods in Fluid Dynamics. (1970) 3. Winternitz, P: Symmetries of discrete system. Lecture presented at the CIMPA Winter School on Discrete Integrable System. (003) 4. Ran, Z: Lie symmetries preservation and shock capturing methods. SIAM J Numer Anal. 46(1), (007) 5. Ran, Z: Symmetry and non-physical oscillation in shock capturing. Chin Quart Mech 6, 4 (005). (In Chinese) 6. Harten, A: High resolution schemes for hyperbolic conservation laws. J Comput Phys. 49, (1983). doi: / (83) Harten, A, Lax, PD: A random choice finite difference scheme for hyperbolic conservation laws. SIAM J Numer Anal. 18, (1981). doi: / Oleinik, OA: Discontinuous solutions of non-linear differential equations. Usp Mat Nauk 1((3(75))):3 73 (1957). (in Russian) 9. Harten, A, Lax, PD, Leer, BV: On upstream difference and Godunov-type schemes for hyperbolic conservation laws. SIAM Rev. 5(1), (1983). doi: / Harten, A, Hyman, JM, Lax, PD: On finite difference approximations and entropy conditions for shocks. Commun. Pure Appl Math. 9, 97 3 (1976). doi:10.100/cpa Osher, S, Chakravarthy, S: High resolution schemes and the entropy condition. SIAM J Numer Anal. 1(5), (1984). doi: / Olver, PJ: Applications of Lie Group to Differential Equations. GTM 107. Springer-Verlag, New York (1986)

16 Page 16 of Bluman, GW, Kumei, S: Symmetries and Differential Equations. In Applied Mathematical Sciences, vol. 81, pp Springer-Verlag, World Publishing Corp (1991) 14. Marsden, JE, Ratiu, TS: Introduction to Mechanics and Symmetry: A Basic Exposition of Classical Mechanical Systems. In Texts Appl Math, vol. 17,Springer Verlag, New York (1994) 15. Shui, HS: Finite Difference Methods for One Dimensional Fluid Dynamics. National Defence Industrial Press. Beijing (1998) (in Chinese) 16. Lax, PD, Wendroff, B: Systems of conservation laws. Commun Pure Appl Math. 3, (1960) 17. Majda, A, Osher, S: A systematic approach for scalar conservation laws. Numer Math. 30, (1978). doi: / BF Rubin, EL, Burstein, SZ: Difference methods for inviscid and viscous equations of compressible gas. J Comput Phys. l, (1976) doi: / Cite this article as: Ran: Galilean invariance and the conservative difference schemes for scalar laws. Advances in Difference Equations :53. Submit your manuscript to a journal and benefit from: 7 Convenient online submission 7 Rigorous peer review 7 Immediate publication on acceptance 7 Open access: articles freely available online 7 High visibility within the field 7 Retaining the copyright to your article Submit your next manuscript at 7 springeropen.com