Next Article in Journal
A Poisson Algebra for Abelian Yang-Mills Fields on Riemannian Manifolds with Boundary
Next Article in Special Issue
The Solutions of Non-Integer Order Burgers’ Fluid Flowing through a Round Channel with Semi Analytical Technique
Previous Article in Journal
Adaptive Backstepping Sliding Mode Control for the Vertical Launching Barrel-Cover of the Underwater Missile
Previous Article in Special Issue
Numerical Investigation of the Savonius Vertical Axis Wind Turbine and Evaluation of the Effect of the Overlap Parameter in Both Horizontal and Vertical Directions on Its Performance
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Two-Dimensional Advection–Diffusion Process with Memory and Concentrated Source

1
Abdus Salam School of Mathematical Sciences, Government College University, 5400 Lahore, Pakistan
2
Department of Mathematics, Lahore Leads University, 5400 Lahore, Pakistan
3
Department of Theoretical Mechanics, Technical University “Gheorghe Asachi” of Iasi, 6600 Iasi, Romania
*
Author to whom correspondence should be addressed.
Symmetry 2019, 11(7), 879; https://doi.org/10.3390/sym11070879
Submission received: 25 June 2019 / Revised: 2 July 2019 / Accepted: 3 July 2019 / Published: 4 July 2019
(This article belongs to the Special Issue Aero/Hydrodynamics and Symmetry)

Abstract

:
Two-dimensional advection–diffusion processes with memory and a source concentrated in the symmetry center of the domain have been investigated. The differential equation of the studied model is a fractional differential equation with short-tail memory (a differential equation with Caputo–Fabrizio time-fractional derivatives). An analytical solution of the initial-boundary value problem has been determined by employing the Laplace transform and double sine-Fourier transforms. A numerical solution of the studied problem has been determined using finite difference approximations. Numerical simulations for both solutions have been carried out using the software Mathcad.

1. Introduction

Partial differential equations with time-fractional derivatives are suitable to describe complex problems in chaotic dynamics, signal processing, wave propagation, classical and continuum mechanics [1,2,3,4,5,6]. These equations admit the non-local memory effects; therefore, the fractional models can eliminate the disadvantage of classical models that cannot satisfactorily describe some complex processes such as anomalous diffusion.
In recent years, fractional advection–diffusion equations were intensively studied due to their multiple applications in many practical problems such as transport of pollutants in air and water, contaminating streams in groundwater, migration of pollutants in rivers and seawater, and tracer dispersion in porous media.
Bhrawy and Baleanu [7] developed an efficient Legendre–Gauss–Lobatto method to solve the advection–diffusion equation with a space-fractional Caputo derivative and non-homogeneous initial-boundary conditions. The advection–dispersion equation with a fractional order dependent on time and the Coimbra derivative was numerically solved by Abdelkawy et al. [8]. Zhuang et al. [9] numerically solved the advection–diffusion equation with derivatives of variable fractional order and a nonlinear source. They presented the explicit and implicit Euler approximation, fractional method of line, matrix transfer technique and the extrapolation method. Fazio et al. [10] have analyzed the advection–diffusion equation with the time-fractional Caputo derivative by using a finite difference method with non-uniform grids. They obtained an accuracy method compared with the classical discrete fractional formula for the time-fractional Caputo derivative. Time–space fractional advection–diffusion equations with Riesz derivatives have been numerically investigated by Arshad et al. [11]. A fractional advection–diffusion equation with a nonlinear source has been studied by Jannelli et al. [12]. They obtained a numerical solution and an analytical solution based on the Mittag–Leffler functions. Using the finite volume element method, Badr et al. [13] have studied the fractional advection–diffusion equation with the time-fractional Caputo derivative. Stability of the method has been analyzed. Pimenov [14] has numerically investigated the fractional advection–diffusion equation with heredity. Singh et al. [15] have obtained analytical solutions for the fractional advection–dispersion equation by employing the q-fractional homotopy analysis transform method. Using the integral transform method, Avci and Yetim [16] have found the fundamental solutions of the fractional advection–diffusion with fractional derivatives of the Atangana–Baleanu type. Fundamental solutions for the Dirichlet problem of the two-dimensional fractional advection–diffusion equation with short-tail memory have been studied by Mirza and Vieru [17]. Using the Laplace transform and the sine–cosine Fourier transform, Mirza et al. [18] have determined analytical solutions for the fractional advection–diffusion equation in the one-dimensional case and with boundary conditions of Robin-type. Interesting fractional diffusive processes have been investigated by Hristov in [19,20,21]. Povstenko and Klekot [22] and Povstenko and Kyrylych [23] have studied the fractional advection diffusion equation with several boundary conditions. It is important to underline that the time–space fractional advection–diffusion equation could be solved with the Monte–Carlo method as a limit of continuous-time random walks [24].
In this work, the analytical and numerical study of two-dimensional advection–diffusion processes with short-tail memory and with the source concentrated in the symmetry center of the geometric domain is carried out.
The mathematical model is described by a fractional differential equation with the time-fractional Caputo–Fabrizio derivative. An analytical solution of the initial-boundary value problem has been determined by employing the Laplace and finite sine-Fourier transforms. A numerical solution of the studied problem has been determined using finite difference approximations. Numerical simulations for both solutions have been carried out using the software Mathcad and a good agreement is found.

2. The Statement of the Problem

Let us consider the two-dimensional advection–diffusion process with memory and the time-dependent concentrated-source in a finite domain:
D = { ( x , y ) , a x b , c y d , a < b , c < d } with the border D .
The mathematical model is described by the following time-fractional advection–diffusion equation [7,9]:
D t α C F u ( x , y , t ) + k 1 u ( x , y , t ) x + k 2 u ( x , y , t ) y ε 1 2 u ( x , y , t ) x 2 ε 2 2 u ( x , y , t ) y 2 = f ( x , y , t ) , 0 < α 1 ,     k 1 ,    k 2 0 ,     ε 1 ,    ε 2 > 0 ,    ( x , y ) D .
where u ( x , y , t ) is the field variable representing the solute concentration, v = k 1 e x + k 2 e y is the constant fluid velocity (the drift velocity), ε 1 , ε 2 are the dispersion coefficients in the x-direction and y-direction, respectively and f ( x , y , t ) represents the sources/sinks.
The operator D t α C F ( ) denotes the time-fractional Caputo–Fabrizio derivative defined as [25]:
D t α C F u ( x , y , t ) = h ( α , t ) u ˙ ( x , y , t ) , h ( α , t ) = { M ( α ) 1 α exp ( α t 1 α ) δ ( t ) , α = 1 , , 0 α < 1 , u ˙ ( x , y , t ) = u ( x , y , t ) t ,
where “ “ denotes the convolution product and δ ( t ) is Dirac’s distribution. In the above definition, the function M ( α ) is so called the normalization function, which should satisfy conditions M ( 0 ) = M ( 1 ) = 1 . In the present paper, for simplicity, we will consider for the normalization function M ( α ) 1 , α [ 0 , 1 ] .
The following properties of the Caputo–Fabrizio operator immediately result from Equation (2):
D t 0 C F u ( x , y , t ) = u ( x , y , t ) u ( x , y , 0 ) , D t 1 C F u ( x , y , t ) = u ( x , y , t ) t , D t α C F C = 0 , C = c o n s t .
The Laplace transform of the Caputo–Fabrizio derivative is given by:
L { D t α C F u ( x , y , t ) } = L { h ( α , t ) u ˙ ( x , y , t ) } = L { h ( α , t ) } L { u ˙ ( x , y , t ) } = s L { u ( x , t , t ) } u ( x , y , 0 ) ( 1 α ) s + α .
The fractional advection–diffusion equation (Equation (1)), is studied along with following initial and boundary conditions:
u ( x , y , 0 ) = ψ 0 ( x , y ) ,    ( x , y ) D ,
u ( x , y , t ) = ψ 1 ( x , y , t ) ,    t > 0 ,    ( x , y ) D , t > 0 ,
Making the transformation:
u ( x , y , t ) = e f 0 ( x , y ) φ ( x , y , t ) + ψ 0 ( x , y ) ,     f 0 ( x , y ) = k 1 x 2 ε 1 + k 2 y 2 ε 2 ,
Equation (1) becomes:
D t α C F φ ( x , y , t ) + 1 2 ( k 1 2 ε 1 + k 2 2 ε 2 ) φ ( x , y , t ) ε 1 2 φ ( x , y , t ) x 2 ε 2 2 φ ( x , y , t ) y 2 = F ( x , y , t ) + G ( x , y ) ,
where
F ( x , y , t ) = f ( x , y , t ) exp ( f 0 ( x , y ) ) G ( x , y ) = [ ε 1 2 ψ 0 ( x , y ) x 2 + ε 2 2 ψ 0 ( x , y ) y 2 k 1 ψ 0 ( x , y ) x k 2 ψ 0 ( x , y ) y ] exp ( f 0 ( x , y ) ) .
After transformation (7), the initial and boundary conditions become:
φ ( x , y , 0 ) = 0 ,    ( x , y ) D ,
φ ( x , y , t ) = F 1 ( x , y , t ) + G 1 ( x , y ) , ( x , y ) D , t > 0 ,
where
F 1 ( x , y , t ) = ψ 1 ( x , y , t ) exp ( f 0 ( x , y ) ) , G 1 ( x , y ) = ψ 0 ( x , y ) exp ( f 0 ( x , y ) ) .

3. Analytical Solution of the Problem

The solution of the problem described by Equation (8) and conditions (10) and (11) will be determined by using the integral transform method.
Applying the Laplace transform to Equation (8), using Equation (4) and Equation (10), we obtain the transformed equation:
[ 1 + β ( 1 α ) ] s + α β ( 1 α ) s + α φ ¯ ( x , y , s ) ε 1 2 φ ¯ ( x , y , s ) x 2 ε 2 2 φ ¯ ( x , y , s ) y 2 = F ¯ ( x , y , s ) + 1 s G ( x , y ) ,
where β = 1 2 ( k 1 2 ε 1 + k 2 2 ε 2 ) and φ ¯ ( x , y , s ) = 0 φ ( x , y , t ) e s t dt , F ¯ ( x , y , s ) = 0 F ( x , y , t ) e s t dt are the Laplace transforms of functions φ ( x , y , t ) , F ( x , y , t ) , respectively.
Equation (13) has to satisfy the boundary condition:
φ ¯ ( x , y , s ) = φ ¯ 1 ( x , y , s ) , ( x , y ) D , φ ¯ 1 ( x , y , s ) = F ¯ 1 ( x , y , s ) + 1 s G 1 ( x , y ) .
Now, we apply this to Equation (13), the double sine-Fourier transform with respect to variables x and y defined as [26]:
φ ¯ ˜ ( n , k , s ) = c d a b sin [ n ( x a ) π b a ] sin [ k ( y c ) π d c ] φ ¯ ( x , y , s ) d x d y ,     n ,    k = 1 , 2 ,
A direct calculation leads to the following relationships:
c d a b 2 φ ¯ ( x , y , s ) x 2 sin [ n ( x a ) π b a ] sin [ k ( y c ) π d c ] d x d y = H ¯ ˜ 1 ( n , k , s ) ( n π b a ) 2 φ ¯ ˜ ( n , k , s ) ,
where
H ¯ ˜ 1 ( n , k , s ) = n π b a c d [ φ ¯ 1 ( a , y , s ) ( 1 ) n φ ¯ 1 ( b , y , s ) ] sin [ k ( y c ) π d c ] d y
a b c d 2 φ ¯ ( x , y , s ) y 2 sin [ n ( x a ) π b a ] sin [ k ( y c ) π d c ] d y d x = H ¯ ˜ 2 ( n , k , s ) ( k π d c ) 2 φ ¯ ˜ ( n , k , s )
where
H ¯ ˜ 2 ( n , k , s ) = k π d c a b [ φ ¯ 1 ( x , c , s ) ( 1 ) k φ ¯ 1 ( x , d , s ) ] sin [ n ( x a ) π b a ] d x .
Applying the double sine-Fourier transform to Equation (13), using the relationships (16)–(19) and notations
a n k = 1 + ( 1 α ) [ β + ε 1 ( n π b a ) 2 + ε 2 ( k π d c ) 2 ] b n k = α [ β + ε 1 ( n π b a ) 2 + ε 2 ( k π d c ) 2 ] , n , k = 1 , 2 , , Ψ ¯ ˜ ( n , k , s ) = F ¯ ˜ ( n , k , s ) + 1 s G ˜ ( n , k ) + ε 1 H ¯ ˜ 1 ( n , k , s ) + ε 2 H ¯ ˜ 2 ( n , k , s ) ,
the transform φ ¯ ˜ ( n , k , s ) is written in the following form:
φ ¯ ˜ ( n , k , s ) = ( 1 α ) s + α a n k s + b n k Ψ ¯ ˜ ( n , k , s ) .
Now, we will write Equation (21) in a suitable form which, after inversion, satisfies the initial and boundary conditions (10) and (11). To do this, we consider the following auxiliary functions:
v i ( y , t ) = [ g i ( y , t ) g i ( d , t ) ] H ( y c ) + [ g i ( y , t ) g i ( c , t ) ] H ( d y ) g i ( y , t ) H ( y c ) H ( d y ) ;      i = 1 , 2 ,     y [ c , d ] ,
w i ( x , t ) = [ h i ( x , t ) h i ( b , t ) ] H ( x a ) + [ h i ( x , t ) h i ( a , t ) ] H ( b x )                h i ( x , t ) H ( x a ) H ( b x ) ;      i = 1 , 2 ,     x [ a , b ] ,
where
g 1 ( y , t ) = F 1 ( a , y , t ) + G 1 ( a , y ) ,    g 2 ( y , t ) = F 1 ( b , y , t ) + G 1 ( b , y ) , h 1 ( x , t ) = F 1 ( x , c , t ) + G 1 ( x , c ) ,    h 2 ( x , t ) = F 1 ( x , d , t ) + G 1 ( x , d ) ,
and
H ( z ) = sign ( z ) [ 1 + sign ( z ) ] 2 = { 0 ,        z 0 1 ,        z > 0 ,
is the Heaviside step unit function.
With the above functions, we define:
θ 1 ( x , y , t ) = b x b a v 1 ( y , t ) + x a b a v 2 ( y , t ) = x [ v 2 ( y , t ) v 1 ( y , t ) ] b a + b v 1 ( y , t ) a v 2 ( y , t ) b a ,
θ 2 ( x , y , t ) = d y d c w 1 ( x , t ) + y c d c w 2 ( x , t ) = y [ w 2 ( x , t ) w 1 ( x , t ) ] d c + d w 1 ( x , t ) c w 2 ( x , t ) d c ,
θ 3 ( x , y , t ) = h 1 ( a , t ) H ( b x ) H ( d y ) + h 1 ( b , t ) H ( x a ) H ( d y ) + h 2 ( a , t ) H ( b x ) H ( y c ) + h 2 ( b , t ) H ( x a ) H ( y c ) ,
The function θ ( x , y , t ) defined by
θ ( x , y , t ) = θ 1 ( x , y , t ) + θ 2 ( x , y , t ) + θ 3 ( x , y , t ) ,
satisfies the following properties:
θ ( a , y , t ) = g 1 ( y , t ) ,     θ ( b , y , t ) = g 2 ( y , t ) ,     θ ( x , c , t ) = h 1 ( x , t ) ,     θ ( x , d , t ) = h 2 ( x , t )
Applying the Laplace transform and double sine-Fourier transform to functions (26)–(28), we obtain:
θ ¯ ˜ 1 ( n , k , s ) = a ( 1 ) n b n π P 1 k ( s ) + 1 ( 1 ) n n π P 2 k ( s ) , P 1 k ( s ) = c d [ v ¯ 2 ( y , s ) v ¯ 1 ( y , s ) ] sin [ k ( y c ) π d c ] d y , P 2 k = c d [ b v ¯ 1 ( y . s ) a v ¯ 2 ( y , s ) ] sin [ k ( y c ) π d c ] d y ,
θ ¯ ˜ 2 ( n , k , s ) = c ( 1 ) k d k π Q 1 n ( s ) + 1 ( 1 ) k k π Q 2 n ( s ) , Q 1 n ( s ) = a b [ w ¯ 2 ( x , s ) w ¯ 1 ( x , s ) ] sin [ n ( x a ) π b a ] d x , Q 2 n ( s ) = a b [ d w ¯ 1 ( x , s ) c w ¯ 2 ( x , s ) ] sin [ n ( x a ) π b a ] d x ,
θ ¯ ˜ 3 ( n , k , s ) = [ h ¯ 1 ( a , s ) + h ¯ 1 ( b , s ) + h ¯ 2 ( a , s ) + h ¯ 2 ( b , s ) ] ( 1 ( 1 ) n n π ) ( 1 ( 1 ) k k π ) ( b a ) ( d c )
respectively,
θ ¯ ˜ ( n , k , s ) = a ( 1 ) n b n π P 1 k ( s ) + 1 ( 1 ) n n π P 2 k ( s ) + c ( 1 ) k d k π Q 1 n ( s ) + 1 ( 1 ) k k π Q 2 n ( s ) +                    [ h ¯ 1 ( a , s ) + h ¯ 1 ( b , s ) + h ¯ 2 ( a , s ) + h ¯ 2 ( b , s ) ] ( 1 ( 1 ) n n π ) ( 1 ( 1 ) k k π ) ( b a ) ( d c ) .
The function φ ¯ ˜ , given by Equation (21), can be written as:
φ ¯ ˜ ( n , k , s ) = θ ¯ ˜ ( n , k , s ) + 1 a n k s + b n k Ψ ¯ ˜ 1 ( n , k , s ) ,
where
Ψ ¯ ˜ 1 ( n , k , s ) = [ ( 1 α ) s + α ] Ψ ¯ ˜ ( n , k , s ) ( a n k s + b n k ) θ ¯ ˜ ( n , k , s ) .
Applying the inverse Laplace and Fourier transforms, we obtain the solution
φ ( x , y , t ) = H ( t ) θ ( x , y , t ) + 4 ( b a ) ( d c ) n = 1 k = 1 sin [ n ( x a ) π b a ] sin [ k ( y c ) π d c ] 0 t exp ( ( t τ ) b n k / a n k ) a n k Ψ ˜ 1 ( n , k , τ ) d τ .
Using the property (30), it is easy to see that function (37) satisfies boundary condition (11). The solution for the classical advection equation is obtained by making α = 1 into Equations (20) and (35)–(37).

4. Particular Case: A Time-Dependent Concentrated Source in the Center of Domain D

In this section, for the source f ( x , y , t ) and for the initial and boundary conditions, we consider:
f ( x , y , t ) = γ 0 δ ( x a + b 2 ) δ ( y c + d 2 ) exp ( p t ) , p 0 , u ( x , y , 0 ) = ψ 0 ( x , y ) = γ 0 exp ( f 0 ( x , y ) ) , ( x , y ) D , γ 0 = c o n s t . , u ( x , y , t ) = ψ 1 ( x , y , t ) = 2 γ 0 exp ( f 0 ( x , y ) ) , ( x , y ) D , t > 0 .
The physical significance of the above source f ( x , y , t ) is that of a source of intensity γ 0 exp ( p t ) concentrated in the center ( a + b 2 , c + d 2 ) of the rectangular domain D.
Using (38), Equations (9), (11), (12) and (14) lead to:
F ( x , y , t ) = γ 0 δ ( x a + b 2 ) δ ( y c + d 2 ) exp ( p t f 0 ( x , y ) ) , G ( x , y ) = β γ 0 2 , F 1 ( x , y , t ) = 2 γ 0 ,    G 1 ( x , y ) = γ 0 ,    φ 1 ( x , y , t ) = γ 0 .
By using φ ¯ 1 ( x , y , s ) = γ 0 s in Equations (17) and (19), we have:
H ¯ ˜ 1 ( n , k , s ) = γ 0 s d c b a n [ 1 ( 1 ) n ] [ 1 ( 1 ) k ] k , H ¯ ˜ 2 ( n , k , s ) = γ 0 s b a d c k [ 1 ( 1 ) n ] [ 1 ( 1 ) k ] n .
Using the “sifting” property of the Dirac distribution [27]
a b f ( x ) δ ( x c ) d x = { 1 2 f ( c + ) , c = a , 1 2 [ f ( c ) + f ( c + ) ] , c ( a , b ) , 1 2 f ( c ) , c = b , 0 , c ( a , b ) .
we find that the double sine-Fourier transform of the functions
F ¯ ( x , y , s ) = γ 0 s + p δ ( x a + b 2 ) δ ( y c + d 2 ) exp ( f 0 ( x , y ) ) , G ( x , y ) = β γ 0 2 ,
are given by
F ¯ ˜ ( n , k , s ) = γ 0 s + p sin ( n π 2 ) sin ( k π 2 ) exp ( k 1 ( a + b ) 2 ε 1 2 k 2 ( c + d ) 2 ε 2 2 ) , G ˜ ( n , k ) = β γ 0 2 ( b a ) ( d c ) k n π 2 [ 1 ( 1 ) n ] [ 1 ( 1 ) k ] .
Since sin ( m π 2 ) = ( 1 ) m 1 2 i m + 1 , i 2 = 1 , m , the function F ¯ ˜ ( n , k , s ) is written as
F ¯ ˜ ( n , k , s ) = γ 0 s + p exp ( k 1 ( a + b ) 2 ε 1 2 k 2 ( c + d ) 2 ε 2 2 ) [ 1 ( 1 ) n ] [ 1 ( 1 ) k ] 4 i n + k .
By using Equations (40) and (43) in the third equation (20), we obtain:
Ψ ¯ ˜ ( n , k , s ) = [ 1 ( 1 ) n [ [ 1 ( 1 ) k ] 4 × [ 2 γ 0 s 2 ε 1 ( n π ) 2 ( d c ) 2 + 2 ε 2 ( k π ) 2 ( b a ) 2 β ( b a ) 2 ( d c ) 2 ( k π ) ( n π ) ( b a ) ( d c ) γ 0 γ 1 i n + k s + p ] ,
where
γ 1 = exp ( k 1 ( a + b ) 4 ε 1 k 2 ( c + d ) 4 ε 2 )
In the considered particular case, Equations (22)–(24) and (31)–(34) become:
g 1 ( y , t ) = g 2 ( y , t ) = γ 0 , v 1 ( y , t ) = v 2 ( y , t ) = γ 0 H ( y c ) H ( d y ) , h 1 ( x , t ) = h 2 ( x , t ) = γ 0 , w 1 ( x , t ) = w 2 ( x , t ) = γ 0 H ( x a ) H ( b x ) ,
θ 1 ( x , y , t ) = γ 0 H ( y c ) H ( d y ) , θ 2 ( x , y , t ) = γ 0 H ( x a ) H ( b x ) , θ 3 ( x , y , t ) = γ 0 H ( b x ) H ( d y ) + γ 0 H ( x a ) H ( d y ) + γ 0 H ( b x ) H ( y c ) + γ 0 H ( x a ) H ( y c ) .
The transformed Laplace–Fourier of the function (47) are given by:
θ ¯ ˜ 1 ( n , k , s ) = θ ¯ ˜ 2 ( n , k , s ) = 1 s γ 0 ( b a ) ( d c ) 1 ( 1 ) n n π 1 ( 1 ) k k π , θ ¯ ˜ 3 ( n , k , s ) = 4 s γ 0 ( b a ) ( d c ) 1 ( 1 ) n n π 1 ( 1 ) k k π , θ ¯ ˜ ( n , k , s ) = 2 s γ 0 ( b a ) ( d c ) 1 ( 1 ) n n π 1 ( 1 ) k k π .
Function Ψ ¯ ˜ 1 ( n , k , s ) given by Equation (36) can be written as:
Ψ ¯ ˜ 1 ( n , k , s ) = C n k 1 s + p + D n k 1 s + E n k ,
where
A n k = [ 1 ( 1 ) n ] [ 1 ( 1 ) k ] B n k = 1 ( b a ) ( d c ) [ 2 ε 1 ( n π ) 2 ( d c ) 2 + 2 ε 2 ( k π ) 2 ( b a ) 2 β ( b a ) 2 ( d c ) 2 ] , C n k = 1 4 γ 0 γ 1 [ ( 1 α ) p α ] i n + k A n k , D n k = γ 0 A n k 2 ( n π ) ( k π ) [ α B n k 4 ( b a ) ( d c ) b n k ] , E n k = A n k 4 ( n π ) ( k π ) [ 2 ( 1 α ) γ 0 B n k 8 γ 0 0 ( b a ) ( d c ) a n k ( 1 α ) γ 0 γ 1 ( n π ) ( k π ) i n + k ] ,
and transform (35) becomes:
φ ¯ ˜ ( n , k , s ) = θ ¯ ˜ ( n , k , s ) + 1 a n k 1 s + b n k / a n k [ C n k s + p + D n k s + E n k ]
Applying the inverse Laplace transform and the inverse Fourier transform to function (51), we obtain:
φ ( x , y , t ) = H ( t ) θ ( x , y , t ) + 4 ( b a ) ( d c ) n = 1 k = 1 1 a n k { E n k exp ( b n k t a n k ) +       0 t exp ( b n k ( t τ ) a n k ) ( D n k + C n k e p τ ) d τ } sin ( n ( x a ) π b a ) sin ( k ( y c ) π d c ) ,
where the function θ ( x , y , t ) is given by
θ ( x , y , t ) = γ 0 [ H ( x a ) + H ( b x ) ] [ H ( y c ) + H ( d y ) ] γ 0 H ( x a ) H ( b x ) γ 0 H ( y c ) H ( d y ) .
Using the property
ϑ ( t ) = H ( t α 1 ) + H ( α 2 t ) H ( t α 1 ) H ( α 2 t ) 1 , t [ α 1 , α 2 ] ,
it is found that θ ( a , y , t ) = θ ( b , y , t ) = γ 0 for y [ c , d ] , respectively θ ( x , c , t ) = θ ( b , d , t ) = γ 0 for x [ a , b ] ; therefore, function (53) satisfied boundary condition (14). Now, to return to the concentration u ( x , y , t ) , the function φ ( x , y , t ) given by Equation (52) will be replaced into transformation (7).
Figure 1, Figure 2, Figure 3 and Figure 4 have been plotted to highlight the variation of the non-dimensional concentration φ ( x , y , t ) versus the fractional parameter α .
To draw these graphs we have used the following values of the parameters: a = 0.5 , b = 1.0 , c = 1.0 , d = 1.5 , p = 0.25 , γ 0 = 4 , ε 1 = ε 2 = 0.2 , k 1 = 0.25 , k 2 = 0.5 .
In Figure 1, we show the plotted curves φ ( x , 1.2 , 0.5 ) and φ ( x , 1.2 , 1.0 ) for different values of the fractional parameter α { 0.3 , 0.5 , 0.7 , 1.0 } . It is observed from Figure 1 that the values of concentration φ decrease with the fractional parameter and have a peak (maximum/minimum) in the central area of the geometrical domain D = [ a , b ] × [ c , d ] .
In Figure 2 and Figure 3, we show the plotted surfaces φ ( x , y , 1 ) for the fractional parameters α = 0.3 and α = 0.5 , respectively, along with the sections through planes y = 1.2 , x = 0.74 .
For small values of the fractional parameter, the concentration φ has maximum values in the central area of the rectangular D = [ a , b ] × [ c , d ] . If the fractional parameter increases, the concentration has a minimum in the same area. This behavior is due to the kernel of the time-fractional Caputo–Fabrizio derivative (see Figure 4) which, for t = 0.5 , increases for α 0.55 and decreases for α > 0.55 . Therefore, at small values of the fractional parameter, the weight function in the fractional derivative has bigger values and the influence on the concentration is stronger.

5. Numerical Solution

In this section, based on the finite difference approximation of the derivatives, we develop a numerical scheme to the problem formulated in Section 2 of the paper.
We consider Equation (8) along with the initial and boundary conditions (10) and (11), with ( x , y , t ) [ a , b ] × [ c , d ] × [ 0 , T ] , T > 0 . The discrete sets of nodes are defined as:
x i = a + i h x , i = 0 , 1 , , N 1 , h x = ( b a ) / N 1 , y j = c + j h y , j = 0 , 1 , , N 2 , h y = ( d c ) / N 2 , t n = n h t , n = 0 , 1 , , N , h t = T / N .
The numerical estimation of the solution φ ( x , y , t ) of Equation (8) in the point ( x i , y j , t n ) is denoted by φ i , j n = φ ( x i , y j , t n ) . The initial and boundary conditions (10) and (11) are written as:
φ i , j 0 = φ ( x i , y j , t 0 ) = 0 , i = 0 , 1 , , N 1 , j = 0 , 1 , , N 2 , φ 0 , j n = φ ( a , y j , t n ) = F 1 ( a , y j , t n ) + G 1 ( a , y j ) = P 0 , j n , j = 0 , 1 , , N 2 , n = 1 , 2 , , N , φ N 1 , j n = φ ( b , y j , t n ) = F 1 ( b , y j , t n ) + G 1 ( b , y j ) = Q N 1 , j n , j = 0 , 1 , , N 2 , n = 1 , 2 , , N , φ i , 0 n = φ ( x i , c , t n ) = F 1 ( x i , c , t n ) + G 1 ( x i , c ) = R i , 0 n , i = 0 , 1 , , N 1 , n = 1 , 2 , , N , φ i , N 2 n = φ ( x i , d , t n ) = F 1 ( x i , d , t n ) + G 1 ( x i , d ) = S i , N 2 n , i = 0 , 1 , , N 1 , n = 1 , 2 , , N ,
Using the following approximation for the first time-derivative:
φ ( x i , y j , t ) t | t = t k = 1 h t [ φ ( x i , y j , t k + 1 ) φ ( x i , y j , t k ) ] = 1 h t ( φ i , j k + 1 φ i , j k ) , k = 0 , 1 , , N 1 ,
the time-fractional Caputo–Fabrizio derivative, in the point t = t n , can be approximated by:
D t α C F φ ( x , y , t ) | ( x i , y j , t n ) = 1 1 α 0 t n φ ( x i , y j , t ) t | t = τ exp ( α ( t n τ ) 1 α ) d τ = k = 0 n 1 1 1 α t k t k + 1 φ ( x i , y j , t ) t | t = τ exp ( α ( t n τ ) 1 α ) d τ = k = 0 n 1 1 ( 1 α ) h t ( φ i , j k + 1 φ i , j k ) t k t k + 1 exp ( α ( t n τ ) 1 α ) d τ = k = 0 n 1 a n , k ( φ i , j k + 1 φ i , j k ) ,    n = 1 , 2 , N ,
where
a n , k = 1 α h t [ exp ( α ( t n t k + 1 ) 1 α ) exp ( α ( t n t k ) 1 α ) ] , α ( 0 , 1 ) , n = 1 , 2 , , N , k = 0 , 1 , , ( n 1 )
It is known [28,29] that Equation (58) approximates the time-fractional Caputo–Fabrizio derivative with an error of order O ( h t 2 ) .
The space derivatives of the second order are approximated by:
2 φ ( x , y , t ) x 2 | ( x i , y j , t n ) = 1 h x 2 ( φ i + 1 , j n 2 φ i , j n + φ i 1 , j n ) , i = 1 , 2 , , N 1 1 , j = 0 , 1 , , N 2 , n = 1 , 2 , N , 2 φ ( x , y , t ) y 2 | ( x i , y j , t n ) = 1 h y 2 ( φ i , j + 1 n 2 φ i , j n + φ i , j 1 n ) , i = 0 , 1 , , N 1 , j = 1 , 2 , , N 2 1 , n = 1 , 2 , N .
Using the approximation formulae (58) and (60), Equation (8) is written as:
k = 0 n 1 a n , k ( φ i , j k + 1 φ i , j k ) + β φ i , j n + β 1 ( φ i + 1 , j n 2 φ i , j n + φ i 1 , j n ) + β 2 ( φ i , j + 1 n 2 φ i , j n + φ i , j 1 n ) = H i , j n , i = 1 , 2 , , N 1 1 , j = 1 , 2 , , N 2 1 , n = 1 , 2 , N ,
where
β 1 = ε 1 h x 2 , β 2 = ε 2 h y 2 , H i , j n = F ( x i , y j , t n ) + G ( x i , y j )
For n = 1 , Equation (61) becomes:
d 1 , 0 φ i , j 1 + β 1 ( φ 1 + 1 , j 1 + φ i 1 , j 1 ) + β 2 ( φ i , j + 1 1 + φ i , j 1 1 ) = H i , j 1 ,
where d 1 , 0 = a 1 , 0 + β 2 ( β 1 + β 2 ) .
Making j = 1 and i = 1 , 2 , , N 1 1 into Equation (63), we obtain the following algebraic system:
M X 1 + M 0 X 2 = D 1 1 ,
where M is a square matrix of order ( N 1 1 ) whose elements are given by:
M 0 , k = d 1 , 0 δ 0 , k + β 1 δ 1 , k , k = 0 , 1 , , N 1 2 , M i , k = β 1 δ i , k + 1 + d 1 , 0 δ i + 1 , k + 1 + β 1 δ i + 2 , k + 1 , i = 1 , 2 , , N 1 2 , k = k = 0 , 1 , , N 1 2 ,
so,
M = ( d 1 , 0 β 1 0 0 0 0 0 0 β 1 d 1 , 0 β 1 0 0 0 0 0 0 β 1 d 1 , 0 β 1 0 0 0 0 0 0 0 0 0 β 1 d 1 , 0 β 1 0 0 0 0 0 0 β 1 d 1 , 0 )
the matrix M 2 is the diagonal square matrix of order ( N 1 1 ) with elements M 1 i , j = β 2 δ i , j , i , j = 0 , 1 , , N 1 2 , δ i , j = { 0 , i j , 1 , i = j , is the Kronecker delta and X 1 , X 2 D 1 1 are the column-matrices.
X 1 = ( φ 1 , 1 1 φ 2 , 1 1 φ N 1 1 , 1 1 ) ,    X 2 = ( φ 1 , 2 1 φ 2 , 2 1 φ N 1 1 , 2 1 ) ,    D 1 1 = ( D 11 1 D 12 1 D 1 N 1 1 1 ) , D 11 1 = H 1 , 1 1 β 1 P 0 , 1 1 β 2 R 1 , 0 1 , D 1 i 1 = H i , 1 1 β 2 R i , 0 1 , i = 2 , 3 , , N 1 2 , D 1 N 1 1 1 = H N 1 1 , 1 1 β 1 Q N 1 , 1 1 β 2 R N 1 1 , 0 1 .
The algebraic system (64) has ( N 1 1 ) equations with 2 ( N 1 1 ) unknown functions.
Making j = 2 , 3 , , N 2 2 and i = 1 , 2 , , N 1 1 into Equation (63), we obtain the algebraic system:
M 0 X j 1 + M X j + M 0 X j + 1 = D j 1 ,
where
D j 1 1 = H 1 , j 1 β 1 P 0 , j 1 , D j i 1 = H i , j 1 , i = 2 , 3 , , N 1 2 , D j N 1 1 1 = H N 1 1 . j 1 β 1 Q N 1 , j 1 .
Making j = N 2 1 and i = 1 , 2 , , N 1 1 into Equation (63), we obtain the algebraic system:
M 0 X N 2 2 + M X N 2 1 = D N 2 1 1 ,
where
D N 2 11 1 = H 1 , N 2 1 1 β 1 P 0 , N 2 1 1 β 2 S 1 , N 2 1 , D N 2 1 i 1 = H i , N 2 1 1 β 2 S i , N 2 1 , i = 2 , 3 , , N 1 2 , D N 2 1 N 1 1 1 = H N 1 1 . N 2 1 1 β 1 Q N 1 , N 2 1 1 β 2 S N 1 1 , N 2 1 .
The algebraic systems (64), (68) and (70) can be written as a unique algebraic system of ( N 1 1 ) ( N 2 1 ) equations with ( N 1 1 ) ( N 2 1 ) unknown functions φ i , j 1 , i = 1 , 2 , , N 1 1 , j = 1 , 2 , , N 2 1 , namely
( M M 0 0 0 0 0 0 0 M 0 M M 0 0 0 0 0 0 0 M 0 M M 0 0 0 0 0 0 0 0 0 0 M 0 M M 0 0 0 0 0 0 0 M 0 M ) ( X 1 X 2 X 3 X N 2 2 X N 2 1 ) = ( D 1 1 D 2 1 D 3 1 D N 2 2 1 D N 2 1 1 )
The solution of system (72) gives the values of function φ ( x , y , t ) in points ( x i , y j , t 1 ) for i = 1 , 2 , , N 1 1 , j = 1 , 2 , , N 2 1 .
The numerical procedure is continued in a similar way for n = 2 , 3 , , N when Equation (61) is written in the equivalent form:
d n , n 1 φ i , j n + β 1 ( φ i + 1 , j n + φ i 1 , j n ) + β 2 ( φ i , j + 1 n + φ i , j 1 n ) = G i , j n , i = 1 , 2 , , N 1 1 , j = 1 , 2 , , N 2 1 ,
where
d n , n 1 = a n , n 1 + β 2 ( β 1 + β 2 ) , n = 2 , 3 , , N , G i , j n = H i , j n + a n , n 1 φ i , j n 1 k = 0 n 2 a n , k ( φ i , j k + 1 φ i , j k ) , i = 1 , 2 , , N 1 1 , j = 1 , 2 , , N 2 1 .
Finally, the proposed numerical procedure gives the values of the solution φ ( x , y , t ) in grid points ( x i , y j , t n ) , i = 1 , 2 , , N 1 1 , j = 1 , 2 , , N 2 1 , n = 1 , 2 , , N .
To compare the numerical values obtained with the analytical solution, respectively, by the numerical method, we have considered the following particular problem:
F ( x , y , t ) = 5000 exp ( k 1 x 2 ε 1 k 2 y 2 ε 2 ) , G ( x , y ) = 0 , ( x , y , t ) [ 0 , 0.05 ] × [ 0 , 0.05 ] × [ 0 , 0.1 ] , F 1 ( x , y , t ) = 0 , G 1 ( x , y ) = 0 , φ ( x , y , 0 ) = 0 , k 1 = k 2 = 0.5 , ε 1 = ε 2 = 2 ,
In this case, the double Fourier transform of the function F ( x , y , t ) is:
F ˜ ( n , k ) = n π ( b a ) [ exp ( a k 1 2 ε 1 ) ( 1 ) n exp ( b k 1 2 ε 1 ) ] ( k 1 2 ε 1 ) 2 ( b a ) 2 + ( n π ) 2 × k π ( b a ) [ exp ( c k 2 2 ε 2 ) ( 1 ) k exp ( d k 2 2 ε 2 ) ] ( k 2 2 ε 2 ) 2 ( d c ) 2 + ( k π ) 2 ,
and the analytical solution is given by:
φ ( x , y , t ) = 4 α ( b a ) ( d c ) n = 1 k = 1 F ˜ ( n , k ) b n k sin [ n ( x a ) π b a ] sin [ k ( y c ) π d c ] + 4 ( b a ) ( d c ) n = 1 k = 1 F ˜ ( n , k ) [ ( 1 α ) b n k α a n k ] a n k b n k sin [ n ( x a ) π b a ] sin [ k ( y c ) π d c ] e b n k t a n k .
The values of function φ ( x , y , t ) at the instant t = 0.01 , obtained with the numerical method and with expression (77) given in the Table 1, are in good agreement.

6. Discussion and Conclusions

A two-dimensional advection–diffusion process with short-tail memory and concentrated source situated in the symmetry center of the geometrical domain has been investigated.
The studied mathematical model is described by a fractional differential equation with a time-fractional Caputo–Fabrizio derivative with exponential kernel.
Analytical solutions of the proposed initial-boundary value problem were obtained using the integral transform method (Laplace and finite Fourier transforms are employed).
Numerical solutions of the studied problem have been determined using finite difference approximations.
Numerical simulations for both analytical and numerical solutions have been carried out using the software Mathcad.
The fractional order of the Caputo–Fabrizio derivative has a significant influence on the pollutant concentration. For small values of the fractional parameter, the concentration has maximum values in the central area of the geometrical domain. If the fractional parameter increases, the concentration has a minimum in the same area. This behavior is due to the memory kernel of the time-fractional Caputo–Fabrizio derivative.

Author Contributions

The authors have equally contributed to this paper.

Funding

This research received no external funding.

Conflicts of Interest

On behalf of all authors, the corresponding author states that there is no conflict of interest.

References

  1. Mainardi, F. Fractals and Fractional Calculus in Continuum Mechanics; Springer: Wien, Austria, 1997. [Google Scholar]
  2. Laskin, N. Fractional Quantum Mechanics; World Scientific: Hackensack, NJ, USA, 2018. [Google Scholar]
  3. Laskin, N. Generalized classical mechanics. Eur. Phys. J. Spec. Top. 2013, 222, 1929–1938. [Google Scholar] [CrossRef]
  4. Hristov, J. Fractional derivative with non-singular kernels: From the Caputo-Fabrizio definition and beyond: Appraising analysis with emphasis on diffusion models. In Frontiers in Fractional Calculus; Bhalekar, S., Ed.; Bentham Publication: Kolhapur, India, 2017; Chapter 10; pp. 243–269. [Google Scholar]
  5. Baleanu, D.; Golmankhaneh, A.K.; Nigmatullin, R.; Golmankhaneh, A.K. Fractional Newtonian mechanics. Cent. Eur. J. Phys. 2010, 8, 120–125. [Google Scholar] [CrossRef]
  6. Povstenko, Y. Fractional Thermoelasticity, Solid Mechanics and Its Applications; Springer International Publishing: Cham, Switzerland, 2015; p. 219. [Google Scholar]
  7. Bhrawy, A.H.; Baleanu, D. A spectral Legendre-Gauss-Lobatto collocation method for a space-fractional advection diffusion equation with variable coefficients. Rep. Math. Phys. 2013, 72, 219–233. [Google Scholar] [CrossRef]
  8. Abdelkawy, M.A.; Zaky, M.A.; Bhrawy, A.H.; Baleanu, D. Numerical simulation of time variable fractional order mobile-immobile advection-dispersion model. Rom. Rep. Phys. 2015, 67, 773–791. [Google Scholar]
  9. Zhuang, P.; Liu, F.; Anh, V.; Turner, I. Numerical methods for the variable-order fractional advection-diffusion equation with a nonlinear source term. SIAM J. Num. Anal. 2009, 47, 1760–1781. [Google Scholar] [CrossRef]
  10. Fazio, R.; Jannelli, A.; Agreste, S. A finit difference method on non-uniform meshes for time-fractional advection-diffusion equations with a source term. Appl. Sci. 2018, 8, 960. [Google Scholar] [CrossRef]
  11. Arshad, S.; Baleanu, D.; Huang, J.; Al Qurashi, M.M.; Tang, Y.; Zhao, Y. Finite difference method for time-space fractional advection-diffusion equations with Riesz derivative. Entropy 2018, 20, 321. [Google Scholar] [CrossRef]
  12. Jannelli, A.; Ruggieri, M.; Speciale, M.P. Exact and numerical solutions of time-fractional advection-diffusion equation with a nonlinear source term by means of the Lie symmetries. Nonlinear Dyn. 2018, 92, 543–555. [Google Scholar] [CrossRef]
  13. Badr, M.; Yazdani, A.; Jafari, H. Stability of a finite volume element method for the time-fractional advection-diffusion equation. Numer. Methods Partial Differ. Eqs. 2018, 34, 1459–1471. [Google Scholar] [CrossRef]
  14. Pimenov, V.G. Numerical method for advection-diffusion equation with heredity. Itogi Naukii Tekniki Ser. Sovrem. Mat. Pril. Temat. Obz. 2017, 132, 86–90. [Google Scholar] [CrossRef]
  15. Singh, J.; Secer, A.; Swroop, R.; Kumar, D. A reliable analytical approach for a fractional model of advection-dispersion equation. Nonlinear Eng. 2018, in press. [Google Scholar] [CrossRef]
  16. Avci, D.; Yetim, A. Analytical solution to the advection-diffusion equation with the Atangana-Baleanu derivative over a finite domain. J. BAUN Inst. Sci. Technol. 2018, 20, 382–395. [Google Scholar] [CrossRef]
  17. Mirza, I.A.; Vieru, D. Fundamental solutions to advection-diffusion equation with time-fractional Caputo-Fabrizio derivative. Comput. Math. Appl. 2017, 73, 1–10. [Google Scholar] [CrossRef]
  18. Mirza, I.A.; Vieru, D.; Ahmed, N. Fractional advection-diffusion equation with memory and Robin-type boundary condition. Math. Model. Nat. Phenom. 2019, 14, 306. [Google Scholar] [CrossRef]
  19. Hristov, J. Integral balance approach to 1-D space-fractional diffusion models. In Mathematical Methods in Engineering Applications in Dynamics of Complex Systems; Tas, K., Baleanu, D., Machado, J.A.T., Eds.; Springer: Cham, Switzerland, 2019; pp. 1–21. [Google Scholar]
  20. Hristov, J. The heat radiation diffusion equation with memory: Constitutive approach and approximate integral-balance solutions. In Heat Conduction, Methods, Applications and Research; Hristov, J., Bennacer, R., Eds.; NOVA Science Publishers: New York, NY, USA, 2019. [Google Scholar]
  21. Hristov, J. Response functions in linear viscoelastic constitutive equations and related fractional operators. Math. Model. Nat. Phenom. 2018, 14, 305. [Google Scholar] [CrossRef]
  22. Povstenko, Y.; Klekot, J. The Dirichlet problem for the time-fractional advection-diffusion equation in a line segment. Bound. Value Probl. 2016, 2016, 89. [Google Scholar] [CrossRef]
  23. Povstenko, Y.; Kyrylych, T. Time-fractional diffusion with mass absorption in a half-space line domain due to boundary value of concentration varying harmonically in time. Entropy 2018, 19, 346. [Google Scholar] [CrossRef]
  24. Fulger, D.; Scalas, E.; Germano, G. Monte Carlo simulation of uncoupled continuous-time random walks yielding a stochastic solution of the space-time fractional diffusion equation. Phys. Rev. E 2008, 77, 021122. [Google Scholar] [CrossRef] [Green Version]
  25. Caputo, M.; Fabrizio, M. A new definition of fractional derivative without singular kernel. Progr. Fract. Differ. Appl. 2015, 1, 73–85. [Google Scholar]
  26. Andrews, L.C.; Shivamoggi, B.K. Integral Transforms for Engineers; SPIE Press: Washington, DC, USA, 1999. [Google Scholar]
  27. Bracewell, R. The Shifting Property. In The Fourier Transform and Its Applications, 3rd ed.; McGraw-Hill Science/Engineering/Math.: New York, NY, USA, 1999; pp. 74–77. [Google Scholar]
  28. Atangana, A.; Alqahtani, R.T. Numerical approximation of the space-time Caputo-Fabrizio fractional derivative and application to groundwater pollution equation. Adv. Diff. Eqs. 2016, 2016, 156. [Google Scholar] [CrossRef]
  29. Rangaig, N.A. Finite difference approximation for Caputo-Fabrizio time fractional derivative on non-uniform mesh and some applications. Phys. J. 2018, 3, 255–263. [Google Scholar]
Figure 1. Profiles of function φ ( x , 1.2 , t ) for different values of the fractional parameter α, y = 1.2, t = 0.5 and t = 1.0.
Figure 1. Profiles of function φ ( x , 1.2 , t ) for different values of the fractional parameter α, y = 1.2, t = 0.5 and t = 1.0.
Symmetry 11 00879 g001
Figure 2. Profiles of surface φ ( x , y , 1 ) for α = 0.3 and their sections with the planes y = 1.2 and x = 0.74 , respectively.
Figure 2. Profiles of surface φ ( x , y , 1 ) for α = 0.3 and their sections with the planes y = 1.2 and x = 0.74 , respectively.
Symmetry 11 00879 g002
Figure 3. Profiles of surface φ ( x , y , 1 ) for α = 0.5 and their sections with the planes y = 1.2 and x = 0.74 , respectively.
Figure 3. Profiles of surface φ ( x , y , 1 ) for α = 0.5 and their sections with the planes y = 1.2 and x = 0.74 , respectively.
Symmetry 11 00879 g003
Figure 4. Profile of kernel h ( α , t ) for two values of the time t.
Figure 4. Profile of kernel h ( α , t ) for two values of the time t.
Symmetry 11 00879 g004
Table 1. Values of function φ ( x , y , 0.01 ) .
Table 1. Values of function φ ( x , y , 0.01 ) .
x y 1 = 0.005 y 2 = 0.010 y 3 = 0.015
Numerical
Solution
Equation (77)Numerical
Solution
Equation (77)Numerical
Solution
Equation (77)
0.0050.0800.0810.1290.1300.1580.159
0.0100.1290.1300.2140.2160.2670.270
0.0150.1580.1590.2670.2700.3380.340
0.0200.1740.1750.2970.2990.3790.381
0.0250.1790.1810.3070.3090.3910.394
0.0300.1740.1750.2970.2990.3780.381
0.0350.1580.1590.2680.2690.3380.340
0.0400.1250.1200.2130.2150.2670.270
0.0450.0800.0810.1290.1300.1580.159
x y 4 = 0.020 y 5 = 0.025 y 6 = 0.030
Numerical
Solution
Equation (77)Numerical
Solution
Equation (77)Numerical
Solution
Equation (77)
0.0050.1740.1750.1790.1810.1740.175
0.0100.2970.2990.3060.3090.2970.299
0.0150.3780.3800.3910.3940.3790.381
0.0200.4240.4270.4390.4410.4240.427
0.0250.4390.4410.4550.4570.4380.442
0.0300.4230.4260.4380.4410.4230.426
0.0350.3770.3800.3910.3930.3770.380
0.0400.2960.2980.3060.3080.2960.299
0.0450.1740.1750.1790.1810.1740.175
x y 7 = 0.035 y 8 = 0.040 y 9 = 0.045
Numerical
Solution
Equation (77)Numerical
Solution
Equation (77)Numerical
Solution
Equation (77)
0.0050.1580.1590.1280.1300.0800.081
0.0100.2670.2690.2130.2150.1280.130
0.0150.3380.3400.2670.2690.1580.159
0.0200.3770.3800.2960.2980.1740.175
0.0250.3900.3930.3060.3080.1790.180
0.0300.3770.3800.2960.2990.1740.175
0.0350.3370.3400.2660.2690.1570.159
0.0400.2660.2690.2120.2140.1280.129
0.0450.1570.1590.1280.1290.0790.081

Share and Cite

MDPI and ACS Style

Ahmed, N.; Shah, N.A.; Vieru, D. Two-Dimensional Advection–Diffusion Process with Memory and Concentrated Source. Symmetry 2019, 11, 879. https://doi.org/10.3390/sym11070879

AMA Style

Ahmed N, Shah NA, Vieru D. Two-Dimensional Advection–Diffusion Process with Memory and Concentrated Source. Symmetry. 2019; 11(7):879. https://doi.org/10.3390/sym11070879

Chicago/Turabian Style

Ahmed, Najma, Nehad Ali Shah, and Dumitru Vieru. 2019. "Two-Dimensional Advection–Diffusion Process with Memory and Concentrated Source" Symmetry 11, no. 7: 879. https://doi.org/10.3390/sym11070879

APA Style

Ahmed, N., Shah, N. A., & Vieru, D. (2019). Two-Dimensional Advection–Diffusion Process with Memory and Concentrated Source. Symmetry, 11(7), 879. https://doi.org/10.3390/sym11070879

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