Next Article in Journal
Non-Hermitian Chiral Magnetic Effect in Equilibrium
Next Article in Special Issue
Method of Third-Order Convergence with Approximation of Inverse Operator for Large Scale Systems
Previous Article in Journal
On Two-Dimensional Fractional Chaotic Maps with Symmetries
Previous Article in Special Issue
A Self-Adaptive Extra-Gradient Methods for a Family of Pseudomonotone Equilibrium Programming with Application in Different Classes of Variational Inequality Problems
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Fast Convergence Methods for Hyperbolic Systems of Balance Laws with Riemann Conditions

by
Kamel Al-Khaled
* and
Nid’a M. Rababah
Department of Mathematics and Statistics, Faculty of Science, Jordan University of Science and Technology, Irbid 22110, Jordan
*
Author to whom correspondence should be addressed.
Symmetry 2020, 12(5), 757; https://doi.org/10.3390/sym12050757
Submission received: 10 March 2020 / Revised: 30 March 2020 / Accepted: 9 April 2020 / Published: 6 May 2020
(This article belongs to the Special Issue Iterative Numerical Functional Analysis with Applications)

Abstract

:
In this paper, we develop an accurate technique via the use of the Adomian decomposition method (ADM) to solve analytically a 2 × 2 systems of partial differential equation that represent balance laws of hyperbolic-elliptic type. We prove that the sequence of iteration obtained by ADM converges strongly to the exact solution by establishing a construction of fixed points. For comparison purposes, we also use the Sinc function methodology to establish a new procedure to solve numerically the same system. It is shown that approximation by Sinc function converges to the exact solution exponentially, also handles changes in type. A numerical example is presented to demonstrate the theoretical results. It is noted that the two methods show the symmetry in the approximate solution. The results obtained by both methods reveal that they are reliable and convenient for solving balance laws where the initial conditions are of the Riemann type.

1. Introduction

Many mathematical and practical problems in physics can be modeled by balance laws, which are systems of nonlinear partial differential equations, usually of hyperbolic type. Balance laws describe very large branches of scientific fields, specifically in fluid and gas dynamics, quantum mechanics, and astrophysics.
In this paper, we present the mathematical formulation for finding numerical solution using Sinc basis functions for 2 × 2 systems of balance laws, while for comparison purposes, we also use the Adomian series technique to find approximate analytical solutions for the same model system that which can be written as
t U + x F ( U ) = G ( U )
The solution U, which is a vector of two components u ( x , t ) that might represent the length of a wave, and another component, say v ( x , t ) which is the velocity of that wave. i.e., U ( x , t ) = [ u ( x , t ) , v ( x , t ) ] T . On the other hand the flux function can be represented as F ( U ) = [ f ( u , v ) , g ( u , v ) ] T , and finally, the source term has the form G ( U ) = [ h 1 ( x , t ) , h 2 ( x , t ) ] T . where F is a smooth function on I R 2 of U, and the derivative of F possesses strictly nonlinear or linearly degenerate eigenfunctions and eigenvalues. The system (1) is of mixed type, it is of elliptic type, for all values of u , v that satisfy ε = { ( u , v ) R 2 : ( f u + g v ) 2 < ( f u + g v ) } and of hyperbolic type for data that lying in H = { ( u , v ) R 2 : ( f u + g v ) 2 > ( f u + g v ) } . We consider in Equation (1) the simplest piece-wise constant initial data, i.e., for certain dispositions of the Riemann initial data
u | t = 0 = u , x 0 u + , x > 0 , v | t = 0 = v , x 0 v + , x > 0
We solve Equation (1) subject to the initial conditions
u ( x , 0 ) = u 0 ( x ) , v ( x , 0 ) = v 0 ( x )
that turn (1) and (3) into a system that is of mixed type (hyperbolic-elliptic). In recent years, the numerical (approximate) solutions for (1) was considered by many authors. In [1] they investigates the Haar-Sinc collocation method for the solution of hyperbolic systems of partial differential equations with applications to the telegraph equation. Existence and uniqueness for a 2 × 2 system of conservation laws, called Brio system, that arising as a model used in magnetohydrodynamics is developed in [2]. In [3] two mixed systems of conservation laws are considered, where numerical and analytical solutions for these two of ill-posed systems are presented. The Adomian decomposition method was used in [4] to develop a fast and accurate algorithm for systems of conservation laws of mixed hyperbolic elliptic type. A survey of several numerical methods for solving nonlinear systems of conservation laws are presented in [5]. A numerical method based on the use of Sinc-collocation method for solving parabolic problems was proposed in [6]. For the construction of numerical solutions; see [7,8,9]. In [10], the Sumudu decomposition method was developed to solve general form of fractional partial differential equations. A new convergence proof for the solution conservation laws of mixed type using the ADM was proposed in [11]. This is with the knowledge that there are many previous studies dealing with this topic, and some of them gave examples from physics and other sciences, to demonstrate the importance of system appeared in (1), and how it is used in practical life, for example we refer to [3,12,13,14].
In this paper, we change the given partial differential equation into a Volterra integral equation, where we discretized in space direction by means of Sinc methodology. After that we apply the indefinite integral formula in time direction which is based on Sinc basis functions. As a result, we are replacing integrals and derivatives with their discrete Sinc approximations, to produce a discrete system of matrices. There are several reasons to approximate by Sinc functions, we mentioned, that they are easily performed and give good accuracy for mixed type problems, also approximation by Sinc are typified by errors of the form O ( exp ( c / h ) ) where c > 0 is a constant and h is a step size.
This research is written in the following order. Initially, in Section 2, we review the important theories and definitions that were mentioned in Sinc theory that would be helpful in writing and building the approximate solution to the issue under study, while Section 3 is devoted to employing what was mentioned in the second section to write the solution using the Sinc method and converting differential equations into a system of algebraic equations, for which solutions can be found in traditional methods such as iterative techniques. As for the Section 4, it is devoted to finding the approximate solution by using the Adomian series, in adding to that, we state and proof a theorem that related to the convergence of the solution. In addition, some relationships related to the order of convergence of the approximate solution are also mentioned. In Section 5, the two methods under discussion are discussed in one example, and all the observations resulting from the solution of the example were mentioned as benefits of the two methods and summarized in tables and graphs. Finally, we will introduce a new issue called the Telegraph system, which is still new and has not been solved by Sinc methodology, and we hope to refer to it in upcoming research.

