Next Article in Journal
Interfacial Dynamics of Miscible Displacement of Shear-Thinning Fluid in a Vertical Channel
Next Article in Special Issue
Investigation of Flow-Induced Instabilities in a Francis Turbine Operating in Non-Cavitating and Cavitating Part-Load Conditions
Previous Article in Journal
Acknowledgment to the Reviewers of Fluids in 2022
Previous Article in Special Issue
Assessment of a RANS Transition Model with Flapping Foils at Moderate Reynolds Numbers
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Generalized Diffusion Equation: Solutions and Anomalous Diffusion

1
Departamento de Física, Universidade Estadual de Ponta Grossa, Av. Gen. Carlos Cavalcanti 4748, Ponta Grossa 84030-900, PR, Brazil
2
National Institute of Science and Technology for Complex Systems, Centro Brasileiro de Pesquisas Físicas, Rio de Janeiro 22290-180, RJ, Brazil
3
Departamento de Física, Universidade Tecnológica Federal do Paraná, R. Marcílio Dias, 635, Apucarana 86812-460, PR, Brazil
4
Departamento de Física, Universidade Federal do Rio Grande do Norte, Natal 59078-900, RN, Brazil
5
Departamento de Engenharia Química, Universidade Federal do Paraná, Curitiba 80060-000, PR, Brazil
*
Author to whom correspondence should be addressed.
Fluids 2023, 8(2), 34; https://doi.org/10.3390/fluids8020034
Submission received: 30 November 2022 / Revised: 7 January 2023 / Accepted: 10 January 2023 / Published: 17 January 2023
(This article belongs to the Special Issue Recent Advances in Fluid Mechanics: Feature Papers, 2022)

Abstract

:
We investigate the solutions of a generalized diffusion-like equation by considering a spatial and time fractional derivative and the presence of non-local terms, which can be related to reaction or adsorption–desorption processes. We use the Green function approach to obtain solutions and evaluate the spreading of the system to show a rich class of behaviors. We also connect the results obtained with the anomalous diffusion processes.

Graphical Abstract

1. Introduction

