The arbitrary-order fractional hyperbolic nonlinear scalar conservation law

In this paper, we use a new powerful technique of arbitrary-order fractional (AOF) characteristic method (CM) to solve the AOF hyperbolic nonlinear scalar conservation law (HNSCL) of time and space. We present the existence and uniqueness of this class of equations in time and one-dimensional space of fractional arbitrary order. We extend Jumarie’s modification of Riemann–Liouville and Caputo’s definition of the fractional arbitrary order to introduce some formulae (Appl. Math. Lett. 22:378–385, 2009; Appl. Math. Lett. 18:739–748, 2005). Then, we use these formulae to prove the main theorem. In the application section, we use the analytical technique that is presented in the theorem to solve examples that are given.

In recent years, considerable attention to time and spatial fractional differential equations (TSFDEs) has been growing rapidly. These problems are deduced by replacing the standard time derivative with time-fractional derivative, and they can be used to describe some physical and industrial processes such as super diffusion and sub-diffusion phenomena [16][17][18][19][20][21][22][23].
In the last two decades, variable-order fractional calculus (VOFC) has attracted many researchers in different fields of science. Ross and Samko introduced the idea of FVOC in 1993 [24,25]. The mathematicians in pure and applied mathematics as well as researchers in physics, chemistry, biology, and engineering, are pursuing this topic. As it is well known, when the order of the fractional operator is variable, some phenomena in physics can be described better than in the case of constant order; for instance, in the diffusion process in an inhomogeneous or heterogeneous medium, or in the processes where the changes in the environment modify the dynamic of the particle [26][27][28][29]. Researchers have considered fractional derivatives of variable order, with α depending on variable t, which was introduced in [28]. We present one of the three types of Caputo fractional derivatives that are defined in the next section. We are going to extend some formulae from constant-order fractional derivative that were introduced by Jumarie [1] for the arbitrary-order fractional derivative. The order of the derivative is considered as a function α(t) taking values on the open interval (0, 1).
The theory of hyperbolic conservation laws arose almost fifty years ago [30,31]. The unique features of this class of systems of partial differential equations (PDEs) had been identified long before. There are many phenomena in mathematical physics that are arising, and their mathematical models are in the form of hyperbolic conservation law. Conservation law problems were presented in the early book by Tyn Myint-U and Lokenath Debnath [30] in Chap. 13. Afterward, some papers have been published. Guo-chang Wu [31] studied a fractional linear conservation law problem that used a fractional characteristic method. Our aim in this paper is the development of this method to address the arbitrary-order fractional nonlinear hyperbolic conservation law. The term time homogeneous hyperbolic conservation law refers to first-order systems of PDEs in the divergence form (1.1) The state vector U = (u 1 , . . . , u n ) T , with values in R n , is to be determined, and u i is a function of the spatial variables (x 1 , . . . , x m ) and time t. The given functions H, c i , and G i , where i = 1, . . . , m, are the smooth maps from R n to R n . Also α : R → (0, 1) and β : R → (0, 1) where α(τ ) and β(τ ) are continuous. The symbol ∂ t stands for ∂ ∂t , ∂ x stands for ∂ ∂x , and ∂ α(τ ) t is an arbitrary-order fractional derivative that will be defined in Sect. 3. The considered problem (1.1) when α = 1 reduces to the classical conservation law, numerical approximations of which have been intensively studied. However, to the best of our knowledge, the analytical solution of time fractional conservation law has not been addressed yet. This article aims to fill the gap and investigates the analytical solution of (1.1). In the present paper, we adopt the fractional characteristic method (FCM), which is a very powerful technique that converts an FPDE to a system of FODEs, which makes it possible to solve (1.1). The FCM method was introduced by Guo-chang Wu [31], and it has been further developed to address the AOF hyperbolic conservation laws. Due to its efficiency in obtaining the exact solution, it becomes a very attractive method for seeking answers to differential equations. The feature of this technique in comparison to the other analytical solution is that it gives us the ability to check if the obtained solution is exact, by substituting the answer in the FPDE and showing it satisfies the differential equation.
The remainder of this paper is organized as follows. In Sect. 2, the definitions of VOF derivative and integral are given. In Sect. 3, we use the definition from the previous section and Jumarie's paper to present the AOF derivative and integral formulae. In Sect. 4, we prove the existence and uniqueness of AOF NHSCL, and in this prosses, the analytic method to solve the AOF FHCL is presented. In Sect. 5, we implemented the analytical approach that is introduced in Theorem 1 to solve a few examples. Also, as a benchmark in Examples 1, 2, and 5, we show that the obtained solution satisfies the AOF NHCL. This test can be done for other cases too. Finally, a summary is given in the last section.