2. Sinc-Approximation Formula

The Sinc Galerkin method was developed by Stenger [15]. For the equation given in (1), together with the initial condition in (3), the Sinc–Galerkin scheme can be developed in both space and time directions as follows: Since our space domain is I R and time is ( 0 , T ) , we define conformal maps ϕ ( z ) = z which mapping from D E , the eye-shaped domain in the z plane onto the infinite strip, D S , where
D E = { z = x + i y : | arg z a b z | < d π / 2 } .
For the sinc method, the basis functions on the interval ( , ) for z D E , and for the time interval ( 0 , T ) are derived from the composite translated sinc functions
S ( m , h x ) ϕ ( x ) = sinc ϕ ( x ) m h x h x , m = N x , , N x
S ( m , h t ) Υ ( t ) = sinc Υ ( t ) m h t h t , k = N t , , N t
where
sinc ( z ) = sin ( π z ) π z , z 0 1 , z = 0
and,
sinc ( k , h ) ( z ) = sinc z k h h = sin ( π ( ( z k h ) / h ) ) π ( ( z k h ) / h ) , z k h 1 , z = k h , k = 0 , 1 , 2 ,
Let f be a function defined on I R , then for h > 0 we define the Whittaker cardinal series
C ( f ; h ) = k = f ( k h ) s i n c x k h h ,
The properties of (6) were studied and surveyed in depth in Stenger and in Lund [15,16,17]. They are based on the infinite strip D d defined in C as
D d = { w = t + i s : | s | < d π 2 , t R }
To approximate by the Sinc method, we will refer to an important class of functions, see [17] called Play-Wiener class, denoted by W ( π / h ) which is the family of all analytic functions that are in C , and satisfy some decaying conditions. For a function f and its derivative f to be in the class W ( π / h ) , the k t h derivative of the function f can be approximated by
f ( k ) ( x ) = h k n = i = f ( i h ) δ i n ( k ) S ( n , h ) ( x )
where δ i n ( k ) denotes the k t h derivative of the cardinal function. Now we present a fundamental definition that has an important role in the development of Sinc methods on arc Γ in C . For the simply connected domain D, let ϕ be a map that take D onto D d with ϕ ( a ) = and ϕ ( b ) = . If the inverse map of ϕ is denoted by ψ , define
Γ = { z C : z = ψ ( u ) , u R } .
When we evaluate a function and its derivative with respect to ϕ , we obtain
δ i n ( 0 ) = S ( n , h ) ϕ ( x ) x = x i = 1 , i = n 0 , i n
δ i n ( 1 ) = h d d ϕ S ( n , h ) ϕ ( x ) x = x i = 0 , i = n , ( 1 ) n i ( n i ) , i n ,
Since we will replace the first derivative by it matrix approximations, here we define the m × m matrix I ( 1 ) whose i n t h entry is given by δ i n ( 1 ) . An important class of functions denoted by L α ( D ) play crucial rule in approximating derivatives will be needed in next formula (see, [17]). Let F L α ( D ) , α is a positive constant, then for h = π d ( α N ) we have
sup x Γ | F ( n ) ( x ) ( d d x ) n k = N N F ( z k ) q ( z k ) q ( x ) S ( k , h ) ϕ ( x ) | C 3 N n + 1 2 e x p ( π d α N )
An application to the above formula, and in order to approximate the first derivative with respect to x for the function u ( x , t ) in the domain ( , ) leads to the approximation
u x ( t ) 1 h I m x ( 1 ) u ( t ) A u ( t )
where A = 1 h I m x ( 1 ) and u ( t ) = u ( x N x , t ) , , u ( x N x , t ) T . With m x = 2 N + 1 and the skew-symmetric matrix I m x ( 1 ) is defined by (9). Finally, we shall indicate a general formula for approximating the integral a ν F ( t ) d t , ν Γ . If F W ( π / h ) , then for all k Z ,
a k h F ( t ) d t = h n = δ k n ( 1 ) F ( n h ) , δ k ( 1 ) = 1 2 + o k sin ( π t ) π t d t .
To simplify the situation, we define the new matrix I ( 1 ) , as, for N positive integer, define the ( 2 N + 1 ) × ( 2 N + 1 ) matrix I ( 1 ) by
I ( 1 ) = δ k n ( 1 ) , k , n = N , , N .