Fractional calculus has quickly become a new efficient mathematical tool to analyze different properties of a given system and connect them with experimental results. A simple extension of the differential operators incorporating non-integer indexes has serious consequences, connecting the formalism with memory effects, long-range correlations, and many other features characterizing complex systems. In this manner, it has brought great insights into many fields of science [1,2,3,4,5,6]. One of them occurs for the diffusion processes, which have been found through fractional calculus a suitable approach to incorporate several effects which are not suitably described in terms of the classical integer order calculus. For instance, infiltration in porous building materials [7], the electrical response of electrolytic cells [6,8], amorphous semiconductors [9], and micellar solutions [10]. In these situations, there is a nonlinear time dependence exhibited by the mean-square displacement, which, in general, is characterized by x x 2 t α , where α characterized the diffusion, e.g., α < 1 , α = 1 , and α > 1 correspond to the sub-, usual, and superdiffusion, respectively. This behavior of the mean square displacement and effects related, e.g., non-Markovian processes and fractal structure has motivated the analysis of different approaches, which extend the usual approach, such as fractional diffusion equations [11,12,13], master equation [14,15], generalized Langevin equations [16], and random walks [17]. One noticeable point regarding these extensions is that the behavior of the solutions is characterized by power laws and stretched exponential for fractional differential operators. In particular, the previous scenarios concerning the diffusion on fractals have indicated that the asymptotic form of the propagator, such as the Sierpinski gasket, is essentially in this form [11,18]. A similar situation is also found in fluid flow through porous media [19,20,21]. Other applications can be found in transport in the porous pellet [22], transport of chloride in concrete [23], and oxidation behaviors of needle-punched carbon/carbon composites [24].
Here, we consider the following extension of the diffusion equation:
t ρ ( x , t ) = 0 D t γ 0 μ d μ ¯ p ( μ ¯ ) D x μ ¯ , η ρ ( x , t ) + 0 t d t Λ ( t t ) ρ ( x , t ) ,
with 1 μ 2 , 1 < η , p ( μ ¯ ) is a distribution, 0 D t γ is a fractional time operator, and D x μ ¯ , η ( ) is a spatial fractional operator [12,25]. The last term can be related to different processes, such as absorption and adsorption–desorption, depending on the choice of the kernel Λ ( t ) . It can also be related to reaction processes of the first order, i.e., irreversible reaction processes or reversible processes depending on the choice of Λ ( t ) . The fractional time operator 0 D t γ ( ) is defined in terms of a generalized kernel K γ ( t ) as follows:
0 D t γ ρ ( x , t ) = t 0 t d t K γ ( t t ) ρ ( x , t ) .
Note that depending on the choice of K γ ( t ) in the previous equation, we may connect it to different integrodifferential operators with singular or non-singular kernels. One of them is the Riemann–Liouville fractional operator, i.e., K γ ¯ ( t ) = t γ ¯ / Γ 1 γ ¯ [26], which implies
γ ¯ t γ ¯ ρ ( x , t ) = 1 Γ 1 γ ¯ t 0 t d t 1 ( t t ) γ ¯ ρ ( x , t ) ,
another one is the Fabrizio–Caputo fractional operator, for K γ ¯ ( t ) = N γ ¯ e γ ¯ 1 γ ¯ t / ( 1 γ ¯ ) , i.e.,
γ ¯ t γ ¯ ρ ( x , t ) = N γ ¯ 1 γ ¯ t 0 t d t e γ ¯ 1 γ ¯ ( t t ) ρ ( x , t ) ,
or the Atangana–Baleanu fractional operator, for K γ ¯ ( t ) = N γ ¯ E γ ¯ γ ¯ t γ ¯ / ( 1 γ ¯ ) / ( 1 γ ¯ ) , given by
γ ¯ t γ ¯ ρ ( x , t ) = N γ ¯ 1 γ ¯ t 0 t d t E γ ¯ γ ¯ 1 γ ¯ ( t t ) γ ¯ ρ ( x , t ) ,
where N γ ¯ is a normalization constant [27,28,29]. These operators may be related to different scenarios in connection with anomalous diffusion, which implies memory effects, long-range correlation, and fractal structures, among others. In particular, these fractional operators have been used in many contexts such as boundary value problems [30,31], electric circuits [32,33], and electrical impedance [34,35] (see also Refs. [36,37,38]). It is also possible to consider other kernels for Equation (2), with different implications for the relaxation processes (see, e.g., Refs. [39,40]).
Following the developments performed in Ref [12,25], the spatial operator is defined as follows:
1 2 d x ψ ± , η ( x , k ) D x μ , η ρ ( x , t ) | k | μ + η ρ ˜ ± ( k , t ) ,
with the integral transform given by:
1 2 d x ψ ± , η ( x , k ) ρ ( x , t ) = ρ ˜ ± ( k , t ) ,
1 2 d k ψ ± , η ( x , k ) ρ ˜ ± ( k , t ) = ρ ( x , t ) ,
where
ψ + , η ( x , k ) = | k | | x | 1 2 ( 1 + η ) J ν 2 | k | | x | 1 2 ( 2 + η ) / ( 2 + η ) and
ψ , η ( x , k ) = x k | k | | x | 1 2 ( 1 + η ) 1 J ν 2 | k | | x | 1 2 ( 2 + η ) / ( 2 + η ) ,
where ψ + , η ( x , k ) and ψ , η ( x , k ) refer to the odd and even functions, ν = ( 1 + η ) / ( 2 + η ) , and J ν ( x ) is the Bessel function [6]. Equations (7) and (8) may be related to a generalized Hankel transform [41,42,43,44] and for μ = 2 , we can directly relate the fractional operator present in Equation (6) with standard differential operators as follows:
D x 2 , η ( ) x | x | η x ( ) .
This case allows us to relate Equation (6) with a diffusion process in heterogeneous media. Such behaviors are also exhibited by diffusion-related problems, such as diffusion on fractals [45,46], turbulence [47,48], diffusion and reaction on fractals [49], and solute transport in fractal porous media [50], where the properties of the media promote an anomalous diffusion. In these scenarios, we have non-Gaussian distributions related to these processes and nonlinear behavior of the mean square displacement. Equation (11) also allows us to connect Equation (1), for p ( μ ) = δ ( μ 2 ) , directly with the continuity equation with an additional term as follows:
t ρ ( x , t ) + x J ( x , t ) = 0 t d t Λ ( t t ) ρ ( x , t ) ,
with the current density given by
J ( x , t ) = 0 D t γ | x | η x ρ ( x , t ) .
Notice that | x | η x ( ) ( 1 + η ) | x | 1 + η ( ) , i.e., it corresponds to a fractal derivative [51,52,53], respectively. Thus, the spatial fractional operator defined above by Equation (6) can be considered a mixing between the Riesz–Wely operator [11] and fractal operators [54]. This feature implies that the solutions of Equation (6) can be related to Lévy distributions and/or distributions with characteristics of stretched exponential. From the above discussion, Equation (1) has a particular case of several situations analyzed in different scenarios and allows analyzing the mixing between different effects on the diffusion process connected to these fractional operators. Further, the reaction term may be connected to different processes, particularly the stochastic resetting [55,56].
We aim to analyze Equation (1) by considering different scenarios. The first considers the absence of the non-local term, i.e., Λ ( t ) = 0 . For this case, we obtain solutions by considering different spatial and time fractional operator choices. In particular, we also consider the case p ( μ ¯ ) = χ μ δ ( μ ¯ μ ) + χ 2 δ ( μ ¯ 2 ) , where χ μ and χ 2 are constants. After we incorporate the non-local term, i.e., Λ ( t ) 0 , in our analysis. We obtain exact solutions in the framework of Green’s function approach for all cases. These formal developments are shown in Section 2. Section 3 discusses the main results and concludes with some remarks.

2. Fractional Dynamics and Diffusion