Variable-order Caputo derivatives for functions of one variable
The generalization of the Caputo derivative from constant order to variable order of fractional differentiation is defined in [28]. Given α(t) ∈ (0, 1), the left and right Caputo fractional derivatives of order α(t) of a function x : indicate the left and right Riemann-Liouville fractional derivatives of variable order α(t).

A definition for Riemann-Liouville and Caputo fractional arbitrary-order derivative
In the above definitions the variable t in α(t) and x(t) is the same, which leads to different types of definitions. However, now we would like to present a definition that has different variables for α and x. In this case, we will have only one definition for the Riemann-Liouville and Caputo variable-order fractional derivatives, so it is appropriate to name it arbitrary-order instead of variable-order fractional derivative.

Definition 3 Riemann-Liouville fractional derivatives of arbitrary order α(t):
Given a function f : [a, b] → R and α : R → (0, 1), where f (x) and α(t) are continuous, then: 1. The left Riemann-Liouville fractional derivative of arbitrary order α(t) is defined by 2. The right Riemann-Liouville fractional derivative of arbitrary order α(t) is defined by

The right Caputo derivative of arbitrary order α(t) is defined by
(3.4)

Some results based on the above definition
Considering Definitions 3 and 4, we can extend all the results from Jumarie's paper 32 from the constant-order fractional derivative to the arbitrary-order fractional derivative by replacing α with α(t). Therefore, we present the following results.
Definition 5 Let f : R → R denote a continuous (but not necessarily differentiable) function, and let h > 0 indicate a constant discretization span. The forward operator FW (h) is defined by the equality (the symbol := means that the left-hand side is defined by the right-hand side) Assume α : R → (0, 1) and continuous, then the fractional difference of arbitrary order α(t) for f (x) is defined by the expression and its fractional arbitrary-order derivative is the limit This definition is close to the standard definition of the derivative, and as a direct result, the α(t)th derivative of a constant, 0 < α(t) < 1, is zero.

Modified Riemann-Liouville derivative
Definition 6 ([32]) Revised Riemann-Liouville definition: Based on Definition 4, its fractional derivative of arbitrary order α(t) < 0 is defined by and we can write Also, for n ∈ N, The difference between (3.8) and (3.9) is that the latter contains the constant f (0), while the first one does not. Equation refers to the modified Riemann-Liouville derivative, the constant of which was introduced by Jumarie.
Caputo's definition can be presented as instead of Definition 2, thus, it assumes explicitly that f (x) is differentiable. With this definition, Laplace's transform L{·} in [32] can be presented in the form of the arbitrary-order fractional derivative as follows:

The basic formula for one-variable functions
A generalized Taylor expansion of constant-order fractional, which applies to nondifferentiable functions (F-Taylor series) in [32], can be generalized to the Taylor expansion of arbitrary-order fractional applicable to nondifferentiable functions as follows.

Proposition 1
Assume that the continuous function f : R → R has a fractional derivative of order kα(t) for any positive integer k and α : R → (0, 1) and continuous, then the following equality holds: is the derivative of variable order α(t)koff (x), and with the notation where Γ (·) denotes the Euler gamma function.
Let E α(t) (u) represent the Mittag-Leffler function defined by the expression and D x be the derivative operator with respect to x, then the above series can be written This arbitrary-order fractional Taylor's series does not hold with the standard Riemann-Liouville derivative, and it only applies to nondifferentiable functions.
Then, integrating this series with respect to h provides The arbitrary order of the derivation in f [kβ(t)+m] (x) is very crucial and should be understood as D kβ(t) f (m) (x), since we start with the fractional Taylor's series of f (m) (x).