3. The Sinc–Galerkin Method: Balance Laws

Among the benefits of the Sinc method, is that the approximate solution is approaching the exact solution exponentially, as well as the Sinc method, it carefully handles systems in which the state changes from elliptic to hyperbolic and vis versa, as this depends on the change of the eigen-values of the system. To find an approximate solution for the balance law appeared in (1), without loss of generality, we begin with a 1-dimensional balance law defined in the domain R × ( 0 , T 0 )
u t ( x , t ) + F ( u ( x , t ) ) x = H ( x , t ) ,
where the initial condition is taken to be
u ( x , 0 ) = k ( x ) , x R
Suppose k ( x ) L α ( D ) , and for the step-size h = π d α N .
Integrate Equation (12) with respect to time, and collocation with respect to the variable x will produce a system of Volterra integral equations of the form
u ( t ) = 0 t [ H F ( u ( τ ) ) A u ( τ ) ] d τ + k
where the column vector u ( t ) represent u ( t ) = u N x ( t ) , , u N x ( t ) T , here we use the notation u i ( t ) = u i ( x i , t ) , k = k ( z N x ) , , k ( z N x ) T and H = H ( z N x ) , , H ( z N x ) T , while the square matrix A = 1 h x I m x ( 1 ) represents the discrete Sinc approximation for the first derivative. Subsequently an application of the use of the indefinite integration formula with respect to time variable t, and since our t ( 0 , T 0 , then the domain D E = { t = x + i y : | a r g ( t / ( T 0 t ) ) | < d π 2 } is conformally mapped onto D d via χ ( t ) = ln t / ( T 0 t ) . Define S ( k , h t ) χ ( t ) , k = N t , , N t , to be the basis function for the interval ( 0 , T 0 ) . Now, we define a matrix E = h t I m t ( 1 ) D 1 χ ( t j ) with the nodes t j = χ 1 ( j h t ) , j = N t , , N t , where
D 1 χ ( t j ) = d i a g D 1 χ ( t N t ) , , D 1 χ ( t N t )
h t = π d α N t , and I m t ( 1 ) as defined in (11), with m t = 2 N t + 1 . Then the solution of Equation (12) is in matrix form is given by the rectangular m x × m t matrix U = [ u i j ] :
U = H F ( U ) A U E T + K
we use ∘ to represents the Hadamard multiplication. For convergence purposes for the obtained solution we state two Theorems, in which the proof can be done by just resemble the theorems in [16].
Theorem 1
([16]). Let the matrix U be as defined in (15). Then for the number of points in both space and time N x , N t to be bigger than some constant that depends on d , α , we have
sup ( x i , t j ) u ( x i , t j ) U C N exp ( π d α N )
where N = min { N x , N t }
To guarantee that our approximate solution as obtained in the discrete system (15) converges uniquely to the true solution, we state the following theorem.
Theorem 2
([16,18]). For any positive constant R, we can find another positive constant T 0 so that whenever U 1 U 0 < R 2 , then the solution of (15) is unique. Also, the iteration system
U n + 1 = H F ( U n ) A U n E T + K
leads to a solution that converges to the unique solution.
Next, we state and prove a fact regarding the behavior of the solution, whenever the initial condition satisfies a decay condition. Consider the initial value problem for a balance law of the form
t u ( x , t ) + x f ( u ( x , t ) ) = h ( x , t , u ) , ( x , t ) I R × ( 0 , T )
with initial condition u ( x , 0 ) = u 0 ( x ) , x I R . We need to assume some decaying conditions on the right hand side of Equation (16).
Definition 1.
We say a function u ( x , t ) decays exponentially in x and uniformly for t ( 0 , T ) if | u ( x , t ) | C exp ( α | x | ) for some positive constants a l p h a , C independent of t.
Theorem 3.
Suppose u 0 ( x ) L α ( D ) , and can be written as u 0 ( x ) = sech ( α x ) g ( x ) , where the function g ( x ) is analytic in the strip D d and grows no more than a polynomial in the strip D d , and suppose the source term h ( x , t , u ) and its derivatives with respect to x , y , u are all exponentially decaying in x and uniformly for small t. Then the solution u ( x , t ) of (16) is in the class L α ( D ) .
Proof. 
In order to find the solution of (16), we consider the system of first order ODEs.
d t 1 = d x f ( u ) = d u h
or,
d x d t = f ( u ( x , t , u ) ) , d u d t = h ( x , t , u )
such that
x ( 0 ) = ξ , u ( 0 ) = u 0 ( ξ )
The system (17) with condition (18) implies that
x u = ξ u 0 ( ξ ) + 0 t f ( u ) h d τ .
However, h ( x , t , u ) and its derivatives with respect to the variables x , and t are all exponentially decaying in x and uniformly for small t. This implies that
f ( u ) h
is Lipschitz in the two variables u , x . In addition, so, by the Picard theorem, the solution u ( x , t ) is unique. □

Treatment of Non-Zero Boundary Conditions

The obtained solution in (15) was valid if the initial condition satisfy some certain decaying conditions, i.e, the initial condition has to belong to the class L α . In Section 5, we will discuss problems which have non-zero boundary conditions with the Riemann-type initial condition, which means that the initial condition is no longer in L α . To illustrate the situation, we consider the problem in (1) with the initial condition (3), and since the Sinc functions composed with the various conformal maps, S ( m , h x ) ϕ , are zero at the end-points of the domain of the problem, and since the boundary conditions are non-homogeneous Dirichlet conditions, then the change of variables
u ˜ ( x , t ) = u ( x , t ) u exp ( x ) + u + exp ( x ) exp ( x ) + exp ( x )
(similarly for v ( x , t ) ), will convert system (1) to another one with homogeneous Dirichlet conditions, with smooth initial condition. We substitute the transformation (19) into Equation (1) that leads to a new balance law for the new unknown u ˜ ( x , t ) with an initial condition that is in the class L α for | x | large. Therefore, to find the solution for u ˜ ( x , t ) we can proceed as mentioned above.

4. The Adomian Decomposition Method (ADM)

A lot of published research papers talked about the method of Adomian, and presented it in a detailed way. Here, for the sake of non-repetition, we will suffice to mention the references that contributed to the development of the method. In addition, because of its abundance and where we cannot limit it, we mention the following [19,20,21,22]. We start by rewriting the system in (1) as follows:
u t + f u u x + f v v x = h 1 ( x , t ) , u ( x , 0 ) = u 0 ( x ) v t + g u u x + g v v x = h 2 ( x , t ) , v ( x , 0 ) = v 0 ( x )
or, we may express Equations (20) using linear differential operators as
L t u + f u L x u + f v L x v = h 1 ( x , t ) , u ( x , 0 ) = u 0 ( x ) L t v + g u L x u + g v L x v = h 2 ( x , t ) , v ( x , 0 ) = v 0 ( x )
In the above, the operators L t , L x are defined to be L t = t and L x = x . Performing the inverse operator L t 1 = 0 t d t PDE system in (21) yields
u ( x , t ) = u 0 ( x ) L t 1 h 1 ( x , t ) ϕ 1 ( u , v ) + ϕ 2 ( u , v ) v ( x , t ) = v 0 ( x ) L t 1 h 2 ( x , t ) ϕ 3 ( u , v ) + ϕ 4 ( u , v )
where ϕ 1 ( u , v ) = f u u x , ϕ 2 ( u , v ) = f v v x , ϕ 3 ( u , v ) = g u u x and ϕ 4 ( u , v ) = g v v x . The ADM [19,20] suggests that our approximate solution for unknown functions u ( x , t ) and v ( x , t ) are given by the infinite series
u ( x , t ) = n = 0 u n ( x , t ) , v ( x , t ) = n = 0 v n ( x , t )
while the four nonlinear operators ϕ i ( u , v ) , i = 1 , 2 , 3 , 4 can be expressed as an infinite series of Adomian polynomials as ϕ ( u , v ) = n = 0 A n i , where A n i are what is called Adomian’s polynomials that can be computed according to specific formula set in [20]. Here, the general formulas for Adomian’s polynomials A n i , i = 1 , 2 , 3 , 4 , has the form
A n i ( u 0 , , u n , v 0 , , v n ) = 1 n ! d n d λ n ϕ i k = 0 λ k u k , k = 0 λ k v k λ = 0 , n 0 .
These formulas are smoothly found using coding programs to get as many terms as we want. Plug the nonlinear system (23) into Equation (22) and we arrive at the recursive relation
u 0 ( x , t ) = u 0 ( x ) , u n + 1 ( x , t ) = L t 1 h 1 ( x , t ) A n 1 + A n 2 , n 1 v 0 ( x , t ) = v 0 ( x ) , v n + 1 ( x , t ) = L t 1 h 2 ( x , t ) A n 3 + A n 4 , n 1
The construction of the solutions of u ( x , t ) and v ( x , t ) are given by
u ( x , t ) = lim n k = 0 n u k ( x , t ) , v ( x , t ) = lim n k = 0 n v k ( x , t )
Regarding the convergence of the above two series, there are previous studies that dealt with the topic in detail and in general, we mention [21].

Convergence of the ADM Approximation

Here we give a proof to show that our obtained solution using ADM converges to the exacts solution via the use of the fixed-point theory. The partial differential equation that we are encountered to solve is
u t ( x , t ) + f ( u ( x , t ) ) x = h ( x , t ) , u ( x , 0 ) = g ( x ) , x I R .
which can be written as
u N ( u ) = F ( x , t ) , u H
where H is the Hilbert space, and N = L t 1 [ f ( u ) x ] is the nonlinear operator N : H H and the source term is F ( x , t ) = L t 1 [ h ( x , t ) ] + g ( x ) that is also in H. For the solution using ADM, we assume the solution u and the nonlinear functions are decomposed into an infinite series of the form
u = n = 0 u n , N ( u ) = n = 0 A n .
Substituting these series back into Equation (26) yields
n = 0 u n n = 0 A n = F ( x , t ) .
Comparing terms, we end up with the algorithm u 0 = F ( x , t ) , u n + 1 = A n ( u 0 , u 1 , , u n ) . Then the approximate solution is given by a finite sum of the series S n = n = 1 n u n = u 1 + u 2 + + u n . If S 0 = 0 , S n + 1 = N ( S n + u 0 ) = k = 0 n A k . If the limit exist, then S = lim n S n is the solution of the fixed-point functional equation S = N ( u 0 + S ) H . Need to show that the sequence { S n } is Cauchy. The distance between two consecutive iterations S n and S n + 1 is a ball S r of radius r is given by
S n + 1 S n = u 0 + u 1 + + u n + 1 u 0 u 1 u n = u n + 1 r u n r 2 u n 1 r n + 1 u 0
To prove that { S n } is Cauchy sequence, then for all positive integers m , n where m > n , we have
S m S n = ( S m S m 1 ) + ( S m 1 S m 2 ) + + ( S n + 1 S n ) r m u 0 + r m 1 u 0 + + r n + 1 u 0 = ( r m + r m 1 + + r n + 1 ) u 0 ( r n + 1 + r n + 2 + ) u 0 = r n + 1 ( 1 + r + r 2 + ) u 0 = 1 1 r u 0
Choose r ( 0 , 1 ) such that, S m S n can be made arbitrarily small provided n is large enough for all m , n I N . With this choice of r and u 0 , u 1 , we see that all iterations will remain in the ball
S r = { u ˜ : u ˜ < 1 1 r u 0 } .
There is an integer N such that u m u n < ϵ , n > N , for any m. So the sequence is Cauchy in the Hilbert space. Hence converges to some S H . , i.e., lim n S n = S , where S = n = 0 u n . Solving Equation (26) is the same as solving the functional equation N ( S + u 0 ) = S , by assuming that N is a continuous operator, we get
N ( S + u 0 ) = N ( lim n ( S + u 0 ) ) = lim n N ( S n + u 0 ) = lim n S n + 1 = S
So S is the solution of Equation (26).
When viewing numerical calculations in the last section, and for purposes of comparison we will calculate what is called the order of convergence, we will provide the following definition
Definition 2.
A sequence S n is said to converge to S of order p I N , if there is a positive constant α such that
lim n | S n + 1 S | | S n S | p = α .
Let us examine the order of convergence for the Cauchy sequence S n . We write Taylor’s series for N ( S n + u 0 ) about ( S + u 0 ) as
N ( S n + u 0 ) = N ( N + u 0 ) + N ( S + u 0 ) ( S n S ) + 1 2 ! N ( S + u 0 ) ( S n S ) 2 + + 1 m ! N ( m ) ( S + u 0 ) ( S n S ) m +
However, it is known that N ( S + u 0 ) = S and N ( S n + u 0 ) = S n + 1 , we may rewrite the above equation as
S n S = N ( S + u 0 ) ( S n S ) + + 1 m ! N ( m ) ( S + u 0 ) ( S n S ) m .
Assume our operator N C p [ a , b ] with the property that N ( m ) ( S + u 0 ) = 0 , for all m = 1 , 2 , , p 1 , while it satisfy N p ( S + u 0 ) 0 , therefore we have
S n + 1 S = 1 p ! N ( p ) ( S + u 0 ) ( S n S ) p + 1 ( p + 1 ) ! N ( p + 1 ) ( S + u 0 ) ( S n S ) p + 1 +
divide by ( S n S ) p and take limit as n , we obtain
lim n | S n + 1 S | | S n S | p = lim n 1 p ! N ( p ) ( S + u 0 ) + 1 ( p + 1 ) ! N ( p + 1 ) ( S + u 0 ) ( S n S ) + .
Using the fact that S n is Cauchy that converges to S, the above equation reduces to
lim n | S n + 1 S | | S n S | p = lim n 1 p ! N ( p ) ( S + u 0 ) = α 0 .
Therefore, the order of convergence of S n is p.
For the purposes of calculating the order of convergence, we assume that the first four iterations of the computed solution are given by x k 2 , x k 1 , x k and x k + 1 so the following formula will be used to calculate the order of convergence
q = log ( x k + 1 x k ) log ( x k x k 1 ) log ( x k x k 1 ) log ( x k 1 x k 2 ) .
As mentioned in [23], the first study comparing ADM and Picard’s method was in 1987 by Rach [24]. The study was done by reviewing several examples. Rach showed that solving differential equations of linear type using ADM was equivalent to the classical method of successive approximations, which is known as Picard iteration. Also it was shown that the Picard method gives more accurate solutions, but it needs more time than ADM; this happens because the main advantage of Adomian’s method over the Picard method is the ease of computation of successive terms.

5. Applications: Riemann Type

To investigate the accuracy of the ADM compared with the Sinc method, we choose examples with known solutions that allows for a more complete error analysis. So a 2 × 2 system of balance laws with Riemann type conditions is tested numerically by using the Sinc function methodology, and for comparison purposes, we also solve the same model by ADM. The example reported here is selected to show the convergence of the two schemes. Consider the system
t U + x F ( U ) = H ( x , t ) , ( x , t ) I R × ( 0 , T )
where T is some small constant. In addition,
U ( x , t ) = u ( x , t ) v ( x , t ) , H ( x , t ) = h 1 ( x , t ) h 2 ( x , t ) .
with flux function F ; I R 2 I R 2 , such that
F ( u , v ) ) = f ( u , v ) g ( u , v ) = 1 2 v 2 1 2 u 2 .
with the approximations to Riemann type condition
U ( x , 0 ) = u ( x , 0 ) v ( x , 0 ) = 1 1 , x < 0 1 1 , x > 0
given by
U ( x , 0 ) = u ( x , 0 ) 2 v ( x , 0 ) = tanh ( 3 x / 2 ) tanh ( 3 x / 2 ) .
The left and right boundary conditions are
U ( , t ) = u ( , t ) v ( , t ) = 1 1 + t 2 exp ( 5 t ) , U ( , t ) = u ( , t ) v ( , t ) = 1 1 + t 2 exp ( 5 t ) .
Here the true solution given by
u ( x , t ) = 1 1 + t 2 tanh ( 3 x / 2 ) , v ( x , t ) = exp ( 5 t ) tanh ( 3 x / 2 ) .
Thus, we can calculate the appropriate source term by applying the differential operator to the solution. The source term H ( x , t ) in (28) turns out to be
h 1 ( x , t ) = sinh ( 3 x / 2 ) ( 4 t cosh 2 ( 3 x / 2 ) + 3 exp ( 10 t ) + 6 t 2 exp ( 10 t ) + 3 t 4 exp ( 10 t ) ) 2 cosh 3 ( 3 x / 2 ) ( 1 + 2 t 2 + t 4 )
and,
h 2 ( x , t ) = 5 sinh ( 3 x / 2 ) ( exp ( 5 t ) cosh 2 ( 3 x / 2 ) + 2 exp ( 5 t ) cosh 2 ( 3 x / 2 ) t 2 + exp ( 5 t ) cosh 2 ( 3 x / 2 ) t 4 3 10 ) cosh 3 ( 3 x / 2 ) ( 1 + 2 t 2 + t 4 ) .
In Section 3 when we construct the approximate solution by means of the Sinc methodology, it was necessary that the initial condition belonged to the class of functions L α ( D ) previously mentioned in Section 3, but we know that the Riemann condition does not belong to this family, so a transformation had to be made to make the primitive condition belong to the family. Since the boundary conditions are non-homogeneous, then the transformation
U ˜ ( x , t ) = u ˜ ( x , t ) v ˜ ( x , t ) = u ( x , t ) v ( x , t ) w 1 ( x , t ) w 2 ( x , t )
where w 1 ( x , t ) = 1 1 + t 2 exp ( x ) + 1 1 + t 2 exp ( x ) exp ( t ) + exp ( z ) , and w 2 ( x , t ) = exp ( 5 t ) exp ( x ) + exp ( 5 t ) exp ( x ) exp ( t ) + exp ( z ) will convert the given boundary conditions to homogeneous conditions, provided that | x | . Now, after substituting the transformation (29) into Equation (28), we get a new system with the unknown U ˜ ( x , t ) = ( u ˜ ( x , t ) , v ˜ ( x , t ) ) T given by
t U ˜ + x F ( U ˜ ) = H ( x , t ) t W
with initial condition
U ˜ ( x , 0 ) = U ( x , 0 ) W ( x , 0 )
where W ( x , 0 ) = ( w 1 ( x , 0 ) , w 2 ( x , 0 ) ) T . Equation (30) can be solved for U ˜ ( x , t ) , and then, using the transformation in (29), for U ( x , t ) . To proceed, since the space and the time domains are ( , ) and ( 0 , T ) , respectively, choose the conformal maps ϕ ( x ) = x and Υ ( t ) = log ( t T t ) . We solve the above system using ADM subject to the initial conditions u ( x , 0 ) , v ( x , 0 ) . Integrate the system in (28) with respect to t and using the initial conditions u ( x , 0 ) , v ( x , 0 ) , we arrive at
u ( x , t ) = u ( x , 0 ) + 0 t [ h 1 ( x , τ ) ( v ( x , τ ) v ( x , τ ) ) x ] d t
and,
v ( x , t ) = v ( x , 0 ) + 0 t [ h 2 ( x , τ ) ( u ( x , τ ) u ( x , τ ) ) x ] d t
In the above two integrals, set u ( x , t ) = n = 0 u n ( x , t ) , v ( x , t ) = n = 0 v n ( x , t ) , and for the nonlinear terms, we find Adomian polynomials as
u u x = n = 0 A n , v v x = n = 0 B n .
The first few terms of Adomian polynomials are given by
A 0 = u 0 u 0 x , A 1 = u 0 x u 1 + u 0 u 1 x , A 2 = u 0 x u 2 + u 1 x u 1 + u 0 u 2 x ,
B 0 = v 0 v 0 x , B 1 = v 0 x v 1 + v 0 v 1 x , B 2 = v 0 x v 2 + v 1 x v 1 + v 0 v 2 x ,
Substitute the above infinite sum for Adomian polynomials, together with the approximate solutions into Equations (32) and (33), exactly as mentioned in Section 5. Balancing terms in Equations we obtain an approximate solution of the form
u A ( x , t ) = u 0 ( x , t ) + u 1 ( x , t ) + u 2 ( x , t ) + u 3 ( x , t ) + = tanh ( 3 x / 2 ) ( 1 t 2 + t 4 t 6 + t 8 ) +
v A ( x , t ) = v 0 ( x , t ) + v 1 ( x , t ) + v 2 ( x , t ) + v 3 ( x , t ) + = tanh ( 3 x / 2 ) ( 1 5 t + 25 t 2 2 125 t 3 6 + 625 t 4 24 ) +
Table 1 shows that the method converges for T = 2 using Sinc methodology, where the second column reports the supremum norm of the error between the exact solution and the Sinc approximate u solution, while the third column reports the error between the exact solution and the Sinc approximate v solution. The error in the approximate u-solution using ADM is reported in Figure 1. Figure 2 reported the v ( x ) -solution using ADM for different values of time t. The error in approximating u , v show that most of the error concentrates at origin, which caused by approximating Riemann type condition by a smooth function. Through our study, and upon consulting the Table 1 and Table 2, we are able to find the values for the order of convergence using Equation (27), it was found that the order of convergence for the ADM is q = 1.0007 , while it is approximately q = 2.843 using the Sinc method. This is in accordance with the theories formulated in Section 4, which show that the Sinc methodology is faster, better, and more comprehensive in finding an approximate solution. Because of the nature of the issues under discussion, it is noted that the two methods show the symmetry of the approximate solution. A quick look at in calculated approximate solutions for both u ( x , t ) and v ( x , t ) shows the symmetry of the solutions where the error value at the point ( x , t ) is the same as at ( x , t ) . The main goal of this research paper is to compare the two methods together, but if the reader wants more ways to solve systems from the same context, we cite the following [22]. It was shown in [25] that traditional methods, like finite difference method is unstable when the system is of mixed-type when Riemann conditions are involved, which means that the two methods used in this research are an improvement over the previous traditional methods. To present a physical and engineering application for the purposes of explaining the importance of the topic, please review the research paper [18] where the p system was solved. Finally, what we intend to refer to in a future research in the hope of discussing it from all sides is the importance of the Telegraph system. As an illustration, we shall consider conductors (such as telephone wires or submarine cables) in which the current I = I ( x , t ) may leak to ground. The resultant decrease in current is governed by I x = G V C V t , where V = V ( x , t ) is voltage, G = G ( x ) is the conductance to ground, and C = C ( x ) is the capacitance with the ground. The change in voltage is governed by V x = R I L I t , where R = R ( x ) is the resistance and L = L ( x ) is the inductance of the cable.

