^{1}

^{2}

In this paper we investigate anomalous diffusion coupled with linear convection, using fractional calculus to describe the anomalous associated memory effects in diffusive term. We get an explicit travelling wave solution, wavefront, with finite propagation. We comment the properties of the solution, including the stationary case.

The fractional calculus is almost as old as the calculus of integer order, but only in 1974 was realized the first Conference on Fractional Calculus its Applications. In the decade 90s fractional calculus was recognized and specific journals and textbooks were being published [

We study the diffusion with linear convection, that is, a nonlinear convection-diffusion problem. Appearing in several physical, biological and chemical applications, B. H. Gilding and R. Kersner [

Our purpose is to find an explicit solution, so that we can investigate its properties in the fractional model. In Section 2 we define the Riesz fractional derivative. In Section 3 the fractional model used to describe the convection-diffusion equation using Riesz fractional derivative. In Section 4 we use the similarity reduction method in fractional equation. In Section 5 we calculate the particular travelling wave, called a wavefront. Section 6 makes our concluding remarks.

The Riesz fractional derivative of order α , with 0 < α < 2 and α ≠ 1 is defined by:

D x α f ( x ) = − D + α f ( x ) + D − α f ( x ) 2 cos ( α π / 2 ) , (1)

where D ± α f ( x ) are Weyl fractional derivatives [

Theorem 1. Let be h ( x ) = | x | − α − 1 , with 1 < α < 2 . We describe the Riesz fractional derivative of order α for an appropriate Fourier convolution product is given by:

D x α f ( x ) = d α ( f ∗ h ) ( x ) , (2)

with d α = − 1 2 Γ ( − α ) cos ( α π 2 ) .

This result is an improvement of the theorem developed by E. C. Grigoletto and E. C. de Oliveira [

Analyze a fractional space nonlinear convection-diffusion equation, in which we apply the Riesz fractional derivative in diffusive term. If we consider only diffusion, Equation (3) defines a fractional porous medium equation [

∂ ∂ t u ( x , t ) = D x α u n ( x , t ) + ∂ ∂ x u ( x , t ) , n > 1 , t > 0 and x ∈ ℝ , (3)

where D x α is Riesz fractional derivative, with 1 < α < 2 . In the case α → 2 , we get the integer classical convection-diffusion equation. The parameter n appears, when we consider the diffusion coefficient n u n − 1 ( x , t ) , that is, density dependent diffusivity, in that n is determined by properties of the medium and phenomenon investigated.

An important application arises in a model describing the unsaturated flow of a fluid through a homogeneous porous column under the influence of capillary pressure and gravity. As example, we cite water movement in a vertical column of the medium, if water movement in a horizontal column of the medium, the problem reduces to the flow porous medium.

Use the similarity reduction method in Equation (3), to find the travelling wave

solutions g ( x ± c t ) , where g ( x ± c t ) = x d g ( 1 ± c t x ) , with d is related to the homogeneity degree.

u ( x , t ) = x a U ( η ) , η = 1 ± c t x , (4)

where t > 0 , x ∈ ℝ − { 0 } and a is unknown.

We apply Equation (4) on the left-side hand of Equation (3):

∂ ∂ t u ( x , t ) = ∂ ∂ t [ x a U ( 1 ± c t x ) ] = ( ± c ) x a − 1 ∂ ∂ η U ( η ) . (5)

In relation to the right-side hand of Equation (3). We calculate Riesz fractional derivative, using Theorem 1:

D x α u n ( x , t ) = D x α [ x a n U n ( 1 ± c t x ) ] = d α [ x a n U n ( 1 ± c t x ) ] ∗ | x | − α − 1 = d α ∫ − ∞ + ∞ s a n U n ( 1 ± c t s ) | x − s | − α − 1 d s = d α [ ∫ x + ∞ ( s − x ) − α − 1 s a n U n ( 1 ± c t s ) d s + ∫ − ∞ x ( x − s ) − α − 1 s a n U n ( 1 ± c t s ) d s ] . (6)

The integrals are given follows:

∫ − ∞ x ( x − s ) − α − 1 s a n U n ( 1 ± c t s ) d s = x a n − α ( 1 − η ) 1 + a n ∫ 1 η ( η − τ ) − 1 − α ( 1 − τ ) − 1 + α − a n U n ( τ ) d τ (7)

and

∫ x + ∞ ( s − x ) − α − 1 s a n U n ( 1 ± c t s ) d s = ( − 1 ) − α x a n − α ( 1 − η ) 1 + a n ∫ 1 η ( η − ψ ) − 1 − α ( 1 − ψ ) − 1 + α − a n U n ( ψ ) d ψ (8)

The Riesz fractional derivative is given by:

D x α u n ( x , t ) = d α [ 1 + ( − 1 ) − α ] x a n − α ( 1 − η ) 1 + a n ∫ 1 η ( η − ψ ) − 1 − α ( 1 − ψ ) − 1 + α − a n U n ( ψ ) d ψ . (9)

Replacing, Equation (5) and Equation (9) in Equation (3), we get:

( ± c ) x a − 1 ∂ ∂ η U ( η ) = d α [ 1 + ( − 1 ) − α ] x a n − α ( 1 − η ) 1 + a n ∫ 1 η ( η − ψ ) − 1 − α ( 1 − ψ ) − 1 + α − a n U n ( ψ ) d ψ + x a − 1 [ a + ( 1 − η ) ] ∂ ∂ η U ( η ) (10)

Imposing invariance on the variable x for Equation (10):

a − 1 = a n − α ⇒ a = α − 1 n − 1 (11)

We obtain the equation:

( ± c ) ∂ ∂ η U ( η ) = − d α [ 1 + ( − 1 ) − α ] ( 1 − η ) 1 + a n ∫ η 1 ( η − ψ ) − 1 − α ( 1 − ψ ) − 1 + α − a n U n ( ψ ) d ψ + [ a + ( 1 − η ) ∂ ∂ η ] U ( η ) (12)

Look for a travelling wave, called wavefront [

U ( η ) = { A η k , η < 0 0 , η ≥ 0 (13)

where A and k are unknowns. To determine these variables, we replace this equation in Equation (12). The calculus of the integral is given by:

d α [ 1 + ( − 1 ) − α ] ( 1 − η ) 1 + a n ∫ η 1 ( η − ψ ) − 1 − α ( 1 − ψ ) − 1 + α − a n U n ( ψ ) d ψ = = d α [ 1 + ( − 1 ) − α ] ( 1 − η ) 1 + a n ∫ η 0 ( η − ψ ) − 1 − α ( 1 − ψ ) − 1 + α − a n A n ψ κ n d ψ = d α [ 1 + ( − 1 ) − α ] η κ n − α ( 1 − η ) 1 + a n A n ∫ 1 0 ( 1 − μ ) − 1 − α ( 1 − η μ ) − 1 + α − a n μ κ n d μ = d α [ 1 + ( − 1 ) − α ] A n η k n − α × ( 1 − η ) 1 + a n ( − 1 ) Γ ( − α ) Γ ( 1 + κ n ) Γ ( 1 − α + κ n ) 2 F 1 ( 1 − α + a n , 1 + κ n ; 1 − α + κ n ; η ) , (14)

where 2 F 1 is hypergeometric function. Thus, we obtain the relation:

( ± c ) A k η k − 1 = − d α [ 1 + ( − 1 ) − α ] A n η k n − α × ( 1 − η ) 1 + a n ( − 1 ) Γ ( − α ) Γ ( 1 + k n ) Γ ( 1 − α + k n ) 2 F 1 ( 1 − α + a n , 1 + k n ; 1 − α + k n ; η ) + A [ ( a − k ) η k + k η k − 1 ] . (15)

Once again, considering k = a , we impose invariance on the variable η for Equation (15), that gives us:

( ± c ) k A = − d α [ 1 + ( − 1 ) − α ] A n ( − 1 ) Γ ( − α ) Γ ( 1 + k n ) Γ ( 1 − α + k n ) + k A . (16)

Thus,

A n − 1 = − exp ( i α π / 2 ) Γ ( 1 − α + k n ) Γ ( 1 + k n ) k [ ( ± c ) − 1 ] , (17)

and using Equation (11), we have:

A = { − exp ( i α π / 2 ) Γ ( 1 + k ) Γ ( 1 + k n ) [ ( ± c ) − 1 ] } 1 n − 1 . (18)

The values of A and a in Equation (4),

u ( x , t ) = { { − exp ( i α π / 2 ) Γ ( 1 + k ) Γ ( 1 + k n ) [ ( ± c ) − 1 ] ( x ± c t ) α − 1 } 1 n − 1 , | x | < c t 0 , | x | ≥ c t . (19)

The integer classical equation is recovered by taking the limit α → 2 in Equation (3):

∂ ∂ t u ( x , t ) = ∂ 2 ∂ x 2 u n ( x , t ) + ∂ ∂ x u ( x , t ) , (20)

and its solution is given by:

u ( x , t ) = { { ( n − 1 ) n [ ( ± c ) − 1 ] ( x ± c t ) } 1 n − 1 , | x | < c t 0 , | x | ≥ c t , (21)

which is same result found by B.H. Gilding and R. Kersner in [

The condition η < 0 gives us the relations: x ∈ ( 0 , ∞ ) associated to “ − c ” and x ∈ ( − ∞ , 0 ) associated to“ + c ”, that is, x → − x implies the change of velocity − c → + c . If we regard u ( x , t ) in Equation (21) nonnegative, we have c ≥ 0 x ∈ ( 0 , ∞ ) and c > 1 x ∈ ( − ∞ , 0 ) . Therefore, the fractional equation in Equation (3) has only stationary solution in ( 0 , ∞ ) . The relation between the values A and c is establishes in Equation (18).

We introduce the space-fractional model based on Riesz fractional derivative with nonlinear diffusion in conjunction linear convection. We obtain the explicit travelling wave solution, with finite propagation, wavefront; that is, Equation (3) admits solution in the region x < c t such that u ( x , t ) vanishes on the front x = c t . This property and the behaviour near the interface are in general described by wavefront; it may give rise to interfaces separating the regions that we have just mentioned. As a continuation of this paper, we can suppose fractional derivatives in time and nonlinear convection for Equation (3) and we search its travelling waves.

The authors are grateful to the referee and Dr. A. J. Silva for their useful suggestions and discussions.

Costa, F.S. and Pereira, M.R.A. (2017) Travelling Waves in Space-Fractional Nonlinear Diffusion with Linear Convection. Journal of Applied Mathematics and Physics, 5, 462-468. https://doi.org/10.4236/jamp.2017.52041