Fractional Taylor's series for two variable functions
The fractional Taylor's series can be generalized in a straightforward way to two-variable functions to yield Therefore, its total deferential for f (x, y) would be

Basic formulae for fractional derivative and integral
The formulae that we present will be used later to prove the theorem and to solve the problems that are given at the end.

A fractional derivative of compounded functions
Equation (3.14) provides a useful differential relation

Corollary 2
The following equalities hold: , then the following properties are satisfied: Remark 1 We denote that u(x) is nondifferentiable in (3.26) and (3.27) and differentiable in (3.28). v(x) is nondifferentiable, and f (u) is differentiable in (3.27) and nondifferentiable in (3.28).
The concept of differentiability and nondifferentiability is described by Jumarie for fractional constant order as follows [33]: "Loosely speaking, there are two main trends in defining the framework of fractional calculus: on the one hand, the formal definition of the fractional derivative as an antiintegral; and on the other hand, the approach via fractional difference which works exactly like the standard (Leibniz) classical differential calculus. "

Integration with respect to (dx) α(t)
The integral with respect to (dx) α(t) is defined as the solution of the fractional differential equation which is provided by the following result.

Lemma 1 ([2]) Let f (x) denote a continuous function, then the solution of equation (3.31)
is defined by the equality We can also derive the following integrating formulae using (3.29) and (3.30):

Existence and uniqueness
We consider AOF HNSCL in time and one-dimensional space, which is defined as follows: Here, we also use the following modified Riemann-Liouville derivative of the arbitraryorder fractional with parameters α(τ ) and β(τ ) that may not be the same: The solutions of the hyperbolic conservation law for α = β = 1 may be visualized as propagating waves. In the case of a nonlinear system, the profiles of compression waves progressively are getting steeper and finally break, generating jump discontinuities or shocks. Therefore, the theory applies only to weak solutions. This problem is complicated further by the fact that, due to weak solutions, uniqueness is lost. On the other hand, the classical characteristic method is an efficient and powerful technique that has been used to solve the initial value problem of FPDEs analytically by converting to a system of FODEs. It is well known that the method of characteristic has played an important role in mathematical physics. To the best of our knowledge, there has been no published work on proving the existence and uniqueness of equation (4.1); in this paper, we attempt to address this matter.
Theorem 1 Let us consider (4.1) where H(u(x, t) with the initial condition Then this nonlinear initial value problem has a unique solution for all 0 < α(τ ), β(τ ) ≤ 1 such that the following conditions are satisfied: (f (ξ )). F and f are differentiable with respect to x and t; (iii) ξ is fractionally differentiable with respect to x and t; Proof The first part of the proof is to show existence. As the main tool, we use the method of arbitrary-order fractional characteristics.
The total differential of u(x, t) for (4.3) is given by These curves Υ are called the characteristic curves of equation (4.3). Thus the solution of (4.3) is reduced to the solution of a pair of simultaneous FODEs (4.6) and (4.7). According to (4.6), which implies u is constant along each characteristic curve, each c(u) remains constant on Υ , and the solutions of (4.7) form the family of the characteristic curves for (4.3) in the (x, t)-plane. It means that if the family of the curves Υ can be obtained, then the general solution of (4.3) is obtained too. If we assume the initial point on the characteristic curve Υ is given by ξ , it means the curve Υ intersects t = 0 at x = ξ , then u(ξ , 0) = f (ξ ) on the entire curve Υ as shown in Fig. 1.
Thus, the family of the characteristic curves Υ is the solution of the following FODEs: From equation (4.11) and its integration, we have which represents the characteristic curves (they are straight lines when α = β = 1). Therefore equations (4.10) and (4.13) present the solution of the initial-value problem (4.3) in a parametric form where F(ξ ) = c f (ξ ) . Now we show that solution (4.14) satisfies (4.3). We obtain the derivative of the 1st equation in (4.14) regarding x with an order of β(τ ) and t with an order of α(τ ): Then we find the derivative of the 2nd equation in (4.14) regarding x and t respectively with an order of β(τ ) and α(τ ): Solution (4.14) also satisfies the initial condition at t = 0 since ξ = x.
The second part of the proof is the uniqueness. Assume that u(x, t) and v(x, t) are two solutions of equation (4.3); therefore, they should satisfy solution (4.14), that is: Or we can say that on Solution: According to the analytical method introduced in Theorem 1, we can obtain the following FODEs: Γ (1+α(τ )) . Using integrating formula (3.33), we have By the initial condition, the constants C 1 and C 2 are obtained. Hence the parametric solution is which leads to the solution Concerning the initial condition (5.2), we have Benchmark 1 (Example 1) We show that (5.6) satisfies (5.1). We obtain the derivative of the 1st equation in (5.6) regarding x with an order of β(τ ) and t with an order of α(τ ): (5.10) where f (ξ ) = f ξ (ξ ). Then the derivative of the 2nd equation in (5.6) regarding x and t respectively with an order of β(τ ) and α(τ ) is (5.11) By formula (3.24), we obtain (5.13) Substituting (5.13) in (5.10) and (5.14) in (5.1) Solution (5.6) also satisfies the initial condition at t = 0 since ξ = x, and solution (5.6) is unique.
Suppose that u(x, t) and v(x, t) are two solutions. Then, on therefore solution (5.6) is unique.
Its solution is given by The graph of u is given in Fig. 3.
Example 2 Let us consider space-time arbitrary-order fractional inviscid Burger equation, which can be defined as follows: with the initial condition u(x, 0) = f (x).
Remark 4 The classic form of Example 2, where α(τ ) = β(τ ) = 1, is the inviscid Burger equation, which is The parametric solution based on Theorem 1 is wave.
Let f (x) = x 2 + 3x be the initial condition for (5.31), then and the graphs of u are given in Fig. 4.
We are investigating a more complicated form of Example 1, the space-time of arbitrary-order fractional, which is u(x, 0) = f (x).
(5.41) Therefore, u(x, t) = f (ln(e -te -x )) and the graph of the solution u with the initial condition f (x) = x 2 + 3x for different values of α(τ ) and β(τ ) is given in Fig. 5.
Example 4 Let us consider the arbitrary-order fractional differential equation of the form

Solution:
We obtain the following FODEs based on the method that was introduced in Theorem 1: Integrating (5.45) by using formula (3.35), the parametric solution of (5.42) is Considering the initial condition f (x) = x 2 + 3x, the solution u can be written as follows: (5.48) Remark 6 The classic case of Example 4, where α(τ ) = β(τ ) = 1, is and its solution is its graph is given in Fig. 6. The graph of u with another value for α(τ ) and β(τ ) is given in Fig. 7.
Example 5 Consider a bit more complicated case of the arbitrary-order fractional differential equation: And from (4.4) we have , ) .

Summary
We have proved the existence and uniqueness of the arbitrary-order fractional hyperbolic nonlinear scalar conservation law in time and space under certain conditions. We have used the generalization of the classical characteristic method and the generalization of some formulae from the constant-order fractional to the arbitrary-order fractional. And finally, we presented a few physical examples in which we have implemented the analytical technique that was introduced in Theorem 1 to find the solutions, and then we showed the graphs for different values of α(τ ) and β(τ ). In the last example that is more complicated, and 1st and 2nd examples as a benchmark, we showed that the solution satisfies the differential equation. This benchmark can be performed for other cases too.
quite a bit of editing and revising too. Mr. MAOA: He had spent quite a bit of time discussing different subjects and problems related to this paper, which caused the paper's improvement; he also read, revised, and corrected the content of the article. All authors read and approved the final manuscript.