Author Contributions

Formal analysis, investigation, and writing-original draft preparation K.A.-K. and N.M.R. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Conflicts of Interest

The authors declare that there is no conflict of interests regarding the publication of this manuscript. The authors declare that they have no competing interests.

References

  1. Pirkhedri, A.; Javdi, H.H.S.; Navidi, H.R. Numerical algorithm based on Haar-Sinc collocation method for solving the hyperbolic PDEs. Sci. World J. 2014, 2014, 340752. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  2. Kalisch, H.; Mitrovic, D.; Teyekpiti, V. Existence and uniqueness of singular solutions for a conservation law arising in magnetohydrodynamics. Nonlinearity 2018, 31, 5463–5483. [Google Scholar] [CrossRef]
  3. Fitt, A.D. The numerical and analytical systems of conservation laws. Appl. Math. Model. 1989, 13, 618–631. [Google Scholar] [CrossRef]
  4. Az-Zo’bi, E.A. Construction of solutions for mixed hyperbolic elliptic Riemann initial value system of conservation laws. Appl. Math. Model. 2013, 37, 6018–6024. [Google Scholar] [CrossRef]
  5. Sod, G.A. A survey of several finite difference methods for systems of nonlinear hyperbolic conservation laws. J. Comput. Phys. 1978, 27, 1–31. [Google Scholar] [CrossRef] [Green Version]
  6. Rashidinia, J.; Barati, A. Numerical solutions of one-dimensional non-linear parabolic equation using Sinc collocation method. Ain Shams Eng. J. 2015, 6, 381–389. [Google Scholar] [CrossRef] [Green Version]
  7. Al-Khaled, K.; Allan, F. Decomposition method for solving nonlinear integro-differential equations. J. Appl. Math. Comput. 2005, 19, 415. [Google Scholar] [CrossRef]
  8. Al-Khaled, K. Cardinal-type approximations for conservation laws of mixed type. Nonlinear Stud. 2014, 21, 423–433. [Google Scholar]
  9. Alquran, M.; Ali, M.; Al-Khaled, K. Solitary wave solutions to shallow water waves arising in fluid dynamics. Nonlinear Stud. 2012, 19, 555–562. [Google Scholar]
  10. Al-Khaled, K. Numerical solutions of time-fractional partial differential equation using Sumudu decomposition method. Rom. J. Phys. 2015, 60, 99–110. [Google Scholar]
  11. Az-Zo’bi, E.A.; Al-Khaled, K. A new convergence proof of the Adomian decomposition method for a mixed hyperbolic elliptic system of conservation laws. Appl. Math. Comput. 2010, 217, 4248–4256. [Google Scholar] [CrossRef]
  12. Fjordholm, U.S.; Lye, K.; Mishra, S.; Weber, F. Statistical solutions of hyperbolic systems of conservation laws: Numerical approximation. arXiv 2019, arXiv:1906.02536. [Google Scholar] [CrossRef] [Green Version]
  13. Hussain, A.; Bano, S.; Khan, I.; Baleanu, D.; Sooppy Nisar, K. Lie Symmetry Analysis, Explicit Solutions and Conservation Laws of a Spatially Two-Dimensional Burgers-Huxley Equation. Symmetry 2020, 12, 170. [Google Scholar] [CrossRef] [Green Version]
  14. Randall, J.L. Numerical Methods for Conservation Laws, Lectures in Mathematics, 2nd ed.; ETH Zürich: Zürich, Switzerland; Birkhäuser: Basel, Switzerland, 1992. [Google Scholar]
  15. Stenger, F. Numerical Methods Based on Sinc and Analytical Functions; Springer: New York, NY, USA, 1993. [Google Scholar]
  16. Al-Khaled, K. Theory and Computation in Hyperbolic Model Problems. Ph.D Thesis, University of Nebraska, Lincoln, NE, USA, 1996. [Google Scholar]
  17. Lund, J.; Rowers, K.L. Sinc Methods for Quadrature and Differential Equations; SIAM: Philadelphia, PA, USA, 1992. [Google Scholar]
  18. Al-Khaled, K. Kamel Al-Khaled; A Sinc-Galerkin approach to the p-system. Saitama Math. J. 1998, 16, 1–13. [Google Scholar]
  19. Adomian, G. A review of the decomposition method in applied Mathematics. J. Math. Anal. Appl. 1988, 135, 501–544. [Google Scholar] [CrossRef] [Green Version]
  20. Adomian, G. Solving Frontier Problems of Physics: The Decomposition Method; Kluwer: Boston, UK, 1994. [Google Scholar]
  21. Cherrualt, Y.; Adomian, G. Decomposition methods: A new proof of convergence. Math. Comput. Model. 1993, 18, 103–106. [Google Scholar] [CrossRef]
  22. Alquran, M.T.; Al-Khaled, K.M. Numerical Comparison of Methods for Solving Systems of Conservation Laws of Mixed Type. Int. J. Math. Anal. 2011, 5, 35–47. [Google Scholar]
  23. El-Sayed, A.M.A.; Hashem, H.H.G.; Ziada, E.A.A. Picard and Adomian decomposition methods for a quadratic integral equation of fractional order. Comput. Appl. Math. 2014, 33, 95–109. [Google Scholar] [CrossRef]
  24. Rach, R. On the Adomian (decomposition) method and comparisons with Picard’s method. J. Math. Anal. Appl. 1987, 128, 480–483. [Google Scholar] [CrossRef] [Green Version]
  25. Holden, H.; Holden, L.; Risebro, N.H. Some qualitative properties of 2 × 2 systems of conservation laws of mixed type. In Nonlinear Evolution Equations That Change Type; Springer: New York, NY, USA, 1990. [Google Scholar]
