Geometrical and Numerical Approach to Solve Transonic Gas Equation

Abstract

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]2 with Dirichlet boundary conditions.

Share and Cite:

Altoum, S. (2018) Geometrical and Numerical Approach to Solve Transonic Gas Equation. Journal of Applied Mathematics and Physics, 6, 1659-1674. doi: 10.4236/jamp.2018.68142.

1. Introduction

Consider the transonic gas equation

u x u x x + u y y = 0 (1)

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 are continuous. 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.

2. Problem

Solve

u x u x x + u y y = 0 in Ω = [ 1 , 1 ] 2 (2)

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.

3. Preliminary

3.1. Definition of a Group

Consider invertible transformations in an n-dimensional Euclidean space n defined by equations of the form

x ¯ i = f i ( x ) , i = 1 , , n ,

where x = ( x 1 , , x n ) n and x ¯ = ( x ¯ 1 , , x ¯ n ) n . It is assumed that the vector function f = ( f 1 , , f n ) is continuous, together with its derivatives. It is also assumed that the coordinates x i and x ¯ i of points x and x ¯ are referred to the same coordinate system. Since the above transformation is invertible, there exists the inverse transformation:

x i = ( f 1 ) i ( x ¯ ) , i = 1 , , n .

We denote T to above transformation and T 1 the inverse transformation. Accordingly, T carries any point x n into a new position x ¯ n , and T 1 returns x ¯ into the original position x. The identical transformation, x ¯ i = x i , i = 1 , , n denoted by I. A function F ( x ) is called an invariant of group G of transformations x ¯ i = f i ( x , a ) , i = 1 , , n if F remains unaltered when one moves along any path curve of the group G. In other words, F is an invariant if F ( f ( x , a ) ) = F ( x ) identically in x and a in a neighborhood of a = 0 .

3.2. Infinitesimal Transformation

Given 1) parameter group G of transformation

x ¯ i = f i ( x , a ) , i = 1 , , n (3)

Let us expand f i ( x , a ) into Taylor series in the parameter (a) at a = 0 , then invoking the initial condition ( f i ( x , 0 ) = x i , i = 1 , , n ) , we arrive to infinitesimal transformation of G [10] .

x ¯ i = x i + a ξ i ( x ) = f i ( x , 0 ) + ( a 0 ) f ( x , 0 ) a + O ( a 2 ) , i = 1 , , n (4)

where

ξ i ( x ) = a [ f i ( x , a ) ] a = 0 (5)

Geometrically, the infinitesimal transformation (4) and (5) defines the tangent vector ξ i ( x ) = ( ξ 1 ( x ) , , ξ n ( x ) ) at the point ( 0, f ( x ,0 ) ) therefore the tangent vector field of the group G denoted by ξ i and written as a first order linear differential operator

X = ξ 1 ( x ) x 1 + + ξ n ( x ) x n

X = ξ i ( x ) x i , i = 1 , , n (6)

3.3. System of Differential Equation

Consider a system of equations in n

F σ ( x ) = 0 , σ = 1 , , s , (7)

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

X F σ ( x ) | ( M ) = 0 , σ = 1 , , s (8)

where the symbol F σ | M 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 x ¯ i = f i ( x , a ) , or that Equation (6) admits G as a symmetry group if F σ ( x ) = 0 , σ = 1 , , s , 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.

3.4. Determining Equations

Theorem (1) provides the infinitesimal criterion of symmetry group of differential equations. The system of DE is,

F σ ( x , u , u 1 , , u k ) = 0 , σ = 1 , , s (9)

is invariant under the group with an infinitesimal generator X if an only if

X F ( σ ) | M = 0 , σ = 1 , , s (10)

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.

3.5. 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] . f ( x ) = f ( x 1 , , x p ) of p independent variables, there are given by:

p k = [ k p + k + 1 ] ,

order partial derivatives of f. We derive kth employ the multi-index notation

k f x k = k f x 1 j 1 x 2 j 2 x n j n ; j 1 + j 2 + + j n = k

Consider the case P = 2 , q = 1 then X 2 has coordinates ( x 1 , x 2 ) = ( x , y ) and U has the single coordinate u. The space U 1 is isomorphic to 2 with coordinates ( u x , u y ) , since these represent all the first order partial derivatives of u with respect to x and y. Similarly U 2 3 has the coordinates ( u x x , u x y , u y y ) , representing the second order partial derivatives of u, and in general U k k + 1 , since there are k + 1 distinct k 1 order partial derivatives of u. We have the heat equation u t u x x = 0 . Here, note that there are two independent variables x , t , one dependent variable, u, so we have the underlying space X × U = 3 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 x , u t , u x x , u x t , u t t ) .

3.6. Lie Algebra of Operator

A Lie algebras of operators, is a vector space L of operator X = ξ i ( x ) x i with the following property. If the operators

X 1 = ξ 1 i ( x ) x i , X 2 = ξ 2 i ( x ) x i

are elements of L then their commutator is

[ X 1 , X 2 ] X 1 X 2 X 2 X 1 = ( X 1 ( ξ 2 i ) X 2 ( ξ 1 i ) ) x i (11)

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 L r to denote an r-dimensional Lie algebra. X L r is their linear combination, X = c α X α = α = 1 r c α X α with constant coefficients c α . In particular, [ X α , X β ] L r , hence

[ X α , X β ] = c α β γ X γ , α , β = 1 , , r (12)

3.7. Main Variables and Total Derivatives

In differential algebra we deal with an infinite derivatives

x = { x i } , u = { u α } , u ( 1 ) = { u i α } , u ( 2 ) = { u i 1 i 2 α } , u ( 3 ) = { u i 1 i 2 i 3 α } ,

where α = 1 , , m ; i , i 1 , = 1 , , n . The variables u i 1 i 2 α , u i 1 i 2 i 3 α , are systematic in the subscripts i 1 i 2 , i 1 i 2 i 3 , . The main operation in the calculus of differential algebra is the total differentiation given by the following

D i = x i + u i α u α + u i i 1 α u i 1 α + u i i 1 i 2 α u i 1 i 2 α + , i = 1 , , n

The total derivatives D i act on functions involving any finite number of variables.

D i ( f ( x , u , u ( 1 ) ) ) = f x i + u i α f u α + u i i 1 α f u i 1 α

In particular, letting f = u α , f = u j α , one obtain:

u i α = D i ( u α ) , u i j α = D i ( u j α ) = D i D j ( u α ) (13)

Consequently, x i are called independent variables and u α differentiable variables with the successive u ( 1 ) , u ( 2 ) ,

Hence the total derivatives written by:

D x = x + u x u + u x x u x + u x y u y + ,

D y = y + u y u + u x y u x + u y y u y +

We will use a generator of a point transformation group

X 2 = ξ 1 ( x , y , u ) x + ξ 2 ( x , y , u ) y + η ( x , y , u ) u . (14)

Extended Equation (14) then X 2 becomes

X 2 = ξ 1 x + ξ 2 y + η u + ς 1 u x + ς 2 u y + ς 11 u x x + ς 12 u x y + ς 22 u y y (15)

ς i α = D i ( η α ) u j α D i ( ξ j ) (16)

and

ς i 1 i 2 = D i 2 ( ς i 1 α ) u j i 1 α D i 2 ( ξ j ) = D i 2 D i 1 ( η α ) u j α D i 2 D i 1 ( ξ j ) u j i 1 D i 2 ( ξ j ) (17)

then

ς 1 = D x ( η ) u x D x ( ξ 1 ) u y D x ( ξ 2 ) (18)

ς 2 = D y ( η ) u x D y ( ξ 1 ) u y D y ( ξ 2 ) (19)

also we can deduce

ς 11 = D x ( ς 1 ) u x x D x ( ξ 1 ) u x y D x ( ξ 2 ) (20)

ς 12 = D y ( ς 1 ) u x x D y ( ξ 1 ) u x y D y ( ξ 2 ) (21)

ς 22 = D y ( ς 2 ) u x y D y ( ξ 1 ) u y y D y ( ξ 2 ) (22)

then

ς 1 = η x + u x η u u x ξ x 1 ( u x ) 2 ξ u 1 u y ξ x 2 u x u y ξ u 2 (23)

ς 2 = η y + u y η u u x ξ y 1 u x u y ξ u 1 u y ξ 2 ( u y ) 2 ξ u 2 (24)

ς 11 = η x x + 2 u x η x u + u x x η u + ( u x ) 2 η u u 2 u x x ξ x 1 u x ξ x x 1 2 ( u x ) 2 ξ x u 1 3 u x u x x ξ u 1 ( u x ) 3 ξ u u 1 2 u x y ξ x 2 u y ξ x x 2 2 u x u y ξ x u 2 ( u y u x x + 2 u x u x y ) ξ u 2 ( u x ) 2 u y ξ u u 2 (25)

ς 12 = η x y + u y η x u + u x η y u + u x y η u + u x u y η u u u x y ( ξ x 1 + ξ y 2 ) u x ξ x y 1 u x x ξ y 1 u x u y ( ξ x u 1 + ξ y u 2 ) ( u x ) 2 ξ y u 1 ( 2 u x u u x y + u y u x x ) ξ u 1 ( u x ) 2 u y ξ u u 1 u y ξ x y 2 u y y ξ x 2 ( u y ) 2 ξ x u 2 ( 2 u y u x y + u x u y y ) ξ u 2 u x ( u y ) 2 ξ u u 2 (26)

ς 22 = η y y + 2 u y η y u + u y y η u + ( u y ) 2 η u u 2 u y y ξ y 2 u y ξ y y 2 2 ( u y ) 2 ξ y u 2 3 u y u y y ξ u 2 ( u y ) 3 ξ u u 2 2 u x y ξ y 1 u x ξ y y 1 2 u x u y ξ y u 1 ( u x u y y + 2 u y u x y ) ξ u 1 u x ( u y ) 2 ξ u u 1 (27)

4. Lie Group

We introduce a geometrical method based on Lie group technique to solve transonic gas equation, namely problem (2).

Considering symmetry operator in the form (15) with unknown coefficients ξ i and η to be found the determining Equation (10)

X F = 0

X ( u x u x x + u y y ) = 0 (28)

u x x ζ 1 + u x ζ 11 + ζ 22 = 0 , (29)

where we substitute ς 1 , ς 11 and ς 22 see (23), (25), (27) from extended generators calculation of symmetries of PDEs, with two independent variables and set u y y = u x u x x . Then Equation (29) contain the variables x , y , u , u x , u y , u x x u y y whereas ξ 1 , ξ 2 and η depend only upon x, y and u. Accordingly we isolate the terms containing u x y , u x x , u x , u y and those free of these variables, and set each term equal to zero. So, the terms containing u x y , lead to following equation

ξ y 1 + u x ξ x 2 + u y ξ u 1 ( u x ) 2 ξ u 2 = 0 ,

whence

ξ y 1 = 0 , ξ u 1 = 0 , ξ x 2 = 0 , ξ u 2 = 0. (30)

The same argument applied to the terms containing u x x yields

η x = 0 , η u 3 ξ x 1 + 2 ξ y 2 = 0 , (31)

then Equation (29) reduces to

η y y = 0 , 2 η y u ξ y y 2 = 0. (32)

Thus Equation (29) is split into the over determined system of linear PDEs (30) & (32). The computation shows that the general solution of this system is given by:

ξ 1 = C 1 x + C 2

ξ 2 = C 3 y + C 4

η = ( 3 C 1 2 C 3 ) u + C 5 y + C 6 .

Finally we can deduce the solution u from Equation (31)

( 3 C 1 2 C 3 ) u + C 5 y + C 6 = 3 ( C 1 x + C 2 ) 2 ( C 3 y + C 4 )

u ( x , y ) = ( C 1 x + C 2 ) 2 ( C 3 y + C 4 ) C 5 y C 6 ( 3 C 1 2 C 3 )

we can write u after simplification

u ( x , y ) = α x + β y + γ

where

α = C 1 3 C 1 2 C 3

β = 2 C 3 + C 5 3 C 1 2 C 3

γ = C 2 3 C 1 2 C 3 + C 4 3 C 1 2 C 3 C 6 3 C 1 2 C 3

So we have six arbitrary constants C i . Thus we arrive at the six-dimensional Lie algebra L 6 spanned by the operators.

X 1 = x , X 2 = y , X 3 = u , X 4 = y u , X 5 = x x + 3 u u , X 6 = y y 2 u u . (33)

Thus, L 6 is the linear span of the operators (33), L 6 = X 1 , , X 6 , from (11), the commutators [ X α , X β ] , 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.

5. Analytical Solution

The exact solution of non-linear PDEs, (2), given by see [10] .

u x u x x + u y y = 0. (34)

Table 1. Lie algebra operators.

The solution will be:

1) Let u ( x , t ) u is the solution of (2). Then the function.

u 1 = c 1 3 c 2 2 u ( c 1 x + c 3 , c 2 y + c 4 ) + c 5 y + c 6 ,

where c 1 , , c 6 are arbitrary constant, is also a solution.

2) Solutions

u ( x , y ) = c 1 x y + c 2 x + c 3 y + c 4

u ( x , y ) = ( x + c 1 ) 2 3 a ( x + c 1 ) 3 + c 3 y + c 4

u ( x , y ) = a 2 c 1 3 39 ( y + A ) 13 + 2 3 a ( y + A ) 8 ( x + B ) + 3 c 1 ( y + A ) 3 ( x + B ) 2 ( x + B ) 3 3 a ( y + A ) 2

u ( x , t ) = a 2 c 1 3 39 ( y + A ) 13 + 2 3 a ( y + A ) 8 ( x + B ) + 3 c 1 ( y + A ) 3 ( x + B ) 2 ( x + B ) 3 3 a ( y + A ) 2

u ( x , y ) = a A 3 y 2 B 2 a A 2 x + c 1 + c 2 ± 4 3 ( A x + B y + c 3 ) 3 2

u ( x , y ) = 1 3 ( A y + B ) ( 2 c 1 x + c 2 ) 3 2 a c 1 3 12 A 2 ( A y + B ) 4 + c 3 y + c 4

u ( x , y ) = 9 a A 2 y + c 1 + 4 A ( x + c 2 y + c 1 ) 3 2 ( x + c 2 ) 3 3 a ( y + c 1 ) 2 + c 3 y + c 4

u ( x , y ) = 3 7 a A 2 ( y + c 1 ) 7 + 4 A ( x + c 2 ) 3 2 ( y + c 1 ) 5 2 ( x + c 2 ) 3 3 a ( y + c 1 ) 2 + c 3 y + c 4

6. Pdsolve Maple Method

The first and second method proposed to study more precisely to find the solution of transonic gas equation. 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 C 1 , C 2 and are global general solution of the second order PDE (transonic gas equation) is

u ( x , y ) = 1 2 c 2 y 2 + C 1 y + C 2 2 3 2 C 3 c 2 x c 2 C 3 2 3 2 C 3 c 2 x c 2 x + C 4 . (35)

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 computation. At first step: Let us, use Dirichlet homogeneous boundary condition and Ω = [ 0 , 1 ] 2 . The system of equation obtained after applying the boundary conditions is:

C 2 2 3 2 C 3 c 2 x c 2 C 3 2 3 2 C 3 c 2 x c 2 x + C 4 = 0

1 2 c 2 + C 1 + C 2 2 3 2 C 3 c 2 x c 2 C 3 2 3 2 C 3 c 2 x c 2 x + C 4 = 0

1 2 c 2 y 2 + C 1 y + C 2 2 3 2 C 3 c 2 C 3 + C 4 = 0

1 2 c 2 y 2 + C 1 y + C 2 2 3 2 C 3 c 2 c 2 C 3 2 3 2 C 3 c 2 c 2 + C 4 = 0

this system is solved using Maple 18, we get the following result:

l Trivial Solution: in case

{ C 1 = 0 C 2 = C 4 C 3 = C 3 C 4 = 0

we remark C 3 we have null solution u = 0 .

l Non-Trivial Solution:

C 1 = 4 9 2 3 i 2 ( γ 1 ) ( 1 / 2 + i / 6 3 ) 12 ( 1 2 + i / 6 3 ) x 2 4 x 3 + 12 ( 1 2 + i / 6 3 ) x + 4 x y 2 ( y 2 2 y + 1 )

where

γ 1 = i 2 + 2 3 2 9 x 3 3 i x 3 3 + 1 + i 3 9 x 2 + 9 i x 2 3 6 i x 3

C 2 = β 1 ( ( β 5 1 / 2 + i / 6 3 ) 48 ( 1 2 + i / 6 3 ) x 2 + 9 C 4 y 2 16 x 3 + 48 ( 1 / 2 + i / 6 3 ) x 9 C 4 y + 16 x )

where

β 1 = 1 9 1 y ( y 1 )

β 2 = 2 9 x 3 3 i x 3 3 + 1 + i 3 9 x 2 + 9 i x 2 3 6 i x 3 ( 1 2 + i / 6 3 ) 12 ( 1 / 2 + i / 6 3 ) x 2

β 3 = ( 1 + i / 3 3 ) ( 2 3 i 2 ( i 2 + 2 3 β 2 ) 4 x 3 + 12 ( 1 2 + i / 6 3 ) x + 4 x )

β 4 = 4 β 3 y 2 ( y 2 2 y + 1 ) 2 ( 1 / 2 + i / 6 3 ) y 2 + 4 2 β 2 y 2 ( y 2 2 y + 1 ) ( 1 2 + i / 6 3 ) y

β 5 = β 4 + 8 3 i 2 ( i 2 + 2 3 2 9 x 3 3 i x 3 3 + 1 + i 3 9 x 2 + 9 i x 2 3 6 i x 3 )

c 2 = 16 3 i 2 ( i 2 + 2 3 2 9 x 3 3 i x 3 3 + 1 + i 3 9 x 2 + 9 i x 2 3 6 i x 3 ) α 1 9 y 2 ( y 2 2 y + 1 )

where

α 1 = ( 1 2 + i / 6 3 ) 96 ( 1 2 + i / 6 3 ) x 2 32 x 3 + 96 ( 1 2 + i / 6 3 ) x + 32 x

C 3 = 1 2 + i 1 6 3 .

Then the nontrivial solution is

u ( x , y ) = 2 2 27 ( δ 1 ( 1 3 2 ( i + 1 3 3 ) δ 2 + ( 4 9 i + i x ( x 1 ) ) 3 + 2 x 3 3 x 2 + x ) ( i 3 + 6 x 3 ) y 2 ( y 1 ) 2 + δ 3 )

where

δ 1 = 2 9 ( 3 / 2 i 3 9 x + 9 2 ) 2 3

δ 2 = ( 3 i x 3 + i + 9 i x 2 6 i x ) 3 + 9 x 3 9 x 2 + 1

δ 3 = 1 3 ( i 3 3 ) 2 3 ( δ 4 + ( 4 9 i + i x ( x 1 ) ) 3 + 2 x 3 3 x 2 + x ) ( i 3 3 ) y 2 ( y 1 ) 2

δ 4 = 1 3 2 ( i + 1 3 3 ) ( 3 i x 3 + i + 9 i x 2 6 i x ) 3 + 9 x 3 9 x 2 + 1

Figure 1 presents the shape of curve of transonic gas equation.

7. Central Finite Difference Method

Let us consider the two-dimensional nonlinear PDE equations:

u x u x x + u y y = 0 ,

and the boundary conditions being:

u ( x , 0 ) = 0 , u ( x , 1 ) = 0 , u ( 0 , y ) = 0 , u ( 1 , y ) = 0.

The finite difference schemes take the form

Figure 1. Pdsolve curve of transonic gas equation.

[ u ( x i + 1 , y i ) u ( x i , y j ) Δ x ] [ u ( x i + 1 , y j ) 2 u ( x i , y j ) + u ( x i + y j ) ( Δ x ) 2 ] + [ u ( x i , y j + 1 ) 2 u ( x i , y j ) + u ( x i , y j 1 ) ] ( Δ y ) 2 = 0

Let Δ x = Δ y = h

( u ( x i + 1 , y j ) u ( x i , y j ) ) ( u ( x i + 1 , y j ) 2 u ( x i , y j ) + u ( x i 1 , y j 1 ) ) + h ( u i , j + 1 2 u i , j + u i , j 1 ) = 0

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

Figure 3 presents the interior points of solution of transonic gas equation.

8. 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:

2 x 2 u ( x , y ) = i = 2 2 β i u ( x + i h , y ) + E r r o r

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 O ( h 6 ) . 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. E r r o r = O ( h 5 ) if and only if:

Figure 2. Central finite difference solution.

Figure 3. Interior points solution.

{ 2 β 2 β 1 + β 1 2 β 2 = 0 1 6 β 1 4 3 β 2 + 4 3 β 2 1 6 β 1 = 0 1 24 β 1 1 24 β 1 2 3 β 2 2 3 β 2 = 0 2 β 2 2 β 2 1 2 β 1 1 2 β 1 = 0 β 2 β 2 β 1 β 1 β 0 = 0

this lead to:

β 2 = 1 12 h 2 , β 1 = 4 3 h 2 , β 0 = 5 2 h 2 , β 1 = 4 3 h 2 , β 2 = 1 12 h 2

2 x 2 u ( t , x ) = 1 12 u ( t , 2 h + x ) h 2 + 4 3 u ( t , h + x ) h 2 5 2 u ( t , x ) h 2 + 4 3 u ( t , h + x ) h 2 1 12 u ( t , 2 h + x ) h 2 + E r r o r

E r r = ( D 2 , 2 , 2 , 2 , 2 , 2 ) ( u ) ( t , x ) h 4 90

2 x 2 u ( t , x ) = 1 12 u ( t , 2 h + x ) h 2 + 4 3 u ( t , h + x ) h 2 5 2 u ( t , x ) h 2 + 4 3 u ( t , h + x ) h 2 1 12 u ( t , 2 h + x ) h 2 + ( 6 x 6 u ( t , x ) ) h 4 90

Note that the last term on the RHS is just the leading order term in the error; terms with higher powers of have been omitted [14] . The power of h in the leading error term is called the order of the stencil.

Figure 4 presents the lattice points of transonic gas equation.

Now we approximate the first derivative

t u ( t , x ) = β 0 u ( t , x ) + β 1 u ( s + t , x ) + β 2 u ( 2 s + t , x ) + E r r o r

the same method as second derivative can be done we obtain

t u ( t , x ) = 3 2 u ( t , x ) s + 2 u ( s + t , x ) s 1 12 u ( 2 s + t , x ) s + 1 3 ( 3 t 3 u ( t , x ) ) s 2 .

We will seek a solution on a rectangular region of the plane

( x , y ) [ L , L ] × [ L , L ]

subject to Dirichlet boundary conditions:

u ( x , y = L ) = g 1 ( x ) , u ( x = + L , y ) = g 2 ( y ) , u ( x , y = + L ) = g 3 ( x ) , u ( x = L , y ) = g 4 (y)

( 2 x 2 u ( x , y ) ) x u ( x , y ) + 2 y 2 u ( x , y ) = 0

( 1 12 u ( x 2 h , y ) h 2 + 4 3 u ( x h , y ) h 2 5 2 u ( x , y ) h 2 + 4 3 u ( x + h , y ) h 2 1 12 u ( x + 2 h , y ) h 2 ) . ( u ( x + h , y ) h u ( x , y ) h ) 1 12 u ( x , y 2 h ) h 2 + 4 3 u ( x , y h ) h 2 5 2 u ( x , y ) h 2 + 4 3 u ( x , y + h ) h 2 1 12 u ( x , y + 2 h ) h 2 = 0

[ α 1 ( u 3 , 2 h u 2 , 2 h ) 1 12 g 1 ( x 2 ) h 2 + 4 3 u 2 , 1 h 2 5 2 u 2 , 2 h 2 + 4 3 u 2 , 3 h 2 1 12 u 2 , 4 h 2 α 2 ( u 4 , 2 h u 3 , 2 h ) 1 12 g 1 ( x 3 ) h 2 + 4 3 u 3 , 1 h 2 5 2 u 3 , 2 h 2 + 4 3 u 3 , 3 h 2 1 12 u 3 , 4 h 2 α 3 ( u 3 , 3 h u 2 , 3 h ) 1 12 u 2 , 1 h 2 + 4 3 u 2 , 2 h 2 5 2 u 2 , 3 h 2 + 4 3 u 2 , 4 h 2 1 12 g 3 ( x 2 ) h 2 α 4 ( u 4 , 3 h u 3 , 3 h ) 1 12 u 3 , 1 h 2 + 4 3 u 3 , 2 h 2 5 2 u 3 , 3 h 2 + 4 3 u 3 , 4 h 2 1 12 g 3 ( x 3 ) h 2 ]

where

Figure 4. Lattice points solution.

α 1 = ( 1 12 g 4 ( y 2 ) h 2 + 4 3 u 1 , 2 h 2 5 2 u 2 , 2 h 2 + 4 3 u 3 , 2 h 2 1 12 u 4 , 2 h 2 )

α 2 = ( 1 12 u 1 , 2 h 2 + 4 3 u 2 , 2 h 2 5 2 u 3 , 2 h 2 + 4 3 u 4 , 2 h 2 1 12 g 2 ( y 2 ) h 2 )

α 3 = ( 1 12 g 4 ( y 3 ) h 2 + 4 3 u 1 , 3 h 2 5 2 u 2 , 3 h 2 4 3 u 3 , 3 h 2 1 12 u 4 , 3 h 2 )

α 4 = ( 1 12 u 1 , 3 h 2 + 4 3 u 2 , 3 h 2 5 2 u 3 , 3 h 2 + 4 3 u 4 , 3 h 2 1 12 g 2 ( y 3 ) h 2 )

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.

9. 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.

Conflicts of Interest

The author declares no conflicts of interest regarding the publication of this paper.

Conflicts of Interest

The authors declare no conflicts of interest.

References

[1] Ames, W.F., Lohner, R.J. and Adams, E. (1981) Group Properties of utt = [f(u)ux]x. International Journal of Non-Linear Mechanics, 16, 439-447.
https://doi.org/10.1016/0020-7462(81)90018-4
[2] Torissi, M. and Valenti, A. (1985) Group Properties and Invariant Solutions for Infinitesimal Transformations of a Non-Linear Wave Equation. International Journal of Non-Linear Mechanics, 20, 135-144.
https://doi.org/10.1016/0020-7462(85)90007-1
[3] Bluman, G. and Kumei, S. (1987) Properties of the Wave Equation. Journal of Mathematical Physics, 28, 307-318.
https://doi.org/10.1063/1.527659
[4] Peters, J.E. and Ames, W.F. (1990) Group Properties of the Non-Linear Dynamic Equations of Elastic String. International Journal of Non-Linear Mechanics, 25, 107-115.
https://doi.org/10.1016/0020-7462(90)90042-8
[5] Fung, R.F., Wang, Y.C. and Wu, J.W. (1999) Group Properties and Group Invariant Solutions for Infinitesimal Transformations of the Non-Linearly Traveling String. International Journal of Non-Linear Mechanics, 34, 693-698.
https://doi.org/10.1016/S0020-7462(98)00047-X
[6] Ozkaya, E. and Pakdemirli, M. (2000) Lie Group Theory and Analytical Solutions for the Axially Accelerating String Problem. Journal of Sound and Vibration, 230, 729-742.
https://doi.org/10.1006/jsvi.1999.2651
[7] Rega, G., Srinil, N., Lacarbonara, W. and Chucheepsakul, S. (2007) Resonant Nonlinear Normal Modes of Inclined Sagged Cables. Journal of Computational and Nonlinear Dynamics, 2, 324-336.
[8] Srinil, N., Rega, G. and Chucheepsakul, S. (2004) Three-Dimensional Non-Linear Coupling and Dynamic Tension in the Large—Amplitude Free Vibrations of Arbitrarily Sagged Cables. Journal of Sound and Vibration, 269, 823-852.
https://doi.org/10.1016/S0022-460X(03)00137-8
[9] Srinil, N., Chucheepsakul, S. and Rega, G. (2005) Internally Resonant Nonlinear Free Vibrations of Horizontal/Inclined Sagged Cables. ASME 2005 International Design Engineering Technical Conferences and Computers and Information in Engineering Conference, Long Beach, 24-28 September 2005, 2021-2029.
[10] Altoum, S.H. (2016) Lie Group and RK4 for Solving Nonlinear First Order ODEs. International Journal of Applied Mathematical Research, 5, 117-122.
https://doi.org/10.14419/ijamr.v5i2.6033
[11] Bertschinger, E. (2002) Introduction to Tenser Calculus for General Relativity.
[12] Ibragimov, N.H. (1996) Elementary Lie Group Analysis and Ordinary Differential Equations. John Wiley Sons, New York.
[13] John, M.L. (2002) Introduction to Smooth Manifolds. Springer Verlag, Berlin.
[14] Chnitia, C., Alhazmic, S.E.A., Altouma, S.H. and Toujanid, M. (2016) DtN and NtD Surface Radiation Conditions for Two-Dimensional Acoustic Scattering: Formal Derivation and Numerical Validation. Applied Numerical Mathematics, 101, 53-70.
https://doi.org/10.1016/j.apnum.2015.08.013

Copyright © 2024 by authors and Scientific Research Publishing Inc.

Creative Commons License

This work and the related PDF file are licensed under a Creative Commons Attribution 4.0 International License.