Let us start our discussion about Equation (1) by establishing the boundary conditions, which are lim x ± ρ ( x , t ) = 0 . An arbitrary function represents the initial condition, i.e., ρ ( x , 0 ) = φ ( x ) . After establishing the boundary and the initial condition, Equation (1) in the absence of the non-local term is given by
t ρ ( x , t ) = 0 D t γ 0 μ d μ ¯ p ( μ ¯ ) D x μ ¯ , η ρ ( x , t ) .
Formally, we can write the solution for this equation as follows:
ρ ( x , t ) = d x G ( x , x , t ) φ ( x ) ,
where the Green’s function [57], G ( x , x , t ) , satisfies the following equation:
t G ( x , x , t ) 0 D t γ 0 μ d μ ¯ p ( μ ¯ ) D x μ ¯ , η G ( x , x , t ) = δ ( x x ) δ ( t ) .
The Green’s function is subjected to the following conditions: lim x ± G ( x , x , t ) = 0 and G ( x , x , t ) = 0 for t < 0 . In terms of the Equations (9) and (10), it is possible to write Green’s function as
G ( x , x , t ) = 1 2 d k ψ + , η ( x , k ) G ˜ + ( k , x , t ) + ψ , η ( x , k ) G ˜ ( k , x , t ) .
Equation (16) can be simplified by using the integral transform defined by Equations (7) and (8), yielding
t G ˜ ± ( k , x , t ) + 0 D t γ 0 μ d μ ¯ p ( μ ¯ ) | k | μ + η G ˜ ± ( k , x , t ) = 1 2 ψ ± , η ( x , k ) δ ( t ) ,
which, after applying the Laplace transform, has the solution given by
G ˜ ^ ± ( k , x , s ) = ψ ± , η ( x , k ) 2 s 1 + K γ ( s ) 0 μ d μ ¯ p ( μ ¯ ) | k | μ ¯ + η .
The inverse Laplace transform of Equation (18) depends on the choice of the fractional time operator, i.e., the option performed to K γ ( t ) .
Let us start with the case K γ ( s ) = 1 / s ( K γ ( t ) = c o n s t ) with p ( μ ¯ ) = χ μ δ ( μ ¯ μ ) . Applying these conditions in Equation (19), we have that
G ˜ ± ( k , x , t ) = 1 2 ψ ± , η ( x , k ) e χ μ | k | μ + η t .
By applying the inverse of the integral transform, it is possible to show that
G ( x , x , t ) = G + ( x , x , t ) + x x | x | | x | G ( x , x , t ) ,
with
G ± ( x , x , t ) = ( 2 + η ) 2 2 ( μ + η ) | x | 2 + η | x | | x | 1 2 1 + η × H 2 , [ 0 : 1 ] , 0 , [ 0 : 2 ] 1 , 0 , 1 , 1 , 1 | x | | x | 2 + η Υ ¯ μ , η t 2 + η μ + η | x | 2 + η 2 ν 2 , 1 ; 2 ± ν 2 , 1 ; ( 0 , 2 + η μ + η ) ; ν 2 , 1 ; ± ν 2 , 1 ; ( 0 , 2 + η μ + η ) , ( 0 , 2 + η μ + η ) ,
where Υ ¯ μ , η = ( 2 + η ) 2 χ μ 2 + η μ + η and
H E , [ A : C ] , F , [ B , D ] L , M , M 1 , N , N 1 x y ( ε 1 , ω 1 ) , , ( ε E , ω E ) ( a 1 , α 1 ) , , ( a A , α A ) ; ( c 1 , β 1 ) , , ( c C , β C ) ( ξ 1 , ϖ 1 ) , , ( ξ F , ϖ F ) ( b 1 , β 1 ) , , ( b B , β B ) ; ( d 1 , δ 1 ) , , ( d C , δ D )
is the generalized H function of Fox [58,59,60].
For the particular case μ = 2 , it is possible to simplify Equation (22) and to show that
G ( x , x , t ) = 1 ( 2 + η ) χ μ t | x | | x | 1 2 1 + η e 1 ( 2 + η ) 2 χ μ t | x | 2 + η + | x | 2 + η × I ν 2 | x | | x | 1 2 ( 2 + η ) ( 2 + η ) 2 χ μ t + sgn ( x x ) I ν 2 | x | | x | 1 2 ( 2 + η ) ( 2 + η ) 2 χ μ t ,
where I ν ( x ) is the Bessel function of order ν of modified argument [57]. Note that the mean square displacement for this particular case is given by σ x 2 = x x 2 t 2 / ( 2 + η ) evidencing the anomalous behavior of the system with the time evolution.
Now, we extend the previous case for K γ ( s ) = 1 / s γ , which corresponds to considering the Riemann–Liouville fractional time derivative. For this case, with p ( μ ¯ ) = χ μ δ ( μ ¯ μ ) , we have that
G ˜ ± ( k , x , t ) = 1 2 ψ ± , η ( x , k ) E γ χ μ | k | μ + η t γ ,
where E γ ( x ) is the Mittag–Leffler function, defined as follows:
E γ ( x ) = n = 0 x n Γ 1 + γ n .
An interesting point concerning the Mittag–Leffler function is the asymptotic behavior governed by a power law instead of an exponential. This feature has implications for the random process connected to these functions, characterized by waiting time distributions with long-tailed behavior. This point can be verified by using the continuous time random walk approach, for example, in Ref. [11], which connects the fractional diffusion equations with a random process. The solution for this case is also given by Equation (21) with
G ± ( x , x , t ) = ( 2 + η ) 2 2 ( μ + η ) | x | 2 + η | x | | x | 1 2 1 + η × H 2 , [ 0 : 1 ] , 0 , [ 0 : 2 ] 1 , 0 , 1 , 1 , 1 | x | | x | 2 + η Υ ¯ μ , η t γ 2 + η μ + η | x | 2 + η 2 ν 2 , 1 ; 2 ± ν 2 , 1 ; ( 0 , 2 + η μ + η ) ; ν 2 , 1 ; ± ν 2 , 1 ; ( 0 , 2 + η μ + η ) , ( 0 , γ 2 + η μ + η ) .
Equation (27) mixing three different parameters, which in connection with a random walk, can be connected to the waiting time and jumping distributions. The parameter γ related to the fractional time derivative is connected to the waiting time distribution, and the other parameters η and μ are connected to the jumping probability. We can have a short or a long-tailed behavior for the spatial distribution depending on the choice of the parameter η and μ . Figure 1 and Figure 2 show the behavior of G + ( x , x , t ) and G ( x , x , t ) for different values of the parameters γ , η , and μ . In particular, it is possible to observe that depending on the choice of the parameters, the previous Green functions can exhibit a unimodal or a bimodal behavior. Figure 3 shows the behavior of the Green function G ( x , x , t ) given by Equations (21) and (27). It is worth mentioning that different choices for the parameters imply different behaviors obtained with the mixing of different behaviors.
In addition, by using the scaling arguments [13], it is possible to show that the solution can be written as ρ ( x , t ) = ρ ξ / t γ / ( μ + η ) with ξ = | x | / t γ / ( μ + η ) and, consequently, σ x 2 t 2 γ / ( μ + η ) for d ξ ξ 2 ρ ( ξ ) finite. This result implies that for 2 γ / ( μ + η ) less, equal, or greater than one, we have subdiffusion, usual diffusion, or superdiffusion, respectively (see Figure 4).
Now, we consider the case p ( μ ¯ ) = χ 2 δ ( μ ¯ 2 ) + χ μ δ ( μ ¯ μ ) . This case can be connected to the mixing of two different behaviors. For this case, we have that Equation (19) can be written as follows:
G ˜ ^ ± ( k , x , s ) = 1 2 ψ ± , η ( x , k ) s + s 1 γ χ 2 | k | 2 + η + χ μ | k | μ ¯ + η .
Applying the inverse of the Laplace transform, we obtain that
G ˜ ± ( k , x , t ) = n = 0 ( 1 ) n Γ ( 1 + n ) χ μ | k | μ + η t γ n E γ ( n ) χ 2 | k | 2 + η t γ = n = 0 ( 1 ) n Γ ( 1 + n ) χ μ | k | μ + η t γ n H 1 , 2 1 , 1 χ 2 | k | 2 + η t γ ( 0 , 1 ) , ( γ n , γ ) 0 , 1 .
Performing the inverse of the integral transform, we obtain that
G ± ( x , x , t ) = 2 + η 2 | x | 2 + η | x | | x | 1 2 1 + η n = 0 ( 1 ) n Γ ( 1 + n ) χ ¯ 2 , μ t 2 μ μ + η γ n × H 2 , [ 0 : 1 ] , 0 , [ 0 : 2 ] 1 , 0 , 1 , 1 , 1 | x | | x | 2 + η χ ¯ 2 t γ | x | 2 + η 2 ν 2 , 1 ; 2 ± ν 2 , 1 ; ( 2 + η μ + η n , 1 ) ; ν 2 , 1 ; ± ν 2 , 1 ; 2 + η μ + η n , 1 , 2 μ 2 + η n γ , γ ,
where χ ¯ 2 , μ = χ μ / χ 2 ( μ + η ) / ( 2 + η ) and χ ¯ 2 = ( 2 + η ) 2 χ 2 .
Figure 5 shows the behavior of the Green function for the previous case, which considers the mixing between two different differential operators.
For the initial condition given by ρ ( x , 0 ) = δ ( x ) , the solution is illustrated in Figure 6 for different values of the parameters γ , μ , and η . We also illustrate the behavior of 1 / ρ ( 0 , t ) 2 , which is connected to the mean square displacement related to this case. In particular, from Figure 6b, it is possible to observe the presence of different diffusion regimes depending on the choice of the parameters.
We can consider a different fractional time operator for the previous scenarios defined by Equation (14). One of them is the Caputo–Fabrizio fractional operator [28]. In this case, Equation (14) can be written as follows:
t ρ ( x , t ) = 0 μ d μ ¯ p ( μ ¯ ) D x μ ¯ , η ρ ( x , t ) + r δ ( x x ) ρ ( x , t ) ,
with r = γ / ( 1 γ ) , for the initial condition ρ ( x , 0 ) = δ ( x x ) . It is worth mentioning that Equation (31) corresponds to a system subjected to a stochastic resetting [61]. In particular, it extends the processes described in Ref. [56]. By using the previous approach, it is possible to show that the solution for this case is given by
ρ ( x , t ) = e r t G ( x , x , t ) + r 0 t d t e r t G ( x , x , t )
with the Green function given by Equations (21) and (27). For the particular case, p ( μ ¯ ) = χ μ δ ( μ ¯ μ ) with x = 0 , we have that the Green function is given by
G ( x , 0 , t ) = 2 + η 2 ( μ + η ) | x | Γ 3 + η 2 + η H 2 , 3 2 , 1 [ | x | 2 + η ( 2 + η ) 2 χ μ t γ 2 + μ 2 + η 1 2 + η , 1 , 1 , 2 + η μ + η , 1 , 1 1 , 2 + η μ + η , 1 , 2 + η μ + η ] .
Figure 7a shows the behavior of Equation (32) and Figure 7b shows the behavior of 1 / [ ρ ( 0 , t ) ] 2 to illustrate the spreading of the system for different values of the parameters η and μ with γ = 1 / 2 . From this figure, we observe that the system has long been stationary. It is also interesting to mention that Equation (32) with the Green function given by Equation (33) allows us to investigate a diffusive process in heterogeneous media with stochastic resetting.
Now, we consider the presence of a non-local term in the diffusion equation. We also consider, for simplicity, the initial condition ρ ( x , 0 ) = δ ( x ) , p ( μ ¯ ) = χ μ δ ( μ ¯ μ ) , and the Riemann–Liouville fractional time operator. For this case, we can apply the previous procedure based on integral transforms, yielding
ρ ˜ ^ + ( k , s ) = 1 2 ψ + ( 0 , k ) s + s 1 γ χ μ | k | μ + Λ ( s ) .
After performing a series of expansions, we have that
ρ ˜ ^ + ( k , s ) = 1 2 ψ + ( 0 , k ) s + s 1 γ χ μ | k | μ n = 0 ( 1 ) n Λ ( s ) s + s 1 γ χ μ | k | μ n ,
which, after performing the inverse integral transforms, yields
ρ ( x , t ) = G + ( x , 0 , t ) + n = 1 1 2 n 0 t d t n d x n Υ ( x , x n , t t n ) × 0 t 2 d t 1 d x 2 Υ ( x 2 , x 1 , t 2 t 1 ) Υ ( x 1 , 0 , t 1 )
where Υ ( x , ξ , t ) = 0 t d t G + ( x , ξ , t t ) Λ ( t ) .

3. Discussion and Conclusions

We have investigated a generalized diffusion equation which has, in particular cases, several situations. We have started our analysis by considering the fractional spatial operator and analyzed the influence of the fractional time operators on the solutions. In this scenario, we have obtained the time behavior of the mean square displacement by using scaling arguments when fractional space and time derivatives are present in the diffusion equation. In this case, we consider a singular kernel for the fractional time derivative that allows a connection with the Riemann–Liouville fractional time derivative. For the spatial fractional operator, we have also considered an operator of distributed order. In particular, we analyzed the mixing between two cases, i.e., p ( μ ¯ ) = χ μ δ ( μ ¯ μ ) + χ 2 δ ( μ 2 ) . In each case, the solutions can be directly connected to the stretched exponential or power laws, depending on the choice of the parameters characterizing the spatial fractional operator. We have also considered Fabrizio–Caputo fractional time operator. For this case, we have related this case with a stochastic resetting process following the approach presented in Ref. [28] and analyzed the behavior of the solutions. For each case, the solutions were obtained using the Green function approach. In addition, we have also considered the solutions for an arbitrary non-local term in the generalized diffusion equation. Finally, we hope that the results found here can be helpful in the discussion of different scenarios in connection with diffusion and anomalous diffusion processes.

Author Contributions

Conceptualization, E.K.L., A.S., R.S.Z., L.R.d.S. and M.K.L.; methodology, E.K.L., A.S., R.S.Z., L.R.d.S. and M.K.L.; formal analysis, E.K.L., A.S., R.S.Z., L.R.d.S. and M.K.L.; investigation, E.K.L., A.S., R.S.Z., L.R.d.S. and M.K.L.; writing—original draft preparation, E.K.L., A.S., R.S.Z., L.R.d.S. and M.K.L.; writing—review and editing, E.K.L., A.S., R.S.Z., L.R.d.S. and M.K.L. All authors have read and agreed to the published version of the manuscript.

Funding

E.K.L. acknowledges the support of the CNPq (Grant No. 302983/2018-0) and the National Institute of Science and Technology of Complex Systems-INCT-SC. R.S.Z thanks CNPq (304634/2020-4) and the National Institute of Science and Technology of Complex Fluids-INCT-FCx. Research developed with the support of LAMAP-UTFPR.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

We acknowledge the support given by CNPq.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Hilfer, R. Applications of Fractional Calculus in Physics; World Scientific: Singapore, 2000. [Google Scholar]
  2. Meerschaert, M.M.; Sikorskii, A. Stochastic Models for Fractional Calculus; de Gruyter: Berlin, Germany, 2019. [Google Scholar]
  3. Tarasov, V.E. Fractional Dynamics: Applications of Fractional Calculus to Dynamics of Particles, Fields and Media; Springer Science & Business Media: Berlin, Germany, 2011. [Google Scholar]
  4. Magin, R. Fractional calculus in bioengineering, part 1. Crit. Rev. Biomed. Eng. 2004, 32, 1–104. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  5. Herrmann, R. Fractional Calculus: An Introduction for Physicists; World Scientific: Singapore, 2011. [Google Scholar]
  6. Evangelista, L.R.; Lenzi, E.K. Fractional Diffusion Equations and Anomalous Diffusion; Cambridge University Press: Cambridge, UK, 2018. [Google Scholar]
  7. Kuntz, M.; Lavallée, P. Experimental evidence and theoretical analysis of anomalous diffusion during water infiltration in porous building materials. J. Phys. D Appl. Phys. 2001, 34, 2547. [Google Scholar] [CrossRef]
  8. Rosseto, M.P.; Evangelista, L.R.; Lenzi, E.K.; Zola, R.S.; Ribeiro de Almeida, R.R. Frequency-Dependent Dielectric Permittivity in Poisson–Nernst–Planck Model. J. Phys. Chem. B 2022, 126, 6446–6453. [Google Scholar] [CrossRef]
  9. Scher, H.; Montroll, E.W. Anomalous transit-time dispersion in amorphous solids. Phys. Rev. B 1975, 12, 2455. [Google Scholar] [CrossRef]
  10. Jeon, J.H.; Leijnse, N.; Oddershede, L.B.; Metzler, R. Anomalous diffusion and power-law relaxation of the time averaged mean squared displacement in worm-like micellar solutions. New J. Phys. 2013, 15, 045011. [Google Scholar] [CrossRef]
  11. Metzler, R.; Klafter, J. The random walk’s guide to anomalous diffusion: A fractional dynamics approach. Phys. Rep. 2000, 339, 1–77. [Google Scholar] [CrossRef]
  12. Lenzi, E.K.; Evangelista, L.; Zola, R.; Scarfone, A. Fractional Schrödinger equation for heterogeneous media and Lévy-like distributions. Chaos Solitons Fractals 2022, 163, 112564. [Google Scholar] [CrossRef]
  13. Magin, R.L.; Lenzi, E.K. Slices of the Anomalous Phase Cube Depict Regions of Sub-and Super-Diffusion in the Fractional Diffusion Equation. Mathematics 2021, 9, 1481. [Google Scholar] [CrossRef]
  14. Kenkre, V.; Montroll, E.; Shlesinger, M. Generalized master equations for continuous-time random walks. J. Stat. Phys. 1973, 9, 45–50. [Google Scholar] [CrossRef]
  15. Swenson, R.J. Derivation of generalized master equations. J. Math. Phys. 1962, 3, 1017–1022. [Google Scholar] [CrossRef]
  16. Cortes, E.; West, B.J.; Lindenberg, K. On the generalized Langevin equation: Classical and quantum mechanicala. J. Chem. Phys. 1985, 82, 2708–2717. [Google Scholar] [CrossRef]
  17. Klafter, J.; Sokolov, I.M. First Steps in Random Walks: From Tools to Applications; OUP Oxford: Oxford, UK, 2011. [Google Scholar]
  18. Giona, M.; Roman, H.E. Fractional diffusion equation on fractals: One-dimensional case and asymptotic behaviour. J. Phys. A Math. Gen. 1992, 25, 2093. [Google Scholar] [CrossRef]
  19. Hashan, M.; Jahan, L.N.; Imtiaz, S.; Hossain, M.E. Modelling of fluid flow through porous media using memory approach: A review. Math. Comput. Simul. 2020, 177, 643–673. [Google Scholar] [CrossRef]
  20. Razminia, K.; Razminia, A.; Baleanu, D. Fractal-fractional modelling of partially penetrating wells. Chaos Solitons Fractals 2019, 119, 135–142. [Google Scholar] [CrossRef]
  21. Raghavan, R.; Chen, C. The Theis solution for subdiffusive flow in rocks. Oil Gas Sci. Technol. Rev. D’Ifp Energies Nouv. 2019, 74, 6. [Google Scholar] [CrossRef]
  22. Zhokh, A.; Strizhak, P. Macroscale modeling the methanol anomalous transport in the porous pellet using the time-fractional diffusion and fractional Brownian motion: A model comparison. Commun. Nonlinear Sci. Numer. Simul. 2019, 79, 104922. [Google Scholar] [CrossRef]
  23. Feng, C.; Si, X.; Li, B.; Cao, L.; Zhu, J. An inverse problem to simulate the transport of chloride in concrete by time–space fractional diffusion model. Inverse Probl. Sci. Eng. 2021, 29, 2429–2445. [Google Scholar] [CrossRef]
  24. Han, M.; Zhou, C.; Silberschmidt, V.V.; Bi, Q. Multiscale heat conduction and fractal oxidation behaviors of needle-punched carbon/carbon composites. Sci. Eng. Compos. Mater. 2022, 29, 508–515. [Google Scholar] [CrossRef]
  25. Lenzi, E.K.; Evangelista, L. Space–time fractional diffusion equations in d-dimensions. J. Math. Phys. 2021, 62, 083304. [Google Scholar] [CrossRef]
  26. Podlubny, I. Fractional Differential Equations; Academic Press: Cambridge, MA, USA, 1999. [Google Scholar]
  27. Atangana, A.; Baleanu, D. New Fractional Derivatives with non-local and non-singular kernel. Therm. Sci. 2016, 20, 763–769. [Google Scholar] [CrossRef]
  28. Tateishi, A.A.; Ribeiro, H.V.; Lenzi, E.K. The role of fractional time-derivative operators on anomalous diffusion. Front. Phys. 2017, 5, 52. [Google Scholar] [CrossRef] [Green Version]
  29. Fernandez, A.; Baleanu, D. Classes of operators in fractional calculus: A case study. Math. Methods Appl. Sci. 2021, 44, 9143–9162. [Google Scholar] [CrossRef]
  30. Singh, H. Chebyshev spectral method for solving a class of local and nonlocal elliptic boundary value problems. Int. J. Nonlinear Sci. Numer. Simul. 2021. [Google Scholar] [CrossRef]
  31. Singh, H. Solving a class of local and nonlocal elliptic boundary value problems arising in heat transfer. Heat Transf. 2022, 51, 1524–1542. [Google Scholar] [CrossRef]
  32. Singh, H. An efficient computational method for non-linear fractional Lienard equation arising in oscillating circuits. In Methods of Mathematical Modelling; CRC Press: Boca Raton, FL, USA, 2019; pp. 39–50. [Google Scholar]
  33. Singh, H.; Srivastava, H. Numerical investigation of the fractional-order Liénard and Duffing equations arising in oscillating circuit theory. Front. Phys. 2020, 8, 120. [Google Scholar] [CrossRef]
  34. Scarfone, A.M.; Barbero, G.; Evangelista, L.R.; Lenzi, E.K. Anomalous Diffusion and Surface Effects on the Electric Response of Electrolytic Cells. Physchem 2022, 2, 163–178. [Google Scholar] [CrossRef]
  35. Barbero, G.; Evangelista, L.; Lenzi, E.K. Time-fractional approach to the electrochemical impedance: The Displacement current. J. Electroanal. Chem. 2022, 920, 116588. [Google Scholar] [CrossRef]
  36. Singh, H.; Srivastava, H.; Nieto, J.J. Handbook of Fractional Calculus for Engineering and Science; CRC Press: Boca Raton, FL, USA, 2022. [Google Scholar]
  37. Singh, H.; Kumar, D.; Baleanu, D. Methods of Mathematical Modelling: Fractional Differential Equations; CRC Press: Boca Raton, FL, USA, 2019. [Google Scholar]
  38. Singh, H.; Kumar, D.; Baleanu, D. Methods of Mathematical Modelling: Infectious Disease; Elsevier Science: Amsterdam, The Netherlands, 2022. [Google Scholar]
  39. Gómez-Aguilar, J.; Atangana, A. Fractional Hunter-Saxton equation involving partial operators with bi-order in Riemann–Liouville and Liouville–Caputo sense. Eur. Phys. J. Plus 2017, 132, 100. [Google Scholar] [CrossRef]
  40. Evangelista, L.R.; Lenzi, E.K. An Introduction to Anomalous Diffusion and Relaxation; Springer Nature: Berlin, Germany, 2023. [Google Scholar]
  41. Ali, I.; Kalla, S. A generalized Hankel transform and its use for solving certain partial differential equations. Anziam J. 1999, 41, 105–117. [Google Scholar] [CrossRef] [Green Version]
  42. Garg, M.; Rao, A.; Kalla, S.L. On a generalized finite Hankel transform. Appl. Math. Comput. 2007, 190, 705–711. [Google Scholar] [CrossRef]
  43. Nakhi, Y.B.; Kalla, S.L. Some boundary value problems of temperature fields in oil strata. Appl. Math. Comput. 2003, 146, 105–119. [Google Scholar] [CrossRef]
  44. Xie, K.; Wang, Y.; Wang, K.; Cai, X. Application of Hankel transforms to boundary value problems of water flow due to a circular source. Appl. Math. Comput. 2010, 216, 1469–1477. [Google Scholar] [CrossRef]
  45. O’Shaughnessy, B.; Procaccia, I. Diffusion on fractals. Phys. Rev. A 1985, 32, 3073–3083. [Google Scholar] [CrossRef]
  46. O’Shaughnessy, B.; Procaccia, I. Analytical Solutions for Diffusion on Fractal Objects. Phys. Rev. Lett. 1985, 54, 455–458. [Google Scholar] [CrossRef]
  47. Richardson, L.F. Atmospheric diffusion shown on a distance-neighbour graph. Proc. Math. Phys. Eng. Sci. 1926, 110, 709–737. [Google Scholar]
  48. Boffetta, G.; Sokolov, I.M. Relative Dispersion in Fully Developed Turbulence: The Richardson’s Law and Intermittency Corrections. Phys. Rev. Lett. 2002, 88, 094501. [Google Scholar] [CrossRef] [Green Version]
  49. Ben Avraham, D.; Havlin, S. Diffusion and Reactions in Fractals and Disordered Systems; CUP: Cambridge, UK, 2000. [Google Scholar]
  50. Su, N.; Sander, G.; Liu, F.; Anh, V.; Barry, D. Similarity solutions for solute transport in fractal porous media using a time- and scale-dependent dispersivity. App. Math. Model. 2005, 29, 852–870. [Google Scholar] [CrossRef] [Green Version]
  51. He, J.H. Fractal calculus and its geometrical explanation. Results Phys. 2018, 10, 272–276. [Google Scholar] [CrossRef]
  52. Cai, W.; Chen, W.; Xu, W. The fractal derivative wave equation: Application to clinical amplitude/velocity reconstruction imaging. J. Acoust. Soc. Am. 2018, 143, 1559–1566. [Google Scholar] [CrossRef]
  53. Chen, W.; Liang, Y. New methodologies in fractional and fractal derivatives modeling. Chaos Solitons Fractals 2017, 102, 72–77. [Google Scholar] [CrossRef]
  54. Liang, Y.; Chen, W.; Cai, W. Hausdorff Calculus: Applications to Fractal Systems; Walter de Gruyter GmbH & Co KG: Berlin, Germany, 2019; Volume 6. [Google Scholar]
  55. Evans, M.R.; Majumdar, S.N. Diffusion with stochastic resetting. Phys. Rev. Lett. 2011, 106, 160601. [Google Scholar] [CrossRef] [PubMed]
  56. Lenzi, M.K.; Lenzi, E.K.; Guilherme, L.; Evangelista, L.R.; Ribeiro, H.V. Transient anomalous diffusion in heterogeneous media with stochastic resetting. Phys. A Stat. Mech. Appl. 2022, 588, 126560. [Google Scholar] [CrossRef]
  57. Wyld, H.W.; Powell, G. Mathematical Methods for Physics; CRC Press: Boca Raton, FL, USA, 2020. [Google Scholar]
  58. Mathai, A.M.; Saxena, R.K.; Haubold, H.J. The H-Function: Theory and Applications; Springer Science & Business Media: Berlin, Germany, 2009. [Google Scholar]
  59. Lenzi, E.K.; Evangelista, L.; Lenzi, M.K.; Ribeiro, H.V.; de Oliveira, E.C. Solutions for a non-Markovian diffusion equation. Phys. Lett. A 2010, 374, 4193–4198. [Google Scholar] [CrossRef]
  60. Jiang, X.; Xu, M. The time fractional heat conduction equation in the general orthogonal curvilinear coordinate and the cylindrical coordinate systems. Phys. A Stat. Mech. Appl. 2010, 389, 3368–3374. [Google Scholar] [CrossRef]
  61. Evans, M.R.; Majumdar, S.N.; Schehr, G. Stochastic resetting and applications. J. Phys. A 2020, 53, 193001. [Google Scholar] [CrossRef]
Figure 1. Trend of G + ( x , x , t ) and G ( x , x , t ) obtained from Equation (22) for μ = 2 and different values of γ and η . We consider, for illustrative purposes, χ 2 = 1 , x = 1 , and t = 1 . Note that in (a,b) show that G + ( x , x , t ) and G ( x , x , t ) have different behavior, in particular, near the origin.
Figure 1. Trend of G + ( x , x , t ) and G ( x , x , t ) obtained from Equation (22) for μ = 2 and different values of γ and η . We consider, for illustrative purposes, χ 2 = 1 , x = 1 , and t = 1 . Note that in (a,b) show that G + ( x , x , t ) and G ( x , x , t ) have different behavior, in particular, near the origin.
Fluids 08 00034 g001
Figure 2. Trend of G + ( x , x , t ) and G ( x , x , t ) obtained from Equation (27) for different values of μ , γ and η , with p ( μ ¯ ) = χ μ δ ( μ ¯ μ ) + χ 2 δ ( μ ¯ 2 ) . We consider, for illustrative purposes, χ 2 = 1 , x = 1 , and t = 1 . Note that in (a,b) show that G + ( x , x , t ) and G ( x , x , t ) have different behavior, in particular, near the origin.
Figure 2. Trend of G + ( x , x , t ) and G ( x , x , t ) obtained from Equation (27) for different values of μ , γ and η , with p ( μ ¯ ) = χ μ δ ( μ ¯ μ ) + χ 2 δ ( μ ¯ 2 ) . We consider, for illustrative purposes, χ 2 = 1 , x = 1 , and t = 1 . Note that in (a,b) show that G + ( x , x , t ) and G ( x , x , t ) have different behavior, in particular, near the origin.
Fluids 08 00034 g002
Figure 3. Trend of G ( x , x , t ) obtained from Equations (21) and (22) for different values of μ , γ and η , with p ( μ ¯ ) = χ μ δ ( μ ¯ μ ) . We consider, for illustrative purposes, χ μ = 1 , x = 1 , and t = 1 .
Figure 3. Trend of G ( x , x , t ) obtained from Equations (21) and (22) for different values of μ , γ and η , with p ( μ ¯ ) = χ μ δ ( μ ¯ μ ) . We consider, for illustrative purposes, χ μ = 1 , x = 1 , and t = 1 .
Fluids 08 00034 g003
Figure 4. Trend of the mean square displacement σ x 2 versus t for different values of γ , μ , and η . Note that depending on the value of the parameter, different behaviors can be obtained.
Figure 4. Trend of the mean square displacement σ x 2 versus t for different values of γ , μ , and η . Note that depending on the value of the parameter, different behaviors can be obtained.
Fluids 08 00034 g004
Figure 5. Trend of G ( x , x , t ) obtained from Equations (21) and (22) for different values of μ , γ and η , with p ( μ ¯ ) = χ μ δ ( μ ¯ μ ) + χ 2 δ ( μ ¯ 2 ) . We consider, for illustrative purposes, χ μ = 1 , χ 2 = 0.6 , x = 1 , and t = 1 .
Figure 5. Trend of G ( x , x , t ) obtained from Equations (21) and (22) for different values of μ , γ and η , with p ( μ ¯ ) = χ μ δ ( μ ¯ μ ) + χ 2 δ ( μ ¯ 2 ) . We consider, for illustrative purposes, χ μ = 1 , χ 2 = 0.6 , x = 1 , and t = 1 .
Fluids 08 00034 g005
Figure 6. (a) shows the profile of ρ ( x , t ) obtained from Equations (15), (21) and (22) for different values of μ , γ and η , with p ( μ ¯ ) = χ μ δ ( μ ¯ μ ) , t = 1 , and ρ ( x , 0 ) = δ ( x ) . (b) shows the behavior of 1 / [ ρ ( 0 , t ) ] 2 to illustrate the spreading of the system. We consider, for illustrative purposes, χ μ = 1 and t = 1 .
Figure 6. (a) shows the profile of ρ ( x , t ) obtained from Equations (15), (21) and (22) for different values of μ , γ and η , with p ( μ ¯ ) = χ μ δ ( μ ¯ μ ) , t = 1 , and ρ ( x , 0 ) = δ ( x ) . (b) shows the behavior of 1 / [ ρ ( 0 , t ) ] 2 to illustrate the spreading of the system. We consider, for illustrative purposes, χ μ = 1 and t = 1 .
Fluids 08 00034 g006
Figure 7. (a) shows the profile of ρ ( x , t ) obtained from Equations (31) and (33) for different values of μ , γ and η , with p ( μ ¯ ) = χ μ δ ( μ ¯ μ ) , t = 1 , and ρ ( x , 0 ) = δ ( x ) . (b) shows the behavior of 1 / [ ρ ( 0 , t ) ] 2 to illustrate the spreading of the system. We consider, for illustrative purposes, χ μ = 1 and r = 1 .
Figure 7. (a) shows the profile of ρ ( x , t ) obtained from Equations (31) and (33) for different values of μ , γ and η , with p ( μ ¯ ) = χ μ δ ( μ ¯ μ ) , t = 1 , and ρ ( x , 0 ) = δ ( x ) . (b) shows the behavior of 1 / [ ρ ( 0 , t ) ] 2 to illustrate the spreading of the system. We consider, for illustrative purposes, χ μ = 1 and r = 1 .
Fluids 08 00034 g007
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Lenzi, E.K.; Somer, A.; Zola, R.S.; da Silva, L.R.; Lenzi, M.K. A Generalized Diffusion Equation: Solutions and Anomalous Diffusion. Fluids 2023, 8, 34. https://doi.org/10.3390/fluids8020034

AMA Style

Lenzi EK, Somer A, Zola RS, da Silva LR, Lenzi MK. A Generalized Diffusion Equation: Solutions and Anomalous Diffusion. Fluids. 2023; 8(2):34. https://doi.org/10.3390/fluids8020034

Chicago/Turabian Style

Lenzi, Ervin K., Aloisi Somer, Rafael S. Zola, Luciano R. da Silva, and Marcelo K. Lenzi. 2023. "A Generalized Diffusion Equation: Solutions and Anomalous Diffusion" Fluids 8, no. 2: 34. https://doi.org/10.3390/fluids8020034

APA Style

Lenzi, E. K., Somer, A., Zola, R. S., da Silva, L. R., & Lenzi, M. K. (2023). A Generalized Diffusion Equation: Solutions and Anomalous Diffusion. Fluids, 8(2), 34. https://doi.org/10.3390/fluids8020034

Article Metrics

Back to TopTop