Figure 1. The error for the u solution using ADM.
Figure 1. The error for the u solution using ADM.
Symmetry 12 00757 g001
Figure 2. The approximate solution v ( x , t ) using ADM for different values of time t.
Figure 2. The approximate solution v ( x , t ) using ADM for different values of time t.
Symmetry 12 00757 g002
Table 1. The absolute error using Sinc method for different values of N x = N t .
Table 1. The absolute error using Sinc method for different values of N x = N t .
N x = N t u v
4 8.5450 × 10 03 4.1314 × 10 02
8 5.7183 × 10 03 1.6575 × 10 02
16 1.3617 × 10 03 3.5789 × 10 03
32 8.8632 × 10 05 2.4625 × 10 04
64 5.3823 × 10 06 9.6498 × 10 05
Table 2. The iterations that use ADM for u ( x , t ) , v ( x , t ) at ( 3 , 0.25 ) .
Table 2. The iterations that use ADM for u ( x , t ) , v ( x , t ) at ( 3 , 0.25 ) .
Iteration No. u ( x , t ) v ( x , t )
First 0.999753 0.289434
Second 0.937269 0.284561
Third 0.941174 0.281119
Fourth 0.904930 0.273621

Share and Cite

MDPI and ACS Style

Al-Khaled, K.; Rababah, N.M. Fast Convergence Methods for Hyperbolic Systems of Balance Laws with Riemann Conditions. Symmetry 2020, 12, 757. https://doi.org/10.3390/sym12050757

AMA Style

Al-Khaled K, Rababah NM. Fast Convergence Methods for Hyperbolic Systems of Balance Laws with Riemann Conditions. Symmetry. 2020; 12(5):757. https://doi.org/10.3390/sym12050757

Chicago/Turabian Style

Al-Khaled, Kamel, and Nid’a M. Rababah. 2020. "Fast Convergence Methods for Hyperbolic Systems of Balance Laws with Riemann Conditions" Symmetry 12, no. 5: 757. https://doi.org/10.3390/sym12050757

APA Style

Al-Khaled, K., & Rababah, N. M. (2020). Fast Convergence Methods for Hyperbolic Systems of Balance Laws with Riemann Conditions. Symmetry, 12(5), 757. https://doi.org/10.3390/sym12050757

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop