Geometrical and Numerical Approach to Solve Transonic Gas Equation

In this paper, the transonic gas equation will be considered. Different methods will be presented and discussed. We will show that there is no unique solution in the computational domain [ ] 1,1 − with Dirichlet boundary conditions.


Introduction
Consider the transonic gas equation The study of gas dynamics is often associated with the flight of modern high speed aircraft and atmospheric reentry of space exploration vehicles; however, its origins lie with simpler machines.At the beginning of the 19th century, investigation into the behavior of fired bullets led to improvement in the accuracy and capabilities of guns and artillery.As the century progressed, inventors such as Gustaf de Laval advanced the field, while researchers such as Ernst Mach sought to understand the physical phenomenon involved through experimentation the model equation gas equation.Lie groups lie at the intersection of two fundamental fields of mathematics: algebra and geometry.A Lie group is first of all a group.Secondly it is a smooth manifold which is a specific kind of geometric object.The circle and the sphere are examples of smooth manifolds.Finally the algebraic structure and the geometric structure must be compatible in a precise way.Informally, a Lie group is a group of symmetries where the symmetries A circle has a continuous group of symmetries, we can rotate the circle an arbitrarily small amount and it looks the same.This is in contrast to the hexagon, for example.If you rotate the hexagon by a small amount then it will look different.Only rotations that are multiples of one-sixth of a full turn are symmetries of a hexagon.Lie groups were studied by the Norwegian mathematician Sophus Lie at the end of the 19th century.Lie was interested in solving equations.At that time techniques for solving equations were basically a bag of tricks.A typical tool was to make a clever change of variables which would make one of the variables drop out of the equations.Lie's basic insight was that when this happened it was due to an underlying symmetric of the equations, and that underlying this symmetry was what is now called a Lie group.Lie groups are ubiquitous in mathematics and all areas of science.Associated to any system which has a continuous group of symmetries is a Lie group.The symmetry Lie groups and invariant solutions for the wave diffusion in gas and the transonic equation were obtained by Ames, R. J. Lohner [1].After their work, Torissi and Valenti [2] applied the infinitesimal group analysis to a second-order nonlinear wave equation telling a non-homogeneous manner.Bluman and Kumei [3] deliberate the wave equations with two and four-parameter symmetry groups and the corresponding invariant (similarity) solutions.Peters and Ames [4] applied Lie group method to the nonlinear dynamic equations of elastic string.Group invariant solutions for a second-order hyperbolic wave equation including an axially transported speed and nonlinear terms were constructed by Fung et al. [5].Ozkaya and Pakdemirli [6] accurate exact solutions of accidental vibrations of a fiber moving with time-dependent velocity.Simple examples of Lie groups contain translation, rotation and scaling.Lie groups and infinitesimal generators can be logically prolonged to act on the space of independent variables [7] [8] [9].Sami H. Altoum, studied Lie group and RK4 for solving nonlinear first order ODEs [10].This article is organized as follows: In Section 2 we state the main problem.In Section 3 introduced some preliminary definitions related to Lie group.In Section 4 A Lie group method approach is purposed to obtain exact solution.In Section 5 we recall a solution obtained see [10].In Section 6 a pdsolve maple method is proposed as technique to get an analytical solution.In Section 7 a central finite difference method proposed to solve transonic gas equation.In Section 8 a high order finite difference method is proposed as technique to show that there is no unique solution of transonic gas equation.Finally we conclude.

Problem
Solve we focus in our study to the homogeneous Dirichlet boundary condition on ∂Ω , when ∂Ω denotes the boundary of the domain Ω .For simplicity we will use homogeneous Dirichlet boundary condition.Up to now there is now way to solve equation manually.So that Lie group technique will be introduced and used as tools to compute a nonzero solution of (2).Furthermore, numerical approach based on central finite difference method will be introduced and a high order accurate finite difference method will be used to show that no unique solution can be found.

Definition of a Group
Consider invertible transformations in an n-dimensional Euclidean space n  defined by equations of the form ( ) It is assumed that the vector function is continuous, together with its derivatives.It is also assumed that the coordinates i x and i x of points x and x are referred to the same coordinate system.Since the above transformation is invertible, there exists the inverse transformation: We denote T to above transformation and 1 T − the inverse transformation.Accordingly, T carries any point identically in x and a in a neighborhood of 0 a = .

Infinitesimal Transformation
Given 1) parameter group G of transformation ( ) Let us expand ( ) x a into Taylor series in the parameter (a) at 0 a = , then invoking the initial condition where Geometrically, the infinitesimal transformation (4) and (5) defines the tangent vector therefore the tangent vector field of the group G denoted by i ξ and written as a first order li- near differential operator ( ) ( )

System of Differential Equation
Consider a system of equations in n  ( ) 0, 1, , , F x s where x in n  and s n < .The system ( 7) is invariant with respect to a group G of transformations (3) [11].Theorem 1 (cf [11] [12]) The system (7) invariant under group G with infinitesimal generator X if and only if where the symbol M F σ means the value F σ of on the manifold M. The system of ( 7) is said to be invariant with respect to group G of transformation ( ) , or that Equation ( 6) admits G as a symmetry group if  , whenever x solve Equation (7).Geometrically it means that transformations of the group G carry any point of the variety M along this variety.In other words, 1, , s σ =  , the path curve of the group G passing through any point x M ∈ lies on M. Consequently, M is an invariant manifold for G.

Determining Equations
Theorem (1) provides the infinitesimal criterion of symmetry group of differential equations.The system of DE is, ( ) is invariant under the group with an infinitesimal generator X if an only if where X is extended to all derivatives involved in (9).We give three methods to solve transonic gas equation to get exact solutions, the first method using Lie groups, second method using classical solution of PDEs and third method using finite difference.

Prolongation
We need to supplant the slightly hazy notion of a system of D.E. by stiff geometric target specified by evanescence of particular functions.To do this we need to prolong the basic space, representing the independent and dependent variables under consideration, to a space which also represents the various partial derivatives occurring in the system [13].
( ) ( ) riables, there are given by: order partial derivatives of f.We derive k th employ the multi-index notation ; , , x x x y = and U ≈  has the single coordinate u.The space 1 U is isomorphic to 2  with coordinates ( ) x y u u , since these represent all the first order partial derivatives of u with respect to x and y.Similarly Here, note that there are two independent variables , x t , one dependent variable, u, so we have the underlying space

3
X U × =  with coordinates , , x u t .Since the heat equation contains a second order partial derivative, we consider the prolongation of the 2-jet, namely to the set X U × which has coordinates ( ) , , , , , , , x t u u u u u u .

Lie Algebra of Operator
A Lie algebras of operators, is a vector space L of operator ( ) ∂ with the following property.If the operators ( ) ( ) is bilinear and skew-symmetric also an element of L. The dimension (dim L) of the Lie algebra is the dimension of the vector space L. We use the symbol r L to denote an r-dimensional Lie algebra.

Main Variables and Total Derivatives
In differential algebra we deal with an infinite derivatives  are systematic in the subscripts 1 2 1 2 3 , , i i i i i  .The main operation in the calculus of differential algebra is the total differentiation given by the following The total derivatives i D act on functions involving any finite number of va- riables.
( ) ( ) ( ) In particular, letting , , Consequently, i x are called independent variables and u α differentiable variables with the successive ( ) ( ) , , u u  Hence the total derivatives written by: , We will use a generator of a point transformation group Extended Equation ( 14) then ( ) ( ) and also we can deduce

Lie Group
We introduce a geometrical method based on Lie group technique to solve transonic gas equation, namely problem (2).
The same argument applied to the terms containing xx u yields & (32).The computation shows that the general solution of this system is given by: Finally we can deduce the solution u from Equation (31) ( ) ( ) ( ) ) we can write u after simplification ( ) So we have six arbitrary constants i C .Thus we arrive at the six-dimensional Lie algebra 6 L spanned by the operators.
Thus, 6 L is the linear span of the operators (33),  , of the basic operators and satisfy (12) for a Lie algebra.The result becomes directly visual if the commutators are disposed as in the following Table 1.

0.
x xx yy The solution will be: 1) Let ( ) , u x t u is the solution of (2).Then the function.
( ) u c c u c x c c y c c y c are arbitrary constant, is also a solution.
2) Solutions ( ) ) x c u x y c y c a x c ) ) 12 x c x c aA u x y A c y c y c y c a y c 3

Pdsolve Maple Method
The first and second method proposed to study more precisely to find the solution of transonic gas equation., C C and are global general solution of the second order PDE (transonic gas equation) is ( ) The following is the solution computed by pdsolve automatically.The approach used in this example serves also to address nonlinear equations of the form (2) with homogenous boundary conditions.This document describes how pdsolve can automatically adjust the arbitrary functions and constants entering the solution.In order, to determine the coefficient involved in (35) we will introduce some boundary conditions on ∂Ω , where Ω is the domain of com- putation.At first step: Let us, use Dirichlet homogeneous boundary condition and [ ] 2 0,1 Ω = . The system of equation obtained after applying the boundary conditions is: this system is solved using Maple 18, we get the following result:  Trivial Solution: in case where Figure 1 presents the shape of curve of transonic gas equation.

Central Finite Difference Method
Let us consider the two-dimensional nonlinear PDE equations:

u x u x u y u y = = = =
The finite difference schemes take the form   Figure 3 presents the interior points of solution of transonic gas equation.

High Finite Difference Method
The key step in solving our PDE numerically using finite difference methods is to replace the derivatives with so-called "finite difference Method".Here is an example of the centered finite difference stencil for the second derivative: , , The stencil is "centered" because u is evaluated at an equal number of points to the right and left of the point where we want to approximate the derivative.
The coefficients i β are specified such that error is ( ) 6 O h .That is, we solve for error, expand in a Taylor series in h and then we choose the coefficients such that the first 5 terms in the Taylor series vanish, i.e. Figure 3. Interior points solution.

h y u x h y u x y u x h y u x h y h h h h h u x h y u x y u x y h u x y h h h h h u x y u x y h u x y h h h h
where ( ) Since, the number of equations is less the number of unknowns parameters, for this reason we cannot have a unique solution for our differential equation.This leads to an open question, what can we do in our case or how can we change completely the PDEs in order to get a solvable system with unique solution.Our code, will be included in the appendix, and it's available for any boundary conditions that can be introduced.

Conclusion
In this paper different methods were discussed.More precisely, we have considered a Lie group approach to get an analytical solution; this solution can be viewed as old approach that gives us a solution.To understand our problem we introduced different numerical schemes were proposed and we arrive to no unique solution that can be founded. and

1 T
− returns x into the original position x.The identical transformation, unaltered when one moves along any path curve of the group G.In other words, F is an invariant if ( )

1 k
− order partial derivatives of u.We have the heat equation 0 η depend only upon x, y and u.Accordingly we isolate the terms containing u and those free of these variables, and set each term equal to zero.So, the terms containing xy u , lead to following equa- 29) is split into the over determined system of linear PDEs (30) The pdsolve routine usually introduces new functions to express the solution for the indeterminate function.The examples below illustrate the use of pdsolve in solving a single PDE.For examples related to solve the transonic gas equation.Any arbitrary constants introduced while separating the variables are represented as 1 2

Figure 2
Figure 2 presents central finite difference solution of transonic gas equation.